Search

We analyzed >700,000 single-nucleus RNA sequencing profiles from 106 donors during prenatal and postnatal developmental stages and identified lineage-specific programs that underlie the development of specific subtypes of excitatory cortical neurons, interneurons, glial cell types, and brain vasculature. By leveraging single-nucleus chromatin accessibility data, we delineated enhancer gene regulatory networks and transcription factors that control commitment of specific cortical lineages. By intersecting our results with genetic risk factors for human brain diseases, we identified the cortical cell types and lineages most vulnerable to genetic insults of different brain disorders, especially autism. We find that lineage-specific gene expression programs up-regulated in female cells are especially enriched for the genetic risk factors of autism. Our study captures the molecular progression of cortical lineages across human development.

The molecular mechanisms and evolutionary changes accompanying synapse development are still poorly understood1,2. Here we generate a cross-species proteomic map of synapse development in the human, macaque and mouse neocortex. By tracking the changes of more than 1,000 postsynaptic density (PSD) proteins from midgestation to young adulthood, we find that PSD maturation in humans separates into three major phases that are dominated by distinct pathways. Cross-species comparisons reveal that human PSDs mature about two to three times slower than those of other species and contain higher levels of Rho guanine nucleotide exchange factors (RhoGEFs) in the perinatal period. Enhancement of RhoGEF signalling in human neurons delays morphological maturation of dendritic spines and functional maturation of synapses, potentially contributing to the neotenic traits of human brain development. In addition, PSD proteins can be divided into four modules that exert stage- and cell-type-specific functions, possibly explaining their differential associations with cognitive functions and diseases. Our proteomic map of synapse development provides a blueprint for studying the molecular basis and evolutionary changes of synapse maturation.

Chronic neurodegeneration and acute injuries lead to neuron losses via diverse processes. We compared retinal ganglion cell (RGC) responses between chronic glaucomatous conditions and the acute injury model. Among major RGC subclasses, αRGCs and intrinsically photosensitive RGCs (ipRGCs) preferentially survive glaucomatous conditions, similar to findings in the retina subject to axotomy. Focusing on an αRGC intrinsic factor, Osteopontin (secreted phosphoprotein 1 [Spp1]), we found an ectopic neuronal expression of Osteopontin (Spp1) in other RGCs subject to glaucomatous conditions. This contrasted with the Spp1 downregulation subject to axotomy. αRGC-specific Spp1 elimination led to significant αRGC loss, diminishing their resiliency. Spp1 overexpression led to robust neuroprotection of susceptible RGC subclasses under glaucomatous conditions. In contrast, Spp1 overexpression did not significantly protect RGCs subject to axotomy. Additionally, SPP1 marked adult human RGC subsets with large somata and SPP1 expression in the aqueous humor correlated with glaucoma severity. Our study reveals Spp1's role in mediating neuronal resiliency in glaucoma.

Radial glial cells, the stem cells of the cerebral cortex, extend a long basal fiber that ends in basal endfeet. A new study in PLOS Biology found that non-muscle myosins control basal endfoot integrity to regulate interneuron organization.

Zika virus (ZIKV) is a flavivirus transmitted via mosquitoes and sex to cause congenital neurodevelopmental defects, including microcephaly. Inherited forms of microcephaly (MCPH) are associated with disrupted centrosome organization. Similarly, we found that ZIKV infection disrupted centrosome organization. ZIKV infection disrupted the organization of centrosomal proteins including CEP63, a MCPH-associated protein. The ZIKV nonstructural protein NS3 bound CEP63, and expression of NS3 was sufficient to alter centrosome architecture and CEP63 localization. Loss of CEP63 suppressed ZIKV-induced centrosome disorganization, indicating that ZIKV requires CEP63 to disrupt centrosome organization. ZIKV infection or CEP63 loss decreased the centrosomal localization and stability of TANK-binding kinase 1 (TBK1), a regulator of the innate immune response. ZIKV infection also increased the centrosomal accumulation of the CEP63 interactor DTX4, a ubiquitin ligase that degrades TBK1. Therefore, we propose that ZIKV disrupts CEP63 function to increase centrosomal DTX4 localization and destabilization of TBK1, thereby tempering the innate immune response.

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 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) readily infects a variety of cell types impacting the function of vital organ systems, with particularly severe impact on respiratory function. Neurological symptoms, which range in severity, accompany as many as one-third of COVID-19 cases, indicating a potential vulnerability of neural cell types. To assess whether human cortical cells can be directly infected by SARS-CoV-2, we utilized stem-cell-derived cortical organoids as well as primary human cortical tissue, both from developmental and adult stages. We find significant and predominant infection in cortical astrocytes in both primary tissue and organoid cultures, with minimal infection of other cortical populations. Infected and bystander astrocytes have a corresponding increase in inflammatory gene expression, reactivity characteristics, increased cytokine and growth factor signaling, and cellular stress. Although human cortical cells, particularly astrocytes, have no observable ACE2 expression, we find high levels of coronavirus coreceptors in infected astrocytes, including CD147 and DPP4. Decreasing coreceptor abundance and activity reduces overall infection rate, and increasing expression is sufficient to promote infection. Thus, we find tropism of SARS-CoV-2 for human astrocytes resulting in inflammatory gliosis-type injury that is dependent on coronavirus coreceptors.

The human cortex contains inhibitory interneurons derived from the medial ganglionic eminence (MGE), a germinal zone in the embryonic ventral forebrain. How this germinal zone generates sufficient interneurons for the human brain remains unclear. We found that the human MGE (hMGE) contains nests of proliferative neuroblasts with ultrastructural and transcriptomic features that distinguish them from other progenitors in the hMGE. When dissociated hMGE cells are transplanted into the neonatal mouse brain, they reform into nests containing proliferating neuroblasts that generate young neurons that migrate extensively into the mouse forebrain and mature into different subtypes of functional interneurons. Together, these results indicate that the nest organization and sustained proliferation of neuroblasts in the hMGE provide a mechanism for the extended production of interneurons for the human forebrain.

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.

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.

Nuclear compartments are thought to play a role in three-dimensional genome organization and gene expression. In mammalian brain, the architecture and dynamics of nuclear compartment-associated genome organization is not known. In this study, we developed Genome Organization using CUT and RUN Technology (GO-CaRT) to map genomic interactions with two nuclear compartments-the nuclear lamina and nuclear speckles-from different regions of the developing mouse, macaque and human brain. Lamina-associated domain (LAD) architecture in cells in vivo is distinct from that of cultured cells, including major differences in LADs previously considered to be cell type invariant. In the mouse and human forebrain, dorsal and ventral neural precursor cells have differences in LAD architecture that correspond to their regional identity. LADs in the human and mouse cortex contain transcriptionally highly active sub-domains characterized by broad depletion of histone-3-lysine-9 dimethylation. Evolutionarily conserved LADs in human, macaque and mouse brain are enriched for transcriptionally active neural genes associated with synapse function. By integrating GO-CaRT maps with genome-wide association study data, we found speckle-associated domains to be enriched for schizophrenia risk loci, indicating a physical relationship between these disease-associated genetic variants and a specific nuclear structure. Our work provides a framework for understanding the relationship between distinct nuclear compartments and genome function in brain development and disease.

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.

Aberrant aggregation of the RNA-binding protein TDP-43 in neurons is a hallmark of frontotemporal lobar degeneration caused by haploinsufficiency in the gene encoding progranulin1,2. However, the mechanism leading to TDP-43 proteinopathy remains unclear. Here we use single-nucleus RNA sequencing to show that progranulin deficiency promotes microglial transition from a homeostatic to a disease-specific state that causes endolysosomal dysfunction and neurodegeneration in mice. These defects persist even when Grn-/- microglia are cultured ex vivo. In addition, single-nucleus RNA sequencing reveals selective loss of excitatory neurons at disease end-stage, which is characterized by prominent nuclear and cytoplasmic TDP-43 granules and nuclear pore defects. Remarkably, conditioned media from Grn-/- microglia are sufficient to promote TDP-43 granule formation, nuclear pore defects and cell death in excitatory neurons via the complement activation pathway. Consistent with these results, deletion of the genes encoding C1qa and C3 mitigates microglial toxicity and rescues TDP-43 proteinopathy and neurodegeneration. These results uncover previously unappreciated contributions of chronic microglial toxicity to TDP-43 proteinopathy during neurodegeneration.

BackgroundStudies of individuals with autism spectrum disorder (ASD) have revealed a strong multigenic basis with the identification of hundreds of ASD susceptibility genes. ASD is characterized by social deficits and a range of other phenotypes, implicating complex genetics and involvement of a variety of brain regions. However, how mutations and mis-expression of select gene sets are associated with the behavioral components of ASD remains unknown. We reasoned that for genes to be associated with ASD core behaviors they must be: (1) expressed in brain regions relevant to ASD social behaviors and (2) expressed during the ASD susceptible window of brain development.MethodsFocusing on the amygdala, a brain region whose dysfunction has been highly implicated in the social component of ASD, we mined publicly available gene expression databases to identify ASD-susceptibility genes expressed during human and mouse amygdala development. We found that a large cohort of known ASD susceptibility genes is expressed in the developing human and mouse amygdala. We further performed analysis of single-nucleus RNA-seq (snRNA-seq) data from microdissected amygdala tissue from five ASD and five control human postmortem brains ranging in age from 4 to 20 years to elucidate cell type specificity of amygdala-expressed genes and their dysregulation in ASD.ResultsOur analyses revealed that of the high-ranking ASD susceptibility genes, 80 are expressed in both human and mouse amygdala during fetal to early postnatal stages of development. Our human snRNA-seq analyses revealed cohorts of genes with altered expression in the ASD amygdala postnatally, especially within excitatory neurons, with dysregulated expression of seven genes predicted from our datamining pipeline.LimitationsWe were limited by the ages for which we were able to obtain human tissue; therefore, the results from our datamining pipeline approach will require validation, to the extent possible, in human tissue from earlier developmental stages.ConclusionsOur pipeline narrows down the number of amygdala-expressed genes possibly involved in the social pathophysiology of ASD. Our human single-nucleus gene expression analyses revealed that ASD is characterized by changes in gene expression in specific cell types in the early postnatal amygdala.

Lineage-specific epigenomic changes during human corticogenesis have been difficult to study owing to challenges with sample availability and tissue heterogeneity. For example, previous studies using single-cell RNA sequencing identified at least 9 major cell types and up to 26 distinct subtypes in the dorsal cortex alone1,2. Here we characterize cell-type-specific cis-regulatory chromatin interactions, open chromatin peaks, and transcriptomes for radial glia, intermediate progenitor cells, excitatory neurons, and interneurons isolated from mid-gestational samples of the human cortex. We show that chromatin interactions underlie several aspects of gene regulation, with transposable elements and disease-associated variants enriched at distal interacting regions in a cell-type-specific manner. In addition, promoters with increased levels of chromatin interactivity-termed super-interactive promoters-are enriched for lineage-specific genes, suggesting that interactions at these loci contribute to the fine-tuning of transcription. Finally, we develop CRISPRview, a technique that integrates immunostaining, CRISPR interference, RNAscope, and image analysis to validate cell-type-specific cis-regulatory elements in heterogeneous populations of primary cells. Our findings provide insights into cell-type-specific gene expression patterns in the developing human cortex and advance our understanding of gene regulation and lineage specification during this crucial developmental window.

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.

Outer radial glial (oRG) cells are a population of neural stem cells prevalent in the developing human cortex that contribute to its cellular diversity and evolutionary expansion. The mammalian Target of Rapamycin (mTOR) signaling pathway is active in human oRG cells. Mutations in mTOR pathway genes are linked to a variety of neurodevelopmental disorders and malformations of cortical development. We find that dysregulation of mTOR signaling specifically affects oRG cells, but not other progenitor types, by changing the actin cytoskeleton through the activity of the Rho-GTPase, CDC42. These effects change oRG cellular morphology, migration, and mitotic behavior, but do not affect proliferation or cell fate. Thus, mTOR signaling can regulate the architecture of the developing human cortex by maintaining the cytoskeletal organization of oRG cells and the radial glia scaffold. Our study provides insight into how mTOR dysregulation may contribute to neurodevelopmental disease.

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.

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.

To discover regulatory elements driving the specificity of gene expression in different cell types and regions of the developing human brain, we generated an atlas of open chromatin from nine dissected regions of the mid-gestation human telencephalon, as well as microdissected upper and deep layers of the prefrontal cortex. We identified a subset of open chromatin regions (OCRs), termed predicted regulatory elements (pREs), that are likely to function as developmental brain enhancers. pREs showed temporal, regional, and laminar differences in chromatin accessibility and were correlated with gene expression differences across regions and gestational ages. We identified two functional de novo variants in a pRE for autism risk gene SLC6A1, and using CRISPRa, demonstrated that this pRE regulates SCL6A1. Additionally, mouse transgenic experiments validated enhancer activity for pREs proximal to FEZF2 and BCL11A. Thus, this atlas serves as a resource for decoding neurodevelopmental gene regulation in health and disease.

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.

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.

Accurate RNA quantification at the single-cell level is critical for understanding the dynamics of gene expression and regulation across space and time. Single molecule FISH (smFISH), such as RNAscope, provides spatial and quantitative measurements of individual transcripts, therefore, can be used to explore differential gene expression among a heterogeneous cell population if combined with cell identify information. However, such analysis is not straightforward, and existing image analysis pipelines cannot integrate both RNA transcripts and cellular staining information to automatically output cell type-specific gene expression. We developed an efficient and customizable analysis method, Single-Molecule Automatic RNA Transcription Quantification (SMART-Q), to enable the analysis of gene transcripts in a cell type-specific manner. SMART-Q efficiently infers cell identity information from multiplexed immuno-staining and quantifies cell type-specific transcripts using a 3D Gaussian fitting algorithm. Furthermore, we have optimized SMART-Q for user experiences, such as flexible parameters specification, batch data outputs, and visualization of analysis results. SMART-Q meets the demands for efficient quantification of single-molecule RNA and can be widely used for cell type-specific RNA transcript analysis.

Pelizaeus-Merzbacher disease (PMD) is an X-linked leukodystrophy caused by mutations in Proteolipid Protein 1 (PLP1), encoding a major myelin protein, resulting in profound developmental delay and early lethality. Previous work showed involvement of unfolded protein response (UPR) and endoplasmic reticulum (ER) stress pathways, but poor PLP1 genotype-phenotype associations suggest additional pathogenetic mechanisms. Using induced pluripotent stem cell (iPSC) and gene-correction, we show that patient-derived oligodendrocytes can develop to the pre-myelinating stage, but subsequently undergo cell death. Mutant oligodendrocytes demonstrated key hallmarks of ferroptosis including lipid peroxidation, abnormal iron metabolism, and hypersensitivity to free iron. Iron chelation rescued mutant oligodendrocyte apoptosis, survival, and differentiationin vitro, and post-transplantation in vivo. Finally, systemic treatment of Plp1 mutant Jimpy mice with deferiprone, a small molecule iron chelator, reduced oligodendrocyte apoptosis and enabled myelin formation. Thus, oligodendrocyte iron-induced cell death and myelination is rescued by iron chelation in PMD pre-clinical models.

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.

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 human amygdala grows during childhood, and its abnormal development is linked to mood disorders. The primate amygdala contains a large population of immature neurons in the paralaminar nuclei (PL), suggesting protracted development and possibly neurogenesis. Here we studied human PL development from embryonic stages to adulthood. The PL develops next to the caudal ganglionic eminence, which generates inhibitory interneurons, yet most PL neurons express excitatory markers. In children, most PL cells are immature (DCX+PSA-NCAM+), and during adolescence many transition into mature (TBR1+VGLUT2+) neurons. Immature PL neurons persist into old age, yet local progenitor proliferation sharply decreases in infants. Using single nuclei RNA sequencing, we identify the transcriptional profile of immature excitatory neurons in the human amygdala between 4-15 years. We conclude that the human PL contains excitatory neurons that remain immature for decades, a possible substrate for persistent plasticity at the interface of the hippocampus and amygdala.

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.

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 (

MicroRNAs (miRNAs) regulate many cellular events during brain development by interacting with hundreds of mRNA transcripts. However, miRNAs operate nonuniformly upon the transcriptional profile with an as yet unknown logic. Shortcomings in defining miRNA-mRNA networks include limited knowledge of in vivo miRNA targets and their abundance in single cells. By combining multiple complementary approaches, high-throughput sequencing of RNA isolated by cross-linking immunoprecipitation with an antibody to AGO2 (AGO2-HITS-CLIP), single-cell profiling and computational analyses using bipartite and coexpression networks, we show that miRNA-mRNA interactions operate as functional modules that often correspond to cell-type identities and undergo dynamic transitions during brain development. These networks are highly dynamic during development and over the course of evolution. One such interaction is between radial-glia-enriched ORC4 and miR-2115, a great-ape-specific miRNA, which appears to control radial glia proliferation rates during human brain development.

The neocortex of primates, including humans, contains more abundant and diverse inhibitory neurons compared with rodents, but the molecular foundations of these observations are unknown. Through integrative gene coexpression analysis, we determined a consensus transcriptional profile of GABAergic neurons in mid-gestation human neocortex. By comparing this profile to genes expressed in GABAergic neurons purified from neonatal mouse neocortex, we identified conserved and distinct aspects of gene expression in these cells between the species. We show here that the calcium-binding protein secretagogin (SCGN) is robustly expressed by neocortical GABAergic neurons derived from caudal ganglionic eminences (CGE) and lateral ganglionic eminences during human but not mouse brain development. Through electrophysiological and morphometric analyses, we examined the effects of SCGN expression on GABAergic neuron function and form. Forced expression of SCGN in CGE-derived mouse GABAergic neurons significantly increased total neurite length and arbor complexity following transplantation into mouse neocortex, revealing a molecular pathway that contributes to morphological differences in these cells between rodents and primates.

Genomic association studies of common or rare protein-coding variation have established robust statistical approaches to account for multiple testing. Here we present a comparable framework to evaluate rare and de novo noncoding single-nucleotide variants, insertion/deletions, and all classes of structural variation from whole-genome sequencing (WGS). Integrating genomic annotations at the level of nucleotides, genes, and regulatory regions, we define 51,801 annotation categories. Analyses of 519 autism spectrum disorder families did not identify association with any categories after correction for 4,123 effective tests. Without appropriate correction, biologically plausible associations are observed in both cases and controls. Despite excluding previously identified gene-disrupting mutations, coding regions still exhibited the strongest associations. Thus, in autism, the contribution of de novo noncoding variation is probably modest in comparison to that of de novo coding variants. Robust results from future WGS studies will require large cohorts and comprehensive analytical strategies that consider the substantial multiple-testing burden.

New neurons continue to be generated in the subgranular zone of the dentate gyrus of the adult mammalian hippocampus. This process has been linked to learning and memory, stress and exercise, and is thought to be altered in neurological disease. In humans, some studies have suggested that hundreds of new neurons are added to the adult dentate gyrus every day, whereas other studies find many fewer putative new neurons. Despite these discrepancies, it is generally believed that the adult human hippocampus continues to generate new neurons. Here we show that a defined population of progenitor cells does not coalesce in the subgranular zone during human fetal or postnatal development. We also find that the number of proliferating progenitors and young neurons in the dentate gyrus declines sharply during the first year of life and only a few isolated young neurons are observed by 7 and 13 years of age. In adult patients with epilepsy and healthy adults (18-77 years; n = 17 post-mortem samples from controls; n = 12 surgical resection samples from patients with epilepsy), young neurons were not detected in the dentate gyrus. In the monkey (Macaca mulatta) hippocampus, proliferation of neurons in the subgranular zone was found in early postnatal life, but this diminished during juvenile development as neurogenesis decreased. We conclude that recruitment of young neurons to the primate hippocampus decreases rapidly during the first years of life, and that neurogenesis in the dentate gyrus does not continue, or is extremely rare, in adult humans. The early decline in hippocampal neurogenesis raises questions about how the function of the dentate gyrus differs between humans and other species in which adult hippocampal neurogenesis is preserved.

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.

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.

A comprehensive characterization of neuronal cell types, their distributions, and patterns of connectivity is critical for understanding the properties of neural circuits and how they generate behaviors. Here we review the experiences of the BRAIN Initiative Cell Census Consortium, ten pilot projects funded by the U.S. BRAIN Initiative, in developing, validating, and scaling up emerging genomic and anatomical mapping technologies for creating a complete inventory of neuronal cell types and their connections in multiple species and during development. These projects lay the foundation for a larger and longer-term effort to generate whole-brain cell atlases in species including mice and humans.

Understanding the development and dysfunction of the human brain is a major goal of neurobiology. Much of our current understanding of human brain development has been derived from the examination of post-mortem and pathological specimens, bolstered by observations of developing non-human primates and experimental studies focused largely on mouse models. However, these tissue specimens and model systems cannot fully capture the unique and dynamic features of human brain development. Recent advances in stem cell technologies that enable the generation of human brain organoids from pluripotent stem cells (PSCs) promise to profoundly change our understanding of the development of the human brain and enable a detailed study of the pathogenesis of inherited and acquired brain diseases.

Classical lissencephaly is a genetic neurological disorder associated with mental retardation and intractable epilepsy, and Miller-Dieker syndrome (MDS) is the most severe form of the disease. In this study, to investigate the effects of MDS on human progenitor subtypes that control neuronal output and influence brain topology, we analyzed cerebral organoids derived from control and MDS-induced pluripotent stem cells (iPSCs) using time-lapse imaging, immunostaining, and single-cell RNA sequencing. We saw a cell migration defect that was rescued when we corrected the MDS causative chromosomal deletion and severe apoptosis of the founder neuroepithelial stem cells, accompanied by increased horizontal cell divisions. We also identified a mitotic defect in outer radial glia, a progenitor subtype that is largely absent from lissencephalic rodents but critical for human neocortical expansion. Our study, therefore, deepens our understanding of MDS cellular pathogenesis and highlights the broad utility of cerebral organoids for modeling human neurodevelopmental disorders.

To understand how diverse progenitor cells contribute to human neocortex development, we examined forebrain progenitor behaviour using timelapse imaging. Here we find that cell cycle dynamics of human neuroepithelial (NE) cells differ from radial glial (RG) cells in both primary tissue and in stem cell-derived organoids. NE cells undergoing proliferative, symmetric divisions retract their basal processes, and both daughter cells regrow a new process following cytokinesis. The mitotic retraction of the basal process is recapitulated by NE cells in cerebral organoids generated from human-induced pluripotent stem cells. In contrast, RG cells undergoing vertical cleavage retain their basal fibres throughout mitosis, both in primary tissue and in older organoids. Our findings highlight developmentally regulated changes in mitotic behaviour that may relate to the role of RG cells to provide a stable scaffold for neuronal migration, and suggest that the transition in mitotic dynamics can be studied in organoid models.

The rapid spread of Zika virus (ZIKV) and its association with abnormal brain development constitute a global health emergency. Congenital ZIKV infection produces a range of mild to severe pathologies, including microcephaly. To understand the pathophysiology of ZIKV infection, we used models of the developing brain that faithfully recapitulate the tissue architecture in early to midgestation. We identify the brain cell populations that are most susceptible to ZIKV infection in primary human tissue, provide evidence for a mechanism of viral entry, and show that a commonly used antibiotic protects cultured brain cells by reducing viral proliferation. In the brain, ZIKV preferentially infected neural stem cells, astrocytes, oligodendrocyte precursor cells, and microglia, whereas neurons were less susceptible to infection. These findings suggest mechanisms for microcephaly and other pathologic features of infants with congenital ZIKV infection that are not explained by neural stem cell infection alone, such as calcifications in the cortical plate. Furthermore, we find that blocking the glia-enriched putative viral entry receptor AXL reduced ZIKV infection of astrocytes in vitro, and genetic knockdown of AXL in a glial cell line nearly abolished infection. Finally, we evaluate 2,177 compounds, focusing on drugs safe in pregnancy. We show that the macrolide antibiotic azithromycin reduced viral proliferation and virus-induced cytopathic effects in glial cell lines and human astrocytes. Our characterization of infection in the developing human brain clarifies the pathogenesis of congenital ZIKV infection and provides the basis for investigating possible therapeutic strategies to safely alleviate or prevent the most severe consequences of the epidemic.

The classic view of cortical development, embodied in the radial unit hypothesis, highlights the ventricular radial glia (vRG) scaffold as a key architectonic feature of the developing neocortex. The scaffold includes continuous fibers spanning the thickness of the developing cortex during neurogenesis across mammals. However, we find that in humans, the scaffold transforms into a physically discontinuous structure during the transition from infragranular to supragranular neuron production. As a consequence of this transformation, supragranular layer neurons arrive at their terminal positions in the cortical plate along outer radial glia (oRG) cell fibers. In parallel, the radial glia that contact the ventricle develop distinct gene expression profile and "truncated" morphology. We propose a supragranular layer expansion hypothesis that posits a deterministic role of oRG cells in the radial and tangential expansion of supragranular layers in primates, with implications for patterns of neuronal migration, area patterning, and cortical folding.

Long non-coding RNAs (lncRNAs) are a diverse and poorly conserved category of transcripts that have expanded greatly in primates, particularly in the brain. We identified an lncRNA, which has acquired 16 microRNA response elements for miR-143-3p in the Catarrhini branch of primates. This lncRNA, termed LncND (neurodevelopment), is expressed in neural progenitor cells and then declines in neurons. Binding and release of miR-143-3p by LncND control the expression of Notch receptors. LncND expression is enriched in radial glia cells (RGCs) in the ventricular and subventricular zones of developing human brain. Downregulation in neuroblastoma cells reduced cell proliferation and induced neuronal differentiation, an effect phenocopied by miR-143-3p overexpression. Gain of function of LncND in developing mouse cortex led to an expansion of PAX6+ RGCs. These findings support a role for LncND in miRNA-mediated regulation of Notch signaling within the neural progenitor pool in primates that may have contributed to the expansion of cerebral cortex.

The recent outbreak of Zika virus (ZIKV) in Brazil has been linked to substantial increases in fetal abnormalities and microcephaly. However, information about the underlying molecular and cellular mechanisms connecting viral infection to these defects remains limited. In this study we have examined the expression of receptors implicated in cell entry of several enveloped viruses including ZIKV across diverse cell types in the developing brain. Using single-cell RNA-seq and immunohistochemistry, we found that the candidate viral entry receptor AXL is highly expressed by human radial glial cells, astrocytes, endothelial cells, and microglia in developing human cortex and by progenitor cells in developing retina. We also show that AXL expression in radial glia is conserved in developing mouse and ferret cortex and in human stem cell-derived cerebral organoids, highlighting multiple experimental systems that could be applied to study mechanisms of ZIKV infectivity and effects on brain development.

γ-Aminobutyric acid (GABA) is the major inhibitory transmitter in the mature brain but is excitatory in the developing cortex. We found that mouse zona incerta (ZI) projection neurons form a GABAergic axon plexus in neonatal cortical layer 1, making synapses with neurons in both deep and superficial layers. A similar depolarizing GABAergic plexus exists in the developing human cortex. Selectively silencing mouse ZI GABAergic neurons at birth decreased synaptic activity and apical dendritic complexity of cortical neurons. The ZI GABAergic projection becomes inhibitory with maturation and can block epileptiform activity in the adult brain. These data reveal an early-developing GABAergic projection from the ZI to cortical layer 1 that is essential for proper development of cortical neurons and balances excitation with inhibition in the adult cortex.

During mammalian neocortical development, newborn excitatory and inhibitory neurons must migrate over long distances to reach their final positions within the cortical plate. In the lissencephalic rodent brain, pyramidal neurons are born in the ventricular and subventricular zones of the pallium and migrate along radial glia fibers to reach the appropriate cortical layer. Although much less is known about neuronal migration in species with a gyrencephalic cortex, retroviral studies in the ferret and primate suggest that, unlike the rodent, pyramidal neurons do not follow strict radial pathways and instead can disperse horizontally. However, the means by which pyramidal neurons laterally disperse remain unknown. In this study, we identified a viral labeling technique for visualizing neuronal migration in the ferret, a gyrencephalic carnivore, and found that migration was predominantly radial at early postnatal ages. In contrast, neurons displayed more tortuous migration routes with a decreased frequency of cortical plate-directed migration at later stages of neurogenesis concomitant with the start of brain folding. This was accompanied by neurons migrating sequentially along several different radial glial fibers, suggesting a mode by which pyramidal neurons may laterally disperse in a folded cortex. These findings provide insight into the migratory behavior of neurons in gyrencephalic species and provide a framework for using nonrodent model systems for studying neuronal migration disorders.Significance statementElucidating neuronal migration dynamics in the gyrencephalic, or folded, cortex is important for understanding neurodevelopmental disorders. Similar to the rodent, we found that neuronal migration was predominantly radial at early postnatal ages in the gyrencephalic ferret cortex. Interestingly, ferret neurons displayed more tortuous migration routes and a decreased frequency of radial migration at later ages coincident with the start of cortical folding. We found that ferret neurons use several different radial glial fibers as migratory guides, including those belonging to the recently described outer radial glia, suggesting a mechanism by which ferret neurons disperse laterally. It is likely that excitatory neurons horizontally disperse in other gyrencephalic mammals, including the primate, suggesting an important modification to the current model deduced primarily from the rodent.

The mammalian neocortex is composed of two major neuronal cell types with distinct origins: excitatory pyramidal neurons and inhibitory interneurons, generated in dorsal and ventral progenitor zones of the embryonic telencephalon, respectively. Thus, inhibitory neurons migrate relatively long distances to reach their destination in the developing forebrain. The role of lineage in the organization and circuitry of interneurons is still not well understood. Utilizing a combination of genetics, retroviral fate mapping, and lineage-specific retroviral barcode labeling, we find that clonally related interneurons can be widely dispersed while unrelated interneurons can be closely clustered. These data suggest that migratory mechanisms related to the clustering of interneurons occur largely independent of their clonal origin.