Molecularly defined and spatially resolved cell atlas of the whole mouse brain

A platform for multimodal in vivo pooled genetic screens reveals regulators of liver function

Organ function requires coordinated activities of thousands of genes in distinct, spatially organized cell types. Understanding the basis of emergent tissue function requires approaches to dissect the genetic control of diverse cellular and tissue phenotypes in vivo. Here, we develop paired imaging and sequencing methods to construct large-scale, multi-modal genotype-phenotypes maps in tissue with pooled genetic perturbations. Using imaging, we identify genetic perturbations in individual cells while simultaneously measuring their gene expression and subcellular morphology. Using single-cell sequencing, we measure transcriptomic responses to the same genetic perturbations. We apply this approach to study hundreds of genetic perturbations in the mouse liver. Our study reveals regulators of hepatocyte zonation and liver unfolded protein response, as well as distinct pathways that cause hepatocyte steatosis. Our approach enables new ways of interrogating the genetic basis of complex cellular and organismal physiology and provides crucial training data for emerging machine-learning models of cellular function.

Identification and characterization of cell niches in tissue from spatial omics data at single-cell resolution

Deciphering the features, structure, and functions of the cell niche in tissues remains a major challenge. Here, we present scNiche, a computational framework to identify and characterize cell niches from spatial omics data at single-cell resolution. We benchmark scNiche with both simulated and biological datasets, and demonstrate that scNiche can effectively and robustly identify cell niches while outperforming other existing methods. In spatial proteomics data from human triple-negative breast cancer, scNiche reveals the influence of the microenvironment on cellular phenotypes, and further dissects patient-specific niches with distinct cellular compositions or phenotypic characteristics. By analyzing mouse liver spatial transcriptomics data across normal and early-onset liver failure donors, scNiche uncovers disease-specific liver injury niches, and further delineates the niche remodeling from normal liver to liver failure. Overall, scNiche enables decoding the cellular microenvironment in tissues from single-cell spatial omics data.

Rarely categorical, always high-dimensional: how the neural code changes along the cortical hierarchy

A long-standing debate in neuroscience concerns whether individual neurons are organized into functionally distinct populations that encode information differently ("categorical" representations [1-3]) and the implications for neural computation. Here, we systematically analyzed how cortical neurons encode cognitive, sensory, and movement variables across 43 cortical regions during a complex task (14,000+ units from the International Brain Laboratory public Brain-wide Map data set [4]) and studied how these properties change across the sensory-cognitive cortical hierarchy [5]. We found that the structure of the neural code was scale-dependent: on a whole-cortex scale, neural selectivity was categorical and organized across regions in a way that reflected their anatomical connectivity. However, within individual regions, categorical representations were rare and limited to primary sensory areas. Remarkably, the degree of categorical clustering of neural selectivity was inversely correlated to the dimensionality of neural representations, suggesting a link between single-neuron selectivity and computational properties of population codes that we explained in a mathematical model. Finally, we found that the fraction of linearly separable combinations of experimental conditions ("Shattering Dimensionality" [6]) was near maximal across all areas, indicating a robust and uniform ability for flexible information encoding throughout the cortex. In conclusion, our results provide systematic evidence for a non-categorical, high-dimensional neural code in all but the lower levels of the cortical hierarchy.

Atypical hippocampal excitatory neurons express and govern object memory

Classically, pyramidal cells of the hippocampus are viewed as flexibly representing spatial and non-spatial information. Recent work has illustrated distinct types of hippocampal excitatory neurons, suggesting that hippocampal representations and functions may be constrained and interpreted by these underlying cell-type identities. In mice, here we reveal a non-pyramidal excitatory neuron type - the "ovoid" neuron - that is spatially adjacent to subiculum pyramidal cells but differs in gene expression, electrophysiology, morphology, and connectivity. Functionally, novel object encounters drive sustained ovoid neuron activity, whereas familiar objects fail to drive activity even months after single-trial learning. Silencing ovoid neurons prevents non-spatial object learning but leaves spatial learning intact, and activating ovoid neurons toggles novel-object seeking to familiar-object seeking. Such function is doubly dissociable from pyramidal neurons, wherein manipulation of pyramidal cells affects spatial assays but not non-spatial learning. Ovoid neurons of the subiculum thus illustrate selective cell-type-specific control of non-spatial memory and behavioral preference.

A Complete Spatial Map of Mouse Retinal Ganglion Cells Reveals Density and Gene Expression Specializations

Retinal ganglion cells (RGCs) transmit visual information from the eye to the brain. In mice, several RGC subtypes show nonuniform spatial distributions, potentially mediating specific visual functions. However, the full extent of RGC specialization remains unknown. Here, we used en-face cryosectioning, spatial transcriptomics, and machine learning to map the spatial distribution of all RGC subtypes identified in previous single-cell studies. While two-thirds of RGC subtypes were evenly distributed, others showed strong biases toward ventral or dorso-temporal regions associated with sky vision and the area retinae temporalis (ART), the predicted homolog of the area centralis. Additionally, we observed unexpected spatial variation in gene expression within several subtypes along the dorso-ventral axis or within vs. outside the ART, independent of RGC density profiles. Finally, we found limited correlations between the gene profiles of the ART and the primate macula, suggesting divergent specialization between the mouse and primate central vision.

Annotation Comparison Explorer (ACE): connecting brain cell types across studies of health and Alzheimer's Disease

Single-cell multiomic technologies have allowed unprecedented access to gene profiles of individual cells across species and organ systems, including >1000 papers focused on brain cell types alone. The Allen Institute has created foundational atlases characterizing mammalian brain cell types in the adult mouse brain and the neocortex of aged humans with and without Alzheimer's disease (AD). With so many public cell type classifications (or 'taxonomies') available and many groups choosing to define their own, linking cell types and associated knowledge between studies remains a major challenge. Here, we introduce Annotation Comparison Explorer (ACE), a web application for comparing cell type assignments and other cell-based annotations (e.g., donor demographics, anatomic locations, batch variables, and quality control metrics). ACE allows filtering of cells and includes an interactive set of tools for comparing two or more taxonomy annotations alongside collected knowledge (e.g., increased abundance in disease conditions, cell type aliases, or other information about a specific cell type). We present three primary use cases for ACE. First, we demonstrate how a user can assign cell type labels from the Seattle Alzheimer's Disease Brain Cell Atlas (SEA-AD) taxonomy to cells from their own study and compare these cell type mappings to existing cell type assignments and cell metadata. Second, we extend this approach to ten published human AD studies which we previously reprocessed through a common data analysis pipeline. This allowed us to compare brain taxonomies across otherwise incomparable studies and identify congruent cell type abundance changes in AD, including a decrease in abundance of subsets of somatostatin interneurons. Finally, ACE includes translation tables between different mouse and human brain cell type taxonomies publicly accessible on Allen Brain Map, from initial studies in individual neocortical areas to more recent studies spanning the whole brain. ACE can be freely and publicly accessed as a web application (https://sea-ad.shinyapps.io/ACEapp/) and on GitHub (github.com/AllenInstitute/ACE).

Spatial transcriptomic clocks reveal cell proximity effects in brain ageing

Old age is associated with a decline in cognitive function and an increase in neurodegenerative disease risk1. Brain ageing is complex and is accompanied by many cellular changes2. Furthermore, the influence that aged cells have on neighbouring cells and how this contributes to tissue decline is unknown. More generally, the tools to systematically address this question in ageing tissues have not yet been developed. Here we generate a spatially resolved single-cell transcriptomics brain atlas of 4.2 million cells from 20 distinct ages across the adult lifespan and across two rejuvenating interventions-exercise and partial reprogramming. We build spatial ageing clocks, machine learning models trained on this spatial transcriptomics atlas, to identify spatial and cell-type-specific transcriptomic fingerprints of ageing, rejuvenation and disease, including for rare cell types. Using spatial ageing clocks and deep learning, we find that T cells, which increasingly infiltrate the brain with age, have a marked pro-ageing proximity effect on neighbouring cells. Surprisingly, neural stem cells have a strong pro-rejuvenating proximity effect on neighbouring cells. We also identify potential mediators of the pro-ageing effect of T cells and the pro-rejuvenating effect of neural stem cells on their neighbours. These results suggest that rare cell types can have a potent influence on their neighbours and could be targeted to counter tissue ageing. Spatial ageing clocks represent a useful tool for studying cell-cell interactions in spatial contexts and should allow scalable assessment of the efficacy of interventions for ageing and disease.

Single-cell spatial transcriptomics reveals distinct patterns of dysregulation in non-neuronal and neuronal cells induced by the Trem2R47H Alzheimer's risk gene mutation

The R47H missense mutation of the TREM2 gene is a known risk factor for development of Alzheimer's Disease. In this study, we analyze the impact of the Trem2R47H mutation on specific cell types in multiple cortical and subcortical brain regions in the context of wild-type and 5xFAD mouse background. We profile 19 mouse brain sections consisting of wild-type, Trem2R47H, 5xFAD and Trem2R47H; 5xFAD genotypes using MERFISH spatial transcriptomics, a technique that enables subcellular profiling of spatial gene expression. Spatial transcriptomics and neuropathology data are analyzed using our custom pipeline to identify plaque and Trem2R47H-induced transcriptomic dysregulation. We initially analyze cell type-specific transcriptomic alterations induced by plaque proximity. Next, we analyze spatial distributions of disease associated microglia and astrocytes, and how they vary between 5xFAD and Trem2R47H; 5xFAD mouse models. Finally, we analyze the impact of the Trem2R47H mutation on neuronal transcriptomes. The Trem2R47H mutation induces consistent upregulation of Bdnf and Ntrk2 across many cortical excitatory neuron types, independent of amyloid pathology. Spatial investigation of genotype enriched subclusters identified spatially localized neuronal subpopulations reduced in 5xFAD and Trem2R47H; 5xFAD mice. Overall, our MERFISH spatial transcriptomics analysis identifies glial and neuronal transcriptomic alterations induced independently by 5xFAD and Trem2R47H mutations, impacting inflammatory responses in microglia and astrocytes, and activity and BDNF signaling in neurons.

Brain-wide cell-type-specific transcriptomic signatures of healthy ageing in mice

Biological ageing can be defined as a gradual loss of homeostasis across various aspects of molecular and cellular function1,2. Mammalian brains consist of thousands of cell types3, which may be differentially susceptible or resilient to ageing. Here we present a comprehensive single-cell RNA sequencing dataset containing roughly 1.2 million 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.

Categorization of 34 computational methods to detect spatially variable genes from spatially resolved transcriptomics data

In the analysis of spatially resolved transcriptomics data, detecting spatially variable genes (SVGs) is crucial. Numerous computational methods exist, but varying SVG definitions and methodologies lead to incomparable results. We review 34 state-of-the-art methods, classifying SVGs into three categories: overall, cell-type-specific, and spatial-domain-marker SVGs. Our review explains the intuitions underlying these methods, summarizes their applications, and categorizes the hypothesis tests they use in the trade-off between generality and specificity for SVG detection. We discuss challenges in SVG detection and propose future directions for improvement. Our review offers insights for method developers and users, advocating for category-specific benchmarking.

Characterizing and targeting glioblastoma neuron-tumor networks with retrograde tracing

Glioblastomas are invasive brain tumors with high therapeutic resistance. Neuron-to-glioma synapses have been shown to promote glioblastoma progression. However, a characterization of tumor-connected neurons has been hampered by a lack of technologies. Here, we adapted retrograde tracing using rabies viruses to investigate and manipulate neuron-tumor networks. Glioblastoma rapidly integrated into neural circuits across the brain, engaging in widespread functional communication, with cholinergic neurons driving glioblastoma invasion. We uncovered patient-specific and tumor-cell-state-dependent differences in synaptogenic gene expression associated with neuron-tumor connectivity and subsequent invasiveness. Importantly, radiotherapy enhanced neuron-tumor connectivity by increased neuronal activity. In turn, simultaneous neuronal activity inhibition and radiotherapy showed increased therapeutic effects, indicative of a role for neuron-to-glioma synapses in contributing to therapeutic resistance. Lastly, rabies-mediated genetic ablation of tumor-connected neurons halted glioblastoma progression, offering a viral strategy to tackle glioblastoma. Together, this study provides a framework to comprehensively characterize neuron-tumor networks and target glioblastoma.

Technical and biological sources of noise confound multiplexed enhancer AAV screening

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 .

Calcineurin: An essential regulator of sleep revealed by biochemical, chemical biological, and genetic approaches

Research into mechanisms underlying sleep traditionally relies on electrophysiology and genetics. Because sleep can only be measured on whole animals by behavioral observations and physical means, no sleep research was initiated by biochemical and chemical biological approaches. We used phosphorylation sites of kinases important for sleep as targets for biochemical and chemical biological approaches. Sleep was increased in mice carrying a threonine-to-alanine substitution at residue T469 of salt-inducible kinase 3 (SIK3). Our biochemical purification and photo-crosslinking revealed calcineurin (CaN) dephosphorylation, both in vitro and in vivo, of SIK3 at T469 and S551, but not T221. Knocking down CaN regulatory subunit reduced daily sleep by more than 5 h, exceeding all known mouse mutants. Our work uncovered a critical physiological role for CaN in sleep and pioneered biochemical purification and chemical biology as effective approaches to study sleep.

InTiCAR: Network-based identification of significant inter-tissue communicators for autoimmune diseases

Inter-tissue communicators (ITCs) are intricate and essential aspects of our body, as they are the keepers of homeostatic equilibrium. It is no surprise that the dysregulation of the exchange between tissues are at the core of various disorders. Among such conditions, autoimmune diseases (AIDs) refer to a collection of pathological conditions where the miscommunication drives the immune system to mistakenly attack one's own body. Due to their myriad and diverse pathophysiologies, AIDs cannot be easily diagnosed or treated, and continuous efforts are required to seek for potential diagnostic markers or therapeutic targets. The identification of ITCs with significant involvement in the disease states is therefore crucial. Here, we present InTiCAR, Inter-Tissue Communicators for Autoimmune diseases by Random walk with restart, which is a network exploration-based analysis method that suggests disease-specific ITCs based on prior knowledge of disease genes, without the need for the external expression data. We first show that distinct ITC profile s can be acquired for various diseases by InTiCAR. We further illustrate that, for autoimmune diseases (AIDs) specifically, the disease-specific ITCs outperform disease genes in diagnosing patients using the UK Biobank plasma proteome dataset. Also, through CMap LINCS dataset, we find that high perturbation on the AIDs genes can be observed by the disease-specific ITCs. Our results provide and highlight unique perspectives on biological network analysis by focusing on the entities of extracellular communications.

Astrocytes in aging

The mammalian nervous system is impacted by aging. Aging alters brain architecture, is associated with molecular damage, and can manifest with cognitive and motor deficits that diminish the quality of life. Astrocytes are glial cells of the CNS that regulate the development, function, and repair of neural circuits during development and adulthood; however, their functions in aging are less understood. Astrocytes change their transcriptome during aging, with astrocytes in areas such as the cerebellum, the hypothalamus, and white matter-rich regions being the most affected. While numerous studies describe astrocyte transcriptional changes in aging, many questions still remain. For example, how is astrocyte function altered by transcriptional changes that occur during aging? What are the mechanisms promoting astrocyte aged states? How do aged astrocytes impact brain function? This review discusses features of aged astrocytes and their potential triggers and proposes ways in which they may impact brain function and health span.

From Physiology to Pathology of Astrocytes: Highlighting Their Potential as Therapeutic Targets for CNS Injury

In the mammalian central nervous system (CNS), astrocytes are the ubiquitous glial cells that have complex morphological and molecular characteristics. These fascinating cells play essential neurosupportive and homeostatic roles in the healthy CNS and undergo morphological, molecular, and functional changes to adopt so-called 'reactive' states in response to CNS injury or disease. In recent years, interest in astrocyte research has increased dramatically and some new biological features and roles of astrocytes in physiological and pathological conditions have been discovered thanks to technological advances. Here, we will review and discuss the well-established and emerging astroglial biology and functions, with emphasis on their potential as therapeutic targets for CNS injury, including traumatic and ischemic injury. This review article will highlight the importance of astrocytes in the neuropathological process and repair of CNS injury.

The neurobiology of parenting and infant-evoked aggression

Parenting behavior comprises a variety of adult-infant and adult-adult interactions across multiple timescales. The state transition from nonparent to parent requires an extensive reorganization of individual priorities and physiology and is facilitated by combinatorial hormone action on specific cell types that are integrated throughout interconnected and brainwide neuronal circuits. In this review, we take a comprehensive approach to integrate historical and current literature on each of these topics across multiple species, with a focus on rodents. New and emerging molecular, circuit-based, and computational technologies have recently been used to address outstanding gaps in our current framework of knowledge on infant-directed behavior. This work is raising fundamental questions about the interplay between instinctive and learned components of parenting and the mutual regulation of affiliative versus agonistic infant-directed behaviors in health and disease. Whenever possible, we point to how these technologies have helped gain novel insights and opened new avenues of research into the neurobiology of parenting. We hope this review will serve as an introduction for those new to the field, a comprehensive resource for those already studying parenting, and a guidepost for designing future studies.

Human-unique brain cell clusters are associated with learning disorders and human episodic memory activity

The advanced evolution of the human cerebral cortex forms the basis for our high-level cognitive functions. Through a comparative analysis of single-nucleus transcriptome data from the human neocortex and that of chimpanzees, macaques, and marmosets, we discovered 20 subgroups of cell types unique to the human brain, which include 11 types of excitatory neurons. Many of these human-unique cell clusters exhibit significant overexpression of genes regulated by human-specific enhancers. Notably, these specific cell clusters also express genes associated with disease risk, particularly those related to brain dysfunctions like learning disorders. Furthermore, genes linked to cortical thickness and human episodic memory encoding activities show heightened expression within these cell subgroups. These findings underscore the critical role of human brain-unique cell clusters in the evolution of human brain functions.

Considerations for building and using integrated single-cell atlases

The rapid adoption of single-cell technologies has created an opportunity to build single-cell 'atlases' integrating diverse datasets across many laboratories. Such atlases can serve as a reference for analyzing and interpreting current and future data. However, it has become apparent that atlasing approaches differ, and the impact of these differences are often unclear. Here we review the current atlasing literature and present considerations for building and using atlases. Importantly, we find that no one-size-fits-all protocol for atlas building exists, but rather we discuss context-specific considerations and workflows, including atlas conceptualization, data collection, curation and integration, atlas evaluation and atlas sharing. We further highlight the benefits of integrated atlases for analyses of new datasets and deriving biological insights beyond what is possible from individual datasets. Our overview of current practices and associated recommendations will improve the quality of atlases to come, facilitating the shift to a unified, reference-based understanding of single-cell biology.

Profiling cell identity and tissue architecture with single-cell and spatial transcriptomics

Single-cell transcriptomics has broadened our understanding of cellular diversity and gene expression dynamics in healthy and diseased tissues. Recently, spatial transcriptomics has emerged as a tool to contextualize single cells in multicellular neighbourhoods and to identify spatially recurrent phenotypes, or ecotypes. These technologies have generated vast datasets with targeted-transcriptome and whole-transcriptome profiles of hundreds to millions of cells. Such data have provided new insights into developmental hierarchies, cellular plasticity and diverse tissue microenvironments, and spurred a burst of innovation in computational methods for single-cell analysis. In this Review, we discuss recent advancements, ongoing challenges and prospects in identifying and characterizing cell states and multicellular neighbourhoods. We discuss recent progress in sample processing, data integration, identification of subtle cell states, trajectory modelling, deconvolution and spatial analysis. Furthermore, we discuss the increasing application of deep learning, including foundation models, in analysing single-cell and spatial transcriptomics data. Finally, we discuss recent applications of these tools in the fields of stem cell biology, immunology, and tumour biology, and the future of single-cell and spatial transcriptomics in biological research and its translation to the clinic.

Degenerate neuronal and circuit mechanisms important for generating rhythmic motor patterns

In 1996, we published a review article (Marder E, Calabrese RL. Physiol Rev 76: 687-717, 1996) describing the state of knowledge about the structure and function of the central pattern-generating circuits important for producing rhythmic behaviors. Although many of the core questions persist, much has changed since 1996. Here, we focus on newer studies that reveal ambiguities that complicate understanding circuit dynamics, despite the enormous technical advances of the recent past. In particular, we highlight recent studies of animal-to-animal variability and our understanding that circuit rhythmicity may be supported by multiple state-dependent mechanisms within the same animal and that robustness and resilience in the face of perturbation may depend critically on the presence of modulators and degenerate circuit mechanisms. Additionally, we highlight the use of computational models to ask whether there are generalizable principles about circuit motifs that can be found across rhythmic motor systems in different animal species.

PCA-based spatial domain identification with state-of-the-art performance

MOTIVATION

The identification of biologically meaningful domains is a central step in the analysis of spatial transcriptomic data.

RESULTS

Following Occam's razor, we show that a simple PCA-based algorithm for unsupervised spatial domain identification rivals the performance of ten competing state-of-the-art methods across six single-cell spatial transcriptomic datasets. Our reductionist approach, NichePCA, provides researchers with intuitive domain interpretation and excels in execution speed, robustness, and scalability.

AVAILABILITY AND IMPLEMENTATION

The code is available at https://github.com/imsb-uke/nichepca.

Spatiotemporal modeling of molecular holograms

Quantifying spatiotemporal dynamics during embryogenesis is crucial for understanding congenital diseases. We developed Spateo (https://github.com/aristoteleo/spateo-release), a 3D spatiotemporal modeling framework, and applied it to a 3D mouse embryogenesis atlas at E9.5 and E11.5, capturing eight million cells. Spateo enables scalable, partial, non-rigid alignment, multi-slice refinement, and mesh correction to create molecular holograms of whole embryos. It introduces digitization methods to uncover multi-level biology from subcellular to whole organ, identifying expression gradients along orthogonal axes of emergent 3D structures, e.g., secondary organizers such as midbrain-hindbrain boundary (MHB). Spateo further jointly models intercellular and intracellular interaction to dissect signaling landscapes in 3D structures, including the zona limitans intrathalamica (ZLI). Lastly, Spateo introduces "morphometric vector fields" of cell migration and integrates spatial differential geometry to unveil molecular programs underlying asymmetrical murine heart organogenesis and others, bridging macroscopic changes with molecular dynamics. Thus, Spateo enables the study of organ ecology at a molecular level in 3D space over time.

Molecular and spatial transcriptomic classification of midbrain dopamine neurons and their alterations in a LRRK2G2019S model of Parkinson's disease

Several studies have revealed that midbrain dopamine (DA) neurons, even within a single neuroanatomical area, display heterogeneous properties. In parallel, studies using single cell profiling techniques have begun to cluster DA neurons into subtypes based on their molecular signatures. Recent work has shown that molecularly defined DA subtypes within the substantia nigra (SNc) display distinctive anatomic and functional properties, and differential vulnerability in Parkinson's disease (PD). Based on these provocative results, a granular understanding of these putative subtypes and their alterations in PD models, is imperative. We developed an optimized pipeline for single-nuclear RNA sequencing (snRNA-seq) and generated a high-resolution hierarchically organized map revealing 20 molecularly distinct DA neuron subtypes belonging to three main families. We integrated this data with spatial MERFISH technology to map, with high definition, the location of these subtypes in the mouse midbrain, revealing heterogeneity even within neuroanatomical sub-structures. Finally, we demonstrate that in the preclinical LRRK2G2019S knock-in mouse model of PD, subtype organization and proportions are preserved. Transcriptional alterations occur in many subtypes including those localized to the ventral tier SNc, where differential expression is observed in synaptic pathways, which might account for previously described DA release deficits in this model. Our work provides an advancement of current taxonomic schemes of the mouse midbrain DA neuron subtypes, a high-resolution view of their spatial locations, and their alterations in a prodromal mouse model of PD.

BDNF augmentation reverses cranial radiation therapy-induced cognitive decline and neurodegenerative consequences

Cranial radiation therapy (RT) for brain cancers is often associated with the development of radiation-induced cognitive dysfunction (RICD). RICD significantly impacts the quality of life for cancer survivors, highlighting an unmet medical need. Previous human studies revealed a marked reduction in plasma brain-derived neurotrophic factor (BDNF) post-chronic chemotherapy, linking this decline to a substantial cognitive dysfunction among cancer survivors. Moreover, riluzole (RZ)-mediated increased BDNF in vivo in the chemotherapy-exposed mice reversed cognitive decline. RZ is an FDA-approved medication for ALS known to increase BDNF in vivo. In an effort to mitigate the detrimental effects of RT-induced BDNF decline in RICD, we tested the efficacy of RZ in a cranially irradiated (9 Gy) adult mouse model. Notably, RT-exposed mice exhibited significantly reduced hippocampal BDNF, accompanied by increased neuroinflammation, loss of neuronal plasticity-related immediate early gene product, cFos, and synaptic density. Spatial transcriptomic profiling comparing the RT + Vehicle with the RT + RZ group showed gene expression signatures of neuroprotection of hippocampal excitatory neurons post-RZ. RT-exposed mice performed poorly on learning and memory, and memory consolidation tasks. However, irradiated mice receiving RZ (13 mg/kg, drinking water) for 6-7 weeks showed a significant improvement in cognitive function compared to RT-exposed mice receiving vehicle. Dual-immunofluorescence staining, spatial transcriptomics, and biochemical assessment of RZ-treated irradiated brains demonstrated preservation of synaptic integrity and mature neuronal plasticity but not neurogenesis and reduced neuroinflammation concurrent with elevated BDNF levels and transcripts compared to vehicle-treated irradiated brains. In summary, oral administration of RZ represents a viable and translationally feasible neuroprotective approach against RICD.

Omics Approaches to Investigate the Pathogenesis of Suicide

Suicide is the second leading cause of death in U.S. adolescents and young adults and is generally associated with a psychiatric disorder. Suicidal behavior has a complex etiology and pathogenesis. Moderate heritability suggests genetic causes. Associations between childhood and recent life adversity indicate contributions from epigenetic factors. Genomic contributions to suicide pathogenesis remain largely unknown. This article is based on a workshop held to design strategies to identify molecular drivers of suicide neurobiology that would be putative new treatment targets. The panel determined that while bulk tissue studies provide comprehensive information, single-nucleus approaches that identify cell type-specific changes are needed. While single-nuclei techniques lack information on cytoplasm, processes, spines, and synapses, spatial multiomic technologies on intact tissue detect cell alterations specific to brain tissue layers and subregions. Because suicide has genetic and environmental drivers, multiomic approaches that combine cell type-specific epigenome, transcriptome, and proteome provide a more complete picture of pathogenesis. To determine the direction of effect of suicide risk gene variants on RNA and protein expression and how these interact with epigenetic marks, single-nuclei and spatial multiomics quantitative trait loci maps should be integrated with whole-genome sequencing and genome-wide association databases. The workshop concluded with a recommendation for the formation of an international suicide biology consortium that will bring together brain banks and investigators with expertise in cutting-edge omics technologies to delineate the biology of suicide and identify novel potential treatment targets to be tested in cellular and animal models for drug and biomarker discovery to guide suicide prevention.

Transcriptomic landscape of mammalian ventral pallidum at single-cell resolution

The ventral pallidum (VP) is critical for motivated behaviors. While contemporary work has begun to elucidate the functional diversity of VP neurons, the molecular heterogeneity underlying this functional diversity remains incompletely understood. We used single-nucleus RNA sequencing and in situ hybridization to define the transcriptional taxonomy of VP cell types in mice, macaques, and baboons. We found transcriptional conservation between all three species, within the broader neurochemical cell types. Unique dopaminoceptive and cholinergic subclusters were identified and conserved across both primate species but had no homolog in mice. This harmonized consensus VP cellular atlas will pave the way for understanding the structure and function of the VP and identified key neuropeptides, neurotransmitters, and neurotransmitter receptors that could be targeted within specific VP cell types for functional investigations.

An integrated single-nucleus and spatial transcriptomics atlas reveals the molecular landscape of the human hippocampus

The hippocampus contains many unique cell types, which serve the structure's specialized functions, including learning, memory and cognition. These cells have distinct spatial organization, morphology, physiology, and connectivity, highlighting the importance of transcriptome-wide profiling strategies that retain cytoarchitectural organization. Here, we generated spatially-resolved transcriptomics (SRT) and single-nucleus RNA-sequencing (snRNA-seq) data from adjacent tissue sections of the anterior human hippocampus in ten adult neurotypical donors to define molecular profiles for hippocampal cell types and spatial domains. Using non-negative matrix factorization (NMF) and label transfer, we integrated these data by defining gene expression patterns within the snRNA-seq data and inferring their expression in the SRT data. We identified NMF patterns that captured transcriptional variation across neuronal cell types and indicated that the response of excitatory and inhibitory postsynaptic specializations were prioritized in different SRT spatial domains. We used the NMF and label transfer approach to leverage existing rodent datasets, identifying patterns of activity-dependent transcription and subpopulations of dentate gyrus granule cells in our SRT dataset that may be predisposed to participate in learning and memory ensembles. Finally, we characterized the spatial organization of NMF patterns corresponding to non-cornu ammonis pyramidal neurons and identified snRNA-seq clusters mapping to distinct regions of the retrohippocampus, to three subiculum layers, and to a population of presubiculum neurons. To make this comprehensive molecular atlas accessible to the scientific community, both raw and processed data are freely available, including through interactive web applications.

MerQuaCo: a computational tool for quality control in image-based spatial transcriptomics

Image-based spatial transcriptomics platforms are powerful tools often used to identify cell populations and describe gene expression in intact tissue. Spatial experiments return large, high-dimension datasets and several open-source software packages are available to facilitate analysis and visualization. Spatial results are typically imperfect. For example, local variations in transcript detection probability are common. Software tools to characterize imperfections and their impact on downstream analyses are lacking so the data quality is assessed manually, a laborious and often a subjective process. Here we describe imperfections in a dataset of 641 fresh-frozen adult mouse brain sections collected using the Vizgen MERSCOPE. Common imperfections included the local loss of tissue from the section, tissue outside the imaging volume due to detachment from the coverslip, transcripts missing due to dropped images, varying detection probability through space, and differences in transcript detection probability between experiments. We describe the incidence of each imperfection and the likely impact on the accuracy of cell type labels. We develop MerQuaCo, open-source code that detects and quantifies imperfections without user input, facilitating the selection of sections for further analysis with existing packages. Together, our results and MerQuaCo facilitate rigorous, objective assessment of the quality of spatial transcriptomics results.

The single-cell opioid responses in the context of HIV (SCORCH) consortium

Substance use disorders (SUD) and drug addiction are major threats to public health, impacting not only the millions of individuals struggling with SUD, but also surrounding families and communities. One of the seminal challenges in treating and studying addiction in human populations is the high prevalence of co-morbid conditions, including an increased risk of contracting a human immunodeficiency virus (HIV) infection. Of the ~15 million people who inject drugs globally, 17% are persons with HIV. Conversely, HIV is a risk factor for SUD because chronic pain syndromes, often encountered in persons with HIV, can lead to an increased use of opioid pain medications that in turn can increase the risk for opioid addiction. We hypothesize that SUD and HIV exert shared effects on brain cell types, including adaptations related to neuroplasticity, neurodegeneration, and neuroinflammation. Basic research is needed to refine our understanding of these affected cell types and adaptations. Studying the effects of SUD in the context of HIV at the single-cell level represents a compelling strategy to understand the reciprocal interactions among both conditions, made feasible by the availability of large, extensively-phenotyped human brain tissue collections that have been amassed by the Neuro-HIV research community. In addition, sophisticated animal models that have been developed for both conditions provide a means to precisely evaluate specific exposures and stages of disease. We propose that single-cell genomics is a uniquely powerful technology to characterize the effects of SUD and HIV in the brain, integrating data from human cohorts and animal models. We have formed the Single-Cell Opioid Responses in the Context of HIV (SCORCH) consortium to carry out this strategy.

Opportunities and challenges of single-cell and spatially resolved genomics methods for neuroscience discovery

Over the past decade, single-cell genomics technologies have allowed scalable profiling of cell-type-specific features, which has substantially increased our ability to study cellular diversity and transcriptional programs in heterogeneous tissues. Yet our understanding of mechanisms of gene regulation or the rules that govern interactions between cell types is still limited. The advent of new computational pipelines and technologies, such as single-cell epigenomics and spatially resolved transcriptomics, has created opportunities to explore two new axes of biological variation: cell-intrinsic regulation of cell states and expression programs and interactions between cells. Here, we summarize the most promising and robust technologies in these areas, discuss their strengths and limitations and discuss key computational approaches for analysis of these complex datasets. We highlight how data sharing and integration, documentation, visualization and benchmarking of results contribute to transparency, reproducibility, collaboration and democratization in neuroscience, and discuss needs and opportunities for future technology development and analysis.

Neonatal Brain Injury Triggers Niche-Specific Changes to Cellular Biogeography

Preterm infants are at risk for brain injury and neurodevelopmental impairment due, in part, to white matter injury following chronic hypoxia exposure. However, the precise molecular mechanisms by which neonatal hypoxia disrupts early neurodevelopment are poorly understood. Here, we constructed a brain-wide map of the regenerative response to newborn brain injury using high-resolution imaging-based spatial transcriptomics to analyze over 800,000 cells in a mouse model of chronic neonatal hypoxia. Additionally, we developed a new method for inferring condition-associated differences in cell type spatial proximity, enabling the identification of niche-specific changes in cellular architecture. We observed hypoxia-associated changes in region-specific cell states, cell type composition, and spatial organization. Importantly, our analysis revealed mechanisms underlying reparative neurogenesis and gliogenesis, while also nominating pathways that may impede circuit rewiring following neonatal hypoxia. Altogether, our work provides a comprehensive description of the molecular response to newborn brain injury.

MYT1L deficiency impairs excitatory neuron trajectory during cortical development

Mutations reducing the function of MYT1L, a neuron-specific transcription factor, are associated with a syndromic neurodevelopmental disorder. MYT1L is used as a pro-neural factor in fibroblast-to-neuron transdifferentiation and is hypothesized to influence neuronal specification and maturation, but it is not clear which neuron types are most impacted by MYT1L loss. In this study, we profile 412,132 nuclei from the forebrains of wild-type and MYT1L-deficient mice at three developmental stages: E14 at the peak of neurogenesis, P1 when cortical neurons have been born, and P21 when neurons are maturing, to examine the role of MYT1L levels on neuronal development. MYT1L deficiency disrupts cortical neuron proportions and gene expression, primarily affecting neuronal maturation programs. Effects are mostly cell autonomous and persistent through development. While MYT1L can both activate and repress gene expression, the repressive effects are most sensitive to haploinsufficiency, likely mediating MYT1L syndrome. These findings illuminate MYT1L's role in orchestrating gene expression during neuronal development, providing insights into the molecular underpinnings of MYT1L syndrome.

Post-Critical Period Transcriptional and Physiological Adaptations of Thalamocortical Connections after Sensory Loss

Unilateral whisker denervation activates plasticity mechanisms and circuit adaptations in adults. Single nucleus RNA sequencing and multiplex fluorescence in situ hybridization revealed differentially expressed genes related to altered glutamate receptor distributions and synaptogenesis in thalamocortical (TC) recipient layer 4 (L4) neurons of the sensory cortex, specifically those receiving input from the intact whiskers after whisker denervation. Electrophysiology detected increased spontaneous excitatory events at L4 neurons, confirming an increase in synaptic connections. Elevated expression levels of Gria2 mRNA and functional GluA2 subunit of AMPA receptors at the TC synapse indicate the presence of stabilized and potentiated TC synapses to L4 excitatory neurons along the intact pathway after unilateral whisker denervation. These adaptations likely underlie the increased cortical activity observed in rodents during intact whisker sensation after unilateral whisker denervation. Our findings provide new insights into the mechanisms by which the adult brain supports recovery after unilateral sensory loss.

Characterization of direct Purkinje cell outputs to the brainstem

Purkinje cells (PCs) primarily project to cerebellar nuclei but also directly innervate the brainstem. Some PC-brainstem projections have been described previously, but most have not been thoroughly characterized. Here we use a PC-specific cre line to anatomically and electrophysiologically characterize PC projections to the brainstem. PC synapses are surprisingly widespread, with the highest densities found in the vestibular and parabrachial nuclei. However, there are pronounced regional differences in synaptic densities within both the vestibular and parabrachial nuclei. Large optogenetically-evoked PC-IPSCs are preferentially observed in subregions with the highest densities of PC synapses, suggesting that PCs selectively influence these areas and the behaviors they regulate. Unexpectedly, the pontine central gray and nearby subnuclei also contained a low density of PC synapses, and large PC-IPSCs are observed in a small fraction of cells. We combined electrophysiological recordings with immunohistochemistry to assess the molecular identities of two putative PC targets: PC synapses onto mesencephalic trigeminal neurons were not observed even though these cells are in close proximity to PC boutons. PC synapses onto locus coeruleus neurons are exceedingly rare or absent, even though previous studies concluded that PCs are a major input to these neurons. The availability of a highly selective cre line for PCs allowed us to study functional synapses, while avoiding complications that can accompany the use of viral approaches. We conclude that PCs directly innervate numerous brainstem nuclei, but only inhibit a small fraction of cells in many nuclei. This suggests that PCs target cell types with specific behavioral roles in brainstem regions.

Widespread co-release of glutamate and GABA throughout the mouse brain

Several brain neuronal populations transmit both the excitatory and inhibitory neurotransmitters, glutamate, and GABA. However, it remains largely unknown whether these opposing neurotransmitters are co-released simultaneously or are independently transmitted at different times and locations. By recording from acute mouse brain slices, we observed biphasic miniature postsynaptic currents, i.e., minis with time-locked excitatory and inhibitory currents, in striatal spiny projection neurons. This observation cannot be explained by accidental coincidence of monophasic excitatory and inhibitory minis. Interestingly, these biphasic minis could either be an excitatory current leading an inhibitory current or vice versa. Deletion of dopaminergic neurons did not eliminate biphasic minis, indicating that they originate from another source. Importantly, we found that both types of biphasic minis were present in multiple striatal neuronal types and in nine out of ten other brain regions. Overall, co-release of glutamate and GABA appears to be a widespread mode of neurotransmission in the brain.

Dual decoding of cell types and gene expression in spatial transcriptomics with PANDA

Sequencing-based spatial transcriptomics technologies have revolutionized our understanding of complex biological systems by enabling transcriptome profiling while preserving spatial context. However, spot-level expression measurements often amalgamate signals from diverse cells, obscuring potential heterogeneity. Existing methods aim to deconvolute spatial transcriptomics data into cell type proportions for each spot using single-cell RNA sequencing references but overlook cell-type-specific gene expression, essential for uncovering intra-type heterogeneity. We present PANDA (ProbAbilistic-based decoNvolution with spot-aDaptive cell type signAtures), a novel method that concurrently deciphers spot-level gene expression into both cell type proportions and cell-type-specific gene expression. PANDA integrates archetypal analysis to capture within-cell-type heterogeneity and dynamically learns cell type signatures for each spot during deconvolution. Simulations demonstrate PANDA's superior performance. Applied to real spatial transcriptomics data from diverse tissues, including tumor, brain, and developing heart, PANDA reconstructs spatial structures and reveals subtle transcriptional variations within specific cell types, offering a comprehensive understanding of tissue dynamics.

Generation of Mammalian Astrocyte Functional Heterogeneity

Mammalian astrocytes have regional roles within the brain parenchyma. Indeed, the notion that astrocytes are molecularly heterogeneous could help explain how the central nervous system (CNS) retains embryonic positional information through development into specialized regions into adulthood. A growing body of evidence supports the concept of morphological and molecular differences between astrocytes in different brain regions, which might relate to their derivation from regionally patterned radial glia and/or local neuron inductive cues. Here, we review evidence for regionally encoded functions of astrocytes to provide an integrated concept on lineage origins and heterogeneity to understand regional brain organization, as well as emerging technologies to identify and further investigate novel roles for astrocytes.

Reimagining Cortical Connectivity by Deconstructing Its Molecular Logic into Building Blocks

Comprehensive maps of neuronal connectivity provide a foundation for understanding the structure of neural circuits. In a circuit, neurons are diverse in morphology, electrophysiology, gene expression, activity, and other neuronal properties. Thus, constructing a comprehensive connectivity map requires associating various properties of neurons, including their connectivity, at cellular resolution. A commonly used approach is to use the gene expression profiles as an anchor to which all other neuronal properties are associated. Recent advances in genomics and anatomical techniques dramatically improved the ability to determine and associate the long-range projections of neurons with their gene expression profiles. These studies revealed unprecedented details of the gene-projection relationship, but also highlighted conceptual challenges in understanding this relationship. In this article, I delve into the findings and the challenges revealed by recent studies using state-of-the-art neuroanatomical and transcriptomic techniques. Building upon these insights, I propose an approach that focuses on understanding the gene-projection relationship through basic features in gene expression profiles and projections, respectively, that associate with underlying cellular processes. I then discuss how the developmental trajectories of projections and gene expression profiles create additional challenges and necessitate interrogating the gene-projection relationship across time. Finally, I explore complementary strategies that, together, can provide a comprehensive view of the gene-projection relationship.

Decoder-seq enhances mRNA capture efficiency in spatial RNA sequencing

Spatial transcriptomics technologies with high resolution often lack high sensitivity in mRNA detection. Here we report a dendrimeric DNA coordinate barcoding design for spatial RNA sequencing (Decoder-seq), which offers both high sensitivity and high resolution. Decoder-seq combines dendrimeric nanosubstrates with microfluidic coordinate barcoding to generate spatial arrays with a DNA density approximately ten times higher than previously reported methods while maintaining flexibility in resolution. We show that the high RNA capture efficiency of Decoder-seq improved the detection of lowly expressed olfactory receptor (Olfr) genes in mouse olfactory bulbs and contributed to the discovery of a unique layer enrichment pattern for two Olfr genes. The near-cellular resolution provided by Decoder-seq has enabled the construction of a spatial single-cell atlas of the mouse hippocampus, revealing dendrite-enriched mRNAs in neurons. When applying Decoder-seq to human renal cell carcinomas, we dissected the heterogeneous tumor microenvironment across different cancer subtypes and identified spatial gradient-expressed genes related to epithelial-mesenchymal transition with the potential to predict tumor prognosis and progression.

VAPOR: Variational autoencoder with transport operators decouples co-occurring biological processes in development

Background

Emerging single-cell and spatial transcriptomic data enable the investigation of gene expression dynamics of various biological processes, especially for development. To this end, existing computational methods typically infer trajectories that sequentially order cells for revealing gene expression changes in development, e.g., to assign a pseudotime to each cell indicating the ordering. However, these trajectories can aggregate different biological processes that cells undergo simultaneously-such as maturation for specialized function and differentiation into specific cell types-which do not occur on the same timescale. Therefore, a single pseudotime axis may not distinguish gene expression dynamics from co-occurring processes.

Methods

We introduce a method, VAPOR (variational autoencoder with transport operators), to decouple dynamic patterns from developmental gene expression data. Particularly, VAPOR learns a latent space for gene expression dynamics and decomposes the space into multiple subspaces. The dynamics on each subspace are governed by an ordinary differential equation model, attempting to recapitulate specific biological processes. Furthermore, we can infer the process-specific pseudotimes, revealing multifaceted timescales of distinct processes in which cells may simultaneously be involved during development.

Results

Initially tested on simulated datasets, VAPOR effectively recovered the topology and decoupled distinct dynamic patterns in the data. We then applied VAPOR to a developmental human brain scRNA-seq dataset across postconceptional weeks and identified gene expression dynamics for several key processes, such as differentiation and maturation. Moreover, our benchmarking analyses also demonstrated the outperformance of VAPOR over other methods. Additionally, we applied VAPOR to spatial transcriptomics data in the human dorsolateral prefrontal cortex. VAPOR captured the 'inside-out' pattern across cortical layers, potentially revealing how layers were formed, characterized by their gene expression dynamics.

Conclusion

VAPOR is open source for general use ( https://github.com/daifengwanglab/VAPOR ) to parameterize and infer developmental gene expression dynamics. It can be further extended for other single-cell and spatial omics such as chromatin accessibility to reveal developmental epigenomic dynamics.

Control of innate olfactory valence by segregated cortical amygdala circuits

Animals exhibit innate behaviors that are stereotyped responses to specific evolutionarily relevant stimuli in the absence of prior learning or experience. These behaviors can be reduced to an axis of valence, whereby specific odors evoke approach or avoidance responses. The posterolateral cortical amygdala (plCoA) mediates innate attraction and aversion to odor. However, little is known about how this brain area gives rise to behaviors of opposing motivational valence. Here, we sought to define the circuit features of plCoA that give rise to innate attraction and aversion to odor. We characterized the physiology, gene expression, and projections of this structure, identifying a divergent, topographic organization that selectively controls innate attraction and avoidance to odor. First, we examined odor-evoked responses in these areas and found sparse encoding of odor identity, but not valence. We next considered a topographic organization and found that optogenetic stimulation of the anterior and posterior domains of plCoA elicits attraction and avoidance, respectively, suggesting a functional axis for valence. Using single cell and spatial RNA sequencing, we identified the molecular cell types in plCoA, revealing an anteroposterior gradient in cell types, whereby anterior glutamatergic neurons preferentially express VGluT2 and posterior neurons express VGluT1. Activation of these respective cell types recapitulates appetitive and aversive behaviors, and chemogenetic inhibition reveals partial necessity for responses to innate appetitive or aversive odors. Finally, we identified topographically organized circuits defined by projections, whereby anterior neurons preferentially project to medial amygdala, and posterior neurons preferentially project to nucleus accumbens, which are respectively sufficient and necessary for innate attraction and aversion. Together, these data advance our understanding of how the olfactory system generates stereotypic, hardwired attraction and avoidance, and supports a model whereby distinct, topographically distributed plCoA populations direct innate olfactory responses by signaling to divergent valence-specific targets, linking upstream olfactory identity to downstream valence behaviors, through a population code. This suggests a novel amygdala circuit motif in which valence encoding is represented not by the firing properties of individual neurons, but by population level identity encoding that is routed through divergent targets to mediate distinct behaviors of opposing appetitive and aversive responses.

Network analysis of α-synuclein pathology progression reveals p21-activated kinases as regulators of vulnerability

α-Synuclein misfolding and progressive accumulation drives a pathogenic process in Parkinson's disease. To understand cellular and network vulnerability to α-synuclein pathology, we developed a framework to quantify network-level vulnerability and identify new therapeutic targets at the cellular level. Full brain α-synuclein pathology was mapped in mice over 9 months. Empirical pathology data was compared to theoretical pathology estimates from a diffusion model of pathology progression along anatomical connections. Unexplained variance in the model enabled us to derive regional vulnerability that we compared to regional gene expression. We identified gene expression patterns that relate to regional vulnerability, including 12 kinases that were enriched in vulnerable regions. Among these, an inhibitor of group II PAKs demonstrated protection from neuron death and α-synuclein pathology, even after delayed compound treatment. This study provides a framework for the derivation of cellular vulnerability from network-based studies and identifies a promising therapeutic pathway for Parkinson's disease.

Developmental mouse brain common coordinate framework

3D brain atlases are key resources to understand the brain's spatial organization and promote interoperability across different studies. However, unlike the adult mouse brain, the lack of developing mouse brain 3D reference atlases hinders advancements in understanding brain development. Here, we present a 3D developmental common coordinate framework (DevCCF) spanning embryonic day (E)11.5, E13.5, E15.5, E18.5, and postnatal day (P)4, P14, and P56, featuring undistorted morphologically averaged atlas templates created from magnetic resonance imaging and co-registered high-resolution light sheet fluorescence microscopy templates. The DevCCF with 3D anatomical segmentations can be downloaded or explored via an interactive 3D web-visualizer. As a use case, we utilize the DevCCF to unveil GABAergic neuron emergence in embryonic brains. Moreover, we map the Allen CCFv3 and spatial transcriptome cell-type data to our stereotaxic P56 atlas. In summary, the DevCCF is an openly accessible resource for multi-study data integration to advance our understanding of brain development.

Transcriptomic mapping of the 5-HT receptor landscape

Serotonin (5-HT) is crucial for regulating brain functions such as mood, sleep, and cognition. This study presents a comprehensive transcriptomic analysis of 5-HT receptors (Htrs) across ≈4 million cells in the adult mouse brain using single-cell RNA sequencing (scRNA-seq) data from the Allen Institute. We observed differential transcription patterns of all 14 Htr subtypes, revealing diverse prevalence and distribution across cell classes. Remarkably, we found that 65.84% of cells transcribe RNA of at least one Htr, with frequent co-transcription of multiple Htrs, underscoring the complexity of the 5-HT system even at the single-cell dimension. Leveraging a multiplexed error-robust fluorescence in situ hybridization (MERFISH) dataset provided by Harvard University of ≈10 million cells, we analyzed the spatial distribution of each Htr, confirming previous findings and uncovering novel transcription patterns. To aid in exploring Htr transcription, we provide an online interactive visualizer.

Cell class-specific long-range axonal projections of neurons in mouse whisker-related somatosensory cortices

Long-range axonal projections of diverse classes of neocortical excitatory neurons likely contribute to brain-wide interactions processing sensory, cognitive and motor signals. Here, we performed light-sheet imaging of fluorescently labeled axons from genetically defined neurons located in posterior primary somatosensory barrel cortex and supplemental somatosensory cortex. We used convolutional networks to segment axon-containing voxels and quantified their distribution within the Allen Mouse Brain Atlas Common Coordinate Framework. Axonal density was analyzed for different classes of glutamatergic neurons using transgenic mouse lines selectively expressing Cre recombinase in layer 2/3 intratelencephalic projection neurons (Rasgrf2-dCre), layer 4 intratelencephalic projection neurons (Scnn1a-Cre), layer 5 intratelencephalic projection neurons (Tlx3-Cre), layer 5 pyramidal tract projection neurons (Sim1-Cre), layer 5 projection neurons (Rbp4-Cre), and layer 6 corticothalamic neurons (Ntsr1-Cre). We found distinct axonal projections from the different neuronal classes to many downstream brain areas, which were largely similar for primary and supplementary somatosensory cortices. Functional connectivity maps obtained from optogenetic activation of sensory cortex and wide-field imaging revealed topographically organized evoked activity in frontal cortex with neurons located more laterally in somatosensory cortex signaling to more anteriorly located regions in motor cortex, consistent with the anatomical projections. The current methodology therefore appears to quantify brain-wide axonal innervation patterns supporting brain-wide signaling.

Spatiotemporal modeling reveals high-resolution invasion states in glioblastoma

BACKGROUND

Diffuse invasion of glioblastoma cells through normal brain tissue is a key contributor to tumor aggressiveness, resistance to conventional therapies, and dismal prognosis in patients. A deeper understanding of how components of the tumor microenvironment (TME) contribute to overall tumor organization and to programs of invasion may reveal opportunities for improved therapeutic strategies.

RESULTS

Towards this goal, we apply a novel computational workflow to a spatiotemporally profiled GBM xenograft cohort, leveraging the ability to distinguish human tumor from mouse TME to overcome previous limitations in the analysis of diffuse invasion. Our analytic approach, based on unsupervised deconvolution, performs reference-free discovery of cell types and cell activities within the complete GBM ecosystem. We present a comprehensive catalogue of 15 tumor cell programs set within the spatiotemporal context of 90 mouse brain and TME cell types, cell activities, and anatomic structures. Distinct tumor programs related to invasion align with routes of perivascular, white matter, and parenchymal invasion. Furthermore, sub-modules of genes serving as program network hubs are highly prognostic in GBM patients.

CONCLUSION

The compendium of programs presented here provides a basis for rational targeting of tumor and/or TME components. We anticipate that our approach will facilitate an ecosystem-level understanding of the immediate and long-term consequences of such perturbations, including the identification of compensatory programs that will inform improved combinatorial therapies.

Spatial Single-Cell Mapping of Transcriptional Differences Across Genetic Backgrounds in Mouse Brains

Genetic variation can alter brain structure and, consequently, function. Comparative statistical analysis of mouse brains across genetic backgrounds requires spatial, single-cell, atlas-scale data, in replicates-a challenge for current technologies. We introduce Atlas-scale Transcriptome Localization using Aggregate Signatures (ATLAS), a scalable tissue mapping method. ATLAS learns transcriptional signatures from scRNAseq data, encodes them in situ with tens of thousands of oligonucleotide probes, and decodes them to infer cell types and imputed transcriptomes. We validated ATLAS by comparing its cell type inferences with direct MERFISH measurements of marker genes and quantitative comparisons to four other technologies. Using ATLAS, we mapped the central brains of five male and five female C57BL/6J (B6) mice and five male BTBR T+ tf/J (BTBR) mice, an idiopathic model of autism, collectively profiling over 40 million cells across over 400 coronal sections. Our analysis revealed over 40 significant differences in cell type distributions and identified 16 regional composition changes across male-female and B6-BTBR comparisons. ATLAS thus enables systematic comparative studies, facilitating organ-level structure-function analysis of disease models.

Molecular states underlying neuronal cell type development and plasticity in the whisker cortex

Mouse whisker somatosensory cortex (wS1) is a major model system to study the experience-dependent plasticity of cortical neuron physiology, morphology, and sensory coding. However, the role of sensory experience in regulating neuronal cell type development and gene expression in wS1 remains poorly understood. We assembled and annotated a transcriptomic atlas of wS1 during postnatal development comprising 45 molecularly distinct neuronal types that can be grouped into eight excitatory and four inhibitory neuron subclasses. Using this atlas, we examined the influence of whisker experience from postnatal day (P) 12, the onset of active whisking, to P22, on the maturation of molecularly distinct cell types. During this developmental period, when whisker experience was normal, ~250 genes were regulated in a neuronal subclass-specific fashion. At the resolution of neuronal types, we found that only the composition of layer (L) 2/3 glutamatergic neuronal types, but not other neuronal types, changed substantially between P12 and P22. These compositional changes resemble those observed previously in the primary visual cortex (V1), and the temporal gene expression changes were also highly conserved between the two regions. In contrast to V1, however, cell type maturation in wS1 is not substantially dependent on sensory experience, as 10-day full-face whisker deprivation did not influence the transcriptomic identity and composition of L2/3 neuronal types. A one-day competitive whisker deprivation protocol also did not affect cell type identity but induced moderate changes in plasticity-related gene expression. Thus, developmental maturation of cell types is similar in V1 and wS1, but sensory deprivation minimally affects cell type development in wS1.

Control of sodium appetite by hindbrain aldosterone-sensitive neurons

Mineralocorticoids play a key role in hydromineral balance by regulating sodium retention and potassium wasting. Through favoring sodium, mineralocorticoids can cause hypertension from fluid overload under conditions of hyperaldosteronism, such as aldosterone-secreting tumors. An often-overlooked mechanism by which aldosterone functions to increase sodium is through stimulation of salt appetite. To drive sodium intake, aldosterone targets neurons in the hindbrain which uniquely express 11β-hydroxysteroid dehydrogenase type 2 (HSD2). This enzyme is a necessary precondition for aldosterone-sensing cells as it metabolizes glucocorticoids - preventing their activation of the mineralocorticoid receptor. In this review, we will consider the role of hindbrain HSD2 neurons in regulating sodium appetite by discussing HSD2 expression in the brain, regulation of hindbrain HSD2 neuron activity, and the circuitry mediating the effects of these aldosterone-sensitive neurons. Reducing the activity of hindbrain HSD2 neurons may be a viable strategy to reduce sodium intake and cardiovascular risk, particularly for conditions of hyperaldosteronism.