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.

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 longstanding debate in gliomagenesis. In contrast to the similarity in gl...

Gliomas with histone H3 lysine27-to-methionine mutations (H3K27M-glioma) arise primarily in the midline of the central nervous system of young children, suggesting a cooperation between genetics and cellular context in tumorigenesis. Although the genetics of H3K27M-glioma are well characterized, their cellular architecture remains uncharted. We performed single-cell RNA sequencing in 3321 cells from six primary H3K27M-glioma and matched models. We found that H3K27M-glioma primarily contain cells that resemble oligodendrocyte precursor cells (OPC-like), whereas more differentiated malignant cells are a minority. OPC-like cells exhibit greater proliferation and tumor-propagating potential than their more differentiated counterparts and are at least in part sustained by signaling. Our study characterizes oncogenic and developmental programs in H3K27M-glioma at single-cell resolution and across genetic subclones, suggesting potential therapeutic targets in this disease.

Summary Glioblastoma, the most common and aggressive malignant brain tumor, is propagated by stem-like cancer cells refractory to existing therapies. Understanding the molecular mechanisms that control glioblastoma stem cell (GSC) proliferation and drug resistance may reveal opportunities for therapeutic interventions. Here we show GSCs can reversibly transition to a slow-cycling, persistent state in response to targeted kinase inhibitors. In this state, GSCs upregulate primitive developmental programs and are dependent upon Notch signaling. This transition is accompanied by widespread redistribution of repressive histone methylation. Accordingly, persister GSCs upregulate, and are dependent on, the histone demethylases KDM6A/B. Importantly, slow-cycling cells with high Notch activity and histone demethylase expression are present in primary glioblastomas before treatment, potentially contributing to relapse. Our findings illustrate how cancer cells may hijack aspects of native developmental programs for deranged proliferation, adaptation, and tolerance. They also suggest strategies for eliminating refractory tumor cells by targeting epigenetic and developmental pathways.

No abstract

The work by Tang et al. provides a comprehensive, single-cell, transcriptomic analysis of kidney and blood cells from the adult zebrafish, identifying novel cell types, including two classes of NK immune cells, classically defined and erythroid-primed hematopoietic stem and progenitor cells, mucin-secreting kidney cells, and kidney stem/progenitor cells.

Melanoma-derived brain metastases (MBM) represent an unmet clinical need because central nervous system (CNS) progression is frequently an end-stage of the disease. Immune checkpoint inhibitors (ICI) provide a clinical opportunity against MBM; however, the MBM tumor microenvironment (TME) has not been fully elucidated in the context of ICI. To dissect unique elements of the MBM TME and correlates of MBM response to ICI, we collected 32 fresh MBM and performed single-cell RNA sequencing of the MBM TME and T-cell receptor clonotyping on T cells from MBM and matched blood and extracranial lesions. We observed myeloid phenotypic heterogeneity in the MBM TME, most notably multiple distinct neutrophil states, including an IL8-expressing population that correlated with malignant cell epithelial-to-mesenchymal transition. Additionally, we observed significant relationships between intracranial T-cell phenotypes and the distribution of T-cell clonotypes intracranially and peripherally. We found that the phenotype, clonotype, and overall number of MBM-infiltrating T cells were associated with response to ICI, suggesting that ICI-responsive MBMs interact with peripheral blood in a manner similar to extracranial lesions. These data identify unique features of the MBM TME that may represent potential targets to improve clinical outcomes for patients with MBM.