Molecular Architecture of the Mouse Nervous System.

The mammalian nervous system executes complex behaviors controlled by specialized, precisely positioned, and interacting cell types. Here, we used RNA sequencing of half a million single cells to create a detailed census of cell types in the mouse nervous system. We mapped cell types spatially and derived a hierarchical, data-driven taxonomy. Neurons were the most diverse and were grouped by developmental anatomical units and by the expression of neurotransmitters and neuropeptides. Neuronal diversity was driven by genes encoding cell identity, synaptic connectivity, neurotransmission, and membrane conductance. We discovered seven distinct, regionally restricted astrocyte types that obeyed developmental boundaries and correlated with the spatial distribution of key glutamate and glycine neurotransmitters. In contrast, oligodendrocytes showed a loss of regional identity followed by a secondary diversification. The resource presented here lays a solid foundation for understanding the molecular architecture of the mammalian nervous system and enables genetic manipulation of specific cell types.

Origin, fate and dynamics of macrophages at central nervous system interfaces.

Perivascular, subdural meningeal and choroid plexus macrophages are non-parenchymal macrophages that mediate immune responses at brain boundaries. Although the origin of parenchymal microglia has recently been elucidated, much less is known about the precursors, the underlying transcriptional program and the dynamics of the other macrophages in the central nervous system (CNS). It was assumed that they have a high turnover from blood-borne monocytes. However, using parabiosis and fate-mapping approaches in mice, we found that CNS macrophages arose from hematopoietic precursors during embryonic development and established stable populations, with the notable exception of choroid plexus macrophages, which had dual origins and a shorter life span. The generation of CNS macrophages relied on the transcription factor PU.1, whereas the MYB, BATF3 and NR4A1 transcription factors were not required.

Single-Cell Profiling Uncovers Evolutionary Divergence of Hypocretin/Orexin Neuronal Subpopulations.

Brain nuclei are traditionally defined by their anatomy, activity, and expression of specific markers. The hypothalamus contains discrete neuronal populations that coordinate fundamental behavioral functions, including sleep and wakefulness, in all vertebrates. Particularly, the diverse roles of hypocretin/orexin (Hcrt)-releasing neurons suggest functional heterogeneity among Hcrt neurons. Using single-cell RNA sequencing (scRNA-seq) and high-resolution imaging of the adult male and female zebrafish hypothalamic periventricular zone, we identified 21 glutamatergic and 28 GABAergic cell types. Integration of zebrafish and mouse scRNA-seq revealed evolutionary conserved and divergent hypothalamic cell types. The expression of specific genes, including npvf, which encodes a sleep-regulating neuropeptide, was enriched in subsets of glutamatergic Hcrt neurons in both larval and adult zebrafish. The genetic profile, activity, and neurite processing of the neuronal subpopulation that coexpresses both Hcrt and Npvf (Hcrt+Npvf+) differ from other Hcrt neurons. These interspecies findings provide a unified annotation of hypothalamic cell types and suggest that the heterogeneity of Hcrt neurons enables multifunctionality, such as consolidation of both wake and sleep by the Hcrt- and Npvf-releasing neuronal subpopulation.

RNA velocity of single cells.

RNA abundance is a powerful indicator of the state of individual cells. Single-cell RNA sequencing can reveal RNA abundance with high quantitative accuracy, sensitivity and throughput1. However, this approach captures only a static snapshot at a point in time, posing a challenge for the analysis of time-resolved phenomena such as embryogenesis or tissue regeneration. Here we show that RNA velocity-the time derivative of the gene expression state-can be directly estimated by distinguishing between unspliced and spliced mRNAs in common single-cell RNA sequencing protocols. RNA velocity is a high-dimensional vector that predicts the future state of individual cells on a timescale of hours. We validate its accuracy in the neural crest lineage, demonstrate its use on multiple published datasets and technical platforms, reveal the branching lineage tree of the developing mouse hippocampus, and examine the kinetics of transcription in human embryonic brain. We expect RNA velocity to greatly aid the analysis of developmental lineages and cellular dynamics, particularly in humans.

Classes and continua of hippocampal CA1 inhibitory neurons revealed by single-cell transcriptomics.

Understanding any brain circuit will require a categorization of its constituent neurons. In hippocampal area CA1, at least 23 classes of GABAergic neuron have been proposed to date. However, this list may be incomplete; additionally, it is unclear whether discrete classes are sufficient to describe the diversity of cortical inhibitory neurons or whether continuous modes of variability are also required. We studied the transcriptomes of 3,663 CA1 inhibitory cells, revealing 10 major GABAergic groups that divided into 49 fine-scale clusters. All previously described and several novel cell classes were identified, with three previously described classes unexpectedly found to be identical. A division into discrete classes, however, was not sufficient to describe the diversity of these cells, as continuous variation also occurred between and within classes. Latent factor analysis revealed that a single continuous variable could predict the expression levels of several genes, which correlated similarly with it across multiple cell types. Analysis of the genes correlating with this variable suggested it reflects a range from metabolically highly active faster-spiking cells that proximally target pyramidal cells to slower-spiking cells targeting distal dendrites or interneurons. These results elucidate the complexity of inhibitory neurons in one of the simplest cortical structures and show that characterizing these cells requires continuous modes of variation as well as discrete cell classes.

Sexual dimorphism in synaptic inputs to the mouse amygdala and orbital cortex.

The medial amygdala (MeA) is a sexually dimorphic brain region that regulates fear responses, emotional memories, and social behaviors. It is known to be larger and contains more cells in males. The MeA integrates information through input connections from olfactory regions, bed nucleus of the stria terminalis, ventral hippocampus, and thalamic and hypothalamic structures. We hypothesize that in addition to the size differences, there are differences in regional connectivity between the sexes. In this study, we utilized G-deleted rabies monosynaptic retrograde tracing to compare amygdala presynaptic cells in male and female whole mouse brains. We report differences in connection patterns to the amygdala, with higher overall connectivity (presynaptic per starter) in males and a larger fraction of inputs originating from the bed nucleus of the stria terminalis, lateral septum, and medial preoptic area. Furthermore, we examined input connections to the orbital cortex (ORB), a brain region shown to be larger in volume in females, and found the opposite trend, where females had more total inputs. Together, our findings extend the evidence for sexual dimorphism in the brain to the neuronal wiring pattern, with likely impacts on behavior and disease susceptibility.

Sex, strain, and lateral differences in brain cytoarchitecture across a large mouse population.

The mouse brain is by far the most intensively studied among mammalian brains, yet basic measures of its cytoarchitecture remain obscure. For example, quantifying cell numbers, and the interplay of sex, strain, and individual variability in cell density and volume is out of reach for many regions. The Allen Mouse Brain Connectivity project produces high-resolution full brain images of hundreds of brains. Although these were created for a different purpose, they reveal details of neuroanatomy and cytoarchitecture. Here, we used this population to systematically characterize cell density and volume for each anatomical unit in the mouse brain. We developed a DNN-based segmentation pipeline that uses the autofluorescence intensities of images to segment cell nuclei even within the densest regions, such as the dentate gyrus. We applied our pipeline to 507 brains of males and females from C57BL/6J and FVB.CD1 strains. Globally, we found that increased overall brain volume does not result in uniform expansion across all regions. Moreover, region-specific density changes are often negatively correlated with the volume of the region; therefore, cell count does not scale linearly with volume. Many regions, including layer 2/3 across several cortical areas, showed distinct lateral bias. We identified strain-specific or sex-specific differences. For example, males tended to have more cells in extended amygdala and hypothalamic regions (MEA, BST, BLA, BMA, and LPO, AHN) while females had more cells in the orbital cortex (ORB). Yet, inter-individual variability was always greater than the effect size of a single qualifier. We provide the results of this analysis as an accessible resource for the community.

Two distinct molecular faces of preeclampsia revealed by single-cell transcriptomics.

INTRODUCTION: Preeclampsia is a multisystemic, pregnancy-specific disorder united by new-onset hypertension but with considerable variation in clinical manifestation, onset, and severity. For symptoms to regress, delivery of the placenta is required. For symptoms to regress, delivery of the placenta is required, making the placenta central to preeclampsia pathophysiology. To dissect which placental functions were impacted in two forms of preeclampsia, we studied molecular changes across the cell types of the placenta. METHODS: We performed a transcriptomic survey of single-cells and single-nuclei on cases of early- and late-onset preeclampsia with gestation-matched controls. FINDINGS: Our data revealed massive dysregulation of gene expression in all cell classes that was almost exclusive to early preeclampsia. For example, an important known receptor/ligand imbalance hallmarking angiogenic disfunction, sFLT1/placental growth factor (PGF), was reflected in striking, cell-autonomous dysregulation of FLT1 and PGF transcription in the syncytium in early preeclampsia only. Stromal cells and vasculature echoed an inflamed, stressed, anti-angiogenic environment. Finally, the placental immune niche set the tone for inflammation in early but not late preeclampsia. Here, fetal-origin Hofbauer and maternal-origin TREM2 macrophages were revealed as surprising main actors, while local cells of the adaptive immune system were largely unaffected. Late preeclampsia showed minimal cellular impact on the placenta. CONCLUSIONS: Our survey provides systematic molecular evidence for two distinct diseases. We resolved systematic molecular dysregulation to individual cell types with strong implications for definition, early detection, diagnosis, and treatment. FUNDING: Funded by the Preeclampsia Foundation through the Peter Joseph Pappas Research Grant.

A telencephalon cell type atlas for goldfish reveals diversity in the evolution of spatial structure and cell types.

Teleost fish form the largest group of vertebrates and show a tremendous variety of adaptive behaviors, making them critically important for the study of brain evolution and cognition. The neural basis mediating these behaviors remains elusive. We performed a systematic comparative survey of the goldfish telencephalon. We mapped cell types using single-cell RNA sequencing and spatial transcriptomics, resulting in de novo molecular neuroanatomy parcellation. Glial cells were highly conserved across 450 million years of evolution separating mouse and goldfish, while neurons showed diversity and modularity in gene expression. Specifically, somatostatin interneurons, famously interspersed in the mammalian isocortex for local inhibitory input, were curiously aggregated in a single goldfish telencephalon nucleus but molecularly conserved. Cerebral nuclei including the striatum, a hub for motivated behavior in amniotes, had molecularly conserved goldfish homologs. We suggest elements of a hippocampal formation across the goldfish pallium. Last, aiding study of the teleostan everted telencephalon, we describe substantial molecular similarities between goldfish and zebrafish neuronal taxonomies.

Neuronal types in the mouse amygdala and their transcriptional response to fear conditioning.

The amygdala is a brain region primarily associated with emotional response. The use of genetic markers and single-cell transcriptomics can provide insights into behavior-associated cell state changes. Here we present a detailed cell-type taxonomy of the adult mouse amygdala during fear learning and memory consolidation. We perform single-cell RNA sequencing on naïve and fear-conditioned mice, identify 130 neuronal cell types and validate their spatial distributions. A subset of all neuronal types is transcriptionally responsive to fear learning and memory retrieval. The activated engram cells upregulate activity-response genes and coordinate the expression of genes associated with neurite outgrowth, synaptic signaling, plasticity and development. We identify known and previously undescribed candidate genes responsive to fear learning. Our molecular atlas may be used to generate hypotheses to unveil the neuron types and neural circuits regulating the emotional component of learning and memory.

Transcriptional maintenance of cortical somatostatin interneuron subtype identity during migration.

Although cardinal cortical interneuron identity is established upon cell-cycle exit, it remains unclear whether specific interneuron subtypes are pre-established, and if so, how their identity is maintained prior to circuit integration. We conditionally removed Sox6 (Sox6-cKO) in migrating somatostatin (Sst+) interneurons and assessed the effects on their mature identity. In adolescent mice, five of eight molecular Sst+ subtypes were nearly absent in the Sox6-cKO cortex without a reduction in cell number. Sox6-cKO cells displayed electrophysiological maturity and expressed genes enriched within the broad class of Sst+ interneurons. Furthermore, we could infer subtype identity prior to cortical integration (embryonic day 18.5), suggesting that the loss in subtype was due to disrupted subtype maintenance. Conversely, Sox6 removal at postnatal day 7 did not disrupt marker expression in the mature cortex. Therefore, Sox6 is necessary during migration for maintenance of Sst+ subtype identity, indicating that subtype maintenance requires active transcriptional programs.

Conserved properties of dentate gyrus neurogenesis across postnatal development revealed by single-cell RNA sequencing.

The dentate gyrus of the hippocampus is a brain region in which neurogenesis persists into adulthood; however, the relationship between developmental and adult dentate gyrus neurogenesis has not been examined in detail. Here we used single-cell RNA sequencing to reveal the molecular dynamics and diversity of dentate gyrus cell types in perinatal, juvenile, and adult mice. We found distinct quiescent and proliferating progenitor cell types, linked by transient intermediate states to neuroblast stages and fully mature granule cells. We observed shifts in the molecular identity of quiescent and proliferating radial glia and granule cells during the postnatal period that were then maintained through adult stages. In contrast, intermediate progenitor cells, neuroblasts, and immature granule cells were nearly indistinguishable at all ages. These findings demonstrate the fundamental similarity of postnatal and adult neurogenesis in the hippocampus and pinpoint the early postnatal transformation of radial glia from embryonic progenitors to adult quiescent stem cells.

Neuronal atlas of the dorsal horn defines its architecture and links sensory input to transcriptional cell types.

The dorsal horn of the spinal cord is critical to processing distinct modalities of noxious and innocuous sensation, but little is known of the neuronal subtypes involved, hampering efforts to deduce principles governing somatic sensation. Here we used single-cell RNA sequencing to classify sensory neurons in the mouse dorsal horn. We identified 15 inhibitory and 15 excitatory molecular subtypes of neurons, equaling the complexity in cerebral cortex. Validating our classification scheme in vivo and matching cell types to anatomy of the dorsal horn by spatial transcriptomics reveals laminar enrichment for each of the cell types. Neuron types, when combined, define a multilayered organization with like neurons layered together. Employing our scheme, we find that heat and cold stimuli activate discrete sets of both excitatory and inhibitory neuron types. This work provides a systematic and comprehensive molecular classification of spinal cord sensory neurons, enabling functional interrogation of sensory processing.

STRT-seq-2i: dual-index 5' single cell and nucleus RNA-seq on an addressable microwell array.

Single-cell RNA-seq has become routine for discovering cell types and revealing cellular diversity, but archived human brain samples still pose a challenge to current high-throughput platforms. We present STRT-seq-2i, an addressable 9600-microwell array platform, combining sampling by limiting dilution or FACS, with imaging and high throughput at competitive cost. We applied the platform to fresh single mouse cortical cells and to frozen post-mortem human cortical nuclei, matching the performance of a previous lower-throughput platform while retaining a high degree of flexibility, potentially also for other high-throughput applications.

Visceral motor neuron diversity delineates a cellular basis for nipple- and pilo-erection muscle control.

Despite the variety of physiological and target-related functions, little is known regarding the cellular complexity in the sympathetic ganglion. We explored the heterogeneity of mouse stellate and thoracic ganglia and found an unexpected variety of cell types. We identified specialized populations of nipple- and pilo-erector muscle neurons. These neurons extended axonal projections and were born among other neurons during embryogenesis, but remained unspecialized until target organogenesis occurred postnatally. Target innervation and cell-type specification was coordinated by an intricate acquisition of unique combinations of growth factor receptors and the initiation of expression of concomitant ligands by the nascent erector muscles. Overall, our results provide compelling evidence for a highly sophisticated organization of the sympathetic nervous system into discrete outflow channels that project to well-defined target tissues and offer mechanistic insight into how diversity and connectivity are established during development.

Oligodendrocyte heterogeneity in the mouse juvenile and adult central nervous system.

Oligodendrocytes have been considered as a functionally homogeneous population in the central nervous system (CNS). We performed single-cell RNA sequencing on 5072 cells of the oligodendrocyte lineage from 10 regions of the mouse juvenile and adult CNS. Thirteen distinct populations were identified, 12 of which represent a continuum from Pdgfra(+) oligodendrocyte precursor cells (OPCs) to distinct mature oligodendrocytes. Initial stages of differentiation were similar across the juvenile CNS, whereas subsets of mature oligodendrocytes were enriched in specific regions in the adult brain. Newly formed oligodendrocytes were detected in the adult CNS and were responsive to complex motor learning. A second Pdgfra(+) population, distinct from OPCs, was found along vessels. Our study reveals the dynamics of oligodendrocyte differentiation and maturation, uncoupling them at a transcriptional level and highlighting oligodendrocyte heterogeneity in the CNS.

Base preferences in non-templated nucleotide incorporation by MMLV-derived reverse transcriptases.

Reverse transcriptases derived from Moloney Murine Leukemia Virus (MMLV) have an intrinsic terminal transferase activity, which causes the addition of a few non-templated nucleotides at the 3' end of cDNA, with a preference for cytosine. This mechanism can be exploited to make the reverse transcriptase switch template from the RNA molecule to a secondary oligonucleotide during first-strand cDNA synthesis, and thereby to introduce arbitrary barcode or adaptor sequences in the cDNA. Because the mechanism is relatively efficient and occurs in a single reaction, it has recently found use in several protocols for single-cell RNA sequencing. However, the base preference of the terminal transferase activity is not known in detail, which may lead to inefficiencies in template switching when starting from tiny amounts of mRNA. Here, we used fully degenerate oligos to determine the exact base preference at the template switching site up to a distance of ten nucleotides. We found a strong preference for guanosine at the first non-templated nucleotide, with a greatly reduced bias at progressively more distant positions. Based on this result, and a number of careful optimizations, we report conditions for efficient template switching for cDNA amplification from single cells.