Search

Diffuse midline gliomas driven by lysine27-to-methionine mutations in histone 3 (H3-K27M DMGs) are among the most fatal brain tumors. Molecular studies including single cell RNA-sequencing (scRNA-seq) of pediatric and predominantly pontine H3-K27M DMGs have shown that the H3-K27M oncohistone keeps glioma cells locked in a stem-like oligodendrocyte precursor cell (OPC) state that is capable of self-renewal and tumor-initiation. However, a comprehensive dissection of the cellular architecture of H3-K27M DMGs across different midline regions and age groups is required to better understand the cell-intrinsic and contextual regulation of H3-K27M DMG cell identities. In particular, the more recently described group of adult H3-K27M DMGs remains understudied. Here, we have collected and characterized 45 H3-K27M mutant patient tumors, spanning pontine (n=26), thalamic (n=17), and spinal (n=2) locations. Median age at surgery was 12 (2-68) years, encompassing 21 early childhood (0-10 years), 12 adolescent (11-20 years), and 12 adult (≥ 21 years) tumors. The majority of samples were obtained pre-treatment (n=28), as opposed to post-treatment or at autopsy (n=17). We profiled all 45 tumors by single cell/single nucleus RNA-seq and selected tumors were further characterized by the single cell assay for transposase-accessible chromatin (scATAC-seq). Our integrated analyses highlight the predominance of transcriptionally and epigenetically defined OPC-like tumor cells as the main cell population of H3-K27M DMGs across all age groups and locations. We further identify distinct age- and location-specific OPC-like cell subpopulations. Comparison of pediatric and adult tumors further demonstrates a significant increase of mesenchymal cell states in adult H3-K27M DMGs, which we link to differences in glioma-associated immune cell compartments between age groups. Together, this study sheds light on the effects of age- and region-dependent microenvironments in shaping cellular identities in H3-K27M DMGs.

Melanoma-derived brain metastases (MBM) represent an unmet clinical need due to central nervous system (CNS) progression as a frequent, end-stage site of disease. Immune checkpoint inhibition (ICI) represents a clinical opportunity against MBM; however, the MBM tumor microenvironment (TME) has not been fully elucidated in the context of ICI. To dissect unique MBM-TME elements and correlates of MBM-ICI response, 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, most notably multiple distinct neutrophil states including an IL-8 expressing population that correlated with malignant cell epithelial-to-mesenchymal transition. Additionally, we observe 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 demonstrate unique features of the MBM-TME, which may represent potential targets to improve clinical outcomes for patients with MBM.

Embryonal tumor with multilayered rosettes (ETMR) is a malignant brain tumor that typically occurs in children under the age of three. Most patients die within two years of diagnosis, and more effective, targeted therapies are urgently needed. To better characterize the oncogenic mechanisms of key driver alterations and identify novel therapeutic targets, we studied the cellular heterogeneity of ETMR using single-cell RNA sequencing. Analyses conducted on >3,000 high-quality cells collected from ten primary and relapse specimens revealed a common cellular hierarchy across all tumors: A highly proliferative neural stem cell-like population (SOX2+) gives rise to intermediate progenitors (ASCL1+) and more differentiated neuron-like cells (STMN2/4+). These malignant populations closely match histological patterns of ETMR (i.e. rosettes, neuropil), as observed by immunofluorescence microscopy. Comparison to single-cell datasets from human embryos indicates resemblance to cell populations of the developing brain, but also reveals key ETMR-specific differences, including expression of the chromosome 19 miRNA cluster (C19MC, the presumed genetic driver of most ETMRs), which is restricted to the stem cell-like population. We next investigated if targeting C19MC is a viable strategy to disrupt the cellular hierarchy of ETMR. Silencing with antisense oligonucleotides shows pronounced reduction of cell line growth for a specific subset of the 46 members of C19MC. These miRNAs share seed sequences with evolutionary conserved miRNAs that have been shown to regulate pluripotency and self-renewal of embryonic stem cells. We hypothesize that select C19MC members play similar roles in ETMR and represent bona fide targets for therapeutic targeting using antisense technology.

Ependymoma represents a heterogeneous disease affecting the entire neuraxis. Extensive molecular profiling efforts have identified molecular ependymoma subgroups based on DNA methylation. However, the intratumoral heterogeneity and developmental origins of these groups are only partially understood, and effective treatments are still lacking for about 50% of patients with high-risk tumors. We interrogated the cellular architecture of ependymoma using single cell/single nucleus RNA-sequencing to analyze 24 tumor specimens across major molecular subgroups and anatomic locations. We additionally analyzed ten patient-derived ependymoma cell models and two patient-derived xenografts (PDXs). Interestingly, we identified an analogous cellular hierarchy across all ependymoma groups, originating from undifferentiated neural stem cell-like populations towards different degrees of impaired differentiation states comprising neuronal precursor-like, astro-glial-like, and ependymal-like tumor cells. While prognostically favorable ependymoma groups predominantly harbored differentiated cell populations, aggressive groups were enriched for undifferentiated subpopulations. Projection of transcriptomic signatures onto an independent bulk RNA-seq cohort stratified patient survival even within known molecular groups, thus refining the prognostic power of DNA methylation-based profiling. Furthermore, we identified novel potentially druggable targets such as IGF- and FGF-signaling within poorly prognostic transcriptional programs. Ependymoma-derived cell models/PDXs widely recapitulated the transcriptional programs identified within fresh tumors and are leveraged to validate identified target genes in functional follow-up analyses. Taken together, our analyses reveal a developmental hierarchy and transcriptomic context underlying the biologically and clinically distinct behavior of ependymoma groups. The newly characterized cellular states and underlying regulatory networks could serve as basis for future therapeutic target identification and reveal biomarkers for clinical trials.

Neurofibrillary tangles (NFTs), composed of hyperphosphorylated tau, are a key pathologic feature of Alzheimer's disease (AD). Tau phosphorylation is under the control of multiple kinases and phosphatases, including Fyn. Previously, our group found an association between 2 regulatory single nucleotide polymorphisms in the FYN gene with increased tau levels in the cerebrospinal fluid. In this study, we hypothesized that Fyn expression in the brain is influenced by AD status and genetic content. We found that Fyn protein, but not messenger RNA, levels were increased in AD patients compared to cognitively normal controls and are associated with regulatory region single nucleotide polymorphisms. In addition, the expression of the FYN 3'UTR can decrease expression in multiple cell lines, suggesting this regulatory region plays an important role in FYN expression. Taken together, these data suggest that FYN expression is regulated according to AD status and regulatory region haplotype, and genetic variants may be instrumental in the development of neurofibrillary tangles in AD and other tauopathies.

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., Recent advances in single-cell, transcriptomic profiling have provided unprecedented access to investigate cell heterogeneity during tissue and organ development. In this study, we used massively parallel, single-cell RNA sequencing to define cell heterogeneity within the zebrafish kidney marrow, constructing a comprehensive molecular atlas of definitive hematopoiesis and functionally distinct renal cells found in adult zebrafish. Because our method analyzed blood and kidney cells in an unbiased manner, our approach was useful in characterizing immune-cell deficiencies within DNA–protein kinase catalytic subunit (prkdc), interleukin-2 receptor γ a (il2rga), and double-homozygous–mutant fish, identifying blood cell losses in T, B, and natural killer cells within specific genetic mutants. Our analysis also uncovered novel cell types, including two classes of natural killer immune cells, classically defined and erythroid-primed hematopoietic stem and progenitor cells, mucin-secreting kidney cells, and kidney stem/progenitor cells. In total, our work provides the first, comprehensive, single-cell, transcriptomic analysis of kidney and marrow cells in the adult zebrafish.

Ependymoma represents a heterogeneous disease affecting the entire neuraxis. Extensive molecular profiling efforts have identified molecular ependymoma subgroups based on DNA methylation. However, the intratumoral heterogeneity and developmental origins of these groups are only partially understood, and effective treatments are still lacking for about 50% of patients with high-risk tumors. We interrogated the cellular architecture of ependymoma using single cell/nucleus RNA-sequencing to analyze 24 tumor specimens across major molecular subgroups and anatomic locations. We additionally analyzed ten patient-derived ependymoma cell models and two patient-derived xenografts (PDXs). Interestingly, we identified an analogous cellular hierarchy across all ependymoma groups, originating from undifferentiated neural stem cell-like populations towards different degrees of impaired differentiation states comprising neuronal precursor-like, astro-glial-like, and ependymal-like tumor cells. While prognostically favorable ependymoma groups predominantly harbored differentiated cell populations, aggressive groups were enriched for undifferentiated subpopulations. Projection of transcriptomic signatures onto an independent bulk RNA-seq cohort stratified patient survival even within known molecular groups, thus refining the prognostic power of DNA methylation-based profiling. Furthermore, we identified novel potentially druggable targets including IGF- and FGF-signaling within poorly prognostic transcriptional programs. Ependymoma-derived cell models/PDXs widely recapitulated the transcriptional programs identified within fresh tumors and are leveraged to validate identified target genes in functional follow-up analyses. Taken together, our analyses reveal a developmental hierarchy and transcriptomic context underlying the biologically and clinically distinct behavior of ependymoma groups. The newly characterized cellular states and underlying regulatory networks could serve as basis for future therapeutic target identification and reveal biomarkers for clinical trials.

Medulloblastoma is a malignant childhood cerebellar tumor comprised of distinct molecular subgroups. Whereas genomic characteristics of these subgroups are well defined, the extent to which cellular diversity underlies their divergent biology and clinical behaviour remains largely unexplored. We used single-cell transcriptomics to investigate intra- and inter-tumoral heterogeneity in twenty-five medulloblastomas spanning all molecular subgroups. WNT, SHH, and Group 3 tumors comprised subgroup-specific undifferentiated and differentiated neuronal-like malignant populations, whereas Group 4 tumors were exclusively comprised of differentiated neuronal-like neoplastic cells. SHH tumors closely resembled granule neurons of varying differentiation states that correlated with patient age. Group 3 and Group 4 tumors exhibited a developmental trajectory from primitive progenitor-like to more mature neuronal-like cells, whose relative proportions distinguished these subgroups. Cross-species transcriptomics defined distinct glutamatergic populations as putative cells-of-origin for SHH and Group 4 subtypes. Collectively, these data provide novel insights into the cellular and developmental states underlying subtype-specific medulloblastoma biology.

High-grade gliomas harboring H3 G34R/V mutations exclusively occur in the cerebral hemispheres of adolescents and young adults, suggesting a distinct neurodevelopmental origin. Combining multimodal bulk and single-cell genomics with unbiased genome-scale CRISPR/Cas9 approaches, we here describe a GABAergic interneuron progenitor lineage as the most likely context from which these H3 G34R/V mutations drive gliomagenesis, conferring unique and tumor-selective gene targets essential for glioma cell survival, as validated genetically and pharmacologically. Phenotypically, we demonstrate that while H3 G34R/V glioma cells harbor the neurotransmitter GABA, they are developmentally stalled, and do not induce the neuronal hyperexcitability described in other glioma subtypes. These findings offer a striking counter-example to the prevailing view of glioma origins in glial precursor cells, resulting in distinct cellular, microenvironmental, and therapeutic consequences.

High-grade gliomas (HGG) are among the most prevalent and fatal cancers in pediatric, adolescent, and young adult (AYA) patients. Especially understudied are older children and young adults, aged 16–39 years. Previously, we profiled primary pediatric HGGs through single-cell transcriptomics and identified the genetic, epigenetic and developmental programs that drive their malignant progression. However, the questions of how these programs compare to those in older HGG patients, what the mechanisms are by which these tumors ultimately evolve to drive recurrence and treatment resistance, and how distinct tumor cell subpopulations bidirectionally communicate with their microenvironment remain to be elucidated. Here, we use single-nucleus RNA sequencing to compare 11 paired, matched high-grade gliomas at diagnosis and recurrence and 15 additional H3K27M primary and recurrent DMG samples in pediatric and AYA patients. In all tumors, we find both undifferentiated and differentiated tumor cells recapitulating distinct glial lineages, as well as diverse microenvironmental cell populations. When longitudinally comparing this tumor architecture within matched pairs, we find substantial differences in transcriptional program expressions. Diagnostic samples include more differentiated, astrocyte-like tumor cells, while cells from recurrent samples more highly express ribosomal and heat-shock protein genes, suggesting tumor progression- and treatment-related shifts. Ongoing sequencing and analysis will allow for unprecedented insight into the evolutionary dynamics of pediatric and AYA high-grade gliomas as well as delineate differences in the biology of DMGs occurring in different age groups. This multi-institutional project was funded by the National Institute of Health.

Summary Medulloblastoma is a malignant childhood cerebellar tumour comprised of distinct molecular subgroups. Whereas genomic characteristics of these subgroups are well defined, the extent to which cellular diversity underlies their divergent biology and clinical behaviour remains largely unexplored. We used single-cell transcriptomics to investigate intra- and inter-tumoural heterogeneity in twenty-five medulloblastomas spanning all molecular subgroups. WNT, SHH, and Group 3 tumours comprised subgroup-specific undifferentiated and differentiated neuronal-like malignant populations, whereas Group 4 tumours were exclusively comprised of differentiated neuronal-like neoplastic cells. SHH tumours closely resembled granule neurons of varying differentiation states that correlated with patient age. Group 3 and Group 4 tumours exhibited a developmental trajectory from primitive progenitor-like to more mature neuronal-like cells, whose relative proportions distinguished these subgroups. Cross-species transcriptomics defined distinct glutamatergic populations as putative cells-of-origin for SHH and Group 4 subtypes. Collectively, these data provide novel insights into the cellular and developmental states underlying subtype-specific medulloblastoma biology.

The cellular composition of H3K27M gliomas Diffuse midline gliomas with histone H3 lysine27-to-methionine mutations (H3K27M-glioma) are an aggressive type of childhood cancer with few options for treatment. Filbin et al. used a single-cell sequencing approach to study the oncogenic programs, genetics, and cellular hierarchies of H3K27M-glioma. Tumors were mainly composed of cells resembling oligodendrocyte precursor cells, whereas differentiated malignant cells were a smaller fraction. In comparison with other gliomas, these cancers had distinct oncogenic programs and stem cell–like profiles that contributed to their stable tumor-propagating potential. The analysis also identified a lineage-specific marker that may be useful in developing therapies. Science , this issue p. 331

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., National Cancer Institute (U.S.) (Grant P30-CA14051)

Human cancers are composed of cells with heterogeneous genetic and epigenetic states, resulting in substantial phenotypic diversity. Medulloblastoma (MB), a clinically challenging, malignant childhood brain tumor, is no exception. A better understanding of its cellular heterogeneity is urgently needed in order to develop more efficacious therapeutic strategies that eliminate all tumor subclones underlying MB initiation, maintenance, progression, and relapse. We used single-cell RNA sequencing to assess intratumoral heterogeneity in 20 MB samples representing all four molecular subgroups. Fresh surgical biopsies were dissociated and flow sorted into single cells prior to full-length transcriptome sequencing using the Smart-seq2 protocol. Computational analyses conducted on >6,000 single cells identified cellular hierarchies ranging from proliferative, undifferentiated cell types to more differentiated neural- and astrocyte-like progeny. Transcriptional differences in undifferentiated cell populations between MB subgroups suggested distinct cellular origins, and informed oncogenic pathways believed to promote MB pathogenesis. Identification of transcriptional programs active in single cells aided the deconvolution of published bulk RNA-seq profiles, providing a deeper understanding of the heterogeneity characteristic of Group 3/4 subtypes. Lastly, inference of copy-number variations and probable driver-gene mutations from the single-cell data provided evidence for genetic subclones in a subset of patients. Overlaying transcriptional profiles on MB genetic subclones unveiled how oncogenic programs adapted during tumor evolution. These interim results provide unparalleled insights into the cellular heterogeneity of MB subgroups. Analyses aimed at further resolving molecular substructure in Group 3/4, verifying the developmental origins of MB subgroups, and implicating defined cell populations driving tumorigenesis are ongoing.

Schwannomas are found either sporadically or as part of a syndrome. Genetic alterations found in schwannomas include those in NF2 and SMARCB1 (implicated in schwannomatosis), and BRAF V600E. Reduced copy number of chromosome 22q can be the mechanism for loss of the unaltered allele ofNF2 and/orSMARCB1in some schwannomas. We analyzed a cohort of 57 vestibular and non-vestibular schwannomas from three clinical subtypes (Neurofibromatosis 2, schwannomatosis, and sporadic) for reduced copy number (NF2 and SMARCB1), and mutations (NF2, SMARCB1, LZTR1, NF1 and BRAF). Tumors were histologically classified as conventional/classic, or hybrid. Genomic DNA was extracted from FFPE biopsies. Full coding regions of genes NF1, NF2, LZTR1, SMARCB1 and BRAF were enriched by anchored multiplex-PCR (AMP). Resulting libraries were quantitated and sequenced on the Illumina MiSeq System. Bioinformatics analyses were performed using custom scripts. Copy number alteration was detected using a read depth approach after consolidation of single molecule derived reads. INI1 immunohistochemistry was performed. RNA ISH for SMARCB1 was performed. Five tumors were of hybrid histology, and all five were syndromic. Reduced copy number for NF2 was identified in 33% of vestibular schwannomas and 18% of non-vestibular schwannomas. There was no significant difference in reduced copy number of NF2 between sporadic and syndromic forms. BRAF V600E and NF1 mutations were not identified in any case. For most INI1 mosaics, there was no molecular event explaining the partial loss of INI1 IHC. This study confirms the association of hybrid histology with syndromic schwannomas. NF1 mutations did not explain hybrid histology. We did not find BRAF V600E mutations in this series. INI1 mosaic pattern was not explained by reduced copy number or mutation of SMARCB1 as this was found in only a small number of the cases. These results suggest that loss of INI1 expression could be the result of epigenetic alterations.

While human tumors are shaped by the genetic evolution of cancer cells, evidence also suggests that they display superimposed hierarchies reminiscent of normal development. Yet, relating the genetic and hierarchical models of cancer organization has been limited by technical challenges. Here, studying H3 K27M-mutant diffuse midline gliomas we profiled thousands of single cells from fresh tumors by RNA-seq and reconstructed the genetic and epigenetic cellular architecture of tumors in patients. We identified a subset of cancer cells that are differentiated along specialized glial programs while a distinct subset expresses neural stem/progenitor expression programs. By identifying distinct genetic clones, we show that these cellular architectures evolve during the progression of the tumors, suggesting that developmental programs and genetic influences coordinately determine tumor architecture. These results provide unprecedented insight into the cellular composition of midline brain tumors at single-cell resolution and may help harmonize the cancer stem cell and the genetic models of cancer, with critical implications for disease management.

Summary Ependymoma is a heterogeneous entity of central nervous system tumors with well-established molecular groups. Here, we apply single-cell RNA sequencing to analyze ependymomas across molecular groups and anatomic locations to investigate their intratumoral heterogeneity and developmental origins. Ependymomas are composed of a cellular hierarchy initiating from undifferentiated populations, which undergo impaired differentiation toward three lineages of neuronal-glial fate specification. While prognostically favorable groups of ependymoma predominantly harbor differentiated cells, aggressive groups are enriched for undifferentiated cell populations. The delineated transcriptomic signatures correlate with patient survival and define molecular dependencies for targeted treatment approaches. Taken together, our analyses reveal a developmental hierarchy underlying ependymomas relevant to biological and clinical behavior., Graphical Abstract, Highlights • Cellular hierarchies in ependymoma reflect impaired neurodevelopment • Ependymoma cells follow neuronal-, astrocytic-, and ependymal-like trajectories • Undifferentiated programs confer inferior prognosis and are enriched at recurrence • Preclinical data suggest feasibility of inhibiting undifferentiated subpopulations, Gojo et al. use single-cell RNA sequencing to investigate intratumoral heterogeneity and cellular hierarchy in pediatric ependymoma, identifying impaired neurodevelopmental trajectories. The malignant trajectories reveal therapeutic targets and prognostic signatures in ependymoma.