- 1Department of Neurosurgery, Tianjin 4th Central Hospital, Tianjin, China
- 2Department of Psychiatry, Characteristic Medical Center of the Chinese People’s Armed Police Force, Tianjin, China
- 3Department of Urology, Tianjin Institute of Urology, The Second Hospital of Tianjin Medical University, Tianjin, China
- 4Central Laboratory, Tianjin 4th Central Hospital, Tianjin, China
- 5Department of Laboratory, Jiangsu Provincial Corps Hospital of Chinese People’s Armed Police Force, Yangzhou, China
- 6Postgraduate School, Medical School of Chinese PLA, Beijing, China
- 7Department of Thoracic Surgery, The First Medical Centre, Chinese PLA General Hospital, Beijing, China
Accumulating studies have confirmed the crucial role of long non-coding RNAs (ncRNAs) as favorable biomarkers for cancer diagnosis, therapy, and prognosis prediction. In our recent study, we established a robust model which is based on multi-gene signature to predict the therapeutic efficacy and prognosis in glioblastoma (GBM), based on Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA) databases. lncRNA-seq data of GBM from TCGA and CGGA datasets were used to identify differentially expressed genes (DEGs) compared to normal brain tissues. The DEGs were then used for survival analysis by univariate and multivariate COX regression. Then we established a risk score model, depending on the gene signature of multiple survival-associated DEGs. Subsequently, Kaplan-Meier analysis was used for estimating the prognostic and predictive role of the model. Gene set enrichment analysis (GSEA) was applied to investigate the potential pathways associated to high-risk score by the R package “cluster profile” and Wiki-pathway. And five survival associated lncRNAs of GBM were identified: LNC01545, WDR11-AS1, NDUFA6-DT, FRY-AS1, TBX5-AS1. Then the risk score model was established and shows a desirable function for predicting overall survival (OS) in the GBM patients, which means the high-risk score significantly correlated with lower OS both in TCGA and CGGA cohort. GSEA showed that the high-risk score was enriched with PI3K-Akt, VEGFA-VEGFR2, TGF-beta, Notch, T-Cell pathways. Collectively, the five-lncRNAs signature-derived risk score presented satisfactory efficacies in predicting the therapeutic efficacy and prognosis in GBM and will be significant for guiding therapeutic strategies and research direction for GBM.
Introduction
As one of the most common malignant brain tumor (1, 2), the 5-year survival rate for glioblastoma (GBM) in patients is less than 5%. Although great progress has occurred in the field of chemotherapy, radiotherapy, and surgical resection, the prognosis for GBM patients is still poor (3–5). Long-non-coding RNA (lncRNA) have been a focus in recent years since they play a nonnegligible role in a variety of biological processes and exert significant influence in human diseases. Particularly, the abnormal expression of lncRNA is tightly connected to the occurrence, prognosis, and survival of patients with cancer (6, 7). Numerous studies have illustrated that certain lncRNAs are aberrantly expressed in GBM tissue, and many of them have been confirmed to be involved in tumor invasion, immune escape and radiation resistance. In recent years, specific lncRNAs were also identified as prognostic biomarkers and therapeutic targets for GBM, while some of them were proposed to be novel indicators for survival prediction in GBM patients.
So, can we create a model on the basis of a multiple-gene signature that can advance the effectiveness of treatment evaluation and prognostic prediction for GBM? In the present study, we collected high-throughput data in TCGA and CGGA databases which were generated by microarrays and next-generation sequencing, then we identified some survival-related differentially expressed genes (DEGs) by analyzing the data of lncRNA expression in GBM, subsequently, a risk model for treatment evaluation and prognostic prediction was established on the basis of the identified gene signature. In this study, five lncRNAs in GBM were revealed notably related with survival independently, and the high-risk score was enriched with items of signaling pathways for oncogenesis and tumor progression. The risk score model based on the five-lncRNAs signature exhibited satisfactory efficacies in predicting the therapeutic efficacy and prognosis in GBM.
Materials and Methods
Publicly Available Clinical Data Sets
RNA-seq (Illumina RNASeqV2, Level 3; Illumina, San Diego, CA) FPKM (fragments per kilobase of transcript per million mapped reads) of GBM were downloaded from The Cancer Genome Atlas (TCGA; https://cancergenome.nih.gov) on October 1, 2019 (The database was updated to January 25, 2018 at that time), including 167 samples. RNA-seq expression data from STAR+RSEM (Illumina Hiseq; mRNAseq_693) of Glioma including 140 GBM samples were downloaded from Chinese Glioma Genome Atlas (CGGA; http://www.cgga.org.cn/index.jsp) on October 1, 2019 (8–10).
Statistical Analyses
All statistical analyses were performed by R. For the TCGA and CGGA dataset, we first converted them to the same gene symbols and selected the common lncRNA, performing log2(FPKM/RSEM+1) transformation. Then the univariate cox regression analysis was used to find the correlation between the expression level of each single gene and patients’ OS. The genes were filtered using a p-value below 0.05. Then, we overlapped the common genes with HR>1 or HR<1 in TCGA and CGGA. We got 3 HR>1 genes and 4 HR<1 genes to perform the multivariate cox regression analysis using the TCGA data set as the training cohort. And we got the risk score formula of the prognostic model. The resulting model according to the formula: In the formula, every single lncRNA represented its expression value, the values of coefficients represented the relative contributions of the five lncRNAs in the multivariate cox regression analysis. TBX5-AS1 is a risk factor while the other four genes may be protective factors. In the calculation of risk score, the higher the value is, the higher the risk of the patient may be. The optimal cut-off value was used to calculate the high-risk and low-risk groups of TCGA, so we got a lower 37% to divide the patients into two groups. And for the CGGA data set, we used the same formula and percentile inherited from the TCGA data set. All the Time-dependent ROC were performed by the R package “survival ROC” (11). In this study, we construct a nomogram based on the risk score groups and clinical traits by the multivariable Cox regression analysis. Subsequently, validations were performed by discrimination and calibration using the R package “rms”. Concordance index (C‐index) was applied to calculate the discrimination of the nomogram. Then we evaluated the nomogram between the prediction probabilities and the observed rates using the calibration curves. The GSEA was performed by the R package “cluster profile” [22455463] (12–14).
Results
Selection of Candidate Genes to Build the Predictive Model
To describe our research clearly, we draw a flow chart of the analysis procedure (Figure S1). Firstly, we selected the common genes of the TCGA and CGGA data sets, including 930 common lncRNAs. Then, we overlapped the genes with HR>1 of TCGA and CGGA, or HR<1 of TCGA and CGGA (Figures 1A, B). And we show the result of the 7 lncRNAs univariate cox regression analysis in Figure 1C. As the results show, the 7 lncRNAs were significantly correlated with the patient’s OS in both TCGA and CGGA. And evaluated H19, LNC00968, and TBX5-AS1 expressions were negatively correlated with the patient’s OS, LNC01545, WDR11-AS1, NDUFA6-DT, and FRY-AS1 were in reverse. Then we used the 7 lncRNAs to perform the multivariate cox regression analysis (Figure 2A) and got the 5-lncRNA signature with C-index 0.65 (Figure 2A). We use the 5-lncRNA signature to divide the patients into two groups based on the optimal cut-off value with 37% low-risk group. The K-M OS curve shows that the high-risk score was significantly correlated with lower OS (HR = 2.12, P < 0.0001; Figure 2B). The AUC of a time-dependent ROC curve was calculated to represent the prognostic capacity of the 5-lncRNA signature (Figure 2C). The AUCs of the 5-lncRNA biomarker prognostic model was 0.67, 0.69, 0.74, 0.70, and 0.70 for the 6, 12, 18, 24, and 36 months OS times, respectively, implicating that the model possesses certain accuracy and classification capability (Figure 2C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/590352/fonc-10-590352-HTML/image_m/fonc-10-590352-g001.jpg)
Figure 1 lncRNA screening in the TCGA and CGGA data sets. (A, B) The overlapping genes with HR>1 or HR<1 between TCGA and CGGA data sets. (C) Univariate cox regression analysis of 7 lncRNAs in the overall survival assessment.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/590352/fonc-10-590352-HTML/image_m/fonc-10-590352-g002.jpg)
Figure 2 The selection of 5 candidate lncRNAs and prognostic capacity assessment. (A) 5 lncRNAs were selected via the multivariate cox regression analysis. (B) K-M OS curve plotting with the 5-lncRNA signature. (C) The prognostic capacity evaluation of the model built by 5-lncRNA signature.
Evaluation the Predictive Model in Validation Cohorts
For independent validation, we assessed the predictive model in CGGA including 140 GBM samples. The risk score of CGGA was calculated in the same way as the TCGA data set. For the survival analysis of the CGGA validation cohort, we also calculated the risk score and divided patients into two groups based on a lower 37% same with TCGA. The result shows the high-risk score was significantly correlated with lower OS in the CGGA cohort (HR = 1.93, p = 0.0013; Figure 3B). Furthermore, the AUCs of time-dependent ROC curves for CGGA cohort were 0.61, 0.69, 0.69, 0.70, and 0.67 at 6, 12, 18, 24, and 36 months (Figure 3A), indicating that the predictive model had capacities for predicting OS in the GBM patients (Figure 3A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/590352/fonc-10-590352-HTML/image_m/fonc-10-590352-g003.jpg)
Figure 3 The evaluation of the predictive model in validation cohorts. (A) AUCs of time-dependent ROC curves for CGGA cohort. (B) The survival analysis of the CGGA validation cohort.
The Five-lncRNAs Signature Was Independent of Clinical Factors
The independent prognostic value of the five-lncRNAs signature was assessed by univariate and multivariate cox regression in 167 TCGA GBM with complete clinical information. As the result of Table 1 shows, age, risk score group, radiation therapy, chemotherapy, and IDH status were correlated with the patient’s OS, and the risk score group of C-index value was almost reached to the radiation therapy and chemotherapy. Then we incorporated the risk score group and all clinical factors to perform the multivariate cox regression. And we find age, radiation therapy, and risk score group were the independent prognostic factors correlated with OS (Table 1). And the risk score group combined with the clinical factors model can increase the C-index from 0.604 to 0.736 (Δ = 0.132), compared to a multivariate model which was mainly based on clinical features (Table 1). In conclusion, the risk score was a significant independent predictor of BCR, indicating that it may contribute to the guidance of treatment decisions in the clinical practice.
Build the Risk Score Combined With Clinical Factors Nomogram
To establish a new model for predicting the patients’ OS in GBM, a nomogram was generated to predict the probability of the 6‐, 12‐, and 18‐months OS in the TCGA cohort. As Table 1 has shown that age, sex, risk score, and radiation therapy were independent prognostic factors (Selecting criteria: significance threshold of log-rank P <0.05), they were therefore embodied in the nomogram (Figure 4A). And the nomogram was visualized with calibration plots. The actual OS was represented by the red line and the blue line represented the predictive OS. Calibration plots showed that the nomogram performed well (Figures 4B–D). And C-index was used to evaluate the combined model (C-index: 0.736), compared to the clinical model (C-index: 0.604) (Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/590352/fonc-10-590352-HTML/image_m/fonc-10-590352-g004.jpg)
Figure 4 The building of a nomogram and its performance on the OS prediction. (A) A nomogram was built with four prognostic factors. (B–D) The OS-predicting performance of the nomogram was evaluated by calibration plots.
Interaction Analysis of Target Proteins and Enrichment of Tumor-Related Signaling Pathway for the Five-lncRNAs Signature
In order to confirm the accuracy of gene screening, we predicted the target proteins of the 5 lncRNAs (Tables S1-S3), and conducted an interaction network analysis of the proteins (Figure S2). Finding that TBX5-AS1 may be related to some oncogenic protein (15–18), while the other 4 lncRNAs may be associated to some tumor suppressor proteins (19–22). Indicating that these lncRNAs may be key genes in GBM. We also performed the differential expression genes (DEGs) analysis using the R package “limma”, and selected the DEGs (adjust. p value <0.05) to perform GSEA using the R package “cluster profile” and Wiki-pathway. The four most significant pathways were showed in Figure S3 and Table S4. The risk group system was accompanied by PI3K-Akt, VEGFA-VEGFR2, TGF-beta, Notch, T-Cell Antigen Receptor signaling pathways, and so on (Figure S3 and Table S4). And these pathways are significantly correlated with tumor progression. Although few lncRNAs have been functionally annotated in GBM, we revealed the relevant signaling pathways of the five lncRNAs through GSEA. Which proved the reliability and distinguishing capability of the formula.
Discussion
LncRNA has been widely recognized as regulators for biological processes related to tumorigenesis and progression in GBM. Jiang et al. stated that blocking Lnc00152 can suppresses glioblastoma malignancy by impairing mesenchymal phenotype through the miR-612/AKT2/NF-κB pathway (23). Lnc-TALC has been proved to be associated with temozolomide (TMZ) resistance induced by AKT signaling pathway in GBM (24). Lnc-SChLAP1 can promote the growth of GBM through stabilization of ACTN4 and activation of NF-κB signaling (25). Growing evidence suggests that lncRNAs can be promising biomarkers for GBM diagnosis and treatment.
However, it is unconvincing to use a single lncRNA as a target due to the significant heterogeneity of GBM (26). More recently, the molecular research and gene exploration of the tumor, especially the statistical analysis research based on big data, has made rapid progress (27, 28). Polygenic research has become the focus of this field. In this present study, after comprehensive analysis of gene expression in GBM from TCGA and CGGA databases, five lncRNAs (LNC01545, WDR11-AS1, NDUFA6-DT, FRY-AS1, TBX5-AS1) were noticed abnormally expressed in GBM, which were also found related with OS. The five-lncRNAs signature-derived risk score presented satisfactory efficacies in predicting the therapeutic efficacy and prognosis in GBM. Additionally, high-risk score was enriched with items of signaling pathways for oncogenesis and tumor progression (29, 30).
The discovered DEGs might be more stably and specifically expressed in GBM based on TCGA and CCGA datasets (including 930 common lncRNAs), compared to just one dataset. Then univariate and multivariate Cox regression analysis were conducted on the DEGs, and five important genes were found associated with survival independently. Results show that TBX5-AS1 expressions were negatively correlated with the patient’s OS, which may function as oncogenes, whereas LNC01545, WDR11-AS1, NDUFA6-DT, and FRY-AS1 were in reverse, which may function as suppressor gene. A risk score model was then established on the basis of the signature of these five DEGs. Studies have indicated that gene signature containing multiple gene components can be more robust and convictive in the prediction of prognosis, while the usage of some single gene may result in instability and predictive bias. In the present study, the risk score based on the five-lncRNAs signature was effective and reproducible in predicting prognosis of the GBM patients from both TCGA and CCGA datasets. Additionally, the risk score model could also predict the radiotherapy response of GBM patient. Generally, the multi-gene signature-derived risk score model may be promising and valid in treatment evaluation and prognostic prediction in GBM.
LncRNA can promote or inhibit GBM by regulating signal transduction pathways, forming sponge adsorption effect, regulating the characteristics of glioma stem cells, regulating hypoxia response, and angiogenesis (31–33). The potential functions of the five genes may partially contribute to the prognostic prediction, but the underlying mechanisms upon these genes in GBM remains to be investigated (34–36). Of the five lncRNAs, apart from TBX5-AS1, which has been to reported be related with unfavorable prognosis in non-small-cell lung carcinoma (NSCLC) (37, 38), the other four has never been explored in cancer, and none of them has ever been explored in GBM. In this case, we believe that it is promising to explore the potential mechanisms of these genes (especially the TBX5-AS1 which acts as a risk factor in our model) related to GBM.
Generally speaking, clinical and pathological classification determine the treatment and prognosis of GBM. Independent of the traditional method, the value of this model lies in its clinical guidance for therapeutic strategies, which may help improve the prognosis in GBM patients. Intensive therapy should be applied to the patients with high-risk scores, while those with low risk scores should avoid excessive treatments that may cause therapeutic toxicities and deterioration. Although the values of AUCs for five genes were close to or more than 0.7, the accuracy and classification capability of the model are still not high enough, the interaction analysis of target proteins and enrichment of tumor-related signaling pathway proved the accuracy of gene screening and the reliability and distinguishing capability of the model. This also enlightens us to make some in-depth exploration on the mechanisms in the further study. On the other hand, compared with those genes that have just been screened and have not yet been verified, it may be another reliable choice to build models based on recognized genes. Moreover, except for TCGA, CGGA, and GEO, it is highly recommended to introduced other databases (e.g. NoncoRNA) in our further studies (39). Despite a long way ahead for the clinical application of this model, the prospects are considerable and worthy of further exploration.
Conclusions
In conclusion, the five-lncRNAs signature-derived risk score presented satisfactory efficacies in predicting the therapeutic efficacy and prognosis in GBM and will be significant for guiding therapeutic strategies and research direction for GBM. Since most of the genes we put forward have not been detailedly researched, confirmatory experiments are necessary in the further studies.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: The Cancer Genome Atlas (TCGA, https://cancergenome.nih.gov), and Chinese Glioma Genome Atlas (CGGA, http://www.cgga.org.cn/index.jsp).
Author Contributions
XN, JS, and HL designed the study. LM, TF, and TZ collated the data. XN and JJ conducted data analyses and drafted the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by Tianjin Major Science and Technology Projects (17ZXMFSY00200), Tianjin Science and Technology Development Strategy Research Projects (18ZLZXZF00740).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.590352/full#supplementary-material
Supplementary Figure 1 | The flow chart of the analysis procedure for this study.
Supplementary Figure 2 | Figure S2 Interaction network analysis of the target proteins for the five lncRNAs. The red circles represent the up-regulated proteins and the green represent the down regulated proteins.
Supplementary Figure 3 | Identification of 5-lncRNAs signature associated biological signaling pathway. PI3K-Akt, VEGFA-VEGFR2, TGF-beta, Notch, T-Cell Antigen Receptor signaling pathways were the main enrichment pathways.
Supplementary Table 1 | Prediction of target proteins for FRY-AS1, LNC01545, WDR11-AS1, and NDUFA6-DT.
Supplementary Table 2 | Prediction of target proteins for TBX5-AS1.
Supplementary Table 3 | Target protein prediction and related biological effects of the five lncRNAs.
Supplementary Table 4 | The related signaling pathways of the 5-lncRNAs signature.
Abbreviations
GBM, Glioblastoma; TCGA, The Cancer Genome Atlas; CGGA, Chinese Glioma Genome Atlas; LncRNA, Long non-coding RNA; DEGs, Differentially expressed genes; GSEA, Gene set enrichment analysis; OS, Overall survival.
References
1. Scott JG, Berglund A, Schell MJ, Mihaylov I, Fulp WJ, Yue B, et al. A genome-based model for adjusting radiotherapy dose (GARD): a retrospective, cohort-based study. Lancet Oncol (2017) 18:202–11. doi: 10.1016/S1470-2045(16)30648-9
2. Wang C, Li J, Sinha S, Peterson A, Grant GA, Yang F. Mimicking brain tumor-vasculature microanatomical architecture via co-culture of brain tumor and endothelial cells in 3D hydrogels. Biomaterials (2019) 202:35–44. doi: 10.1016/j.biomaterials.2019.02.024
3. Ostrom QT, Gittleman H, Fulop J, Liu M, Blanda R, Kromer C, et al. CBTRUS Statistical Report: Primary Brain and Central Nervous System Tumors Diagnosed in the United States in 2008-2012. Neuro Oncol (2015) 17(Suppl 4):iv1–1iv62. doi: 10.1093/neuonc/nov189
4. Stupp R, Taillibert S, Kanner AA, Kesari S, Steinberg DM, Toms SA, et al. Maintenance Therapy With Tumor-Treating Fields Plus Temozolomide vs Temozolomide Alone for Glioblastoma: A Randomized Clinical Trial. JAMA (2015) 314:2535–43. doi: 10.1001/jama.2015.16669
5. Tan AC, Ashley DM, López GY, Malinzak M, Friedman HS, Khasraw M. Management of glioblastoma: State of the art and future directions. CA Cancer J Clin (2020) 70:299–312. doi: 10.3322/caac.21613
6. Hu Q, Ye Y, Chan LC, Li Y, Liang K, Lin A, et al. Oncogenic lncRNA downregulates cancer cell antigen presentation and intrinsic tumor suppression. Nat Immunol (2019) 20:835–51. doi: 10.1038/s41590-019-0400-7
7. Peng Z, Liu C, Wu M. New insights into long noncoding RNAs and their roles in glioma. Mol Cancer (2018) 17:61. doi: 10.1186/s12943-018-0812-2
8. Jalili V, Matteucci M, Morelli MJ, Masseroli M. MuSERA: Multiple Sample Enriched Region Assessment. Brief Bioinform (2017) 18:367–81. doi: 10.1093/bib/bbw029
9. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol (2014) 15:550. doi: 10.1186/s13059-014-0550-8
10. Yin T, Cook D, Lawrence M. ggbio: an R package for extending the grammar of graphics for genomic data. Genome Biol (2012) 13:R77. doi: 10.1186/gb-2012-13-8-r77
11. Heagerty PJ, Lumley T, Pepe MS. Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics (2000) 56:337–44. doi: 10.1111/j.0006-341x.2000.00337.x
12. Ellis SE, Collado-Torres L, Jaffe A, Leek JT. Improving the value of public RNA-seq expression data by phenotype prediction. Nucleic Acids Res (2018) 46:e54. doi: 10.1093/nar/gky102
13. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43:e47. doi: 10.1093/nar/gkv007
14. Seyednasrollah F, Laiho A, Elo LL. Comparison of software packages for detecting differential expression in RNA-seq studies. Brief Bioinform (2015) 16:59–70. doi: 10.1093/bib/bbt086
15. Deininger MH, Wybranietz WA, Graepler FT, Lauer UM, Meyermann R, Schluesener HJ. Endothelial endostatin release is induced by general cell stress and modulated by the nitric oxide/cGMP pathway. FASEB J (2003) 17:1267–76. doi: 10.1096/fj.02-1118com
16. Yang J, Fan J, Li Y, Li F, Chen P, Fan Y, et al. Genome-wide RNAi screening identifies genes inhibiting the migration of glioblastoma cells. PloS One (2013) 8:e61915. doi: 10.1371/journal.pone.0061915
17. Jiang R, Choi W, Hu L, Gerner EW, Hamilton SR, Zhang W. Activation of polyamine catabolism by N1, N11-diethylnorspermine alters the cellular localization of mTOR and downregulates mTOR protein level in glioblastoma cells. Cancer Biol Ther (2007) 6:1644–8. doi: 10.4161/cbt.6.10.4800
18. Zhang P, Guo Z, Zhang Y, Gao Z, Ji N, Wang D, et al. A preliminary quantitative proteomic analysis of glioblastoma pseudoprogression. Proteome Sci (2015) 13:12. doi: 10.1186/s12953-015-0066-5
19. Izaguirre DI, Zhu W, Hai T, Cheung HC, Krahe R, Cote GJ. PTBP1-dependent regulation of USP5 alternative RNA splicing plays a role in glioblastoma tumorigenesis. Mol Carcinog (2012) 51:895–906. doi: 10.1002/mc.20859
20. Xu J, Hou X, Pang L, Sun S, He S, Yang Y, et al. Identification of Dysregulated Competitive Endogenous RNA Networks Driven by Copy Number Variations in Malignant Gliomas. Front Genet (2019) 10:1055. doi: 10.3389/fgene.2019.01055
21. Han B, Wang R, Chen Y, Meng X, Wu P, Li Z, et al. QKI deficiency maintains glioma stem cell stemness by activating the SHH/GLI1 signaling pathway. Cell Oncol (Dordr) (2019) 42:801–13. doi: 10.1007/s13402-019-00463-x
22. Gupta MK, Polisetty RV, Sharma R, Ganesh RA, Gowda H, Purohit AK, et al. Altered transcriptional regulatory proteins in glioblastoma and YBX1 as a potential regulator of tumor invasion. Sci Rep (2019) 9:10986. doi: 10.1038/s41598-019-47360-9
23. Cai J, Zhang J, Wu P, Yang W, Ye Q, Chen Q, et al. Blocking LINC00152 suppresses glioblastoma malignancy by impairing mesenchymal phenotype through the miR-612/AKT2/NF-κB pathway. J Neurooncol (2018) 140:225–36. doi: 10.1007/s11060-018-2951-0
24. Wu P, Cai J, Chen Q, Han B, Meng X, Li Y, et al. Lnc-TALC promotes O6-methylguanine-DNA methyltransferase expression via regulating the c-Met pathway by competitively binding with miR-20b-3p. Nat Commun (2019) 10:2045. doi: 10.1038/s41467-019-10025-2
25. Ji J, Xu R, Ding K, Bao G, Zhang X, Huang B, et al. Long Noncoding RNA SChLAP1 Forms a Growth-Promoting Complex with HNRNPL in Human Glioblastoma through Stabilization of ACTN4 and Activation of NF-κB Signaling. Clin Cancer Res (2019) 25:6868–81. doi: 10.1158/1078-0432.CCR-19-0747
26. Gong Z, Hong F, Wang H, Zhang X, Chen J. An eight-mRNA signature outperforms the lncRNA-based signature in predicting prognosis of patients with glioblastoma. BMC Med Genet (2020) 21:56. doi: 10.1186/s12881-020-0992-7
27. Budczies J, Kluck K, Walther W, Stein U. Decoding and targeting the molecular basis of MACC1-driven metastatic spread: Lessons from big data mining and clinical-experimental approaches. Semin Cancer Biol (2020) 60:365–79. doi: 10.1016/j.semcancer.2019.08.010
28. Kuenzi BM, Ideker T. A census of pathway maps in cancer systems biology. Nat Rev Cancer (2020) 20:233–46. doi: 10.1038/s41568-020-0240-7
29. Starmans MH, Lieuwes NG, Span PN, Haider S, Dubois L, Nguyen F, et al. Independent and functional validation of a multi-tumour-type proliferation signature. Br J Cancer (2012) 107:508–15. doi: 10.1038/bjc.2012.269
30. Subramanian J, Simon R. What should physicians look for in evaluating prognostic gene-expression signatures. Nat Rev Clin Oncol (2010) 7:327–34. doi: 10.1038/nrclinonc.2010.60
31. Lin A, Li C, Xing Z, Hu Q, Liang K, Han L, et al. The LINK-A lncRNA activates normoxic HIF1α signalling in triple-negative breast cancer. Nat Cell Biol (2016) 18:213–24. doi: 10.1038/ncb3295
32. Wang Y, Wu S, Zhu X, Zhang L, Deng J, Li F, et al. LncRNA-encoded polypeptide ASRPS inhibits triple-negative breast cancer angiogenesis. J Exp Med (2020) 217:jem.20190950. doi: 10.1084/jem.20190950
33. Yang X, Xiao Z, Du X, Huang L, Du G. Silencing of the long non-coding RNA NEAT1 suppresses glioma stem-like properties through modulation of the miR-107/CDK6 pathway. Oncol Rep (2017) 37:555–62. doi: 10.3892/or.2016.5266
34. Liu C, Zhang Y, She X, Fan L, Li P, Feng J, et al. A cytoplasmic long noncoding RNA LINC00470 as a new AKT activator to mediate glioblastoma cell autophagy. J Hematol Oncol (2018) 11:77. doi: 10.1186/s13045-018-0619-z
35. Pastori C, Kapranov P, Penas C, Peschansky V, Volmar CH, Sarkaria JN, et al. The Bromodomain protein BRD4 controls HOTAIR, a long noncoding RNA essential for glioblastoma proliferation. Proc Natl Acad Sci U S A (2015) 112:8326–31. doi: 10.1073/pnas.1424220112
36. Tan SK, Pastori C, Penas C, Komotar RJ, Ivan ME, Wahlestedt C, et al. Serum long noncoding RNA HOTAIR as a novel diagnostic and prognostic biomarker in glioblastoma multiforme. Mol Cancer (2018) 17:74. doi: 10.1186/s12943-018-0822-0
37. Qu QH, Jiang SZ, Li XY. LncRNA TBX5-AS1 Regulates the Tumor Progression Through the PI3K/AKT Pathway in Non-Small Cell Lung Cancer. Onco Targets Ther (2020) 13:7949–61. doi: 10.2147/OTT.S255195
38. Shih JH, Chen HY, Lin SC, Yeh YC, Shen R, Lang YD, et al. Integrative analyses of noncoding RNAs reveal the potential mechanisms augmenting tumor malignancy in lung adenocarcinoma. Nucleic Acids Res (2020) 48:1175–91. doi: 10.1093/nar/gkz1149
Keywords: lncRNA, glioblastoma, prognosis, The Cancer Genome Atlas, Chinese Glioma Genome Atlas
Citation: Niu X, Sun J, Meng L, Fang T, Zhang T, Jiang J and Li H (2020) A Five-lncRNAs Signature-Derived Risk Score Based on TCGA and CGGA for Glioblastoma: Potential Prospects for Treatment Evaluation and Prognostic Prediction. Front. Oncol. 10:590352. doi: 10.3389/fonc.2020.590352
Received: 01 August 2020; Accepted: 10 November 2020;
Published: 17 December 2020.
Edited by:
Jing Sun, George Washington University, United StatesReviewed by:
Manqi Zhang, Duke University, United StatesJinquan Cai, Harbin Medical University, China
Copyright © 2020 Niu, Sun, Meng, Fang, Zhang, Jiang and Li. 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: Huanming Li, lihuanming0716@163.com; Jipeng Jiang, 13502034360@163.com
†These authors have contributed equally to this work