Search

Analysis of tumor gene expression is an important approach for the classification and identification of therapeutic vulnerabilities. However, targeting glioblastoma (GBM) based on molecular subtyping has not yet translated into successful therapies. Here, we present an integrative approach based on molecular pathways to expose new potentially actionable targets. We used gene set enrichment analysis to conduct an unsupervised clustering analysis to condense the gene expression data from bulk patient samples and patient-derived gliomasphere lines into new gene signatures. We identified key targets that are predicted to be differentially activated between tumors and were functionally validated in a library of gliomasphere cultures. Resultant cluster-specific gene signatures associated not only with hallmarks of cell cycle and stemness gene expression, but also with cell type–specific markers and different cellular states of GBM. Several upstream regulators, such as PIK3R1 and EBF1 were differentially enriched in cells bearing stem cell like signatures and bear further investigation. We identified the transcription factor E2F1 as a key regulator of tumor cell proliferation and self-renewal in only a subset of gliomasphere cultures predicted to be E2F1 signaling dependent. Our in vivo work also validated the functional significance of E2F1 in tumor formation capacity in the predicted samples. E2F1 inhibition also differentially sensitized E2F1-dependent gliomasphere cultures to radiation treatment. Our findings indicate that this novel approach exploring cancer pathways highlights key therapeutic vulnerabilities for targeting GBM. Significance: Molecular classification of GBM has not yet resulted in the development of effective therapies. We have developed an integrative approach to identify molecular targets differentially utilized by individual tumors. This approach could lead to patient- and tumor-specific therapeutics.

TCGA gene ontology analysis of clusters shows differentially enriched terms

GO terms and TF enriched in modules associated with GS clusters

Analysis of tumor gene expression is an important approach for the classification and identification of therapeutic vulnerabilities. However, targeting glioblastoma (GBM) based on molecular subtyping has not yet translated into successful therapies. Here, we present an integrative approach based on molecular pathways to expose new potentially actionable targets. We used gene set enrichment analysis to conduct an unsupervised clustering analysis to condense the gene expression data from bulk patient samples and patient-derived gliomasphere lines into new gene signatures. We identified key targets that are predicted to be differentially activated between tumors and were functionally validated in a library of gliomasphere cultures. Resultant cluster-specific gene signatures associated not only with hallmarks of cell cycle and stemness gene expression, but also with cell type–specific markers and different cellular states of GBM. Several upstream regulators, such as PIK3R1 and EBF1 were differentially enriched in cells bearing stem cell like signatures and bear further investigation. We identified the transcription factor E2F1 as a key regulator of tumor cell proliferation and self-renewal in only a subset of gliomasphere cultures predicted to be E2F1 signaling dependent. Our in vivo work also validated the functional significance of E2F1 in tumor formation capacity in the predicted samples. E2F1 inhibition also differentially sensitized E2F1-dependent gliomasphere cultures to radiation treatment. Our findings indicate that this novel approach exploring cancer pathways highlights key therapeutic vulnerabilities for targeting GBM.Significance:Molecular classification of GBM has not yet resulted in the development of effective therapies. We have developed an integrative approach to identify molecular targets differentially utilized by individual tumors. This approach could lead to patient- and tumor-specific therapeutics.

Canonical (H3.1/H3.2) and noncanonical (H3.3) histone 3 K27M-mutant gliomas have unique spatiotemporal distributions, partner alterations, and molecular profiles. The contribution of the cell-of-origin to these differences has been challenging to uncouple from the oncogenic reprogramming induced by the mutation. Here, we perform an integrated analysis of 116 tumors, including single-cell transcriptome and chromatin accessibility, 3D chromatin architecture and epigenomic profiles, and show that K27M-mutant gliomas faithfully maintain chromatin configuration at developmental genes consistent with anatomically distinct oligodendrocyte-precursor-like cells (OPC). H3.3K27M thalamic gliomas map to prosomere 2-derived lineages. In turn, H3.1K27M ACVR1-mutant pontine gliomas uniformly mirror early ventral NKX6-1+/SHH-dependent brainstem OPCs, while H3.3K27M gliomas frequently resemble dorsal PAX3+/BMP-dependent progenitors. Our data suggest a context-specific vulnerability in H3.1K27M-mutant SHH-dependent ventral OPCs, which rely on acquisition of ACVR1 mutations to drive aberrant BMP signaling required for oncogenesis. The unifying action of K27M mutations is to restrict H3K27me3 at PRC2 landing sites, while other epigenetic changes are mainly contingent on the cell-of-origin chromatin state and cycling rate.

The lipidome is currently understudied but fundamental to life. Within the brain, little is known about cell-type lipid heterogeneity, and even less is known about cell-to-cell lipid diversity because it is difficult to study the lipids within individual cells. Here, we used single-cell mass spectrometry-based protocols to profile the lipidomes of 154 910 single cells across ten individuals consisting of five developmental ages and five brain regions, resulting in a unique lipid atlas available via a web browser of the developing human brain. From these data, we identify differentially expressed lipids across brain structures, cortical areas, and developmental ages. We inferred lipid profiles of several major cell types from this data set and additionally detected putative cell-type specific lipids. This data set will enable further interrogation of the developing human brain lipidome.

Glucose is a universal energy currency for living organisms, however, its non-energetic functions in processes such as differentiation are undefined. In epidermis, differentiating cells exhibit dynamic changes in gene expression1–4driven by specific transcription factors (TFs)5–9. The interplay between such TFs and biomolecules that also change in this process is not understood. Metabolomic analyses revealed that increased intracellular glucose accompanies differentiation of epidermal keratinocytes. This elevation also occurred in differentiating cells from other tissues and was verified in epidermal tissue engineered with glucose sensors, which detected a glucose gradient that peaked in the outermost differentiated layers. Free glucose accumulation, unaccompanied by its increased metabolism, was essential for epidermal differentiation and required GLUT1, GLUT3, and SGLT1 transporters. Glucose affinity chromatography and azido-glucose click chemistry uncovered glucose binding to diverse regulatory proteins, including the IRF6 TF, whose epidermal knockout confirmed its requirement in glucose-dependent differentiation. Direct glucose binding enabled IRF6 dimerization, DNA binding, genomic localization, and induction of IRF6 target genes, including essential pro-differentiation TFsGRHL1, GRHL3, HOPXandPRDM1. The IRF6R84Cmutant found in undifferentiated cancers was unable to bind glucose. These data identify a new role for glucose as a gradient morphogen that modulates protein multimerization in cellular differentiation.

The human cortex comprises diverse cell types that emerge from an initially uniform neuroepithelium that gives rise to radial glia, the neural stem cells of the cortex. To characterize the earliest stages of human brain development, we performed single-cell RNA-sequencing across regions of the developing human brain, including the telencephalon, diencephalon, midbrain, hindbrain and cerebellum. We identify nine progenitor populations physically proximal to the telencephalon, suggesting more heterogeneity than previously described, including a highly prevalent mesenchymal-like population that disappears once neurogenesis begins. Comparison of human and mouse progenitor populations at corresponding stages identifies two progenitor clusters that are enriched in the early stages of human cortical development. We also find that organoid systems display low fidelity to neuroepithelial and early radial glia cell types, but improve as neurogenesis progresses. Overall, we provide a comprehensive molecular and spatial atlas of early stages of human brain and cortical development., Eze et al. use single-cell sequencing and immunohistochemical validation to create an atlas of early human brain development. In the telencephalon, they discover a diversity of progenitor subtypes, including two that are enriched in humans.

In children, brain tumors account for the highest mortality rates of all malignancies, and in contrast to adult brain cancers, pediatric tumors often arise in infratentorial regions. Two tumors in particular, H3K27-altered diffuse midline gliomas (DMGs) and group-A posterior fossa ependymomas (PFAs) arise from structures around the fourth ventricle, and strikingly, both tumors exhibit a global decrease in the histone mark histone 3, lysine 27 tri-methylation (H3K27me3) that is associated with transcriptional repression. H3K27-altered DMGs primarily are found primarily in the pons and exhibit by lysine-to-methionine mutations at H3K27 (H3K27M). This missense mutation results in inhibition of the catalytic activity of the polycomb repressive complex 2 (PRC2) which mediates H3K27me3 deposition. PFAs are found in the lateral recess of the fourth ventricle and the cerebellum and demonstrate overexpression of EZH inhibitory protein (EZHIP). EZHIP contains a domain with striking resemblance to the mutant H3K27M histone tail and similarly blocks PRC2 activity. Although H3K27me3 levels are dramatically reduced in both areas, select genes retain H3K27me3 in both tumor types. Increasingly, data have emerged demonstrating that cases of PFAs actually harbor H3K27M mutations and cases of H3-wildtype DMGs in fact overexpress EZHIP. The neighboring origins and overlapping recurrent H3K27me3-modulating alterations observed in these two pediatric brain tumors suggest certain populations of cells in the developing hindbrain may be particularly prone to malignant transformation in the setting of PRC2 inhibition. A systematic comparison of genomic, transcriptomic, and epigenomic datasets for both tumors identified key molecular similarities between subsets of DMGs and PFAs and highlighted how developmental signatures may be preserved in tumor expression profiles. Specifically, these analyses identify shared copy-number profiles between PFAs and different subtypes of DMGs and demonstrate that the Activin A Receptor Type 1 (ACVR1), a gene with recurrent activating mutations in H3K27-altered DMGs, is overexpressed in a subset PFAs with worse survival outcomes. Additionally, interrogation of the H3K27me3 and enhancer landscapes identified critical genes with similar chromatin profiles and revealed heterogeneity in repression of H3K27me3-marked genes that correlated with patterns of inferred anatomic origins for each tumor. Moreover, one of these genes, cellular retinoic acid binding partner 1 (CRABP1) showed enrichment in progenitor cell types in single-cell transcriptomic datasets. Together, these analyses characterize the common and unique molecular features of H3K27-altered DMGs and PFAs and shed light on the cellular contexts of the developing hindbrain that lead to development of each tumor. Citation Format: Matthew Pun, Drew Pratt, Patricia R. Nano, Piyush K. Joshi, Li Jiang, Bernhard Englinger, Arvind Rao, Marcin Cieslik, Arul M. Chinnaiyan, Kenneth Aldape, Stefan Pfister, Mariella G. Filbin, Aparna Bhaduri, Sriram Venneti. A systematic comparison of molecular features shared by H3K27-altered diffuse midline gliomas and posterior fossa A ependymomas. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3523.

BACKGROUND The Hedgehog pathway directs gene expression programs that are essential for glioblastoma stem cells in vitro, but a clinical trial suggested Hedgehog pathway inhibition does not improve glioblastoma outcomes (ABTC-0904). In pancreatic and bladder cancer, Hedgehog signaling through the tumor microenvironment induces stromal morphogens to inhibit cancer stem cell proliferation. Thus, we hypothesized Hedgehog signaling through the glioblastoma microenvironment may inhibit glioblastoma stem cell proliferation in vivo. METHODS Hedgehog signaling was studied using patient-derived, fluorescently-labeled IDH wild-type glioblastoma cells (GBM6, SF11956, SF11907, GPMP004, GPMP005) in vitro, in vivo after intracranial implantation in mice, or in 3D co-culture with iPSC-derived organoids comprised of human astrocytes or neurons. Genetic loss-of-function experiments were performed using CRISPR interference. Cell proliferation, stem cell marker expression, Sonic Hedgehog (SHH) expression, or primary cilia architecture were assessed using immunofluorescence, confocal microscopy, or live-cell imaging. Hedgehog signaling or morphogen expression were assessed using QPCR or single-cell RNA sequencing. Pharmacologic experiments were performed using vismodegib to inhibit the Hedgehog pathway, or FK506 to induce stromal morphogen expression. RESULTS Vismodegib increased glioblastoma cell proliferation and stem cell marker expression in vivo or in co-culture with astrocyte organoids. Single-cell RNA sequencing demonstrated glioblastoma cell co-culture with astrocyte organoids induced Hedgehog target genes and stem cell markers in glioblastoma cells, or stromal morphogens in astrocytes, which also expressed primary cilia. CRISPRi suppression of SHH in glioblastoma cells increased glioblastoma cell proliferation and stem cell marker expression in co-culture with astrocyte organoids. Neither vismodegib nor SHH suppression induced glioblastoma cell proliferation or stem cell marker expression in vitro or in co-culture with neuron organoids. FK506 induced astrocyte morphogen expression, inhibiting glioblastoma cell proliferation, stem cell marker expression, and the effects of vismodegib in co-culture with astrocyte organoids. CONCLUSIONS SHH from glioblastoma cells signals through astrocytes to inhibit cancer stem cell proliferation.

Innovations in organoid-based models of human tissues have made them an exciting experimental platform for studying development and disease. However, these models require systematic benchmarking against primary tissue to establish their value. We discuss key parameters that impact the utility of organoid models, primarily focusing on cerebral organoids as examples.

Summary As the use of single-cell technologies has grown, so has the need for tools to explore these large, complicated datasets. The UCSC Cell Browser is a tool that allows scientists to visualize gene expression and metadata annotation distribution throughout a single-cell dataset or multiple datasets. Availability and implementation We provide the UCSC Cell Browser as a free website where scientists can explore a growing collection of single-cell datasets and a freely available python package for scientists to create stable, self-contained visualizations for their own single-cell datasets. Learn more at https://cells.ucsc.edu. Supplementary information Supplementary data are available at Bioinformatics online.

The human brain is subdivided into distinct anatomical structures, including the neocortex, which in turn encompasses dozens of distinct specialized cortical areas. Early morphogenetic gradients are known to establish early brain regions and cortical areas, but how early patterns result in finer and more discrete spatial differences remains poorly understood1. Here we use single-cell RNA sequencing to profile ten major brain structures and six neocortical areas during peak neurogenesis and early gliogenesis. Within the neocortex, we find that early in the second trimester, a large number of genes are differentially expressed across distinct cortical areas in all cell types, including radial glia, the neural progenitors of the cortex. However, the abundance of areal transcriptomic signatures increases as radial glia differentiate into intermediate progenitor cells and ultimately give rise to excitatory neurons. Using an automated, multiplexed single-molecule fluorescent in situ hybridization approach, we find that laminar gene-expression patterns are highly dynamic across cortical regions. Together, our data suggest that early cortical areal patterning is defined by strong, mutually exclusive frontal and occipital gene-expression signatures, with resulting gradients giving rise to the specification of areas between these two poles throughout successive developmental timepoints., RNA-sequencing analysis of the prenatal human brain at different stages of development shows that areal transcriptional signatures are dynamic and coexist with developmental and cell-type signatures.

The glioblastoma microenvironment is an immunosuppressive barrier to therapeutic innovation. We hypothesized intratumor convection enhanced delivery (CED) of gene therapy vectors could reprogram the glioblastoma immune microenvironment and elucidate therapeutic vulnerabilities. To test this, SB28 or GL261 glioblastoma allografts were implanted into immunocompetent mice and treated with CED of attenuated adeno-associated virus 9 vectors (AAV9) encoding experimental cytokines (Il1b, Ccl4, or Apoa1) underlying infiltration or activation of anti-tumor immune cells in other intracranial tumors. Serial intracranial bioluminescence was used to assess glioblastoma growth and animals were monitored for survival. The impact of gene therapy perturbations on the glioblastoma immune microenvironment was assessed using histology, immunohistochemistry, single-cell mass cytometry (CyTOF), and multiplexed cytokine assays. Serial body weight and systemic cytokine measurements showed no evidence of treatment toxicity. AAV9-APOA1 or AAV9-IL1B gene therapy CED treatments attenuated SB28 growth and prolonged survival, decreasing immunosuppressive macrophage infiltration and increasing CD8 T cell and microglia infiltration of the glioblastoma microenvironment compared to control AAV9 vectors. Gene therapy CED treatments did not attenuate GL261 growth or prolong survival, but CyTOF of human glioblastomas (n=6) in comparison to preclinical models revealed untreated GL261 glioblastomas were endogenously enriched in CD8 T cells and other lymphoid lineages compared to untreated SB28 or human glioblastomas. Multiplexed cytokine assays demonstrated suppression of intratumor IL6 is a conserved mechanism of action underlying glioblastoma gene therapy responses. Single-cell RNA sequencing analysis of 32,877 cells from human glioblastomas (n=11) showed IL6 is predominantly produced by radial glial like cancer stem cells or endothelial cells in the tumor microenvironment. In support of these findings, survival from intracranial SB28 glioblastomas was prolonged in Il6 knockout C57BL/6J mice compared to wildtype mice. In summary, we report a novel strategy using gene therapy and CED to reprogram the glioblastoma immune microenvironment, revealing IL6 drives glioblastoma immunosuppression and growth.

There are no approved targeted therapies for meningiomas and the cell types underlying meningeal tumorigenesis are incompletely understood. To address these limitations, we performed single-cell RNA sequencing of 57,114 cells from 8 human meningiomas and 54,607 cells from 3 canine meningiomas. Pseudotime, gene ontology, and copy number variant analyses revealed a population of pericyte-like meningioma cells that were conserved across human and canine tumors and were enriched in expression of Notch3 and other cancer stem cell genes. Deconvolution of cell types from bulk RNA sequencing and DNA methylation profiling of 200 human meningiomas integrated with immunohistochemistry (IHC), immunofluorescence (IF), and RNAScope demonstrated Notch3+ pericytes and Notch3 expression were enriched in high grade or Immune-enriched meningiomas, which were distinguished from other meningioma DNA methylation groups by genes driving vasculature development. IHC and IF of human meninges integrated with lineage tracing approaches using Notch3-CreERT2 ROSAmT/mG alleles in mice demonstrated Notch3 expression was restricted to the perivascular stem cell niche during meningeal development and homeostasis. Mice harboring Notch3-CreERT2 Nf2fl/fl alleles developed meningeal hyperproliferation. Overexpression of constitutively activated Notch3 (Notch3AICD) in Immune-enriched human meningioma cells increased the expression of cancer stem cell genes, driving clonogenic growth in vitro, limiting dilution tumor-initiating capacity in vivo, and resistance to radiotherapy in vivo. A selective Notch3 neutralizing antibody (αNRR3) blocked meningioma cell proliferation and expression of Notch3 target genes, inhibiting meningioma xenograft growth and prolonging overall survival. Single-cell RNA sequencing of 187,366 cells from meningioma xenografts after αNRR3 or radiotherapy treatment ± Notch3AICD overexpression revealed distinct meningioma cell-intrinsic or cell-extrinsic mechanisms driving responses to radiotherapy or αNRR3, respectively. Combined treatment with αNRR3 and radiotherapy additively blocked meningioma xenograft growth and extended survival benefit. In sum, these data shed light on a novel cell type, molecular mechanism, and therapeutic vulnerability in the most common primary intracranial tumor.

Interactions between angiogenesis and neurogenesis regulate embryonic brain development. However, a comprehensive understanding of the stages of vascular cell maturation is lacking, especially in the prenatal human brain. Using fluorescence-activated cell sorting, single-cell transcriptomics, and histological and ultrastructural analyses, we show that an ensemble of endothelial and mural cell subtypes tile the brain vasculature during the second trimester. These vascular cells follow distinct developmental trajectories and utilize diverse signaling mechanisms, including collagen, laminin, and midkine, to facilitate cell-cell communication and maturation. Interestingly, our results reveal that tip cells, a subtype of endothelial cells, are highly enriched near the ventricular zone, the site of active neurogenesis. Consistent with these observations, prenatal vascular cells transplanted into cortical organoids exhibit restricted lineage potential that favors tip cells, promotes neurogenesis, and reduces cellular stress. Together, our results uncover important mechanisms into vascular maturation during this critical period of human brain development.

The development and diversity of neuronal subtypes in the human hypothalamus has been insufficiently characterized. We sequenced the transcriptomes of 40,927 cells from the prenatal human hypothalamus spanning from 6 to 25 gestational weeks and 25,424 mature neurons in regions of the adult human hypothalamus, revealing a temporal trajectory from proliferative stem cell populations to mature neurons and glia. Developing hypothalamic neurons followed branching trajectories leading to 170 transcriptionally distinct neuronal subtypes in ten hypothalamic nuclei in the adult. The uniqueness of hypothalamic neuronal lineages was examined developmentally by comparing excitatory lineages present in cortex and inhibitory lineages in ganglionic eminence from the same individuals, revealing both distinct and shared drivers of neuronal maturation across the human forebrain. Cross-species comparisons to the mouse hypothalamus identified human-specificPOMCpopulations expressing unique combinations of transcription factors and neuropeptides. These results provide the first comprehensive transcriptomic view of human hypothalamus development at cellular resolution.One-Sentence SummaryUsing single-cell genomics, we reconstructed the developmental lineages by which precursor populations give rise to 170 distinct neuronal subtypes in the human hypothalamus.

Targeting glioblastoma (GBM) based on molecular subtyping have not yet translated into successful therapies. Here, we used gene set enrichment analysis (GSEA) to conduct an unsupervised clustering analysis to condense the gene expression data from bulk patient samples and patient-derived gliomasphere lines into new gene lists. We then identified key molecular pathways differentially regulated between tumors. These gene lists associated not only with cell cycle and stemness signatures, but also with cell-type specific markers and different cellular states of GBM. We identified the transcription factor E2F1 as a key regulator of tumor cell proliferation and self-renewal in only the subset of proliferating gliomasphere cultures predicted to be E2F1-activated and validated its functional significance in tumor formation capacity. E2F1 inhibition also sensitized E2F1-activated gliomasphere cultures to radiation treatment. Our findings indicate that a pathway-based approach can be leveraged to deconstruct inter-tumoral heterogeneity and uncover key therapeutic vulnerabilities for targeting GBM.

Studies of the spatiotemporal, transcriptomic, and morphological diversity of radial glia (RG) have spurred our current models of human corticogenesis. In the developing cortex, neural intermediate progenitor cells (nIPCs) are a neuron-producing transit-amplifying cell type born in the germinal zones of the cortex from RG. The potential diversity of the nIPC population, that produces a significant portion of excitatory cortical neurons, is understudied, particularly in the developing human brain. Here we explore the spatiotemporal, transcriptomic, and morphological variation that exists within the human nIPC population and provide a resource for future studies. We observe that the spatial distribution of nIPCs in the cortex changes abruptly around gestational week (GW) 19/20, marking a distinct shift in cellular distribution and organization during late neurogenesis. We also identify five transcriptomic subtypes, one of which appears at this spatiotemporal transition. Finally, we observe a diversity of nIPC morphologies that do not correlate with specific transcriptomic subtypes. These results provide an analysis of the spatiotemporal, transcriptional, and morphological diversity of nIPCs in developing brain tissue and provide an atlas of nIPC subtypes in the developing human cortex that can benchmark in vitro models of human development such as cerebral organoids and help inform future studies of how nIPCs contribute to cortical neurogenesis.

The human brain is subdivided into distinct anatomical structures. The neocortex, one of these structures, enables higher-order sensory, associative, and cognitive functions, and in turn encompasses dozens of distinct specialized cortical areas. Early morphogenetic gradients are known to establish an early blueprint for the specification of brain regions and cortical areas. Furthermore, recent studies have uncovered distinct transcriptomic signatures between opposing poles of the developing neocortex1. However, how early, broad developmental patterns result in finer and more discrete spatial differences across the adult human brain remains poorly understood2. Here, we use single-cell RNA-sequencing to profile ten major brain structures and six neocortical areas during peak neurogenesis and early gliogenesis. Our data reveal that distinct cell subtypes are predominantly brain-structure specific. Within the neocortex, we find that even early in the second trimester, a large number of genes are differentially expressed across distinct cortical areas in all cell types, including radial glia, the neural progenitors of the cortex. However, the abundance of areal transcriptomic signatures increases as radial glia differentiate into intermediate progenitor cells and ultimately give rise to excitatory neurons. Using an automated, multiplexed single-molecule fluorescent in situ hybridization (smFISH) approach, we validated the expression pattern of area-specific neuronal genes and also discover that laminar gene expression patterns are highly dynamic across cortical regions. Together, our data suggest that early cortical areal patterning is defined by strong, mutually exclusive frontal and occipital gene expression signatures, with resulting gradients giving rise to the specification of areas between these two poles throughout successive developmental timepoints.

BACKGROUND Meningiomas are the most common primary intracranial tumors in humans and dogs, but biologic drivers and cell types underlying meningeal tumorigenesis are incompletely understood. Here we integrate meningioma single-cell RNA sequencing with stem cell approaches to define a perivascular stem cell underlying vertebrate meningeal tumorigenesis. METHODS Single-cell RNA sequencing was performed on 57,114 cells from 8 human meningiomas, 54,607 cells from 3 dog meningiomas, and human meningioma xenografts in mice. Results were validated using immunofluorescence (IF), immunohistochemistry (IHC), and deconvolution of bulk RNA sequencing of 200 human meningiomas. Mechanistic and functional studies were performed using clonogenic and limiting dilution assays, xenografts, and genetically engineered mouse models. RESULTS Copy number variant identification from human meningioma single cells distinguished tumor cells with loss of chr22q from non-tumor cells with intact chr22q. A single cluster distinguished by expression of Notch3 and other cancer stem cell genes had an intermediate level of loss of chr22q, suggesting this cluster may represent meningioma stem cells. In support of this hypothesis, pseudotime trajectory analysis demonstrated transcriptomic progression starting from Notch3+ cells and encompassing all other meningioma cell types. Notch3+ meningioma cells had transcriptomic concordance to mural pericytes, and IF/IHC of prenatal and adult human meninges, as well as lineage tracing using a Notch3-CreERT2 allele in mice, confirmed Notch3+ cells were restricted to the perivascular stem cell niche in mammalian meningeal development and homeostasis. Integrating human and dog meningioma single cells revealed Notch3+ cells in tumor and non-tumor clusters in dog meningiomas. Notch3 IF/IHC and cell-type deconvolution of bulk RNA sequencing showed Notch3+ cells were enriched in high-grade human meningiomas. Notch3 overexpression in human meningioma cells increased clonogenic growth in vitro, and increased tumorigenesis and tumor growth in vivo, decreasing overall survival. CONCLUSIONS Notch3+ stem cells in the perivascular niche underlie vertebrate meningeal tumorigenesis.

Human embryonic stem cell-derived retinal pigment epithelial cells (hESC-RPE) are a promising cell source to treat age-related macular degeneration (AMD). Despite several ongoing clinical studies, detailed single cell mapping of the transient cellular and molecular dynamics from the pluripotent state to mature RPE has not been performed. Here we conduct single-cell transcriptomic analyses of 25,718 cells during differentiation as well as in embryonic and adult retina references, revealing differentiation progression through an un-expected initial cell diversification recapitulating early embryonic development before converging towards an RPE lineage. We also identified NCAM1 to track and capture an intermediate retinal progenitor with the potential to give rise to multiple neuroepithelial lineages. Finally, we profiled hESC-RPE cells after subretinal transplantation into the rabbit eye, uncovering robustin vivomaturation towards an adult state. Our detailed evaluation of hESC-RPE differentiation supports the development of safe and efficient pluripotent stem cell-based therapies for AMD.

SummaryAs the use of single-cell technologies has grown, so has the need for tools to explore these large, complicated datasets. The UCSC Cell Browser is a tool that allows scientists to visualize gene expression and metadata annotation distribution throughout a single-cell dataset or multiple datasets.Availability and implementationWe provide the UCSC Cell Browser as a free website where users can explore a growing collection of single-cell datasets and a freely available python package for scientists to create stable, self-contained visualizations for their own single-cell datasets. Learn more at https://cells.ucsc.edu.Contactcells@ucsc.edu

Author(s): Crawford, Emily D; Acosta, Irene; Ahyong, Vida; Anderson, Erika C; Arevalo, Shaun; Asarnow, Daniel; Axelrod, Shannon; Ayscue, Patrick; Azimi, Camillia S; Azumaya, Caleigh M; Bachl, Stefanie; Bachmutsky, Iris; Bhaduri, Aparna; Brown, Jeremy Bancroft; Batson, Joshua; Behnert, Astrid; Boileau, Ryan M; Bollam, Saumya R; Bonny, Alain R; Booth, David; Borja, Michael Jerico B; Brown, David; Buie, Bryan; Burnett, Cassandra E; Byrnes, Lauren E; Cabral, Katelyn A; Cabrera, Joana P; Caldera, Saharai; Canales, Gabriela; Castaneda, Gloria R; Chan, Agnes Protacio; Chang, Christopher R; Charles-Orszag, Arthur; Cheung, Carly; Chio, Unseng; Chow, Eric D; Citron, Y Rose; Cohen, Allison; Cohn, Lillian B; Chiu, Charles; Cole, Mitchel A; Conrad, Daniel N; Constantino, Angela; Cote, Andrew; Crayton-Hall, Tre'Jon; Darmanis, Spyros; Detweiler, Angela M; Dial, Rebekah L; Dong, Shen; Duarte, Elias M; Dynerman, David; Egger, Rebecca; Fanton, Alison; Frumm, Stacey M; Fu, Becky Xu Hua; Garcia, Valentina E; Garcia, Julie; Gladkova, Christina; Goldman, Miriam; Gomez-Sjoberg, Rafael; Gordon, M Grace; Grove, James CR; Gupta, Shweta; Haddjeri-Hopkins, Alexis; Hadley, Pierce; Haliburton, John; Hao, Samantha L; Hartoularos, George; Herrera, Nadia; Hilberg, Melissa; Ho, Kit Ying E; Hoppe, Nicholas; Hosseinzadeh, Shayan; Howard, Conor J; Hussmann, Jeffrey A; Hwang, Elizabeth; Ingebrigtsen, Danielle; Jackson, Julia R; Jowhar, Ziad M; Kain, Danielle; Kim, James YS; Kistler, Amy; Kreutzfeld, Oriana; Kulsuptrakul, Jessie; Kung, Andrew F

Human cerebral cortex size and complexity has increased greatly during evolution. While increased progenitor diversity and enhanced proliferative potential play important roles in human neurogenesis and gray matter expansion, the mechanisms of human oligodendrogenesis and white matter expansion remain largely unknown. Here, we identify EGFR-expressing “Pre-OPCs” that originate from outer radial glial cells (oRGs) and undergo mitotic somal translocation (MST) during division. oRG-derived Pre-OPCs provide an additional source of human cortical oligodendrocyte precursor cells (OPCs) and define a lineage trajectory. We further show that human OPCs undergo consecutive symmetric divisions to exponentially increase the progenitor pool size. Additionally, we find that the OPC-enriched gene, PCDH15, mediates daughter cell repulsion and facilitates proliferation. These findings indicate properties of OPC derivation, proliferation, and dispersion important for human white matter expansion and myelination.

The human cortex is comprised of diverse cell types that emerge from an initially uniform neuroepithelium that gives rise to radial glia, the neural stem cells of the cortex. To characterize the earliest stages of human brain development, we performed single-cell RNA-sequencing across regions of the developing human brain, including the telencephalon, diencephalon, midbrain, hindbrain, and cerebellum. We identify nine progenitor populations physically proximal to the telencephalon, suggesting more heterogeneity than previously described, including a highly prevalent mesenchymal-like population that disappears once neurogenesis begins. Comparison of human and mouse progenitor populations at corresponding stages identifies two progenitor clusters that are enriched in the early stages of human cortical development. We also find that organoid systems display low fidelity to neuroepithelial and early radial glia cell types, but improve as neurogenesis progresses. Overall, we provide a comprehensive molecular and spatial atlas of early stages of human brain and cortical development.

SUMMARY Neural stem cells (NSCs) in the prenatal neocortex progressively generate different subtypes of glutamatergic projection neurons. Following that, NSCs have a major switch in their progenitor properties and produce γ-aminobutyric acid (GABAergic) interneurons for the olfactory bulb (OB), cortical oligodendrocytes, and astrocytes. Herein, we provide evidence for the molecular mechanism that underlies this switch in the state of neocortical NSCs. We show that, at around E16.5, mouse neocortical NSCs start to generate GSX2-expressing (GSX2+) intermediate progenitor cells (IPCs). In vivo lineage-tracing study revealed that GSX2+ IPC population gives rise not only to OB interneurons but also to cortical oligodendrocytes and astrocytes, suggesting that they are a tri-potential population. We demonstrated that Sonic hedgehog signaling is both necessary and sufficient for the generation of GSX2+ IPCs by reducing GLI3R protein levels. Using single-cell RNA sequencing, we identify the transcriptional profile of GSX2+ IPCs and the process of the lineage switch of cortical NSCs., Graphical Abstract, In Brief Zhang et al. reveal that cortical radial glia-derived GSX2+ cells at the late embryonic stage are tri-potential intermediate progenitors, which give rise to a subset of cortical oligodendrocytes, astrocytes, and olfactory bulb interneurons. SHH signaling is crucial for the generation of GSX2+ cells by reducing GLI3R protein level.

Cortical organoids are self-organizing three-dimensional cultures that model features of the developing human cerebral cortex1,2. However, the fidelity of organoid models remains unclear3-5. Here we analyse the transcriptomes of individual primary human cortical cells from different developmental periods and cortical areas. We find that cortical development is characterized by progenitor maturation trajectories, the emergence of diverse cell subtypes and areal specification of newborn neurons. By contrast, organoids contain broad cell classes, but do not recapitulate distinct cellular subtype identities and appropriate progenitor maturation. Although the molecular signatures of cortical areas emerge in organoid neurons, they are not spatially segregated. Organoids also ectopically activate cellular stress pathways, which impairs cell-type specification. However, organoid stress and subtype defects are alleviated by transplantation into the mouse cortex. Together, these datasets and analytical tools provide a framework for evaluating and improving the accuracy of cortical organoids as models of human brain development.

International audience; We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

The brain contains immense cellular diversity. Early neuroanatomists identified broad structural brain features and used a variety of stains to characterize, draw, and name a number of constituent cell types. However, our granular understanding of what defines cellular identity, and how this relates to human neurogenesis, has been transformed by the sequencing of the human genome and subsequent technological advances. The genome enabled a translation of single gene level studies to genome-wide, unbiased characterization of gene expression patterns. Rapidly developing gene sequencing technology has revolutionized our understanding of cellular diversity, first based on microarray technology, then RNA-sequencing of 3′ fragments, to more recent whole-transcriptome sequencing, and ultimately single-cell sequencing. Single-cell sequencing allows for a heterogeneous tissue, like the brain, to be broken down into its component cell types based on fundamentally relevant single cell level gene programs. The advent of single-cell sequencing has been transformative for many fields, but using single-cell sequencing in the study of the brain and specifically brain development has enabled a convergence of the understanding of gene expression with features such as area, layer, age, connectivity, morphology, and ultimately cellular function, to achieve a core foundational understanding of cell type identity. While a number of these approaches are now being applied to human development, they have already been widely applied to rodent models, especially the mouse. Due to limitations of accessibility of human brain tissue, in vitro model systems such as the cerebral organoid have been developed to further our understanding of cellular and molecular developmental programs and to functionally test hypotheses gathered from sequencing studies. Here, we focus on established and emerging single-cell technologies, how they have enabled discovery in human brain development, and how organoid models recapitulate cellular developmental trajectories.

Glioblastoma is a devastating form of brain cancer. To identify aspects of tumor heterogeneity that may illuminate drivers of tumor invasion, we created a glioblastoma tumor cell atlas with single-cell transcriptomics of cancer cells mapped onto a reference framework of the developing and adult human brain. We find that multiple GSC subtypes exist within a single tumor. Within these GSCs, we identify an invasive cell population similar to outer radial glia (oRG), a fetal cell typethat expands the stem cell niche in normal human cortex. Using live time-lapse imaging of primary resected tumors, we discover that tumor-derived oRG-like cells undergo characteristic mitotic somal translocation behavior previously only observed inhuman development, suggesting areactivation ofdevelopmental programs. In addition, we showthat PTPRZ1 mediates both mitotic somal translocation and glioblastoma tumor invasion. These data suggest that the presence of heterogeneous GSCs may underlie glioblastoma's rapid progression and invasion.

Genetic changes causing brain size expansion in human evolution have remained elusive. Notch signaling is essential for radial glia stem cell proliferation and is a determinant of neuronal number in the mammalian cortex. We find that three paralogs of human-specific NOTCH2NL are highly expressed in radial glia. Functional analysis reveals that different alleles of NOTCH2NL have varying potencies to enhance Notch signaling by interacting directly with NOTCH receptors. Consistent with a role in Notch signaling, NOTCH2NL ectopic expression delays differentiation of neuronal progenitors, while deletion accelerates differentiation into cortical neurons. Furthermore, NOTCH2NL genes provide the breakpoints in 1q21.1 distal deletion/duplication syndrome, where duplications are associated with macrocephaly and autism and deletions with microcephaly and schizophrenia. Thus, the emergence of human-specific NOTCH2NL genes may have contributed to the rapid evolution of the larger human neocortex, accompanied by loss of genomic stability at the 1q21.1 locus and resulting recurrent neurodevelopmental disorders.

Building a brain The human brain is built in an inside-out manner as a series of layers. Although progenitor cells spin off new neurons in a seemingly organized fashion, the devil is in the details. Nowakowski et al. analyzed the transcriptomes of single cells from the developing brain to elucidate the hidden complexity of brain construction. For each cell, its position within the brain matters, as well as what type of neuron is being made at what point during overall development. These individual expression patterns result in organized diversity in the brain's cortex. Science , this issue p. 1318

Multiple sclerosis (MS) is a neuroinflammatory disease with a relapsing-remitting disease course at early stages, distinct lesion characteristics in cortical grey versus subcortical white matter and neurodegeneration at chronic stages. Here we used single-nucleus RNA sequencing to assess changes in expression in multiple cell lineages in MS lesions and validated the results using multiplex in situ hybridization. We found selective vulnerability and loss of excitatory CUX2-expressing projection neurons in upper-cortical layers underlying meningeal inflammation; such MS neuron populations exhibited upregulation of stress pathway genes and long non-coding RNAs. Signatures of stressed oligodendrocytes, reactive astrocytes and activated microglia mapped most strongly to the rim of MS plaques. Notably, single-nucleus RNA sequencing identified phagocytosing microglia and/or macrophages by their ingestion and perinuclear import of myelin transcripts, confirmed by functional mouse and human culture assays. Our findings indicate lineage- and region-specific transcriptomic changes associated with selective cortical neuron damage and glial activation contributing to progression of MS lesions.

In the developing human neocortex, progenitor cells generate diverse cell types prenatally. Progenitor cells and newborn neurons respond to signaling cues, including neurotransmitters. While single-cell RNA sequencing has revealed cellular diversity, physiological heterogeneity has yet to be mapped onto these developing and diverse cell types. By combining measurements of intracellular Ca2+ elevations in response to neurotransmitter receptor agonists and RNA sequencing of the same single cells, we show that Ca2+ responses are cell-type-specific and change dynamically with lineage progression. Physiological response properties predict molecular cell identity and additionally reveal diversitynot captured by single-cell transcriptomics. We find that the serotonin receptor HTR2A selectively activates radial glia cells in the developing human, butnot mouse, neocortex, and inhibiting HTR2A receptorsin human radial glia disrupts the radial glial scaffold. We show highly specific neurotransmittersignaling during neurogenesis in the developing human neocortex and highlight evolutionarilydivergent mechanisms of physiological signaling.

The cerebral cortex is the hallmark of the mammalian nervous system, and its large size and cellular diversity in humans supports our most sophisticated cognitive abilities. Although the basic cellular organization of the cortex is conserved across mammals, cells have diversified during evolution. An increasingly integrated taxonomy of cell types, especially with the advent of single cell transcriptomic data, has revealed an unprecedented variety of human cortical cell subtypes. Here, we broadly review the cellular composition and diversity of the mammalian brain, and how progenitor pools generate cell subtypes during development. We then discuss human cortical cells that are distinct from rodent cells, as well as the challenges and advantages of using model systems to study human cell types in health and disease.