Molecular Diversity and Specializations among the Cells of the Adult Mouse Brain.

The mammalian brain is composed of diverse, specialized cell populations. To systematically ascertain and learn from these cellular specializations, we used Drop-seq to profile RNA expression in 690,000 individual cells sampled from 9 regions of the adult mouse brain. We identified 565 transcriptionally distinct groups of cells using computational approaches developed to distinguish biological from technical signals. Cross-region analysis of these 565 cell populations revealed features of brain organization, including a gene-expression module for synthesizing axonal and presynaptic components, patterns in the co-deployment of voltage-gated ion channels, functional distinctions among the cells of the vasculature and specialization of glutamatergic neurons across cortical regions. Systematic neuronal classifications for two complex basal ganglia nuclei and the striatum revealed a rare population of spiny projection neurons. This adult mouse brain cell atlas, accessible through interactive online software (DropViz), serves as a reference for development, disease, and evolution.

A molecular and cellular perspective on human brain evolution and tempo.

The evolution of the modern human brain was accompanied by distinct molecular and cellular specializations, which underpin our diverse cognitive abilities but also increase our susceptibility to neurological diseases. These features, some specific to humans and others shared with related species, manifest during different stages of brain development. In this multi-stage process, neural stem cells proliferate to produce a large and diverse progenitor pool, giving rise to excitatory or inhibitory neurons that integrate into circuits during further maturation. This process unfolds over varying time scales across species and has progressively become slower in the human lineage, with differences in tempo correlating with differences in brain size, cell number and diversity, and connectivity. Here we introduce the terms 'bradychrony' and 'tachycrony' to describe slowed and accelerated developmental tempos, respectively. We review how recent technical advances across disciplines, including advanced engineering of in vitro models, functional comparative genetics and high-throughput single-cell profiling, are leading to a deeper understanding of how specializations of the human brain arise during bradychronic neurodevelopment. Emerging insights point to a central role for genetics, gene-regulatory networks, cellular innovations and developmental tempo, which together contribute to the establishment of human specializations during various stages of neurodevelopment and at different points in evolution.

Calcium-permeable AMPA receptors govern PV neuron feature selectivity.

The brain helps us survive by forming internal representations of the external world1,2. Excitatory cortical neurons are often precisely tuned to specific external stimuli3,4. However, inhibitory neurons, such as parvalbumin-positive (PV) interneurons, are generally less selective5. PV interneurons differ from excitatory neurons in their neurotransmitter receptor subtypes, including AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptors (AMPARs)6,7. Excitatory neurons express calcium-impermeable AMPARs that contain the GluA2 subunit (encoded by GRIA2), whereas PV interneurons express receptors that lack the GluA2 subunit and are calcium-permeable (CP-AMPARs). Here we demonstrate a causal relationship between CP-AMPAR expression and the low feature selectivity of PV interneurons. We find low expression stoichiometry of GRIA2 mRNA relative to other subunits in PV interneurons that is conserved across ferrets, rodents, marmosets and humans, and causes abundant CP-AMPAR expression. Replacing CP-AMPARs in PV interneurons with calcium-impermeable AMPARs increased their orientation selectivity in the visual cortex. Manipulations to induce sparse CP-AMPAR expression demonstrated that this increase was cell-autonomous and could occur with changes beyond development. Notably, excitatory-PV interneuron connectivity rates and unitary synaptic strength were unaltered by CP-AMPAR removal, which suggested that the selectivity of PV interneurons can be altered without markedly changing connectivity. In Gria2-knockout mice, in which all AMPARs are calcium-permeable, excitatory neurons showed significantly degraded orientation selectivity, which suggested that CP-AMPARs are sufficient to drive lower selectivity regardless of cell type. Moreover, hippocampal PV interneurons, which usually exhibit low spatial tuning, became more spatially selective after removing CP-AMPARs, which indicated that CP-AMPARs suppress the feature selectivity of PV interneurons independent of modality. These results reveal a new role of CP-AMPARs in maintaining low-selectivity sensory representation in PV interneurons and implicate a conserved molecular mechanism that distinguishes this cell type in the neocortex.

Conserved and divergent gene regulatory programs of the mammalian neocortex.

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

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.

Author Correction: Innovations present in the primate interneuron repertoire.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Innovations present in the primate interneuron repertoire.

Primates and rodents, which descended from a common ancestor around 90 million years ago1, exhibit profound differences in behaviour and cognitive capacity; the cellular basis for these differences is unknown. Here we use single-nucleus RNA sequencing to profile RNA expression in 188,776 individual interneurons across homologous brain regions from three primates (human, macaque and marmoset), a rodent (mouse) and a weasel (ferret). Homologous interneuron types-which were readily identified by their RNA-expression patterns-varied in abundance and RNA expression among ferrets, mice and primates, but varied less among primates. Only a modest fraction of the genes identified as 'markers' of specific interneuron subtypes in any one species had this property in another species. In the primate neocortex, dozens of genes showed spatial expression gradients among interneurons of the same type, which suggests that regional variation in cortical contexts shapes the RNA expression patterns of adult neocortical interneurons. We found that an interneuron type that was previously associated with the mouse hippocampus-the 'ivy cell', which has neurogliaform characteristics-has become abundant across the neocortex of humans, macaques and marmosets but not mice or ferrets. We also found a notable subcortical innovation: an abundant striatal interneuron type in primates that had no molecularly homologous counterpart in mice or ferrets. These interneurons expressed a unique combination of genes that encode transcription factors, receptors and neuropeptides and constituted around 30% of striatal interneurons in marmosets and humans.

Transcriptional profiles of supragranular-enriched genes associate with corticocortical network architecture in the human brain.

The human brain is patterned with disproportionately large, distributed cerebral networks that connect multiple association zones in the frontal, temporal, and parietal lobes. The expansion of the cortical surface, along with the emergence of long-range connectivity networks, may be reflected in changes to the underlying molecular architecture. Using the Allen Institute's human brain transcriptional atlas, we demonstrate that genes particularly enriched in supragranular layers of the human cerebral cortex relative to mouse distinguish major cortical classes. The topography of transcriptional expression reflects large-scale brain network organization consistent with estimates from functional connectivity MRI and anatomical tracing in nonhuman primates. Microarray expression data for genes preferentially expressed in human upper layers (II/III), but enriched only in lower layers (V/VI) of mouse, were cross-correlated to identify molecular profiles across the cerebral cortex of postmortem human brains (n = 6). Unimodal sensory and motor zones have similar molecular profiles, despite being distributed across the cortical mantle. Sensory/motor profiles were anticorrelated with paralimbic and certain distributed association network profiles. Tests of alternative gene sets did not consistently distinguish sensory and motor regions from paralimbic and association regions: (i) genes enriched in supragranular layers in both humans and mice, (ii) genes cortically enriched in humans relative to nonhuman primates, (iii) genes related to connectivity in rodents, (iv) genes associated with human and mouse connectivity, and (v) 1,454 gene sets curated from known gene ontologies. Molecular innovations of upper cortical layers may be an important component in the evolution of long-range corticocortical projections.

Functional Specialization and Flexibility in Human Association Cortex.

The association cortex supports cognitive functions enabling flexible behavior. Here, we explored the organization of human association cortex by mathematically formalizing the notion that a behavioral task engages multiple cognitive components, which are in turn supported by multiple overlapping brain regions. Application of the model to a large data set of neuroimaging experiments (N = 10 449) identified complex zones of frontal and parietal regions that ranged from being highly specialized to highly flexible. The network organization of the specialized and flexible regions was explored with an independent resting-state fMRI data set (N = 1000). Cortical regions specialized for the same components were strongly coupled, suggesting that components function as partially isolated networks. Functionally flexible regions participated in multiple components to different degrees. This heterogeneous selectivity was predicted by the connectivity between flexible and specialized regions. Functionally flexible regions might support binding or integrating specialized brain networks that, in turn, contribute to the ability to execute multiple and varied tasks.

Estimates of segregation and overlap of functional connectivity networks in the human cerebral cortex.

The organization of the human cerebral cortex has recently been explored using techniques for parcellating the cortex into distinct functionally coupled networks. The divergent and convergent nature of cortico-cortical anatomic connections suggests the need to consider the possibility of regions belonging to multiple networks and hierarchies among networks. Here we applied the Latent Dirichlet Allocation (LDA) model and spatial independent component analysis (ICA) to solve for functionally coupled cerebral networks without assuming that cortical regions belong to a single network. Data analyzed included 1000 subjects from the Brain Genomics Superstruct Project (GSP) and 12 high quality individual subjects from the Human Connectome Project (HCP). The organization of the cerebral cortex was similar regardless of whether a winner-take-all approach or the more relaxed constraints of LDA (or ICA) were imposed. This suggests that large-scale networks may function as partially isolated modules. Several notable interactions among networks were uncovered by the LDA analysis. Many association regions belong to at least two networks, while somatomotor and early visual cortices are especially isolated. As examples of interaction, the precuneus, lateral temporal cortex, medial prefrontal cortex and posterior parietal cortex participate in multiple paralimbic networks that together comprise subsystems of the default network. In addition, regions at or near the frontal eye field and human lateral intraparietal area homologue participate in multiple hierarchically organized networks. These observations were replicated in both datasets and could be detected (and replicated) in individual subjects from the HCP.

Conservation, alteration, and redistribution of mammalian striatal interneurons.

Mammalian brains vary in size, structure, and function, but the extent to which evolutionarily novel cell types contribute to this variation remains unresolved1-4. Recent studies suggest there is a primate-specific population of striatal inhibitory interneurons, the TAC3 interneurons5. However, there has not yet been a detailed analysis of the spatial and phylogenetic distribution of this population. Here, we profile single cell gene expression in the developing pig (an ungulate) and ferret (a carnivore), representing 94 million years divergence from primates, and assign newborn inhibitory neurons to initial classes first specified during development6. We find that the initial class of TAC3 interneurons represents an ancestral striatal population that is also deployed towards the cortex in pig and ferret. In adult mouse, we uncover a rare population expressing Tac2, the ortholog of TAC3, in ventromedial striatum, prompting a reexamination of developing mouse striatal interneuron initial classes by targeted enrichment of their precursors. We conclude that the TAC3 interneuron initial class is conserved across Boreoeutherian mammals, with the mouse population representing Th striatal interneurons, a subset of which expresses Tac2. This study suggests that initial classes of telencephalic inhibitory neurons are largely conserved and that during evolution, neuronal types in the mammalian brain change through redistribution and fate refinement, rather than by derivation of novel precursors early in development.

Kynurenic acid inflammatory signaling expands in primates and impairs prefrontal cortical cognition.

Cognitive deficits from dorsolateral prefrontal cortex (dlPFC) dysfunction are common in neuroinflammatory disorders, including long-COVID, schizophrenia and Alzheimer's disease, and have been correlated with kynurenine inflammatory signaling. Kynurenine is further metabolized to kynurenic acid (KYNA) in brain, where it blocks NMDA and α7-nicotinic receptors (nic-α7Rs). These receptors are essential for neurotransmission in dlPFC, suggesting that KYNA may cause higher cognitive deficits in these disorders. The current study found that KYNA and its synthetic enzyme, KAT II, have greatly expanded expression in primate dlPFC in both glia and neurons. Local application of KYNA onto dlPFC neurons markedly reduced the delay-related firing needed for working memory via actions at NMDA and nic-α7Rs, while inhibition of KAT II enhanced neuronal firing in aged macaques. Systemic administration of agents that reduce KYNA production similarly improved cognitive performance in aged monkeys, suggesting a therapeutic avenue for the treatment of cognitive deficits in neuroinflammatory disorders.

Astrocyte regional specialization is shaped by postnatal development.

Astrocytes are an abundant class of glial cells with critical roles in neural circuit assembly and function. Though many studies have uncovered significant molecular distinctions between astrocytes from different brain regions, how this regionalization unfolds over development is not fully understood. We used single-nucleus RNA sequencing to characterize the molecular diversity of brain cells across six developmental stages and four brain regions in the mouse and marmoset brain. Our analysis of over 170,000 single astrocyte nuclei revealed striking regional heterogeneity among astrocytes, particularly between telencephalic and diencephalic regions, at all developmental time points surveyed in both species. At the stages sampled, most of the region patterning was private to astrocytes and not shared with neurons or other glial types. Though astrocytes were already regionally patterned in late embryonic stages, this region-specific astrocyte gene expression signature changed dramatically over postnatal development, and its composition suggests that regional astrocytes further specialize postnatally to support their local neuronal circuits. Comparing across species, we found divergence in the expression of astrocytic region- and age-differentially expressed genes and the timing of astrocyte maturation relative to birth between mouse and marmoset, as well as hundreds of species differentially expressed genes. Finally, we used expansion microscopy to show that astrocyte morphology is largely conserved across gray matter regions of prefrontal cortex, striatum, and thalamus in the mouse, despite substantial molecular divergence.

Key Roles of CACNA1C/Cav1.2 and CALB1/Calbindin in Prefrontal Neurons Altered in Cognitive Disorders.

Importance: The risk of mental disorders is consistently associated with variants in CACNA1C (L-type calcium channel Cav1.2) but it is not known why these channels are critical to cognition, and whether they affect the layer III pyramidal cells in the dorsolateral prefrontal cortex that are especially vulnerable in cognitive disorders. Objective: To examine the molecular mechanisms expressed in layer III pyramidal cells in primate dorsolateral prefrontal cortices. Design, Setting, and Participants: The design included transcriptomic analyses from human and macaque dorsolateral prefrontal cortex, and connectivity, protein expression, physiology, and cognitive behavior in macaques. The research was performed in academic laboratories at Yale, Harvard, Princeton, and the University of Pittsburgh. As dorsolateral prefrontal cortex only exists in primates, the work evaluated humans and macaques. Main Outcomes and Measures: Outcome measures included transcriptomic signatures of human and macaque pyramidal cells, protein expression and interactions in layer III macaque pyramidal cells using light and electron microscopy, changes in neuronal firing during spatial working memory, and working memory performance following pharmacological treatments. Results: Layer III pyramidal cells in dorsolateral prefrontal cortex coexpress a constellation of calcium-related proteins, delineated by CALB1 (calbindin), and high levels of CACNA1C (Cav1.2), GRIN2B (NMDA receptor GluN2B), and KCNN3 (SK3 potassium channel), concentrated in dendritic spines near the calcium-storing smooth endoplasmic reticulum. L-type calcium channels influenced neuronal firing needed for working memory, where either blockade or increased drive by ß1-adrenoceptors, reduced neuronal firing by a mean (SD) 37.3% (5.5%) or 40% (6.3%), respectively, the latter via SK potassium channel opening. An L-type calcium channel blocker or ß1-adrenoceptor antagonist protected working memory from stress. Conclusions and Relevance: The layer III pyramidal cells in the dorsolateral prefrontal cortex especially vulnerable in cognitive disorders differentially express calbindin and a constellation of calcium-related proteins including L-type calcium channels Cav1.2 (CACNA1C), GluN2B-NMDA receptors (GRIN2B), and SK3 potassium channels (KCNN3), which influence memory-related neuronal firing. The finding that either inadequate or excessive L-type calcium channel activation reduced neuronal firing explains why either loss- or gain-of-function variants in CACNA1C were associated with increased risk of cognitive disorders. The selective expression of calbindin in these pyramidal cells highlights the importance of regulatory mechanisms in neurons with high calcium signaling, consistent with Alzheimer tau pathology emerging when calbindin is lost with age and/or inflammation.

Sibling chimerism among microglia in marmosets.

Chimerism happens rarely among most mammals but is common in marmosets and tamarins, a result of fraternal twin or triplet birth patterns in which in utero connected circulatory systems (through which stem cells transit) lead to persistent blood chimerism (12-80%) throughout life. The presence of Y-chromosome DNA sequences in other organs of female marmosets has long suggested that chimerism might also affect these organs. However, a longstanding question is whether this chimerism is driven by blood-derived cells or involves contributions from other cell types. To address this question, we analyzed single-cell RNA-seq data from blood, liver, kidney and multiple brain regions across a number of marmosets, using transcribed single nucleotide polymorphisms (SNPs) to identify cells with the sibling's genome in various cell types within these tissues. Sibling-derived chimerism in all tissues arose entirely from cells of hematopoietic origin (i.e., myeloid and lymphoid lineages). In brain tissue this was reflected as sibling-derived chimerism among microglia (20-52%) and macrophages (18-64%) but not among other resident cell types (i.e., neurons, glia or ependymal cells). The percentage of microglia that were sibling-derived showed significant variation across brain regions, even within individual animals, likely reflecting distinct responses by siblings' microglia to local recruitment or proliferation cues or, potentially, distinct clonal expansion histories in different brain areas. In the animals and tissues we analyzed, microglial gene expression profiles bore a much stronger relationship to local/host context than to sibling genetic differences. Naturally occurring marmoset chimerism will provide new ways to understand the effects of genes, mutations and brain contexts on microglial biology and to distinguish between effects of microglia and other cell types on brain phenotypes.

A marmoset brain cell census reveals regional specialization of cellular identities.

The mammalian brain is composed of many brain structures, each with its own ontogenetic and developmental history. We used single-nucleus RNA sequencing to sample over 2.4 million brain cells across 18 locations in the common marmoset, a New World monkey primed for genetic engineering, and examined gene expression patterns of cell types within and across brain structures. The adult transcriptomic identity of most neuronal types is shaped more by developmental origin than by neurotransmitter signaling repertoire. Quantitative mapping of GABAergic types with single-molecule FISH (smFISH) reveals that interneurons in the striatum and neocortex follow distinct spatial principles, and that lateral prefrontal and other higher-order cortical association areas are distinguished by high proportions of VIP+ neurons. We use cell type-specific enhancers to drive AAV-GFP and reconstruct the morphologies of molecularly resolved interneuron types in neocortex and striatum. Our analyses highlight how lineage, local context, and functional class contribute to the transcriptional identity and biodistribution of primate brain cell types.

Comparative single cell epigenomic analysis of gene regulatory programs in the rodent and primate neocortex.

Sequence divergence of cis- regulatory elements drives species-specific traits, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains to be elucidated. We investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset, and mouse with single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome, and chromosomal conformation profiles from a total of over 180,000 cells. For each modality, we determined species-specific, divergent, and conserved gene expression and epigenetic features at multiple levels. We find that cell type-specific gene expression evolves more rapidly than broadly expressed genes and that epigenetic status at distal candidate cis -regulatory elements (cCREs) evolves faster than promoters. Strikingly, transposable elements (TEs) contribute to nearly 80% of the human-specific cCREs in cortical cells. Through machine learning, we develop sequence-based predictors of cCREs in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Lastly, we show that epigenetic conservation combined with sequence similarity helps uncover functional cis -regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.