A glycolysis-based ten-gene signature correlates with the clinical outcome,molecular subtype and IDH1 mutation in glioblastoma

Reprogrammed metabolism is a hallmark of cancer. Glioblastoma (GBM) tumor cells predominantly utilize aerobic glycolysis for the biogenesis of energy and intermediate nutrients. However, in GBM, the clinical significance of glycolysis and its underlying relations with the molecular features such as IDH1 mutation and subtype have not been elucidated yet. Herein, based on glioma datasets including TCGA (The Cancer Genome Atlas), REMBRANDT (Repository for Molecular Brain Neoplasia Data) and GSE16011, we established a glycolytic gene expression signature score (GGESS) by incorporating ten glycolytic genes. Then we performed survival analyses and investigated the correlations between GGESS and IDH1 mutation as well as the molecular subtypes in GBMs. The results showed that GGESS independently predicted unfavorable prognosis and poor response to chemotherapy of GBM patients. Notably, GGESS was high in GBMs of mesenchymal subtype but low in IDH1-mutant GBMs. Furthermore, we found that the promoter regions of tumor-promoting glycolytic genes were hypermethylated in IDH1-mutant GBMs. Finally, we found that high GGESS also predicted poor prognosis and poor response to chemotherapy when investigating IDH1-wildtype GBM patients only. Collectively, glycolysis represented by GGESS predicts unfavorable clinical outcome of GBM patients and is closely associated with mesenchymal subtype and IDH1 mutation in GBMs.     

Integrative subtype discovery in glioblastoma using iCluster

Large-scale cancer genome projects, such as the Cancer Genome Atlas (TCGA) project, are comprehensive molecular characterization efforts to accelerate our understanding of cancer biology and the discovery of new therapeutic targets. The accumulating wealth of multidimensional data provides a new paradigm for important research problems including cancer subtype discovery. The current standard approach relies on separate clustering analyses followed by manual integration. Results can be highly data type dependent, restricting the ability to discover new insights from multidimensional data. In this study, we present an integrative subtype analysis of the TCGA glioblastoma (GBM) data set. Our analysis revealed new insights through integrated subtype characterization. We found three distinct integrated tumor subtypes. Subtype 1 lacks the classical GBM events of chr 7 gain and chr 10 loss. This subclass is enriched for the G-CIMP phenotype and shows hypermethylation of genes involved in brain development and neuronal differentiation. The tumors in this subclass display a Proneural expression profile. Subtype 2 is characterized by a near complete association with EGFR amplification, overrepresentation of promoter methylation of homeobox and G-protein signaling genes, and a Classical expression profile. Subtype 3 is characterized by NF1 and PTEN alterations and exhibits a Mesenchymal-like expression profile. The data analysis workflow we propose provides a unified and computationally scalable framework to harness the full potential of large-scale integrated cancer genomic data for integrative subtype discovery.     

A Fourteen Gene GBM Prognostic Signature Identifies Association of Immune Response Pathway and Mesenchymal Subtype with High Risk Group

Background

Recent research on glioblastoma (GBM) has focused on deducing gene signatures predicting prognosis. The present study evaluated the mRNA expression of selected genes and correlated with outcome to arrive at a prognostic gene signature.

Methods

Patients with GBM (n = 123) were prospectively recruited, treated with a uniform protocol and followed up. Expression of 175 genes in GBM tissue was determined using qRT-PCR. A supervised principal component analysis followed by derivation of gene signature was performed. Independent validation of the signature was done using TCGA data. Gene Ontology and KEGG pathway analysis was carried out among patients from TCGA cohort.

Results

A 14 gene signature was identified that predicted outcome in GBM. A weighted gene (WG) score was found to be an independent predictor of survival in multivariate analysis in the present cohort (HR = 2^.507; B = 0^.919; p<0^.001) and in TCGA cohort. Risk stratification by standardized WG score classified patients into low and high risk predicting survival both in our cohort (p = <0^.001) and TCGA cohort (p = 0^.001). Pathway analysis using the most differentially regulated genes (n = 76) between the low and high risk groups revealed association of activated inflammatory/immune response pathways and mesenchymal subtype in the high risk group.

Conclusion

We have identified a 14 gene expression signature that can predict survival in GBM patients. A network analysis revealed activation of inflammatory response pathway specifically in high risk group. These findings may have implications in understanding of gliomagenesis, development of targeted therapies and selection of high risk cancer patients for alternate adjuvant therapies.     

Increased microglia/macrophage gene expression in a subset of adult and pediatric astrocytomas

Glioblastoma (GBM) is a highly malignant brain tumor with a dismal prognosis. Gene expression profiling of GBM has revealed clinically relevant tumor subtypes, and this provides exciting opportunities to better understand disease pathogenesis. Results from an increasing number of studies demonstrate a role for the immune response in cancer progression, yet it is unclear how the immune response differs across tumor subtypes and how it affects outcome. Utilizing gene expression data from The Cancer Genome Atlas Project and the Gene Expression Omnibus database, we demonstrate an enrichment of immune response-related gene expression in the mesenchymal subtype of adult GBM (n = 173) and pediatric high-grade gliomas (n = 53). In an independent cohort of pediatric astrocytomas (n = 24) from UCSF, we stratified tumors into subtypes and confirmed these findings. Using novel immune cell-specific gene signatures we demonstrate selective enrichment of microglia/macrophage-related genes in adult and pediatric GBM tumors of the mesenchymal subtype. Furthermore, immunostaining of adult GBM tumors showed significantly higher cell numbers of microglia/macrophages in mesenchymal versus non-mesenchymal tumors (p = 0.04). Interestingly, adult GBM tumors with the shortest survival had significant enrichment of microglia/macrophage-related genes but this was not true for pediatric GBMs. Consistent with an association with poor outcome, immune response-related genes were highly represented in an adult poor prognosis gene signature, with the expression of genes related to macrophage recruitment and activation being most strongly associated with survival (p<0.05) using CoxBoost multivariate modeling. Using a microglia/macrophage high gene signature derived from quantification of tumor-infiltrating cells in adult GBM, we identified enrichment of genes characteristic of CD4 T cells, granulocytes, and microglia/macrophages (n = 573). These studies support a role for the immune response, particularly the microglia/macrophage response, in the biology of an important subset of GBM. Identification of this subset may be important for future therapeutic stratification.     

Network Signatures of Survival in Glioblastoma Multiforme

To determine a molecular basis for prognostic differences in glioblastoma multiforme (GBM), we employed a combinatorial network analysis framework to exhaustively search for molecular patterns in protein-protein interaction (PPI) networks. We identified a dysregulated molecular signature distinguishing short-term (survival<225 days) from long-term (survival>635 days) survivors of GBM using whole genome expression data from The Cancer Genome Atlas (TCGA). A 50-gene subnetwork signature achieved 80% prediction accuracy when tested against an independent gene expression dataset. Functional annotations for the subnetwork signature included “protein kinase cascade,” “IκB kinase/NFκB cascade,” and “regulation of programmed cell death” – all of which were not significant in signatures of existing subtypes. Finally, we used label-free proteomics to examine how our subnetwork signature predicted protein level expression differences in an independent GBM cohort of 16 patients. We found that the genes discovered using network biology had a higher probability of dysregulated protein expression than either genes exhibiting individual differential expression or genes derived from known GBM subtypes. In particular, the long-term survivor subtype was characterized by increased protein expression of DNM1 and MAPK1 and decreased expression of HSPA9, PSMD3, and CANX. Overall, we demonstrate that the combinatorial analysis of gene expression data constrained by PPIs outlines an approach for the discovery of robust and translatable molecular signatures in GBM.     

The Proneural Molecular Signature Is Enriched in Oligodendrogliomas and Predicts Improved Survival among Diffuse Gliomas

The Cancer Genome Atlas Project (TCGA) has produced an extensive collection of ‘-omic’ data on glioblastoma (GBM), resulting in several key insights on expression signatures. Despite the richness of TCGA GBM data, the absence of lower grade gliomas in this data set prevents analysis genes related to progression and the uncovering of predictive signatures. A complementary dataset exists in the form of the NCI Repository for Molecular Brain Neoplasia Data (Rembrandt), which contains molecular and clinical data for diffuse gliomas across the full spectrum of histologic class and grade. Here we present an investigation of the significance of the TCGA consortium''s expression classification when applied to Rembrandt gliomas. We demonstrate that the proneural signature predicts improved clinical outcome among 176 Rembrandt gliomas that includes all histologies and grades, including GBMs (log rank test p = 1.16e-6), but also among 75 grade II and grade III samples (p = 2.65e-4). This gene expression signature was enriched in tumors with oligodendroglioma histology and also predicted improved survival in this tumor type (n = 43, p = 1.25e-4). Thus, expression signatures identified in the TCGA analysis of GBMs also have intrinsic prognostic value for lower grade oligodendrogliomas, and likely represent important differences in tumor biology with implications for treatment and therapy. Integrated DNA and RNA analysis of low-grade and high-grade proneural gliomas identified increased expression and gene amplification of several genes including GLIS3, TGFB2, TNC, AURKA, and VEGFA in proneural GBMs, with corresponding loss of DLL3 and HEY2. Pathway analysis highlights the importance of the Notch and Hedgehog pathways in the proneural subtype. This demonstrates that the expression signatures identified in the TCGA analysis of GBMs also have intrinsic prognostic value for low-grade oligodendrogliomas, and likely represent important differences in tumor biology with implications for treatment and therapy.     

Novel gene signatures for prognosis prediction in ovarian cancer

Ovarian cancer (OV) is one of the leading causes of cancer deaths in women worldwide. Late diagnosis and heterogeneous treatment result to poor survival outcomes for patients with OV. Therefore, we aimed to develop novel biomarkers for prognosis prediction from the potential molecular mechanism of tumorigenesis. Eight eligible data sets related to OV in GEO database were integrated to identify differential expression genes (DEGs) between tumour tissues and normal. Enrichment analyses discovered DEGs were most significantly enriched in G2/M checkpoint signalling pathway. Subsequently, we constructed a multi‐gene signature based on the LASSO Cox regression model in the TCGA database and time‐dependent ROC curves showed good predictive accuracy for 1‐, 3‐ and 5‐year overall survival. Utility in various types of OV was validated through subgroup survival analysis. Risk scores formulated by the multi‐gene signature stratified patients into high‐risk and low‐risk, and the former inclined worse overall survival than the latter. By incorporating this signature with age and pathological tumour stage, a visual predictive nomogram was established, which was useful for clinicians to predict survival outcome of patients. Furthermore, SNRPD1 and EFNA5 were selected from the multi‐gene signature as simplified prognostic indicators. Higher EFNA5 expression or lower SNRPD1 indicated poorer outcome. The correlation between signature gene expression and clinical characteristics was observed through WGCNA. Drug‐gene interaction was used to identify 16 potentially targeted drugs for OV treatment. In conclusion, we established novel gene signatures as independent prognostic factors to stratify the risk of OV patients and facilitate the implementation of personalized therapies.     

Glioblastoma models reveal the connection between adult glial progenitors and the proneural phenotype

Background

Tumor heterogeneity is a major obstacle for finding effective treatment of Glioblastoma (GBM). Based on global expression analysis, GBM can be classified into distinct subtypes: Proneural, Neural, Classical and Mesenchymal. The signatures of these different tumor subtypes may reflect the phenotypes of cells giving rise to them. However, the experimental evidence connecting any specific subtype of GBM to particular cells of origin is lacking. In addition, it is unclear how different genetic alterations interact with cells of origin in determining tumor heterogeneity. This issue cannot be addressed by studying end-stage human tumors.

Methodology/Principal Findings

To address this issue, we used retroviruses to deliver transforming genetic lesions to glial progenitors in adult mouse brain. We compared the resulting tumors to human GBM. We found that different initiating genetic lesions gave rise to tumors with different growth rates. However all mouse tumors closely resembled the human Proneural GBM. Comparative analysis of these mouse tumors allowed us to identify a set of genes whose expression in humans with Proneural GBM correlates with survival.

Conclusions/Significance

This study offers insights into the relationship between adult glial progenitors and Proneural GBM, and allows us to identify molecular alterations that lead to more aggressive tumor growth. In addition, we present a new preclinical model that can be used to test treatments directed at a specific type of GBM in future studies.     

Radiogenomic mapping of edema/cellular invasion MRI-phenotypes in glioblastoma multiforme

Background

Despite recent discoveries of new molecular targets and pathways, the search for an effective therapy for Glioblastoma Multiforme (GBM) continues. A newly emerged field, radiogenomics, links gene expression profiles with MRI phenotypes. MRI-FLAIR is a noninvasive diagnostic modality and was previously found to correlate with cellular invasion in GBM. Thus, our radiogenomic screen has the potential to reveal novel molecular determinants of invasion. Here, we present the first comprehensive radiogenomic analysis using quantitative MRI volumetrics and large-scale gene- and microRNA expression profiling in GBM.

Methods

Based on The Cancer Genome Atlas (TCGA), discovery and validation sets with gene, microRNA, and quantitative MR-imaging data were created. Top concordant genes and microRNAs correlated with high FLAIR volumes from both sets were further characterized by Kaplan Meier survival statistics, microRNA-gene correlation analyses, and GBM molecular subtype-specific distribution.

Results

The top upregulated gene in both the discovery (4 fold) and validation (11 fold) sets was PERIOSTIN (POSTN). The top downregulated microRNA in both sets was miR-219, which is predicted to bind to POSTN. Kaplan Meier analysis demonstrated that above median expression of POSTN resulted in significantly decreased survival and shorter time to disease progression (P<0.001). High POSTN and low miR-219 expression were significantly associated with the mesenchymal GBM subtype (P<0.0001).

Conclusion

Here, we propose a novel diagnostic method to screen for molecular cancer subtypes and genomic correlates of cellular invasion. Our findings also have potential therapeutic significance since successful molecular inhibition of invasion will improve therapy and patient survival in GBM.     

A ten-microRNA expression signature predicts survival in glioblastoma

Glioblastoma (GBM) is the most common and aggressive primary brain tumor with very poor patient median survival. To identify a microRNA (miRNA) expression signature that can predict GBM patient survival, we analyzed the miRNA expression data of GBM patients (n=222) derived from The Cancer Genome Atlas (TCGA) dataset. We divided the patients randomly into training and testing sets with equal number in each group. We identified 10 significant miRNAs using Cox regression analysis on the training set and formulated a risk score based on the expression signature of these miRNAs that segregated the patients into high and low risk groups with significantly different survival times (hazard ratio [HR]=2.4; 95% CI=1.4-3.8; p<0.0001). Of these 10 miRNAs, 7 were found to be risky miRNAs and 3 were found to be protective. This signature was independently validated in the testing set (HR=1.7; 95% CI=1.1-2.8; p=0.002). GBM patients with high risk scores had overall poor survival compared to the patients with low risk scores. Overall survival among the entire patient set was 35.0% at 2 years, 21.5% at 3 years, 18.5% at 4 years and 11.8% at 5 years in the low risk group, versus 11.0%, 5.5%, 0.0 and 0.0% respectively in the high risk group (HR=2.0; 95% CI=1.4-2.8; p<0.0001). Cox multivariate analysis with patient age as a covariate on the entire patient set identified risk score based on the 10 miRNA expression signature to be an independent predictor of patient survival (HR=1.120; 95% CI=1.04-1.20; p=0.003). Thus we have identified a miRNA expression signature that can predict GBM patient survival. These findings may have implications in the understanding of gliomagenesis, development of targeted therapy and selection of high risk cancer patients for adjuvant therapy.     

Comparative analyses identify molecular signature of MRI-classified SVZ-associated glioblastoma

Glioblastoma (GBM) is a highly aggressive brain cancer with limited therapeutic options. While efforts to identify genes responsible for GBM have revealed mutations and aberrant gene expression associated with distinct types of GBM, patients with GBM are often diagnosed and classified based on MRI features. Therefore, we seek to identify molecular representatives in parallel with MRI classification for group I and group II primary GBM associated with the subventricular zone (SVZ). As group I and II GBM contain stem-like signature, we compared gene expression profiles between these 2 groups of primary GBM and endogenous neural stem progenitor cells to reveal dysregulation of cell cycle, chromatin status, cellular morphogenesis, and signaling pathways in these 2 types of MRI-classified GBM. In the absence of IDH mutation, several genes associated with metabolism are differentially expressed in these subtypes of primary GBM, implicating metabolic reprogramming occurs in tumor microenvironment. Furthermore, histone lysine methyltransferase EZH2 was upregulated while histone lysine demethylases KDM2 and KDM4 were downregulated in both group I and II primary GBM. Lastly, we identified 9 common genes across large data sets of gene expression profiles among MRI-classified group I/II GBM, a large cohort of GBM subtypes from TCGA, and glioma stem cells by unsupervised clustering comparison. These commonly upregulated genes have known functions in cell cycle, centromere assembly, chromosome segregation, and mitotic progression. Our findings highlight altered expression of genes important in chromosome integrity across all GBM, suggesting a common mechanism of disrupted fidelity of chromosome structure in GBM.     

急性髓系白血病铁死亡相关基因筛选及预后模型的建立

基于急性髓系白血病(Acute Myeloid Leukemia,AML)临床大数据及多组学数据库探讨铁死亡相关基因在AML中的作用,并建立铁死亡基因表达相关预后模型。整合TCGA数据库中151例AML患者和GTEx数据库中337例正常人外周血的临床和转录组数据。将Wilcoxon检验和单因素Cox分析结果取交集,筛选出预后相关差异表达基因(Differential Expression Genes, DEGs),使用Lasso回归建立基因标志物预后模型,利用受试者工作特征曲线(Receiver Operating Characteristic Curve,ROC曲线)评价预测价值,Kaplan-Meier法进行生存分析,对AML患者临床数据进行单因素和多因素Cox回归分析,使用差异基因表达分析等方法比较高、低风险患者间的组学差异,最后,利用BeatAML数据库对基因标志物进行验证。将差异基因表达分析和单因素分析结果取交集,得到13个预后相关DEGs。构建了8个基因标志物的预后评分模型,并将患者分为高、低风险两组;ROC曲线分析证实了模型良好的预测性能;生存分析提示高、低风险组患者的生存率具有显著差异;单因素分析显示年龄和风险评分与患者整体生存显著相关,多因素分析显示,年龄和风险评分是独立预后指标。在2个风险组之间筛选出384个DEGs,GO富集分析结果显示,富集的基因大多与中性粒细胞和白细胞的趋化与迁移等免疫相关分子和通路显著相关,KEGG富集通路主要与TNF信号通路、细胞因子与细胞因子受体相互作用相关。BeatAML数据库验证结果显示,5个基因与预后显著相关。铁死亡相关基因在AML中显著表达,且高风险患者预后较差,该研究对AML铁死亡相关潜在生物标志物的发现和应用奠定了一定的基础。