Search

Background: High 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., Results: To 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., Conclusions: Together, 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 (< 1%) may be identifiable only with much higher cell numbers.

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.

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

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.

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.

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.

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.

Neuroserpin is a serine-protease inhibitor mainly expressed in the CNS and involved in the inhibition of the proteolytic cascade. Animal models confirmed its neuroprotective role in perinatal hypoxia-ischaemia and adult stroke. Although neuroserpin may be a potential therapeutic target in the treatment of the aforementioned conditions, there is still no information in the literature on its distribution during human brain development. The present study provides a detailed description of the changing spatiotemporal patterns of neuroserpin focusing on physiological human brain development. Five stages were distinguished within our examined age range which spanned from the 7th gestational week until adulthood. In particular, subplate and deep cortical plate neurons were identified as the main sources of neuroserpin production between the 25th gestational week and the first postnatal month. Our immunohistochemical findings were substantiated by single cell RNA sequencing data showing specific neuronal and glial cell types expressing neuroserpin. The characterization of neuroserpin expression during physiological human brain development is essential for forthcoming studies which will explore its involvement in pathological conditions, such as perinatal hypoxia-ischaemia and adult stroke in human.

Adult cortical areas consist of specialized cell types and circuits that support unique higher-order cognitive functions. How this regional diversity develops from an initially uniform neuroepithelium has been the subject of decades of seminal research, and emerging technologies, including single-cell transcriptomics, provide a new perspective on area-specific molecular diversity. Here, we review the early developmental processes that underlie cortical arealization, including both cortex intrinsic and extrinsic mechanisms as embodied by the protomap and protocortex hypotheses, respectively. We propose an integrated model of serial homology whereby intrinsic genetic programs and local factors establish early transcriptomic differences between excitatory neurons destined to give rise to broad "proto-regions," and activity-dependent mechanisms lead to progressive refinement and formation of sharp boundaries between functional areas. Finally, we explore the potential of these basic developmental processes to inform our understanding of the emergence of functional neural networks and circuit abnormalities in neurodevelopmental disorders.

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.

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 diversity not captured by single-cell transcriptomics. We find that the serotonin receptor HTR2A selectively activates radial glia cells in the developing human, but not mouse, neocortex, and inhibiting HTR2A receptors in human radial glia disrupts the radial glial scaffold. We show highly specific neurotransmitter signaling during neurogenesis in the developing human neocortex and highlight evolutionarily divergent mechanisms of physiological signaling.

Direct comparisons of human and non-human primate brains can reveal molecular pathways underlying remarkable specializations of the human brain. However, chimpanzee tissue is inaccessible during neocortical neurogenesis when differences in brain size first appear. To identify human-specific features of cortical development, we leveraged recent innovations that permit generating pluripotent stem cell-derived cerebral organoids from chimpanzee. Despite metabolic differences, organoid models preserve gene regulatory networks related to primary cell types and developmental processes. We further identified 261 differentially expressed genes in human compared to both chimpanzee organoids and macaque cortex, enriched for recent gene duplications, and including multiple regulators of PI3K-AKT-mTOR signaling. We observed increased activation of this pathway in human radial glia, dependent on two receptors upregulated specifically in human: INSR and ITGB8. Our findings establish a platform for systematic analysis of molecular changes contributing to human brain development and evolution.

BACKGROUND:Tumor-associated macrophages (TAMs) are abundant in gliomas and immunosuppressive TAMs are a barrier to emerging immunotherapies. It is unknown to what extent macrophages derived from peripheral blood adopt the phenotype of brain-resident microglia in pre-treatment gliomas. The relative proportions of blood-derived macrophages and microglia have been poorly quantified in clinical samples due to a paucity of markers that distinguish these cell types in malignant tissue. RESULTS:We perform single-cell RNA-sequencing of human gliomas and identify phenotypic differences in TAMs of distinct lineages. We isolate TAMs from patient biopsies and compare them with macrophages from non-malignant human tissue, glioma atlases, and murine glioma models. We present a novel signature that distinguishes TAMs by ontogeny in human gliomas. Blood-derived TAMs upregulate immunosuppressive cytokines and show an altered metabolism compared to microglial TAMs. They are also enriched in perivascular and necrotic regions. The gene signature of blood-derived TAMs, but not microglial TAMs, correlates with significantly inferior survival in low-grade glioma. Surprisingly, TAMs frequently co-express canonical pro-inflammatory (M1) and alternatively activated (M2) genes in individual cells. CONCLUSIONS:We conclude that blood-derived TAMs significantly infiltrate pre-treatment gliomas, to a degree that varies by glioma subtype and tumor compartment. Blood-derived TAMs do not universally conform to the phenotype of microglia, but preferentially express immunosuppressive cytokines and show an altered metabolism. Our results argue against status quo therapeutic strategies that target TAMs indiscriminately and in favor of strategies that specifically target immunosuppressive blood-derived TAMs.