Integrated Morphoelectric and Transcriptomic Classification of Cortical GABAergic Cells.

Neurons are frequently classified into distinct types on the basis of structural, physiological, or genetic attributes. To better constrain the definition of neuronal cell types, we characterized the transcriptomes and intrinsic physiological properties of over 4,200 mouse visual cortical GABAergic interneurons and reconstructed the local morphologies of 517 of those neurons. We find that most transcriptomic types (t-types) occupy specific laminar positions within visual cortex, and, for most types, the cells mapping to a t-type exhibit consistent electrophysiological and morphological properties. These properties display both discrete and continuous variation among t-types. Through multimodal integrated analysis, we define 28 met-types that have congruent morphological, electrophysiological, and transcriptomic properties and robust mutual predictability. We identify layer-specific axon innervation pattern as a defining feature distinguishing different met-types. These met-types represent a unified definition of cortical GABAergic interneuron types, providing a systematic framework to capture existing knowledge and bridge future analyses across different modalities.

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.

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.

Human neocortical expansion involves glutamatergic neuron diversification.

The neocortex is disproportionately expanded in human compared with mouse1,2, both in its total volume relative to subcortical structures and in the proportion occupied by supragranular layers composed of neurons that selectively make connections within the neocortex and with other telencephalic structures. Single-cell transcriptomic analyses of human and mouse neocortex show an increased diversity of glutamatergic neuron types in supragranular layers in human neocortex and pronounced gradients as a function of cortical depth3. Here, to probe the functional and anatomical correlates of this transcriptomic diversity, we developed a robust platform combining patch clamp recording, biocytin staining and single-cell RNA-sequencing (Patch-seq) to examine neurosurgically resected human tissues. We demonstrate a strong correspondence between morphological, physiological and transcriptomic phenotypes of five human glutamatergic supragranular neuron types. These were enriched in but not restricted to layers, with one type varying continuously in all phenotypes across layers 2 and 3. The deep portion of layer 3 contained highly distinctive cell types, two of which express a neurofilament protein that labels long-range projection neurons in primates that are selectively depleted in Alzheimer's disease4,5. Together, these results demonstrate the explanatory power of transcriptomic cell-type classification, provide a structural underpinning for increased complexity of cortical function in humans, and implicate discrete transcriptomic neuron types as selectively vulnerable in disease.

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.

A comprehensive transcriptional map of primate brain development.

The transcriptional underpinnings of brain development remain poorly understood, particularly in humans and closely related non-human primates. We describe a high-resolution transcriptional atlas of rhesus monkey (Macaca mulatta) brain development that combines dense temporal sampling of prenatal and postnatal periods with fine anatomical division of cortical and subcortical regions associated with human neuropsychiatric disease. Gene expression changes more rapidly before birth, both in progenitor cells and maturing neurons. Cortical layers and areas acquire adult-like molecular profiles surprisingly late in postnatal development. Disparate cell populations exhibit distinct developmental timing of gene expression, but also unexpected synchrony of processes underlying neural circuit construction including cell projection and adhesion. Candidate risk genes for neurodevelopmental disorders including primary microcephaly, autism spectrum disorder, intellectual disability, and schizophrenia show disease-specific spatiotemporal enrichment within developing neocortex. Human developmental expression trajectories are more similar to monkey than rodent, although approximately 9% of genes show human-specific regulation with evidence for prolonged maturation or neoteny compared to monkey.

Transcriptional landscape of the prenatal human brain.

The anatomical and functional architecture of the human brain is mainly determined by prenatal transcriptional processes. We describe an anatomically comprehensive atlas of the mid-gestational human brain, including de novo reference atlases, in situ hybridization, ultra-high-resolution magnetic resonance imaging (MRI) and microarray analysis on highly discrete laser-microdissected brain regions. In developing cerebral cortex, transcriptional differences are found between different proliferative and post-mitotic layers, wherein laminar signatures reflect cellular composition and developmental processes. Cytoarchitectural differences between human and mouse have molecular correlates, including species differences in gene expression in subplate, although surprisingly we find minimal differences between the inner and outer subventricular zones even though the outer zone is expanded in humans. Both germinal and post-mitotic cortical layers exhibit fronto-temporal gradients, with particular enrichment in the frontal lobe. Finally, many neurodevelopmental disorder and human-evolution-related genes show patterned expression, potentially underlying unique features of human cortical formation. These data provide a rich, freely-accessible resource for understanding human brain development.

An anatomically comprehensive atlas of the adult human brain transcriptome.

Neuroanatomically precise, genome-wide maps of transcript distributions are critical resources to complement genomic sequence data and to correlate functional and genetic brain architecture. Here we describe the generation and analysis of a transcriptional atlas of the adult human brain, comprising extensive histological analysis and comprehensive microarray profiling of ∼900 neuroanatomically precise subdivisions in two individuals. Transcriptional regulation varies enormously by anatomical location, with different regions and their constituent cell types displaying robust molecular signatures that are highly conserved between individuals. Analysis of differential gene expression and gene co-expression relationships demonstrates that brain-wide variation strongly reflects the distributions of major cell classes such as neurons, oligodendrocytes, astrocytes and microglia. Local neighbourhood relationships between fine anatomical subdivisions are associated with discrete neuronal subtypes and genes involved with synaptic transmission. The neocortex displays a relatively homogeneous transcriptional pattern, but with distinct features associated selectively with primary sensorimotor cortices and with enriched frontal lobe expression. Notably, the spatial topography of the neocortex is strongly reflected in its molecular topography-the closer two cortical regions, the more similar their transcriptomes. This freely accessible online data resource forms a high-resolution transcriptional baseline for neurogenetic studies of normal and abnormal human brain function.

Conserved molecular signatures of neurogenesis in the hippocampal subgranular zone of rodents and primates.

The neurogenic potential of the subgranular zone (SGZ) of the hippocampal dentate gyrus is likely to be regulated by molecular cues arising from its complex heterogeneous cellular environment. Through transcriptome analysis using laser microdissection coupled with DNA microarrays, in combination with analysis of genome-wide in situ hybridization data, we identified 363 genes selectively enriched in adult mouse SGZ. These genes reflect expression in the different constituent cell types, including progenitor and dividing cells, immature granule cells, astrocytes, oligodendrocytes and GABAergic interneurons. Similar transcriptional profiling in the rhesus monkey dentate gyrus across postnatal development identified a highly overlapping set of SGZ-enriched genes, which can be divided based on temporal profiles to reflect maturation of glia versus granule neurons. Furthermore, we identified a neurogenesis-related gene network with decreasing postnatal expression that is highly correlated with the declining number of proliferating cells in dentate gyrus over postnatal development. Many of the genes in this network showed similar postnatal downregulation in mouse, suggesting a conservation of molecular mechanisms underlying developmental and adult neurogenesis in rodents and primates. Conditional deletion of Sox4 and Sox11, encoding two neurogenesis-related transcription factors central in this network, produces a mouse with no hippocampus, confirming the crucial role for these genes in regulating hippocampal neurogenesis.

High-throughput analysis of dendrite and axonal arbors reveals transcriptomic correlates of neuroanatomy.

Neuronal anatomy is central to the organization and function of brain cell types. However, anatomical variability within apparently homogeneous populations of cells can obscure such insights. Here, we report large-scale automation of neuronal morphology reconstruction and analysis on a dataset of 813 inhibitory neurons characterized using the Patch-seq method, which enables measurement of multiple properties from individual neurons, including local morphology and transcriptional signature. We demonstrate that these automated reconstructions can be used in the same manner as manual reconstructions to understand the relationship between some, but not all, cellular properties used to define cell types. We uncover gene expression correlates of laminar innervation on multiple transcriptomically defined neuronal subclasses and types. In particular, our results reveal correlates of the variability in Layer 1 (L1) axonal innervation in a transcriptomically defined subpopulation of Martinotti cells in the adult mouse neocortex.

The ANTsX Ecosystem for Mapping the Mouse Brain.

Precision mapping techniques coupled with high resolution image acquisition of the mouse brain permit the study of the spatial organization of gene expression and their mutual interaction for a comprehensive view of salient structural/functional relationships. Such research is facilitated by standardized anatomical coordinate systems, such as the well-known Allen Common Coordinate Framework (AllenCCFv3), and the ability to spatially map to such standardized spaces. The Advanced Normalization Tools Ecosystem is a comprehensive open-source software toolkit for generalized quantitative imaging with applicability to multiple organ systems, modalities, and animal species. Herein, we illustrate the utility of ANTsX for generating precision spatial mappings of the mouse brain and potential subsequent quantitation. We describe ANTsX-based workflows for mapping domain-specific image data to AllenCCFv3 accounting for common artefacts and other confounds. Novel contributions include ANTsX functionality for velocity flow-based mapping spanning the spatiotemporal domain of a longitudinal trajectory which we apply to the Developmental Common Coordinate Framework. Additionally, we present an automated structural morphological pipeline for determining volumetric and cortical thickness measurements analogous to the well-utilized ANTsX pipeline for human neuroanatomical structural morphology which illustrates a general open-source framework for tailored brain parcellations.

Areal specializations in the morpho-electric and transcriptomic properties of primate layer 5 extratelencephalic projection neurons.

Large-scale analysis of single-cell gene expression has revealed transcriptomically defined cell subclasses present throughout the primate neocortex with gene expression profiles that differ depending upon neocortical region. Here, we test whether the interareal differences in gene expression translate to regional specializations in the physiology and morphology of infragranular glutamatergic neurons by performing Patch-seq experiments in brain slices from the temporal cortex (TCx) and motor cortex (MCx) of the macaque. We confirm that transcriptomically defined extratelencephalically projecting neurons of layer 5 (L5 ET neurons) include retrogradely labeled corticospinal neurons in the MCx and find multiple physiological properties and ion channel genes that distinguish L5 ET from non-ET neurons in both areas. Additionally, while infragranular ET and non-ET neurons retain distinct neuronal properties across multiple regions, there are regional morpho-electric and gene expression specializations in the L5 ET subclass, providing mechanistic insights into the specialized functional architecture of the primate neocortex.

Integrated multimodal cell atlas of Alzheimer's disease.

Alzheimer's disease (AD) is the leading cause of dementia in older adults. Although AD progression is characterized by stereotyped accumulation of proteinopathies, the affected cellular populations remain understudied. Here we use multiomics, spatial genomics and reference atlases from the BRAIN Initiative to study middle temporal gyrus cell types in 84 donors with varying AD pathologies. This cohort includes 33 male donors and 51 female donors, with an average age at time of death of 88 years. We used quantitative neuropathology to place donors along a disease pseudoprogression score. Pseudoprogression analysis revealed two disease phases: an early phase with a slow increase in pathology, presence of inflammatory microglia, reactive astrocytes, loss of somatostatin+ inhibitory neurons, and a remyelination response by oligodendrocyte precursor cells; and a later phase with exponential increase in pathology, loss of excitatory neurons and Pvalb+ and Vip+ inhibitory neuron subtypes. These findings were replicated in other major AD studies.

Trypanosoma brucei mitochondrial respiratome: composition and organization in procyclic form.

The mitochondrial respiratory chain is comprised of four different protein complexes (I-IV), which are responsible for electron transport and generation of proton gradient in the mitochondrial intermembrane space. This proton gradient is then used by F0F1-ATP synthase (complex V) to produce ATP by oxidative phosphorylation. In this study, the respiratory complexes I, II, and III were affinity purified from Trypanosoma brucei procyclic form cells and their composition was determined by mass spectrometry. The results along with those that we previously reported for complexes IV and V showed that the respiratome of Trypanosoma is divergent because many of its proteins are unique to this group of organisms. The studies also identified two mitochondrial subunit proteins of respiratory complex IV that are encoded by edited RNAs. Proteomics data from analyses of complexes purified using numerous tagged component proteins in each of the five complexes were used to generate the first predicted protein-protein interaction network of the Trypanosoma brucei respiratory chain. These results provide the first comprehensive insight into the unique composition of the respiratory complexes in Trypanosoma brucei, an early diverged eukaryotic pathogen.

Morphoelectric and transcriptomic divergence of the layer 1 interneuron repertoire in human versus mouse neocortex.

Neocortical layer 1 (L1) is a site of convergence between pyramidal-neuron dendrites and feedback axons where local inhibitory signaling can profoundly shape cortical processing. Evolutionary expansion of human neocortex is marked by distinctive pyramidal neurons with extensive L1 branching, but whether L1 interneurons are similarly diverse is underexplored. Using Patch-seq recordings from human neurosurgical tissue, we identified four transcriptomic subclasses with mouse L1 homologs, along with distinct subtypes and types unmatched in mouse L1. Subclass and subtype comparisons showed stronger transcriptomic differences in human L1 and were correlated with strong morphoelectric variability along dimensions distinct from mouse L1 variability. Accompanied by greater layer thickness and other cytoarchitecture changes, these findings suggest that L1 has diverged in evolution, reflecting the demands of regulating the expanded human neocortical circuit.

Signature morphoelectric properties of diverse GABAergic interneurons in the human neocortex.

Human cortex transcriptomic studies have revealed a hierarchical organization of γ-aminobutyric acid-producing (GABAergic) neurons from subclasses to a high diversity of more granular types. Rapid GABAergic neuron viral genetic labeling plus Patch-seq (patch-clamp electrophysiology plus single-cell RNA sequencing) sampling in human brain slices was used to reliably target and analyze GABAergic neuron subclasses and individual transcriptomic types. This characterization elucidated transitions between PVALB and SST subclasses, revealed morphological heterogeneity within an abundant transcriptomic type, identified multiple spatially distinct types of the primate-specialized double bouquet cells (DBCs), and shed light on cellular differences between homologous mouse and human neocortical GABAergic neuron types. These results highlight the importance of multimodal phenotypic characterization for refinement of emerging transcriptomic cell type taxonomies and for understanding conserved and specialized cellular properties of human brain cell types.

Integrated multimodal cell atlas of Alzheimer's disease.

Alzheimer's disease (AD) is the most common cause of dementia in older adults. Neuropathological and imaging studies have demonstrated a progressive and stereotyped accumulation of protein aggregates, but the underlying molecular and cellular mechanisms driving AD progression and vulnerable cell populations affected by disease remain coarsely understood. The current study harnesses single cell and spatial genomics tools and knowledge from the BRAIN Initiative Cell Census Network to understand the impact of disease progression on middle temporal gyrus cell types. We used image-based quantitative neuropathology to place 84 donors spanning the spectrum of AD pathology along a continuous disease pseudoprogression score and multiomic technologies to profile single nuclei from each donor, mapping their transcriptomes, epigenomes, and spatial coordinates to a common cell type reference with unprecedented resolution. Temporal analysis of cell-type proportions indicated an early reduction of Somatostatin-expressing neuronal subtypes and a late decrease of supragranular intratelencephalic-projecting excitatory and Parvalbumin-expressing neurons, with increases in disease-associated microglial and astrocytic states. We found complex gene expression differences, ranging from global to cell type-specific effects. These effects showed different temporal patterns indicating diverse cellular perturbations as a function of disease progression. A subset of donors showed a particularly severe cellular and molecular phenotype, which correlated with steeper cognitive decline. We have created a freely available public resource to explore these data and to accelerate progress in AD research at SEA-AD.org.

Connecting single-cell transcriptomes to projectomes in mouse visual cortex.

The mammalian brain is composed of diverse neuron types that play different functional roles. Recent single-cell RNA sequencing approaches have led to a whole brain taxonomy of transcriptomically-defined cell types, yet cell type definitions that include multiple cellular properties can offer additional insights into a neuron's role in brain circuits. While the Patch-seq method can investigate how transcriptomic properties relate to the local morphological and electrophysiological properties of cell types, linking transcriptomic identities to long-range projections is a major unresolved challenge. To address this, we collected coordinated Patch-seq and whole brain morphology data sets of excitatory neurons in mouse visual cortex. From the Patch-seq data, we defined 16 integrated morpho-electric-transcriptomic (MET)-types; in parallel, we reconstructed the complete morphologies of 300 neurons. We unified the two data sets with a multi-step classifier, to integrate cell type assignments and interrogate cross-modality relationships. We find that transcriptomic variations within and across MET-types correspond with morphological and electrophysiological phenotypes. In addition, this variation, along with the anatomical location of the cell, can be used to predict the projection targets of individual neurons. We also shed new light on infragranular cell types and circuits, including cell-type-specific, interhemispheric projections. With this approach, we establish a comprehensive, integrated taxonomy of excitatory neuron types in mouse visual cortex and create a system for integrated, high-dimensional cell type classification that can be extended to the whole brain and potentially across species.