- Open Access
Identification of intrinsic subtype-specific prognostic microRNAs in primary glioblastoma
- Rui Li†1,
- Kaiming Gao†1,
- Hui Luo1,
- Xiefeng Wang1,
- Yan Shi1,
- Qingsheng Dong1,
- WenKang Luan1 and
- Yongping You1Email author
© Li et al.; licensee BioMed Central Ltd. 2014
Received: 23 December 2013
Accepted: 13 January 2014
Published: 19 January 2014
Glioblastoma multiforme (GBM) is the most malignant type of glioma. Integrated classification based on mRNA expression microarrays and whole–genome methylation subdivides GBM into five subtypes: Classical, Mesenchymal, Neural, Proneural-CpG island methylator phenotype (G-CIMP) and Proneural-non G-CIMP. Biomarkers that can be used to predict prognosis in each subtype have not been systematically investigated.
In the present study, we used Cox regression and risk-score analysis to construct respective prognostic microRNA (miRNA) signatures in the five intrinsic subtypes of primary glioblastoma in The Cancer Genome Atlas (TCGA) dataset.
Patients who had high-risk scores had poor overall survival compared with patients who had low-risk scores. The prognostic miRNA signature for the Mesenchymal subtype (four risky miRNAs: miR-373, miR-296, miR-191, miR-602; one protective miRNA: miR-223) was further validated in an independent cohort containing 41 samples.
We report novel diagnostic tools for deeper prognostic sub-stratification in GBM intrinsic subtypes based upon miRNA expression profiles and believe that such signature could lead to more individualized therapies to improve survival rates and provide a potential platform for future studies on gene treatment for GBM.
Glioblastoma multiforme (GBM) is the most lethal type of glioma in adults due to its poor prognosis and resistance to clinical therapy . Despite the continuous improvements in radiotherapy, chemotherapy, and surgical treatments, the prognosis remains unsatisfactory . Obstacles result, in part, from the heterogeneous nature and variable genetic aberrations among affected individuals. Developments of molecular biology offer a potential gene therapy for GBM. With the deepening of the research, a growing number of potential gene therapy targets have been identified, which calls for an urgent and objective classification based on molecular signatures rather than the histopathologic classification system . Currently, many databases, such as The Cancer Genome Atlas (TCGA) Research Network, are being continuously improved and provide us with a comprehensive catalogue of genomic information . Several molecular-based classification systems, such as mRNA expression-based and DNA methylation-based systems, have identified abnormalities driving tumorigenesis and correlated clinical data of samples.
The generally accepted molecular classification divides GBM into five subtypes: Proneural-G-CIMP, Proneural-non G-CIMP, Neural, Classical and Mesenchymal. Each subtype exhibits distinct biological behavior and characteristic gene expression [5, 6]. However, biomarkers, which can be used to predict prognosis in each subtype, have not been systematically investigated. In our study, we identified significant miRNAs associated with clinical outcomes in each subtype. The prognostic miRNA signature for the Mesenchymal subtype was validated in an independent cohort containing 41 samples. The findings reveal the potential prognostic meaning of each subtype and provide an alternative individual treatment for patients with GBM.
The miRNA expression microarray data (Level 3) for GBM samples were downloaded from The Cancer Genome Atlas (TCGA) database (http://cancergenome.nih.gov). The corresponding clinical data with the molecular subtype annotations were obtained from Ref . In total, 448 primary GBM cases with molecular annotations (120 cases of Classical subtype, 141 cases of Mesenchymal subtype, 73 cases of Neural subtype, 29 cases Proneural-G-CIMP and 85 cases of Proneural-non G-CIMP) and miRNA expression microarray data were included in our analysis. A total of 41 Mesenchymal primary GBM from the Chinese Glioma Genome Atlas (CGGA) were used as a validation cohort by qRT-PCR assay.
Total RNA (tRNA) was extracted from frozen tissues using the mirVana miRNA Isolation Kit (Ambion, Inc., Austin, Tex), and its concentration and quality were determined with a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, Del) .
Real-time quantification of miRNAs by stem-loop RT-PCR
For the TaqMan-based real-time reverse transcription polymerase chain reaction (RT-PCR) assays, an ABI 7300 HT Sequence Detection system (Applied Biosystems, Foster City, CA) was used. All primers and probes of hsa-miR-373, hsa-miR-296, hsa-miR-191, hsa-miR-602 and hsa-miR-223, and RNU6B endogenous controls for TaqMan miRNA assays were purchased from Applied Biosystems. Real-time PCR was performed as described by Ref. . The relative gene expression was calculated via a 2−ΔΔCt method .
Based on the generally accepted molecular classification of GBM (the classical, mesenchymal, neural, proneural-G-CIMP and proneural-non G-CIMP subtypes), we first analyzed the internal prognosis stratification of the five subtypes based on the MGMT promoter methylator phenotype. Kaplan-Meier survival analysis was used to estimate the survival distributions. The log-rank test was used to assess the statistical significance between the stratified survival groups using GraphPad Prism 6.0 statistical software. P < 0.05 was considered significant.
The expression level of each miRNA (n = 470) was assessed by Cox regression analysis using the BRB array tools package developed by Richard Simon and the BRB-ArrayTools Development Team . Permutation Tests were performed with 10,000 permutations to select genes that were significantly associated with overall survival. For the Classical subtype, seven miRNAs with permutation P-values < 0.05 were selected as the candidate genes; for Proneural-G-CIMP and Proneural-non G-CIMP subtypes, three and ten miRNAs with permutation P-values <0.05 were selected as candidate genes. For Neural subtypes, eight miRNAs with permutation P-values < 0.015 were selected as candidate genes, and for Mesenchymal subtypes, five miRNAs with permutation P-values < 0.01 were selected as candidate genes. The significant miRNAs were divided into risky and protective types. Risky miRNAs were defined as miRNAs with a hazard ratio for death greater than 1. In contrast, protective miRNAs were defined based on a hazard ratio for death less than 1.
For each subtype, patients were divided into high-risk and low-risk groups according to the cutoff value (median risk score); patients in high-risk group are expected to have a poor outcome. The Kaplan-Meier method was used to estimate overall survival time for the two groups. Differences in survival times were analyzed using the two-sided log rank test. The significant miRNAs remained the same in validation set.
SPSS 13.0 for Windows (SPSS, Inc., Chicago, Ill) was used to conduct survival analyses. All tests were 2-tailed, and the significance level was set at P <0 .05.
Identification prognostic microRNA signatures in the five molecular subtypes of primary glioblastoma
In another subtype of GBM, Neural, we filtered 8 miRNAs (one risky miRNAs: hsa-miR-222 and seven protective miRNAs: hsa-miR-422a, hsa-miR-662, hsa-miR-566, hsa-miR-24, hsa-miR-370, hsa-miR-492, hsa-miR-629) as signature genes for predicting the patient outcomes (p<0.015) (Figure 1B). Our result revealed that the prognoses of patients belonging to the high-risk group (n = 36) were worse than those of the low-risk group (n = 37) (Figure 2B).
A total of 5 miRNAs (four risky miRNAs: miR-373, miR-296, miR-191, miR-602; one protective miRNA: miR-223) were selected to identify the overall survival of patients in Mesenchymal subtype GBM (Figure 1C). Compared with the high-risk group, the low-risk group exhibited a better prognosis (p<0.001) (Figure 2C).
As for Proneural-G-CIMP and Proneural-non G-CIMP subtype GBM, three (one risky miRNAs: hsa-miR-582 and two protective miRNAs: hsa-miR-130a, hsa-miR-195) and ten miRNAs (four risky miRNAs: hsa-miR-335, hsa-miR-34a, hsa-miR-581, hsa-miR-21 and six protective miRNAs: hsa-miR-361, hsa-miR-145, hsa-miR-143, hsa-miR-378, hsa-miR-182, hsa-miR-183) were filtered for intrinsic prognostic analysis, respectively (Figure 1D-E). In both subtypes, the patients in the high-risk group exhibited shorter survival than patients in the low-risk group, as shown in Figure 2D-E.
Validation of the prognostic value of the miRNA signatures of mesenchymal subtype in an independent cohort
Glioblastoma multiforme (GBM) is the most malignant type of glioma. People who suffer from GBM exhibit poor prognosis. Survival ranges from one week to longer than three years, but most patients only survive approximately one year [13, 14]. Decades of effort in clinical treatment have not satisfactorily improved prognosis. To improve the situation, researchers have begun to use molecular biological technologies to identify new approaches. Data on some important biomarkers in the progression of GBM, such as loss of heterozygosity (LOH) 10q, EGFR amplification, TP53 mutations, and PTEN mutations, have been collected and considered as certain gene targets [15, 16]. To improve target efficiency, classifications of GBM have been established based on different gene signatures, including mRNA expression [17, 18], microRNA expression, and methylation [19, 20]. The generally accepted gene expression-based molecular classification by Phillips, H.S. et al. divides GBM into the Proneural, Neural, Classical and Mesenchymal subtypes . Further studies divide the Proneural subtype of GBM into Proneural-G-CIMP and Proneural-non-G-CIMP subtypes based on the CPG Island methylator phenotype [22, 23].
However, due to the complex mechanism of development and progression of GBM and the diversity of gene mutation, prognoses vary in samples of each subtype. The isocitrate dehydrogenase 1 (IDH1) mutation was reported to be a biomarker in the prediction of clinical outcomes for patients with GBM, and patients with IDH1 mutations exhibit better outcome than those with wild-type IDH1 in gliomas. However, the genomic alteration was frequent in grade II and grade III glial tumors but rare in primary glioblastoma (pGBM) [24, 25]. Another prognostic biomarker, MGMT DNA methylation, was reported by Brennan CW et al. . Only in classical subtype GBM, MGMT DNA methylation may be a predictive biomarker for treatment response. Patients with MGMT DNA methylation exhibit significantly better outcomes than atients with MGMT DNA unmethylation. However, the other four subtypes, Proneural-G-CIMP, Proneural-non G-CIMP, Neural, and Mesenchymal, exhibit no significant difference in clinical prognosis based on the MGMT DNA methylation biomarker. This prompted us to identify molecular prognostic markers that exist in all subtypes of GBM. In present study, we choose miRNAs as predictive biomarkers for intrinsic prognostic stratification. By using permutation tests and Cox regression analyses, significant genes were filtered and risk scores were calculated to help us divide patients into low- and high-risk groups. Our results revealed that in all five subtypes, patients belonging to the low-risk groups had significant longer overall survival than those in the high-risk group. Additionally, we validated the results in an independent database (CGGA) and obtained similar results. Although significant miRNAs were selected as predictive biomarkers for each subtype, there are many other factors, including age, sex, race and different treatment strategies, which may affect clinical outcomes and should also be taken into consideration in the treatment of glioblastoma.
In summary, our results indicate that miRNAs can be considered biomarkers for prognosis stratification in each subtype of glioblastoma, which could provide a novel approach in the evaluation of the prognosis and selection of the gene therapy targets.
This work was supported by National High Technology Research and Development Program 863 (2012AA02A508), Jiangsu Provincial Special Program of Medical Science (BL2012028), China National Natural Scientific Fund (91229121 and 81172389).
- Penas-Prado M, Armstrong TS, Gilbert MR: Glioblastoma. Handb Clin Neurol. 2012, 105: 485-506.View ArticlePubMedGoogle Scholar
- Omuro A, DeAngelis LM: Glioblastoma and other malignant gliomas: a clinical review. JAMA. 2013, 310: 1842-1850. 10.1001/jama.2013.280319.View ArticlePubMedGoogle Scholar
- Belden CJ, Valdes PA, Ran C, Pastel DA, Harris BT, Fadul CE, Israel MA, Paulsen K, Roberts DW: Genetics of glioblastoma: a window into its imaging and histopathologic variability. Radiographics. 2011, 31: 1717-1740. 10.1148/rg.316115512.View ArticlePubMedGoogle Scholar
- Rao SA, Santosh V, Somasundaram K: Genome-wide expression profiling identifies deregulated miRNAs in malignant astrocytoma. Mod Pathol. 2010, 23: 1404-1417. 10.1038/modpathol.2010.135.View ArticlePubMedGoogle Scholar
- Ashton CH, Rawlins MD, Tyrer SP: Buspirone in benzodiazepine withdrawal. Br J Psychiatry. 1991, 158: 283-284.PubMedGoogle Scholar
- Brennan CW, Verhaak RG, McKenna A, Campos B, Noushmehr H, Salama SR, Zheng S, Chakravarty D, Sanborn JZ, Berman SH, et al: The somatic genomic landscape of glioblastoma. Cell. 2013, 155: 462-477. 10.1016/j.cell.2013.09.034.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang W, Zhang J, Yan W, You G, Bao Z, Li S, Kang C, Jiang C, You Y, Zhang Y, et al: Whole-genome microRNA expression profiling identifies a 5-microRNA signature as a prognostic biomarker in Chinese patients with primary glioblastoma multiforme. Cancer. 2013, 119: 814-824. 10.1002/cncr.27826.View ArticlePubMedGoogle Scholar
- Chen SY, Su YH, Wu SF, Sha T, Zhang YP: Mitochondrial diversity and phylogeographic structure of Chinese domestic goats. Mol Phylogenet Evol. 2005, 37: 804-814. 10.1016/j.ympev.2005.06.014.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- Simon R, Lam A, Li MC, Ngan M, Menenzes S, Zhao Y: Analysis of gene expression data using BRB-ArrayTools. Canc Informat. 2007, 3: 11-17.Google Scholar
- Dave SS, Wright G, Tan B, Rosenwald A, Gascoyne RD, Chan WC, Fisher RI, Braziel RM, Rimsza LM, Grogan TM, et al: Prediction of survival in follicular lymphoma based on molecular features of tumor-infiltrating immune cells. N Engl J Med. 2004, 351: 2159-2169. 10.1056/NEJMoa041869.View ArticlePubMedGoogle Scholar
- Zhao Q, Sun J: Cox survival analysis of microarray gene expression data using correlation principal component regression. Stat Appl Genet Mol Biol. 2007, 6: Article16Google Scholar
- Ohgaki H, Kleihues P: Population-based studies on incidence, survival rates, and genetic alterations in astrocytic and oligodendroglial gliomas. J Neuropathol Exp Neurol. 2005, 64: 479-489.PubMedGoogle Scholar
- Ohgaki H, Dessen P, Jourde B, Horstmann S, Nishikawa T, Di Patre PL, Burkhard C, Schuler D, Probst-Hensch NM, Maiorka PC, et al: Genetic pathways to glioblastoma: a population-based study. Cancer Res. 2004, 64: 6892-6899. 10.1158/0008-5472.CAN-04-1337.View ArticlePubMedGoogle Scholar
- Verhaak RG, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, Miller CR, Ding L, Golub T, Mesirov JP, et al: Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer cell. 2010, 17: 98-110. 10.1016/j.ccr.2009.12.020.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan W, Zhang W, Jiang T: Oncogene addiction in gliomas: implications for molecular targeted therapy. J Exp Clin Cancer Res. 2011, 30: 58-10.1186/1756-9966-30-58.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan W, Zhang W, You G, Zhang J, Han L, Bao Z, Wang Y, Liu Y, Jiang C, Kang C, et al: Molecular classification of gliomas based on whole genome gene expression: a systematic report of 225 samples from the Chinese Glioma Cooperative Group. Neuro Oncol. 2012, 14: 1432-1440. 10.1093/neuonc/nos263.PubMed CentralView ArticlePubMedGoogle Scholar
- Li A, Walling J, Ahn S, Kotliarov Y, Su Q, Quezado M, Oberholtzer JC, Park J, Zenklusen JC, Fine HA: Unsupervised analysis of transcriptomic profiles reveals six glioma subtypes. Cancer Res. 2009, 69: 2091-2099.PubMed CentralView ArticlePubMedGoogle Scholar
- Noushmehr H, Weisenberger DJ, Diefes K, Phillips HS, Pujara K, Berman BP, Pan F, Pelloski CE, Sulman EP, Bhat KP, et al: Identification of a CpG island methylator phenotype that defines a distinct subgroup of glioma. Cancer cell. 2010, 17: 510-522. 10.1016/j.ccr.2010.03.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang W, Yan W, You G, Bao Z, Wang Y, Liu Y, You Y, Jiang T: Genome-wide DNA methylation profiling identifies ALDH1A3 promoter methylation as a prognostic predictor in G-CIMP- primary glioblastoma. Cancer Lett. 2013, 328: 120-125. 10.1016/j.canlet.2012.08.033.View ArticlePubMedGoogle Scholar
- Phillips HS, Kharbanda S, Chen R, Forrest WF, Soriano RH, Wu TD, Misra A, Nigro JM, Colman H, Soroceanu L, et al: Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer cell. 2006, 9: 157-173. 10.1016/j.ccr.2006.02.019.View ArticlePubMedGoogle Scholar
- Theeler BJ, Yung WK, Fuller GN, De Groot JF: Moving toward molecular classification of diffuse gliomas in adults. Neurology. 2012, 79: 1917-1926. 10.1212/WNL.0b013e318271f7cb.PubMed CentralView ArticlePubMedGoogle Scholar
- Li B, Senbabaoglu Y, Peng W, Yang ML, Xu J, Li JZ: Genomic estimates of aneuploid content in glioblastoma multiforme and improved classification. Clin Cancer Res. 2012, 18: 5595-5605. 10.1158/1078-0432.CCR-12-1427.PubMed CentralView ArticlePubMedGoogle Scholar
- Cohen AL, Holmen SL, Colman H: IDH1 and IDH2 mutations in gliomas. Curr Neurol Neurosci Rep. 2013, 13: 345-PubMed CentralView ArticlePubMedGoogle Scholar
- Ducray F, Idbaih A, Wang XW, Cheneau C, Labussiere M, Sanson M: Predictive and prognostic factors for gliomas. Expert Rev Anticancer Ther. 2011, 11: 781-789. 10.1586/era.10.202.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.