A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain.

The mammalian brain consists of millions to billions of cells that are organized into many cell types with specific spatial distribution patterns and structural and functional properties1-3. Here we report a comprehensive and high-resolution transcriptomic and spatial cell-type atlas for the whole adult mouse brain. The cell-type atlas was created by combining a single-cell RNA-sequencing (scRNA-seq) dataset of around 7 million cells profiled (approximately 4.0 million cells passing quality control), and a spatial transcriptomic dataset of approximately 4.3 million cells using multiplexed error-robust fluorescence in situ hybridization (MERFISH). The atlas is hierarchically organized into 4 nested levels of classification: 34 classes, 338 subclasses, 1,201 supertypes and 5,322 clusters. We present an online platform, Allen Brain Cell Atlas, to visualize the mouse whole-brain cell-type atlas along with the single-cell RNA-sequencing and MERFISH datasets. We systematically analysed the neuronal and non-neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell-type organization in different brain regions-in particular, a dichotomy between the dorsal and ventral parts of the brain. The dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. Our study also uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types. Finally, we found that transcription factors are major determinants of cell-type classification and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole mouse brain transcriptomic and spatial cell-type atlas establishes a benchmark reference atlas and a foundational resource for integrative investigations of cellular and circuit function, development and evolution of the mammalian brain.

Transcriptomic cytoarchitecture reveals principles of human neocortex organization.

Variation in cytoarchitecture is the basis for the histological definition of cortical areas. We used single cell transcriptomics and performed cellular characterization of the human cortex to better understand cortical areal specialization. Single-nucleus RNA-sequencing of 8 areas spanning cortical structural variation showed a highly consistent cellular makeup for 24 cell subclasses. However, proportions of excitatory neuron subclasses varied substantially, likely reflecting differences in connectivity across primary sensorimotor and association cortices. Laminar organization of astrocytes and oligodendrocytes also differed across areas. Primary visual cortex showed characteristic organization with major changes in the excitatory to inhibitory neuron ratio, expansion of layer 4 excitatory neurons, and specialized inhibitory neurons. These results lay the groundwork for a refined cellular and molecular characterization of human cortical cytoarchitecture and areal specialization.

A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain.

The mammalian brain is composed of millions to billions of cells that are organized into numerous cell types with specific spatial distribution patterns and structural and functional properties. An essential step towards understanding brain function is to obtain a parts list, i.e., a catalog of cell types, of the brain. Here, we report a comprehensive and high-resolution transcriptomic and spatial cell type atlas for the whole adult mouse brain. The cell type atlas was created based on the combination of two single-cell-level, whole-brain-scale datasets: a single-cell RNA-sequencing (scRNA-seq) dataset of ~7 million cells profiled, and a spatially resolved transcriptomic dataset of ~4.3 million cells using MERFISH. The atlas is hierarchically organized into five nested levels of classification: 7 divisions, 32 classes, 306 subclasses, 1,045 supertypes and 5,200 clusters. We systematically analyzed the neuronal, non-neuronal, and immature neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell type organization in different brain regions, in particular, a dichotomy between the dorsal and ventral parts of the brain: the dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. We also systematically characterized cell-type specific expression of neurotransmitters, neuropeptides, and transcription factors. The study uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types across the brain, suggesting they mediate a myriad of modes of intercellular communications. Finally, we found that transcription factors are major determinants of cell type classification in the adult mouse brain and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole-mouse-brain transcriptomic and spatial cell type atlas establishes a benchmark reference atlas and a foundational resource for deep and integrative investigations of cell type and circuit function, development, and evolution of the mammalian brain.

Comparative cellular analysis of motor cortex in human, marmoset and mouse.

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch-seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

A transcriptomic and epigenomic cell atlas of the mouse primary motor cortex.

Single-cell transcriptomics can provide quantitative molecular signatures for large, unbiased samples of the diverse cell types in the brain1-3. With the proliferation of multi-omics datasets, a major challenge is to validate and integrate results into a biological understanding of cell-type organization. Here we generated transcriptomes and epigenomes from more than 500,000 individual cells in the mouse primary motor cortex, a structure that has an evolutionarily conserved role in locomotion. We developed computational and statistical methods to integrate multimodal data and quantitatively validate cell-type reproducibility. The resulting reference atlas-containing over 56 neuronal cell types that are highly replicable across analysis methods, sequencing technologies and modalities-is a comprehensive molecular and genomic account of the diverse neuronal and non-neuronal cell types in the mouse primary motor cortex. The atlas includes a population of excitatory neurons that resemble pyramidal cells in layer 4 in other cortical regions4. We further discovered thousands of concordant marker genes and gene regulatory elements for these cell types. Our results highlight the complex molecular regulation of cell types in the brain and will directly enable the design of reagents to target specific cell types in the mouse primary motor cortex for functional analysis.

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.

Conserved cell types with divergent features in human versus mouse cortex.

Elucidating the cellular architecture of the human cerebral cortex is central to understanding our cognitive abilities and susceptibility to disease. Here we used single-nucleus RNA-sequencing analysis to perform a comprehensive study of cell types in the middle temporal gyrus of human cortex. We identified a highly diverse set of excitatory and inhibitory neuron types that are mostly sparse, with excitatory types being less layer-restricted than expected. Comparison to similar mouse cortex single-cell RNA-sequencing datasets revealed a surprisingly well-conserved cellular architecture that enables matching of homologous types and predictions of properties of human cell types. Despite this general conservation, we also found extensive differences between homologous human and mouse cell types, including marked alterations in proportions, laminar distributions, gene expression and morphology. These species-specific features emphasize the importance of directly studying human brain.

Single-nucleus and single-cell transcriptomes compared in matched cortical cell types.

Transcriptomic profiling of complex tissues by single-nucleus RNA-sequencing (snRNA-seq) affords some advantages over single-cell RNA-sequencing (scRNA-seq). snRNA-seq provides less biased cellular coverage, does not appear to suffer cell isolation-based transcriptional artifacts, and can be applied to archived frozen specimens. We used well-matched snRNA-seq and scRNA-seq datasets from mouse visual cortex to compare cell type detection. Although more transcripts are detected in individual whole cells (~11,000 genes) than nuclei (~7,000 genes), we demonstrate that closely related neuronal cell types can be similarly discriminated with both methods if intronic sequences are included in snRNA-seq analysis. We estimate that the nuclear proportion of total cellular mRNA varies from 20% to over 50% for large and small pyramidal neurons, respectively. Together, these results illustrate the high information content of nuclear RNA for characterization of cellular diversity in brain tissues.

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.

B cell clonal lineage alterations upon recombinant HIV-1 envelope immunization of rhesus macaques.

Broadly neutralizing HIV-1 antibodies (bNAbs) isolated from infected subjects display protective potential in animal models. Their elicitation by immunization is thus highly desirable. The HIV-1 envelope glycoprotein (Env) is the sole viral target of bnAbs, but is also targeted by binding, non-neutralizing antibodies. Env-based immunogens tested so far in various animal species and humans have elicited binding and autologous neutralizing antibodies but not bNAbs (with a few notable exceptions). The underlying reasons for this are not well understood despite intensive efforts to characterize the binding specificities of the elicited antibodies; mostly by employing serologic methodologies and monoclonal antibody isolation and characterization. These approaches provide limited information on the ontogenies and clonal B cell lineages that expand following Env-immunization. Thus, our current understanding on how the expansion of particular B cell lineages by Env may be linked to the development of non-neutralizing antibodies is limited. Here, in addition to serological analysis, we employed high-throughput BCR sequence analysis from the periphery, lymph nodes and bone marrow, as well as B cell- and antibody-isolation and characterization methods, to compare in great detail the B cell and antibody responses elicited in non-human primates by two forms of the clade C HIV Env 426c: one representing the full length extracellular portion of Env while the other lacking the variable domains 1, 2 and 3 and three conserved N-linked glycosylation sites. The two forms were equally immunogenic, but only the latter elicited neutralizing antibodies by stimulating a more restricted expansion of B cells to a narrower set of IGH/IGK/IGL-V genes that represented a small fraction (0.003-0.02%) of total B cells. Our study provides new information on how Env antigenic differences drastically affect the expansion of particular B cell lineages and supports immunogen-design efforts aiming at stimulating the expansion of cells expressing particular B cell receptors.

Murine spinal cord transcriptome analysis following reduction of prevalent myelin cDNA sequences.

From 1,000 randomly selected colonies from cDNA libraries derived from murine spinal cord subtracted against white matter by means of suppression subtractive hybridization, 220 clones were identified as differentially expressed by dot blot analysis. Sequence analysis by the BLAST programming identified 140 unique genes. (1) The percentage of known sequences from myelin and other glial sources was reduced by approximately 75% over previous, similar subtractions employing visual cortex as the driver. (2) Differentially expressed genes tended to reflect existing expectations concerning structure and function of the spinal cord. (3) About 35% of all genes differentially expressed in the spinal cord in this study are also known to be differentially expressed for this structure as tabulated in the UniGene database. (4) About 33% of all genes differentially expressed in the present study are recorded as not present when measured in the spinal cord according to the UniGene database indicating that present techniques are not recording about a third of differentially expressed genes in this structure. (5) About 15% of all differentially expressed genes are for unknown, putative or hypothetical protein products. (6) About 4% of all differentially expressed genes are novel expressed sequence tags for the mouse. The current study demonstrates the importance of reducing the presence of glial associated sequences when comparing brain regions. It is concluded that the persistence of some myelin sequences in the spinal cord when white matter is used as the driver indicates that myelination is more active in this structure than for those areas represented by white matter and corpus callosum.

Spinal cord transcriptome analysis using suppression subtractive hybridization and mirror orientation selection.

Comparison of cDNA libraries derived from the spinal cord with those derived from the visual cortex by means of forward and reverse subtractive hybridization resulted in the cataloguing of 60 genes differentially expressed in the spinal cord. 1. The differentially expressed genes represent a mixture of novel and known sequences with known and unknown protein products. 2. The possibility that the subtraction process was simply overwhelmed by background sequences was significantly reduced by several observations including comparisons between suppression subtractive hybridization (SSH) and mirror orientation selection (MOS). 3. Nearly half of all genes up-regulated in the spinal cord are of myelin origin. 4. Twenty-five percent of all up-regulated clones in the spinal cord versus the visual cortex are for proteolipid protein. 5. Ten percent of all up-regulated clones in spinal cord versus visual cortex are for ferretin heavy chain, which is known to be produced in oligodendroglial cells in the CNS. 6. Two of the up-regulated sequences, proteolipid protein and N-myc down-regulated gene 4, are identified with genes known to directly affect neuron survival. 7. Two of the up-regulated genes, ferritin and transferrin, are indirectly associated with apoptosis through their ability to sequester iron and reduce free radical formation.