Mario L. Suvà

Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma

Cancer at single-cell resolution Single-cell sequencing can illuminate the genetic properties of brain cancers and reveal heterogeneity within a tumor. Patel et al. examined the genome sequence of single cells isolated from brain glioblastomas. The findings revealed shared chromosomal changes but also extensive transcription variation, including genes related to signaling, which represent potential therapeutic targets. The authors suggest that the variation in tumor cells reflects neural development and that such variation among cancer cells may prove to have clinical significance. Science, this issue p. 1396 Screening individual cancer cells within a brain tumor may help to guide treatment and predict prognosis. Human cancers are complex ecosystems composed of cells with distinct phenotypes, genotypes, and epigenetic states, but current models do not adequately reflect tumor composition in patients. We used single-cell RNA sequencing (RNA-seq) to profile 430 cells from five primary glioblastomas, which we found to be inherently variable in their expression of diverse transcriptional programs related to oncogenic signaling, proliferation, complement/immune response, and hypoxia. We also observed a continuum of stemness-related expression states that enabled us to identify putative regulators of stemness in vivo. Finally, we show that established glioblastoma subtype classifiers are variably expressed across individual cells within a tumor and demonstrate the potential prognostic implications of such intratumoral heterogeneity. Thus, we reveal previously unappreciated heterogeneity in diverse regulatory programs central to glioblastoma biology, prognosis, and therapy.

Epigenetic Reprogramming in Cancer

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 interrelationships among transcription factors, chromatin regulators, and preexisting epigenetic states. The shared mechanisms provide insights into oncogenic transformation, tumor heterogeneity, and cancer stem cell models.

Insulator dysfunction and oncogene activation in IDH mutant gliomas

Gain-of-function IDH mutations are initiating events that define major clinical and prognostic classes of gliomas. Mutant IDH protein produces a new onco-metabolite, 2-hydroxyglutarate, which interferes with iron-dependent hydroxylases, including the TET family of 5′-methylcytosine hydroxylases. TET enzymes catalyse a key step in the removal of DNA methylation. IDH mutant gliomas thus manifest a CpG island methylator phenotype (G-CIMP), although the functional importance of this altered epigenetic state remains unclear. Here we show that human IDH mutant gliomas exhibit hypermethylation at cohesin and CCCTC-binding factor (CTCF)-binding sites, compromising binding of this methylation-sensitive insulator protein. Reduced CTCF binding is associated with loss of insulation between topological domains and aberrant gene activation. We specifically demonstrate that loss of CTCF at a domain boundary permits a constitutive enhancer to interact aberrantly with the receptor tyrosine kinase gene PDGFRA, a prominent glioma oncogene. Treatment of IDH mutant gliomaspheres with a demethylating agent partially restores insulator function and downregulates PDGFRA. Conversely, CRISPR-mediated disruption of the CTCF motif in IDH wild-type gliomaspheres upregulates PDGFRA and increases proliferation. Our study suggests that IDH mutations promote gliomagenesis by disrupting chromosomal topology and allowing aberrant regulatory interactions that induce oncogene expression.

Alternative Lengthening of Telomeres Renders Cancer Cells Hypersensitive to ATR Inhibitors

Cancers alternative means to an end To stay alive and proliferating, tumor cells must maintain their telomeres: the DNA sequences at the ends of chromosomes. The majority accomplish this by activating the enzyme telomerase. However, certain tumor types favor a different mechanism called alternative lengthening of telomeres (ALT), which involves DNA recombination. Flynn et al. delineated the molecular events that occur at the telomeres of ALT-proficient tumor cells by studying the function of a protein that is altered by mutation in these tumors. The analysis revealed a specific protein kinase that is essential for ALT, which could in principle be targeted to halt tumor growth. Science, this issue p. 273 A potential therapeutic strategy is identified for tumor cells that maintain their telomeres by an unusual mechanism. Cancer cells rely on telomerase or the alternative lengthening of telomeres (ALT) pathway to overcome replicative mortality. ALT is mediated by recombination and is prevalent in a subset of human cancers, yet whether it can be exploited therapeutically remains unknown. Loss of the chromatin-remodeling protein ATRX associates with ALT in cancers. Here, we show that ATRX loss compromises cell-cycle regulation of the telomeric noncoding RNA TERRA and leads to persistent association of replication protein A (RPA) with telomeres after DNA replication, creating a recombinogenic nucleoprotein structure. Inhibition of the protein kinase ATR, a critical regulator of recombination recruited by RPA, disrupts ALT and triggers chromosome fragmentation and apoptosis in ALT cells. The cell death induced by ATR inhibitors is highly selective for cancer cells that rely on ALT, suggesting that such inhibitors may be useful for treatment of ALT-positive cancers.

Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma

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 progeny. 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.

Imp2 controls oxidative phosphorylation and is crucial for preserving glioblastoma cancer stem cells

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.

EWS-FLI1 Utilizes Divergent Chromatin Remodeling Mechanisms to Directly Activate or Repress Enhancer Elements in Ewing Sarcoma

The aberrant transcription factor EWS-FLI1 drives Ewing sarcoma, but its molecular function is not completely 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 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.

Decoupling genetics, lineages, and microenvironment in IDH-mutant gliomas by single-cell RNA-seq.

Single-cell RNA sequencing identifies a common origin for specific types of human glioma brain tumors. Effects of the tumor microenvironment Glioma brain tumors that carry mutant copies of the IDH gene can be subdivided into two major classes. However, the development of and differences between these two classes are not well characterized. Venteicher et al. coupled bulk sequencing and publicly available data with single-cell RNA sequencing data on glioma patient tissue samples. They identified a common lineage program that is shared between glioma subtypes. This suggests that the observed differences between the two glioma classes originate from lineage-specific genetic changes and the tumor microenvironment. Science, this issue p. eaai8478 INTRODUCTION Tumor fitness, evolution, and resistance to therapy are governed by selection of malignant cells with specific genotypes, by expression programs related to cellular phenotypes, and by influences of the tumor microenvironment (TME). Although bulk tumor analysis can interrogate the genetic state of tumor cells with high precision, bulk expression profiles average the diverse cells within each tumor, thereby masking critical differences and providing limited insight into cancer cell programs and TME influences. Single-cell RNA sequencing (scRNA-seq) can help to address those challenges but incurs financial and logistic considerations, including the time required to accrue large cohorts of fresh tumor specimen for single-cell analysis. RATIONALE We reasoned that scRNA-seq of a limited number of representative tumors could be combined with bulk data from large cohorts to decipher differences between tumor subclasses. In this approach, bulk samples collected for large cohorts, such as from The Cancer Genome Atlas (TCGA), are first used to define the combined effects of differences in cancer cell genotypes, phenotypes, and the composition of the TME. Single-cell analysis of a limited set of representative tumors is then used to distinguish those effects. We applied this approach to understand the differences between two types of isocitrate dehydrogenase (IDH)–mutant gliomas: astrocytoma (IDH-A) and oligodendroglioma (IDH-O). IDH-A and IDH-O are distinguished by co-occurring signature genetic events and by histopathology and are thought to recapitulate distinct glial lineages. By combining 9879 scRNA-seq profiles from 10 IDH-A tumors, 4347 scRNA-seq profiles from six IDH-O tumors, and 165 TCGA bulk RNA profiles, we could decipher differences between these two tumor types at single-cell resolution. RESULTS We find that differences in bulk expression profiles between IDH-A and IDH-O are primarily explained by the impact of signature genetic events and TME composition, but not by distinct expression programs of glial lineages in the malignant cells. We infer that both IDH-A and IDH-O share the same developmental hierarchy, consisting in each case of three subpopulations of malignant cells: nonproliferating cells differentiated along the astrocytic and oligodendrocytic lineages, and proliferative undifferentiated cells that resemble neural stem/progenitor cells. By analyzing tumors of different clinical grades, we observe that higher-grade tumors present enhanced proliferation, larger pools of undifferentiated glioma cells, and an increase in macrophage over microglia programs in the TME. CONCLUSION Our approach provides a general framework to decipher differences between classes of human tumors by decoupling cancer cell genotypes, phenotypes, and the composition of the TME. The shared glial lineages and developmental hierarchies observed in IDH-A and IDH-O suggest a common progenitor for all IDH-mutant gliomas, shedding light on a long-standing debate in gliomagenesis. In contrast to the similarity in glial lineages, IDH-A and IDH-O differ significantly in their TME, and in particular in the abundance of microglia/macrophage cells. Microglia and macrophages also differ between IDH-A tumors of different grades. Our study redefines the cellular composition of human IDH-mutant gliomas, with important implications for disease management. Single-cell RNA-seq of IDH-mutant gliomas reveals tumor architecture. (Top) Human samples were dissociated and analyzed by scRNA-seq. (Bottom) IDH-O and IDH-A differ in genetics and TME but are both primarily composed of three main types of malignant cells: cycling stem-like cells and noncycling astrocyte-like and oligodendrocyte-like cells. Tumor progression is associated with increased proliferation, decreased differentiation, and increase in macrophages over microglia in the TME. Tumor subclasses differ according to the genotypes and phenotypes of malignant cells as well as the composition of the tumor microenvironment (TME). We dissected these influences in isocitrate dehydrogenase (IDH)–mutant gliomas by combining 14,226 single-cell RNA sequencing (RNA-seq) profiles from 16 patient samples with bulk RNA-seq profiles from 165 patient samples. Differences in bulk profiles between IDH-mutant astrocytoma and oligodendroglioma can be primarily explained by distinct TME and signature genetic events, whereas both tumor types share similar developmental hierarchies and lineages of glial differentiation. As tumor grade increases, we find enhanced proliferation of malignant cells, larger pools of undifferentiated glioma cells, and an increase in macrophage over microglia expression programs in TME. Our work provides a unifying model for IDH-mutant gliomas and a general framework for dissecting the differences among human tumor subclasses.

A Tell-Tail Sign of Chromatin: Histone Mutations Drive Pediatric Glioblastoma

Recent genomic analyses of pediatric glioblastoma, a poorly understood tumor with dismal outcome, have identified mutations in histone H3 variants that affect critical amino acids in the tail. The findings extend discoveries of chromatin regulator inactivation and gain-of-function mutations by documenting alteration of a modifiable histone residue in human cancer.

Multiplexed RNAi therapy against brain tumor-initiating cells via lipopolymeric nanoparticle infusion delays glioblastoma progression

Significance Glioblastoma is a deadly brain tumor with no cure. Brain tumor-initiating cells (BTICs) have been recognized as the key driver behind the unstoppable malignant growth, therapy resistance, and recurrence. BTICs are exceptionally difficult to target because of heterogeneous genetic and epigenetic aberrations that are challenging to reverse therapeutically using conventional pharmaceuticals or biologics. Here we report a lipopolymeric nanoparticle (LPNP) formulation that demonstrates a surprisingly high affinity for BTICs and the capacity to encapsulate multiple siRNAs for potent and targeted anti-BTIC therapy. We show that direct infusion of LPNP siRNAs to brain tumors effectively impedes tumor growth in mouse and provides encouraging survival benefits. This multiplexed nanomedicine platform carries strong potential for personalized anti-BTIC therapies. Brain tumor-initiating cells (BTICs) have been identified as key contributors to therapy resistance, recurrence, and progression of diffuse gliomas, particularly glioblastoma (GBM). BTICs are elusive therapeutic targets that reside across the blood–brain barrier, underscoring the urgent need to develop novel therapeutic strategies. Additionally, intratumoral heterogeneity and adaptations to therapeutic pressure by BTICs impede the discovery of effective anti-BTIC therapies and limit the efficacy of individual gene targeting. Recent discoveries in the genetic and epigenetic determinants of BTIC tumorigenesis offer novel opportunities for RNAi-mediated targeting of BTICs. Here we show that BTIC growth arrest in vitro and in vivo is accomplished via concurrent siRNA knockdown of four transcription factors (SOX2, OLIG2, SALL2, and POU3F2) that drive the proneural BTIC phenotype delivered by multiplexed siRNA encapsulation in the lipopolymeric nanoparticle 7C1. Importantly, we demonstrate that 7C1 nano-encapsulation of multiplexed RNAi is a viable BTIC-targeting strategy when delivered directly in vivo in an established mouse brain tumor. Therapeutic potential was most evident via a convection-enhanced delivery method, which shows significant extension of median survival in two patient-derived BTIC xenograft mouse models of GBM. Our study suggests that there is potential advantage in multiplexed targeting strategies for BTICs and establishes a flexible nonviral gene therapy platform with the capacity to channel multiplexed RNAi schemes to address the challenges posed by tumor heterogeneity.