Genome-wide atlas of gene expression in the adult mouse brain

Claustrum and dorsal endopiriform cortex complex cell-identity is determined by Nurr1 and regulates hallucinogenic-like states in mice

The Claustrum/dorsal endopiriform cortex complex (CLA) is an enigmatic brain region with extensive glutamatergic projections to multiple cortical areas. The transcription factor Nurr1 is highly expressed in the CLA, but its role in this region is not understood. By using conditional gene-targeted mice, we show that Nurr1 is a crucial regulator of CLA neuron identity. Although CLA neurons remain intact in the absence of Nurr1, the distinctive gene expression pattern in the CLA is abolished. CLA has been hypothesized to control hallucinations, but little is known of how the CLA responds to hallucinogens. After the deletion of Nurr1 in the CLA, both hallucinogen receptor expression and signaling are lost. Furthermore, functional ultrasound and Neuropixel electrophysiological recordings revealed that the hallucinogenic-receptor agonists' effects on functional connectivity between prefrontal and sensorimotor cortices are altered in Nurr1-ablated mice. Our findings suggest that Nurr1-targeted strategies provide additional avenues for functional studies of the CLA.

Concerning neuromodulation as treatment of neurological and neuropsychiatric disorder: Insights gained from selective targeting of the subthalamic nucleus, para-subthalamic nucleus and zona incerta in rodents

Neuromodulation such as deep brain stimulation (DBS) is advancing as a clinical intervention in several neurological and neuropsychiatric disorders, including Parkinson's disease, dystonia, tremor, and obsessive-compulsive disorder (OCD) for which DBS is already applied to alleviate severely afflicted individuals of symptoms. Tourette syndrome and drug addiction are two additional disorders for which DBS is in trial or proposed as treatment. However, some major remaining obstacles prevent this intervention from reaching its full therapeutic potential. Side-effects have been reported, and not all DBS-treated individuals are relieved of their symptoms. One major target area for DBS electrodes is the subthalamic nucleus (STN) which plays important roles in motor, affective and associative functions, with impact on for example movement, motivation, impulsivity, compulsivity, as well as both reward and aversion. The multifunctionality of the STN is complex. Decoding the anatomical-functional organization of the STN could enhance strategic targeting in human patients. The STN is located in close proximity to zona incerta (ZI) and the para-subthalamic nucleus (pSTN). Together, the STN, pSTN and ZI form a highly heterogeneous and clinically important brain area. Rodent-based experimental studies, including opto- and chemogenetics as well as viral-genetic tract tracings, provide unique insight into complex neuronal circuitries and their impact on behavior with high spatial and temporal precision. This research field has advanced tremendously over the past few years. Here, we provide an inclusive review of current literature in the pre-clinical research fields centered around STN, pSTN and ZI in laboratory mice and rats; the three highly heterogeneous and enigmatic structures brought together in the context of relevance for treatment strategies. Specific emphasis is placed on methods of manipulation and behavioral impact.

A graph self-supervised residual learning framework for domain identification and data integration of spatial transcriptomics

Spatial transcriptomics (ST) technologies allow for comprehensive characterization of gene expression patterns in the context of tissue microenvironment. However, accurately identifying domains with spatial coherence in both gene expression and histology in situ and effectively integrating data from multi-sample remains challenging. Here, we propose ResST, a graph self-supervised residual learning model based on graph neural network and Margin Disparity Discrepancy (MDD) theory. ResST aggregates gene expression, biological effects, spatial location, and morphological information to capture nonlinear relationships between a cell and surrounding cells for spatial domain identification. Also, ResST integrates multiple ST datasets and aligns latent embeddings based on MDD theory for correcting batch effects. Results show that ResST identifies continuous spatial domains at a finer scale in ten ST datasets acquired with different technologies. Moreover, ResST efficiently integrated data from multiple tissue sections vertically or horizontally while correcting batch effects. Overall, ResST demonstrates exceptional performance in analyzing ST datasets.

Unsupervised pattern identification in spatial gene expression atlas reveals mouse brain regions beyond established ontology

The rapid growth of large-scale spatial gene expression data demands efficient and reliable computational tools to extract major trends of gene expression in their native spatial context. Here, we used stability-driven unsupervised learning (i.e., staNMF) to identify principal patterns (PPs) of 3D gene expression profiles and understand spatial gene distribution and anatomical localization at the whole mouse brain level. Our subsequent spatial correlation analysis systematically compared the PPs to known anatomical regions and ontology from the Allen Mouse Brain Atlas using spatial neighborhoods. We demonstrate that our stable and spatially coherent PPs, whose linear combinations accurately approximate the spatial gene data, are highly correlated with combinations of expert-annotated brain regions. These PPs yield a brain ontology based purely on spatial gene expression. Our PP identification approach outperforms principal component analysis and typical clustering algorithms on the same task. Moreover, we show that the stable PPs reveal marked regional imbalance of brainwide genetic architecture, leading to region-specific marker genes and gene coexpression networks. Our findings highlight the advantages of stability-driven machine learning for plausible biological discovery from dense spatial gene expression data, streamlining tasks that are infeasible by conventional manual approaches.

Ultrasound pulse repetition frequency preferentially activates different neuron populations independent of cell type

Objective. Transcranial ultrasound (US) stimulation serves as an external input to a neuron, and thus the evoked response relies on neurons' intrinsic properties. Neural activity is limited to a couple hundred hertz and often exhibits preference to input frequencies. Accordingly, US pulsed at specific physiologic pulse repetition frequencies (PRFs) may selectively engage neurons with the corresponding input frequency preference. However, most US parametric studies examine the effects of supraphysiologic PRFs. It remains unclear whether pulsing US at different physiologic PRFs could activate distinct neurons in the awake mammalian brain.Approach. We recorded cellular calcium responses of individual motor cortex neurons to US pulsed at PRFs of 10, 40, and 140 Hz in awake mice. We compared the evoked responses across these PRFs in the same neurons. To further understand the cell-type dependent effects, we categorized the recorded neurons as parvalbumin positive fast spiking interneurons or putative excitatory neurons and analyzed single-cell mechanosensitive channel expression in mice and humans using the Allen Brain Institute's RNA-sequencing databases.Main results. We discovered that many neurons were preferentially activated by only one PRF and different PRFs selectively engaged distinct neuronal populations. US-evoked cellular calcium responses exhibited the same characteristics as those naturally occurring during spiking, suggesting that US increases intrinsic neuronal activity. Furthermore, evoked responses were similar between fast-spiking inhibitory neurons and putative excitatory neurons. Thus, variation in individual neuron's cellular properties dominates US-evoked response heterogeneity, consistent with our observed cell-type independent expression patterns of mechanosensitive channels across individual neurons in mice and humans. Finally, US transiently increased network synchrony without producing prolonged over-synchronization that could be detrimental to neural circuit functions.Significance. These results highlight the feasibility of activating distinct neuronal subgroups by varying PRF and the potential to improve neuromodulation effects by combining physiologic PRFs.

Contributions of mirror-image hair cell orientation to mouse otolith organ and zebrafish neuromast function

Otolith organs in the inner ear and neuromasts in the fish lateral-line harbor two populations of hair cells oriented to detect stimuli in opposing directions. The underlying mechanism is highly conserved: the transcription factor EMX2 is regionally expressed in just one hair cell population and acts through the receptor GPR156 to reverse cell orientation relative to the other population. In mouse and zebrafish, loss of Emx2 results in sensory organs that harbor only one hair cell orientation and are not innervated properly. In zebrafish, Emx2 also confers hair cells with reduced mechanosensory properties. Here, we leverage mouse and zebrafish models lacking GPR156 to determine how detecting stimuli of opposing directions serves vestibular function, and whether GPR156 has other roles besides orienting hair cells. We find that otolith organs in Gpr156 mouse mutants have normal zonal organization and normal type I-II hair cell distribution and mechano-electrical transduction properties. In contrast, gpr156 zebrafish mutants lack the smaller mechanically-evoked signals that characterize Emx2-positive hair cells. Loss of GPR156 does not affect orientation-selectivity of afferents in mouse utricle or zebrafish neuromasts. Consistent with normal otolith organ anatomy and afferent selectivity, Gpr156 mutant mice do not show overt vestibular dysfunction. Instead, performance on two tests that engage otolith organs is significantly altered - swimming and off-vertical-axis rotation. We conclude that GPR156 relays hair cell orientation and transduction information downstream of EMX2, but not selectivity for direction-specific afferents. These results clarify how molecular mechanisms that confer bi-directionality to sensory organs contribute to function, from single hair cell physiology to animal behavior.

Brain morphological analysis in mice with hyperactivation of the hedgehog signaling pathway

Hedgehog signaling is a highly conserved pathway that plays pivotal roles in morphogenesis, tumorigenesis, osteogenesis, and wound healing. Previous investigations in patients with Gorlin syndrome found low harm avoidance traits, and increased volumes in the cerebrum, cerebellum, and cerebral ventricles, suggesting the association between brain morphology and the constitutive hyperactivation of hedgehog signaling, while the changes of regional brain volumes in upregulated hedgehog signaling pathway remains unclear so far. Herein, we investigated comprehensive brain regional volumes using quantitative structural brain MRI, and identified increased volumes of amygdala, striatum, and pallidum on the global segmentation, and increased volumes of the lateral and medial parts of the central nucleus of the amygdala on the detail segmentation in Ptch heterozygous deletion mice. Our data may enhance comprehension of the association between brain morphogenic changes and hyperactivity in hedgehog signaling.

Sex chromosomes and hormones independently influence healthy brain development but act similarly after cranial radiation

The course of normal development and response to pathology are strongly influenced by biological sex. For instance, female childhood cancer survivors who have undergone cranial radiation therapy (CRT) tend to display more pronounced cognitive deficits than their male counterparts. Sex effects can be the result of sex chromosome complement (XX vs. XY) and/or gonadal hormone influence. The contributions of each can be separated using the four-core genotype mouse model (FCG), where sex chromosome complement and gonadal sex are decoupled. While studies of FCG mice have evaluated brain differences in adulthood, it is still unclear how sex chromosome and sex hormone effects emerge through development in both healthy and pathological contexts. Our study utilizes longitudinal MRI with the FCG model to investigate sex effects in healthy development and after CRT in wildtype and immune-modified Ccl2-knockout mice. Our findings in normally developing mice reveal a relatively prominent chromosome effect prepubertally, compared to sex hormone effects which largely emerge later. Spatially, sex chromosome and hormone influences were independent of one another. After CRT in Ccl2-knockout mice, both male chromosomes and male hormones similarly improved brain outcomes but did so more separately than in combination. Our findings highlight the crucial role of sex chromosomes in early development and identify roles for sex chromosomes and hormones after CRT-induced inflammation, highlighting the influences of biological sex in both normal brain development and pathology.

Systematic comparison of sequencing-based spatial transcriptomic methods

Recent developments of sequencing-based spatial transcriptomics (sST) have catalyzed important advancements by facilitating transcriptome-scale spatial gene expression measurement. Despite this progress, efforts to comprehensively benchmark different platforms are currently lacking. The extant variability across technologies and datasets poses challenges in formulating standardized evaluation metrics. In this study, we established a collection of reference tissues and regions characterized by well-defined histological architectures, and used them to generate data to compare 11 sST methods. We highlighted molecular diffusion as a variable parameter across different methods and tissues, significantly affecting the effective resolutions. Furthermore, we observed that spatial transcriptomic data demonstrate unique attributes beyond merely adding a spatial axis to single-cell data, including an enhanced ability to capture patterned rare cell states along with specific markers, albeit being influenced by multiple factors including sequencing depth and resolution. Our study assists biologists in sST platform selection, and helps foster a consensus on evaluation standards and establish a framework for future benchmarking efforts that can be used as a gold standard for the development and benchmarking of computational tools for spatial transcriptomic analysis.

Sustained fentanyl exposure inhibits neuronal activity in dissociated striatal neuronal-glial cocultures through actions independent of opioid receptors

Besides having high potency and efficacy at the µ-opioid (MOR) and other opioid receptor types, fentanyl has some affinity for some adrenergic receptor types, which may underlie its unique pathophysiological differences from typical opioids. To better understand the unique actions of fentanyl, we assessed the extent to which fentanyl alters striatal medium spiny neuron (MSN) activity via opioid receptors or α1-adrenoceptors in dopamine type 1 or type 2 receptor (D1 or D2)-expressing MSNs. In neuronal and mixed-glial cocultures from the striatum, acute fentanyl (100 nM) exposure decreased the frequency of spontaneous action potentials. Overnight exposure of cocultures to 100 nM fentanyl severely reduced the proportion of MSNs with spontaneous action potentials, which was unaffected by coexposure to the opioid receptor antagonist naloxone (10 µM) but fully negated by coadministering the pan-α1-adrenoceptor inverse agonist prazosin (100 nM) and partially reversed by the selective α1A-adrenoceptor antagonist RS 100329 (300 nM). Acute fentanyl (100 nM) exposure modestly reduced the frequency of action potentials and caused firing rate adaptations in D2, but not D1, MSNs. Prolonged (2-5 h) fentanyl (100 nM) application dramatically attenuated firing rates in both D1 and D2 MSNs. To identify possible cellular sites of α1-adrenoceptor action, α1-adrenoceptors were localized in subpopulations of striatal astroglia and neurons by immunocytochemistry and Adra1a mRNA by in situ hybridization in astrocytes. Thus, sustained fentanyl exposure can inhibit striatal MSN activity via a nonopioid receptor-dependent pathway, which may be modulated via complex actions in α1-adrenoceptor-expressing striatal neurons and/or glia.NEW & NOTEWORTHY Acute fentanyl exposure attenuated the activity of striatal medium spiny neurons (MSNs) in vitro and in dopamine D2, but not D1, receptor-expressing MSNs in ex vivo slices. By contrast, sustained fentanyl exposure suppressed the spontaneous activity of MSNs cocultured with glia through a nonopioid receptor-dependent mechanism modulated, in part, by α1-adrenoceptors. Fentanyl exposure can affect striatal function via a nonopioid receptor mechanism of action that appears mediated by α1-adrenoreceptor-expressing striatal neurons and/or astroglia.

Association of medullary reticular formation ventral part with spasticity in mice suffering from photothrombotic stroke

Strokes cause spasticity via stretch reflex hyperexcitability in the spinal cord, and spastic paralysis due to involuntary muscle contraction in the hands and fingers can severely restrict skilled hand movements. However, the underlying neurological mechanisms remain unknown. Using a mouse model of spasticity after stroke, we demonstrate changes in neuronal activity with and without electrostimulation of the afferent nerve to induce the stretch reflex, measured using quantitative activation-induced manganese-enhanced magnetic resonance imaging. Neuronal activity increased within the ventral medullary reticular formation (MdV) in the contralesional brainstem during the acute post-stroke phase, and this increase was characterised by activation of circuits involved in spasticity. Interestingly, ascending electrostimulation inhibited the MdV activity on the stimulation side in normal conditions. Moreover, immunohistochemical staining showed that, in the acute phase, the density of GluA1, one of the α-amino-3 hydroxy‑5 methyl -4 isoxazolepropionic acid receptor (AMPAR) subunits, at the synapses of MdV neurons was significantly increased. In addition, the GluA1/GluA2 ratio in these receptors was altered at 2 weeks post-stroke, confirming homeostatic plasticity as the underlying mechanisms of spasticity. These results provide new insights into the relationship between impaired skilled movements and spasticity at the acute post-stroke phase.

Interaction between the TBC1D24 TLDc domain and the KIBRA C2 domain is disrupted by two epilepsy-associated TBC1D24 missense variants

Mutations of human TBC1D24 are associated with deafness, epilepsy, or DOORS syndrome (deafness, onychodystrophy, osteodystrophy, cognitive disability, and seizures). The causal relationships between TBC1D24 variants and the different clinical phenotypes are not understood. Our hypothesis is that phenotypic heterogeneity of missense mutations of TBC1D24 results, in part, from perturbed binding of different protein partners. To discover novel protein partners of TBC1D24, we conducted yeast two-hybrid (Y2H) screen using mouse full-length TBC1D24 as bait. Kidney and brain protein (KIBRA), a scaffold protein encoded by Wwc1, was identified as a partner of TBC1D24. KIBRA functions in the Hippo signaling pathway and is important for human cognition and memory. The TBC1D24 TLDc domain binds to KIBRA full-length and to its C2 domain, confirmed by Y2H assays. No interaction was detected with Y2H assays between the KIBRA C2 domain and TLDc domains of NCOA7, MEAK7, and OXR1. Moreover, the C2 domains of other WWC family proteins do not interact with the TLDc domain of TBC1D24, demonstrating specificity. The mRNAs encoding TBC1D24 and KIBRA proteins in mouse are coexpressed at least in a subset of hippocampal cells indicating availability to interact in vivo. As two epilepsy-associated recessive variants (Gly511Arg and Ala515Val) in the TLDc domain of human TBC1D24 disrupt the interaction with the human KIBRA C2 domain, this study reveals a pathogenic mechanism of TBC1D24-associated epilepsy, linking the TBC1D24 and KIBRA pathways. The interaction of TBC1D24-KIBRA is physiologically meaningful and necessary to reduce the risk of epilepsy.

STIE: Single-cell level deconvolution, convolution, and clustering in in situ capturing-based spatial transcriptomics

In in situ capturing-based spatial transcriptomics, spots of the same size and printed at fixed locations cannot precisely capture the randomly-located single cells, therefore inherently failing to profile transcriptome at the single-cell level. To this end, we present STIE, an Expectation Maximization algorithm that aligns the spatial transcriptome to its matched histology image-based nuclear morphology and recovers missing cells from ~70% gap area, thereby achieving the real single-cell level and whole-slide scale deconvolution, convolution, and clustering for both low- and high-resolution spots. STIE characterizes cell-type-specific gene expression and demonstrates outperforming concordance with true cell-type-specific transcriptomic signatures than the other spot- and subspot-level methods. Furthermore, STIE reveals the single-cell level insights, for instance, lower actual spot resolution than its reported spot size, unbiased evaluation of cell type colocalization, superior power of high-resolution spot in distinguishing nuanced cell types, and spatial cell-cell interactions at the single-cell level other than spot level.

The short-term plasticity of VIP interneurons in motor cortex

Short-term plasticity is an important feature in the brain for shaping neural dynamics and for information processing. Short-term plasticity is known to depend on many factors including brain region, cortical layer, and cell type. Here we focus on vasoactive-intestinal peptide (VIP) interneurons (INs). VIP INs play a key disinhibitory role in cortical circuits by inhibiting other IN types, including Martinotti cells (MCs) and basket cells (BCs). Despite this prominent role, short-term plasticity at synapses to and from VIP INs is not well described. In this study, we therefore characterized the short-term plasticity at inputs and outputs of genetically targeted VIP INs in mouse motor cortex. To explore inhibitory to inhibitory (I → I) short-term plasticity at layer 2/3 (L2/3) VIP IN outputs onto L5 MCs and BCs, we relied on a combination of whole-cell recording, 2-photon microscopy, and optogenetics, which revealed that VIP IN→MC/BC synapses were consistently short-term depressing. To explore excitatory (E) → I short-term plasticity at inputs to VIP INs, we used extracellular stimulation. Surprisingly, unlike VIP IN outputs, E → VIP IN synapses exhibited heterogeneous short-term dynamics, which we attributed to the target VIP IN cell rather than the input. Computational modeling furthermore linked the diversity in short-term dynamics at VIP IN inputs to a wide variability in probability of release. Taken together, our findings highlight how short-term plasticity at VIP IN inputs and outputs is specific to synapse type. We propose that the broad diversity in short-term plasticity of VIP IN inputs forms a basis to code for a broad range of contrasting signal dynamics.

The influence of extended fasting on thyroid hormone: local and differentiated regulatory mechanisms

The hypometabolism induced by fasting has great potential in maintaining health and improving survival in extreme environments, among which thyroid hormone (TH) plays an important role in the adaptation and the formation of new energy metabolism homeostasis during long-term fasting. In the present review, we emphasize the potential of long-term fasting to improve physical health and emergency rescue in extreme environments, introduce the concept and pattern of fasting and its impact on the body's energy metabolism consumption. Prolonged fasting has more application potential in emergency rescue in special environments. The changes of THs caused by fasting, including serum biochemical characteristics, responsiveness of the peripheral and central hypothalamus-pituitary-thyroid (HPT) axis, and differential changes of TH metabolism, are emphasized in particular. It was proposed that the variability between brain and liver tissues in THs uptake, deiodination activation and inactivation is the key regulatory mechanism for the cause of peripheral THs decline and central homeostasis. While hypothalamic tanycytes play a pivotal role in the fine regulation of the HPT negative feedback regulation during long-term fasting. The study progress of tanycytes on thyrotropin-releasing hormone (TRH) release and deiodination is described in detail. In conclusion, the combination of the decrease of TH metabolism in peripheral tissues and stability in the central HPT axis maintains the basal physiological requirement and new energy metabolism homeostasis to adapt to long-term food scarcity. The molecular mechanisms of this localized and differential regulation will be a key research direction for developing measures for hypometabolic applications in extreme environment.

Dendritically localized RNAs are packaged as diversely composed ribonucleoprotein particles with heterogeneous copy number states

Localization of mRNAs to dendrites is a fundamental mechanism by which neurons achieve spatiotemporal control of gene expression. Translationally repressed neuronal mRNA transport granules, also referred to as ribonucleoprotein particles (RNPs), have been shown to be trafficked as single or low copy number RNPs and as larger complexes with multiple copies and/or species of mRNAs. However, there is little evidence of either population in intact neuronal circuits. Using single molecule fluorescence in situ hybridization studies in the dendrites of adult rat and mouse hippocampus, we provide evidence that supports the existence of multi-transcript RNPs with the constituents varying in amounts for each RNA species. By competing-off fluorescently labeled probe with serial increases of unlabeled probe, we detected stepwise decreases in Arc RNP number and fluorescence intensity, suggesting Arc RNAs localize to dendrites in both low- and multiple-copy number RNPs. When probing for multiple mRNAs, we find that localized RNPs are heterogeneous in size and colocalization patterns that vary per RNA. Further, localized RNAs that are targeted by the same trans-acting element (FMRP) display greater levels of colocalization compared to an RNA not targeted by FMRP. Simultaneous visualization of a dozen FMRP-targeted mRNA species using highly multiplexed imaging demonstrates that dendritic RNAs are mostly trafficked as heteromeric cargoes of multiple types of RNAs (at least one or more RNAs). Moreover, the composition of these RNA cargoes, as assessed by colocalization, correlates with the abundance of the transcripts even after accounting for the expected differences in colocalization based on expression. Collectively, these results suggest that dendritic RNPs are packaged as heterogeneous co-assemblies of different mRNAs and that RNP contents may be driven, at least partially, by highly abundant dendritic RNAs; a model that favors efficiency over fine-tuned control for sustaining long-distance trafficking of thousands of messenger molecules.

Opto-seq reveals input-specific immediate-early gene induction in ventral tegmental area cell types

The ventral tegmental area (VTA) is a critical node in circuits governing motivated behavior and is home to diverse populations of neurons that release dopamine, gamma-aminobutyric acid (GABA), glutamate, or combinations of these neurotransmitters. The VTA receives inputs from many brain regions, but a comprehensive understanding of input-specific activation of VTA neuronal subpopulations is lacking. To address this, we combined optogenetic stimulation of select VTA inputs with single-nucleus RNA sequencing (snRNA-seq) and highly multiplexed in situ hybridization to identify distinct neuronal clusters and characterize their spatial distribution and activation patterns. Quantification of immediate-early gene (IEG) expression revealed that different inputs activated select VTA subpopulations, which demonstrated cell-type-specific transcriptional programs. Within dopaminergic subpopulations, IEG induction levels correlated with differential expression of ion channel genes. This new transcriptomics-guided circuit analysis reveals the diversity of VTA activation driven by distinct inputs and provides a resource for future analysis of VTA cell types.

Nova-ST: Nano-patterned ultra-dense platform for spatial transcriptomics

Spatial transcriptomics workflows using barcoded capture arrays are commonly used for resolving gene expression in tissues. However, existing techniques are either limited by capture array density or are cost prohibitive for large-scale atlasing. We present Nova-ST, a dense nano-patterned spatial transcriptomics technique derived from randomly barcoded Illumina sequencing flow cells. Nova-ST enables customized, low-cost, flexible, and high-resolution spatial profiling of large tissue sections. Benchmarking on mouse brain sections demonstrates significantly higher sensitivity compared to existing methods at a reduced cost.

Cyclase-associated protein (CAP) inhibits inverted formin 2 (INF2) to induce dendritic spine maturation

The morphology of dendritic spines, the postsynaptic compartment of most excitatory synapses, decisively modulates the function of neuronal circuits as also evident from human brain disorders associated with altered spine density or morphology. Actin filaments (F-actin) form the backbone of spines, and a number of actin-binding proteins (ABP) have been implicated in shaping the cytoskeleton in mature spines. Instead, only little is known about the mechanisms that control the reorganization from unbranched F-actin of immature spines to the complex, highly branched cytoskeleton of mature spines. Here, we demonstrate impaired spine maturation in hippocampal neurons upon genetic inactivation of cyclase-associated protein 1 (CAP1) and CAP2, but not of CAP1 or CAP2 alone. We found a similar spine maturation defect upon overactivation of inverted formin 2 (INF2), a nucleator of unbranched F-actin with hitherto unknown synaptic function. While INF2 overactivation failed in altering spine density or morphology in CAP-deficient neurons, INF2 inactivation largely rescued their spine defects. From our data we conclude that CAPs inhibit INF2 to induce spine maturation. Since we previously showed that CAPs promote cofilin1-mediated cytoskeletal remodeling in mature spines, we identified them as a molecular switch that control transition from filopodia-like to mature spines.

Mesenchymal stem cell-derived small extracellular vesicles alleviate the immunometabolic dysfunction in murine septic encephalopathy

Sepsis is a life-threatening organ dysfunction caused by a dysregulated host response to infection that results in high mortality and long-term sequela. The central nervous system (CNS) is susceptible to injury from infectious processes, which can lead to clinical symptoms of septic encephalopathy (SE). SE is linked to a profound energetic deficit associated with immune dysregulation. Here, we show that intravenous administration of adipose tissue mesenchymal stem cell (MSC)-derived small extracellular vesicles (sEVs) in septic mice improved disease outcomes by reducing SE clinical severity, restoring aerobic metabolism, and lowering pro-inflammatory cytokines in the cerebellum, a key region affected by SE. Our high throughput analysis showed that MSC-derived sEVs partially reversed sepsis-induced transcriptomic changes, highlighting the potential association of miRNA regulators in the cerebellum of MSC-derived sEV-treated mice with miRNAs identified in sEV cargo. MSC-derived sEVs could serve as a promising therapeutic agent in SE through their favorable immunometabolic properties.

Benchmarking computational methods for single-cell chromatin data analysis

BACKGROUND

Single-cell chromatin accessibility assays, such as scATAC-seq, are increasingly employed in individual and joint multi-omic profiling of single cells. As the accumulation of scATAC-seq and multi-omics datasets continue, challenges in analyzing such sparse, noisy, and high-dimensional data become pressing. Specifically, one challenge relates to optimizing the processing of chromatin-level measurements and efficiently extracting information to discern cellular heterogeneity. This is of critical importance, since the identification of cell types is a fundamental step in current single-cell data analysis practices.

RESULTS

We benchmark 8 feature engineering pipelines derived from 5 recent methods to assess their ability to discover and discriminate cell types. By using 10 metrics calculated at the cell embedding, shared nearest neighbor graph, or partition levels, we evaluate the performance of each method at different data processing stages. This comprehensive approach allows us to thoroughly understand the strengths and weaknesses of each method and the influence of parameter selection.

CONCLUSIONS

Our analysis provides guidelines for choosing analysis methods for different datasets. Overall, feature aggregation, SnapATAC, and SnapATAC2 outperform latent semantic indexing-based methods. For datasets with complex cell-type structures, SnapATAC and SnapATAC2 are preferred. With large datasets, SnapATAC2 and ArchR are most scalable.

Cell type-specific dysregulation of gene expression due to Chd8 haploinsufficiency during mouse cortical development

Disruptive variants in the chromodomain helicase CHD8, which acts as a transcriptional regulator during neurodevelopment, are strongly associated with risk for autism spectrum disorder (ASD). Loss of CHD8 function is hypothesized to perturb gene regulatory networks in the developing brain, thereby contributing to ASD etiology. However, insight into the cell type-specific transcriptional effects of CHD8 loss of function remains limited. We used single-cell and single-nucleus RNA-sequencing to globally profile gene expression and identify dysregulated genes in the embryonic and juvenile wild type and Chd8 +/- mouse cortex, respectively. Chd8 and other ASD risk-associated genes showed a convergent expression trajectory that was largely conserved between the mouse and human developing cortex, increasing from the progenitor zones to the cortical plate. Genes associated with risk for neurodevelopmental disorders and genes involved in neuron projection development, chromatin remodeling, signaling, and migration were dysregulated in Chd8 +/- embryonic day (E) 12.5 radial glia. Genes implicated in synaptic organization and activity were dysregulated in Chd8 +/- postnatal day (P) 25 deep- and upper-layer excitatory cortical neurons, suggesting a delay in synaptic maturation or impaired synaptogenesis due to CHD8 loss of function. Our findings reveal a complex pattern of transcriptional dysregulation in Chd8 +/- developing cortex, potentially with distinct biological impacts on progenitors and maturing neurons in the excitatory neuronal lineage.

High-resolution diffusion magnetic resonance imaging and spatial-transcriptomic in developing mouse brain

Brain development is a highly complex process regulated by numerous genes at the molecular and cellular levels. Brain tissue exhibits serial microstructural changes during the development process. High-resolution diffusion magnetic resonance imaging (dMRI) affords a unique opportunity to probe these changes in the developing brain non-destructively. In this study, we acquired multi-shell dMRI datasets at 32 µm isotropic resolution to investigate the tissue microstructure alterations, which we believe to be the highest spatial resolution dMRI datasets obtained for postnatal mouse brains. We adapted the Allen Developing Mouse Brain Atlas (ADMBA) to integrate quantitative MRI metrics and spatial transcriptomics. Diffusion tensor imaging (DTI), diffusion kurtosis imaging (DKI), and neurite orientation dispersion and density imaging (NODDI) metrics were used to quantify brain development at different postnatal days. We demonstrated that the differential evolutions of fiber orientation distributions contribute to the distinct development patterns in white matter (WM) and gray matter (GM). Furthermore, the genes enriched in the nervous system that regulate brain structure and function were expressed in spatial correlation with age-matched dMRI. This study is the first one providing high-resolution dMRI, including DTI, DKI, and NODDI models, to trace mouse brain microstructural changes in WM and GM during postnatal development. This study also highlighted the genotype-phenotype correlation of spatial transcriptomics and dMRI, which may improve our understanding of brain microstructure changes at the molecular level.

Current insights and assumptions on α-synuclein in Lewy body disease

Lewy body disorders are heterogeneous neurological conditions defined by intracellular inclusions composed of misshapen α-synuclein protein aggregates. Although α-synuclein aggregates are only one component of inclusions and not strictly coupled to neurodegeneration, evidence suggests they seed the propagation of Lewy pathology within and across cells. Genetic mutations, genomic multiplications, and sequence polymorphisms of the gene encoding α-synuclein are also causally linked to Lewy body disease. In nonfamilial cases of Lewy body disease, the disease trigger remains unidentified but may range from industrial/agricultural toxicants and natural sources of poisons to microbial pathogens. Perhaps due to these peripheral exposures, Lewy inclusions appear at early disease stages in brain regions connected with cranial nerves I and X, which interface with inhaled and ingested environmental elements in the nasal or gastrointestinal cavities. Irrespective of its identity, a stealthy disease trigger most likely shifts soluble α-synuclein (directly or indirectly) into insoluble, cross-β-sheet aggregates. Indeed, β-sheet-rich self-replicating α-synuclein multimers reside in patient plasma, cerebrospinal fluid, and other tissues, and can be subjected to α-synuclein seed amplification assays. Thus, clinicians should be able to capitalize on α-synuclein seed amplification assays to stratify patients into potential responders versus non-responders in future clinical trials of α-synuclein targeted therapies. Here, we briefly review the current understanding of α-synuclein in Lewy body disease and speculate on pathophysiological processes underlying the potential transmission of α-synucleinopathy across the neuraxis.

Dorsoventral Heterogeneity of Synaptic Connectivity in Hippocampal CA3 Pyramidal Neurons

The hippocampal CA3 region plays an important role in learning and memory. CA3 pyramidal neurons (PNs) receive two prominent excitatory inputs-mossy fibers (MFs) from dentate gyrus (DG) and recurrent collaterals (RCs) from CA3 PNs-that play opposing roles in pattern separation and pattern completion, respectively. Although the dorsoventral heterogeneity of the hippocampal anatomy, physiology, and behavior has been well established, nothing is known about the dorsoventral heterogeneity of synaptic connectivity in CA3 PNs. In this study, we performed Timm's sulfide silver staining, dendritic and spine morphological analyses, and ex vivo electrophysiology in mice of both sexes to investigate the heterogeneity of MF and RC pathways along the CA3 dorsoventral axis. Our morphological analyses demonstrate that ventral CA3 (vCA3) PNs possess greater dendritic lengths and more complex dendritic arborization, compared with dorsal CA3 (dCA3) PNs. Moreover, using ChannelRhodopsin2 (ChR2)-assisted patch-clamp recording, we found that the ratio of the RC-to-MF excitatory drive onto CA3 PNs increases substantially from dCA3 to vCA3, with vCA3 PNs receiving significantly weaker MFs, but stronger RCs, excitation than dCA3 PNs. Given the distinct roles of MF versus RC inputs in pattern separation versus completion, our findings of the significant dorsoventral variations of MF and RC excitation in CA3 PNs may have important functional implications for the contribution of CA3 circuit to the dorsoventral difference in hippocampal function.

Prefrontal Cortex subregions provide distinct visual and behavioral feedback modulation to the Primary Visual Cortex

The mammalian Prefrontal Cortex (PFC) has been suggested to modulate sensory information processing across multiple cortical regions via long-range axonal projections. These axonal projections arise from PFC subregions with unique brain-wide connectivity and functional repertoires, which may provide the architecture for modular feedback intended to shape sensory processing. Here, we used axonal tracing, axonal and somatic 2-photon calcium imaging, and chemogenetic manipulations in mice to delineate how projections from the Anterior Cingulate Cortex (ACA) and ventrolateral Orbitofrontal Cortex (ORB) of the PFC modulate sensory processing in the primary Visual Cortex (VISp) across behavioral states. Structurally, we found that ACA and ORB have distinct patterning of projections across both cortical regions and layers. ACA axons in VISp had a stronger representation of visual stimulus information than ORB axons, but both projections showed non-visual, behavior-dependent activity. ACA input to VISp enhanced the encoding of visual stimuli by VISp neurons, and modulation of visual responses scaled with arousal. On the other hand, ORB input shaped movement and arousal related modulation of VISp visual responses, but specifically reduced the encoding of high-contrast visual stimuli. Thus, ACA and ORB feedback have separable projection patterns and encode distinct visual and behavioral information, putatively providing the substrate for their unique effects on visual representations and behavioral modulation in VISp. Our results offer a refined model of cortical hierarchy and its impact on sensory information processing, whereby distinct as opposed to generalized properties of PFC projections contribute to VISp activity during discrete behavioral states.

Benchmarking clustering, alignment, and integration methods for spatial transcriptomics

BACKGROUND

Spatial transcriptomics (ST) is advancing our understanding of complex tissues and organisms. However, building a robust clustering algorithm to define spatially coherent regions in a single tissue slice and aligning or integrating multiple tissue slices originating from diverse sources for essential downstream analyses remains challenging. Numerous clustering, alignment, and integration methods have been specifically designed for ST data by leveraging its spatial information. The absence of comprehensive benchmark studies complicates the selection of methods and future method development.

RESULTS

In this study, we systematically benchmark a variety of state-of-the-art algorithms with a wide range of real and simulated datasets of varying sizes, technologies, species, and complexity. We analyze the strengths and weaknesses of each method using diverse quantitative and qualitative metrics and analyses, including eight metrics for spatial clustering accuracy and contiguity, uniform manifold approximation and projection visualization, layer-wise and spot-to-spot alignment accuracy, and 3D reconstruction, which are designed to assess method performance as well as data quality. The code used for evaluation is available on our GitHub. Additionally, we provide online notebook tutorials and documentation to facilitate the reproduction of all benchmarking results and to support the study of new methods and new datasets.

CONCLUSIONS

Our analyses lead to comprehensive recommendations that cover multiple aspects, helping users to select optimal tools for their specific needs and guide future method development.

Haploinsufficiency of GABAA Receptor-Associated Clptm1 Enhances Phasic and Tonic Inhibitory Neurotransmission, Suppresses Excitatory Synaptic Plasticity, and Impairs Memory

The mechanisms utilized by neurons to regulate the efficacy of phasic and tonic inhibition and their impacts on synaptic plasticity and behavior are incompletely understood. Cleft lip and palate transmembrane protein 1 (Clptm1) is a membrane-spanning protein that interacts with multiple γ-aminobutyric acid type A receptor (GABAAR) subunits, trapping them in the endoplasmic reticulum and Golgi network. Overexpression and knock-down studies suggest that Clptm1 modulates GABAAR-mediated phasic inhibition and tonic inhibition as well as activity-induced inhibitory synaptic homeostasis in cultured hippocampal neurons. To investigate the role of Clptm1 in the modulation of GABAARs in vivo, we generated Clptm1 knock-out (KO) mice. Here, we show that genetic KO of Clptm1 elevated phasic and tonic inhibitory transmission in both male and female heterozygous mice. Although basal excitatory synaptic transmission was not affected, Clptm1 haploinsufficiency significantly blocked high-frequency stimulation-induced long-term potentiation (LTP) in hippocampal CA3→CA1 synapses. In the hippocampus-dependent contextual fear-conditioning behavior task, both male and female Clptm1 heterozygous KO mice exhibited impairment in contextual fear memory. In addition, LTP and contextual fear memory were rescued by application of L-655,708, a negative allosteric modulator of the extrasynaptic GABAAR α5 subunit. These results suggest that haploinsufficiency of Clptm1 contributes to cognitive deficits through altered synaptic transmission and plasticity by elevation of inhibitory neurotransmission, with tonic inhibition playing a major role.

The loss of βΙ spectrin alters synaptic size and composition in the ja/ja mouse

Introduction

Deletion or mutation of members of the spectrin gene family contributes to many neurologic and neuropsychiatric disorders. While each spectrinopathy may generate distinct neuropathology, the study of βΙ spectrin's role (Sptb) in the brain has been hampered by the hematologic consequences of its loss.

Methods

Jaundiced mice (ja/ja) that lack βΙ spectrin suffer a rapidly fatal hemolytic anemia. We have used exchange transfusion of newborn ja/ja mice to blunt their hemolytic pathology, enabling an examination of βΙ spectrin deficiency in the mature mouse brain by ultrastructural and biochemical analysis.

Results

βΙ spectrin is widely utilized throughout the brain as the βΙΣ2 isoform; it appears by postnatal day 8, and concentrates in the CA1,3 region of the hippocampus, dentate gyrus, cerebellar granule layer, cortical layer 2, medial habenula, and ventral thalamus. It is present in a subset of dendrites and absent in white matter. Without βΙ spectrin there is a 20% reduction in postsynaptic density size in the granule layer of the cerebellum, a selective loss of ankyrinR in cerebellar granule neurons, and a reduction in the level of the postsynaptic adhesion molecule NCAM. While we find no substitution of another spectrin for βΙ at dendrites or synapses, there is curiously enhanced βΙV spectrin expression in the ja/ja brain.

Discussion

βΙΣ2 spectrin appears to be essential for refining postsynaptic structures through interactions with ankyrinR and NCAM. We speculate that it may play additional roles yet to be discovered.

Fasting induces metabolic switches and spatial redistributions of lipid processing and neuronal interactions in tanycytes

The ependyma lining the third ventricle (3V) in the mediobasal hypothalamus plays a crucial role in energy balance and glucose homeostasis. It is characterized by a high functional heterogeneity and plasticity, but the underlying molecular mechanisms governing its features are not fully understood. Here, 5481 hypothalamic ependymocytes were cataloged using FACS-assisted scRNAseq from fed, 12h-fasted, and 24h-fasted adult male mice. With standard clustering analysis, typical ependymal cells and β2-tanycytes appear sharply defined, but other subpopulations, β1- and α-tanycytes, display fuzzy boundaries with few or no specific markers. Pseudospatial approaches, based on the 3V neuroanatomical distribution, enable the identification of specific versus shared tanycyte markers and subgroup-specific versus general tanycyte functions. We show that fasting dynamically shifts gene expression patterns along the 3V, leading to a spatial redistribution of cell type-specific responses. Altogether, we show that changes in energy status induce metabolic and functional switches in tanycyte subpopulations, providing insights into molecular and functional diversity and plasticity within the tanycyte population.

Dependency-aware deep generative models for multitasking analysis of spatial omics data

Spatially resolved transcriptomics (SRT) technologies have significantly advanced biomedical research, but their data analysis remains challenging due to the discrete nature of the data and the high levels of noise, compounded by complex spatial dependencies. Here, we propose spaVAE, a dependency-aware, deep generative spatial variational autoencoder model that probabilistically characterizes count data while capturing spatial correlations. spaVAE introduces a hybrid embedding combining a Gaussian process prior with a Gaussian prior to explicitly capture spatial correlations among spots. It then optimizes the parameters of deep neural networks to approximate the distributions underlying the SRT data. With the approximated distributions, spaVAE can contribute to several analytical tasks that are essential for SRT data analysis, including dimensionality reduction, visualization, clustering, batch integration, denoising, differential expression, spatial interpolation, resolution enhancement and identification of spatially variable genes. Moreover, we have extended spaVAE to spaPeakVAE and spaMultiVAE to characterize spatial ATAC-seq (assay for transposase-accessible chromatin using sequencing) data and spatial multi-omics data, respectively.

Enhanced In Situ Spatial Proteomics by Effective Combination of MALDI Imaging and LC-MS/MS

Highly specialized cells are fundamental for the proper functioning of complex organs. Variations in cell-type-specific gene expression and protein composition have been linked to a variety of diseases. Investigation of the distinctive molecular makeup of these cells within tissues is therefore critical in biomedical research. Although several technologies have emerged as valuable tools to address this cellular heterogeneity, most workflows lack sufficient in situ resolution and are associated with high costs and extremely long analysis times. Here, we present a combination of experimental and computational approaches that allows a more comprehensive investigation of molecular heterogeneity within tissues than by either shotgun LC-MS/MS or MALDI imaging alone. We applied our pipeline to the mouse brain, which contains a wide variety of cell types that not only perform unique functions but also exhibit varying sensitivities to insults. We explored the distinct neuronal populations within the hippocampus, a brain region crucial for learning and memory that is involved in various neurological disorders. As an example, we identified the groups of proteins distinguishing the neuronal populations of the dentate gyrus (DG) and the cornu ammonis (CA) in the same brain section. Most of the annotated proteins matched the regional enrichment of their transcripts, thereby validating the method. As the method is highly reproducible, the identification of individual masses through the combination of MALDI-IMS and LC-MS/MS methods can be used for the much faster and more precise interpretation of MALDI-IMS measurements only. This greatly speeds up spatial proteomic analyses and allows the detection of local protein variations within the same population of cells. The method's general applicability has the potential to be used to investigate different biological conditions and tissues and a much higher throughput than other techniques making it a promising approach for clinical routine applications.

Neutral sphingomyelinase controls acute and chronic alcohol effects on brain activity

Alcohol consumption is a widespread phenomenon throughout the world. However, how recreational alcohol use evolves into alcohol use disorder (AUD) remains poorly understood. The Smpd3 gene and its coded protein neutral sphingomyelinase (NSM) are associated with alcohol consumption in humans and alcohol-related behaviors in mice, suggesting a potential role in this transition. Using multiparametric magnetic resonance imaging, we characterized the role of NSM in acute and chronic effects of alcohol on brain anatomy and function in female mice. Chronic voluntary alcohol consumption (16 vol% for at least 6 days) affected brain anatomy in WT mice, reducing regional structure volume predominantly in cortical regions. Attenuated NSM activity prevented these anatomical changes. Functional MRI linked these anatomical adaptations to functional changes: Chronic alcohol consumption in mice significantly modulated resting state functional connectivity (RS FC) in response to an acute ethanol challenge (i.p. bolus of 2 g kg-1) in heterozygous NSM knockout (Fro), but not in WT mice. Acute ethanol administration in alcohol-naïve WT mice significantly decreased RS FC in cortical and brainstem regions, a key finding that was amplified in Fro mice. Regarding direct pharmacological effects, acute ethanol administration increased the regional cerebral blood volume (rCBV) in many brain areas. Here, chronic alcohol consumption otherwise attenuated the acute rCBV response in WT mice but enhanced it in Fro mice. Altogether, these findings suggest a differential role for NSM in acute and chronic functional brain responses to alcohol. Therefore, targeting NSM may be useful in the prevention or treatment of AUD.

Genetic control of MRI contrast using the manganese transporter Zip14

PURPOSE

Gene-expression reporter systems, such as green fluorescent protein, have been instrumental to understanding biological processes in living organisms at organ system, tissue, cell, and molecular scales. More than 30 years of work on developing MRI-visible gene-expression reporter systems has resulted in a variety of clever application-specific methods. However, these techniques have not yet been widely adopted, so a general-purpose expression reporter is still required. Here, we demonstrate that the manganese ion transporter Zip14 is an in vivo MRI-visible, flexible, and robust gene-expression reporter to meet this need.

METHODS

Plasmid constructs consisting of a cell type-specific promoter, gene coding for human Zip14, and a histology-visible tag were packaged into adeno-associated viruses. These viruses were intracranially injected into the mouse brain. Serial in vivo MRI was performed using a vendor-supplied 3D-MPRAGE sequence. No additional contrast agents were administered. Animals were sacrificed after the last imaging timepoint for immunohistological validation.

RESULTS

Neuron-specific overexpression of Zip14 produced substantial and long-lasting changes in MRI contrast. Using appropriate viruses enabled both anterograde and retrograde neural tracing. Expression of Zip14 in astrocytes also enabled MRI of glia populations in the living mammalian brain.

CONCLUSIONS

The flexibility of this system as an MRI-visible gene-expression reporter will enable many applications of serial, high-resolution imaging of gene expression for basic science and therapy development.

Neuroanatomical and cognitive biomarkers of alpha-synuclein propagation in a mouse model of synucleinopathy prior to onset of motor symptoms

Significant evidence suggests that misfolded alpha-synuclein (aSyn), a major component of Lewy bodies, propagates in a prion-like manner contributing to disease progression in Parkinson's disease (PD) and other synucleinopathies. In fact, timed inoculation of M83 hemizygous mice with recombinant human aSyn preformed fibrils (PFF) has shown symptomatic deficits after substantial spreading of pathogenic alpha-synuclein, as detected by markers for the phosphorylation of S129 of aSyn. However, whether accumulated toxicity impact human-relevant cognitive and structural neuroanatomical measures is not fully understood. Here we performed a single unilateral striatal PFF injection in M83 hemizygous mice, and using two assays with translational potential, ex vivo magnetic resonance imaging (MRI) and touchscreen testing, we examined the combined neuroanatomical and behavioral impact of aSyn propagation. In PFF-injected mice, we observed widespread atrophy in bilateral regions that project to or receive input from the injection site using MRI. We also identified early deficits in reversal learning prior to the emergence of motor symptoms. Our findings highlight a network of regions with related cellular correlates of pathology that follow the progression of aSyn spreading, and that affect brain areas relevant for reversal learning. Our experiments suggest that M83 hemizygous mice injected with human PFF provides a model to understand how misfolded aSyn affects human-relevant pre-clinical measures and suggest that these pre-clinical biomarkers could be used to detect early toxicity of aSyn and provide better translational measures between mice and human disease.

Molecular Organization of Autonomic, Respiratory, and Spinally-Projecting Neurons in the Mouse Ventrolateral Medulla

The ventrolateral medulla (VLM) is a crucial region in the brain for visceral and somatic control, serving as a significant source of synaptic input to the spinal cord. Experimental studies have shown that gene expression in individual VLM neurons is predictive of their function. However, the molecular and cellular organization of the VLM has remained uncertain. This study aimed to create a comprehensive dataset of VLM cells using single-cell RNA sequencing in male and female mice. The dataset was enriched with targeted sequencing of spinally-projecting and adrenergic/noradrenergic VLM neurons. Based on differentially expressed genes, the resulting dataset of 114,805 VLM cells identifies 23 subtypes of neurons, excluding those in the inferior olive, and five subtypes of astrocytes. Spinally-projecting neurons were found to be abundant in seven subtypes of neurons, which were validated through in situ hybridization. These subtypes included adrenergic/noradrenergic neurons, serotonergic neurons, and neurons expressing gene markers associated with premotor neurons in the ventromedial medulla. Further analysis of adrenergic/noradrenergic neurons and serotonergic neurons identified nine and six subtypes, respectively, within each class of monoaminergic neurons. Marker genes that identify the neural network responsible for breathing were concentrated in two subtypes of neurons, delineated from each other by markers for excitatory and inhibitory neurons. These datasets are available for public download and for analysis with a user-friendly interface. Collectively, this study provides a fine-scale molecular identification of cells in the VLM, forming the foundation for a better understanding of the VLM's role in vital functions and motor control.

Hierarchical Gradients of Encoded Spatial and Sensory Information in the Neocortex Are Attenuated by Dorsal Hippocampal Lesions

During navigation, the neocortex actively integrates learned spatial context with current sensory experience to guide behaviors. However, the relative encoding of spatial and sensorimotor information among cortical cells, and whether hippocampal feedback continues to modify these properties after learning, remains poorly understood. Thus, two-photon microscopy of male and female Thy1-GCaMP6s mice was used to longitudinally image neurons spanning superficial retrosplenial cortex and layers II-Va of primary and secondary motor cortices before and after bilateral dorsal hippocampal lesions. During behavior on a familiar cued treadmill, the locations of two obstacles were interchanged to decouple place-tuning from cue-tuning among position-correlated cells with fields at those locations. Subpopulations of place and cue cells each formed interareal gradients such that higher-level cortical regions exhibited higher fractions of place cells, whereas lower-level regions exhibited higher fractions of cue cells. Position-correlated cells in the motor cortex also formed translaminar gradients; more superficial cells were more likely to exhibit fields and were more sparsely and precisely tuned than deeper cells. After dorsal hippocampal lesions, a neural representation of the learned environment persisted, but retrosplenial cortex exhibited significantly increased cue-tuning, and, in motor cortices, both position-correlated cell recruitment and population activity at the unstable obstacle locations became more homogeneously elevated across laminae. Altogether, these results support that the hippocampus continues to modulate cortical responses in familiar environments, and the relative impact of descending feedback obeys hierarchical interareal and interlaminar gradients opposite to the flow of ascending sensory inputs.

Environmental enrichment reduces restricted repetitive behavior by altering gray matter microstructure

Restricted, repetitive behaviors are common symptoms in neurodevelopmental disorders including autism spectrum disorder. Despite being associated with poor developmental outcomes, repetitive behaviors remain poorly understood and have limited treatment options. Environmental enrichment attenuates the development of repetitive behaviors, but the exact mechanisms remain obscure. Using the C58 mouse model of repetitive behavior, we performed diffusion tensor imaging to examine microstructural alterations associated with the development of repetitive behavior and its attenuation by environmental enrichment. The C57BL/6 mouse strain, which displays little or no repetitive behavior, was used as a control group. We observed widespread differences in diffusion metrics between C58 mice and C57BL/6 mice. In juvenile C58 mice, repetitive motor behavior displayed strong negative correlations with fractional anisotropy in multiple gray matter regions, whereas in young adult C58 mice, high repetitive motor behavior was most strongly associated with lower fractional anisotropy and higher radial diffusivity in the striatum. Environmental enrichment increased fractional anisotropy and axial diffusivity throughout gray matter regions in the brains of juvenile C58 mice and overlapped predominantly with cerebellar and sensory regions associated with repetitive behavior. Our results suggest environmental enrichment reduces repetitive behavior development by altering gray matter microstructure in the cerebellum, medial entorhinal cortex, and sensory processing regions in juvenile C58 mice. Under standard laboratory conditions, early pathology in these regions appears to contribute to later striatal and white matter dysfunction in adult C58 mice. Future studies should examine the role these regions play in the development of repetitive behavior and the relationship between sensory processing and cerebellar deficits and repetitive behavior.

The nucleus accumbens shell: a neural hub at the interface of homeostatic and hedonic feeding

Feeding behavior is a complex physiological process regulated by the interplay between homeostatic and hedonic feeding circuits. Among the neural structures involved, the nucleus accumbens (NAc) has emerged as a pivotal region at the interface of these two circuits. The NAc comprises distinct subregions and in this review, we focus mainly on the NAc shell (NAcSh). Homeostatic feeding circuits, primarily found in the hypothalamus, ensure the organism's balance in energy and nutrient requirements. These circuits monitor peripheral signals, such as insulin, leptin, and ghrelin, and modulate satiety and hunger states. The NAcSh receives input from these homeostatic circuits, integrating information regarding the organism's metabolic needs. Conversely, so-called hedonic feeding circuits involve all other non-hunger and -satiety processes, i.e., the sensory information, associative learning, reward, motivation and pleasure associated with food consumption. The NAcSh is interconnected with hedonics-related structures like the ventral tegmental area and prefrontal cortex and plays a key role in encoding hedonic information related to palatable food seeking or consumption. In sum, the NAcSh acts as a crucial hub in feeding behavior, integrating signals from both homeostatic and hedonic circuits, to facilitate behavioral output via its downstream projections. Moreover, the NAcSh's involvement extends beyond simple integration, as it directly impacts actions related to food consumption. In this review, we first focus on delineating the inputs targeting the NAcSh; we then present NAcSh output projections to downstream structures. Finally we discuss how the NAcSh regulates feeding behavior and can be seen as a neural hub integrating homeostatic and hedonic feeding signals, via a functionally diverse set of projection neuron subpopulations.

Whole brain alignment of spatial transcriptomics between humans and mice with BrainAlign

The increasing utilization of mouse models in human neuroscience research places higher demands on computational methods to translate findings from the mouse brain to the human one. In this study, we develop BrainAlign, a self-supervised learning approach, for the whole brain alignment of spatial transcriptomics (ST) between humans and mice. BrainAlign encodes spots and genes simultaneously in two separated shared embedding spaces by a heterogeneous graph neural network. We demonstrate that BrainAlign could integrate cross-species spots into the embedding space and reveal the conserved brain regions supported by ST information, which facilitates the detection of homologous regions between humans and mice. Genomic analysis further presents gene expression connections between humans and mice and reveals similar expression patterns for marker genes. Moreover, BrainAlign can accurately map spatially similar homologous regions or clusters onto a unified spatial structural domain while preserving their relative positions.

Graded control of Purkinje cell outputs by cAMP through opposing actions on axonal action potential and transmitter release

All-or-none signalling by action potentials (APs) in neuronal axons is pivotal for the precisely timed and identical size of outputs to multiple distant targets. However, technical limitations with respect to measuring the signalling in small intact axons have hindered the evaluation of high-fidelity signal propagation. Here, using direct recordings from axonal trunks and/or terminals of cerebellar Purkinje cells in slice and culture, we demonstrate that the timing and amplitude of axonal outputs are gradually modulated by cAMP depending on the length of axon. During the propagation in long axon, APs were attenuated and slowed in conduction by cAMP via specifically decreasing axonal Na+ currents. Consequently, the Ca2+ influx and transmitter release at distal boutons are reduced by cAMP, counteracting its direct facilitating effect on release machinery as observed at various CNS synapses. Together, our tour de force functional dissection has unveiled the axonal distance-dependent graded control of output timing and strength by intracellular signalling. KEY POINTS: The information processing in the nervous system has been classically thought to rely on the axonal faithful and high-speed conduction of action potentials (APs). We demonstrate that the strength and timing of axonal outputs are weakened and delayed, respectively, by cytoplasmic cAMP depending on the axonal length in cerebellar Purkinje cells (PCs). Direct axonal patch clamp recordings uncovered axon-specific attenuation of APs by cAMP through reduction of axonal Na+ currents. cAMP directly augments transmitter release at PC terminals without changing presynaptic Ca2+ influx or readily releasable pool of vesicles, although the extent is weaker compared to other CNS synapses. Two opposite actions of cAMP on PC axons, AP attenuation and release augmentation, together give rise to graded control of synaptic outputs in a manner dependent on the axonal length.

Oxytocin induces the formation of distinctive cortical representations and cognitions biased toward familiar mice

Social recognition is essential for the formation of social structures. Many times, recognition comes with lesser exploration of familiar animals. This lesser exploration has led to the assumption that recognition may be a habituation memory. The underlying memory mechanisms and the thereby acquired cortical representations of familiar mice have remained largely unknown, however. Here, we introduce an approach directly examining the recognition process from volatile body odors among male mice. We show that volatile body odors emitted by mice are sufficient to identify individuals and that more salience is assigned to familiar mice. Familiarity is encoded by reinforced population responses in two olfactory cortex hubs and communicated to other brain regions. The underlying oxytocin-induced plasticity promotes the separation of the cortical representations of familiar from other mice. In summary, neuronal encoding of familiar animals is distinct and utilizes the cortical representational space more broadly, promoting storage of complex social relationships.

Spatiotemporal transcriptomic profiling and modeling of mouse brain at single-cell resolution reveals cell proximity effects of aging and rejuvenation

Old age is associated with a decline in cognitive function and an increase in neurodegenerative disease risk1. Brain aging is complex and accompanied by many cellular changes2-20. However, the influence that aged cells have on neighboring cells and how this contributes to tissue decline is unknown. More generally, the tools to systematically address this question in aging tissues have not yet been developed. Here, we generate spatiotemporal data at single-cell resolution for the mouse brain across lifespan, and we develop the first machine learning models based on spatial transcriptomics ('spatial aging clocks') to reveal cell proximity effects during brain aging and rejuvenation. We collect a single-cell spatial transcriptomics brain atlas of 4.2 million cells from 20 distinct ages and across two rejuvenating interventions-exercise and partial reprogramming. We identify spatial and cell type-specific transcriptomic fingerprints of aging, rejuvenation, and disease, including for rare cell types. Using spatial aging clocks and deep learning models, we find that T cells, which infiltrate the brain with age, have a striking pro-aging proximity effect on neighboring cells. Surprisingly, neural stem cells have a strong pro-rejuvenating effect on neighboring cells. By developing computational tools to identify mediators of these proximity effects, we find that pro-aging T cells trigger a local inflammatory response likely via interferon-γ whereas pro-rejuvenating neural stem cells impact the metabolism of neighboring cells possibly via growth factors (e.g. vascular endothelial growth factor) and extracellular vesicles, and we experimentally validate some of these predictions. These results suggest that rare cells can have a drastic influence on their neighbors and could be targeted to counter tissue aging. We anticipate that these spatial aging clocks will not only allow scalable assessment of the efficacy of interventions for aging and disease but also represent a new tool for studying cell-cell interactions in many spatial contexts.

Medial preoptic circuits governing instinctive social behaviors

The medial preoptic area (MPOA) has long been implicated in maternal and male sexual behavior. Modern neuroscience methods have begun to reveal the cellular networks responsible, while also implicating the MPOA in other social behaviors, affiliative social touch, and aggression. The social interactions rely on input from conspecifics whose most important modalities in rodents are olfaction and somatosensation. These inputs bypass the cerebral cortex to reach the MPOA to influence the social function. Hormonal inputs also directly act on MPOA neurons. In turn, the MPOA controls social responses via various projections for reward and motor output. The MPOA thus emerges as one of the major brain centers for instinctive social behavior. While key elements of MPOA circuits have been identified, a synthesis of these new data is now provided for further studies to reveal the mechanisms by which the area controls social interactions.

Antibody-assisted selective isolation of Purkinje cell nuclei from mouse cerebellar tissue

We developed a method that utilizes fluorescent labeling of nuclear envelopes alongside cytometry sorting for the selective isolation of Purkinje cell (PC) nuclei. Beginning with SUN1 reporter mice, we GFP-tagged envelopes to confirm that PC nuclei could be accurately separated from other cell types. We then developed an antibody-based protocol to make PC nuclear isolation more robust and adaptable to cerebellar tissues of any genotypic background. Immunofluorescent labeling of the nuclear membrane protein RanBP2 enabled the isolation of PC nuclei from C57BL/6 cerebellum. By analyzing the expression of PC markers, nuclear size, and nucleoli number, we confirmed that our method delivers a pure fraction of PC nuclei. To demonstrate its applicability, we isolated PC nuclei from spinocerebellar ataxia type 7 (SCA7) mice and identified transcriptional changes in known and new disease-associated genes. Access to pure PC nuclei offers insights into PC biology and pathology, including the nature of selective neuronal vulnerability.

Rem2 interacts with CaMKII at synapses and restricts long-term potentiation in hippocampus

Synaptic plasticity, the process whereby neuronal connections are either strengthened or weakened in response to stereotyped forms of stimulation, is widely believed to represent the molecular mechanism that underlies learning and memory. The holoenzyme calcium/calmodulin-dependent protein kinase II (CaMKII) plays a well-established and critical role in the induction of a variety of forms of synaptic plasticity such as long-term potentiation (LTP), long-term depression (LTD) and depotentiation. Previously, we identified the GTPase Rem2 as a potent, endogenous inhibitor of CaMKII. Here, we report that knock out of Rem2 enhances LTP at the Schaffer collateral to CA1 synapse in hippocampus, consistent with an inhibitory action of Rem2 on CaMKII in vivo. Further, re-expression of WT Rem2 rescues the enhanced LTP observed in slices obtained from Rem2 conditional knock out (cKO) mice, while expression of a mutant Rem2 construct that is unable to inhibit CaMKII in vitro fails to rescue increased LTP. In addition, we demonstrate that CaMKII and Rem2 interact in dendritic spines using a 2pFLIM-FRET approach. Taken together, our data lead us to propose that Rem2 serves as a brake on synaptic potentiation via inhibition of CaMKII activity. Further, the enhanced LTP phenotype we observe in Rem2 cKO slices reveals a previously unknown role for Rem2 in the negative regulation of CaMKII function.

An Improved Postprocessing Method to Mitigate the Macroscopic Cross-Slice B0 Field Effect on R2* Measurements in the Mouse Brain at 7T

The MR transverse relaxation rate, R2*, has been widely used to detect iron and myelin content in tissue. However, it is also sensitive to macroscopic B0 inhomogeneities. One approach to correct for the B0 effect is to fit gradient-echo signals with the three-parameter model, a sinc function-weighted monoexponential decay. However, such three-parameter models are subject to increased noise sensitivity. To address this issue, this study presents a two-stage fitting procedure based on the three-parameter model to mitigate the B0 effect and reduce the noise sensitivity of R2* measurement in the mouse brain at 7T. MRI scans were performed on eight healthy mice. The gradient-echo signals were fitted with the two-stage fitting procedure to generate R2corr_t*. The signals were also fitted with the monoexponential and three-parameter models to generate R2nocorr* and R2corr*, respectively. Regions of interest (ROIs), including the corpus callosum, internal capsule, somatosensory cortex, caudo-putamen, thalamus, and lateral ventricle, were selected to evaluate the within-ROI mean and standard deviation (SD) of the R2* measurements. The results showed that the Akaike information criterion of the monoexponential model was significantly reduced by using the three-parameter model in the selected ROIs (p = 0.0039-0.0078). However, the within-ROI SD of R2corr* using the three-parameter model was significantly higher than that of the R2nocorr* in the internal capsule, caudo-putamen, and thalamus regions (p = 0.0039), a consequence partially due to the increased noise sensitivity of the three-parameter model. With the two-stage fitting procedure, the within-ROI SD of R2corr* was significantly reduced by 7.7-30.2% in all ROIs, except for the somatosensory cortex region with a fast in-plane variation of the B0 gradient field (p = 0.0039-0.0078). These results support the utilization of the two-stage fitting procedure to mitigate the B0 effect and reduce noise sensitivity for R2* measurement in the mouse brain.

GeneSurfer Enables Transcriptome-wide Exploration and Functional Annotation of Gene Co-expression Modules in 3D Spatial Transcriptomics Data

Gene co-expression provides crucial insights into biological functions, however, there is a lack of exploratory analysis tools for localized gene co-expression in large-scale datasets. We present GeneSurfer, an interactive interface designed to explore localized transcriptome-wide gene co-expression patterns in the 3D spatial domain. Key features of GeneSurfer include transcriptome-wide gene filtering and gene clustering based on spatial local co-expression within transcriptomically similar cells, multi-slice 3D rendering of average expression of gene clusters, and on-the-fly Gene Ontology term annotation of co-expressed gene sets. Additionally, GeneSurfer offers multiple linked views for investigating individual genes or gene co-expression in the spatial domain at each exploration stage. Demonstrating its utility with both spatial transcriptomics and single-cell RNA sequencing data from the Allen Brain Cell Atlas, GeneSurfer effectively identifies and annotates localized transcriptome-wide co-expression, providing biological insights and facilitating hypothesis generation and validation.

In vivo CRISPR base editing for treatment of Huntington's disease

Huntington's disease (HD) is an inherited and ultimately fatal neurodegenerative disorder caused by an expanded polyglutamine-encoding CAG repeat within exon 1 of the huntingtin (HTT) gene, which produces a mutant protein that destroys striatal and cortical neurons. Importantly, a critical event in the pathogenesis of HD is the proteolytic cleavage of the mutant HTT protein by caspase-6, which generates fragments of the N-terminal domain of the protein that form highly toxic aggregates. Given the role that proteolysis of the mutant HTT protein plays in HD, strategies for preventing this process hold potential for treating the disorder. By screening 141 CRISPR base editor variants targeting splice elements in the HTT gene, we identified platforms capable of producing HTT protein isoforms resistant to caspase-6-mediated proteolysis via editing of the splice acceptor sequence for exon 13. When delivered to the striatum of a rodent HD model, these base editors induced efficient exon skipping and decreased the formation of the N-terminal fragments, which in turn reduced HTT protein aggregation and attenuated striatal and cortical atrophy. Collectively, these results illustrate the potential for CRISPR base editing to decrease the toxicity of the mutant HTT protein for HD.

Phospholipid Scramblase 1 Controls Efficient Neurotransmission and Synaptic Vesicle Retrieval at Cerebellar Synapses

Phospholipids (PLs) are asymmetrically distributed at the plasma membrane. This asymmetric lipid distribution is transiently altered during calcium-regulated exocytosis, but the impact of this transient remodeling on presynaptic function is currently unknown. As phospholipid scramblase 1 (PLSCR1) randomizes PL distribution between the two leaflets of the plasma membrane in response to calcium activation, we set out to determine its role in neurotransmission. We report here that PLSCR1 is expressed in cerebellar granule cells (GrCs) and that PLSCR1-dependent phosphatidylserine egress occurred at synapses in response to neuron stimulation. Synaptic transmission is impaired at GrC Plscr1 -/- synapses, and both PS egress and synaptic vesicle (SV) endocytosis are inhibited in Plscr1 -/- cultured neurons from male and female mice, demonstrating that PLSCR1 controls PL asymmetry remodeling and SV retrieval following neurotransmitter release. Altogether, our data reveal a novel key role for PLSCR1 in SV recycling and provide the first evidence that PL scrambling at the plasma membrane is a prerequisite for optimal presynaptic performance.