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.

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.

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 transcriptomics reveals human-specific cortical features.

The cognitive abilities of humans are distinctive among primates, but their molecular and cellular substrates are poorly understood. We used comparative single-nucleus transcriptomics to analyze samples of the middle temporal gyrus (MTG) from adult humans, chimpanzees, gorillas, rhesus macaques, and common marmosets to understand human-specific features of the neocortex. Human, chimpanzee, and gorilla MTG showed highly similar cell-type composition and laminar organization as well as a large shift in proportions of deep-layer intratelencephalic-projecting neurons compared with macaque and marmoset MTG. Microglia, astrocytes, and oligodendrocytes had more-divergent expression across species compared with neurons or oligodendrocyte precursor cells, and neuronal expression diverged more rapidly on the human lineage. Only a few hundred genes showed human-specific patterning, suggesting that relatively few cellular and molecular changes distinctively define adult human cortical structure.

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.

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.

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.

Gene expression links functional networks across cortex and striatum.

The human brain is comprised of a complex web of functional networks that link anatomically distinct regions. However, the biological mechanisms supporting network organization remain elusive, particularly across cortical and subcortical territories with vastly divergent cellular and molecular properties. Here, using human and primate brain transcriptional atlases, we demonstrate that spatial patterns of gene expression show strong correspondence with limbic and somato/motor cortico-striatal functional networks. Network-associated expression is consistent across independent human datasets and evolutionarily conserved in non-human primates. Genes preferentially expressed within the limbic network (encompassing nucleus accumbens, orbital/ventromedial prefrontal cortex, and temporal pole) relate to risk for psychiatric illness, chloride channel complexes, and markers of somatostatin neurons. Somato/motor associated genes are enriched for oligodendrocytes and markers of parvalbumin neurons. These analyses indicate that parallel cortico-striatal processing channels possess dissociable genetic signatures that recapitulate distributed functional networks, and nominate molecular mechanisms supporting cortico-striatal circuitry in health and disease.

Gradients in cytoarchitectural landscapes of the isocortex: Diprotodont marsupials in comparison to eutherian mammals.

Although it has been claimed that marsupials possess a lower density of isocortical neurons compared with other mammals, little is known about cross-cortical variation in neuron distributions in this diverse taxonomic group. We quantified upper-layer (layers II-IV) and lower-layer (layers V-VI) neuron numbers per unit of cortical surface area in three diprotodont marsupial species (two macropodiformes, the red kangaroo and the parma wallaby, and a vombatiform, the koala) and compared these results to eutherian mammals (e.g., xenarthrans, rodents, primates). In contrast to the notion that the marsupial isocortex contains a low density of neurons, we found that neuron numbers per unit of cortical surface area in several marsupial species overlap with those found in eutherian mammals. Furthermore, neuron numbers vary systematically across the isocortex of the marsupial mammals examined. Neuron numbers under a unit of cortical surface area are low toward the frontal cortex and high toward the caudo-medial (occipital) pole. Upper-layer neurons (i.e., layers II-IV) account for most of the variation in neuron numbers across the isocortex. The variation in neuron numbers across the rostral to the caudal pole resembles primates. These findings suggest that diprotodont marsupials and eutherian mammals share a similar cortical architecture despite their distant evolutionary divergence.

Gradients of Connectivity in the Cerebral Cortex.

The human neocortex is organized with distributed networks that connect distant regions together, but what determines their spatial layout? A recent study sheds light on the topological placement of regions along the cortical surface in relation to gradients of connectivity in both humans and macaques.

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.

Reconfigurable task-dependent functional coupling modes cluster around a core functional architecture.

Functional coupling across distributed brain regions varies across task contexts, yet there are stable features. To better understand the range and central tendencies of network configurations, coupling patterns were explored using functional MRI (fMRI) across 14 distinct continuously performed task states ranging from passive fixation to increasingly demanding classification tasks. Mean global correlation profiles across the cortex ranged from 0.69 to 0.82 between task states. Network configurations from both passive fixation and classification tasks similarly predicted task coactivation patterns estimated from meta-analysis of the literature. Thus, even across markedly different task states, central tendencies dominate the coupling configurations. Beyond these shared components, distinct task states displayed significant differences in coupling patterns in response to their varied demands. One possibility is that anatomical connectivity provides constraints that act as attractors pulling network configurations towards a limited number of robust states. Reconfigurable coupling modes emerge as significant modifications to a core functional architecture.

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.