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.
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.
Proper brain function requires the assembly and function of diverse populations of neurons and glia. Single cell gene expression studies have mostly focused on characterization of neuronal cell diversity; however, recent studies have revealed substantial diversity of glial cells, particularly astrocytes. To better understand glial cell types and their roles in neurobiology, we built a new suite of adeno-associated viral (AAV)-based genetic tools to enable genetic access to astrocytes and oligodendrocytes. These oligodendrocyte and astrocyte enhancer-AAVs are highly specific (usually > 95% cell type specificity) with variable expression levels, and our astrocyte enhancer-AAVs show multiple distinct expression patterns reflecting the spatial distribution of astrocyte cell types. To provide the best glial-specific functional tools, several enhancer-AAVs were: optimized for higher expression levels, shown to be functional and specific in rat and macaque, shown to maintain specific activity in epilepsy where traditional promoters changed activity, and used to drive functional transgenes in astrocytes including Cre recombinase and acetylcholine-responsive sensor iAChSnFR. The astrocyte-specific iAChSnFR revealed a clear reward-dependent acetylcholine response in astrocytes of the nucleus accumbens during reinforcement learning. Together, this collection of glial enhancer-AAVs will enable characterization of astrocyte and oligodendrocyte populations and their roles across species, disease states, and behavioral epochs.
Experimental access to cell types within the mammalian spinal cord is severely limited by the availability of genetic tools. To enable access to spinal motor neurons (SMNs) and SMN subtypes, we generated single-cell multiome datasets from mouse and macaque spinal cords and discovered putative enhancers for each neuronal population. We cloned these enhancers into adeno-associated viral vectors driving a reporter fluorophore and functionally screened them in the mouse. We extensively characterized the most promising candidate enhancers in rat and macaque and developed an optimized pan-SMN enhancer virus. Additionally, we generated derivative viruses expressing iCre297T recombinase or ChR2-EYFP for labeling and functional studies, and we created a single vector with combined enhancer elements to achieve simultaneous labeling of layer 5 extratelencephalic projecting neurons and SMNs. This unprecedented SMN toolkit will enable future investigations of cell type function across species and potential therapeutic interventions for human neurodegenerative diseases.
Cis -acting regulatory enhancer elements are valuable tools for gaining cell type-specific genetic access. Leveraging large chromatin accessibility atlases, putative enhancer sequences can be identified and deployed in adeno-associated virus (AAV) delivery platforms. However, a significant bottleneck in enhancer AAV discovery is charting their detailed expression patterns in vivo , a process that currently requires gold-standard one-by-one testing. Here we present a barcoded multiplex strategy for screening enhancer AAVs at cell type resolution using single cell RNA sequencing and taxonomy mapping. We executed a proof-of-concept study using a small pool of validated enhancer AAVs expressing in a variety of neuronal and non-neuronal cell types across the mouse brain. Unexpectedly, we encountered substantial technical and biological noise including chimeric packaging products, necessitating development of novel techniques to accurately deconvolve enhancer expression patterns. These results underscore the need for improved methods to mitigate noise and highlight the complexity of enhancer AAV biology in vivo .
Biological ageing can be defined as a gradual loss of homeostasis across various aspects of molecular and cellular function(1,2). Mammalian brains consist of thousands of cell types(3), which may be differentially susceptible or resilient to ageing. Here we present a comprehensive single-cell RNA sequencing dataset containing roughly 1.2million high-quality single-cell transcriptomes of brain cells from young adult and aged mice of both sexes, from regions spanning the forebrain, midbrain and hindbrain. High-resolution clustering of all cells results in 847 cell clusters and reveals at least 14 age-biased clusters that are mostly glial types. At the broader cell subclass and supertype levels, we find age-associated gene expression signatures and provide a list of 2,449 unique differentially expressed genes (age-DE genes) for many neuronal and non-neuronal cell types. Whereas most age-DE genes are unique to specific cell types, we observe common signatures with ageing across cell types, including a decrease in expression of genes related to neuronal structure and function in many neuron types, major astrocyte types and mature oligodendrocytes, and an increase in expression of genes related to immune function, antigen presentation, inflammation, and cell motility in immune cell types and some vascular cell types. Finally, we observe that some of the cell types that demonstrate the greatest sensitivity to ageing are concentrated around the third ventricle in the hypothalamus, including tanycytes, ependymal cells, and certain neuron types in the arcuate nucleus, dorsomedial nucleus and paraventricular nucleus that express genes canonically related to energy homeostasis. Many of these types demonstrate both a decrease in neuronal function and an increase in immune response. These findings suggest that the third ventricle in the hypothalamus may be a hub for ageing in the mouse brain. Overall, this study systematically delineates a dynamic landscape of cell-type-specific transcriptomic changes in the brain associated with normal ageing that will serve as a foundation for the investigation of functional changes in ageing and the interaction of ageing and disease.
The telencephalon of the mammalian brain comprises multiple regions and circuit pathways that play adaptive and integrative roles in a variety of brain functions. There is a wide array of GABAergic neurons in the telencephalon; they play a multitude of circuit functions, and dysfunction of these neurons has been implicated in diverse brain disorders. In this study, we conducted a systematic and in-depth analysis of the transcriptomic and spatial organization of GABAergic neuronal types in all regions of the mouse telencephalon and their developmental origins. This was accomplished by utilizing 611,423 single-cell transcriptomes from the comprehensive and high-resolution transcriptomic and spatial cell type atlas for the adult whole mouse brain we have generated, supplemented with an additional single-cell RNA-sequencing dataset containing 99,438 high-quality single-cell transcriptomes collected from the pre- and postnatal developing mouse brain. We present a hierarchically organized adult telencephalic GABAergic neuronal cell type taxonomy of 7 classes, 52 subclasses, 284 supertypes, and 1,051 clusters, as well as a corresponding developmental taxonomy of 450 clusters across different ages. Detailed charting efforts reveal extraordinary complexity where relationships among cell types reflect both spatial locations and developmental origins. Transcriptomically and developmentally related cell types can often be found in distant and diverse brain regions indicating that long-distance migration and dispersion is a common characteristic of nearly all classes of telencephalic GABAergic neurons. Additionally, we find various spatial dimensions of both discrete and continuous variations among related cell types that are correlated with gene expression gradients. Lastly, we find that cortical, striatal and some pallidal GABAergic neurons undergo extensive postnatal diversification, whereas septal and most pallidal GABAergic neuronal types emerge simultaneously during the embryonic stage with limited postnatal diversification. Overall, the telencephalic GABAergic cell type taxonomy can serve as a foundational reference for molecular, structural and functional studies of cell types and circuits by the entire community.
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.
The distinctive physiology of striatal medium spiny neurons (MSNs) underlies their ability to integrate sensory and motor input. In rodents, MSNs have a hyperpolarized resting potential and low input resistance. When activated, they have a delayed onset of spiking and regular spike rate. Here, we show that in the macaque putamen, latency to spike is reduced and spike rate adaptation is increased relative to mouse. We use whole-cell brain slice recordings and recover single-cell gene expression using Patch-seq to distinguish macaque MSN cell types. Species differences in the expression of ion channel genes including the calcium-activated chloride channel, ANO2, and an auxiliary subunit of the A-type potassium channel, DPP10, are correlated with species differences in spike rate adaptation and latency to the first spike, respectively. These surprising divergences in physiology better define the strengths and limitations of mouse models for understanding neuronal and circuit function in the primate basal ganglia.
The mammalian cortex is composed of a highly diverse set of cell types and develops through a series of temporally regulated events that build out the cell type and circuit foundation for cortical function. The mechanisms underlying the development of different cell types remain elusive. Single-cell transcriptomics provides the capacity to systematically study cell types across the entire temporal range of cortical development. Here, we present a comprehensive and high-resolution transcriptomic and epigenomic cell type atlas of the developing mouse visual cortex. The atlas was built from a single-cell RNA-sequencing dataset of 568,674 high-quality single-cell transcriptomes and a single-nucleus Multiome dataset of 194,545 high-quality nuclei providing both transcriptomic and chromatin accessibility profiles, densely sampled throughout the embryonic and postnatal developmental stages from E11.5 to P56. We computationally reconstructed a transcriptomic developmental trajectory map of all excitatory, inhibitory, and non-neuronal cell types in the visual cortex, identifying branching points marking the emergence of new cell types at specific developmental ages and defining molecular signatures of cellular diversification. In addition to neurogenesis, gliogenesis and early postmitotic maturation in the embryonic stage which gives rise to all the cell classes and nearly all subclasses, we find that increasingly refined cell types emerge throughout the postnatal differentiation process, including the late emergence of many cell types during the eye-opening stage (P11-P14) and the onset of critical period (P21), suggesting continuous cell type diversification at different stages of cortical development. Throughout development, we find cooperative dynamic changes in gene expression and chromatin accessibility in specific cell types, identifying both chromatin peaks potentially regulating the expression of specific genes and transcription factors potentially regulating specific peaks. Furthermore, a single gene can be regulated by multiple peaks associated with different cell types and/or different developmental stages. Collectively, our study provides the most detailed dynamic molecular map directly associated with individual cell types and specific developmental events that reveals the molecular logic underlying the continuous refinement of cell type identities in the developing visual cortex.
Large-scale analysis of single-cell gene expression has revealed transcriptomically defined cell subclasses present throughout the primate neocortex with gene expression profiles that differ depending upon neocortical region. Here, we test whether the interareal differences in gene expression translate to regional specializations in the physiology and morphology of infragranular glutamatergic neurons by performing Patch-seq experiments in brain slices from the temporal cortex (TCx) and motor cortex (MCx) of the macaque. We confirm that transcriptomically defined extratelencephalically projecting neurons of layer 5 (L5 ET neurons) include retrogradely labeled corticospinal neurons in the MCx and find multiple physiological properties and ion channel genes that distinguish L5 ET from non-ET neurons in both areas. Additionally, while infragranular ET and non-ET neurons retain distinct neuronal properties across multiple regions, there are regional morpho-electric and gene expression specializations in the L5 ET subclass, providing mechanistic insights into the specialized functional architecture of the primate neocortex.
Alzheimer's disease (AD) is the most common cause of dementia in older adults. Neuropathological and imaging studies have demonstrated a progressive and stereotyped accumulation of protein aggregates, but the underlying molecular and cellular mechanisms driving AD progression and vulnerable cell populations affected by disease remain coarsely understood. The current study harnesses single cell and spatial genomics tools and knowledge from the BRAIN Initiative Cell Census Network to understand the impact of disease progression on middle temporal gyrus cell types. We used image-based quantitative neuropathology to place 84 donors spanning the spectrum of AD pathology along a continuous disease pseudoprogression score and multiomic technologies to profile single nuclei from each donor, mapping their transcriptomes, epigenomes, and spatial coordinates to a common cell type reference with unprecedented resolution. Temporal analysis of cell-type proportions indicated an early reduction of Somatostatin-expressing neuronal subtypes and a late decrease of supragranular intratelencephalic-projecting excitatory and Parvalbumin-expressing neurons, with increases in disease-associated microglial and astrocytic states. We found complex gene expression differences, ranging from global to cell type-specific effects. These effects showed different temporal patterns indicating diverse cellular perturbations as a function of disease progression. A subset of donors showed a particularly severe cellular and molecular phenotype, which correlated with steeper cognitive decline. We have created a freely available public resource to explore these data and to accelerate progress in AD research at SEA-AD.org.
Single-cell analyses parse the brain's billions of neurons into thousands of 'cell-type' clusters residing in different brain structures1. Many cell types mediate their functions through targeted long-distance projections allowing interactions between specific cell types. Here we used epi-retro-seq2 to link single-cell epigenomes and cell types to long-distance projections for 33,034 neurons dissected from 32 different regions projecting to 24 different targets (225 source-to-target combinations) across the whole mouse brain. We highlight uses of these data for interrogating principles relating projection types to transcriptomics and epigenomics, and for addressing hypotheses about cell types and connections related to genetics. We provide an overall synthesis with 926 statistical comparisons of discriminability of neurons projecting to each target for every source. We integrate this dataset into the larger BRAIN Initiative Cell Census Network atlas, composed of millions of neurons, to link projection cell types to consensus clusters. Integration with spatial transcriptomics further assigns projection-enriched clusters to smaller source regions than the original dissections. We exemplify this by presenting in-depth analyses of projection neurons from the hypothalamus, thalamus, hindbrain, amygdala and midbrain to provide insights into properties of those cell types, including differentially expressed genes, their associated cis-regulatory elements and transcription-factor-binding motifs, and neurotransmitter use.
Volume 66, no. 1, e0164321, 2022, https://doi.org/10.1128/AAC.01643-21. Page 12, “Data availability” paragraph, line 2: The reported GenBank accession numbers (MZ568467 to MZ568625) are incorrect. The GenBank accession numbers should be MZ568548 to MZ568625.
Alzheimer’s disease (AD) is the most common form of dementia. The progression of AD throughout the brain follows a stereotypical pattern that can be described by the quantification of histopathological changes in several brain regions. As part of The Seattle Alzheimer’s Disease Cell Atlas (SEA-AD, https://sea-ad.org ), we previously quantified pathological proteins in the middle temporal gyrus (MTG) and used them to describe disease progression, vulnerable cell types, and the MTG molecular changes related to disease In this work, we expanded our results and undertook a comparative approach. We measured a battery of neuropathological proteins (staining NeuN, AT8, IBA1, GFAP, 6e10, pTDP43, and a-Syn immunoreactivity, and analyzing them with machine learning) in the middle temporal gyrus, the middle frontal gyrus, and medial entorhinal cortex in 84 donors. In addition, we undertook a multimodal single-cell analysis of such regions, profiling cells using single-cell RNA-seq, ATAC-seq, Multiome, and spatial transcriptomic. Borrowing from our hierarchically resolved cell types in MTG, based on the BRAIN initiative cell type reference, we used machine learning algorithms to create multimodal cell type maps across regions. We used Bayesian statistics to identify differentially accessible regulatory elements, link them to the affected genes and the regulatory transcriptional machinery binding in such elements. We identified regional vulnerable neuronal populations and a core set of transcriptionally similar cell types susceptible to disease across the cortex. We created a cortex-wide disease progression map combining neuropathological data from our three profiled regions with latent Bayesian statistical algorithms. In each area, we developed an AD pseudo-progression timescale that describes the burden of disease in the area and, next, combined all areas to create a cortex-wide holistic pseudo-progression. We used the mapped cell types and our developed scale to synchronize changes occurring in each patient, describing abundances, gene expression, and chromatin accessibility modifications. The progression scale developed in this study serves as a reference to align changes across brain regions and, in the future, to reconcile differences across cohorts.
The mammalian brain consists of millions to billions of cells that are organized into many cell types with specific spatial distribution patterns and structural and functional properties1-3. Here we report a comprehensive and high-resolution transcriptomic and spatial cell-type atlas for the whole adult mouse brain. The cell-type atlas was created by combining a single-cell RNA-sequencing (scRNA-seq) dataset of around 7 million cells profiled (approximately 4.0 million cells passing quality control), and a spatial transcriptomic dataset of approximately 4.3 million cells using multiplexed error-robust fluorescence in situ hybridization (MERFISH). The atlas is hierarchically organized into 4 nested levels of classification: 34 classes, 338 subclasses, 1,201 supertypes and 5,322 clusters. We present an online platform, Allen Brain Cell Atlas, to visualize the mouse whole-brain cell-type atlas along with the single-cell RNA-sequencing and MERFISH datasets. We systematically analysed the neuronal and non-neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell-type organization in different brain regions-in particular, a dichotomy between the dorsal and ventral parts of the brain. The dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. Our study also uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types. Finally, we found that transcription factors are major determinants of cell-type classification and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole mouse brain transcriptomic and spatial cell-type atlas establishes a benchmark reference atlas and a foundational resource for integrative investigations of cellular and circuit function, development and evolution of the mammalian brain.