The human neocortex underlies higher cognition and is the engine of complex thought. Yet our understanding of its neuronal diversity is limited by sparse access to tissue, inconsistent sampling across studies, and a lack of multiple modality data. Although single-cell transcriptomic taxonomies are an important framework for characterizing cell type diversity, transcriptomic information alone cannot reveal the cellular properties that define neuronal computations. To address this, we performed Patch-seq, a method for collecting Morphology, Electrophysiology, and Transcriptomic data from a single neuron. We focused on glutamatergic, neocortical, excitatory neurons, the principal long-range projecting neurons of the cortex, and systematically integrated their morphoelectric features with transcriptomic identity. In combination with spatial transcriptomic data, we interrogated 39 of 42 transcriptomically-defined neuron types with a layer-centric perspective. Morphoelectric properties, such as cortical depth, apical dendrite structure, and excitability clearly distinguish transcriptomic subclasses and support many finer transcriptomic types. Morphoelectric properties are influenced by spatial location in supragranular layers, while deeper layers exhibit greater heterogeneity. Cross-species comparisons reveal conserved subclass organization but pronounced differences in apical dendrite arborization between mouse and human, and surprising similarities between human and macaque. Together, these datasets provide a unified multimodal reference that advances our understanding of human cortical circuitry and establishes a foundation for experimental and computational studies of human brain function and disease.

Functional genomics studies have provided critical insights into cell type-specific gene regulatory programs, but to date most have been conducted in wild-type tissues or cell cultures. Here, we present a gene expression functional atlas across the mouse brain. We use an enhanced in vivo Perturb-seq platform to analyze transcriptome-wide responses to loss of 1,947 disease-associated genes, profiling over 7.7 million cells spanning major brain regions and neuronal populations. We find striking cell-type-specific essentiality and transcriptional programs and show that closely related disease genes such as two NMDA receptor subunits can drive opposing transcriptional programs. Together, this work reveals insights into the genetics and mechanisms of neurodevelopmental, psychiatric, and neurodegenerative diseases in vivo, paving the way for the design of future genetic medicine.

The external globus pallidus (GPe) is traditionally viewed as a homogeneous relay in the indirect basal ganglia pathway that broadly suppresses movement. Using whole-brain anterograde axon mapping, rabies tracing, single-neuron reconstruction, and single-nucleus RNA sequencing, we reveal that the GPe is a major basal ganglia output nucleus, composed of anatomically and molecularly distinct populations. Besides the neurons with canonical projections, and projections to striatum and cortex, the GPe contains specific neuronal populations that target the thalamus and brainstem directly, including the parafascicular thalamus (GPePf) and the pedunculopontine nucleus (GPePPN). These projection-defined subpopulations exhibit distinct behavioral functions. GPe-PPN neurons are selectively suppressed at locomotor onset, and bidirectional manipulation of this pathway is sufficient to promote or suppress locomotion. Notably, D2-MSN stimulation upstream of these neurons evokes locomotion, and co-activation of GPePPN pathway blocks this D2-MSN-evoked locomotion, demonstrating a disinhibitory circuit for action stemming from D2-MSNs. In contrast, GPePf neurons are not engaged during locomotion but are selectively recruited during skilled forelimb actions, and their activation disrupts forelimb movements without affecting locomotion. Together, these findings establish the GPe as a basal ganglia output hub composed of modules that mediate distinct behaviors. This organization revises canonical models of the indirect pathway by demonstrating that D2-MSNs can also facilitate, rather than only suppress, movement depending on the downstream GPe output channels they engage.

The basal ganglia comprise interconnected subcortical nuclei essential for motor control, learning, emotion and cognition, yet how cell types relate to their circuit organization remains elusive. Here we show that spatial patterns of transcriptomic cell types in the basal ganglia and thalamic parafascicular nucleus relate to module-organized cortical inputs in the mouse. By co-registering genoarchitectural and connectivity datasets to the common coordinate framework, we delineate distinct three-dimensional subdivisions of these nuclei. Analyses of anterograde and retrograde viral tracing data and single-neuron reconstructions reveal that each subdivision receives convergent, modular and cell-type-specific cortical and subcortical inputs. Caudoputamen subdivisions primarily receive layer 5 intratelencephalic inputs from the entire isocortex, whereas smaller basal ganglia subdivisions receive ipsilateral layer 5 extratelencephalic inputs from limited cortical areas. The parafascicular nucleus subdivisions are connected with specific cortical modules and basal ganglia subdivisions. Our results suggest that combinatorial gene expression underlies a global map of parallel, modular and cell-type-specific basal ganglia and parafascicular nucleus subnetworks.

Precise manipulation of neurons across cortical layers requires optical approaches that maintain single-cell specificity deep within scattering brain tissue. Two-photon optogenetics provides targeted photostimulation in vivo and has expanded experimental access to defined neuronal populations beyond the limits of single-photon excitation. However, its application remains largely confined to the upper cortical layers (∼250–300 µm) in mice due to scattering. Here we introduce three-photon temporally focused computer-generated holography (3P TF-CGH) for precise optogenetic activation beyond this depth limit. We characterize 3P excitation of multiple excitatory and inhibitory opsins in organotypic slices, demonstrating cubic power dependence, efficient photocurrents and preserved temporal fidelity. We demonstrated that temporally focused holography maintains tight axial confinement in scattering tissue and minimizes superficial out-of-focus excitation under conditions required for deep targeting. In vivo, we combined 3P holographic stimulation at 1700 nm with simultaneous 3P calcium imaging at 1300 nm to achieve reliable neuronal activation across cortical layers down to 800 µm. Photostimulation remained stable across repeated trials without detectable physiological perturbation. By extending cell-resolved optogenetic control well beyond the established depth limits of 2P approaches, 3P holographic optogenetics enables non-invasive, all-optical interrogation of deep cortical circuits.

Understanding how cell types are organized across the central nervous system (CNS) is key to uncovering neural function. Here, we integrate single-nucleus Multiome (RNA+ATAC) sequencing, spatial transcriptomics, and computational analyses to map conserved cell type signatures in the adult mouse brainstem and spinal cord. We identify a shared core of neuronal and non-neuronal cell types, alongside region-specific specializations reflecting distinct functions. Spatial data reveal conserved cellular niches across the brainstem-spinal cord boundary, indicating a continuous organizational logic. Cross-region comparisons uncover recurrent gene expression modules and signaling programs that may support shared circuit features. Chromatin accessibility profiling highlights cell-type-specific regulatory programs and implicates Hox transcription factors in positional identity. Notably, cell-type and positional identities are largely orthogonal, with varying regional influence across neuronal classes: motor neurons show strong positional coupling, whereas glutamatergic and GABAergic interneurons show minimal entrainment. This work provides a reference for the shared molecular architecture of these CNS regions.

Patterns of brain activity moving in waves occur across brain regions and species, yet their spatial organization, anatomical basis, and brain-wide distribution remain unclear. Using cortex-wide imaging and electrophysiology in awake mice, we revealed a prominent wave motif across spatial scales. Waves frequently formed rotational patterns centered on somatosensory cortex and sweeping across somatotopic maps. Axonal architecture within sensory cortex exhibited a matching circular arrangement. Rotating waves were mirrored between hemispheres and between sensory and motor cortex and were coordinated with subcortical spiking. Bilaterally cutting the circular circuitry diminished rotating waves. Rotating waves were modulated across behavioral states, evoked by sensory inputs, and recruited during correct visuomotor performance. These results establish that rotating waves are sculpted by axonal architecture across diverse brain systems and behavioral contexts.

The basal ganglia (BG) are a set of topographically organized, interconnected structures that are pivotal for regulating volitional movement and other aspects of cognitive, motivational, and affective behavior. Recently generated taxonomies of transcriptomically-defined cell types (T-types) have revealed both fine-grained distinctions in gene expression between neurons in these structures as well as continuous transcriptomic variation across similar T-types1-9, which are both related to location within a structure. However, it remains unclear to what extent these and other cellular properties co-vary with each other. Therefore, we performed Patch-seq experiments10,11 on over 900 neurons in mouse brain slices from BG to provide an integrated view of the co-variation between gene expression, location, physiology, and morphology measured from the same neurons. Medium spiny neurons (MSNs) from both the direct and indirect pathways across the dorsal and ventral striatum follow a gene expression gradient that varies in a dorsolateral to ventromedial direction; we find that this gradient also corresponds with systematic differences in action potential kinetics and dendritic arborization. Our analysis also characterizes additional multimodal dimensions of MSN variation, such as those between direct and indirect pathway neurons and between matrix and striosome neurons. Furthermore, through comparison with Patch-seq data from macaque, we demonstrate that the relationship between the transcriptomic/spatial gradient and electrophysiological and morphological properties is conserved across these two species. We also find that properties of striatal interneurons, such as action potential kinetics, vary across the striatum in a manner consistent with the MSN gradient. Outside the striatum, our multimodal Patch-seq dataset from the globus pallidus, subthalamic nucleus, and substantia nigra enabled us to characterize transcriptomically-defined types and link them to prior descriptions of cell types in these structures. Finally, we examined to what extent the MSN transcriptomic/spatial gradient persisted across different stages of the BG circuit by comparing Patch-seq neurons to reconstructed whole-neuron morphologies and the topography of their interareal projections, finding that the gradient is better preserved in GPe and GPi compared to SNr. Our study links transcriptomic variation across T-types in mouse BG to spatial localization and phenotypic differences at the level of individual cells, improving our understanding of cell type architecture in topographically organized circuits of the brain.

Single-cell mapping methods convert raw, heterogeneous single-cell datasets into interpretable and comparable representations of biological identity. As reference cell-type taxonomies mature, mapping new datasets to shared references has become a central strategy for enabling cross-study integration, reproducible annotation, and cumulative biological knowledge. Here we present MapMyCells, an open-source framework designed to align diverse single-cell omics datasets to hierarchical reference taxonomies with minimal preprocessing. MapMyCells provides out-of-the-box support for an expanding set of high-quality brain cell-type references generated by the Allen Institute for Brain Science, the BRAIN Initiative, and the Seattle Alzheimer's Disease Brain Cell Atlas, including whole-brain mouse and human atlases, aging and Alzheimer's disease cohorts, and a cross-species consensus taxonomy initially focused on the basal ganglia. MapMyCells enables efficient mapping of hundreds of thousands of cells on standard workstations without specialized hardware, providing a deterministic, scalable, and modality-agnostic approach that is robust across species and molecular assays. The framework produces interpretable confidence metrics and quantitative summaries of mapping performance, allowing users to evaluate assignment precision and accuracy. We demonstrate the mapping of unlabeled transcriptomic, epigenomic, and spatial datasets to reference taxonomies and describe a general workflow for preparing arbitrary hierarchical taxonomies for reference-based mapping. As the ecosystem of single-cell reference atlases expands, MapMyCells offers a practical and reproducible solution for community-scale cell-type annotation and cross-dataset integration, supporting the development of unified and extensible brain cell atlases.

A major goal of the BRAIN Initiative Cell Atlas Network (BICAN) is to create a suite of foundational reference cell atlases and associated standards for human and non-human primate brains. Central to this goal is the creation of cross-species harmonized cellular taxonomies and structural parcellations with formal ontologies that can be mapped into 3D reference frameworks bridging neuroimaging and cellular and histological resolutions. We describe here an iterative approach, focused initially on the basal ganglia, to co-create structural and cellular ontologies in human, macaque and marmoset brains, including a Harmonized Ontology of Mammalian Brain Anatomy (HOMBA), and to map and refine structural parcellations into neuroimaging-based common coordinate frameworks. These references provide the framework for documenting and mapping all experimental sampling in BICAN, allowing analyses of cellular and molecular variation as a function of topographic position, and enabling comparisons of cellular, molecular and neuroimaging-based functional variation within and between primate species.

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.

Spatial transcriptomics enables the precise mapping of gene expression patterns within tissue architecture, offering unprecedented insights into cellular interactions, tissue heterogeneity, and disease pathology that are unattainable with traditional transcriptomic approaches. We present a tool for processing spatial transcriptomics data, SCALPEL (Spatial Cell Analysis, Labeling, Processing, and Expression Linking). SCALPEL is specifically designed to support the analysis of large, atlas-level datasets. Our new workflow features advanced 3D segmentation optimized for dense and heterogeneous tissues, refined filtering criteria, and transcriptome-based doublet detection to remove low-quality or artifactual cells. Cell type label transfer from existing taxonomies is further improved through updated filtering thresholds. Spatial domain detection is incorporated to capture local transcriptomic organization, and tissue sections are registered to the Allen Mouse Brain Common Coordinate Framework version 3 (CCFv3) for precise anatomical alignment. Genome-wide expression imputation from single-cell RNA-sequencing (scRNAseq) further enriches the dataset. Crucially, we benchmark the performance of this updated pipeline against a previously published version of our whole-mouse-brain (WMB) dataset (Yao et al., 2023b), demonstrating substantial improvements in cell number, expression profile clarity, and spatial registration. These advances provide a robust foundation for downstream spatial analyses and set a new standard for large-scale spatial transcriptomics studies.

The mammalian basal ganglia (BG) orchestrate motor, cognitive, and affective functions, yet cell type-specific genetic access remains limited, especially beyond rodents. Key structures implicated in movement and psychiatric disorders, including pallidum, subthalamic nucleus, and dopaminergic midbrain, lack scalable tools for cross-species targeting. Here, we present a comprehensive enhancer-AAV library enabling selective labeling and manipulation of major BG neuronal populations: striatal projection neuron subtypes, pallidal and subthalamic neurons, and midbrain dopaminergic and GABAergic populations. Using an evolutionarily informed discovery pipeline, we identified enhancers targeting canonical, non-canonical, and disease-relevant cell types, with validation demonstrating robust cross-species conservation of specificity between mouse and macaque. Computational modeling revealed sequence features predictive of in vivo performance, including motif grammar, chromatin accessibility, and evolutionary conservation, and identified distinct regulatory architectures across glial, projection, and interneuron lineages. This work establishes a comprehensive cross-species viral toolkit for the BG, unlocking previously inaccessible cell types for circuit dissection.

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 mammalian brain consists of diverse neuron types with various functions. Recent single-cell RNA sequencing approaches have led to a whole-brain taxonomy of transcriptomically defined cell types1. Patch-seq experiments augment these cell-type descriptions by linking transcriptomic profiles with local morphological and electrophysiological properties2-7. However, linking transcriptomic identities to long-range axonal projections remains a major unresolved challenge. Here, to address this, we collected two datasets from the mouse visual cortex consisting of: (1) 1,528 excitatory Patch-seq neurons, with local morphological, electrophysiological and transcriptomic data collected from each cell, and (2) 341 excitatory, whole-neuron morphologies. From the Patch-seq data, we defined 17 morphoelectric-transcriptomic types and built a multistep classifier to integrate cell-type assignments with whole-neuron morphology and interrogate cross-modality relationships. We find that transcriptomic variation within and across morphoelectric-transcriptomic types corresponds with morphological and electrophysiological phenotypes. In addition, these gene expression patterns, along with the anatomical location of the cell, can be used to predict projection targets of individual neurons. We observed novel multimodal cell-type signatures for layer 5 intratelencephalic and extratelencephalic neurons and shed new light on their axonal circuitry, including interhemispheric intratelencephalic projections. With this approach, we establish a comprehensive, integrated taxonomy of cortical, excitatory neuron types, and create a system for high-dimensional cell-type classification that can be extended to the whole brain and potentially across species.

The complexity of the mammalian brain’s vast population of interconnected neurons poses a formidable challenge to elucidate its underlying mechanisms of coordination and computation. A key step forward will be technologies that can perform large-scale, cellular-resolution monitoring and interrogation of distributed brain circuit activity in behaving animals. Here, we present an all-optical strategy for precise optogenetic activity control of ∼103 neurons and simultaneous activity monitoring of ∼104 neurons within and across areas of mouse cortex—an order-of-magnitude leap beyond previous capabilities. Tracking population responses following delivery of precisely-defined widely-distributed activity patterns to the visual cortex of awake mice, we were surprised to identify neurons robustly responsive to stimulation of diverse ensembles, defying conventional like-to-like wiring rules. These cells were primarily deep L2/3 somatostatin-positive (SST) interneurons with functional properties distinct from other SST neurons, and appeared to play a role in brain dynamics that could only have been identified through broad cellular-resolution circuit interrogation. Our work reveals the value of measuring large-scale circuit-dynamical properties of functionally-resolved single cells, beyond genetic and anatomical classification, to define and explore the roles of cell types in brain function. ### Competing Interest Statement K.D. is a founder and scientific advisor for Maplight Therapeutics and Stellaromics and a scientific advisor to RedTree LLC and Modulight.

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.

Task learning involves learning associations between stimuli and outcomes and storing these relationships in memory. While this information can be reliably decoded from population activity, individual neurons encoding this representation can drift over time. The circuit or molecular mechanisms underlying this drift and its role in learning are unclear. We performed two-photon calcium imaging in the perirhinal cortex during task training. Using post hoc spatial transcriptomics, we measured immediate-early gene (IEG) expression and assigned monitored neurons to excitatory or inhibitory subtypes. We discovered an IEG-defined network spanning multiple subtypes that form stimulus-outcome associations. Targeted deletion of brain-derived neurotrophic factor in the perirhinal cortex disrupted IEG expression and impaired task learning. Representational drift slowed with prolonged training. Pre-existing representations were strengthened while stimulus-reward associations failed to form. Our findings reveal the cell types and molecules regulating long-term network stability that is permissive for task learning and memory allocation.