The basal ganglia are evolutionary ancient subcortical nuclei that form interconnected loops with the neocortex and limbic system to regulate movement, learning, habit formation, emotion, and motivation. Their dysfunction contributes to major neurological and psychiatric disorders, yet most cellular-level insights derive from rodent studies, leaving knowledge gaps in humans and translationally relevant primate species. To address this, we generated multi-modal Patch-seq data linking transcriptomic identity with morphological and electrophysiological properties in macaque striatum, the input nucleus of the basal ganglia. We found underappreciated diversity among medium spiny neurons, including non-canonical types, and variation aligned with functional gradients. Interneurons also exhibited spatial variation and even greater morphoelectric diversity, highlighting their functional modularity. Despite broad evolutionary conservation, we identified primate-specific features and key differences from rodent striatal neurons. By integrating molecular classification with cellular properties that shape network function, our findings provide insights into the functional organization of the primate striatum.

The spinal cord contains evolutionarily conserved cell types critical for motor function, sensory processing, and autonomic regulation, many of which are implicated in diverse neurological diseases and injuries. Yet the field lacks a comprehensive molecular characterization of cellular diversity in human, macaque, and mouse spinal cord. Here, we present a unified, cross-species cell type atlas based on the integration of single-nucleus gene expression, chromatin accessibility, and spatial transcriptomic data from segments within cervical, thoracic, lumbar, and sacral regions, including motor neurons (MNs) sampled across the entire rostro-caudal axis of the macaque spinal cord. Leveraging the spatial distributions of our molecularly defined cell types, we generated a cell type-guided anatomical map of spinal cord laminae and nuclei. We identified both conserved and species-specific cellular features, including gene expression patterns across distinct MN subtypes in the primate spinal cord. Cross-species cis-regulatory analysis and deep learning sequence models dissected the enhancer logic underlying viral targeting, uncovering conserved transcription factor grammar encoding cellular identity. Together, these results establish a unifying molecular and anatomical taxonomy of spinal cord cell types across species.

The basal ganglia (BG) are conserved brain regions essential for motor control, learning, emotion, and cognition, and are implicated in neurological and psychiatric disease. Yet a unified cross-species taxonomy of BG cell types is lacking, limiting translation of BG circuit mechanisms, interpretation of human genetic risk, and development of cell type-targeted tools. We present a multiomic consensus atlas of 1.8 million nuclei from human, macaque, and marmoset spanning eight BG structures. Integrating cross-species gene expression, open chromatin, and spatial profiling enables definition of conserved and divergent cell types. Alignment to existing mouse and human atlases identifies 61 homologous cell types conserved over 80 million years. We identify a STRd D2 StrioMat Hybrid medium spiny neuron (MSN) type with molecular, electrophysiological, and morphological features that clarify hybrid MSN identities. Comparative cis-regulatory analysis reveals conserved sequence grammars that encode cell identity and inform viral targeting strategies, providing a foundational resource for BG evolution, function, and disease.

We present an enhancer-AAV toolbox for accessing and perturbing striatal cell types and circuits. Best-in-class vectors were curated for accessing major striatal neuron populations including medium spiny neurons (MSNs), direct- and indirect-pathway MSNs, Sst-Chodl, Pvalb-Pthlh, and cholinergic interneurons. Specificity was evaluated by multiple modes of molecular validation, by three different routes of virus delivery, and with diverse transgene cargos. Importantly, we provide detailed information necessary to achieve reliable cell-type-specific labeling under different experimental contexts. We demonstrate direct pathway circuit-selective optogenetic perturbation of behavior and multiplex labeling of striatal interneuron types for targeted analysis of cellular features. Lastly, we show conserved in vivo activity for exemplary MSN enhancers in rats and macaques. This collection of striatal enhancer AAVs offers greater versatility compared to available transgenic lines and can readily be applied for cell type and circuit studies in diverse mammalian species beyond the mouse model.

Single-cell transcriptomics has enhanced our understanding of the brain's cellular composition. Biologists now analyze complex datasets to explore how marker genes influence biological processes, genetic variations, and phenotypic traits. A challenge is comparing these datasets across species to detect subtle differences or similarities to evolutionary development. Here, we present Cytosplore EvoViewer to facilitate examining relationships between transcriptomic datasets across species, simplifying the analysis of marker gene regulation and its impact on biological functions and integrating these findings with prior evolutionary knowledge or species-specific traits. We conducted a design study, including domain analysis, implementation of the results into Cytosplore EvoViewer, and an expert evaluation. Cytosplore EvoViewer offers valuable insights into genetic variations and evolutionary dynamics, helping to understand the diversity and the unity within diversity across species and their evolutionary development. The Cytosplore EvoViewer installer application can be downloaded from the Cytosplore Viewer website(1), and its source code is available on the ManiVault Studio GitHub(2).

The mammalian cortex is comprised of cells classified into types according to shared properties. Defining the contribution of each cell type to the processes guided by the cortex is essential for understanding its function in health and disease. We use transcriptomic and epigenomic cortical cell-type taxonomies from mouse and human to define marker genes and putative enhancers and create a large toolkit of transgenic lines and enhancer adeno-associated viruses (AAVs) for selective targeting of cortical cell populations. We report creation and evaluation of fifteen transgenic driver lines, two reporter lines, and >1,000 different enhancer AAV vectors covering most subclasses of cortical cells. The tools reported here have been made publicly available, and along with the scaled process of tool creation, evaluation, and modification, they will enable diverse experimental strategies toward understanding mammalian cortex and brain function.

The mammalian cortex is comprised of cells classified into types according to shared properties. Defining the contribution of each cell type to the processes guided by the cortex is essential for understanding its function in health and disease. We used transcriptomic and epigenomic cortical cell type taxonomies from mouse and human to define marker genes and putative enhancers and created a large toolkit of transgenic lines and enhancer AAVs for selective targeting of cortical cell populations. We report evaluation of fifteen new transgenic driver lines, two new reporter lines, and >800 different enhancer AAVs covering most subclasses of cortical cells. The tools reported here as well as the scaled process of tool creation and modification enable diverse experimental strategies towards understanding mammalian cortex and brain function.

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Evolutionary constraint and acceleration are powerful, cell-type agnostic measures of functional importance. Previous studies in mammals were limited by species number and reliance on human-referenced alignments. We explore the evolution of placental mammals, including humans, through reference-free whole-genome alignment of 240 species and protein-coding alignments for 428 species. We estimate 10.7% of the human genome is evolutionarily constrained. We resolve constraint to single nucleotides, pinpointing functional positions, and refine and expand by over seven-fold the catalog of ultraconserved elements. Overall, 48.5% of constrained bases are as yet unannotated, suggesting yet-to-be-discovered functional importance. Using species-level phenotypes and an updated phylogeny, we associate coding and regulatory variation with olfaction and hibernation. Focusing on biodiversity conservation, we identify genomic metrics that predict species at risk of extinction.

Organismal behavior, with its tremendous complexity and diversity, is generated by numerous physiological systems acting in coordination. Understanding how these systems evolve to support differences in behavior within and among species is a longstanding goal in biology that has captured the imagination of researchers who work on a multitude of taxa, including humans. Of particular importance are the physiological determinants of behavioral evolution, which are sometimes overlooked because we lack a robust conceptual framework to study mechanisms underlying adaptation and diversification of behavior. Here, we discuss a framework for such an analysis that applies a "systems view" to our understanding of behavioral control. This approach involves linking separate models that consider behavior and physiology as their own networks into a singular vertically integrated behavioral control system. In doing so, hormones commonly stand out as the links, or edges, among nodes within this system. To ground our discussion, we focus on studies of manakins (Pipridae), a family of Neotropical birds. These species have numerous physiological and endocrine specializations that support their elaborate reproductive displays. As a result, manakins provide a useful example to help imagine and visualize the way systems concepts can inform our appreciation of behavioral evolution. In particular, manakins help clarify how connectedness among physiological systems-which is maintained through endocrine signaling-potentiate and/or constrain the evolution of complex behavior to yield behavioral differences across taxa. Ultimately, we hope this review will continue to stimulate thought, discussion, and the emergence of research focused on integrated phenotypes in behavioral ecology and endocrinology.

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth’s vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Background Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high. Results We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin. Conclusion The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

Genetic studies are rapidly identifying non-protein-coding human disease-associated loci. Understanding the regulatory mechanisms underlying these loci remains a challenge because the causal variants and the tissues in which they act are often unclear. Massively parallel reporter assays (MPRAs) have the potential to link differences in genome sequence, including genetic variants, to tissue-specific regulatory function. Although MPRA and similar technologies have been widely adopted in cell culture, there have been several barriers to widespread use in animals. We overcome these challenges with a new whole-animal MPRA (WhAMPRA), where systemic intravenous AAV effectively transduces the plasmid MPRA library to mouse tissues. Our WhAMPRA approach revealed models of tissue-specific regulation that generally match machine learning model predictions. In addition, we measured the regulatory effects of disrupting MEF2C transcription factor binding sites and impacts of late onset Alzheimer’s disease-associated genetic variations. Overall, our WhAMPRA technology simultaneously determines the transcriptional functions of hundreds of enhancers in vivo across multiple tissues.### Competing Interest StatementThe authors have declared no competing interest.

Vocal learning is a skilled motor behavior observed in several mammalian and avian species and is critical for human speech. While convergent gene expression patterns have highlighted similar primary motor and striatal pathways for vocal imitation in songbirds and humans, the extent of molecular and circuit convergence remains unresolved. Here we profiled the four principal song nuclei of the zebra finch (HVC, LMAN, RA, Area X) and their surrounding brain regions using RNA-Seq and compared them with specialized markers we identified for human speech brain regions. Expanding previous work, both songbird RA and HVC exhibited convergent specialized gene expression of ∼350 genes with human laryngeal sensorimotor cortex. The songbird HVC RA intratelencephalic (IT) neurons were the predominant cell type that was convergent with human, specifically layer 2/3 IT neurons, while the songbird RA extratelencephalic (ET) projection neurons exhibited convergent expression with human layer 5 ET projection neurons. The molecular specializations of both songbird LMAN and human Broca’s area were more unique to each species. These findings demonstrate the extent of convergent molecular specializations in distantly related species for vocal imitation and emphasize important evolutionary constraints for this complex trait. One-Sentence Summary Our data provide hundreds of candidate genes to study the molecular basis and evolution of song and speech across species.

Recent discoveries of extreme cellular diversity in the brain warrant rapid development of technologies to access specific cell populations within heterogeneous tissue. Available approaches for engineering-targeted technologies for new neuron subtypes are low yield, involving intensive transgenic strain or virus screening. Here, we present Specific Nuclear-Anchored Independent Labeling (SNAIL), an improved virus-based strategy for cell labeling and nuclear isolation from heterogeneous tissue. SNAIL works by leveraging machine learning and other computational approaches to identify DNA sequence features that confer cell type-specific gene activation and then make a probe that drives an affinity purification-compatible reporter gene. As a proof of concept, we designed and validated two novel SNAIL probes that target parvalbumin-expressing (PV+) neurons. Nuclear isolation using SNAIL in wild-type mice is sufficient to capture characteristic open chromatin features of PV+ neurons in the cortex, striatum, and external globus pallidus. The SNAIL framework also has high utility for multispecies cell probe engineering; expression from a mouse PV+ SNAIL enhancer sequence was enriched in PV+ neurons of the macaque cortex. Expansion of this technology has broad applications in cell type-specific observation, manipulation, and therapeutics across species and disease models.

Recent large genome-wide association studies have identified multiple confident risk loci linked to addiction-associated behavioral traits. Most genetic variants linked to addiction-associated traits lie in noncoding regions of the genome, likely disruptingcis-regulatory element (CRE) function. CREs tend to be highly cell type-specific and may contribute to the functional development of the neural circuits underlying addiction. Yet, a systematic approach for predicting the impact of risk variants on the CREs of specific cell populations is lacking. To dissect the cell types and brain regions underlying addiction-associated traits, we applied stratified linkage disequilibrium score regression to compare genome-wide association studies to genomic regions collected from human and mouse assays for open chromatin, which is associated with CRE activity. We found enrichment of addiction-associated variants in putative CREs marked by open chromatin in neuronal (NeuN+) nuclei collected from multiple prefrontal cortical areas and striatal regions known to play major roles in reward and addiction. To further dissect the cell type-specific basis of addiction-associated traits, we also identified enrichments in human orthologs of open chromatin regions of female and male mouse neuronal subtypes: cortical excitatory, D1, D2, and PV. Last, we developed machine learning models to predict mouse cell type-specific open chromatin, enabling us to further categorize human NeuN+open chromatin regions into cortical excitatory or striatal D1 and D2 neurons and predict the functional impact of addiction-associated genetic variants. Our results suggest that different neuronal subtypes within the reward system play distinct roles in the variety of traits that contribute to addiction.SIGNIFICANCE STATEMENTWe combine statistical genetic and machine learning techniques to find that the predisposition to for nicotine, alcohol, and cannabis use behaviors can be partially explained by genetic variants in conserved regulatory elements within specific brain regions and neuronal subtypes of the reward system. Our computational framework can flexibly integrate open chromatin data across species to screen for putative causal variants in a cell type- and tissue-specific manner for numerous complex traits.