Novel LncRNA OXCT1-AS1 indicates poor prognosis and contributes to tumorigenesis by regulating miR-195/CDC25A axis in glioblastoma
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 123 (2021)
Long noncoding RNAs (lncRNAs) contribute to multiple biological processes in human glioblastoma (GBM). However, identifying a specific lncRNA target remains a challenge. In this study, bioinformatics methods and competing endogenous RNA (ceRNA) network regulatory rules were used to identify GBM-related lncRNAs and revealed that OXCT1 antisense RNA 1 (OXCT1-AS1) is a potential therapeutic target for the treatment of glioma.
Based on the Gene Expression Omnibus (GEO) dataset, we identified differential lncRNAs, microRNAs and mRNAs and constructed an lncRNA-associated ceRNA network.
The novel lncRNA OXCT1-AS1 was proposed to function as a ceRNA, and its potential target miRNAs were predicted through the database LncBase Predicted v.2. The expression patterns of OXCT1-AS1 in glioma and normal tissue samples were measured. The effect of OXCT1-AS1 on glioma cells was checked using the Cell Counting Kit 8 assay, cell colony formation assay, Transwell assay and flow cytometry in vitro. The dual-luciferase activity assay was performed to investigate the potential mechanism of the ceRNA network. Finally, orthotopic mouse models of glioma were created to evaluate the influence of OXCT1-AS1 on tumour growth in vivo.
In this study, it was found that the expression of lncRNA OXCT1-AS1 was upregulated in both The Cancer Genome Atlas (TCGA) GBM patients and GBM tissue samples, and high expression of OXCT1-AS1 predicted a poor prognosis. Suppressing OXCT1-AS1 expression significantly decreased GBM cell proliferation and inhibited cell migration and invasion. We further investigated the potential mechanism and found that OXCT1-AS1 may act as a ceRNA of miR-195 to enhance CDC25A expression and promote glioma cell progression. Finally, knocking down OXCT1-AS1 notably attenuated the severity of glioma in vivo.
OXCT1-AS1 inhibits glioma progression by regulating the miR-195-5p/CDC25A axis and is a specific tumour marker and a novel potential therapeutic target for glioma treatment.
Glioma is an aggressive subtype of primary brain tumours with an extremely poor prognosis; GBM is the most severe type of glioma and accounts for approximately 1/5–1/4 of primary intracranial malignancies [1, 2]. Patients with GBM have a poor prognosis, and the 5-year survival rate is less than 10% [3, 4]. The prognosis of patients diagnosed with malignant gliomas remains dismal, although treatments such as radical surgery, radiotherapy, and chemotherapy are valuable in managing these tumours . Additionally, prognostic biomarkers and therapeutic targets for gliomas have not been fully characterised [6, 7]. Therefore, it is necessary to identify novel biomarkers in glioma and reveal the molecular mechanism underlying glioma progression, to improve the early diagnosis and effective treatment of the disease.
In the eukaryotic cell genome, more than 90% of human transcripts have limited protein-coding capacity but encode noncoding RNAs, including microRNAs (miRNAs), lncRNAs, and circular RNAs (circRNAs) [8,9,10]. LncRNAs are noncoding RNA molecules of more than 200 nucleotides in length that are transcribed by RNA polymerase II and exert several regulatory functions at both the transcriptional and post-transcriptional levels . Although lncRNAs do not encode proteins, they regulate gene expression in various ways, such as genome modification, transcriptional activation, transcriptional interference and chromosome sedimentation [12, 13]. Accumulating evidence has shown that the abnormal expression of lncRNAs is closely related to the pathogenesis, progression and prognosis of malignant tumours, including GBM . Mechanistically, the ceRNA hypothesis suggests that lncRNAs function as competing endogenous RNAs in multiple human malignancies to regulate the miRNA-mRNA axis . Increasing data suggest that lncRNA/miRNA/target gene axis play important roles in GBM. For example, LncRNA HOXA-AS3/miR-455-5p/USP3 axis promotes the malignancy of glioblastoma  and Linc00152/miR-103a-3p/FEZF1 axis promotes the malignant progression of glioma stem cells .
The present study aimed to better understand the pathological process of GBM at the genome level and identify new and specific lncRNA targets. We constructed an lncRNA-related ceRNA network in GBM and identified lncRNA OXCT1-AS1 as a potential specific prognostic biomarker and therapeutic target in GBM. TCGA and quantitative real-time PCR (qRT-PCR) assays revealed that the expression of OXCT1-AS1 was significantly increased in GBM tissue samples and cell lines and that patients with higher OXCT1-AS1 expression had a worse prognosis. Additionally, we demonstrated that OXCT1-AS1 promotes GBM proliferation by regulating the miR-195/CDC25A axis in GBM. Our study is the first to report the expression pattern, biological function and potential regulatory mechanism of lncRNA OXCT1-AS1 in GBM and may provide a novel diagnostic biomarker and therapeutic target for GBM.
Materials and methods
Series matrix files of the GSE4290 dataset containing the mRNA microarray data of tissues samples from 23 epileptic and 81 GBM patients , the GSE90603 dataset containing the miRNA microarray data of 16 fresh-frozen GBM multiforme and 7 healthy brain tissue samples, and the GSE104267 dataset containing the mRNA microarray data of 9 tumour and 3 healthy tissue samples were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/) . All 26 glioma specimens and paired non-tumour tissues used in this study were collected from glioma patients who had undergone surgery. All the GBM samples were immediately frozen in liquid nitrogen until RNA was extracted.
Patient tissue preparation
A total of 26 pairs of GBM samples and their adjacent normal tissues were obtained from the First Affiliated Hospital of Harbin Medical University between 2018 and 2020. For qRT-PCR and western blot analysis, tissues were immediately frozen in liquid nitrogen. All experimental ethics were approved by the Institutional Animal Care and Use Committee at Harbin Medical University (No. HMUIRB-2008-06).
ceRNA network construction
The differentially expressed mRNAs (DEmRNAs), differentially expressed lncRNAs (DElncRNAs), and differentially expressed miRNAs (DEmiRNAs) in the GSE4290, GSE90603 and GSE104267 datasets were identified using the limma package in R . Differential expression analyses were performed using the limma settings of an FDR threshold of 0.05 and an FC lower threshold of |log(FC)| ≥ 1. We used online tools (http://bioinfogp.cnb.csic.es/tools/venny/) to conduct integrated bioinformatics analyses [21, 22]. The miRcode database (http://www.mircode.org) was used to define the relationships between the DElncRNAs and DEmiRNAs. Target genes were predicted by miRTarBase, TargetScan, and miRDB. The DElncRNA-DEmiRNA-DEmRNA network was constructed using Cytoscape (version 3.6.1) .
Functional annotation and pathway analyses
We conducted Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway analyses for the mRNAs involved in the constructed lncRNA-related network using the R package “clusterprofiler”, with a set of cut-off criteria at P < 0.05 . The online database Retrieval of Interacting Gene (STRING) was used to construct a protein-protein interaction (PPI) network . DEmRNAs were incorporated into the PPI network when they had a combined score greater than 0.4 .
Cell lines and cell culture
All human GBM cell lines (A172, LN229, U87, U251, and U373) and a normal human astrocyte cell line (NHA) were obtained from the China Infrastructure of Cell Line Resource. The cells were cultured and preserved in DMEM (GIBCO-BRL, 11965092) supplemented with 10% FBS, 100 U/mL penicillin and 100 mg/mL streptomycin in humidified air at 37 °C with 5% CO2.
OXCT1-AS1 inhibitor shRNA and an empty vector were commercially synthesised by General Biology (Anhui, China). The plasmids were transfected at 2.5 μg/well in a 6-well plate. All the transfections were performed using Lipofectamine 2000 (Invitrogen, Carlsbad, USA, 11668019). The transfection efficiency was assessed by qRT-PCR analysis. Subsequent experiments were performed at 48 h post transfection.
Total RNA was extracted from glioma tissues and cell lines using TRIzol reagent (Ambion Life Technologies, USA, A33254). Reverse transcription was performed using a First Strand cDNA Synthesis Kit (TOYOBO Life Science, Shanghai, China, FSQ-101). Next, qRT-PCR analyses were performed using the Universal SYBR-Green Master Mix (Roche, Germany, 04707516001). GAPDH was used as an endogenous control for lncRNAs and mRNAs. U6 was used as endogenous controls for miRNAs. The 2− ▵▵Ct method was used to analyse the results.
Cell proliferation assay
The transfected cells were cultured in 96-well plates at 1 × 104 per well and incubated for 24, 48, and 72 h. Cell proliferation was assayed using the Cell Counting Kit 8 assay (CCK-8; MedChem Express, China, HY-K0301) according to the manufacturer’s protocol.
Colony formation assay
The cells were plated at a density of 5 × 103 cells/well in 96-well plates. After 24 h of transfection with sh-OXCT1-AS1 and the empty vector, the cells were seeded into six-well plates (500 cells per well) and cultured for 2 weeks. Next, 0.1% crystal violet was used to stain clones, and cells were photographed using a ChemiDoc™ MP system (Bio-Rad, USA). The number of colonies was counted using ImageJ.
After transfection, the cells were fixed with 4% paraformaldehyde, permeabilised with Triton X-100 and blocked with 5% bovine serum albumin (BSA; BOSTER, USA, AR0004). Next, the cells were incubated with the Ki-67 primary antibody overnight at 4 °C, followed by Alexa Fluor 594 AffiniPure goat anti-rabbit IgG secondary antibody (ZSGB-BIO, China, ZF-0516) and DAPI (Beyotime, China, C1005). The results were determined by fluorescence microscopy the next day.
Flow cytometry analysis of the cell cycle
After transfection, the cells were harvested and stained using the Cycletest™ PLUS DNA Reagent Kit (BD Biosciences, 340242). Next, the samples were analysed using an Accuri C5 flow cytometer (BD Biosciences) to determine the distribution of the cells in the G0-G1, S, and G2-M phases.
Total cell protein was lysed using RIPA buffer (Beyotime Institute of Biotechnology, Beijing, China, P0013B) containing protease and phosphatase inhibitors. Equal amounts of lysates were separated on 12.5% SDS-PAGE gels and transferred to PVDF membranes (Millipore, Billerica, MA). After blocking in 5% skim milk with TBST, the membranes were incubated with primary antibodies overnight at 4 °C. The membranes were incubated with secondary antibodies at room temperature for 1 h. The protein bands were visualised by ChemiDoc™ MP System.
Eight-micrometre Transwell chambers (Corning Company, NY) were used to conduct migration assays. After transfection, the cells were planted into the upper chambers (1 × 105 cells/well) and cultured for 24 h. The chamber was fixed with 4% paraformaldehyde and stained with 0.1% crystal violet for 20 ~ 30 min. After washing out the crystal violet, the stained cells were counted under a microscope. Next, 100 μl of Matrigel (BD Biosciences, 354230) was added to the upper chambers, and the assay was performed as previously described to detect cell invasion ability.
Luciferase reporter assay
Fragments of OXCT1-AS1 and CDC25A 3’UTR containing miR-195 binding sites were amplified and cloned into the psiCheck2 reporter vector (Promega, Shanghai, China). Next, HEK293 cells were cotransfected with reporter vector and miR-195 mimics for 48 h. Luciferase activity was measured using the Dual Luciferase Reporter Assay (Promega, E1910), normalising firefly (experimental group) luciferase to Renilla (control group) activity.
In vivo tumour formation assay
Four-week-old BALB/C nude mice (15–20 g) were obtained from the Vital River Animal Center (Beijing, China) and randomly divided into two groups. Mice anaesthetised with isoflurane were placed in a stereotaxic frame, and the skull was exposed for intracranial injections or infusions. After stable transfection with the empty vector or shRNA, LN229 cells were collected, and 5 × 106 resuspended cells were used for each intracranial injection. Magnetic resonance imaging was performed to evaluate the intracranial lesions every 15 days after the vaccinations. The tumour volume was assessed and calculated (volume = 4/3π(a*b*c), a,b,c are the 1/2 maximum diameter at coronal, sagittal and axial position) every 3 days. The survival time of the mice was recorded, and Kaplan-Meier survival curves were plotted for each group.
Formalin-fixed tumour tissues were embedded in paraffin and sliced into 5-μm-thick sections. After blocking with 5% BSA, the sections were incubated with the CDC25A primary antibody (1:500; Affinity Biosciences, AF6252) at 4 °C overnight and secondary antibodies at 37 °C for 30 min. Next, the samples were visualised using the diaminobenzidine substrate kit (Abcam, ab64238) for 10 min. After intensive washing, the samples were counterstained with haematoxylin, dehydrated and covered with a coverslip.
The data from three independent experiments were shown as means ± standard deviations (SD). Two-group statistical analysis was performed using Student’s t test (two-tailed). For multiple group statistical analysis, we used one-way ANOVA with SPSS 19.0 (IBM, USA). A P-value< 0.05 was considered statistically significant.
Construction of the ceRNA network in GBM
From the GEO database, we downloaded the microarray data of GSE4290, which included tissue samples from 23 epileptic and 81 GBM patients, and the microarray data of GSE90603, which included 16 fresh-frozen GBM multiforme and 7 healthy brain tissue samples, and the microarray data of GSE104267, including 9 GBM and 3 epileptic tissue samples. We identified 1292 upregulated mRNAs and 1709 downregulated mRNAs in GSE4290 (Fig. 1a, S1A), 170 upregulated miRNAs and 81 downregulated miRNAs (Fig. 1b, S1B) in GSE90603, and 21 upregulated lncRNAs and 69 downregulated lncRNAs (Fig. 1c, S1C) in GSE104267, with an FC lower threshold of |log (FC)| ≥ 1 and P < 0.05. We selected the interaction pairs of lncRNAs and miRNAs based on the miRcode database and interaction pairs of miRNAs and mRNAs based on the miRDB database, TargetScan and miRTarBase. We constructed an lncRNA-miRNA-mRNA competing network. Finally, 12 DElncRNAs, 3 DEmiRNAs, and 116 DEmRNAs were employed to construct the ceRNA network using Cytoscape (Fig. 1d). We found that lncRNA OXCT1-AS1 has a connection with two of the three DEmiRNAs; therefore, we focused on OXCT1-AS1 in the subsequent experiments. According to the ceRNA network hypothesis, we found that OXCT1-AS1 expression was only inversely related to the miR-195 level. Hence, we hypothesised that miR-195 might be the downstream target of OXCT1-AS1 in the ceRNA network.
GO functional annotation and KEGG pathway enrichment analyses
We Used the ClusterProfile package in R software to perform the GO functional enrichment and the KEGG pathway enrichment analyses. GO analysis showed that the DEmRNAs we obtained were associated with “cell cycle,” “membrane enclosed lumen,” “transcription factor activity,” and so on (Fig. 2a-d, Table S1, S2 and S3). KEGG pathway analysis showed that the DEmRNAs were enriched in cancer associated pathways, such as MicroRNA in cancer, cell cycle and so on (Fig. 2f, Table 1). Both GO and KEGG analysis all enriched in cell cycle process. From the ceRNA network, we searched the cell cycle related DEmRNAs can be regulated by miR-195 and found cell division cycle 25A (CDC25A) was the only candidate. We also used the STRING database to explore the hub genes from the DEmRNAs and constructed the PPI network (Fig. 2g). We identified the top 10 DEmRNAs as the hub genes based on their linkage degree, including CDC25A (Fig. 2f). Thus, CDC25A might be the component of OXCT1-AS1/miR-195 axis regulated ceRNA network.
LncRNA OXCT1-AS1 is highly expressed in GBM and associated with prognosis
We first examined the expression pattern of lncRNA OXCT1-AS1 in various common solid cancers based TCGA database analysis and found that OXCT1-AS1 was aberrantly expressed in GBM and low-grade glioma (LGG) (Fig. 3a). In GBM, patients with high OXCT1-AS1 expression had a lower median survival (37.75 months) than those with low OXCT1-AS1 expression (53.52 months). The median survival status significantly differed between the high OXCT1-AS1 and low OXCT1-AS1 expression groups (HR = 1.52; 95% CI = 1.094–2.163; P value = 0.0132; Fig. 3b). We also analysed 26 pairs of GBM samples and their adjacent normal tissues, and the expression of OXCT1-AS1 was remarkably increased in the GBM samples (Fig. 3c). Kaplan-Meier analysis was used to evaluate OXCT1-AS1-related patient survival and revealed that higher OXCT1-AS1 expression was associated with poor survival (HR = 1.922; 95% CI = 1.154–3.203; P value = 0.0245; Fig. 3e). We wondered whether OXCT1-AS1 could be an indicator for the diagnosis of GBM. We used the 26 paired adjacent normal tissues as controls to build the ROC curve (Fig. 3d). The sensitivity and specificity were 0.769 and 0.808, respectively, and the cut-off value was 1.1545. The area under the curve was 0.817 (95% CI = 0.692–0.941; P < 0.000), and the Youden index was 0.577. In the univariable and multivariable Cox regression models, we found that OXCT1-AS1 expression was an indicator of GBM (HR = 0.468; 95% CI = 0.094–0.361; P value = 0.014; Table 2). qRT-PCR assays also demonstrated that CDC25A was upregulated in glioma tissues (Fig. 3f), in contrast to miR-195 expression but consistent with OXCT1-AS1 expression (Fig. 3g-i).
LncRNA OXCT1-AS1 has a higher expression level in GBM cell lines, and OXCT1-AS1 suppression arrests GBM cell proliferation
We performed qRT-PCR assays to assess the OXCT1-AS1 expression level in GBM cell lines, and the expression levels of OXCT1-AS1 were upregulated in GBM cells compared with those in NHA cells (Fig. 4a). Because A172 and LN229 have the highest expression of OXCT1-AS1 among the GBM cell lines, we selected them for this study. We transfected shRNA into A172 and LN229 cells to suppress OXCT1-AS1 expression. Forty-eight hours post-transfection, qRT-PCR analysis was performed, and the 1st and 2nd shRNAs significantly downregulated OXCT1-AS1 expression (Fig. 4b). We selected the 1st shRNA for further experiments. The CCK-8 assay results showed that the viability of A172 and LN229 cells was obviously decreased by shRNA transfection (Fig. 4c, f). Similarly, reduced OXCT1-AS1 expression impaired the colony formation capacities of GBM cells (Fig. 4d, e). Finally, we used Ki-67 staining to measure cell proliferation. The proliferative marker Ki67 was decreased after OXCT1-AS1 knockdown (Fig. 4g, h).
LncRNA OXCT1-AS1 promotes the migration and invasion of GBM cells in vitro
We also studied whether OXCT1-AS1 affected the migration and invasion of GBM cells. Cell migration and invasion potential were measured by Transwell assays. Both the migration and invasion of cells were significantly decreased after silencing OXCT1-AS1 (Fig. 5a-d). Additionally, knocking down OXCT1-AS1 increased the protein level of E-cadherin while decreasing the protein levels of N-cadherin and snail, further confirming our Transwell results (Fig. 5e and f).
Reduced OXCT1-AS1 levels induce cell cycle arrest in GBM cells
Our KEGG pathway enrichment analyses indicated that the obtained hub genes may regulate cell cycle processes. Therefore, we conducted flow cytometry assays to evaluate the impact of OXCT1-AS1 on the GBM cell cycle. After silencing OXCT1-AS1 in GBM cells, the number of cells in the G0/G1 phase was increased and the number of cells in the G2/M phase was decreased (Fig. 6a-d), confirming that OXCT1-AS1 promotes GBM malignant progression, partly by regulating the cell cycle process. Additionally, by constructing a ceRNA network, we found that OXCT1-AS1 may influence the GBM cell cycle by regulating CDC25A. Furthermore, we measured the levels of proteins associated with the cell cycle, such as CDC25A, CCNA1, CCNE1 and CDK2 (the downstream tyrosine dephosphorylating target of CDC25A). Knockdown of OXCT1-AS1 significantly decreased the protein levels of CDC25A, CCNA1, and CCNE1 (Fig. 5e) and upregulated CDK2 phosphorylation levels (Fig. 5f). All these experiments further supported our hypothesis.
OXCT1-AS1 enhances the expression level of CDC25A by competitively binding miR-195 in GBM cells
Initially, we conducted microarray analysis and constructed a ceRNA network to explore the underlying molecular mechanism of OXCT1-AS1. To examine whether OXCT1-AS1 functions as a ceRNA in GBM cells, we used the online software LncBase V2.0 (http://carolina.imis.athena-innovation.gr) to predict the binding sites of OXCT1-AS1 and miR-195 (Fig. 7a). qRT-PCR showed that knocking down OXCT1-AS1 increased the level of miR-195 in GBM cells (Fig. 7b). To confirm the competing sponging mechanism, we constructed wild-type and mutant (mut) OXCT1-AS1 luciferase reporters separately using the psiCheck2 vector. The dual luciferase reporter assay revealed that miR-195 mimics reduced the luciferase activity of the wild-type OXCT1-AS1 reporter, but no significant change was observed in the OXCT1-AS1 mutant reporter (Fig. 7d). These results demonstrated that OXCT1-AS1 directly targeted miR-195 and that the expression of the latter was inhibited. Based on our previous microarray analysis results, CDC25A was predicted to function downstream of miR-195 in the ceRNA network. The binding sequence between miR-195 and CDC25A was predicted and is shown in Fig. 7c. Similarly, the luciferase activity of wild-type CDC25A but not mutant-type CDC25A was reduced by miR-195 mimics (Fig. 7e). We also conducted western blotting, which revealed that CDC25A expression was decreased in the OXCT1-AS1 knockdown group compared with that in the vector group, whereas the inhibitory effect could be partially reversed by adding the miR-195 inhibitor AMO-195 (Fig. 7f). Functional colony formation assays also revealed that the miR-195 inhibitor partly reversed OXCT1-AS1 knockdown-induced GBM cell growth arrest (Fig. 7g). In summary, OXCT1-AS1 promotes GBM cell proliferation by competitively binding miR-195 and negatively regulating the miR-195/CDC25A axis.
Knocking down OXCT1-AS1 inhibits GBM growth in vivo
To investigate the impact of OXCT1-AS1 on tumorigenesis, an in vivo experiment was conducted. Primary GBM cell-derived mouse intracranial tumours were obviously smaller after knocking down OXCT1-AS1 compared with those derived from cells transfected with the empty vector (Fig. 8a and b). Consistently, the survival time in the OXCT1-AS1 silencing group was longer than that in the vector group (Fig. 8c). Next, qRT-PCR determined an obvious increase in miR-195 expression in the tumour tissue derived from sh-OXCT1-AS1-transfected cells (Fig. 8d). Through an immunohistochemistry assay, we found that the expression level of CDC25A was decreased in the OXCT1-AS1 blocking group compared with that in the vector group (Fig. 8e-f). These findings indicated that OXCT1-AS1 promotes GBM growth in vivo.
Because of the lack of effective therapeutic targets, patients with GBM usually have the worst survival rate [7, 27]. Thus, identifying new, suitable and effective therapeutic targets is critical . In recent studies, deregulation of lncRNAs was revealed in various malignancies and implicated potential biomarkers and therapeutic targets [29,30,31]. For example, lncRNA LINC00899 suppresses breast cancer progression by inhibiting miR-425 , lncRNA ZFAS1 regulates oesophageal squamous cell carcinoma cell malignant behaviours via the miR-124/STAT3 axis  and lncRNA MALAT1 promotes GBM proliferation and progression by targeting the miR-199a/ZHX1 axis . These studies indicated that lncRNAs play an important role in diagnosing and treating human cancers. In our present study, we identified DElncRNAs, DEmiRNAs and DEmRNAs in GBM and normal brain samples from the GEO database. We further constructed a putative ceRNA network and performed GO and KEGG pathway enrichment analyses. From the ceRNA network, we identified the core lncRNA OXCT1-AS1. In our subsequent studies, we focused on verifying that OXCT1-AS1 is a specific biomarker and potent therapeutic target for GBM. We identified that lncRNA OXCT1-AS1 was aberrantly expressed in GBM tissue and that increased OXCT1-AS1 expression was associated with poor median survival. This result was confirmed in an in vitro experiment, and blocking OXCT1-AS1 efficiently suppressed cell proliferation, induced cell cycle arrest and inhibited the migration and invasion and EMT related gene expression of GBM cells. It has been reported that ectopic miR-195 over-expression increases E-cadherin in prostate cancer and breast cancer cells [35, 36], these evidence may further support OXCT-AS1/miR-195 axis is associated with cell migration and invasion. In vivo experiments demonstrated that knocking down OXCT1-AS1 inhibited GBM growth. Collectively, our data demonstrated that OXCT1-AS1 acts as an oncogene in GBM and can become a specific diagnostic and prognostic biomarker as well as a therapeutic target in GBM.
To further investigate the potential mechanisms, we filtered the downregulated DEmiRNAs from the ceRNA network we constructed. Based on the ceRNA network hypothesis, we conducted bioinformatics analysis using LncBase V2.0 and showed that OXCT1-AS1 contained a conserved target site of miR-195. Previous studies have proven that miR-195 acts as a tumour suppressor in various cancers [37,38,39]. In our study, we explored the potential targets of miR-195. According to KEGG enrichment analysis, we screened the DEmRNAs matched with the cell cycle pathway and used TargetScan to predict conserved target sites. We found that CDC25A can be the direct target of miR-195. CDC25A is a member of the CDC25 phosphatase family and plays an important role in cell cycle regulation by activating cyclins-CDKs through the removal of inhibitory phosphates . Ectopic CDC25A expression accelerates the G1/S phase transition [41, 42]. In our study, knocking down OXCT1-AS1 reduced the CDC25A expression level and increased the CDK2 phosphorylation level. Flow cytometry assays also confirmed that OXCT1-AS1 knockdown caused cell cycle arrest. Additionally, luciferase activity assays and functional experiments demonstrated that OXCT1-AS1 promotes GBM proliferation by acting as a ceRNA and regulating the miR-195/CDC25A axis.
All in all, OXCT1-AS1, a novel potential prognostic factor and therapeutic target, is closely related to the malignant proliferation of GBM. High expression of OXCT1-AS1 is associated with poor prognosis and survival rate, especially in patients with recurrence GBM (Fig. S2). How to apply our gene to clinical treatment is also the focus of our future research. Recent advancements in CRISPR/Cas9 technologies for gene knockout and point mutations may facilitate the development of OXCT1-AS1 targeted cancer therapy.
In the present study, we first identified lncRNA OXCT1-AS1 as a potential predictor and therapeutic target of GBM. Knocking down OXCT1-AS1 significantly inhibited proliferation, migration, invasion and tumorigenesis in GBM. We also investigated the molecular mechanism involved and demonstrated that OXCT1-AS1 might act as a ceRNA and regulate the miR-195/CDC25A axis. Our findings reveal a novel regulatory network in GBM that can help gain insight into the pathogenesis of GBM and improve the treatment of GBM patients.
Availability of data and materials
Because of our internal policy, raw data cannot be shared.
Long noncoding RNAs
Competing endogenous RNA
OXCT1 antisense RNA 1
Gene Expression Omnibus
The Cancer Genome Atlas
Quantitative real-time PCR
Differentially expressed mRNAs
Differentially expressed miRNAs
Differentially expressed lncRNAs
Kyoto Encyclopedia of Genes and Genome
Search Tool for the Retrieval of Interacting Genes
Cell division cycle 25A
Mary D. Glioblastoma: overview of disease and treatment. Clin J Oncol Nurs. 2016;20(5 Suppl):S2–8.
Brandner S, Jaunmuktane Z. Neurological update: gliomas and other primary brain tumours in adults. J Neurol. 2018;265(3):717–27. https://doi.org/10.1007/s00415-017-8652-3.
Chen R, Smith-Cohn M, Cohen AL, Colman H. Glioma subclassifications and their clinical significance. Neurotherapeutics. 2017;14(2):1–14.
Alan O, Telli TA, Tuylu TB, Arikan R, Demircan NC, Ercelep O, et al. Prognostic factors in progressive high-grade glial tumors treated with systemic approach: A single center experience. J Oncol Pharm Pract. 2020;0(0):1078155220920684.
Jiang HH, Lin S. Current status and prospect in the treatment of glioblastoma. Zhonghua wai ke za zhi Chinese J Surg. 2020;58(1):70–4. https://doi.org/10.3760/cma.j.issn.0529-5815.2020.01.015.
Chen RS-CM, Cohen AL, Colman H. Molecular profiling of gliomas: potential therapeutic implications. Expert Rev Anticancer Ther. 2015;15(8):955–62.
Zhou L, Tang H, Wang F, Chen L, Ou S, Wu T, et al. Bioinformatics analyses of significant genes, related pathways and candidate prognostic biomarkers in glioblastoma. Mol Med Rep. 2018;18(5):4185–96. https://doi.org/10.3892/mmr.2018.9411.
Dragomir MP, Kopetz S, Ajani JA, Calin GA. Non-coding RNAs in GI cancers: from cancer hallmarks to clinical utility. Gut. 2020;69(4):748–63. https://doi.org/10.1136/gutjnl-2019-318279.
Grillone K, Riillo C, Scionti F, Rocca R, Tassone P. Non-coding RNAs in cancer: platforms and strategies for investigating the genomic “dark matter”. J Exp Clin Cancer Res. 2020;39(1):117. https://doi.org/10.1186/s13046-020-01622-x.
Liao Y, Jung SH, Kim T. A-to-I RNA editing as a tuner of noncoding RNAs in cancer. Cancer Lett. 2020;494:88–93. https://doi.org/10.1016/j.canlet.2020.08.004.
Costa FF. Non-coding RNAs: meet thy masters. Bioessays. 2010;32(7):599–608. https://doi.org/10.1002/bies.200900112.
Choudhari R, Sedano MJ, Harrison AL, Subramani R, Lin KY, Ramos EI, et al. Long noncoding RNAs in cancer: from discovery to therapeutic targets. Adv Clin Chem. 2020;95:105–47.
Zhang X-Z, Liu H, Chen S-R. Mechanisms of Long non-coding RNAs in cancers and their dynamic regulations. Cancers. 2020;12(5):1245. https://doi.org/10.3390/cancers12051245.
Peng Z, Liu C, Wu M. New insights into long noncoding RNAs and their roles in glioma. Mol Cancer. 2018;17(1):61. https://doi.org/10.1186/s12943-018-0812-2.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505(7483):344–52. https://doi.org/10.1038/nature12986.
Chen W, Li Q, Zhang G, Wang H, Zhu Z, Chen L. LncRNA HOXA-AS3 promotes the malignancy of glioblastoma through regulating miR-455-5p/USP3 axis. J Cell Mol Med. 2020;24(20):11755–67. https://doi.org/10.1111/jcmm.15788.
Yu M, Xue Y, Zheng J, Liu X, Yu H, Liu L, et al. Linc00152 promotes malignant progression of glioma stem cells by regulating miR-103a-3p/FEZF1/CDC25A pathway. Mol Cancer. 2017;16(1):110. https://doi.org/10.1186/s12943-017-0677-9.
Sun L, Hui AM, Su Q, Vortmeyer A, Kotliarov Y, Pastorino S, et al. Neuronal and glioma-derived stem cell factor induces angiogenesis within the brain. Cancer Cell. 2006;9(4):287–300. https://doi.org/10.1016/j.ccr.2006.03.003.
TE Gulluoglu S, Sahin M, Kuskucu A, Kaan Yaltirik C, Ture U, Kucukkaraduman B, et al. Simultaneous miRNA and mRNA Transcriptome profiling of Glioblastoma samples reveals a novel set of OncomiR candidates and their target genes. Brain Res. 2018;1700:199–210.
Feng F, Wu J, Gao Z, Yu S, Cui Y. Screening the key microRNAs and transcription factors in prostate cancer based on microRNA functional synergistic relationships. Medicine (Baltimore). 2017;96(1):e5679. https://doi.org/10.1097/MD.0000000000005679.
Zhang J, Zhou YJ, Yu ZH, Chen AX, Cao XC. Identification of core genes and clinical roles in pregnancy-associated breast cancer based on integrated analysis of different microarray profile datasets. Bioence Reports. 2019;39(6):BSR20190019.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106. https://doi.org/10.1186/gb-2010-11-10-r106.
Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, et al. Integration of biological networks and gene expression data using Cytoscape. Nat Protoc. 2007;2(10):2366–82. https://doi.org/10.1038/nprot.2007.324.
Yin J, Zeng X, Ai Z, Yu M, Yo W, Li S. Construction and analysis of the lncRNA-miRNA-mRNA network based on competitive endogenous RNA reveal functional lncRNAs in oral cancer. BMC Med Genet. 2020;13(1):84.
von Mering CHM, Jaeggi D, Schmidt S, Bork P, Snel B. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003;31(1):258–61. https://doi.org/10.1093/nar/gkg034.
Yi F, Sun M, Dai G, Ramani K. The intrinsic geometric structure of protein-protein interaction networks for protein interaction prediction. IEEE/ACM Trans Comput Biol Bioinform. 2016;13(1):76–85.
Ana M. María, Blanco-Prieto, João, Sousa, et al. breaching barriers in glioblastoma. Part I: molecular pathways and novel treatment approaches. Int J Pharm. 2017;531(1):372–88.
Taylor OG, Brzozowski JS, Skelding KA. Glioblastoma Multiforme: an overview of emerging therapeutic targets. Front Oncol. 2019;9:963. https://doi.org/10.3389/fonc.2019.00963.
Athina K-S, Kleita M, Andreas S, et al. Long Noncoding RNAs in Digestive System Malignancies: A Novel Class of Cancer Biomarkers and Therapeutic Targets? Gastroenterol Res Pract. 2015;2015:319861.
Fatima R, Akhade VS, Pal D, Rao SM. Long noncoding RNAs in development and cancer: potential biomarkers and therapeutic targets. Mol Cell Ther. 2015;3(1):5. https://doi.org/10.1186/s40591-015-0042-6.
Wang L, Cho KB, Li Y, Tao G, Guo B. Long noncoding RNA (lncRNA)-mediated competing endogenous RNA networks provide novel potential biomarkers and therapeutic targets for colorectal Cancer. Int J Mol Ences. 2019;20(22):5758.
Cheng B, Rong A, Zhou Q, Li W. LncRNA LINC00662 promotes colon cancer tumor growth and metastasis by competitively binding with miR-340-5p to regulate CLDN8/IL22 co-expression and activating ERK signaling pathway. J Exp Clin Cancer Res. 2020;39(1):5. https://doi.org/10.1186/s13046-019-1510-7.
Zhang H, Zhu C, He Z, Chen S, Li L, Sun C. LncRNA PSMB8-AS1 contributes to pancreatic cancer progression via modulating miR-382-3p/STAT1/PD-L1 axis. J Exp Clin Cancer Res. 2020;39(1):179. https://doi.org/10.1186/s13046-020-01687-8.
Liao K, Lin Y, Gao W, Xiao Z, Guo L. Blocking lncRNA MALAT1/miR-199a/ZHX1 axis inhibits glioblastoma proliferation and progression. Mol Ther Nucleic Acids. 2019;18:388–99. https://doi.org/10.1016/j.omtn.2019.09.005.
Jin M, Wang L, Zheng T, Yu J, Sheng R, Zhu H. MiR-195-3p inhibits cell proliferation in cervical cancer by targeting BCDIN3D. J Reprod Immunol. 2020;143:103211.
Liu X, Zhou Y, Ning YE, Gu H, Tong Y, Wang N. MiR-195-5p inhibits malignant progression of cervical Cancer by targeting YAP1. Oncotargets Ther. 2020;13:931–44. https://doi.org/10.2147/OTT.S227826.
Long Z, Wang Y. miR-195-5p Suppresses Lung Cancer Cell Proliferation, Migration, and Invasion Via FOXK1. Technol Cancer Res Treatment. 2020;19:1533033820922587.
Cai C, Chen QB, Han ZD, Zhang YQ, He HC, Chen JH, et al. miR-195 inhibits tumor progression by targeting RPS6KB1 in human prostate Cancer. Clin Cancer Res. 2015;21(21):4922–34. https://doi.org/10.1158/1078-0432.CCR-15-0217.
Singh R, Yadav V, Kumar S, Saini N. MicroRNA-195 inhibits proliferation, invasion and metastasis in breast cancer cells by targeting FASN, HMGCR, ACACA and CYP27B1. Sci Rep. 2015;5(1):17454. https://doi.org/10.1038/srep17454.
Sadeghi H, Golalipour M, Yamchi A, Farazmandfar T, Shahbazi M. CDC25A pathway toward tumorigenesis: molecular targets of CDC25A in cell-cycle regulation. J Cell Biochem. 2019;120(3):2919–28. https://doi.org/10.1002/jcb.26838.
Blomberg I, Hoffmann I. Ectopic expression of Cdc25A accelerates the G1/S transition and leads to premature activation of Cyclin E- and Cyclin A-dependent kinases. Mol Cell Biol. 1999;19(9):6183–94. https://doi.org/10.1128/MCB.19.9.6183.
Kabakci Z, Käppeli S, Cantù C, Jensen LD, König C, Toggweiler J, et al. Pharmacophore-guided discovery of CDC25 inhibitors causing cell cycle arrest and tumor regression. Sci Rep. 2019;9(1):1335. https://doi.org/10.1038/s41598-019-38579-7.
We thank Yan Tao for excellent technical assistance and thank Xiaoxiong Wang for professional comments.
This work was supported by the China Postdoctoral Science Foundation (grant no. 2020 M681888), National Key R&D Program of China (grant no. 2018YFC1312600; 2018YFC1312603), Key Research and Development Project of Zhejiang Province (grant no. 2018C03011) and TCM Key Discipline of Zhejiang Province (grant no. 2017-XK-A39).
Ethics approval and consent to participate
All experimental ethics and animal experiments conformed to the European Parliament Directive (2010/63/EU) and were approved by the Institutional Animal Care and Use Committee at Harbin Medical University (No. HMUIRB-2008-06).
Consent for publication
The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Heatmaps of DEmRNAs in the GSE4290 dataset (A), DEmiRNAs in the GSE90603 dataset (B) and DElncRNAs in the GSE104267 dataset.
GO term enrichment for biological processes of altered genes in the ceRNA network.
GO term enrichment for cellular components of altered genes in the ceRNA network.
GO term enrichment for molecular functions of altered genes in the ceRNA network.
Antibodies used in this study.
qRT-PCR primer sequences used in this study.
OXCT1-AS1 is significantly associated with poor prognosis in patients with recurrent GBM. (A) OXCT1-AS1 expression patterns in patients with recurrent GBM and the survival time and status from the TCGA database. (B) Kaplan-Meier curve analysis of OXCT1-AS1 based on TCGA recurrent GBM samples. (C) ROC curve of OXCT1-AS1 and poor prognosis in patients with recurrent GBM.
About this article
Cite this article
Zhong, C., Yu, Q., Peng, Y. et al. Novel LncRNA OXCT1-AS1 indicates poor prognosis and contributes to tumorigenesis by regulating miR-195/CDC25A axis in glioblastoma. J Exp Clin Cancer Res 40, 123 (2021). https://doi.org/10.1186/s13046-021-01928-4