- 1Department of Neurosurgery, Henan Provincial People's Hospital, Zhengzhou, China
- 2Department of Intensive Care Unit, The Third People's Hospital of Zhengzhou, Zhengzhou, China
Background: Isocitrate dehydrogenase (IDH) mutant is one of the most robust and important genetic aberrations in glioma. However, the underlying regulation mechanism of long non-coding RNA (lncRNA) in IDH mutant glioma has not been systematically portrayed.
Methods: In this work, 775 IDH mutant glioma samples with transcriptome data, including 167 samples from the Chinese Glioma Genome Atlas (CGGA) RNAseq dataset, 390 samples from The Cancer Genome Atlas (TCGA) dataset, 79 samples from GSE16011 dataset, and 139 samples from CGGA microarray dataset, were enrolled. R language and GraphPad Prism software were applied for the statistical analysis and graphical work.
Results: By comparing the differentially lncRNA genes between IDH mutant and IDH wild-type glioma samples, a four-lncRNA (JAG1, PVT1, H19, and HAR1A) signature was identified in IDH mutant glioma patients. The signature model was established based on the expression level and the regression coefficient of the four lncRNA genes. IDH mutant glioma samples could be successfully stratified into low-risk and high-risk groups in CGGA RNAseq, TCGA, GSE16011, and CGGA microarray databases. Meanwhile, multivariate Cox analysis showed that the four-lncRNA signature was an independent prognostic biomarker after adjusting for other clinicopathologic factors. Moreover, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the immune response and cellular metabolism were significantly associated with the four-lncRNA risk signature.
Conclusion: Taken together, the four-lncRNA risk signature was identified as a novel prognostic marker for IDH mutant glioma patients and may potentially lead to improvements in the lives of glioma patients.
The new established four-lncRNA risk signature was a prognostic factor for patients with IDH mutant glioma and may provide therapeutic targets in future.
Introduction
Long non-coding RNAs (lncRNAs), which are initially regarded as spurious transcriptional noise, are recently reported to play roles in cellular development, cell metabolism, and posttranscriptional modification (1). The biological functions and specific expression patterns of lncRNA have given rise to novel viewpoints on a variety of diseases, including cancer (2). Many reports were proposed to investigate the molecular mechanism of lncRNA and its impact on cancer development. The lncRNA HOX Transcript Antisense RNA (HOTAIR), which regulates the transcriptional silencing of EZH2 and the epigenetic modification of H3K27, is reported as a negative prognostic factor in breast cancer and colon cancer (3, 4).
Isocitrate dehydrogenase (IDH) gene, encoded IDH enzyme, could catalyze the oxidative decarboxylation of isocitrate to α-ketoglutarate during the production of nicotinamide adenine dinucleotide phosphate (NADPH) (5). IDH mutation, one of the earliest detectable genetic alterations, has been identified as a favorable prognostic marker of glioma patients (6). Several clinical researches reported that IDH mutation was involved in the transcription and posttranslational modification (7, 8). However, the aberrant expressions of lncRNA in IDH mutant glioma remain unknown.
In this study, we investigated the change of lncRNA transcriptional profile and lncRNA-related biological function in IDH mutant gliomas. Moreover, a four-lncRNA risk signature was established as an independent prognostic factor for IDH mutant patients accurately and robustly. The findings in our study may help in understanding the features of the lncRNA in IDH-mutant glioma and potentially contribute to the clinical management of IDH mutant glioma patients.
Materials and Methods
Samples
The RNA sequence and microarray data and corresponding clinical data were downloaded from the Chinese Glioma Genome Atlas (CGGA) database (http://www.cgga.org.cn), The Cancer Genome Atlas (TCGA) database (http://cancergenome.nih.gov/), and the GSE16011 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=gse16011). Four independent datasets were collected in our study, including CGGA RNAseq dataset (n = 167), TCGA RNAseq dataset (n = 390), GSE16011 dataset (n = 79), and CGGA microarray dataset (n = 139). Overall survival (OS) was calculated from the date of diagnosis until death or the end of follow-up. The clinicopathologic characteristics of the patients from three datasets are shown in Table 1.
Gene Ontology Analysis
Gene Ontology (GO) enrichment analysis, annotation, and visualization were performed using Metascape (http://metascape.org/) resource (9). A probability value of <0.05 was considered statistically significant.
Statistical Analysis
Differentially expressed lncRNAs (DELs) between patients in IDH mutant group and IDH wild-type group were calculated with Limma R package. Associations between continuous variables were tested using Pearson correlation analysis. The survival distribution was estimated with Kaplan–Meier survival curves. Univariate Cox and multivariate Cox regression analyses were performed to identify independent prognostic factors. The risk-score formula for predicting survival was developed based on the expression level of lncRNA genes and the regression coefficient derived from the univariate Cox regression analysis. The four-lncRNA risk score for each patient was calculated as follows:
Patients were separated into a low-risk group and a high-risk group based on the cutoff point (median value). The Venn diagram was generated by VennDiagram R package. Heatmap was conducted by pheatmap R package. The volcano plot was plotted by ggplot2 R package. Kaplan–Meier survival curves were conducted by GraphPad Prism (version 7.0). All statistical computations were performed with the statistical software environment R (version 3.6.0) (http://www.r-project.org/) and GraphPad Prism (version 7.0). p < 0.05 was considered statistically significant.
Results
Whole Genome Long Non-coding RNA Profiling in Isocitrate Dehydrogenase Mutant Glioma
We downloaded the list of 56,946 lncRNA genes from LNCipedia database (https://lncipedia.org/) and identified the DELs between IDH1 mutant and IDH1 wild-type samples. The lncRNAs with a minimum of 2-fold difference in expression level and false discovery rate (FDR) of <0.05 were considered to be statistically significant. As shown in Figure 1A, volcano plots were generated to visualize the DELs between IDH mutant and IDH wild-type glioma in CGGA RNAseq, TCGA, and GSE16011 datasets. Red and purple plots represented the downregulated and upregulated lncRNAs, respectively. In total, 215 DELs, 78 DELs, and 51 DELs were significantly upregulated in IDH mutant glioma in CGGA, TCGA, and GSE16011 datasets, respectively (Figure 1B). Meanwhile, 214 DELs, 68 DELs, and 18 DELs were significantly downregulated in IDH mutant glioma in CGGA, TCGA, and GSE16011, respectively (Figure 1C). A total of 33 lncRNAs were found to be differentially expressed in all three independent datasets, including 21 lncRNA genes that were upregulated and 12 lncRNA genes that were downregulated in IDH1 mutant glioma (Figure 1D). The list of DELs is shown in Supplementary Table 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/573264/fneur-11-573264-HTML/image_m/fneur-11-573264-g001.jpg)
Figure 1. The deferentially expressed long non-coding RNA (lncRNA) genes that are significantly expressed between isocitrate dehydrogenase (IDH) mutant and IDH wild-type glioma. Volcano plots show the differentially expressed lncRNA genes in Chinese Glioma Genome Atlas (CGGA) RNAseq, The Cancer Genome Atlas (TCGA), and GSE16011 datasets. Each point represents a different gene, and the red and purple colors of the points categorize the downregulated and upregulated genes, respectively (A). Venn diagram presents the workflow used to identify common upregulated and downregulated lncRNA genes in three datasets (B,C). Expression patterns of 33 differentially expressed lncRNA genes in IDH mutant and IDH wild-type glioma (D).
Identification of the Four-lncRNA Signature for Isocitrate Dehydrogenase Mutant Glioma
We assessed the prognostic performance of 33 DELs in IDH mutant glioma by univariate Cox regression analysis. The DELs were separated into risky lncRNAs [hazard ratio (HR) > 1] group and protective lncRNAs (HR < 1) group. With this definition, we found three risky lncRNAs (JAG1, PVT1, and H19) and one protective lncRNA (HAR1A) in all three databases. We then used four identified lncRNAs to construct a signature by risk score method. The risk score was calculated for each sample in three datasets. And then the IDH mutant patients were divided into the low-risk group and high-risk group based on the cutoff value (median four-lncRNA risk score). The four-lncRNA risk score and OS distribution are also shown in Figures 2A–C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/573264/fneur-11-573264-HTML/image_m/fneur-11-573264-g002.jpg)
Figure 2. Distribution of risk score (A), overall survival (OS) (B), and the expression profile of four long non-coding RNA (lncRNA) genes with prognostic value (C).
Validation of Four-lncRNA Signature in Four Datasets
To evaluate whether the four-lncRNA risk signature could stratify glioma patients according to prognostic outcome, IDH mutant glioma samples in CGGA RNAseq dataset were divided into low-risk and high-risk groups according to the median risk score. Survival curves revealed that patients in the low-risk group had a significantly more favorable outcome than those in the high-risk group (Figure 3A, p < 0.0001). Meanwhile, the patients in the high-risk group had survival advantage over those in IDH wild-type group (Figure 3A, p = 0.0007). Then we evaluated the prognostic value of four-lncRNA risk signature in TCGA and GSE16011 datasets. As expected, the patients in the low-risk group had statically significantly longer OS than those in the high-risk group in TCGA (Figure 3B, p = 0.0013) and GSE16011 (Figure 3C, p < 0.0001) datasets. The CGGA microarray dataset, an independent validation dataset, was used to assess the prognostic value of four-lncRNA signature further. Similarly, the survival curves indicated that OS of high-risk patients was significantly shorter compared with that in the low-risk group (Figure 3D, p = 0.027). In addition, samples in IDH wild type displayed significantly shorter OS than those in the high-risk group in TCGA and CGGA microarray datasets (Figures 3B,D, p < 0.0001). However, the survival was similar in high-risk and IDH wild-type samples in GSE16011 dataset (Figure 3C, p = 0.4269).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/573264/fneur-11-573264-HTML/image_m/fneur-11-573264-g003.jpg)
Figure 3. Based on the value of four-long non-coding RNA (lncRNA) risk signature, the isocitrate dehydrogenase (IDH) mutant glioma samples were divided into the low-risk and high-risk groups with different prognoses in Chinese Glioma Genome Atlas (CGGA) RNAseq (A), The Cancer Genome Atlas (TCGA) (B), GSE16011 (C), and CGGA microarray datasets (D). (E) Clinicopathologic features in different groups stratified by the signature in TCGA dataset. Rows represent corresponding items (DNA methylation cluster, age, KPS, grade, TCGA molecular subtype, MGMT promoter status, TERT promoter status, and ATRX status). *p < 0.05, **p < 0.01, ***p < 0.001 and ****p < 0.0001, respectively.
Then we combined the methylation data, whole genome sequencing data, and clinical data to explore the features of IDH mutant glioma patients in the low-risk group and high-risk group. As shown in Figure 3E, the sample clusters, including codel DNA methylation cluster, the WHO grade II cluster, proneural cluster, MGMT promoter methylated cluster, TERT promoter mutant cluster, and ATRX wild-type cluster, were highly enriched in the low-risk group, whereas those in the high-risk group were moderately enriched for G-CIMP-high and G-CIMP-low clusters, WHO grade III and IV clusters, classical and mesenchymal clusters, MGMT promoter unmethylated cluster, TERT promoter wild-type cluster, and ATRX mutant cluster (p < 0.05, chi-square test). To our surprise, the patients in the high-risk group were significantly younger than in the low-risk group (p < 0.001, unpaired t-test).
These findings confirmed a strong correlation between the four-lncRNA risk signature and malignancy in glioma and indicated the powerful ability of the signature in stratifying the IDH mutant glioma.
The Prognostic Value of Four-Long Non-coding RNA in Different Subtype of Isocitrate Dehydrogenase Mutant Glioma
Then we observed the prognostic value of four-lncRNA risk signature in the different subtypes of glioma. We demonstrated the prognostic value of four-lncRNA signature in both male (Figure 4A, p < 0.0001) and female (Figure 4B, p < 0.0001). In addition, similar trends were observed in younger patients (Figure 4C, p < 0.0001) or older (Figure 4D, p = 0.001) patients. Moreover, we also found that an obvious difference in survival outcome existed between low-risk patients and high-risk patients in grade II (Figure 4E, p = 0.0332), grade III (Figure 4F, p = 0.0059), and grade IV (Figure 4G, p = 0.0368), respectively. Those results indicated that the prediction power of four-lncRNA risk signature was not decreased in the subgroup analysis for gender, age, and grade.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/573264/fneur-11-573264-HTML/image_m/fneur-11-573264-g004.jpg)
Figure 4. The four-long non-coding RNA (lncRNA) risk signature retained its prognostic significance in the subgroups stratified by gender (A,B), age (C,D), and grades (E–G).
The Four-Long Non-coding RNA Signature Is an Independent Factor for Isocitrate Dehydrogenase Mutant Glioma
Furthermore, the univariate and multivariate Cox regression analyses were conducted to evaluate the prognostic independence of the four-lncRNA risk signature in four independent datasets. As shown in Figure 5, after the clinical factors were adjusted, including age, gender, WHO grades, and 1p/19q status, the four-lncRNA signature was an independent prognostic factor in IDH mutant glioma patients in CGGA RNAseq dataset, TCGA dataset, GSE dataset, and CGGA microarray dataset. These results confirmed the independent prognostic value of the risk score.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/573264/fneur-11-573264-HTML/image_m/fneur-11-573264-g005.jpg)
Figure 5. (A-D) Forest plot of hazard ratios for overall survival assessed by four-long non-coding RNA (lncRNA) signature and clinicopathologic factors. The four-lncRNA signature was an independent prognostic factor after adjusting age, gender, WHO grades, and 1p19q status in four independent datasets.
Functional Analysis of Four-Long Non-coding RNA Signature
To investigate the biological functions of IDH mutant gliomas correlated with four-lncRNA risk signature, we screened co-expressed genes using Pearson correlation analysis (Pearson |r| > 0.4) in CGGA RNAseq, TCGA, and GSE16011 datasets. We found that 453 genes were significantly upregulated (Figure 6A) and 268 genes were significantly downregulated (Figure 6D) in the three datasets. IDH mutant glioma samples were arranged in order of increasing risk factor (Figure 6B). We found that patients in the high-risk group have to 1p/19q intact, WHO II grade and WHO III grade, and classical and mesenchymal TCGA subtypes. The GO and KEGG analyses depicted that the upregulated genes were strongly correlated with the immune-related, cellular component-related, and metabolism-related functions, such as myeloid leukocyte activation, immune response, and extracellular matrix organization (Figure 6C). Meanwhile, the downregulated genes were more involved in normal biological progress, such as chemical synaptic transmission and neuronal system (Figure 6E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/573264/fneur-11-573264-HTML/image_m/fneur-11-573264-g006.jpg)
Figure 6. Functional and pathway analyses of four-long non-coding RNA (lncRNA) risk signature in isocitrate dehydrogenase (IDH) mutant glioma. Venn diagram shows the overlapped positive/negative genes of four-lncRNA risk signature in three datasets (A,D). The landscape of clinical and molecular features in associations with four-lncRNA risk score (B). The biological function significantly associated with positive (C) and negative genes (E).
Discussion
Even with numerous attempts to explore the new molecular and clinical markers to improve the accuracy of diagnosis and prognosis in glioma, the IDH mutant is the most robust and widely applied in clinical practices (10). IDH mutant has been found to be an inciting genetic alteration in gliomagenesis and has profound implications in tumor progression and clinical outcome (11). LncRNA genes are functionally defined as transcripts >200 nucleotides in length with no appreciable peptide products. Recent researches have reported that lncRNA exceeded a more exquisitely regulation function than mRNA in many specific organs (12). The transcriptome profile, genetic alteration characteristics, and methylation data were widely used to stratify the patients in clinical practices (13). However, the alteration of the lncRNA profile in IDH mutant glioma patients remains largely unknown.
Recently, numerous studies reported that lncRNAs emerged as regulators of human carcinogenesis by affecting the expression of key tumor suppressor genes and oncogenes. By analyzing the transcriptomics data from 775 IDH mutant glioma patients, we found that 21 lncRNAs (62.6% of 33 lncRNAs) were upregulated and 12 lncRNAs were downregulated, which suggested that increased lncRNAs play a more important role in IDH mutant glioma. The univariate Cox regression method was further performed to screen the lncRNAs, which were significantly associated with OS in IDH mutant glioma samples in three independent datasets. Finally, four lncRNAs (JAG1, PVT1, H19, and HAR1A) with the independent prognostic value were selected. Among the four lncRNAs, JAG1, PVT1, and H19 were risk-associated genes, while HAR1A was a protective gene.
JAG1 (Jagged Canonical Notch Ligand 1) is located on chromosome 20p12.2 and composed of 26 exons. Notch signaling system, which plays important roles in human development, tissue homeostasis, and disease, is a highly conserved and well-characterized pathway. JAG1 is one of five canonical ligands in Notch signaling pathway expressed by mammalian cells (14, 15). Jag1/Notch axis participated in tumorigenesis and in promoting the progression of colorectal cancer (16), ovarian cancer (17), and cancer stem cell (18). Meanwhile, multiple associations between Jag1/Notch signaling and other oncogenic factors in regulating tumor microenvironment were reported, including interleukins, TGF-β, and the Hippo, Wnt/β-catenin, and NF-κB signaling pathways (19–22). Circular RNAs (circRNAs) are covalently closed loop-like structure with no 5′ to 3′ polarity (23). Plasmacytoma variant translocation 1 (PVT1) gene is encoded by chromosome band 8q24.21 and has both circRNAs and linear lncRNA isoforms (24). The genomic amplification, rearrangements, and increased PVT1 expression exhibited oncogenic and epigenetic regulation effects on breast cancer (25), prostate cancer (26), and acute myeloid leukemia (27). Meanwhile, c-Myc and PVT1 were co-amplified in many cancers (28). In addition, PVT1 was involved in the differentiation of immune cells under physiological and pathological conditions (29, 30). H19 Imprinted Maternally Expressed Transcript (H19), a 2.3-kb lncRNA encoded by chromosome 11p15.5, is associated with the regulation of stem cell differentiation (31). The miR-675-5p, derived from H19 (32), is crucial in downregulated ubiquitin ligase E3 family and subsequently influenced the stability and function of the epidermal growth factor receptor (EGFR) and c-Met (33). Meanwhile, Chao Cui et al. found that H19 plays an important role in the promotion of hypoxia/oxygenation (h/R) injury and apoptosis by induction of autophagy via the inhibiting PI3K–Akt–mTOR pathway in hepatocellular cancer (34). Guan et al. confirmed that H19 is not only overexpressed in glioma tissue and cells, but its expression was also associated with the progression of glioma (35). The investigations performed on H19 indicated that H19 acts as an oncogene in tumors and is associated with tumorigenic phenotypes (36). Highly Accelerated Region 1A (HAR1A), a rapidly evolving cis-antisense segment on chromosome 20q13.33, is mainly expressed in the human brain (37, 38). Previous studies reported that the HAR1A expression levels were reduced in tissue specimens coinfected with hepatocellular carcinoma (HCC) and that low expression of HAR1A in patients with HCC was significantly associated with advanced histological grade and TNM stage. The HAR1A was a favorable prognostic marker for patients with HCC (37). Moreover, consistent with our findings, Zou et al. demonstrated that diffuse glioma patients with lower expression of HAR1A had an unfavorable outcome (39). HAR1A could be potentially used as a prognostic biomarker for HCC and glioma.
The lncRNAs can affect the expression of their flanking coding genes and be functionally related to neighboring coding genes (40). Our analysis indicated that JAG1, PVT1, H19, and HAR1A were aberrantly expressed in IDH mutant gliomas and were core lncRNA genes in regulating the expression of protein-encoding gene in IDH mutant gliomas.
Risk score method is widely applied in clinical studies (41). Here, taking consideration of the expression level and the prognostic contribution of lncRNAs, we established a four-lncRNA risk signature to stratify IDH mutant patients. By applying the risk signature, IDH mutant patients could be divided into the low-risk group and high-risk group. The risk signature exhibited its prognostic value in CGGA RNAseq, TCGA, GSE16011, and CGGA microarray datasets. Patients in the low-risk group significantly had a longer survival than those in the high-risk group. In addition, we found that the four-lncRNA risk signature can also identify high-risk patients in age ≤ 40 and age > 40; grade prognostic prediction power of four-lncRNA risk signature was not decreased in molecular or clinical factor subgroups of IDH mutant glioma. This result emphasized that the mechanism of glioma initiation and progression in low-risk patients and high-risk patients has great distinction. Meanwhile, the univariate and multivariate Cox regression models showed that the four-lncRNA risk signature was an independent prognostic factor.
Moreover, the GO and KEGG pathway analyses revealed that the four-lncRNA signature was positively correlated with leukocyte activation, regulation of innate immune response, and regulation of the Rho GTPase signaling. The risk score was negatively correlated with chemical synaptic transmission, neuronal system, behavior, and regulation of iron transport. The biological functional differences partially interpreted the survival differences between low-risk and high-risk group patients. The abnormal expression of lncRNA genes in the high-risk group could probably influence the outcome of IDH mutant patients by disturbing immune environment and cancer metabolism.
Targeting lncRNAs in cancer is a promising novel approach to cancer treatment. It is well-known that Notch signaling plays a key role in the occurrence, progression, and recurrence of cancer (42). As a novel candidate for blocking Notch signaling, JAG1 had unique characteristics in modifying Notch signal. MicroRNA-26b and microRNA-377-3p could suppress tumor development and invasion through sponging JAG1/Notch signaling (43, 44). The role of PVT1 as a cancer biomarker is widely accepted (45). It has been demonstrated that PVT1 was quantifiable in cancer cells and in serum from cancer patients (46, 47). Knockdown of PVT1 inhibited cell proliferation and arrested cell cycle at G1 stage via recruiting EZH2 (zeste homolog 2) and regulating p15 and p16 epigenetically (28). Moreover, some new insights on the role of PVT1 in drug response are emerging, including the elevated expression of PVT1 in the vincristine-resistant AMU-ML2 DLBCL line (48). Recently, the research performed by Li et al. confirmed that the downregulated expression of H19 increased the efficacy of temozolomide (TMZ) chemotherapy in glioblastoma multiforme (GBM) cells (49). Guan et al. indicated that the downregulation of H19 resulted in Wnt/β-catenin pathway inactivation and thereby inhibited glioma cell proliferation, invasion, and migration (35). Zou et al. found that glioma patients with high HAR1A expression could be benefit more from chemotherapy and radiotherapy (39). Therefore, the lncRNAs are potential therapeutic targets because of the various biological functions in cancers. In addition, the combination of TMZ chemotherapy and lncRNA inhibitor might show synergistic effects in increasing the chances to cure or maintain a long-term remission for glioma patients.
The four-lncRNA risk signature could potentially serve as a prognostic marker for IDH mutant glioma. However, this research has a few limitations. Our study was designed on the basis of retrospective analysis. A large prospective study with large sample size is required to assess the risk score model in future.
To sum up, we established a four-lncRNA risk signature that could predict the survival for IDH mutant glioma patients across multi-platform and multi-population. Our findings indicated that the immune environment and metabolism regulation were more chaotic in IDH mutant glioma patients with higher risk score, and patients in that group may suffer from worse outcome. The combination of lncRNAs was more objective and robust in evaluating the prognosis for IDH mutant glioma patients. The four-lncRNA risk signature could potentially serve as a prognostic score and therapeutic target for IDH mutant glioma patients in the future.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
YC: conception, analysis, and design. YG and HC: acquisition of data. FM and HC: reviewed submitted version of manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We appreciate the generosity of TCGA and CGGA network for sharing a huge amount of data.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2020.573264/full#supplementary-material
References
1. Sanchez Calle A, Kawamura Y, Yamamoto Y, Takeshita F, Ochiya T. Emerging roles of long non-coding RNA in cancer. Cancer Sci. (2018) 109:2093–100. doi: 10.1111/cas.13642
2. Schmitt AM, Chang HY. Long noncoding RNAs in cancer pathways. Cancer Cell. (2016) 29:452–63. doi: 10.1016/j.ccell.2016.03.010
3. Pawlowska E, Szczepanska J, Blasiak J. The long noncoding RNA HOTAIR in breast cancer: does autophagy play a role? Int J Mol Sci. (2017) 18:2317. doi: 10.3390/ijms18112317
4. Luo ZF, Zhao D, Li XQ, Cui YX, Ma N, Lu CX, et al. Clinical significance of HOTAIR expression in colon cancer. World J Gastroenterol. (2016) 22:5254–9. doi: 10.3748/wjg.v22.i22.5254
5. Parsons DW, Jones S, Zhang X, Lin JC, Leary RJ, Angenendt P, et al. An integrated genomic analysis of human glioblastoma multiforme. Science. (2008) 321:1807–12. doi: 10.1126/science.1164382
6. Yan H, Parsons DW, Jin G, McLendon R, Rasheed BA, Yuan W, et al. IDH1 and IDH2 mutations in gliomas. N Engl J Med. (2009) 360:765–73. doi: 10.1056/NEJMoa0808710
7. Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, et al. The 2016 World Health Organization classification of tumors of the central nervous system: a summary. Acta Neuropathol. (2016) 131:803–20. doi: 10.1007/s00401-016-1545-1
8. Tesileanu CMS, Dirven L, Wijnenga MMJ, Koekkoek JAF, Vincent A, Dubbink HJ, et al. Survival of diffuse astrocytic glioma, IDH1/2-wildtype, with molecular features of glioblastoma, WHO grade IV: a confirmation of the cIMPACT-NOW criteria. Neuro Oncol. (2019) 22:515–23. doi: 10.1093/neuonc/noz200
9. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. (2019) 10:1523. doi: 10.1038/s41467-019-09234-6
10. Zhang CB, Bao ZS, Wang HJ, Yan W, Liu YW, Li MY, et al. Correlation of IDH1/2 mutation with clinicopathologic factors and prognosis in anaplastic gliomas: a report of 203 patients from China. J Cancer Res Clin Oncol. (2014) 140:45–51. doi: 10.1007/s00432-013-1519-9
11. Yan W, Zhang W, You G, Bao Z, Wang Y, Liu Y, et al. Correlation of IDH1 mutation with clinicopathologic factors and prognosis in primary glioblastoma: a report of 118 patients from China. PLoS ONE. (2012) 7:e30339. doi: 10.1371/journal.pone.0030339
12. Fang Y, Fullwood MJ. Roles, functions, and mechanisms of long non-coding RNAs in cancer. Genomics Proteomics Bioinformatics. (2016) 14:42–54. doi: 10.1016/j.gpb.2015.09.006
13. Molinaro AM, Taylor JW, Wiencke JK, Wrensch MR. Genetic and molecular epidemiology of adult diffuse glioma. Nat Rev Neurol. (2019) 15:405–17. doi: 10.1038/s41582-019-0220-2
14. Weinmaster G, Fischer JA. Notch ligand ubiquitylation: what is it good for? Dev Cell. (2011) 21:134–44. doi: 10.1016/j.devcel.2011.06.006
15. Yamakawa T, Atsumi Y, Kubo S, Yamagishi A, Morita I, Matsuno K. Insight into notch signaling steps that involve pecanex from dominant-modifier screens in drosophila. Genetics. (2018) 209:1099–119. doi: 10.1534/genetics.118.300935
16. Dai Y, Wilson G, Huang B, Peng M, Teng G, Zhang D, et al. Silencing of Jagged1 inhibits cell growth and invasion in colorectal cancer. Cell Death Dis. (2014) 5:e1170. doi: 10.1038/cddis.2014.137
17. Li D, Masiero M, Banham AH, Harris AL. The notch ligand JAGGED1 as a target for anti-tumor therapy. Front Oncol. (2014) 4:254. doi: 10.3389/fonc.2014.00254
18. Sun J, Luo Z, Wang G, Wang Y, Wang Y, Olmedo M, et al. Notch ligand Jagged1 promotes mesenchymal stromal cell-based cartilage repair. Exp Mol Med. (2018) 50:126. doi: 10.1038/s12276-018-0151-9
19. Hu M, Yuan X, Liu Y, Tang S, Miao J, Zhou Q, et al. IL-1beta-induced NF-kappaB activation down-regulates miR-506 expression to promotes osteosarcoma cell growth through JAG1. Biomed Pharmacother. (2017) 95:1147–55. doi: 10.1016/j.biopha.2017.08.120
20. Colombo M, Galletti S, Bulfamante G, Falleni M, Tosi D, Todoerti K, et al. Multiple myeloma-derived Jagged ligands increases autocrine and paracrine interleukin-6 expression in bone marrow niche. Oncotarget. (2016) 7:56013–29. doi: 10.18632/oncotarget.10820
21. Zhou J, Jain S, Azad AK, Xu X, Yu HC, Xu Z, et al. Notch and TGFbeta form a positive regulatory loop and regulate EMT in epithelial ovarian cancer cells. Cell Signal. (2016) 28:838–49. doi: 10.1016/j.cellsig.2016.03.016
22. Tschaharganeh DF, Chen X, Latzko P, Malz M, Gaida MM, Felix K, et al. Yes-associated protein up-regulates Jagged-1 and activates the Notch pathway in human hepatocellular carcinoma. Gastroenterology. (2013) 144:1530–42.e12. doi: 10.1053/j.gastro.2013.02.009
23. Ebbesen KK, Hansen TB, Kjems J. Insights into circular RNA biology. RNA Biol. (2017) 14:1035–45. doi: 10.1080/15476286.2016.1271524
24. Shtivelman E, Bishop JM. Effects of translocations on transcription from PVT. Mol Cell Biol. (1990) 10:1835–9. doi: 10.1128/MCB.10.4.1835
25. Tang J, Li Y, Sang Y, Yu B, Lv D, Zhang W, et al. LncRNA PVT1 regulates triple-negative breast cancer through KLF5/beta-catenin signaling. Oncogene. (2018) 37:4723–34. doi: 10.1038/s41388-018-0310-4
26. Yang J, Li C, Mudd A, Gu X. LncRNA PVT1 predicts prognosis and regulates tumor growth in prostate cancer. Biosci Biotechnol Biochem. (2017) 81:2301–6. doi: 10.1080/09168451.2017.1387048
27. El-Khazragy N, Elayat W, Matbouly S, Seliman S, Sami A, Safwat G, et al. The prognostic significance of the long non-coding RNAs “CCAT1, PVT1” in t(8;21) associated acute myeloid leukemia. Gene. (2019) 707:172–7. doi: 10.1016/j.gene.2019.03.055
28. Salehi M, Sharifi M, Bagheri M. Knockdown of long noncoding RNA plasmacytoma variant translocation 1 with antisense locked nucleic acid GAPMERS exerts tumor-suppressive functions in human acute erythroleukemia cells through downregulation of C-MYC expression. Cancer Biother Radiopharm. (2019) 34:371–9. doi: 10.1089/cbr.2018.2510
29. Zhang S, Zhang G, Liu J. Long noncoding RNA PVT1 promotes cervical cancer progression through epigenetically silencing miR-200b. APMIS. (2016) 124:649–58. doi: 10.1111/apm.12555
30. Shen CJ, Cheng YM, Wang CL. LncRNA PVT1 epigenetically silences miR-195 and modulates EMT and chemoresistance in cervical cancer cells. J Drug Target. (2017) 25:637–44. doi: 10.1080/1061186X.2017.1307379
31. Zhu A, Chu L, Ma Q, Li Y. Long non-coding RNA H19 down-regulates miR-181a to facilitate endothelial angiogenic function. Artif Cells Nanomed Biotechnol. (2019) 47:2698–705. doi: 10.1080/21691401.2019.1634577
32. Dey BK, Pfeifer K, Dutta A. The H19 long noncoding RNA gives rise to microRNAs miR-675-3p and miR-675-5p to promote skeletal muscle differentiation and regeneration. Genes Dev. (2014) 28:491–501. doi: 10.1101/gad.234419.113
33. Vennin C, Spruyt N, Dahmani F, Julien S, Bertucci F, Finetti P, et al. H19 non coding RNA-derived miR-675 enhances tumorigenesis and metastasis of breast cancer cells by downregulating c-Cbl and Cbl-b. Oncotarget. (2015) 6:29209–23. doi: 10.18632/oncotarget.4976
34. Cui C, Li Z, Wu D. The long non-coding RNA H19 induces hypoxia/reoxygenation injury by up-regulating autophagy in the hepatoma carcinoma cells. Biol Res. (2019) 52:32. doi: 10.1186/s40659-019-0239-2
35. Guan N, Wang R, Guo WS, Lai YJ, Zhang YD, Cheng YY. Long non-coding RNA H19 regulates the development of gliomas through the Wnt/beta-catenin signaling pathway. Eur Rev Med Pharmacol Sci. (2019) 23:4243–53. doi: 10.26355/eurrev_201905_17929
36. Gao LM, Xu SF, Zheng Y, Wang P, Zhang L, Shi SS, et al. Long non-coding RNA H19 is responsible for the progression of lung adenocarcinoma by mediating methylation-dependent repression of CDH1 promoter. J Cell Mol Med. (2019) 23:6411–28. doi: 10.1111/jcmm.14533
37. Shi Z, Luo Y, Zhu M, Zhou Y, Zheng B, Wu D, et al. Expression analysis of long non-coding RNA HAR1A and HAR1B in HBV-induced hepatocullular carcinoma in Chinese patients. Lab Med. (2019) 50:150–7. doi: 10.1093/labmed/lmy055
38. Tolosa A, Sanjuan J, Leal C, Costas J, Molto MD, de Frutos R. Rapid evolving RNA gene HAR1A and schizophrenia. Schizophr Res. (2008) 99:370–2. doi: 10.1016/j.schres.2007.10.011
39. Zou H, Wu LX, Yang Y, Li S, Mei Y, Liu YB, et al. lncRNAs PVT1 and HAR1A are prognosis biomarkers and indicate therapy outcome for diffuse glioma patients. Oncotarget. (2017) 8:78767–80. doi: 10.18632/oncotarget.20226
40. Ponjavic J, Oliver PL, Lunter G, Ponting CP. Genomic and transcriptional co-localization of protein-coding and long non-coding RNA pairs in the developing brain. PLoS Genet. (2009) 5:e1000617. doi: 10.1371/journal.pgen.1000617
41. Zhang CB, Zhu P, Yang P, Cai JQ, Wang ZL, Li QB, et al. Identification of high risk anaplastic gliomas by a diagnostic and prognostic signature derived from mRNA expression profiling. Oncotarget. (2015) 6:36643–51. doi: 10.18632/oncotarget.5421
42. Xue TC, Zou JH, Chen RX, Cui JF, Tang ZY, Ye SL. Spatial localization of the JAG1/Notch1/osteopontin cascade modulates extrahepatic metastasis in hepatocellular carcinoma. Int J Oncol. (2014) 45:1883–90. doi: 10.3892/ijo.2014.2630
43. Wang L, Wang W, Wu Y. MicroRNA-26b acts as an antioncogene and prognostic factor in cervical cancer. Oncol Lett. (2019) 17:3418–24. doi: 10.3892/ol.2019.9965
44. Tang L, Yang B, Cao X, Li Q, Jiang L, Wang D. MicroRNA-377-3p inhibits growth and invasion through sponging JAG1 in ovarian cancer. Genes Genomics. (2019) 41:919–26. doi: 10.1007/s13258-019-00822-w
45. Onagoruwa OT, Pal G, Ochu C, Ogunwobi OO. Oncogenic role of PVT1 and therapeutic implications. Front Oncol. (2020) 10:17. doi: 10.3389/fonc.2020.00017
46. Xie Z, Chen X, Li J, Guo Y, Li H, Pan X, et al. Salivary HOTAIR and PVT1 as novel biomarkers for early pancreatic cancer. Oncotarget. (2016) 7:25408–19. doi: 10.18632/oncotarget.8323
47. Huang C, Liu S, Wang H, Zhang Z, Yang Q, Gao F. LncRNA PVT1 overexpression is a poor prognostic biomarker and regulates migration and invasion in small cell lung cancer. Am J Transl Res. (2016) 8:5025–34. Available online at: http://www.ajtr.org/files/ajtr0024362.pdf
48. Mizuno S, Hanamura I, Ota A, Karnan S, Kanasugi J, Nakamura A, et al. Establishment and characterization of a novel vincristine-resistant diffuse large B-cell lymphoma cell line containing the 8q24 homogeneously staining region. FEBS Open Bio. (2018) 8:1977–91. doi: 10.1002/2211-5463.12538
Keywords: lncRNA profile, risk signature, IDH mutant, glioma, prognosis
Citation: Chen Y, Guo Y, Chen H and Ma F (2020) Long Non-coding RNA Expression Profiling Identifies a Four-Long Non-coding RNA Prognostic Signature for Isocitrate Dehydrogenase Mutant Glioma. Front. Neurol. 11:573264. doi: 10.3389/fneur.2020.573264
Received: 17 June 2020; Accepted: 05 October 2020;
Published: 20 November 2020.
Edited by:
Leonora Balaj, Massachusetts General Hospital and Harvard Medical School, United StatesReviewed by:
Christopher James Pirozzi, Duke University, United StatesJinquan Cai, Harbin Medical University, China
Copyright © 2020 Chen, Guo, Chen and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fengjin Ma, cys0367@126.com; Hang Chen, adam7055@sina.com