A taxonomy of transcriptomic cell types across the isocortex and hippocampal formation.

The isocortex and hippocampal formation (HPF) in the mammalian brain play critical roles in perception, cognition, emotion, and learning. We profiled ∼1.3 million cells covering the entire adult mouse isocortex and HPF and derived a transcriptomic cell-type taxonomy revealing a comprehensive repertoire of glutamatergic and GABAergic neuron types. Contrary to the traditional view of HPF as having a simpler cellular organization, we discover a complete set of glutamatergic types in HPF homologous to all major subclasses found in the six-layered isocortex, suggesting that HPF and the isocortex share a common circuit organization. We also identify large-scale continuous and graded variations of cell types along isocortical depth, across the isocortical sheet, and in multiple dimensions in hippocampus and subiculum. Overall, our study establishes a molecular architecture of the mammalian isocortex and hippocampal formation and begins to shed light on its underlying relationship with the development, evolution, connectivity, and function of these two brain structures.

Multimodal Analysis of Cell Types in a Hypothalamic Node Controlling Social Behavior.

The ventrolateral subdivision of the ventromedial hypothalamus (VMHvl) contains ∼4,000 neurons that project to multiple targets and control innate social behaviors including aggression and mounting. However, the number of cell types in VMHvl and their relationship to connectivity and behavioral function are unknown. We performed single-cell RNA sequencing using two independent platforms-SMART-seq (∼4,500 neurons) and 10x (∼78,000 neurons)-and investigated correspondence between transcriptomic identity and axonal projections or behavioral activation, respectively. Canonical correlation analysis (CCA) identified 17 transcriptomic types (T-types), including several sexually dimorphic clusters, the majority of which were validated by seqFISH. Immediate early gene analysis identified T-types exhibiting preferential responses to intruder males versus females but only rare examples of behavior-specific activation. Unexpectedly, many VMHvl T-types comprise a mixed population of neurons with different projection target preferences. Overall our analysis revealed that, surprisingly, few VMHvl T-types exhibit a clear correspondence with behavior-specific activation and connectivity.

A Suite of Transgenic Driver and Reporter Mouse Lines with Enhanced Brain-Cell-Type Targeting and Functionality.

Modern genetic approaches are powerful in providing access to diverse cell types in the brain and facilitating the study of their function. Here, we report a large set of driver and reporter transgenic mouse lines, including 23 new driver lines targeting a variety of cortical and subcortical cell populations and 26 new reporter lines expressing an array of molecular tools. In particular, we describe the TIGRE2.0 transgenic platform and introduce Cre-dependent reporter lines that enable optical physiology, optogenetics, and sparse labeling of genetically defined cell populations. TIGRE2.0 reporters broke the barrier in transgene expression level of single-copy targeted-insertion transgenesis in a wide range of neuronal types, along with additional advantage of a simplified breeding strategy compared to our first-generation TIGRE lines. These novel transgenic lines greatly expand the repertoire of high-precision genetic tools available to effectively identify, monitor, and manipulate distinct cell types in the mouse brain.

Molecularly defined and spatially resolved cell atlas of the whole mouse brain.

In mammalian brains, millions to billions of cells form complex interaction networks to enable a wide range of functions. The enormous diversity and intricate organization of cells have impeded our understanding of the molecular and cellular basis of brain function. Recent advances in spatially resolved single-cell transcriptomics have enabled systematic mapping of the spatial organization of molecularly defined cell types in complex tissues1-3, including several brain regions (for example, refs. 1-11). However, a comprehensive cell atlas of the whole brain is still missing. Here we imaged a panel of more than 1,100 genes in approximately 10 million cells across the entire adult mouse brains using multiplexed error-robust fluorescence in situ hybridization12 and performed spatially resolved, single-cell expression profiling at the whole-transcriptome scale by integrating multiplexed error-robust fluorescence in situ hybridization and single-cell RNA sequencing data. Using this approach, we generated a comprehensive cell atlas of more than 5,000 transcriptionally distinct cell clusters, belonging to more than 300 major cell types, in the whole mouse brain with high molecular and spatial resolution. Registration of this atlas to the mouse brain common coordinate framework allowed systematic quantifications of the cell-type composition and organization in individual brain regions. We further identified spatial modules characterized by distinct cell-type compositions and spatial gradients featuring gradual changes of cells. Finally, this high-resolution spatial map of cells, each with a transcriptome-wide expression profile, allowed us to infer cell-type-specific interactions between hundreds of cell-type pairs and predict molecular (ligand-receptor) basis and functional implications of these cell-cell interactions. These results provide rich insights into the molecular and cellular architecture of the brain and a foundation for functional investigations of neural circuits and their dysfunction in health and disease.

Single-cell DNA methylome and 3D multi-omic atlas of the adult mouse brain.

Cytosine DNA methylation is essential in brain development and is implicated in various neurological disorders. Understanding DNA methylation diversity across the entire brain in a spatial context is fundamental for a complete molecular atlas of brain cell types and their gene regulatory landscapes. Here we used single-nucleus methylome sequencing (snmC-seq3) and multi-omic sequencing (snm3C-seq)1 technologies to generate 301,626 methylomes and 176,003 chromatin conformation-methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell taxonomy with 4,673 cell groups and 274 cross-modality-annotated subclasses. We identified 2.6 million differentially methylated regions across the genome that represent potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide spatial transcriptomics data validated the association of spatial epigenetic diversity with transcription and improved the anatomical mapping of our epigenetic datasets. Furthermore, chromatin conformation diversities occurred in important neuronal genes and were highly associated with DNA methylation and transcription changes. Brain-wide cell-type comparisons enabled the construction of regulatory networks that incorporate transcription factors, regulatory elements and their potential downstream gene targets. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a whole-brain SMART-seq2 dataset. Our study establishes a brain-wide, single-cell DNA methylome and 3D multi-omic atlas and provides a valuable resource for comprehending the cellular-spatial and regulatory genome diversity of the mouse brain.

A transcriptomic taxonomy of mouse brain-wide spinal projecting neurons.

The brain controls nearly all bodily functions via spinal projecting neurons (SPNs) that carry command signals from the brain to the spinal cord. However, a comprehensive molecular characterization of brain-wide SPNs is still lacking. Here we transcriptionally profiled a total of 65,002 SPNs, identified 76 region-specific SPN types, and mapped these types into a companion atlas of the whole mouse brain1. This taxonomy reveals a three-component organization of SPNs: (1) molecularly homogeneous excitatory SPNs from the cortex, red nucleus and cerebellum with somatotopic spinal terminations suitable for point-to-point communication; (2) heterogeneous populations in the reticular formation with broad spinal termination patterns, suitable for relaying commands related to the activities of the entire spinal cord; and (3) modulatory neurons expressing slow-acting neurotransmitters and/or neuropeptides in the hypothalamus, midbrain and reticular formation for 'gain setting' of brain-spinal signals. In addition, this atlas revealed a LIM homeobox transcription factor code that parcellates the reticulospinal neurons into five molecularly distinct and spatially segregated populations. Finally, we found transcriptional signatures of a subset of SPNs with large soma size and correlated these with fast-firing electrophysiological properties. Together, this study establishes a comprehensive taxonomy of brain-wide SPNs and provides insight into the functional organization of SPNs in mediating brain control of bodily functions.

Single-cell analysis of chromatin accessibility in the adult mouse brain.

Recent advances in single-cell technologies have led to the discovery of thousands of brain cell types; however, our understanding of the gene regulatory programs in these cell types is far from complete1-4. Here we report a comprehensive atlas of candidate cis-regulatory DNA elements (cCREs) in the adult mouse brain, generated by analysing chromatin accessibility in 2.3 million individual brain cells from 117 anatomical dissections. The atlas includes approximately 1 million cCREs and their chromatin accessibility across 1,482 distinct brain cell populations, adding over 446,000 cCREs to the most recent such annotation in the mouse genome. The mouse brain cCREs are moderately conserved in the human brain. The mouse-specific cCREs-specifically, those identified from a subset of cortical excitatory neurons-are strongly enriched for transposable elements, suggesting a potential role for transposable elements in the emergence of new regulatory programs and neuronal diversity. Finally, we infer the gene regulatory networks in over 260 subclasses of mouse brain cells and develop deep-learning models to predict the activities of gene regulatory elements in different brain cell types from the DNA sequence alone. Our results provide a resource for the analysis of cell-type-specific gene regulation programs in both mouse and human brains.

Brain-wide correspondence of neuronal epigenomics and distant projections.

Single-cell analyses parse the brain's billions of neurons into thousands of 'cell-type' clusters residing in different brain structures1. Many cell types mediate their functions through targeted long-distance projections allowing interactions between specific cell types. Here we used epi-retro-seq2 to link single-cell epigenomes and cell types to long-distance projections for 33,034 neurons dissected from 32 different regions projecting to 24 different targets (225 source-to-target combinations) across the whole mouse brain. We highlight uses of these data for interrogating principles relating projection types to transcriptomics and epigenomics, and for addressing hypotheses about cell types and connections related to genetics. We provide an overall synthesis with 926 statistical comparisons of discriminability of neurons projecting to each target for every source. We integrate this dataset into the larger BRAIN Initiative Cell Census Network atlas, composed of millions of neurons, to link projection cell types to consensus clusters. Integration with spatial transcriptomics further assigns projection-enriched clusters to smaller source regions than the original dissections. We exemplify this by presenting in-depth analyses of projection neurons from the hypothalamus, thalamus, hindbrain, amygdala and midbrain to provide insights into properties of those cell types, including differentially expressed genes, their associated cis-regulatory elements and transcription-factor-binding motifs, and neurotransmitter use.

Isoform cell-type specificity in the mouse primary motor cortex.

Full-length SMART-seq1 single-cell RNA sequencing can be used to measure gene expression at isoform resolution, making possible the identification of specific isoform markers for different cell types. Used in conjunction with spatial RNA capture and gene-tagging methods, this enables the inference of spatially resolved isoform expression for different cell types. Here, in a comprehensive analysis of 6,160 mouse primary motor cortex cells assayed with SMART-seq, 280,327 cells assayed with MERFISH2 and 94,162 cells assayed with 10x Genomics sequencing3, we find examples of isoform specificity in cell types-including isoform shifts between cell types that are masked in gene-level analysis-as well as examples of transcriptional regulation. Additionally, we show that isoform specificity helps to refine cell types, and that a multi-platform analysis of single-cell transcriptomic data leveraging multiple measurements provides a comprehensive atlas of transcription in the mouse primary motor cortex that improves on the possibilities offered by any single technology.

Morphological diversity of single neurons in molecularly defined cell types.

Dendritic and axonal morphology reflects the input and output of neurons and is a defining feature of neuronal types1,2, yet our knowledge of its diversity remains limited. Here, to systematically examine complete single-neuron morphologies on a brain-wide scale, we established a pipeline encompassing sparse labelling, whole-brain imaging, reconstruction, registration and analysis. We fully reconstructed 1,741 neurons from cortex, claustrum, thalamus, striatum and other brain regions in mice. We identified 11 major projection neuron types with distinct morphological features and corresponding transcriptomic identities. Extensive projectional diversity was found within each of these major types, on the basis of which some types were clustered into more refined subtypes. This diversity follows a set of generalizable principles that govern long-range axonal projections at different levels, including molecular correspondence, divergent or convergent projection, axon termination pattern, regional specificity, topography, and individual cell variability. Although clear concordance with transcriptomic profiles is evident at the level of major projection type, fine-grained morphological diversity often does not readily correlate with transcriptomic subtypes derived from unsupervised clustering, highlighting the need for single-cell cross-modality studies. Overall, our study demonstrates the crucial need for quantitative description of complete single-cell anatomy in cell-type classification, as single-cell morphological diversity reveals a plethora of ways in which different cell types and their individual members may contribute to the configuration and function of their respective circuits.

Unsupervised pattern identification in spatial gene expression atlas reveals mouse brain regions beyond established ontology.

The rapid growth of large-scale spatial gene expression data demands efficient and reliable computational tools to extract major trends of gene expression in their native spatial context. Here, we used stability-driven unsupervised learning (i.e., staNMF) to identify principal patterns (PPs) of 3D gene expression profiles and understand spatial gene distribution and anatomical localization at the whole mouse brain level. Our subsequent spatial correlation analysis systematically compared the PPs to known anatomical regions and ontology from the Allen Mouse Brain Atlas using spatial neighborhoods. We demonstrate that our stable and spatially coherent PPs, whose linear combinations accurately approximate the spatial gene data, are highly correlated with combinations of expert-annotated brain regions. These PPs yield a brain ontology based purely on spatial gene expression. Our PP identification approach outperforms principal component analysis and typical clustering algorithms on the same task. Moreover, we show that the stable PPs reveal marked regional imbalance of brainwide genetic architecture, leading to region-specific marker genes and gene coexpression networks. Our findings highlight the advantages of stability-driven machine learning for plausible biological discovery from dense spatial gene expression data, streamlining tasks that are infeasible by conventional manual approaches.

The Q system: a repressible binary system for transgene expression, lineage tracing, and mosaic analysis.

We describe a new repressible binary expression system based on the regulatory genes from the Neurospora qa gene cluster. This "Q system" offers attractive features for transgene expression in Drosophila and mammalian cells: low basal expression in the absence of the transcriptional activator QF, high QF-induced expression, and QF repression by its repressor QS. Additionally, feeding flies quinic acid can relieve QS repression. The Q system offers many applications, including (1) intersectional "logic gates" with the GAL4 system for manipulating transgene expression patterns, (2) GAL4-independent MARCM analysis, and (3) coupled MARCM analysis to independently visualize and genetically manipulate siblings from any cell division. We demonstrate the utility of the Q system in determining cell division patterns of a neuronal lineage and gene function in cell growth and proliferation, and in dissecting neurons responsible for olfactory attraction. The Q system can be expanded to other uses in Drosophila and to any organism conducive to transgenesis.

Single-cell transcriptomic classification of rabies-infected cortical neurons.

Cortical circuit tracing using modified rabies virus can identify input neurons making direct monosynaptic connections onto neurons of interest. However, challenges remain in our ability to establish the cell type identity of rabies-labeled input neurons. While transcriptomics may offer an avenue to characterize inputs, the extent of rabies-induced transcriptional changes in distinct neuronal cell types remains unclear, and whether these changes preclude characterization of rabies-infected neurons according to established transcriptomic cell types is unknown. We used single-nucleus RNA sequencing to survey the gene expression profiles of rabies-infected neurons and assessed their correspondence with established transcriptomic cell types. We demonstrated that when using transcriptome-wide RNA profiles, rabies-infected cortical neurons can be transcriptomically characterized despite global and cell-type-specific rabies-induced transcriptional changes. Notably, we found differential modulation of neuronal marker gene expression, suggesting that caution should be taken when attempting to characterize rabies-infected cells with single genes or small gene sets.

Distinct descending motor cortex pathways and their roles in movement.

Activity in the motor cortex predicts movements, seconds before they are initiated. This preparatory activity has been observed across cortical layers, including in descending pyramidal tract neurons in layer 5. A key question is how preparatory activity is maintained without causing movement, and is ultimately converted to a motor command to trigger appropriate movements. Here, using single-cell transcriptional profiling and axonal reconstructions, we identify two types of pyramidal tract neuron. Both types project to several targets in the basal ganglia and brainstem. One type projects to thalamic regions that connect back to motor cortex; populations of these neurons produced early preparatory activity that persisted until the movement was initiated. The second type projects to motor centres in the medulla and mainly produced late preparatory activity and motor commands. These results indicate that two types of motor cortex output neurons have specialized roles in motor control.

Shared and distinct transcriptomic cell types across neocortical areas.

The neocortex contains a multitude of cell types that are segregated into layers and functionally distinct areas. To investigate the diversity of cell types across the mouse neocortex, here we analysed 23,822 cells from two areas at distant poles of the mouse neocortex: the primary visual cortex and the anterior lateral motor cortex. We define 133 transcriptomic cell types by deep, single-cell RNA sequencing. Nearly all types of GABA (γ-aminobutyric acid)-containing neurons are shared across both areas, whereas most types of glutamatergic neurons were found in one of the two areas. By combining single-cell RNA sequencing and retrograde labelling, we match transcriptomic types of glutamatergic neurons to their long-range projection specificity. Our study establishes a combined transcriptomic and projectional taxonomy of cortical cell types from functionally distinct areas of the adult mouse cortex.

RecV recombinase system for in vivo targeted optogenomic modifications of single cells or cell populations.

Brain circuits comprise vast numbers of interconnected neurons with diverse molecular, anatomical and physiological properties. To allow targeting of individual neurons for structural and functional studies, we created light-inducible site-specific DNA recombinases based on Cre, Dre and Flp (RecVs). RecVs can induce genomic modifications by one-photon or two-photon light induction in vivo. They can produce targeted, sparse and strong labeling of individual neurons by modifying multiple loci within mouse and zebrafish genomes. In combination with other genetic strategies, they allow intersectional targeting of different neuronal classes. In the mouse cortex they enable sparse labeling and whole-brain morphological reconstructions of individual neurons. Furthermore, these enzymes allow single-cell two-photon targeted genetic modifications and can be used in combination with functional optical indicators with minimal interference. In summary, RecVs enable spatiotemporally precise optogenomic modifications that can facilitate detailed single-cell analysis of neural circuits by linking genetic identity, morphology, connectivity and function.

Ventral tegmental area glutamate neurons mediate nonassociative consequences of stress.

Exposure to trauma is a risk factor for the development of a number of mood disorders, and may enhance vulnerability to future adverse life events. Recent data demonstrate that ventral tegmental area (VTA) neurons expressing the vesicular glutamate transporter 2 (VGluT2) signal and causally contribute to behaviors that involve aversive or threatening stimuli. However, it is unknown whether VTA VGluT2 neurons regulate transsituational outcomes of stress and whether these neurons are sensitive to stressor controllability. This work adapted an operant mouse paradigm to examine the impact of stressor controllability on VTA VGluT2 neuron function as well as the role of VTA VGluT2 neurons in mediating transsituational stressor outcomes. Uncontrollable (inescapable) stress, but not physically identical controllable (escapable) stress, produced social avoidance and exaggerated fear in male mice. Uncontrollable stress in females led to exploratory avoidance of a novel brightly lit environment. Both controllable and uncontrollable stressors increased VTA VGluT2 neuronal activity, and chemogenetic silencing of VTA VGluT2 neurons prevented the behavioral sequelae of uncontrollable stress in male and female mice. Further, we show that stress activates multiple genetically-distinct subtypes of VTA VGluT2 neurons, especially those that are VGluT2+VGaT+, as well as lateral habenula neurons receiving synaptic input from VTA VGluT2 neurons. Our results provide causal evidence that mice can be used for identifying stressor controllability circuitry and that VTA VGluT2 neurons contribute to transsituational stressor outcomes, such as social avoidance, exaggerated fear, or anxiety-like behavior that are observed within trauma-related disorders.

Identification of preoptic sleep neurons using retrograde labelling and gene profiling.

In humans and other mammalian species, lesions in the preoptic area of the hypothalamus cause profound sleep impairment, indicating a crucial role of the preoptic area in sleep generation. However, the underlying circuit mechanism remains poorly understood. Electrophysiological recordings and c-Fos immunohistochemistry have shown the existence of sleep-active neurons in the preoptic area, especially in the ventrolateral preoptic area and median preoptic nucleus. Pharmacogenetic activation of c-Fos-labelled sleep-active neurons has been shown to induce sleep. However, the sleep-active neurons are spatially intermingled with wake-active neurons, making it difficult to target the sleep neurons specifically for circuit analysis. Here we identify a population of preoptic area sleep neurons on the basis of their projection target and discover their molecular markers. Using a lentivirus expressing channelrhodopsin-2 or a light-activated chloride channel for retrograde labelling, bidirectional optogenetic manipulation, and optrode recording, we show that the preoptic area GABAergic neurons projecting to the tuberomammillary nucleus are both sleep active and sleep promoting. Furthermore, translating ribosome affinity purification and single-cell RNA sequencing identify candidate markers for these neurons, and optogenetic and pharmacogenetic manipulations demonstrate that several peptide markers (cholecystokinin, corticotropin-releasing hormone, and tachykinin 1) label sleep-promoting neurons. Together, these findings provide easy genetic access to sleep-promoting preoptic area neurons and a valuable entry point for dissecting the sleep control circuit.

Cortical representations of olfactory input by trans-synaptic tracing.

In the mouse, each class of olfactory receptor neurons expressing a given odorant receptor has convergent axonal projections to two specific glomeruli in the olfactory bulb, thereby creating an odour map. However, it is unclear how this map is represented in the olfactory cortex. Here we combine rabies-virus-dependent retrograde mono-trans-synaptic labelling with genetics to control the location, number and type of 'starter' cortical neurons, from which we trace their presynaptic neurons. We find that individual cortical neurons receive input from multiple mitral cells representing broadly distributed glomeruli. Different cortical areas represent the olfactory bulb input differently. For example, the cortical amygdala preferentially receives dorsal olfactory bulb input, whereas the piriform cortex samples the whole olfactory bulb without obvious bias. These differences probably reflect different functions of these cortical areas in mediating innate odour preference or associative memory. The trans-synaptic labelling method described here should be widely applicable to mapping connections throughout the mouse nervous system.

DICE, an efficient system for iterative genomic editing in human pluripotent stem cells.

To reveal the full potential of human pluripotent stem cells, new methods for rapid, site-specific genomic engineering are needed. Here, we describe a system for precise genetic modification of human embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs). We identified a novel human locus, H11, located in a safe, intergenic, transcriptionally active region of chromosome 22, as the recipient site, to provide robust, ubiquitous expression of inserted genes. Recipient cell lines were established by site-specific placement of a 'landing pad' cassette carrying attP sites for phiC31 and Bxb1 integrases at the H11 locus by spontaneous or TALEN-assisted homologous recombination. Dual integrase cassette exchange (DICE) mediated by phiC31 and Bxb1 integrases was used to insert genes of interest flanked by phiC31 and Bxb1 attB sites at the H11 locus, replacing the landing pad. This system provided complete control over content, direction and copy number of inserted genes, with a specificity of 100%. A series of genes, including mCherry and various combinations of the neural transcription factors LMX1a, FOXA2 and OTX2, were inserted in recipient cell lines derived from H9 ESC, as well as iPSC lines derived from a Parkinson's disease patient and a normal sibling control. The DICE system offers rapid, efficient and precise gene insertion in ESC and iPSC and is particularly well suited for repeated modifications of the same locus.