A taxonomy of transcriptomic cell types across the isocortex and hippocampal formation

A single-nucleus transcriptomic atlas of medium spiny neurons in the rat nucleus accumbens

Neural processing of rewarding stimuli involves several distinct regions, including the nucleus accumbens (NAc). The majority of NAc neurons are GABAergic projection neurons known as medium spiny neurons (MSNs). MSNs are broadly defined by dopamine receptor expression, but evidence suggests that a wider array of subtypes exist. To study MSN heterogeneity, we analyzed single-nucleus RNA sequencing data from the largest available rat NAc dataset. Analysis of 48,040 NAc MSN nuclei identified major populations belonging to the striosome and matrix compartments. Integration with mouse and human data indicated consistency across species and disease-relevance scoring using genome-wide association study results revealed potentially differential roles for MSN populations in substance use disorders. Additional high-resolution clustering identified 34 transcriptomically distinct subtypes of MSNs definable by a limited number of marker genes. Together, these data demonstrate the diversity of MSNs in the NAc and provide a basis for more targeted genetic manipulation of specific populations.

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.

Visualizing scRNA-Seq Data at Population Scale with GloScope

Increasingly, scRNA-Seq studies explore cell populations across different samples and the effect of sample heterogeneity on organism's phenotype. However, relatively few bioinformatic methods have been developed which adequately address the variation between samples for such population-level analyses. We propose a framework for representing the entire single-cell profile of a sample, which we call a GloScope representation. We implement GloScope on scRNA-Seq datasets from study designs ranging from 12 to over 300 samples and demonstrate how GloScope allows researchers to perform essential bioinformatic tasks at the sample-level, in particular visualization and quality control assessment.

Integrative proteomics identifies a conserved Aβ amyloid responsome, novel plaque proteins, and pathology modifiers in Alzheimer's disease

Alzheimer's disease (AD) is a complex neurodegenerative disorder that develops over decades. AD brain proteomics reveals vast alterations in protein levels and numerous altered biologic pathways. Here, we compare AD brain proteome and network changes with the brain proteomes of amyloid β (Aβ)-depositing mice to identify conserved and divergent protein networks with the conserved networks identifying an Aβ amyloid responsome. Proteins in the most conserved network (M42) accumulate in plaques, cerebrovascular amyloid (CAA), and/or dystrophic neuronal processes, and overexpression of two M42 proteins, midkine (Mdk) and pleiotrophin (PTN), increases the accumulation of Aβ in plaques and CAA. M42 proteins bind amyloid fibrils in vitro, and MDK and PTN co-accumulate with cardiac transthyretin amyloid. M42 proteins appear intimately linked to amyloid deposition and can regulate amyloid deposition, suggesting that they are pathology modifiers and thus putative therapeutic targets. We posit that amyloid-scaffolded accumulation of numerous M42+ proteins is a central mechanism mediating downstream pathophysiology in AD.

Opposing effects of the purinergic P2X7 receptor on seizures in neurons and microglia in male mice

BACKGROUND

The purinergic ATP-gated P2X7 receptor (P2X7R) is increasingly recognized to contribute to pathological neuroinflammation and brain hyperexcitability. P2X7R expression has been shown to be increased in the brain, including both microglia and neurons, in experimental models of epilepsy and patients. To date, the cell type-specific downstream effects of P2X7Rs during seizures remain, however, incompletely understood.

METHODS

Effects of P2X7R signaling on seizures and epilepsy were analyzed in induced seizure models using male mice including the kainic acid model of status epilepticus and pentylenetetrazole model and in male and female mice in a genetic model of Dravet syndrome. RNA sequencing was used to analyze P2X7R downstream signaling during seizures. To investigate the cell type-specific role of the P2X7R during seizures and epilepsy, we generated mice lacking exon 2 of the P2rx7 gene in either microglia (P2rx7:Cx3cr1-Cre) or neurons (P2rx7:Thy-1-Cre). To investigate the protective potential of overexpressing P2X7R in GABAergic interneurons, P2X7Rs were overexpressed using adeno-associated virus transduction under the mDlx promoter.

RESULTS

RNA sequencing of hippocampal tissue from wild-type and P2X7R knock-out mice identified both glial and neuronal genes, in particular genes involved in GABAergic signaling, under the control of the P2X7R following seizures. Mice with deleted P2rx7 in microglia displayed less severe acute seizures and developed a milder form of epilepsy, and microglia displayed an anti-inflammatory molecular profile. In contrast, mice lacking P2rx7 in neurons showed a more severe seizure phenotype when compared to epileptic wild-type mice. Analysis of single-cell expression data revealed that human P2RX7 expression is elevated in the hippocampus of patients with temporal lobe epilepsy in excitatory and inhibitory neurons. Functional studies determined that GABAergic interneurons display increased responses to P2X7R activation in experimental epilepsy. Finally, we show that viral transduction of P2X7R in GABAergic interneurons protects against evoked and spontaneous seizures in experimental temporal lobe epilepsy and in mice lacking Scn1a, a model of Dravet syndrome.

CONCLUSIONS

Our results suggest a dual and opposing action of P2X7R in epilepsy and suggest P2X7R overexpression in GABAergic interneurons as a novel therapeutic strategy for acquired and, possibly, genetic forms of epilepsy.

Spatial enrichment and genomic analyses reveal the link of NOMO1 with amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a severe motor neuron disease with uncertain genetic predisposition in most sporadic cases. The spatial architecture of cell types and gene expression are the basis of cell-cell interactions, biological function and disease pathology, but are not well investigated in the human motor cortex, a key ALS-relevant brain region. Recent studies indicated single nucleus transcriptomic features of motor neuron vulnerability in ALS motor cortex. However, the brain regional vulnerability of ALS-associated genes and the genetic link between region-specific genes and ALS risk remain largely unclear. Here, we developed an entropy-weighted differential gene expression matrix-based tool (SpatialE) to identify the spatial enrichment of gene sets in spatial transcriptomics. We benchmarked SpatialE against another enrichment tool (multimodal intersection analysis) using spatial transcriptomics data from both human and mouse brain tissues. To investigate regional vulnerability, we analysed three human motor cortex and two dorsolateral prefrontal cortex tissues for spatial enrichment of ALS-associated genes. We also used Cell2location to estimate the abundance of cell types in ALS-related cortex layers. To dissect the link of regionally expressed genes and ALS risk, we performed burden analyses of rare loss-of-function variants detected by whole-genome sequencing in ALS patients and controls, then analysed differential gene expression in the TargetALS RNA-sequencing dataset. SpatialE showed more accurate and specific spatial enrichment of regional cell type markers than multimodal intersection analysis in both mouse brain and human dorsolateral prefrontal cortex. Spatial transcriptomic analyses of human motor cortex showed heterogeneous cell types and spatial gene expression profiles. We found that 260 manually curated ALS-associated genes are significantly enriched in layer 5 of the motor cortex, with abundant expression of upper motor neurons and layer 5 excitatory neurons. Burden analyses of rare loss-of-function variants in Layer 5-associated genes nominated NOMO1 as a novel ALS-associated gene in a combined sample set of 6814 ALS patients and 3324 controls (P = 0.029). Gene expression analyses in CNS tissues revealed downregulation of NOMO1 in ALS, which is consistent with a loss-of-function disease mechanism. In conclusion, our integrated spatial transcriptomics and genomic analyses identified regional brain vulnerability in ALS and the association of a layer 5 gene (NOMO1) with ALS risk.

Optimized ultrasound neuromodulation for non-invasive control of behavior and physiology

Focused ultrasound can non-invasively modulate neural activity, but whether effective stimulation parameters generalize across brain regions and cell types remains unknown. We used focused ultrasound coupled with fiber photometry to identify optimal neuromodulation parameters for four different arousal centers of the brain in an effort to yield overt changes in behavior. Applying coordinate descent, we found that optimal parameters for excitation or inhibition are highly distinct, the effects of which are generally conserved across brain regions and cell types. Optimized stimulations induced clear, target-specific behavioral effects, whereas non-optimized protocols of equivalent energy resulted in substantially less or no change in behavior. These outcomes were independent of auditory confounds and, contrary to expectation, accompanied by a cyclooxygenase-dependent and prolonged reduction in local blood flow and temperature with brain-region-specific scaling. These findings demonstrate that carefully tuned and targeted ultrasound can exhibit powerful effects on complex behavior and physiology.

Nanoligomers targeting NF-κB and NLRP3 reduce neuroinflammation and improve cognitive function with aging and tauopathy

Neuroinflammation contributes to impaired cognitive function in brain aging and neurodegenerative disorders like Alzheimer's disease, which is characterized by the aggregation of pathological tau. One major driver of both age- and tau-associated neuroinflammation is the NF-κB and NLRP3 signaling axis. However, current treatments targeting NF-κB or NLRP3 may have adverse/systemic effects, and most have not been clinically translatable. In this study, we tested the efficacy of a novel, nucleic acid therapeutic (Nanoligomer) cocktail specifically targeting both NF-κB and NLRP3 in the brain for reducing neuroinflammation and improving cognitive function in old (aged 19 months) wildtype mice, and in rTg4510 tau pathology mice (aged 2 months). We found that 4 weeks of NF-κB/NLRP3-targeting Nanoligomer treatment strongly reduced neuro-inflammatory cytokine profiles in the brain and improved cognitive-behavioral function in both old and rTg4510 mice. These effects of NF-κB/NLRP3-targeting Nanoligomers were also associated with reduced glial cell activation and pathology, favorable changes in transcriptome signatures of glia-associated inflammation (reduced) and neuronal health (increased), and positive systemic effects. Collectively, our results provide a basis for future translational studies targeting both NF-κB and NLRP3 in the brain, perhaps using Nanoligomers, to inhibit neuroinflammation and improve cognitive function with aging and neurodegeneration.

Genetic approaches to elucidating cortical and hippocampal GABAergic interneuron diversity

GABAergic interneurons (INs) in the mammalian forebrain represent a diverse population of cells that provide specialized forms of local inhibition to regulate neural circuit activity. Over the last few decades, the development of a palette of genetic tools along with the generation of single-cell transcriptomic data has begun to reveal the molecular basis of IN diversity, thereby providing deep insights into how different IN subtypes function in the forebrain. In this review, we outline the emerging picture of cortical and hippocampal IN speciation as defined by transcriptomics and developmental origin and summarize the genetic strategies that have been utilized to target specific IN subtypes, along with the technical considerations inherent to each approach. Collectively, these methods have greatly facilitated our understanding of how IN subtypes regulate forebrain circuitry via cell type and compartment-specific inhibition and thus have illuminated a path toward potential therapeutic interventions for a variety of neurocognitive disorders.

Single-cell multiregion dissection of Alzheimer's disease

Alzheimer's disease is the leading cause of dementia worldwide, but the cellular pathways that underlie its pathological progression across brain regions remain poorly understood1-3. Here we report a single-cell transcriptomic atlas of six different brain regions in the aged human brain, covering 1.3 million cells from 283 post-mortem human brain samples across 48 individuals with and without Alzheimer's disease. We identify 76 cell types, including region-specific subtypes of astrocytes and excitatory neurons and an inhibitory interneuron population unique to the thalamus and distinct from canonical inhibitory subclasses. We identify vulnerable populations of excitatory and inhibitory neurons that are depleted in specific brain regions in Alzheimer's disease, and provide evidence that the Reelin signalling pathway is involved in modulating the vulnerability of these neurons. We develop a scalable method for discovering gene modules, which we use to identify cell-type-specific and region-specific modules that are altered in Alzheimer's disease and to annotate transcriptomic differences associated with diverse pathological variables. We identify an astrocyte program that is associated with cognitive resilience to Alzheimer's disease pathology, tying choline metabolism and polyamine biosynthesis in astrocytes to preserved cognitive function late in life. Together, our study develops a regional atlas of the ageing human brain and provides insights into cellular vulnerability, response and resilience to Alzheimer's disease pathology.

Transcriptional determinism and stochasticity contribute to the complexity of autism-associated SHANK family genes

Precision of transcription is critical because transcriptional dysregulation is disease causing. Traditional methods of transcriptional profiling are inadequate to elucidate the full spectrum of the transcriptome, particularly for longer and less abundant mRNAs. SHANK3 is one of the most common autism causative genes. Twenty-four Shank3-mutant animal lines have been developed for autism modeling. However, their preclinical validity has been questioned due to incomplete Shank3 transcript structure. We apply an integrative approach combining cDNA-capture and long-read sequencing to profile the SHANK3 transcriptome in humans and mice. We unexpectedly discover an extremely complex SHANK3 transcriptome. Specific SHANK3 transcripts are altered in Shank3-mutant mice and postmortem brain tissues from individuals with autism spectrum disorder. The enhanced SHANK3 transcriptome significantly improves the detection rate for potential deleterious variants from genomics studies of neuropsychiatric disorders. Our findings suggest that both deterministic and stochastic transcription of the genome is associated with SHANK family genes.

An enhancer-AAV toolbox to target and manipulate distinct interneuron subtypes

In recent years, we and others have identified a number of enhancers that, when incorporated into rAAV vectors, can restrict the transgene expression to particular neuronal populations. Yet, viral tools to access and manipulate fine neuronal subtypes are still limited. Here, we performed systematic analysis of single cell genomic data to identify enhancer candidates for each of the cortical interneuron subtypes. We established a set of enhancer-AAV tools that are highly specific for distinct cortical interneuron populations and striatal cholinergic neurons. These enhancers, when used in the context of different effectors, can target (fluorescent proteins), observe activity (GCaMP) and manipulate (opto- or chemo-genetics) specific neuronal subtypes. We also validated our enhancer-AAV tools across species. Thus, we provide the field with a powerful set of tools to study neural circuits and functions and to develop precise and targeted therapy.

Pyramidal neurons proportionately alter the identity and survival of specific cortical interneuron subtypes

The mammalian cerebral cortex comprises a complex neuronal network that maintains a delicate balance between excitatory neurons and inhibitory interneurons. Previous studies, including our own research, have shown that specific interneuron subtypes are closely associated with particular pyramidal neuron types, forming stereotyped local inhibitory microcircuits. However, the developmental processes that establish these precise networks are not well understood. Here we show that pyramidal neuron types are instrumental in driving the terminal differentiation and maintaining the survival of specific associated interneuron subtypes. In a wild-type cortex, the relative abundance of different interneuron subtypes aligns precisely with the pyramidal neuron types to which they synaptically target. In Fezf2 mutant cortex, characterized by the absence of layer 5 pyramidal tract neurons and an expansion of layer 6 intratelencephalic neurons, we observed a corresponding decrease in associated layer 5b interneurons and an increase in layer 6 subtypes. Interestingly, these shifts in composition are achieved through mechanisms specific to different interneuron types. While SST interneurons adjust their abundance to the change in pyramidal neuron prevalence through the regulation of programmed cell death, parvalbumin interneurons alter their identity. These findings illustrate two key strategies by which the dynamic interplay between pyramidal neurons and interneurons allows local microcircuits to be sculpted precisely. These insights underscore the precise roles of extrinsic signals from pyramidal cells in the establishment of interneuron diversity and their subsequent integration into local cortical microcircuits.

Distinct transcriptional programs define a heterogeneous neuronal ensemble for social interaction

Social interactions are encoded by the coordinated activity of heterogeneous cell types within distributed brain regions including the medial prefrontal cortex (mPFC). However, our understanding of the cell types which comprise the social ensemble has been limited by available mouse lines and reliance on single marker genes. We identified differentially active neuronal populations during social interactions by quantifying immediate-early gene (IEG) expression using snRNA-sequencing. These studies revealed that distinct prefrontal neuron populations composed of heterogeneous cell types are activated by social interaction. Evaluation of IEG expression within these recruited neuronal populations revealed cell-type and region-specific programs, suggesting that reliance on a single molecular marker is insufficient to quantify activation across all cell types. Our findings provide a comprehensive description of cell-type specific transcriptional programs invoked by social interactions and reveal insights into the neuronal populations which compose the social ensemble.

A brief history of somatostatin interneuron taxonomy or: how many somatostatin subtypes are there, really?

We provide a brief (and unabashedly biased) overview of the pre-transcriptomic history of somatostatin interneuron taxonomy, followed by a chronological summary of the large-scale, NIH-supported effort over the last ten years to generate a comprehensive, single-cell RNA-seq-based taxonomy of cortical neurons. Focusing on somatostatin interneurons, we present the perspective of experimental neuroscientists trying to incorporate the new classification schemes into their own research while struggling to keep up with the ever-increasing number of proposed cell types, which seems to double every two years. We suggest that for experimental analysis, the most useful taxonomic level is the subdivision of somatostatin interneurons into ten or so "supertypes," which closely agrees with their more traditional classification by morphological, electrophysiological and neurochemical features. We argue that finer subdivisions ("t-types" or "clusters"), based on slight variations in gene expression profiles but lacking clear phenotypic differences, are less useful to researchers and may actually defeat the purpose of classifying neurons to begin with. We end by stressing the need for generating novel tools (mouse lines, viral vectors) for genetically targeting distinct supertypes for expression of fluorescent reporters, calcium sensors and excitatory or inhibitory opsins, allowing neuroscientists to chart the input and output synaptic connections of each proposed subtype, reveal the position they occupy in the cortical network and examine experimentally their roles in sensorimotor behaviors and cognitive brain functions.

A multiomic atlas of the aging hippocampus reveals molecular changes in response to environmental enrichment

Aging involves the deterioration of organismal function, leading to the emergence of multiple pathologies. Environmental stimuli, including lifestyle, can influence the trajectory of this process and may be used as tools in the pursuit of healthy aging. To evaluate the role of epigenetic mechanisms in this context, we have generated bulk tissue and single cell multi-omic maps of the male mouse dorsal hippocampus in young and old animals exposed to environmental stimulation in the form of enriched environments. We present a molecular atlas of the aging process, highlighting two distinct axes, related to inflammation and to the dysregulation of mRNA metabolism, at the functional RNA and protein level. Additionally, we report the alteration of heterochromatin domains, including the loss of bivalent chromatin and the uncovering of a heterochromatin-switch phenomenon whereby constitutive heterochromatin loss is partially mitigated through gains in facultative heterochromatin. Notably, we observed the multi-omic reversal of a great number of aging-associated alterations in the context of environmental enrichment, which was particularly linked to glial and oligodendrocyte pathways. In conclusion, our work describes the epigenomic landscape of environmental stimulation in the context of aging and reveals how lifestyle intervention can lead to the multi-layered reversal of aging-associated decline.

Sex-biasing influence of autism-associated Ube3a gene overdosage at connectomic, behavioral, and transcriptomic levels

Genomic mechanisms enhancing risk in males may contribute to sex bias in autism. The ubiquitin protein ligase E3A gene (Ube3a) affects cellular homeostasis via control of protein turnover and by acting as transcriptional coactivator with steroid hormone receptors. Overdosage of Ube3a via duplication or triplication of chromosomal region 15q11-13 causes 1 to 2% of autistic cases. Here, we test the hypothesis that increased dosage of Ube3a may influence autism-relevant phenotypes in a sex-biased manner. We show that mice with extra copies of Ube3a exhibit sex-biasing effects on brain connectomics and autism-relevant behaviors. These effects are associated with transcriptional dysregulation of autism-associated genes, as well as genes differentially expressed in 15q duplication and in autistic people. Increased Ube3a dosage also affects expression of genes on the X chromosome, genes influenced by sex steroid hormone, and genes sex-differentially regulated by transcription factors. These results suggest that Ube3a overdosage can contribute to sex bias in neurodevelopmental conditions via influence on sex-differential mechanisms.

Cell type-specific and subcellular expression of phospholipid phosphatase-related proteins to modulate lyso-phosphatidic acid synaptic signaling in the developing and adult CNS

Lysophosphatidic acid (LPA) is a bioactive phospholipid that participates in critical processes in neural development and adult brain function and is implicated in various pathophysiological conditions. Along with its six well-characterized receptors, atypical regulators of LPA signaling have also been suggested, including phospholipid phosphatase-related proteins (PLPPRs). PLPPRs have been mostly studied in the developing brain where they control LPA-dependent axon guidance, cortical network hyperexcitability, and glutamatergic neurotransmission. PLPPR4 and PLPPR3 represent two closely related proteins reported to localize predominantly in dendrites and axons, respectively, and differ in their developmental expression patterns. Herein, we have revised the expression patterns of PLPPRs in the cerebellum, dorsal and ventral hippocampus, prefrontal cortex (PFC), nucleus accumbens, and striatum during development and in the adult using quantitative PCR. Expression patterns of Plppr2,4 and 5 were consistent with previous studies, whereas Plppr3 and Plppr1 exhibited a unique expression profile in nucleus accumbens (NAc) and striatum in later developmental and adult stages, which we verified at the protein level for PLPPR3. To investigate neuron type-specific expression at the single cell level, we developed a bioinformatic tool to analyze recent single-cell RNA-sequencing data in the cerebral cortex and hippocampus of adult mice. Our analysis revealed a widespread but also selective adult neuron-type expression with higher expression levels of Plppr3, Plppr1, and Plppr5 in GABAergic and Plppr4 and Plppr2 in glutamatergic neurons. PLPPR4 has been identified as a post-synaptic modulator of LPA levels in glutamatergic synapses operating via an uptake mechanism, to control LPA-dependent cortical network hyperexcitability. Using subcellular fractionation experiments, we found that both PLPPR4 and PLPPR3 are co-expressed in adult synaptosomal membranes. Furthermore, flow cytometry experiments in HEK293 cells showed comparable LPA uptake by PLPPR4 and PLPPR3, whereas PLPRR3, but not PLPPR4, induced also uptake of monoacylglycerol, the dephosphorylation product of LPA. We propose that synaptic LPA may be subject to both pre-synaptic and post-synaptic mechanisms of regulation by PLPPRs in addition to LPARs in developing and adult synapses.

The Pgb1 locus controls glycogen aggregation in astrocytes of the aged hippocampus without impacting cognitive function

In aged humans and mice, aggregates of hypobranched glycogen molecules called polyglucosan bodies (PGBs) accumulate in hippocampal astrocytes. PGBs are known to drive cognitive decline in neurological diseases but remain largely unstudied in the context of typical brain aging. Here, we show that PGBs arise in autophagy-dysregulated astrocytes of the aged C57BL/6J mouse hippocampus. To map the genetic cause of age-related PGB accumulation, we quantified PGB burden in 32 fully sequenced BXD-recombinant inbred mouse strains, which display a 400-fold variation in hippocampal PGB burden at 16-18 months of age. A major modifier locus was mapped to chromosome 1 at 72-75 Mb, which we defined as the Pgb1 locus. To evaluate candidate genes and downstream mechanisms by which Pgb1 controls the aggregation of glycogen, extensive hippocampal transcriptomic and proteomic datasets were produced for aged mice of the BXD family. We utilized these datasets to identify Smarcal1 and Usp37 as potential regulators of PGB accumulation. To assess the effect of PGB burden on age-related cognitive decline, we performed phenome-wide association scans, transcriptomic analyses as well as conditioned fear memory and Y-maze testing. Importantly, we did not find any evidence suggesting a negative impact of PGBs on cognition. Taken together, our study demonstrates that the Pgb1 locus controls glycogen aggregation in astrocytes of the aged hippocampus without affecting age-related cognitive decline.

Detection of allele-specific expression in spatial transcriptomics with spASE

Spatial transcriptomics technologies permit the study of the spatial distribution of RNA at near-single-cell resolution genome-wide. However, the feasibility of studying spatial allele-specific expression (ASE) from these data remains uncharacterized. Here, we introduce spASE, a computational framework for detecting and estimating spatial ASE. To tackle the challenges presented by cell type mixtures and a low signal to noise ratio, we implement a hierarchical model involving additive mixtures of spatial smoothing splines. We apply our method to allele-resolved Visium and Slide-seq from the mouse cerebellum and hippocampus and report new insight into the landscape of spatial and cell type-specific ASE therein.

Inhibitory neurons defined by a synaptogenic molecule impair memory discrimination

The CA3 region is central to hippocampal function during learning and memory because of its unique connectivity. CA3 pyramidal neurons are the targets of huge, excitatory mossy fiber synapses from DG axons and have an unusually high degree of excitatory recurrent connectivity. Thus, inhibition likely plays an outsized importance in constraining runaway excitation and shaping CA3 ensembles during learning and memory. Here, we investigate the function of a group of dendrite-targeting, hippocampal GABAergic neurons defined by expression of the synaptogenic adhesion molecule, Kirrel3. We discovered that activating Kirrel3-expressing GABAergic neurons impairs memory discrimination by inhibiting CA3 pyramidal neurons in novel contexts. Kirrel3 is required for DG-to-GABA synapse formation and variants in Kirrel3 are strong risk factors for neurodevelopmental disorders. Thus, our work suggests that Kirrel3-GABA neurons are a critical source of feed-forward inhibition from DG to CA3 during contextual memory whose activity may be specifically disrupted in some brain disorders.

Control of innate olfactory valence by segregated cortical amygdala circuits

Animals perform 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. The 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 olfactory behaviors of valence. 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 Slc17a6 and posterior neurons express Slc17a7. Activation of these respective cell types recapitulates appetitive and aversive valence behaviors, and chemogenetic inhibition reveals partial necessity for valence 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 negative and positive olfactory valence. 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 valence responses by signaling to divergent valence-specific targets, linking upstream olfactory identity to downstream valence behaviors, through a population code. This represents a novel 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 valence.

Making Ramón y Cajal proud: Development of cell identity and diversity in the cerebral cortex

Since the beautiful images of Santiago Ramón y Cajal provided a first glimpse into the immense diversity and complexity of cell types found in the cerebral cortex, neuroscience has been challenged and inspired to understand how these diverse cells are generated and how they interact with each other to orchestrate the development of this remarkable tissue. Some fundamental questions drive the field's quest to understand cortical development: what are the mechanistic principles that govern the emergence of neuronal diversity? How do extrinsic and intrinsic signals integrate with physical forces and activity to shape cell identity? How do the diverse populations of neurons and glia influence each other during development to guarantee proper integration and function? The advent of powerful new technologies to profile and perturb cortical development at unprecedented resolution and across a variety of modalities has offered a new opportunity to integrate past knowledge with brand new data. Here, we review some of this progress using cortical excitatory projection neurons as a system to draw out general principles of cell diversification and the role of cell-cell interactions during cortical development.

Complexin-1 enhances ultrasound neurotransmission in the mammalian auditory pathway

Unlike megabats, which rely on well-developed vision, microbats use ultrasonic echolocation to navigate and locate prey. To study ultrasound perception, here we compared the auditory cortices of microbats and megabats by constructing reference genomes and single-nucleus atlases for four species. We found that parvalbumin (PV)+ neurons exhibited evident cross-species differences and could respond to ultrasound signals, whereas their silencing severely affected ultrasound perception in the mouse auditory cortex. Moreover, megabat PV+ neurons expressed low levels of complexins (CPLX1-CPLX4), which can facilitate neurotransmitter release, while microbat PV+ neurons highly expressed CPLX1, which improves neurotransmission efficiency. Further perturbation of Cplx1 in PV+ neurons impaired ultrasound perception in the mouse auditory cortex. In addition, CPLX1 functioned in other parts of the auditory pathway in microbats but not megabats and exhibited convergent evolution between echolocating microbats and whales. Altogether, we conclude that CPLX1 expression throughout the entire auditory pathway can enhance mammalian ultrasound neurotransmission.

Gastrointestinal Incomplete Degradation Exacerbates Neurotoxic Effects of PLA Microplastics via Oligomer Nanoplastics Formation

Biodegradable plastics, hailed for their environmental friendliness, may pose unforeseen risks as they undergo gastrointestinal degradation, forming oligomer nanoplastics. Despite this, the influence of gastrointestinal degradation on the potential human toxicity of biodegradable plastics remains poorly understood. To this end, the impact of the murine in vivo digestive system is investigated on the biotransformation, biodistribution, and toxicity of PLA polymer and PLA oligomer MPs. Through a 28-day repeated oral gavage study in mice, it is revealed that PLA polymer and oligomer microplastics undergo incomplete and complete degradation, respectively, in the gastrointestinal tract. Incompletely degraded PLA polymer microplastics transform into oligomer nanoplastics, heightening bioavailability and toxicity, thereby exacerbating overall toxic effects. Conversely, complete degradation of PLA oligomer microplastics reduces bioavailability and mitigates toxicity, offering a potential avenue for toxicity reduction. Additionally, the study illuminates shared targets and toxicity mechanisms in Parkinson's disease-like neurotoxicity induced by both PLA polymer and PLA oligomer microplastics. This involves the upregulation of MICU3 in midbrains, leading to neuronal mitochondrial calcium overload. Notably, neurotoxicity is mitigated by inhibiting mitochondrial calcium influx with MCU-i4 or facilitating mitochondrial calcium efflux with DBcAMP in mice. These findings enhance the understanding of the toxicological implications of biodegradable microplastics on human health.

SPRITE: improving spatial gene expression imputation with gene and cell networks

MOTIVATION

Spatially resolved single-cell transcriptomics have provided unprecedented insights into gene expression in situ, particularly in the context of cell interactions or organization of tissues. However, current technologies for profiling spatial gene expression at single-cell resolution are generally limited to the measurement of a small number of genes. To address this limitation, several algorithms have been developed to impute or predict the expression of additional genes that were not present in the measured gene panel. Current algorithms do not leverage the rich spatial and gene relational information in spatial transcriptomics. To improve spatial gene expression predictions, we introduce Spatial Propagation and Reinforcement of Imputed Transcript Expression (SPRITE) as a meta-algorithm that processes predictions obtained from existing methods by propagating information across gene correlation networks and spatial neighborhood graphs.

RESULTS

SPRITE improves spatial gene expression predictions across multiple spatial transcriptomics datasets. Furthermore, SPRITE predicted spatial gene expression leads to improved clustering, visualization, and classification of cells. SPRITE can be used in spatial transcriptomics data analysis to improve inferences based on predicted gene expression.

AVAILABILITY AND IMPLEMENTATION

The SPRITE software package is available at https://github.com/sunericd/SPRITE. Code for generating experiments and analyses in the manuscript is available at https://github.com/sunericd/sprite-figures-and-analyses.

Proteomic features of gray matter layers and superficial white matter of the rhesus monkey neocortex: comparison of prefrontal area 46 and occipital area 17

In this novel large-scale multiplexed immunofluorescence study we comprehensively characterized and compared layer-specific proteomic features within regions of interest of the widely divergent dorsolateral prefrontal cortex (A46) and primary visual cortex (A17) of adult rhesus monkeys. Twenty-eight markers were imaged in rounds of sequential staining, and their spatial distribution precisely quantified within gray matter layers and superficial white matter. Cells were classified as neurons, astrocytes, oligodendrocytes, microglia, or endothelial cells. The distribution of fibers and blood vessels were assessed by quantification of staining intensity across regions of interest. This method revealed multivariate similarities and differences between layers and areas. Protein expression in neurons was the strongest determinant of both laminar and regional differences, whereas protein expression in glia was more important for intra-areal laminar distinctions. Among specific results, we observed a lower glia-to-neuron ratio in A17 than in A46 and the pan-neuronal markers HuD and NeuN were differentially distributed in both brain areas with a lower intensity of NeuN in layers 4 and 5 of A17 compared to A46 and other A17 layers. Astrocytes and oligodendrocytes exhibited distinct marker-specific laminar distributions that differed between regions; notably, there was a high proportion of ALDH1L1-expressing astrocytes and of oligodendrocyte markers in layer 4 of A17. The many nuanced differences in protein expression between layers and regions observed here highlight the need for direct assessment of proteins, in addition to RNA expression, and set the stage for future protein-focused studies of these and other brain regions in normal and pathological conditions.

The Brain-Heart Network of Syncope

Observed and recorded in various forms since ancient times, 'syncope' is often popularly called 'fainting', such that the two terms are used synonymously. Syncope/fainting can be caused by a variety of conditions, including but not limited to head injuries, vertigo, and oxygen deficiency. Here, we draw on a large body of literature on syncope, including the role of a recently discovered set of specialized mammalian neurons. Although the etiology of syncope still remains a mystery, we have attempted to provide a comprehensive account of what is known and what still needs to be performed. Much of our understanding of syncope is owing to studies in the laboratory mouse, whereas evidence from human patients remains scarce. Interestingly, the cardioinhibitory Bezold-Jarisch reflex, recognized in the early 1900s, has an intriguing similarity to-and forms the basis of-syncope. In this review, we have integrated this minimal model into the modern view of the brain-neuron-heart signaling loop of syncope, to which several signaling events contribute. Molecular signaling is our major focus here, presented in terms of a normal heart, and thus, syncope due to abnormal or weak heart activity is not discussed in detail. In addition, we have offered possible directions for clinical intervention based on this model. Overall, this article is expected to generate interest in chronic vertigo and syncope/fainting, an enigmatic condition that affects most humans at some point in life; it is also hoped that this may lead to a mechanism-based clinical intervention in the future.

A single-cell transcriptomic atlas of sensory-dependent gene expression in developing mouse visual cortex

Sensory experience drives the refinement and maturation of neural circuits during postnatal brain development through molecular mechanisms that remain to be fully elucidated. One likely mechanism involves the sensory-dependent expression of genes that encode direct mediators of circuit remodeling within developing cells. However, while studies in adult systems have begun to uncover crucial roles for sensory-induced genes in modifying circuit connectivity, the gene programs induced by brain cells in response to sensory experience during development remain to be fully characterized. Here, we present a single-nucleus RNA-sequencing dataset describing the transcriptional responses of cells in mouse visual cortex to sensory deprivation or sensory stimulation during a developmental window when visual input is necessary for circuit refinement. We sequenced 118,529 individual nuclei across sixteen neuronal and non-neuronal cortical cell types isolated from control, sensory deprived, and sensory stimulated mice, identifying 1,268 unique sensory-induced genes within the developing brain. To demonstrate the utility of this resource, we compared the architecture and ontology of sensory-induced gene programs between cell types, annotated transcriptional induction and repression events based upon RNA velocity, and discovered Neurexin and Neuregulin signaling networks that underlie cell-cell interactions via CellChat . We find that excitatory neurons, especially layer 2/3 pyramidal neurons, are highly sensitive to sensory stimulation, and that the sensory-induced genes in these cells are poised to strengthen synapse-to-nucleus crosstalk by heightening protein serine/threonine kinase activity. Altogether, we expect this dataset to significantly broaden our understanding of the molecular mechanisms through which sensory experience shapes neural circuit wiring in the developing brain.

Microglia-astrocyte crosstalk in the amyloid plaque niche of an Alzheimer's disease mouse model, as revealed by spatial transcriptomics

The amyloid plaque niche is a pivotal hallmark of Alzheimer's disease (AD). Here, we employ two high-resolution spatial transcriptomics (ST) platforms, CosMx and Spatial Enhanced Resolution Omics-sequencing (Stereo-seq), to characterize the transcriptomic alterations, cellular compositions, and signaling perturbations in the amyloid plaque niche in an AD mouse model. We discover heterogeneity in the cellular composition of plaque niches, marked by an increase in microglial accumulation. We profile the transcriptomic alterations of glial cells in the vicinity of plaques and conclude that the microglial response to plaques is consistent across different brain regions, while the astrocytic response is more heterogeneous. Meanwhile, as the microglial density of plaque niches increases, astrocytes acquire a more neurotoxic phenotype and play a key role in inducing GABAergic signaling and decreasing glutamatergic signaling in hippocampal neurons. We thus show that the accumulation of microglia around hippocampal plaques disrupts astrocytic signaling, in turn inducing an imbalance in neuronal synaptic signaling.

Physiological regulation of neuronal Wnt activity is essential for TDP-43 localization and function

Nuclear exclusion of the RNA- and DNA-binding protein TDP-43 can induce neurodegeneration in different diseases. Diverse processes have been implicated to influence TDP-43 mislocalization, including disrupted nucleocytoplasmic transport (NCT); however, the physiological pathways that normally ensure TDP-43 nuclear localization are unclear. The six-transmembrane enzyme glycerophosphodiester phosphodiesterase 2 (GDE2 or GDPD5) cleaves the glycosylphosphatidylinositol (GPI) anchor that tethers some proteins to the membrane. Here we show that GDE2 maintains TDP-43 nuclear localization by regulating the dynamics of canonical Wnt signaling. Ablation of GDE2 causes aberrantly sustained Wnt activation in adult neurons, which is sufficient to cause NCT deficits, nuclear pore abnormalities, and TDP-43 nuclear exclusion. Disruption of GDE2 coincides with TDP-43 abnormalities in postmortem tissue from patients with amyotrophic lateral sclerosis (ALS). Further, GDE2 deficits are evident in human neural cell models of ALS, which display erroneous Wnt activation that, when inhibited, increases mRNA levels of genes regulated by TDP-43. Our study identifies GDE2 as a critical physiological regulator of Wnt signaling in adult neurons and highlights Wnt pathway activation as an unappreciated mechanism contributing to nucleocytoplasmic transport and TDP-43 abnormalities in disease.

Hypnotic treatment improves sleep architecture and EEG disruptions and rescues memory deficits in a mouse model of fragile X syndrome

Fragile X syndrome (FXS) is associated with disrupted cognition and sleep abnormalities. Sleep loss negatively impacts cognitive function, and one untested possibility is that disrupted cognition in FXS is exacerbated by abnormal sleep. We tested whether ML297, a hypnotic acting on G-protein-activated inward-rectifying potassium (GIRK) channels, could reverse sleep phenotypes and disrupted memory in Fmr1-/y mice. Fmr1-/y mice exhibit reduced non-rapid eye movement (NREM) sleep and fragmented NREM architecture, altered sleep electroencephalogram (EEG) oscillations, and reduced EEG coherence between cortical areas; these are partially reversed following ML297 administration. Treatment following contextual fear or spatial learning restores disrupted memory consolidation in Fmr1-/y mice. During memory recall, Fmr1-/y mice show an altered balance of activity among hippocampal principal neurons vs. parvalbumin-expressing interneurons; this is partially reversed by ML297. Because sleep disruption could impact neurophysiological phenotypes in FXS, augmenting sleep may improve disrupted cognition in this disorder.

Neocortical long-range inhibition promotes cortical synchrony and sleep

Behavioral states such as sleep and wake are highly correlated with specific patterns of rhythmic activity in the cortex. During low arousal states such as slow wave sleep, the cortex is synchronized and dominated by low frequency rhythms coordinated across multiple regions. Although recent evidence suggests that GABAergic inhibitory neurons are key players in cortical state modulation, the in vivo circuit mechanisms coordinating synchronized activity among local and distant neocortical networks are not well understood. Here, we show that somatostatin and chondrolectin co-expressing cells (Sst-Chodl cells), a sparse and unique class of neocortical inhibitory neurons, are selectively active during low arousal states and are largely silent during periods of high arousal. In contrast to other neocortical inhibitory neurons, we show these neurons have long-range axons that project across neocortical areas. Activation of Sst-Chodl cells is sufficient to promote synchronized cortical states characteristic of low arousal, with increased spike co-firing and low frequency brain rhythms, and to alter behavioral states by promoting sleep. Contrary to the prevailing belief that sleep is exclusively driven by subcortical mechanisms, our findings reveal that these long-range inhibitory neurons not only track changes in behavioral state but are sufficient to induce both sleep-like cortical states and sleep behavior, establishing a crucial circuit component in regulating behavioral states.

Transcriptome analysis identifies an ASD-Like phenotype in oligodendrocytes and microglia from C58/J amygdala that is dependent on sex and sociability

BACKGROUND

Autism Spectrum Disorder (ASD) is a group of neurodevelopmental disorders with higher incidence in males and is characterized by atypical verbal/nonverbal communication, restricted interests that can be accompanied by repetitive behavior, and disturbances in social behavior. This study investigated brain mechanisms that contribute to sociability deficits and sex differences in an ASD animal model.

METHODS

Sociability was measured in C58/J and C57BL/6J mice using the 3-chamber social choice test. Bulk RNA-Seq and snRNA-Seq identified transcriptional changes in C58/J and C57BL/6J amygdala within which DMRseq was used to measure differentially methylated regions in amygdala.

RESULTS

C58/J mice displayed divergent social strata in the 3-chamber test. Transcriptional and pathway signatures revealed immune-related biological processes differ between C58/J and C57BL/6J amygdala. Hypermethylated and hypomethylated genes were identified in C58/J versus C57BL/6J amygdala. snRNA-Seq data in C58/J amygdala identified differential transcriptional signatures within oligodendrocytes and microglia characterized by increased ASD risk gene expression and predicted impaired myelination that was dependent on sex and sociability. RNA velocity, gene regulatory network, and cell communication analysis showed diminished oligodendrocyte/microglia differentiation. Findings were verified using Bulk RNA-Seq and demonstrated oxytocin's beneficial effects on myelin gene expression.

LIMITATIONS

Our findings are significant. However, limitations can be noted. The cellular mechanisms linking reduced oligodendrocyte differentiation and reduced myelination to an ASD phenotype in C58/J mice need further investigation. Additional snRNA-Seq and spatial studies would determine if effects in oligodendrocytes/microglia are unique to amygdala or if this occurs in other brain regions. Oxytocin's effects need further examination to understand its' potential as an ASD therapeutic.

CONCLUSIONS

Our work demonstrates the C58/J mouse model's utility in evaluating the influence of sex and sociability on the transcriptome in concomitant brain regions involved in ASD. Our single-nucleus transcriptome analysis elucidates potential pathological roles of oligodendrocytes and microglia in ASD. This investigation provides details regarding regulatory features disrupted in these cell types, including transcriptional gene dysregulation, aberrant cell differentiation, altered gene regulatory networks, and changes to key pathways that promote microglia/oligodendrocyte differentiation. Our studies provide insight into interactions between genetic risk and epigenetic processes associated with divergent affiliative behavior and lack of positive sociability.

The transcriptomic and spatial organization of telencephalic GABAergic neuronal types

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.

Social Isolation Induces Changes in the Monoaminergic Signalling in the Rat Medial Prefrontal Cortex

(1) Background: The effects of short-term social isolation during adulthood have not yet been fully established in rats behaviourally, and not at all transcriptomically in the medial prefrontal cortex (mPFC). (2) Methods: We measured the behavioural effects of housing adult male rats in pairs or alone for 10 days. We also used RNA sequencing to measure the accompanying gene expression alterations in the mPFC of male rats. (3) Results: The isolated animals exhibited reduced sociability and social novelty preference, but increased social interaction. There was no change in their aggression, anxiety, or depression-like activity. Transcriptomic analysis revealed a differential expression of 46 genes between the groups. The KEGG pathway analysis showed that differentially expressed genes are involved in neuroactive ligand-receptor interactions, particularly in the dopaminergic and peptidergic systems, and addiction. Subsequent validation confirmed the decreased level of three altered genes: regulator of G protein signalling 9 (Rgs9), serotonin receptor 2c (Htr2c), and Prodynorphin (Pdyn), which are involved in dopaminergic, serotonergic, and peptidergic function, respectively. Antagonizing Htr2c confirmed its role in social novelty discrimination. (4) Conclusions: Social homeostatic regulations include monoaminergic and peptidergic systems of the mPFC.

The ENCODE mouse postnatal developmental time course identifies regulatory programs of cell types and cell states

Postnatal genomic regulation significantly influences tissue and organ maturation but is under-studied relative to existing genomic catalogs of adult tissues or prenatal development in mouse. The ENCODE4 consortium generated the first comprehensive single-nucleus resource of postnatal regulatory events across a diverse set of mouse tissues. The collection spans seven postnatal time points, mirroring human development from childhood to adulthood, and encompasses five core tissues. We identified 30 cell types, further subdivided into 69 subtypes and cell states across adrenal gland, left cerebral cortex, hippocampus, heart, and gastrocnemius muscle. Our annotations cover both known and novel cell differentiation dynamics ranging from early hippocampal neurogenesis to a new sex-specific adrenal gland population during puberty. We used an ensemble Latent Dirichlet Allocation strategy with a curated vocabulary of 2,701 regulatory genes to identify regulatory "topics," each of which is a gene vector, linked to cell type differentiation, subtype specialization, and transitions between cell states. We find recurrent regulatory topics in tissue-resident macrophages, neural cell types, endothelial cells across multiple tissues, and cycling cells of the adrenal gland and heart. Cell-type-specific topics are enriched in transcription factors and microRNA host genes, while chromatin regulators dominate mitosis topics. Corresponding chromatin accessibility data reveal dynamic and sex-specific regulatory elements, with enriched motifs matching transcription factors in regulatory topics. Together, these analyses identify both tissue-specific and common regulatory programs in postnatal development across multiple tissues through the lens of the factors regulating transcription.

SIMS: A deep-learning label transfer tool for single-cell RNA sequencing analysis

Cell atlases serve as vital references for automating cell labeling in new samples, yet existing classification algorithms struggle with accuracy. Here we introduce SIMS (scalable, interpretable machine learning for single cell), a low-code data-efficient pipeline for single-cell RNA classification. We benchmark SIMS against datasets from different tissues and species. We demonstrate SIMS's efficacy in classifying cells in the brain, achieving high accuracy even with small training sets (<3,500 cells) and across different samples. SIMS accurately predicts neuronal subtypes in the developing brain, shedding light on genetic changes during neuronal differentiation and postmitotic fate refinement. Finally, we apply SIMS to single-cell RNA datasets of cortical organoids to predict cell identities and uncover genetic variations between cell lines. SIMS identifies cell-line differences and misannotated cell lineages in human cortical organoids derived from different pluripotent stem cell lines. Altogether, we show that SIMS is a versatile and robust tool for cell-type classification from single-cell datasets.

DOT: a flexible multi-objective optimization framework for transferring features across single-cell and spatial omics

Single-cell transcriptomics and spatially-resolved imaging/sequencing technologies have revolutionized biomedical research. However, they suffer from lack of spatial information and a trade-off of resolution and gene coverage, respectively. We propose DOT, a multi-objective optimization framework for transferring cellular features across these data modalities, thus integrating their complementary information. DOT uses genes beyond those common to the data modalities, exploits the local spatial context, transfers spatial features beyond cell-type information, and infers absolute/relative abundance of cell populations at tissue locations. Thus, DOT bridges single-cell transcriptomics data with both high- and low-resolution spatially-resolved data. Moreover, DOT combines practical aspects related to cell composition, heterogeneity, technical effects, and integration of prior knowledge. Our fast implementation based on the Frank-Wolfe algorithm achieves state-of-the-art or improved performance in localizing cell features in high- and low-resolution spatial data and estimating the expression of unmeasured genes in low-coverage spatial data.

A suite of enhancer AAVs and transgenic mouse lines for genetic access to cortical cell types

The mammalian cortex is comprised of cells with different morphological, physiological, and molecular properties that can be classified according to shared properties into cell types. Defining the contribution of each cell type to the computational and cognitive processes that are guided by the cortex is essential for understanding its function in health and disease. We use transcriptomic and epigenomic cortical cell type taxonomies from mice and humans to define marker genes and enhancers, and to build genetic tools for cortical cell types. Here, we present a large toolkit for selective targeting of cortical populations, including mouse transgenic lines and recombinant adeno-associated virus (AAV) vectors containing genomic enhancers. We report evaluation of fifteen new transgenic driver lines and over 680 different enhancer AAVs covering all major subclasses of cortical cells, with many achieving a high degree of specificity, comparable with existing transgenic lines. We find that the transgenic lines based on marker genes can provide exceptional specificity and completeness of cell type labeling, but frequently require generation of a triple-transgenic cross for best usability/specificity. On the other hand, enhancer AAVs are easy to screen and use, and can be easily modified to express diverse cargo, such as recombinases. However, their use depends on many factors, such as viral titer and route of administration. The tools reported here as well as the scaled process of tool creation provide an unprecedented resource that should enable diverse experimental strategies towards understanding mammalian cortex and brain function.

Deciphering the genetic code of neuronal type connectivity through bilinear modeling

Understanding how different neuronal types connect and communicate is critical to interpreting brain function and behavior. However, it has remained a formidable challenge to decipher the genetic underpinnings that dictate the specific connections formed between neuronal types. To address this, we propose a novel bilinear modeling approach that leverages the architecture similar to that of recommendation systems. Our model transforms the gene expressions of presynaptic and postsynaptic neuronal types, obtained from single-cell transcriptomics, into a covariance matrix. The objective is to construct this covariance matrix that closely mirrors a connectivity matrix, derived from connectomic data, reflecting the known anatomical connections between these neuronal types. When tested on a dataset of Caenorhabditis elegans, our model achieved a performance comparable to, if slightly better than, the previously proposed spatial connectome model (SCM) in reconstructing electrical synaptic connectivity based on gene expressions. Through a comparative analysis, our model not only captured all genetic interactions identified by the SCM but also inferred additional ones. Applied to a mouse retinal neuronal dataset, the bilinear model successfully recapitulated recognized connectivity motifs between bipolar cells and retinal ganglion cells, and provided interpretable insights into genetic interactions shaping the connectivity. Specifically, it identified unique genetic signatures associated with different connectivity motifs, including genes important to cell-cell adhesion and synapse formation, highlighting their role in orchestrating specific synaptic connections between these neurons. Our work establishes an innovative computational strategy for decoding the genetic programming of neuronal type connectivity. It not only sets a new benchmark for single-cell transcriptomic analysis of synaptic connections but also paves the way for mechanistic studies of neural circuit assembly and genetic manipulation of circuit wiring.

Cell-type-specific cholinergic control of granular retrosplenial cortex with implications for angular velocity coding across brain states

Cholinergic receptor activation enables the persistent firing of cortical pyramidal neurons, providing a key cellular basis for theories of spatial navigation involving working memory, path integration, and head direction encoding. The granular retrosplenial cortex (RSG) is important for spatially-guided behaviors, but how acetylcholine impacts RSG neurons is unknown. Here, we show that a transcriptomically, morphologically, and biophysically distinct RSG cell-type - the low-rheobase (LR) neuron - has a very distinct expression profile of cholinergic muscarinic receptors compared to all other neighboring excitatory neuronal subtypes. LR neurons do not fire persistently in response to cholinergic agonists, in stark contrast to all other principal neuronal subtypes examined within the RSG and across midline cortex. This lack of persistence allows LR neuron models to rapidly compute angular head velocity (AHV), independent of cholinergic changes seen during navigation. Thus, LR neurons can consistently compute AHV across brain states, highlighting the specialized RSG neural codes supporting navigation.

Haploinsufficiency underlies the neurodevelopmental consequences of SLC6A1 variants

Heterozygous variants in SLC6A1, encoding the GAT-1 GABA transporter, are associated with seizures, developmental delay, and autism. The majority of affected individuals carry missense variants, many of which are recurrent germline de novo mutations, raising the possibility of gain-of-function or dominant-negative effects. To understand the functional consequences, we performed an in vitro GABA uptake assay for 213 unique variants, including 24 control variants. De novo variants consistently resulted in a decrease in GABA uptake, in keeping with haploinsufficiency underlying all neurodevelopmental phenotypes. Where present, ClinVar pathogenicity reports correlated well with GABA uptake data; the functional data can inform future reports for the remaining 72% of unscored variants. Surface localization was assessed for 86 variants; two-thirds of loss-of-function missense variants prevented GAT-1 from being present on the membrane while GAT-1 was on the surface but with reduced activity for the remaining third. Surprisingly, recurrent de novo missense variants showed moderate loss-of-function effects that reduced GABA uptake with no evidence for dominant-negative or gain-of-function effects. Using linear regression across multiple missense severity scores to extrapolate the functional data to all potential SLC6A1 missense variants, we observe an abundance of GAT-1 residues that are sensitive to substitution. The extent of this missense vulnerability accounts for the clinically observed missense enrichment; overlap with hypermutable CpG sites accounts for the recurrent missense variants. Strategies to increase the expression of the wild-type SLC6A1 allele are likely to be beneficial across neurodevelopmental disorders, though the developmental stage and extent of required rescue remain unknown.

Neuronal cell populations in circumoral nerve ring of sea cucumber Apostichopus japonicus: Ultrastructure and transcriptional profile

The echinoderm nervous system has been studied as a model for understanding the evolution of the chordate nervous system. Neuronal cells are essential groups that release a 'cocktail' of messenger molecules providing a spectrum of biological actions in the nervous system. Among echinoderms, most evidence on neuronal cell types has been obtained from starfish and sea urchin. In sea cucumbers, most research has focused on the location of neuronal cells, whereas their transcriptional features have rarely been investigated. Here, we observed the ultrastructure of neuronal cells in the sea cucumber, Apostichopus japonicus. The transcriptional profile of neuronal cells from the circumoral nerve ring (CNR) was investigated using single-cell RNA sequencing (scRNA-seq), and a total of six neuronal cell types were identified. 26 neuropeptide precursor genes (NPPs) and 28 G-protein-coupled receptors (GPCR) were expressed in the six neuronal cell types, comprising five NPP/NP-GPCR pairs. Unsupervised pseudotime analysis of neuronal cells showed their different differentiation status. We also located the neuronal cells in the CNR by immunofluorescence (IF) and identified the potential hub genes of key cell populations. This broad resource serves as a valuable support in the development of cell-specific markers for accurate cell-type identification in sea cucumbers. It also contributes to facilitating comparison across species, providing a deeper understanding of the evolutionary processes of neuronal cells.

Repairing the in situ hybridization missing data in the hippocampus region by using a 3D residual U-Net model

The hippocampus is a critical brain region. Transcriptome data provides valuable insights into the structure and function of the hippocampus at the gene level. However, transcriptome data is often incomplete. To address this issue, we use the convolutional neural network model to repair the missing voxels in the hippocampus region, based on Allen institute coronal slices in situ hybridization (ISH) dataset. Moreover, we analyze the gene expression correlation between coronal and sagittal dataset in the hippocampus region. The results demonstrated that the trend of gene expression correlation between the coronal and sagittal datasets remained consistent following the repair of missing data in the coronal ISH dataset. In the last, we use repaired ISH dataset to identify novel genes specific to hippocampal subregions. Our findings demonstrate the accuracy and effectiveness of using deep learning method to repair ISH missing data. After being repaired, ISH has the potential to improve our comprehension of the hippocampus's structure and function.

Single-cell investigation of lead toxicity from neurodevelopment to neurodegeneration: Current review and future opportunities

Human exposure to the metal lead (Pb) is prevalent and associated with adverse neurodevelopmental and neurodegenerative outcomes. Pb disrupts normal brain function by inducing oxidative stress and neuroinflammation, altering cellular metabolism, and displacing essential metals. Prior studies on the molecular impacts of Pb have examined bulk tissues, which collapse information across all cell types, or in targeted cells, which are limited to cell autonomous effects. These approaches are unable to represent the complete biological implications of Pb exposure because the brain is a cooperative network of highly heterogeneous cells, with cellular diversity and proportions shifting throughout development, by brain region, and with disease. New technologies are necessary to investigate whether Pb and other environmental exposures alter cell composition in the brain and whether they cause molecular changes in a cell-type-specific manner. Cutting-edge, single-cell approaches now enable research resolving cell-type-specific effects from bulk tissues. This article reviews existing Pb neurotoxicology studies with genome-wide molecular signatures and provides a path forward for the field to implement single-cell approaches with practical recommendations.

Compensation between FOXP transcription factors maintains proper striatal function

Spiny projection neurons (SPNs) of the striatum are critical in integrating neurochemical information to coordinate motor and reward-based behavior. Mutations in the regulatory transcription factors expressed in SPNs can result in neurodevelopmental disorders (NDDs). Paralogous transcription factors Foxp1 and Foxp2, which are both expressed in the dopamine receptor 1 (D1) expressing SPNs, are known to have variants implicated in NDDs. Utilizing mice with a D1-SPN-specific loss of Foxp1, Foxp2, or both and a combination of behavior, electrophysiology, and cell-type-specific genomic analysis, loss of both genes results in impaired motor and social behavior as well as increased firing of the D1-SPNs. Differential gene expression analysis implicates genes involved in autism risk, electrophysiological properties, and neuronal development and function. Viral-mediated re-expression of Foxp1 into the double knockouts is sufficient to restore electrophysiological and behavioral deficits. These data indicate complementary roles between Foxp1 and Foxp2 in the D1-SPNs.

GHSR in a Subset of GABA Neurons Controls Food Deprivation-Induced Hyperphagia in Male Mice

The growth hormone secretagogue receptor (GHSR), primarily known as the receptor for the hunger hormone ghrelin, potently controls food intake, yet the specific Ghsr-expressing cells mediating the orexigenic effects of this receptor remain incompletely characterized. Since Ghsr is expressed in gamma-aminobutyric acid (GABA)-producing neurons, we sought to investigate whether the selective expression of Ghsr in a subset of GABA neurons is sufficient to mediate GHSR's effects on feeding. First, we crossed mice that express a tamoxifen-dependent Cre recombinase in the subset of GABA neurons that express glutamic acid decarboxylase 2 (Gad2) enzyme (Gad2-CreER mice) with reporter mice, and found that ghrelin mainly targets a subset of Gad2-expressing neurons located in the hypothalamic arcuate nucleus (ARH) and that is predominantly segregated from Agouti-related protein (AgRP)-expressing neurons. Analysis of various single-cell RNA-sequencing datasets further corroborated that the primary subset of cells coexpressing Gad2 and Ghsr in the mouse brain are non-AgRP ARH neurons. Next, we crossed Gad2-CreER mice with reactivable GHSR-deficient mice to generate mice expressing Ghsr only in Gad2-expressing neurons (Gad2-GHSR mice). We found that ghrelin treatment induced the expression of the marker of transcriptional activation c-Fos in the ARH of Gad2-GHSR mice, yet failed to induce food intake. In contrast, food deprivation-induced refeeding was higher in Gad2-GHSR mice than in GHSR-deficient mice and similar to wild-type mice, suggesting that ghrelin-independent roles of GHSR in a subset of GABA neurons is sufficient for eliciting full compensatory hyperphagia in mice.

Organization of corticocortical and thalamocortical top-down inputs in the primary visual cortex

Unified visual perception requires integration of bottom-up and top-down inputs in the primary visual cortex (V1), yet the organization of top-down inputs in V1 remains unclear. Here, we used optogenetics-assisted circuit mapping to identify how multiple top-down inputs from higher-order cortical and thalamic areas engage V1 excitatory and inhibitory neurons. Top-down inputs overlap in superficial layers yet segregate in deep layers. Inputs from the medial secondary visual cortex (V2M) and anterior cingulate cortex (ACA) converge on L6 Pyrs, whereas ventrolateral orbitofrontal cortex (ORBvl) and lateral posterior thalamic nucleus (LP) inputs are processed in parallel in Pyr-type-specific subnetworks (Pyr←ORBvl and Pyr←LP) and drive mutual inhibition between them via local interneurons. Our study deepens understanding of the top-down modulation mechanisms of visual processing and establishes that V2M and ACA inputs in L6 employ integrated processing distinct from the parallel processing of LP and ORBvl inputs in L5.

A single-nucleus transcriptomic atlas of medium spiny neurons in the rat nucleus accumbens

Neural processing of rewarding stimuli involves several distinct regions, including the nucleus accumbens (NAc). The majority of NAc neurons are GABAergic projection neurons known as medium spiny neurons (MSNs). MSNs are broadly defined by dopamine receptor expression, but evidence suggests that a wider array of subtypes exist. To study MSN heterogeneity, we analyzed single-nucleus RNA sequencing data from the largest available rat NAc dataset. Analysis of 48,040 NAc MSN nuclei identified major populations belonging to the striosome and matrix compartments. Integration with mouse and human data indicated consistency across species and disease-relevance scoring using genome-wide association study results revealed potentially differential roles for MSN populations in substance use disorders. Additional high-resolution clustering identified 34 transcriptomically distinct subtypes of MSNs definable by a limited number of marker genes. Together, these data demonstrate the diversity of MSNs in the NAc and provide a basis for more targeted genetic manipulation of specific populations.