Search

Glioblastoma (GBM) is the deadliest form of primary brain tumor with limited treatment options. Recent studies have profiled GBM tumor heterogeneity, revealing numerous axes of variation that explain the molecular and spatial features of the tumor. Here, we seek to bridge descriptive characterization of GBM cell type heterogeneity with the functional role of individual populations within the tumor. Our lens leverages a gene program-centric meta-atlas of published transcriptomic studies to identify commonalities between diverse tumors and cell types in order to decipher the mechanisms that drive them. This approach led to the discovery of a tumor-derived stem cell population with mixed vascular and neural stem cell features, termed a neurovascular progenitor (NVP). Following in situ validation and molecular characterization of NVP cells in GBM patient samples, we characterized their function in vivo. Genetic depletion of NVP cells resulted in altered tumor cell composition, fewer cycling cells, and extended survival, underscoring their critical functional role. Clonal analysis of primary patient tumors in a human organoid tumor transplantation system demonstrated that the NVP has dual potency, generating both neuronal and vascular tumor cells. Although NVP cells comprise a small fraction of the tumor, these clonal analyses demonstrated that they strongly contribute to the total number of cycling cells in the tumor and generate a defined subset of the whole tumor. This study represents a paradigm by which cell type-specific interrogation of tumor populations can be used to study functional heterogeneity and therapeutically targetable vulnerabilities of GBM.

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

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

CD19-directed immunotherapies are clinically effective for treating B cell malignancies but also cause a high incidence of neurotoxicity. A subset of patients treated with chimeric antigen receptor (CAR) T cells or bispecific T cell engager (BiTE) antibodies display severe neurotoxicity, including fatal cerebral edema associated with T cell infiltration into the brain. Here, we report that mural cells, which surround the endothelium and are critical for blood-brain-barrier integrity, express CD19. We identify CD19 expression in brain mural cells using single-cell RNA sequencing data and confirm perivascular staining at the protein level. CD19 expression in the brain begins early in development alongside the emergence of mural cell lineages and persists throughout adulthood across brain regions. Mouse mural cells demonstrate lower levels of Cd19 expression, suggesting limitations in preclinical animal models of neurotoxicity. These data suggest an on-target mechanism for neurotoxicity in CD19-directed therapies and highlight the utility of human single-cell atlases for designing immunotherapies.

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.

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 type that 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 in human development, suggesting a reactivation of developmental programs. In addition, we show that 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.

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.

The cerebral cortex is the hallmark of the mammalian nervous system, and its large size and cellular diversity in humans support 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.

Despite the clinical and genetic heterogeneity of autism, bulk gene expression studies show that changes in the neocortex of autism patients converge on common genes and pathways. However, direct assessment of specific cell types in the brain affected by autism has not been feasible until recently. We used single-nucleus RNA sequencing of cortical tissue from patients with autism to identify autism-associated transcriptomic changes in specific cell types. We found that synaptic signaling of upper-layer excitatory neurons and the molecular state of microglia are preferentially affected in autism. Moreover, our results show that dysregulation of specific groups of genes in cortico-cortical projection neurons correlates with clinical severity of autism. These findings suggest that molecular changes in upper-layer cortical circuits are linked to behavioral manifestations of autism.

BackgroundHigh throughput methods for profiling the transcriptomes of single cells have recently emerged as transformative approaches for large-scale population surveys of cellular diversity in heterogeneous primary tissues. However, the efficient generation of such atlases will depend on sufficient sampling of diverse cell types while remaining cost-effective to enable a comprehensive examination of organs, developmental stages, and individuals.ResultsTo examine the relationship between sampled cell numbers and transcriptional heterogeneity in the context of unbiased cell type classification, we explored the population structure of a publicly available 1.3 million cell dataset from E18.5 mouse brain and validated our findings in published data from adult mice. We propose a computational framework for inferring the saturation point of cluster discovery in a single-cell mRNA-seq experiment, centered around cluster preservation in downsampled datasets. In addition, we introduce a "complexity index," which characterizes the heterogeneity of cells in a given dataset. Using Cajal-Retzius cells as an example of a limited complexity dataset, we explored whether the detected biological distinctions relate to technical clustering. Surprisingly, we found that clustering distinctions carrying biologically interpretable meaning are achieved with far fewer cells than the originally sampled, though technical saturation of rare populations such as Cajal-Retzius cells is not achieved. We additionally validated these findings with a recently published atlas of cell types across mouse organs and again find using subsampling that a much smaller number of cells recapitulates the cluster distinctions of the complete dataset.ConclusionsTogether, these findings suggest that most of the biologically interpretable cell types from the 1.3 million cell database can be recapitulated by analyzing 50,000 randomly selected cells, indicating that instead of profiling few individuals at high "cellular coverage," cell atlas studies may instead benefit from profiling more individuals, or many time points at lower cellular coverage and then further enriching for populations of interest. This strategy is ideal for scenarios where cost and time are limited, though extremely rare populations of interest (

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.

Systematic analyses of spatiotemporal gene expression trajectories during organogenesis have been challenging because diverse cell types at different stages of maturation and differentiation coexist in the emerging tissues. We identified discrete cell types as well as temporally and spatially restricted trajectories of radial glia maturation and neurogenesis in developing human telencephalon. These lineage-specific trajectories reveal the expression of neurogenic transcription factors in early radial glia and enriched activation of mammalian target of rapamycin signaling in outer radial glia. Across cortical areas, modest transcriptional differences among radial glia cascade into robust typological distinctions among maturing neurons. Together, our results support a mixed model of topographical, typological, and temporal hierarchies governing cell-type diversity in the developing human telencephalon, including distinct excitatory lineages emerging in rostral and caudal cerebral cortex.

Pediatric-type follicular lymphoma and pediatric marginal zone lymphoma are two of the rarest B-cell lymphomas. These lymphomas occur predominantly in the pediatric population and show features distinct from their more common counterparts in adults: adult-type follicular lymphoma and adult-type nodal marginal zone lymphoma. Here we report a detailed whole-exome deep sequencing analysis of a cohort of pediatric-type follicular lymphomas and pediatric marginal zone lymphomas. This analysis revealed a recurrent somatic variant encoding p.Lys66Arg in the transcription factor interferon regulatory factor 8 (IRF8) in 3 of 6 cases (50%) of pediatric-type follicular lymphoma. This specific point mutation was not detected in pediatric marginal zone lymphoma or in adult-type follicular lymphoma. Additional somatic point mutations in pediatric-type follicular lymphoma were observed in genes involved in transcription, intracellular signaling, and cell proliferation. In pediatric marginal zone lymphoma, no recurrent mutation was identified; however, somatic point mutations were observed in genes involved in cellular adhesion, cytokine regulatory elements, and cellular proliferation. A somatic variant in AMOTL1, a recurrently mutated gene in splenic marginal zone lymphoma, was also identified in a case of pediatric marginal zone lymphoma. The overall non-synonymous mutational burden was low in both pediatric-type follicular lymphoma and pediatric marginal zone lymphoma (4.6 mutations per exome). Altogether, these findings support a distinctive genetic basis for pediatric-type follicular lymphoma and pediatric marginal zone lymphoma when compared with adult subtypes and to one another. Moreover, identification of a recurrent point mutation in IRF8 provides insight into a potential driver mutation in the pathogenesis of pediatric-type follicular lymphoma with implications for novel diagnostic or therapeutic strategies.