Long non-coding RNA NCK1-AS1 promotes the tumorigenesis of glioma through sponging microRNA-138-2-3p and activating the TRIM24/Wnt/β-catenin axis
Journal of Experimental & Clinical Cancer Research volume 39, Article number: 63 (2020)
Glioma is a common brain malignancy with high mortality. The competing endogenous RNA (ceRNA) networks may play key roles in cancer progression. This study was conducted to probe the role of long noncoding RNA (lncRNA) NCK1-AS1 in glioma progression and the involved mechanisms.
Microarray analyses were performed to explore the lncRNAs/miRNAs/genes with differential expression in glioma. NCK1-AS1 levels in glioma tissues and normal brain tissues, and in glioma cell lines and normal human glial cells were identified. The interactions among NCK1-AS1, miR-138-2-3p and TRIM24 were validated through luciferase reporter, RNA immunoprecipitation and RNA pull-down assays. Gain- and loss-of functions of NCK1-AS1, miR-138-2-3p and TRIM24 were performed to identify their roles in the behaviors of glioma cells. The activity of the Wnt/β-catenin pathway was measured. In vivo experiments were performed as well.
High expression of NCK1-AS1 was found in glioma tissues and cells, especially in U251 cells. Online predictions and the integrated experiments identified that NCK1-AS1 elevated the TRIM24 expression through sponging miR-138-2-3p, and further activated the Wnt/β-catenin pathway. Artificial silencing of NCK1-AS1 or up-regulation of miR-138-2-3p led to inhibited proliferation, invasion and migration but promoted cell apoptosis of U251 cells, while up-regulation of TRIM24 reversed these changes, and it activated the Wnt/β-catenin pathway. The in vitro results were reproduced in in vivo experiments.
Our study suggested that NCK1-AS1 might elevate TRIM24 expression and further activate the Wnt/β-catenin pathway via acting as a ceRNA for miR-138-2-3p. Silencing of NCK1-AS1 might inhibit the progression of glioma.
Gliomas are tumors that initiate from precursor or glial cells and include astrocytoma (including glioblastoma), oligodendroglioma, ependymoma, oligoastrocytoma (mixed glioma) and malignant glioma and account for about 25.5% of all primary brain and other central nervous system tumors and 80.8% of malignant tumors . The current clinical treatment of glioma includes surgery, chemotherapy, radiotherapy, targeted therapy and immunotherapy; however, the prognosis of glioma remains unfavorable with a high recurrence rate after initial treatment . Diffuse glioma comprises no more than 1% of all newly diagnosed cancers, though, it leads to high mortality and morbidity, especially the most lethal type, glioblastoma, which takes place with 70–75% of all glioma cases with a median overall survival of 14–17 months . Developing novel therapeutic options and identifying new molecular mechanisms are of great importance in glioma control.
Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) are two large classes of non-protein-coding transcripts that participate in diverse essential cellular processes through multiple mechanisms [4, 5]. LncRNAs are larger than 200 nucleotides, and their importance in gene expression, key cellular processes, metastasis and disease prognosis in cancer has been largely studied [6, 7], including in glioma [8, 9]. miRNAs are 20–25 nucleotides in length, and they down-regulate gene expression by binding with target mRNAs post-transcriptionally thus resulting in transcript degradation . miRNAs exert critical functions in fundamental cellular process such as cell proliferation, apoptosis, development and inflammation . Likewise, dysregulation of several miRNAs has been found to participate in human glioma progression [12, 13]. Importantly, a competing endogenous RNA (ceRNA) theory proposed recently, which suggests that non-coding RNAs and protein-coding RNAs work as ceRNAs through competing for miRNAs via shared miRNA recognition elements [14, 15], has aroused wide concerns. Several ceRNA networks have been validated to play key functions in metastasis, proliferation, and apoptosis in human glioma cells [16, 17]. The possible ceRNA networks in glioma and the mechanisms involved remain largely unknown.
To this end, this study first figured out the representative differentially expressed lncRNAs through the glioma microarrays, and lncRNA NCK1-AS1, which has been noted to be aberrantly expressed in several human cancers [18, 19] with its role in glioma remaining unknown, was selected as the subject of this study. We further determined tripartite motif-containing 24 (TRIM24) as an indirect target of NCK1-AS1 through the miR-138-2-3p sponge. TRIM24 is a member of the TRIM family with its oncogenic role being identified in prostate cancer . In addition, TRIM24 has been suggested to promote cancer progression by triggering the Wnt/β-catenin signaling pathway [21, 22]. This pathway has been noted to be closely associated with cellular proliferation, migration, and angiogenesis and the following glioma malignancy . We assumed that NCK1-AS1 could affect glioma progression through the miR-138-2-3p/TRIM24 network and the following Wnt/β-catenin pathway, with gain- and loss-of functions of these molecules performed in both cell and animal experiments to validate this hypothesis.
Materials and methods
The glioma gene microarrays (GSE50161 and GSE35493) and miRNA microarray (GSE65626) acquired from the GEO Database (https://www.ncbi.nlm.nih.gov/geo/) were used for differential expression analysis utilizing an R-Language Limma Package . Differentially expressed lncRNAs/miRNAs/genes were analyzed with |log2FC| > 1.0 and p < 0.05 as the criteria, and the heatmaps for differentially expressed lncRNAs/miRNAs/genes were drawn by an online pheatmap package (https://cran.r-project.org/web/packages/pheatmap/index.html). The binding sites between lncRNAs and miRNAs were predicted on the RNA22 database (https://cm.jefferson.edu/rna22/). Besides, the RNA22, miRDB (http://www.mirdb.org), miRDIP (http://ophid.utoronto.ca/mirDIP/), miRWalk (http://mirwalk.umm.uni-heidelberg.de/) and DIANA (http://diana.imis.athena-innovation.gr/DianaTools/index.php?r=microT_CDS/index) databases were applied to predict the target genes of miRNAs. The prediction outcomes were integrated and analyzed using a Venn diagram produced by online webtools (http://bioinformatics.psb.ugent.be/webtools/Venn/).
Clinical tissue sample collection
A total of 12 pairs of normal brain tissue samples and 32 pairs of glioma tissue samples were acquired from the Zhejiang Provincial Hospital of Traditional Chinese Medicine. Normal brain tissues collected from epileptic patients who subjected to surgery were used as control. The study gained the approval of the Clinical Ethical Committee of the Zhejiang Provincial Hospital of Traditional Chinese Medicine. All procedures were conducted as per the Declaration of Helsinki. All eligible participants signed the informed consent.
RNA-in situ hybridization (RNA-ISH)
RNA-ISH of lncRNA NCK1-AS1 was performed using an ISH assay kit (Boye Biotechnology Co., Ltd., Guangzhou, Guangdong, China). In brief, the tissue samples were fixed, embedded in paraffin, and warm-incubated in gradient alcohol and then in 3% H2O2 for 30 min. Then, the streptavidin-horseradish peroxidase (HRP) conjugate and biotin conjugate probes were introduced into the samples for hybridization. Then samples were then stained with hematoxylin and observed under an optical microscope (Leica, Solms, Germany).
Normal human glial cell line HEB from Yan-Yu Bio-technology Co., Ltd. (Shanghai, China) (http://www.hdbsw.com/) and 4 glioma cell lines U251, SHG-441, U87 and T98 (ATCC, Manassas, USA) were incubated in DMEM (Gibco, NY, USA) supplemented with 10% fetal bovine serum (FBS), penicillin (100 U/mL) and streptomycin (100 mg/mL) (37 °C with 5% CO2). The cells were passaged when they reached an 85% confluence.
Cell treatment and grouping
The glioma cells were trypsinized to 2 × 106 cells/mL. Then the cell suspension was sorted into 12-well plates at 1 mL per well and incubated in 5% CO2 at 37 °C. Next, the cells were transfected with different vectors using the Lipofectamine 2000 (Thermo Fisher, MA, USA) as per the protocols and allocated into the corresponding groups. Three batches of grouping were performed. As for the first batch, cells were allocated into over-expression (oe)-negative control (NC) group (oe-NC group, cells were transfected with NC of NCK1-AS1 oe vector), oe-NCK1-AS1 group (cells were transfected with NCK1-AS1 oe vector), short hairpin (sh)-NC group (cells were transfected with NC of shRNA of NCK1-AS1) and sh-NCK1–AS1 (cells were transfected with sh-NCK1-AS1).
As for the second batch, cells were allocated into mimic NC group (cells were transfected with mimic NC), miR-138-2-3p mimic group (cells were transfected with miR-138-2-3p mimic), inhibitor NC group (cells were transfected with inhibitor NC) and miR-138-2-3p inhibitor group (cells were transfected with miR-138-2-3p inhibitor).
In terms of the third batch, cells were allocated into mimic NC group, miR-138-2-3p mimic group, inhibitor NC group, miR-138-2-3p inhibitor group, miR-138-2-3p mimic + oe-NC group (cells were transfected with miR-138-2-3p mimic and NC of NCK1-AS1 oe vector) and miR-138-2-3p mimic + oe-TRIM24 group (cells were transfected with miR-138-2-3p mimic and TRIM24 oe vector).
After transfection, the cells were incubated in the transfection solution for 4 h and then in normal cell culture medium for following experiments. The sequences for the vectors or mimic are listed in Supplementary Table 1.
Dual luciferase reporter gene assay
The wide type (WT) sequence based on the binding site between NCK1-AS1 and miR-138-2-3p and the corresponding mutant type (MUT) sequence were inserted into the target sequence of the psiCheck2 vector to construct psiChech2-NCK1-AS1-WT vector and psiChech2-NCK1-AS1-MUT vectors. The psiChech2-TRIM24-WT psiChech2-TRIM24-MUT vectors were constructed in a similar manner. Next, well-constructed reporter vectors were co-transfected with miR-138-2-3p mimic or mimic NC into HEK293T cells, and then the luciferase activity was determined in accordance with the instructions of a dual luciferase reporter assay system (Promega, Madison, WI, USA). Three independent experiments were performed.
5-ethynyl-2′-deoxyuridine (EdU) labeling assay
Exponentially growing cells were sorted into 24-well plates, and to each group, 3 duplicated cells were set up. Next, EdU (RiboBio Co., Ltd., Guangdong, China) was filled into the culture medium and adjusted to 10 μmol/mL. After 2 h of incubation in 5% CO2 at 37 °C, the medium was absorbed, and the cells were fixed in 4% paraformaldehyde-supplemented phosphate buffer saline (PBS) for 15 min, washed twice in 3% bovine serum albumin (BSA)-supplemented PBS, and then incubated in PBS supplemented with 0.5% Triton-100 for 20 min. Following two 3%BSA-PBS washes, each well was filled with 100 μL Apollo® 567 (RiboBio) and incubated for 30 min without light exposure at room temperature. Then the cells were washed twice in 3% BSA-PBS again, and stained with 1× Hoechst 33342 for 30 min. Thereafter, the wells were washed 3 more times in PBS, sealed, and observed under a fluorescence microscope (Olympus Optical Co., Ltd., Tokyo, Japan) with 5 views randomly selected. The positive cells, which were stained in red under the scope, were counted and recorded. The procedures were repeated 3 times.
Forty-eight hours after transfection, each group of cells were sorted into 6-well plates at 5 × 105 cells/well. A 200 μL pipette tip (21–0200, Biologix, Shanghai, China) was used to produce a scratch through the midline of the wells when the confluence reached 85%. Then the floating cells were washed away with PBS, while the remaining cells were continually incubated in serum-free medium for 1 h of recovery. At 0 h and 24 h after recovery, the cells were photographed to measure the migration of cells on an Image-Pro Plus Anaysis software (Media Cybernetics, USA). The experiment was performed in triplicate.
Extracellular matrix (ECM) gel was allowed to stand at 4 °C overnight. The next day (all the pipette tips and Transwells were pre-cooled on ice for half an hour), the gel was diluted in serum-free medium at 1:9 till the final concentration reaching 1 mg/mL. Each apical chamber of the 24-well Transwells was loaded with 40 μL ECM gel on the polycarbonate membrane and then placed in a 37 °C incubator with 5% CO2 for 5 h to polymerize the ECM gel. Next, the remaining liquid was removed, and each well was filled with 70 μL DMEM and incubated at 37 °C with 5% CO2 for 30 min to hydrate the gel again with the remaining culture medium discarded. After that, each group of cells were hungered in a serum-free condition for 24 h, detached, centrifuged, and diluted in FBS-free DMEM to 2.5 × 105 cells/mL. Then 0.2 mL suspension was loaded in each apical chamber where the substrate membraned was already hydrated, while each basolateral chamber was loaded with 700 μL pre-cooled 10% FBS-DMEM. The chambers were then incubated in 5% CO2 with saturated humidity at 37 °C for 24 h. Thereafter, the apical chambers were removed, and the cells in apical chambers and on the substrates were discarded using wet cotton swabs. The remaining cells were methanol-fixed, stained by 0.1% crystal violet, air dried, and photographed under a microscope with 5 random fields (200 ×) observed, and the average volume of invaded cells was calculated. The experiments were performed in triplicate.
Cells were washed 3 times in PBS 48 h after transfection and centrifuged at 3000 r/min for 20 min to discard the supernatant, and then diluted in PBS to adjust the concentration to 1 × 105 cells/mL. Then the suspension was successively treated with 1 mL − 20 °C pre-cooled 75% ethanol for 1 h of fixing, centrifuged for 5 min at 1500 r/min, washed in PBS, treated with 100 μL Rnase A (Thermo Fisher) without light exposure, bathed in 37 °C water for 30 min, and stained by 400 μL propidium iodide (PI) (Sigma-Aldrich, Merck KGaA, Darmstadt, Germany) at 4 °C for 30 min. Next, the cell cycle was measured based on the red fluorescence at the 488 nm (excitation wavelength).
In terms of apoptosis detection, the cells were detached in EDTA-free trypsin (Thermo Fisher) 48 h after transfection and filled into the flow cytometer tubes. Following 30 min of centrifugation at 3000 r/min and 3 cold PBS washes, the cells were further centrifuged at 3000 r/min for 20 min to discard the supernatant. Then the Annexin-V-FITC/PI staining solution was compounded by HEPES buffer (Thermo Fisher), Annexin-V-FITC and PI (50:1:2). Then the cells were stained for 15 min at room temperature, and then treated with 1 mL HEPES buffer. After vibration, the cell apoptosis was measured on the flow cytometer at 488 nm.
Fluorescence in situ hybridization (FISH)
The sub-cellular localization of NCK1-AS1 in glioma cell line was assessed according to the protocols of a lncRNA FISH probe Mix (RiboBio). In brief, cover glasses were put into 6-well plates on which glioma cells were seeded. One day later when the cell confluence got to 80%, the cell slides were collected, washed in PBS, and treated with 1 mL 4% paraformaldehyde. Next, the cells were given Protease-K (2 μg/mL), glycine and acetylating agent, and then cultured with prehybridization agent (250 μL) at 42 °C for 1 h, and then with hybridization agent (250 μL) containing probes (300 ng/mL) at 42 °C overnight. After 3 PBST washes, the cells were treated with PBST-DAPI (1:800) for nucleus staining for 5 min. Next, the cell slides were washed in PBST (3 × 3 min), sealed with anti-fluorescence quencher, and then observed under the fluorescence microscope (× 400) with 5 random fields included.
RNA immunoprecipitation (RIP)
A RIP kit (Millipore Corp., Billerica, MA, USA) was applied to measure the binding relation between NCK1-AS1 and AGO2. Glioma cells were washed in PBS. Then the cells were lysed in an equal volume of RIPA cell lysis solution (Beyotime, Shanghai, China) on ice for 5 min, and then centrifuged at 14000 r/min at 4 °C to collect the supernatant. A part of cell lysis was used as Input, and another part was cultured with the antibodies for coprecipitation. In each coprecipitation system, 50 μL magnetic beads were resuspended in RIP wash buffer (100 μL) and then incubated with 5 μg antibodies according the grouping. The magnet bead-antibody compounds were resuspended in RIP wash buffer (900 μL), and further cultured with 100 μL cell lysis overnight at 4 °C. The samples were set on the magnet base to obtain the magnet bead-antibody compounds. The samples and Input were respectively detached with protease K with the RNA extracted for the following reverse transcription quantitative polymerase chain reaction (RT-qPCR). The antibodies are AGO2 (ab32381, 1:50) and immunoglobulin G (IgG, 1:100, A109489) (all provided by Abcam Inc., Cambridge, MA, USA).
Glioma cells were either transfected with MUT-biotinylated or WT-biotinylated miR-138-2-3p (50 nM for each). Cells were harvested 48 h after, washed in PBS, and whirled. Then the cells were cultured for 10 min in specific cell lysis solution (Ambion, Austin, Texas, USA). The cell lysates were co-cultured with RNase-free and yeast tRNA-precoated M-280 streptavidin beads (all from Sigma-Aldrich) for 3 h at 4 °C and then washed twice in cold lysis solution, 3 times in low salt buffer, and once in high salt buffer. The combined RNA was purified with Trizol, and the NCK1-AS1 expression was evaluated using RT-qPCR.
Western blot analysis
Glioma cells and tissues were lysed in protease inhibitor-supplemented enhanced RIPA lysis solution (Boster, Wuhan, China) and then to have the total protein collected, and a bicinchoninic acid kit (Boster) was applied to assess the protein concentration. The proteins were run on 10% SDS-PAGE and transferred onto PVDF membranes. Then the membranes were sealed in 5% BSA for 2 h to block non-specific binding, and then cultured with the primary antibodies (Table 1) at 4 °C overnight. Afterwards, the membranes were washed 3 times in PBST and incubated with horseradish peroxidase-labeled goat-anti-rabbit (ab205719, 1:2000, Abcam) for 1 h at room temperature. Then the membranes were further washed in PBST and measured using the enhanced chemiluminescence reagent (EMD Millipore, USA). The signal intensity of the protein bands was analyzed using Image J software with the value of GAPDH set as the internal reference. The procedures were conducted for 3 times.
The Trizol reagent (15,596,026, Invitrogen Inc., Carlsbad, CA, USA) was utilized to obtain the total RNA of cells and tissues. Reverse transcription was performed using a PrimeScript RT reagent kit (RR047A, Takara Holdings Inc., Tokyo, Japan) to produce cDNA. Quantification was performed as per the instructions of a RT-qPCR assay kit (Thermo Fisher). The relative expression of miR-148a was normalized by a U6 transcript. Three duplicate wells were set for each group. U6 was set as the internal reference for miR-138-2-3p while GAPDH for other genes. Relative gene expression was evaluated using the 2-ΔΔCt method. The primers were synthetized by TransGen Biotech (Beijing, China) and the sequences are presented in Table 2.
The expression of stemness-related proteins SOX2 (1:1000, ab97959, Abcam) and OCT4 (1:1000, ab18976, Abcam) in U251 and SHG441 cells was determined as guided by a previous study .
Tumor sphere formation assay
U251 and SHG-441 cells were seeded into 6-well ultra-low attachment culture dishes (Corning Incorporated, Corning, NU, USA). Then the cells were cultured in serum-free DMEM/F12 supplemented with 1% B27, 20 ng/mL human epidermal growth factor (EGF) and 20 ng/mL human fibroblast growth factor (FGF) at 37 °C for 2 weeks to form tumor spheres. Two weeks later, the number of newly formed tumor spheres were counted under an optical microscope (Nikon Corporation, Tokyo, Japan) at a × 200 magnification.
Xenograft tumors in nude mice
Thirty female BALB/c nude mice (3–4 weeks old, 14 ± 2 g, provided by SJA Laboratory Animal Co., Ltd., Hunan, China) were fed at 25–27 °C with controlled humidity (45–50%) with free access to food and water. Then the mice were randomized into 6 groups: sh-NC group, sh-NCK1-AS1 group, mimic NC group, miR-138-2-3p mimic group, sh-NCK1-AS1 + oe-NC group and sh-NCK1-AS1 + oe-TRIM24 group, 5 mice in each. Cells with stable transfection were resuspended and adjusted to 1 × 107 cells/mL, and each nude mouse was implanted with 20 μL cell suspension. The tumor growth was observed and photographed every 7 d for a total of 35 d, and a tumor growth curve was drawn. The tumor volume was calculated as (a × b2)/2, in which “a” refers to the longest diameter while “b” refers to the shortest. The mice were euthanized on the 35th d via intraperitoneal injection of 150 mg/kg pentobarbital to collect and weigh the tumors. Thereafter, the tumor samples were preserved in liquid nitrogen for following RT-qPCR and western blot analysis. Animal studies were conducted as per the principles and procedures ratified by the Committee on the Ethics of Animal Experiments of the Zhejiang Provincial Hospital of Traditional Chinese Medicine. Great attempts were made to minimize the number and pain of animals.
SPSS 21.0 (IBM Corp. Armonk, NY, USA) was applied for data analysis. Data were in normal distribution. Measurement data were exhibited as mean ± standard derivation (mean ± SD). Differences between every two groups were analyzed using the t test, while those among multiple groups were analyzed using one-way analysis of variance (ANOVA) and Tukey’s multiple comparison test. Differences among multiple groups at different time points were analyzed using repeated measurement ANOVA and Bonferroni test. The Log-rank test was used for single-factor analysis. The p value was acquired from two-tailed tests, and p < 0.05 was considered to show significant difference.
NCK1-AS1 is highly expressed in PC tissues
Differentially expressed lncRNAs in the glioma microarrays GSE50161 and GSE35493 were analyzed, and the heatmap for top 10 highly expressed lncRNAs in the microarray GSE50161 was produced (Fig. 1a). And we chose NCK1-AS1, whose oncogenic roles have been noted in many other cancer types excluding glioma, as the subject of this study. Besides, high expression of NCK1-AS1 was also found in the other microarray GSE35493 (Fig. 1b). The online system GEPIA (http://gepia.cancer-pku.cn/) suggested that NCK1-AS1 is highly expressed in the TCGA-GBM DataBase (Fig. 1c). Moreover, in line with the prediction results, RT-qPCR identified high NCK1-AS1 expression in glioma tissues (p < 0.05) (Fig. 1d). Still, the NCK1-AS1 expression was increased in glioma cell lines U251, SHG-441, U87 and T98 compared to the normal glial cell line HEB, with the highest expression found in the U251 cell line while the lowest in the SHG-441 cell line (all p < 0.05) (Fig. 1e). Then, the U251 and SHG-441 cell lines were used for the subsequent experiments. In addition, RNA-ISH was further performed to validate NCK1-AS1 expression in brain and glioma tissues, and the results also suggested a higher NCK1-AS1 level in glioma tissues than in brain tissues (Fig. 1f).
Artificial silencing of NCK1-AS1 was introduced in U251 cells while over-expression of NCK1-AS1 was introduced in SHG-441 cells. The transfection was successfully performed based on the NCK1-AS1 expression in cells measured by RT-qPCR, and the sh-NCK1-AS1–1 presented the best silencing efficacy (all p < 0.05) (Fig. 2a). The EdU assasy, scratch test and Transwell assay results suggested that silencing of NCK1-AS1 led to inhibited proliferation, migration and invasion of U251 cells, while over-expression of NCK1-AS1 resulted in opposite trends in SHG-441 cells (all p < 0.05) (Fig. 2b-d). Moreover, the flow cytometry results suggested that silencing of NCK1-AS1 led to increased cell ratio in the G0/G1 phase while decreased cells in the S phase, correspondingly the cell apoptosis rate was increased in U251 cells. Conversely, the cell cycle arrest in G0/G1 phase and the apoptosis of SHG-441 cells with over-expressed NCK1-AS1 were decreased (all p < 0.05) (Fig. 2e-f). These results indicated that silencing of NCK1-AS1 might inhibit proliferation, invasion, migration and the resistance to death of glioma cells.
NCK1-AS1 might bind to miR-138-2-3p to regulate TRIM24
A heatmap for top 10 differentially expressed miRNAs according the data of GSE65626 microarray were drawn (Fig. 3a), and it was suggested that Homo sapiens (has)-miR-138-2-3p was poorly expressed in glioma tissues. A total of 230 potential miRNAs that may sponged with NCK1-AS1 were predicted via RNA22 with p < 0.05 as the selection criterion. Meanwhile, 27 miRNAs with low expression in the microarray GSE65626 were found, among which 3 miRNAs were found as the down-stream miRNAs of NCK1-AS1: hsa-miR-487a-5p, hsa-miR-184 and hsa-miR-138-2-3p (Fig. 3b). The binding sites between NCK1-AS1 and has-miR-138-2-3p predicted on RNA22 are shown in Fig. 3c, indicating NCK1-AS1 might sponge miR-138-2-3p. Moreover, the target genes of miR-138-2-3p were acquired in several databases (RNA22, miRDB, miRDIP, miRWalk and DIANA), and then a Venn diagram was produced (Fig. 3d) with 36 intersections found. The potential intersect genes were further compared with the highly expressed genes in the glioma microarrays GSE50161 and GSE35493, after which 2 intersections TRIM24 and GLISE were found (Fig. 3e). As aforementioned, TRIM24 was noted as an oncogene, but its role in glioma remains unknown. High expression of TRIM24 in microarrays GSE50161 and GSE35493 are shown in Fig. 3f and Fig. 3g. TRIM24 is also highly expressed in the TCGA-GBM DataBase (Fig. 3h). The biding site between miR-138-2-3p and TRIM24 predicted on RNA22 are shown in Fig. 3i. As mentioned above, TRIM24 has been documented to regulate the Wnt/β-catenin signaling pathway. Herein, we speculated that NCK1-AS1 might regulate TRIM24 expression via sponging miR-138-2-3p and further regulate the Wnt/β-catenin pathway in glioma.
The expression of miR-138-2-3p, TRIM24 and β-catenin in glioma and normal brain tissues was measured using RT-qPCR and western blot analysis. It was shown that miR-138-2-3p was lowly expressed, while the mRNA and protein expression of TRIM24 was elevated in glioma tissues (all p < 0.05) (Fig. 3j-k). The putative binding sites between NCK1-AS1 and miR-138-2-3p and between miR-138-2-3p and TRIM24 are shown in Fig. 3c and Fig. 3h. The dual luciferase reporter gene assay results suggested that NCK1-AS1 could directly bind to miR-138-2-3p, and miR-138-2-3p could bind to TRIM24 (all p < 0.05) (Fig. 3l-m).
NCK1-AS1 acts as a ceRNA for miR-138-2-3p in glioma cancer
We further detected the sub-cellular localization of NCK1-AS1 in U251 and SHG-441cells using FISH. As shown in Fig. 4a, in which the blue parts indicate nuclei while the red parts indicate NCK1-AS1, it was shown that NCK1-AS1 is mainly sub-localized in cytoplasm, indicating NCK1-AS1 might exert function through the ceRNA network. Following the findings above, RNA pull-down and RIP assays were further performed to identify the interactions among NCK1-AS1, miR-138-2-3p and TRIM24. The RNA pull-down assay found that compared to the MUT-miR-138-2-3p and the Bio-NC groups, the binding with NCK1-AS1 increased in the WT-miR-138-2-3p group (Fig. 4b). Meanwhile, the RIP assay results suggested that NCK1-AS1 could form compound with anti-AGO2 in U251 and SHG-441 cells (Fig. 4c). Since AGO2 is capable of forming compound with microRNAs, the fact that NCK1-AS1 combined with the AGO2-miR-1382-3p compound further identified the binding relationship between NCK1-AS1 and miR-138-2-3p.
Next, U251 cells with silenced NCK1-AS1 were further transfected with miR-138-2-3p inhibitor, while SHG-441 cells with over-expressed NCK1-AS1 were transfected with miR-138-2-3p mimic. RT-qPCR found that silencing of NCK1-AS1 increased miR-138-2-3p expression in U251 cells, leading to decreased TRIM24 expression. Then, miR-138-2-3p inhibitor transfection partly recovered the TRIM24 expression. In SHG-441 cells, overexpression of NCK1-AS1 led to decreased miR-138-2-3p expression and increased TRIM24 expression, but miR-138-2-3p mimic transfection significantly decreased TRIM24 expression (Figure 4d,e).
Up-regulation of miR-138-2-3p inhibits invasion and migration but promotes apoptosis of U251 cells via down-regulating TRIM24
Following the above findings, we further identified that inhibition of miR-138-2-3p or over-expression of TRIM24 aggregated the malignant behaviors of U251 cells with silenced NCK1-AS1, including increased cell proliferation, invasion and migration while reduced cell cycle arrest and cell apoptosis. Consistently, up-regulation of miR-138-2-3p inhibited the above malignant behaviors in SHG-441 cells with over-expressed NCK1-AS1 (Fig. 5a-e).
Meanwhile, western blot analysis was applied to measure the protein levels of Wnt/β-catenin pathway-related factor (Wnt1, Wnt3a and β-catenin), proliferation-related factor (PCDA), invasion and migration-related factors (N-cadherin and MMP), and apoptosis-related factors (Bcl-2 and Bax) in cells. The results told silencing of miR-138-2-3p or over-expression of TRIM24 increased the protein levels of Wnt1, Wnt3a, β-catenin, PCNA, N-cadherin, MMP-9 and Bcl-2 while inhibited the Bax expression. But over-expression of miR-138-2-3p in SHG-441 cells led to opposite trends (all p < 0.05) (Fig. 5f).
Silencing of NCK1-AS1 promotes the stemness of glioma cells
As one of the most malignant tumors, glioma is featured with strong cancer stem cell (CSC) characteristics. Here we determined the CSC-related biomarkers SOX2 and OCT4 in U251 and SHG-441 cells. The results showed that silencing of NCK1-AS1 inhibited the expression of SOX2 and OCT4 in cells, while the further silencing of miR-138-2-3p partly recovered the SOX2 and OCT4 expression. The reverse trends were found in SHG-441 cells with over-expressed NCK1-AS1 and further over-expressed miR-138-2-3p (all p < 0.05) (Fig. 6a). Moreover, the similar results were produced in tumor sphere formation assay (all p < 0.05) (Fig. 6b).
Silencing of NCK1-AS1 or up-regulation of miR-138-2-3p inhibits tumor formation in nude mice
Nude mice were implanted with cells with stable transfections to further investigate the effect of NCK1-AS1 on tumor formation in vivo. The tumor volume and weight were detected and recorded (Fig. 7a-c). It was found that silencing of NCK1-AS1 or up-regulation of miR-138-2-3p led to increased TRIM24 expression and decreased tumor formation speed in mice, while the further over-expression of TRIM24 promoted the tumor formation speed (all p < 0.05) (Figure 7 d,e).
It is quite necessary to develop new potential target for glioma treatment since the overall 5-year survival rate of patients with glioma remains lower than 5% even following comprehensive chemotherapy, radiotherapy and surgery . LncRNAs can sponge miRNAs and block these miRNAs from regulating the target mRNAs, thus mediating target genes post transcriptionally . In light of this theory and several recently identified ceRNA networks in cancers, the study figured out a novel NCK1-AS1/miR-138-2-3p/TRIM24 ceNRA network in glioma.
The evidence concerning the correlation between abnormal lncRNA expression with cancers has been growing [28, 29]. Initially, data on glioma microarrays suggested that NCK1-AS1 is aberrantly highly expressed in glioma, which was further identified by RT-qPCR. As for the other screened out lncRNAs, lncRNA CRNDE, for example, has already been noted to promote glioma cell growth , while LINC00998, has hardly been investigated in human cancers. NCK1-AS1 is a relatively newly recognized lncRNA with its oncogenic role in many several human malignancies [19, 31, 32] reported, thus, was selected as the current study subject. Then, we found that silencing of NCK1-AS1 led to decreased proliferation, invasion and migration abilities while increased cell cycle arrest and apoptosis of glioma cells. Moreover, silencing of NCK1-AS1 was also found to inhibit the expression of CSC biomarkers SOX2 and OCT4 in our study. As CSCs are a class of self-renewal cells with strong tumorigenic potency and are more resistant to conventional therapies than other cancer cells . This finding further evidenced the promoting role of NCK1-AS1 in the malignant behaviors of glioma cells. The same trends were reproduced in in vivo experiments, which showed that silencing of NCK1-AS1 inhibited tumor formation and growth in nude mice.
The findings above triggered us to identify the downstream mechanisms involved in the events. As above mentioned, lncRNAs might exert functions through the crosstalk with miRNAs and mRNAs [14, 15]. Therefore, we analyzed the potential target miRNAs of NCK1-AS1 and the following genes via integrated online bioinformation system, glioma microarrays, dual luciferase reporter gene, RNA pull-down and RIP assays, after which we identified the interactions among NCK1-AS1, miR-138-2-3p and TRIM24. Importantly, we noticed that up-regulation of miR-138-2-3p decreased proliferation, invasion and migration while promoted apoptosis of SHG-441 cells with over-expressed NCK1-AS1, showing as decreased levels of PCNA, N-cadherin, MMP-9 and Bcl-2 while increased level of Bax. But artificial over-expression of TRIM24 led to opposite trends in U251 cells with silenced NCK1-AS1. PCNA is an important replication accessory factor that supports DNA replication, repair and recombination, and cell cycle regulation [34, 35]. Abnormal expression of the N-cadherin is a crucial biomarker of epithelial-to-mesenchymal transition in many cancer types, promoting the aggressiveness of tumors . As one of the most studied MMPs, MMP-9 is well-known for the key roles in invasion of cancer cells and metastasis of tumors . Bcl-2 is an anti-apoptotic protein that mediates apoptosis by regulating the permeability of the mitochondrial membrane, while Bax can break the outer mitochondrial membrane thus promoting apoptosis . In addition, hsa-miR-487a-5p and hsa-miR-184 were also identified as NCK1-AS1 targets after integrated analysis. However, the role of miR-184 in glioma has been largely studied [39, 40], while the role of miR-487a-5p in cancer development has been little concerned. miR-138-2-3p is a rarely mentioned miRNA, though, some same trends were found in a previous study , which suggested that over-expression of miR-138-2-3p decreased proliferation and invasion but promoted apoptosis of human laryngeal CSCs following radiotherapy. Meanwhile, lowly-expressed miR-138-2-3p has been found in drug-resistant non-small cell lung cancer patients . As for TRIM24, it is a member of the tripartite motif (TRIM) family , and the oncogenic role of TRIM24 in carcinogenesis has been demonstrated in several cancer types [43,44,45]. Moreover, our study identified that inhibited miR-138-2-3p or promoted TRIM24 stimulated Wnt/β-catenin activation in glioma cells, which was in line with the previous studies [21, 41]. Wnt/β-catenin controls diverse cellular processes and drives cancer progression , with its activation revealed to promotes tumorigenesis of multiple cancers [47, 48], including glioma . Taken the above discussion together, in can be inferred that NCK1-AS1 could promote the malignant behaviors of glioma cells and tumor metastasis through the crosstalk with miR-138-2-3p and TRIM24, and the following activation of the Wnt/β-catenin pathway.
By the way, GLIS2, the other intersected target gene of miR-138-2-3p, is a member of GLI-similar zinc finger protein family and is closely linked with acute myeloid leukaemia, has been also noted to correlated with tumor progression [49, 50]. We would like to investigate if GLIS2 plays roles in glioma progression in our future experiments.
To sum up, this study identified a novel ceRNA network in glioma that NCK1-AS1 up-regulates TRIM24 expression through sponging miR-138-2-3p. Silencing of NCK1-AS1 might inhibit the proliferation, invasion, migration and resistance to death of glioma cells, and inhibit the metastasis of tumor in vivo through down-regulating TRIM24 and the following Wnt/β-catenin inactivation (Fig. 8). These findings might provide novel insights in glioma treatment, and we hope more studies in this field would be carried out to validate our findings, and to identify more potential mechanisms for better understanding of glioma.
Availability of data and materials
All the data generated or analyzed during this study are included in this published article.
Analysis of variance
Bovine serum albumin
Dulbecco’s modified eagle’s medium
Fetal bovine serum
Fluorescence in situ hybridization;
- mean ± SD:
mean ± standard derivation
Phosphate buffer saline
Proliferating cell nuclear antigen
reverse transcription quantitative polymerase chain reaction
Sodium dodecyl sulfate-polyacrylamide gel electrophoresis
Tripartite motif containing 24
Ostrom QT, Cioffi G, Gittleman H, Patil N, Waite K, Kruchko C, et al. CBTRUS Statistical Report: Primary Brain and Other Central Nervous System Tumors Diagnosed in the United States in 2012–2016. Neuro Oncol. 2019;21(Supplement_5):v1–v100.
Chen W, Lei C, Liu P, Liu Y, Guo X, Kong Z, et al. Progress and prospects of recurrent Glioma: a recent Scientometric analysis of the web of science in 2019. World Neurosurg. 2020;134:e387–e99.
Molinaro AM, Taylor JW, Wiencke JK, Wrensch MR. Genetic and molecular epidemiology of adult diffuse glioma. Nat Rev Neurol. 2019;15(7):405–17.
Wan ZY, Song F, Sun Z, Chen YF, Zhang WL, Samartzis D, et al. Aberrantly expressed long noncoding RNAs in human intervertebral disc degeneration: a microarray related study. Arthritis Res Ther. 2014;16(5):465.
Xue M, Zhuo Y, Shan B. MicroRNAs, long noncoding RNAs, and their functions in human disease. Methods Mol Biol. 2017;1617:1–25.
Feng S, Yao J, Chen Y, Geng P, Zhang H, Ma X, et al. Expression and functional role of reprogramming-related long noncoding RNA (lincRNA-ROR) in Glioma. J Mol Neurosci. 2015;56(3):623–30.
Hauptman N, Glavac D. Long non-coding RNA in cancer. Int J Mol Sci. 2013;14(3):4655–69.
Ma Y, Zhou G, Li M, Hu D, Zhang L, Liu P, et al. Long noncoding RNA DANCR mediates cisplatin resistance in glioma cells via activating AXL/PI3K/Akt/NF-kappaB signaling pathway. Neurochem Int. 2018;118:233–41.
Zhang D, Zhou H, Liu J, Mao J. Long noncoding RNA ASB16-AS1 promotes proliferation, migration, and invasion in Glioma cells. Biomed Res Int. 2019;2019:5437531.
Sun L, Zhao M, Wang Y, Liu A, Lv M, Li Y, et al. Neuroprotective effects of miR-27a against traumatic brain injury via suppressing FoxO3a-mediated neuronal autophagy. Biochem Biophys Res Commun. 2017;482(4):1141–7.
Dong RF, Zhang B, Tai LW, Liu HM, Shi FK, Liu NN. The Neuroprotective role of MiR-124-3p in a 6-Hydroxydopamine-induced cell model of Parkinson's disease via the regulation of ANAX5. J Cell Biochem. 2018;119(1):269–77.
Lv QL, Zhu HT, Li HM, Cheng XH, Zhou HH, Chen SH. Down-regulation of miRNA-320c promotes tumor growth and metastasis and predicts poor prognosis in human glioma. Brain Res Bull. 2018;139:125–32.
Zhou Q, Liu J, Quan J, Liu W, Tan H, Li W. MicroRNAs as potential biomarkers for the diagnosis of glioma: a systematic review and meta-analysis. Cancer Sci. 2018;109(9):2651–9.
Chen Y, Shen Z, Zhi Y, Zhou H, Zhang K, Wang T, et al. Long non-coding RNA ROR promotes radioresistance in hepatocelluar carcinoma cells by acting as a ceRNA for microRNA-145 to regulate RAD18 expression. Arch Biochem Biophys. 2018;645:117–25.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Jiang J, Wang X, Gao G, Liu X, Chang H, Xiong R, et al. Silencing of lncRNA HOXA11-AS inhibits cell migration, invasion, proliferation, and promotes apoptosis in human glioma cells via upregulating microRNA-125a: in vitro and in vivo studies. Am J Transl Res. 2019;11(10):6382–92.
Liang Q, Li X, Guan G, Xu X, Chen C, Cheng P, et al. Long non-coding RNA, HOTAIRM1, promotes glioma malignancy by forming a ceRNA network. Aging (Albany NY). 2019;11(17):6805–38.
Hu H, Li H, Feng X. Downregulation of lncRNA NCK1-AS1 inhibits Cancer cell migration and invasion in nasopharyngeal carcinoma by Upregulating miR-135a. Cancer Manag Res. 2019;11:10531–7.
Huang L, Gan X, He L, Wang L, Yu J. Silencing of long non-coding RNA NCK1-AS1 inhibits cell proliferation and migration via inhibition of microRNA-134 in cervical cancer. Exp Ther Med. 2019;18(3):2314–22.
Offermann A, Roth D, Hupe MC, Hohensteiner S, Becker F, Joerg V, et al. TRIM24 as an independent prognostic biomarker for prostate cancer. Urol Oncol. 2019;37(9):576 e1–e10.
Fang Z, Deng J, Zhang L, Xiang X, Yu F, Chen J, et al. TRIM24 promotes the aggression of gastric cancer via the Wnt/beta-catenin signaling pathway. Oncol Lett. 2017;13(3):1797–806.
Yuan YH, Zhou J, Zhang Y, Xu MD, Wu J, Li W, et al. Identification of key genes and pathways downstream of the beta-catenin-TCF7L1 complex in pancreatic cancer cells using bioinformatics analysis. Oncol Lett. 2019;18(2):1117–32.
He L, Zhou H, Zeng Z, Yao H, Jiang W, Qu H. Wnt/beta-catenin signaling cascade: a promising target for glioma therapy. J Cell Physiol. 2019;234(3):2217–28.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:3.
Su YK, Bamodu OA, Tzeng YM, Hsiao M, Yeh CT, Lin CM. Ovatodiolide inhibits the oncogenicity and cancer stem cell-like phenotype of glioblastoma cells, as well as potentiate the anticancer effect of temozolomide. Phytomedicine. 2019;61:152840.
Binder DC, Davis AA, Wainwright DA. Immunotherapy for cancer in the central nervous system: current and future directions. Oncoimmunology. 2016;5(2):e1082027.
Cheng Y, Geng L, Wang K, Sun J, Xu W, Gong S, et al. Long noncoding RNA expression signatures of Colon Cancer based on the ceRNA network and their prognostic value. Dis Markers. 2019;2019:7636757.
Li CH, Chen Y. Insight into the role of long noncoding RNA in Cancer development and progression. Int Rev Cell Mol Biol. 2016;326:33–65.
Lorenzi L, Avila Cobos F, Decock A, Everaert C, Helsmoortel H, Lefever S, et al. Long noncoding RNA expression profiling in cancer: challenges and opportunities. Genes Chromosomes Cancer. 2019;58(4):191–9.
Kiang KM, Zhang XQ, Zhang GP, Li N, Cheng SY, Poon MW, et al. CRNDE expression positively correlates with EGFR activation and modulates Glioma cell growth. Target Oncol. 2017;12(3):353–63.
Chang H, Li B, Zhang X, Meng X. NCK1-AS1 promotes NCK1 expression to facilitate tumorigenesis and chemo-resistance in ovarian cancer. Biochem Biophys Res Commun. 2019.
Li H, Jia Y, Cheng J, Liu G, Song F. LncRNA NCK1-AS1 promotes proliferation and induces cell cycle progression by crosstalk NCK1-AS1/miR-6857/CDK1 pathway. Cell Death Dis. 2018;9(2):198.
Najafi M, Mortezaee K, Majidpoor J. Cancer stem cell (CSC) resistance drivers. Life Sci. 2019;234:116781.
Park SY, Jeong MS, Han CW, Yu HS, Jang SB. Structural and functional insight into proliferating cell nuclear antigen. J Microbiol Biotechnol. 2016;26(4):637–47.
Strzalka W, Ziemienowicz A. Proliferating cell nuclear antigen (PCNA): a key factor in DNA replication and cell cycle regulation. Ann Bot. 2011;107(7):1127–40.
Mrozik KM, Blaschuk OW, Cheong CM, Zannettino ACW, Vandyke K. N-cadherin in cancer metastasis, its emerging role in haematological malignancies and potential as a therapeutic target in cancer. BMC Cancer. 2018;18(1):939.
Huang H. Matrix Metalloproteinase-9 (MMP-9) as a Cancer Biomarker and MMP-9 Biosensors: Recent Advances. Sensors (Basel). 2018;18(10):3249.
Lv J, Liang Y, Tu Y, Chen J, Xie Y. Hypoxic preconditioning reduces propofol-induced neuroapoptosis via regulation of Bcl-2 and Bax and downregulation of activated caspase-3 in the hippocampus of neonatal rats. Neurol Res. 2018;40(9):767–73.
Cheng Z, Wang HZ, Li X, Wu Z, Han Y, Li Y, et al. MicroRNA-184 inhibits cell proliferation and invasion, and specifically targets TNFAIP2 in Glioma. J Exp Clin Cancer Res. 2015;34:27.
Feng R, Dong L. Inhibitory effect of miR-184 on the potential of proliferation and invasion in human glioma and breast cancer cells in vitro. Int J Clin Exp Pathol. 2015;8(8):9376–82.
Zhu Y, Shi LY, Lei YM, Bao YH, Li ZY, Ding F, et al. Radiosensitization effect of hsa-miR-138-2-3p on human laryngeal cancer stem cells. PeerJ. 2017;5:e3233.
Ma Y, Pan X, Xu P, Mi Y, Wang W, Wu X, et al. Plasma microRNA alterations between EGFR-activating mutational NSCLC patients with and without primary resistance to TKI. Oncotarget. 2017;8(51):88529–36.
Fang Z, Zhang L, Liao Q, Wang Y, Yu F, Feng M, et al. Regulation of TRIM24 by miR-511 modulates cell proliferation in gastric cancer. J Exp Clin Cancer Res. 2017;36(1):17.
Guan Y, Guan X, An H, Baihetiya A, Wang W, Shao W, et al. Epigenetic silencing of miR-137 induces resistance to bicalutamide by targeting TRIM24 in prostate cancer cells. Am J Transl Res. 2019;11(5):3226–37.
Ma L, Yuan L, An J, Barton MC, Zhang Q, Liu Z. Histone H3 lysine 23 acetylation is associated with oncogene TRIM24 expression and a poor prognosis in breast cancer. Tumour Biol. 2016;37(11):14803–12.
Li YF, Zhang J, Yu L. Circular RNAs regulate Cancer onset and progression via Wnt/beta-catenin signaling pathway. Yonsei Med J. 2019;60(12):1117–28.
Jiang J, Yang X, He X, Ma W, Wang J, Zhou Q, et al. MicroRNA-449b-5p suppresses the growth and invasion of breast cancer cells via inhibiting CREPT-mediated Wnt/beta-catenin signaling. Chem Biol Interact. 2019;302:74–82.
Wang F, Zhu W, Yang R, Xie W, Wang D. LncRNA ZEB2-AS1 contributes to the tumorigenesis of gastric cancer via activating the Wnt/beta-catenin pathway. Mol Cell Biochem. 2019;456(1–2):73–83.
Masetti R, Bertuccio SN, Pession A, Locatelli F. CBFA2T3-GLIS2-positive acute myeloid leukaemia. A peculiar paediatric entity. Br J Haematol. 2019;184(3):337–47.
Sun H, Gu J, Li Z, Liu Q, Lin J, Tian Y, et al. Low expression of GLIS2 gene might associate with Radiosensitivity of gastric Cancer. J Oncol. 2019;2019:2934925.
This study was supported by the National Natural Science Foundation of China (81873186).
Ethics approval and consent to participate
The study gained the approval of the Clinical Ethical Committee of the Zhejiang Provincial Hospital of Traditional Chinese Medicine. All procedures were conducted as per the Declaration of Helsinki. All eligible participants signed the informed consent. Animal studies were conducted as per the principles and procedures ratified by the Committee on the Ethics of Animal Experiments of the Zhejiang Provincial Hospital of Traditional Chinese Medicine. Great attempts were made to minimize the number and pain of animals.
Consent for publication
The authors declared no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this article was revised: the authors would like to correct the Funding section from ‘This study was supported by the National Natural Science Foundation of China (81873186).’ to ‘None.’
About this article
Cite this article
Huang, L., Li, X., Ye, H. et al. Long non-coding RNA NCK1-AS1 promotes the tumorigenesis of glioma through sponging microRNA-138-2-3p and activating the TRIM24/Wnt/β-catenin axis. J Exp Clin Cancer Res 39, 63 (2020). https://doi.org/10.1186/s13046-020-01567-1