Gene expression profiling distinguishes proneural glioma stem cells from mesenchymal glioma stem cells

Abstract

Tumor heterogeneity of high-grade glioma (HGG) is recognized by four clinically relevant subtypes based on core gene signatures. However, molecular signaling in glioma stem cells (GSCs) in individual HGG subtypes is poorly characterized. Previously we identified and characterized two mutually exclusive GSC subtypes with distinct activated signaling pathways and biological phenotypes. One GSC subtype presented with a gene signature resembling Proneural (PN) HGG, whereas the other was similar to mesenchymal (Mes) HGG. Classical HGG-derived GSCs were sub-classified as either one of these two subtypes. Differential mRNA expression analysis of PN and Mes GSCs identified 5796 differentially expressed genes, revealing a pronounced correlation with the corresponding PN or Mes HGGs. Mes GSCs displayed more aggressive phenotypes in vitro and as intracranial xenografts in mice. Further, Mes GSCs were markedly resistant to radiation compared with PN GSCs. Expression of ALDH1A3 — one of the most up-regulated Mes representative genes and a universal cancer stem cell marker in non-brain cancers — was associated with self-renewal and a multi-potent stem cell population in Mes but not PN samples. Moreover, inhibition of ALDH1A3 attenuated the growth of Mes but not PN GSCs in vitro. Lastly, radiation treatment of PN GSCs up-regulated Mes-associated markers and down-regulated PN-associated markers, whereas inhibition of ALDH1A3 attenuated an irradiation-induced gain of Mes identity in PN GSCs in vitro. Taken together, our data suggest that two subtypes of GSCs, harboring distinct metabolic signaling pathways, represent intertumoral glioma heterogeneity and highlight previously unidentified roles of ALDH1A3-associated signaling that promotes aberrant proliferation of Mes HGGs and GSCs. Inhibition of ALDH1A3-mediated pathways therefore might provide a promising therapeutic approach for a subset of HGGs with the Mes signature. Here, we describe the gene expression analysis, including pre-processing methods for the data published by Mao and colleagues in PNAS [1], integration of microarray data from this study with The Cancer Genome Atlas (TCGA) glioblastoma data and also with another published study. The raw CEL files and processed data were submitted to Gene Expression Omnibus (GEO) under the accession GSE67089.