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.

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.

Here we report the generation of a multimodal cell census and atlas of the mammalian primary motor cortex as the initial product of the BRAIN Initiative Cell Census Network (BICCN). This was achieved by coordinated large-scale analyses of single-cell transcriptomes, chromatin accessibility, DNA methylomes, spatially resolved single-cell transcriptomes, morphological and electrophysiological properties and cellular resolution input–output mapping, integrated through cross-modal computational analysis. Our results advance the collective knowledge and understanding of brain cell-type organization1–5. First, our study reveals a unified molecular genetic landscape of cortical cell types that integrates their transcriptome, open chromatin and DNA methylation maps. Second, cross-species analysis achieves a consensus taxonomy of transcriptomic types and their hierarchical organization that is conserved from mouse to marmoset and human. Third, in situ single-cell transcriptomics provides a spatially resolved cell-type atlas of the motor cortex. Fourth, cross-modal analysis provides compelling evidence for the transcriptomic, epigenomic and gene regulatory basis of neuronal phenotypes such as their physiological and anatomical properties, demonstrating the biological validity and genomic underpinning of neuron types. We further present an extensive genetic toolset for targeting glutamatergic neuron types towards linking their molecular and developmental identity to their circuit function. Together, our results establish a unifying and mechanistic framework of neuronal cell-type organization that integrates multi-layered molecular genetic and spatial information with multi-faceted phenotypic properties.

Mammalian brain cells show remarkable diversity in gene expression, anatomy and function, yet the regulatory DNA landscape underlying this extensive heterogeneity is poorly understood. Here we carry out a comprehensive assessment of the epigenomes of mouse brain cell types by applying single-nucleus DNA methylation sequencing1,2 to profile 103,982 nuclei (including 95,815 neurons and 8,167 non-neuronal cells) from 45 regions of the mouse cortex, hippocampus, striatum, pallidum and olfactory areas. We identified 161 cell clusters with distinct spatial locations and projection targets. We constructed taxonomies of these epigenetic types, annotated with signature genes, regulatory elements and transcription factors. These features indicate the potential regulatory landscape supporting the assignment of putative cell types and reveal repetitive usage of regulators in excitatory and inhibitory cells for determining subtypes. The DNA methylation landscape of excitatory neurons in the cortex and hippocampus varied continuously along spatial gradients. Using this deep dataset, we constructed an artificial neural network model that precisely predicts single neuron cell-type identity and brain area spatial location. Integration of high-resolution DNA methylomes with single-nucleus chromatin accessibility data3 enabled prediction of high-confidence enhancer–gene interactions for all identified cell types, which were subsequently validated by cell-type-specific chromatin conformation capture experiments4. By combining multi-omic datasets (DNA methylation, chromatin contacts, and open chromatin) from single nuclei and annotating the regulatory genome of hundreds of cell types in the mouse brain, our DNA methylation atlas establishes the epigenetic basis for neuronal diversity and spatial organization throughout the mouse cerebrum.

23The primary motor cortex (M1) is essential for voluntary fine motor control and is functionally conserved 24 across mammals. Using high-throughput transcriptomic and epigenomic profiling of over 450,000 single 25 nuclei in human, marmoset monkey, and mouse, we demonstrate a broadly conserved cellular makeup 26 of this region, whose similarity mirrors evolutionary distance and is consistent between the 27 transcriptome and epigenome. The core conserved molecular identity of neuronal and non-neuronal 28 types allowed the generation of a cross-species consensus cell type classification and inference of 29 conserved cell type properties across species. Despite overall conservation, many species 30 specializations were apparent, including differences in cell type proportions, gene expression, DNA 31 methylation, and chromatin state. Few cell type marker genes were conserved across species, 32 providing a short list of candidate genes and regulatory mechanisms responsible for conserved features 33 of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic 34 classification allowed the Patch-seq identification of layer 5 (L5) corticospinal Betz cells in non-human 35 primate and human and characterization of their highly specialized physiology and anatomy. These 36 findings highlight the robust molecular underpinnings of cell type diversity in M1 across mammals and 37 point to the genes and regulatory pathways responsible for the functional identity of cell types and their 38 species-specific adaptations. 39 40 distinguished on the basis of regions of open chromatin or DNA methylation 5,9,10 . Furthermore, several 48 recent studies have shown that transcriptomically-defined cell types can be aligned across species 2,11-49 13 , indicating that these methods provide a path to quantitatively study evolutionary conservation and 50 divergence at the level of cell types. However, application of these methods has been highly 51 fragmented to date. Human and mouse comparisons have been performed in different cortical regions, 52 using single-cell (with biases in cell proportions) versus single-nucleus (with biases in transcript 53 makeup) analysis, and most single-cell transcriptomic and epigenomic studies have been performed 54 independently. 55 56The primary motor cortex (MOp in mouse, M1 in human and non-human primates, all referred to as M1 57 herein) provides an ideal cortical region to address questions about cellular evolution in rodents and 58 primates by integrating these approaches. Unlike the primary visual cortex (V1), which is highly 59 specialized in primates, or frontal and temporal association areas, whose homologues in rodents 60 remain poorly defined, M1 is essential for fine motor control and is functionally conserved across 61 placental mammals. M1 is an agranular cortex, lacking a defined L4, although neurons with L4-like 62properties have been described 14 . L5 of carnivore and primate M1 contains exceptionally large 63 "giganto-cellular" corticospinal neurons (Betz c...