Diverse genetic, epigenetic, and developmental programs drive glioblastoma, an incurable and poorly understood tumor, but their precise characterization remains challenging. Here, we use an integrative approach spanning single-cell RNA-sequencing of 28 tumors, bulk genetic and expression analysis of 401 specimens from the The Cancer Genome Atlas (TCGA), functional approaches, and single-cell lineage tracing to derive a unified model of cellular states and genetic diversity in glioblastoma. We find that malignant cells in glioblastoma exist in four main cellular states that recapitulate distinct neural cell types, are influenced by the tumor microenvironment, and exhibit plasticity. The relative frequency of cells in each state varies between glioblastoma samples and is influenced by copy number amplifications of the CDK4, EGFR, and PDGFRA loci and by mutations in the NF1 locus, which each favor a defined state. Our work provides a blueprint for glioblastoma, integrating the malignant cell programs, their plasticity, and their modulation by genetic drivers.

Although human tumours are shaped by the genetic evolution of cancer cells, evidence also suggests that they display hierarchies related to developmental pathways and epigenetic programs in which cancer stem cells (CSCs) can drive tumour growth and give rise to differentiated progeny1. Yet, unbiased evidence for CSCs in solid human malignancies remains elusive. Here we profile 4,347 single cells from six IDH1 or IDH2 mutant human oligodendrogliomas by RNA sequencing (RNA-seq) and reconstruct their developmental programs from genome-wide expression signatures. We infer that most cancer cells are differentiated along two specialized glial programs, whereas a rare subpopulation of cells is undifferentiated and associated with a neural stem cell expression program. Cells with expression signatures for proliferation are highly enriched in this rare subpopulation, consistent with a model in which CSCs are primarily responsible for fuelling the growth of oligodendroglioma in humans. Analysis of copy number variation (CNV) shows that distinct CNV sub-clones within tumours display similar cellular hierarchies, suggesting that the architecture of oligodendroglioma is primarily dictated by developmental programs. Subclonal point mutation analysis supports a similar model, although a full phylogenetic tree would be required to definitively determine the effect of genetic evolution on the inferred hierarchies. Our single-cell analyses provide insight into the cellular architecture of oligodendrogliomas at single-cell resolution and support the cancer stem cell model, with substantial implications for disease management.

Summary Developmental fate decisions are dictated by master transcription factors (TFs) that interact with cis-regulatory elements to direct transcriptional programs. Certain malignant tumors may also depend on cellular hierarchies reminiscent of normal development but superimposed on underlying genetic aberrations. In glioblastoma (GBM), a subset of stem-like tumor-propagating cells (TPCs) appears to drive tumor progression and underlie therapeutic resistance, yet remain poorly understood. Here, we identify a core set of neurodevelopmental TFs (POU3F2, SOX2, SALL2, OLIG2) essential for GBM propagation. These TFs coordinately bind and activate TPC-specific regulatory elements, and are sufficient to fully reprogram differentiated GBM cells to ‘induced’ TPCs, recapitulating the epigenetic landscape and phenotype of native TPCs. We reconstruct a network model that highlights critical interactions and identifies novel therapeutic targets for eliminating TPCs. Our study establishes the epigenetic basis of a developmental hierarchy in GBM, provides detailed insight into underlying gene regulatory programs, and suggests attendant therapeutic strategies.

Summary The aberrant transcription factor EWS-FLI1 drives Ewing sarcoma yet its molecular function is incompletely understood. We find that EWS-FLI1 reprograms gene regulatory circuits in Ewing sarcoma by directly inducing or repressing enhancers. At GGAA repeat elements, which lack evolutionary conservation and regulatory potential in other cell types, EWS-FLI1 multimers induce chromatin opening and create de novo enhancers that physically interact with target promoters. Conversely, EWS-FLI1 inactivates conserved enhancers containing canonical ETS motifs by displacing wild type ETS transcription factors. These divergent chromatin-remodeling patterns repress tumor suppressors and mesenchymal lineage regulators, while activating oncogenes and new potential therapeutic targets, such as the kinase VRK1. Our findings demonstrate how EWS-FLI1 establishes an oncogenic regulatory program governing both tumor survival and differentiation.

The demonstration of induced pluripotency and direct lineage conversion has led to remarkable insights regarding the roles of transcription factors and chromatin regulators in mediating cell state transitions. Beyond its considerable implications for regenerative medicine, this body of work is highly relevant to multiple stages of oncogenesis, from the initial cellular transformation to the hierarchical organization of established malignancies. Here we review conceptual parallels between the respective biological phenomena, highlighting important inter-relationships among transcription factors, chromatin regulators and pre-existing epigenetic states. The shared mechanisms provide insights into oncogenic transformation, tumor heterogeneity and cancer stem cell models.

Alterations in transcriptional regulators can orchestrate oncogenic gene expression programs in cancer. Here, we show that the BRG1/BRM-associated factor (BAF) chromatin remodeling complex, which is mutated in over 20% of human tumors, interacts with EWSR1, a member of a family of proteins with prion-like domains (PrLD) that are frequent partners in oncogenic fusions with transcription factors. In Ewing sarcoma, we find that the BAF complex is recruited by the EWS-FLI1 fusion protein to tumor-specific enhancers and contributes to target gene activation. This process is a neomorphic property of EWS-FLI1 compared to wild-type FLI1 and depends on tyrosine residues that are necessary for phase transitions of the EWSR1 prion-like domain. Furthermore, fusion of short fragments of EWSR1 to FLI1 is sufficient to recapitulate BAF complex retargeting and EWS-FLI1 activities. Our studies thus demonstrate that the physical properties of prion-like domains can retarget critical chromatin regulatory complexes to establish and maintain oncogenic gene expression programs.

Overexpression of the polycomb group protein enhancer of zeste homologue 2 (EZH2) occurs in diverse malignancies, including prostate cancer, breast cancer, and glioblastoma multiforme (GBM). Based on its ability to modulate transcription of key genes implicated in cell cycle control, DNA repair, and cell differentiation, EZH2 is believed to play a crucial role in tissue-specific stem cell maintenance and tumor development. Here, we show that targeted pharmacologic disruption of EZH2 by the S-adenosylhomocysteine hydrolase inhibitor 3-deazaneplanocin A (DZNep), or its specific downregulation by short hairpin RNA (shRNA), strongly impairs GBM cancer stem cell (CSC) self-renewal in vitro and tumor-initiating capacity in vivo. Using genome-wide expression analysis of DZNep-treated GBM CSCs, we found the expression of c-myc, recently reported to be essential for GBM CSCs, to be strongly repressed upon EZH2 depletion. Specific shRNA-mediated downregulation of EZH2 in combination with chromatin immunoprecipitation experiments revealed that c-myc is a direct target of EZH2 in GBM CSCs. Taken together, our observations provide evidence that direct transcriptional regulation of c-myc by EZH2 may constitute a novel mechanism underlying GBM CSC maintenance and suggest that EZH2 may be a valuable new therapeutic target for GBM management.

Growth of numerous cancer types is believed to be driven by a subpopulation of poorly differentiated cells, often referred to as cancer stem cells (CSCs), that have the capacity for self-renewal, tumor initiation, and generation of nontumorigenic progeny. Despite their potentially key role in tumor establishment and maintenance, the energy requirements of these cells and the mechanisms that regulate their energy production are unknown. Here, we show that the oncofetal insulin-like growth factor 2 mRNA-binding protein 2 (IMP2, IGF2BP2) regulates oxidative phosphorylation (OXPHOS) in primary glioblastoma (GBM) sphere cultures (gliomaspheres), an established in vitro model for CSC expansion. We demonstrate that IMP2 binds several mRNAs that encode mitochondrial respiratory chain complex subunits and that it interacts with complex I (NADH:ubiquinone oxidoreductase) proteins. Depletion of IMP2 in gliomaspheres decreases their oxygen consumption rate and both complex I and complex IV activity that results in impaired clonogenicity in vitro and tumorigenicity in vivo. Importantly, inhibition of OXPHOS but not of glycolysis abolishes GBM cell clonogenicity. Our observations suggest that gliomaspheres depend on OXPHOS for their energy production and survival and that IMP2 expression provides a key mechanism to ensure OXPHOS maintenance by delivering respiratory chain subunit-encoding mRNAs to mitochondria and contributing to complex I and complex IV assembly.