- Open Access
Identification of intrinsic subtype-specific prognostic microRNAs in primary glioblastoma
Journal of Experimental & Clinical Cancer Researchvolume 33, Article number: 9 (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.
Using these significant miRNAs, a risk-score formula for predicting survival was developed based on a linear combination of the gene expression level (expr) weighted by the regression coefficient derived from the univariate Cox regression analysis (β) [11, 12]. The risk score for each patient was calculated as follows:
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 classical subtype GBM, seven miRNAs (five risky miRNAs: hsa-miR-26a, hsa-miR-767-3p, hsa-miR-153, hsa-miR-31, hsa-miR-222, and two protective miRNAs: hsa-miR-654 and hsa-miR-422b) were found to be significantly correlated with clinical outcomes (p<0.05). The risk score was calculated based on the expression of these genes and was obtained in order to predict patient survival. A total of 120 samples were divided into a high-risk group (n = 60) and a low-risk group (n = 60) according to their risk score. The heatmap shows that, protective miRNAs exhibit high expression in low-risk group, while the risky miRNAs have high expression in high-risk group (Figure 1A). As shown in Figure 2A, the patients in the high-risk group suffered obvious worse overall survival than those in the low-risk group (P < 0.0001).
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
The prognostic miRNA signature for the Mesenchymal subtype was further validated in an independent cohort containing 41 Mesenchymal pGBM samples from CGGA. The expression levels of the five prognostic miRNAs (hsa-miR-373, hsa-miR-296, hsa-miR-191, hsa-miR-602 and hsa-miR-223) in the Mesenchymal subtype were analyzed by qRT-PCR in the independent cohort. The risk-score in each sample of validation cohort were calculated according the above formula. Furthermore, patients with Mesenchymal subtype GBM were divided into two groups based on their risk-scores. Similar results were obtained: the low-risk group exhibited extended survival, whereas the high-risk group exhibited shorter survival (Figure 3).
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.
Penas-Prado M, Armstrong TS, Gilbert MR: Glioblastoma. Handb Clin Neurol. 2012, 105: 485-506.
Omuro A, DeAngelis LM: Glioblastoma and other malignant gliomas: a clinical review. JAMA. 2013, 310: 1842-1850. 10.1001/jama.2013.280319.
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.
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.
Ashton CH, Rawlins MD, Tyrer SP: Buspirone in benzodiazepine withdrawal. Br J Psychiatry. 1991, 158: 283-284.
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.
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.
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.
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.
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.
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.
Zhao Q, Sun J: Cox survival analysis of microarray gene expression data using correlation principal component regression. Stat Appl Genet Mol Biol. 2007, 6: Article16
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Cohen AL, Holmen SL, Colman H: IDH1 and IDH2 mutations in gliomas. Curr Neurol Neurosci Rep. 2013, 13: 345-
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.
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).
The authors declare that they have no competing interests.
LR, GKM, WXF performed the experiments. LH, SY, DQS, LWK analyzed the data. LR Wrote the manuscript. YYP, LR conceived and designed the experiments. All authors read and approved the final manuscript.
Rui Li, Kaiming Gao contributed equally to this work.