A transcriptomic atlas of mouse cerebellar cortex comprehensively defines cell types

Lack of junctional adhesion molecule (JAM)-B traps CD8 T cells in CNS border zones and ameliorates autoimmune neuroinflammation

The endothelial blood-brain barrier (BBB) tightly controls T cell entry into the central nervous system (CNS). T cell extravasation across the BBB involves a multi-step cascade with a predominant role of α4β1-integrins. In contrast to CD4 T cells, α4β1-integrin mediated CD8 T cell interaction with the BBB was proposed to involve the tight junction protein junctional adhesion molecule (JAM)-B. Here, we made use of ODC-OVA mice expressing ovalbumin as neo-self-antigen in oligodendrocytes that is solely visible to CD8 T cells, allowing to investigate CD8 T cell-mediated autoimmune neuroinflammation. We generated JAM-B-deficient ODC-OVA mice (ODC-OVA; JAM-BKO mice) and compared CD8 T cell mediated autoimmune neuroinflammation to their ODC-OVA; JAM-BWT littermates. ODC-OVA; JAM-BKO mice developed ameliorated clinical disease, which was associated with a marked reduction in CD8 T cell infiltration into the CNS parenchyma. Surprisingly, lack of JAM-B did not affect CD8 T cell arrest or extravasation in spinal cord microvessels but rather resulted in CD8 T cell accumulation in the subarachnoid space and perivascular spaces in ODC-OVA; JAM-BKO mice. Detection of Jam-2 RNA expression in cells other than BBB endothelial cells contributing to CNS barriers including astrocytes forming the glia limitans, Bergmann glial cells, meningeal fibroblasts and choroid plexus epithelial cells suggests that JAM-B may regulate CD8 T cell entry into the CNS at barriers other than the BBB, particularly at the glia limitans. Thus, targeting JAM-B could provide a therapeutic strategy for treating neuroinflammation without disrupting T cell-mediated immune surveillance in CNS border compartments.

Functional interaction of electrical coupling and H-current and its putative impact on inhibitory transmission

The flow of information within neural circuits depends on the communication between neurons, primarily taking place at chemical and electrical synapses. The coexistence of these two modalities of synaptic transmission and their dynamical interaction with voltage-gated membrane conductances enables a rich repertoire of complex functional operations. One such operation, coincidence detection, allows electrically coupled neurons to respond more strongly to simultaneous synaptic inputs than to temporally dispersed ones. Using the mesencephalic trigeminal (MesV) nucleus-a structure composed of large, somatically coupled neurons-as an experimental model, we first demonstrate that electrical coupling strength in the hyperpolarized voltage range is highly time-dependent due to the involvement of the IH current. We then show how this property influences the coincidence detection of hyperpolarizing signals. Specifically, simultaneous hyperpolarizing inputs induce larger membrane potential changes, resulting in stronger IH current activation. This, in turn, shortens the temporal window for coincidence detection. We propose that this phenomenon may be crucial for network dynamics in circuits of electrically coupled neurons that receive inhibitory synaptic inputs and express the IH current. In particular, molecular layer interneurons (MLIs) of the cerebellar cortex provide an ideal model for studying coincidence detection of inhibitory synaptic inputs, and how this operation is shaped by the voltage-dependent conductances like the IH current, potentially impacting on motor coordination and learning.

MicroRNA-138-5p suppresses excitatory synaptic strength at the cerebellar input layer

MicroRNAs are small, highly conserved non-coding RNAs that negatively regulate mRNA translation and stability. In the brain, miRNAs contribute to neuronal development, synaptogenesis, and synaptic plasticity. MicroRNA 138-5p (miR-138-5p) controls inhibitory synaptic transmission in the hippocampus and is highly expressed in cerebellar excitatory neurons. However, its specific role in cerebellar synaptic transmission remains unknown. Here, we investigated excitatory transmission in the cerebellum of mice expressing a sponge construct that sequesters endogenous miR-138-5p. Mossy fibre stimulation-evoked EPSCs in granule cells were ∼40% larger in miR-138-5p sponge mice compared to controls. Furthermore, we observed larger miniature EPSC amplitudes, suggesting an increased number of functional postsynaptic AMPA receptors. High-frequency train stimulation revealed enhanced short-term depression following miR-138-5p downregulation. Together with computational modelling, this suggests a negative regulation of presynaptic release probability. Overall, our results demonstrate that miR-138-5p suppresses synaptic strength through pre- and postsynaptic mechanisms, providing a potentially powerful mechanism for tuning excitatory synaptic input into the cerebellum. KEY POINTS: MicroRNAs are powerful regulators of mRNA translation and control key cell biological processes including synaptic transmission, but their role in regulating synaptic function in the cerebellum has remained elusive. In this study, we investigated how microRNA-138-5p (miR-138-5p) modulates excitatory transmission at adult murine cerebellar mossy fibre to granule cell synapses. Downregulation of miR-138-5p enhances excitatory synaptic strength at the cerebellar input layer and increases short-term depression. miR-138-5p exerts its regulatory function through both pre- and postsynaptic mechanisms by negatively regulating release probability at mossy fibre boutons, as well as functional AMPA receptor numbers in granule cells. These findings provide insights into the role of miR-138-5p in the cerebellum and expand our understanding of microRNA-dependent control of excitatory synaptic transmission and short-term plasticity.

A single-cell transcriptomic atlas of all cell types in the brain of 5xFAD Alzheimer mice in response to dietary inulin supplementation

BACKGROUND

Alzheimer's disease (AD) is a progressive neurodegenerative disease that is a major threat to the aging population. Due to lack of effective therapy, preventive treatments are important strategies to limit AD onset and progression, of which dietary regimes have been implicated as a key factor. Diet with high fiber content is known to have beneficial effects on cognitive decline in AD. However, a global survey on microbiome and brain cell dynamics in response to high fiber intake at single-cell resolution in AD mouse models is still missing.

RESULTS

Here, we show that dietary inulin supplementation synergized with AD progression to specifically increase the abundance of Akkermansia muciniphila in gut microbiome of 5 × Familial AD (FAD) mice. By performing single-nucleus RNA sequencing on different regions of the whole brain with three independent biological replicates, we reveal region-specific changes in the proportion of neuron, astrocyte, and granule cell subpopulations upon inulin supplementation in 5xFAD mice. In addition, we find that astrocytes have more pronounced region-specific diversity than microglia. Intriguingly, such dietary change reduces amyloid-β plaque burden and alleviates microgliosis in the forebrain region, without affecting the spatial learning and memory.

CONCLUSIONS

These results provide a comprehensive overview on the transcriptomic changes in individual cells of the entire mouse brain in response to high fiber intake and a resourceful foundation for future mechanistic studies on the influence of diet and gut microbiome on the brain during neurodegeneration.

Genomic evolution reshapes cell-type diversification in the amniote brain

Over 320 million years of evolution, amniotes have developed complex brains and cognition through largely unexplored genetic and gene expression mechanisms. We created a comprehensive single-cell atlas of over 1.3 million cells from the telencephalon and cerebellum of turtles, zebra finches, pigeons, mice, and macaques, employing single-cell resolution spatial transcriptomics to validate gene expression patterns across species. Our study identifies significant species-specific variations in cell types, highlighting their conservation and diversification in evolution. We found pronounced differences in telencephalon excitatory neurons (EXs) and cerebellar cell types between birds and mammals. Birds predominantly express SLC17A6 in EX, whereas mammals express SLC17A7 in the neocortex and SLC17A6 elsewhere, possibly due to loss of function of SLC17A7 in birds. Additionally, we identified a bird-specific Purkinje cell subtype (SVIL+), implicating the lysine-specific demethylase 11 (LSD1)/KDM1A pathway in learning and circadian rhythms and containing numerous positively selected genes, which suggests an evolutionary optimization of cerebellar functions for ecological and behavioral adaptation. Our findings elucidate the complex interplay between genetic evolution and environmental adaptation, underscoring the role of genetic diversification in the development of specialized cell types across amniotes.

DSCT: a novel deep-learning framework for rapid and accurate spatial transcriptomic cell typing

Unraveling complex cell-type-composition and gene-expression patterns at the cellular spatial resolution is crucial for understanding intricate cell functions in the brain. In this study, we developed Deep Neural Network-based Spatial Cell Typing (DSCT)-an innovative framework for spatial cell typing within spatial transcriptomic data sets. This approach utilizes a synergistic integration of an enhanced gene-selection strategy and a lightweight deep neural network for data training, offering a more rapid and accurate solution for the analysis of spatial transcriptomic data. Based on comprehensive analysis, DSCT achieved exceptional accuracy in cell-type identification across various brain regions, species and spatial transcriptomic platforms. It also performed well in mapping finer cell types, thereby showcasing its versatility and adaptability across diverse data sets. Strikingly, DSCT exhibited high efficiency and remarkable processing speed, with fewer computational resource demands. As such, this novel approach opens new avenues for exploring the spatial organization of cell types and gene-expression patterns, advancing our understanding of biological functions and pathologies within the nervous system.

Spatial multi-omics reveals cell-type-specific nuclear compartments

The mammalian nucleus is compartmentalized by diverse subnuclear structures. These subnuclear structures, marked by nuclear bodies and histone modifications, are often cell-type specific and affect gene regulation and 3D genome organization1-3. Understanding their relationships rests on identifying the molecular constituents of subnuclear structures and mapping their associations with specific genomic loci and transcriptional levels in individual cells, all in complex tissues. Here, we introduce two-layer DNA seqFISH+, which enables simultaneous mapping of 100,049 genomic loci, together with the nascent transcriptome for 17,856 genes and subnuclear structures in single cells. These data enable imaging-based chromatin profiling of diverse subnuclear markers and can capture their changes at genomic scales ranging from 100-200 kilobases to approximately 1 megabase, depending on the marker and DNA locus. By using multi-omics datasets in the adult mouse cerebellum, we showed that repressive chromatin regions are more variable by cell type than are active regions across the genome. We also discovered that RNA polymerase II-enriched foci were locally associated with long, cell-type-specific genes (bigger than 200 kilobases) in a manner distinct from that of nuclear speckles. Furthermore, our analysis revealed that cell-type-specific regions of heterochromatin marked by histone H3 trimethylated at lysine 27 (H3K27me3) and histone H4 trimethylated at lysine 20 (H4K20me3) are enriched at specific genes and gene clusters, respectively, and shape radial chromosomal positioning and inter-chromosomal interactions in neurons and glial cells. Together, our results provide a single-cell high-resolution multi-omics view of subnuclear structures, associated genomic loci and their effects on gene regulation, directly within complex tissues.

Non-allometric expansion and enhanced compartmentalization of Purkinje cell dendrites in the human cerebellum

Purkinje cell (PC) dendrites are optimized to integrate the vast cerebellar input array and drive the sole cortical output. PCs are classically seen as stereotypical computational units, yet mouse PCs are morphologically diverse and those with multi-branched structure can receive non-canonical climbing fiber (CF) multi-innervation that confers independent compartment-specific signaling. While otherwise uncharacterized, human PCs are universally multi-branched. Do they exceed allometry to achieve enhanced integrative capacities relative to mouse PCs? To answer this, we used several comparative histology techniques in adult human and mouse to analyze cellular morphology, parallel fiber (PF) and CF input arrangement, and regional PC demographics. Human PCs are substantially larger than previously described; they exceed allometric constraint by cortical thickness and are the largest neuron in the brain with 6-7 cm total dendritic length. Unlike mouse, human PC dendrites ramify horizontally to form a multi-compartment motif that we show can receive multiple CFs. Human spines are denser (6.9 vs 4.9 spines/μm), larger (~0.36 vs 0.29 μm), and include an unreported 'spine cluster' structure-features that may be congruent with enhanced PF association and amplification as human-specific adaptations. By extrapolation, human PCs may receive 500,000 to 1 million synaptic inputs compared with 30-40,000 in mouse. Collectively, human PC morphology and input arrangement is quantitatively and qualitatively distinct from rodent. Multi-branched PCs are more prevalent in posterior and lateral cerebellum, co-varying with functional boundaries, supporting the hypothesis that this morphological motif permits expanded input multiplexing and may subserve task-dependent needs for input association.

Crosstalk Signaling Between the Epithelial and Non-Epithelial Compartments of the Mouse Inner Ear

PURPOSE

The otolith organs of the inner ear consist of the utricle and saccule that detect linear acceleration. These organs rely on mechanosensitive hair cells for transduction of signals to the central nervous system. In the murine utricle, about half of the hair cells are born during the first postnatal week. Here, we wanted to explore the role and interaction of the non-epithelial mesenchymal cells with the sensory epithelium and provide a resource for the auditory neurosciences community.

METHODS

We utilized full-length Smart-seq2 single-cell RNA sequencing at postnatal days 4 and 6 along with a host of computational methods to infer interactions between the epithelial and non-epithelial compartments of the mouse utricle. We validated these findings using a combination of immunohistochemistry and quantitative multiplex in situ hybridization.

RESULTS

We report diverse cell-cell crosstalk among the 12 annotated cell populations (n = 955 cells) in the developing neonatal mouse utricle, including epithelial and non-epithelial cellular signaling. The mesenchymal cells are the dominant signal senders during the postnatal period. Epithelial to mesenchymal signaling, as well as mesenchymal to epithelial signaling, are quantitatively shown through the TGFβ and pleiotrophin pathways.

CONCLUSION

This study highlights the dynamic process of postnatal vestibular organ development that relies not only on epithelial cells, but also on crosstalk between spatial compartments and among different cell groups. We further provide a data-rich resource for the inner ear community.

Deletion of the voltage-gated calcium channel gene, Ca V 1.3, reduces Purkinje cell dendritic complexity without altering cerebellar-mediated eyeblink conditioning

Genetic variation in CACNA1D , the gene that encodes the pore-forming subunit of the L-type calcium channel Ca V 1.3, has been associated with increased risk for neuropsychiatric disorders that display abnormalities in cerebellar structures. We sought to clarify if deletion of Ca V 1.3 in mice would induce abnormalities in cerebellar cortex cytoarchitecture or synapse morphology. Since Ca V 1.3 is highly expressed in cerebellar molecular layer interneurons (MLIs) and L-type channels appear to regulate GABA release from MLIs, we hypothesized that loss of Ca V 1.3 would alter GABAergic synapses between MLIs and Purkinje cells (PCs) without altering MLI numbers or PC structure. As expected, we did not observe changes in the numbers of MLIs or PCs. Surprisingly, Ca V 1.3 KO mice do have decreased complexity of PC dendritic arbors without differences in the number or structure of GABAergic synapses onto PCs. Loss of Ca V 1.3 was not associated with impaired acquisition of delay eyeblink conditioning. Therefore, our data suggest that Ca V 1.3 expression is important for PC structure but does not affect other measures of cerebellar cortex morphology or cerebellar function as assessed by delay eyeblink conditioning.

The cerebellar connectome

The cerebellum, once primarily associated with motor functions, has emerged as a critical component in higher cognitive processes and emotional regulation. This paradigm shift frames the cerebellum as an essential focal point for elucidating sophisticated functional brain circuitry. Network neuroscience often maintains a cortical-centric viewpoint, potentially overlooking the significant contributions of the cerebellum in connectome organization. Enhanced recognition and integration of cerebellar aspects in connectomic analyses hold significant potential for elucidating cerebellar circuitry within comprehensive brain networks and in neuropsychiatric conditions where cerebellar involvement is evident. This review explores the intricate anatomy, connectivity, and functional organization of the cerebellum within the broader context of large-scale brain networks. Cerebellar-specific networks are examined, emphasizing their role in supporting diverse cognitive functions via the cerebellum's hierarchical functional organization. The clinical significance of cerebellar connectomics is then addressed, highlighting the interplay between cerebellar circuitry and neurological and psychiatric conditions. The paper concludes by considering neurostimulation treatments and future directions in the field. This comprehensive review underscores the cerebellum's integral role in the human connectome.

Spatiotemporal diversity in molecular and functional abnormalities in the mdx dystrophic brain

Duchenne muscular dystrophy (DMD) is characterized by progressive muscle degeneration and neuropsychiatric abnormalities. Loss of full-length dystrophins is both necessary and sufficient to initiate DMD. These isoforms are expressed in the hippocampus, cerebral cortex (Dp427c), and cerebellar Purkinje cells (Dp427p). However, our understanding of the consequences of their absence, which is crucial for developing targeted interventions, remains inadequate. We combined RNA sequencing with genome-scale metabolic modelling (GSMM), immunodetection, and mitochondrial assays to investigate dystrophic alterations in the brains of the mdx mouse model of DMD. The cerebra and cerebella were analysed separately to discern the roles of Dp427c and Dp427p, respectively. Investigating these regions at 10 days (10d) and 10 weeks (10w) followed the evolution of abnormalities from development to early adulthood. These time points also encompass periods before onset and during muscle inflammation, enabling assessment of the potential damage caused by inflammatory mediators crossing the dystrophic blood-brain barrier. For the first time, we demonstrated that transcriptomic and functional dystrophic alterations are unique to the cerebra and cerebella and vary substantially between 10d and 10w. The common anomalies involved altered numbers of retained introns and spliced exons across mdx transcripts, corresponding with alterations in the mRNA processing pathways. Abnormalities in the cerebra were significantly more pronounced in younger mice. The top enriched pathways included those related to metabolism, mRNA processing, and neuronal development. GSMM indicated dysregulation of glucose metabolism, which corresponded with GLUT1 protein downregulation. The cerebellar dystrophic transcriptome, while significantly altered, showed an opposite trajectory to that of the cerebra, with few changes identified at 10 days. These late defects are specific and indicate an impact on the functional maturation of the cerebella that occurs postnatally. Although no classical neuroinflammation markers or microglial activation were detected at 10 weeks, specific differences indicate that inflammation impacts DMD brains. Importantly, some dystrophic alterations occur late and may therefore be amenable to therapeutic intervention, offering potential avenues for mitigating DMD-related neuropsychiatric defects.

Adaptive protein synthesis in genetic models of copper deficiency and childhood neurodegeneration

Rare inherited diseases caused by mutations in the copper transporters SLC31A1 (CTR1) or ATP7A induce copper deficiency in the brain, causing seizures and neurodegeneration in infancy through poorly understood mechanisms. Here, we used multiple model systems to characterize the molecular mechanisms by which neuronal cells respond to copper deficiency. Targeted deletion of CTR1 in neuroblastoma cells produced copper deficiency that produced a metabolic shift favoring glycolysis over oxidative phosphorylation. Proteomic and transcriptomic analysis of CTR1 knockout (KO) cells revealed simultaneous up-regulation of mTORC1 and S6K signaling and reduced PERK signaling. Patterns of gene and protein expression and pharmacogenomics show increased activation of the mTORC1-S6K pathway as a prosurvival mechanism, ultimately resulting in increased protein synthesis. Spatial transcriptomic profiling of Atp7aflx/Y :: Vil1Cre/+ mice identified up-regulated protein synthesis machinery and mTORC1-S6K pathway genes in copper-deficient Purkinje neurons in the cerebellum. Genetic epistasis experiments in Drosophila demonstrated that copper deficiency dendritic phenotypes in class IV neurons are improved or rescued by increased S6k expression or 4E-BP1 (Thor) RNAi, while epidermis phenotypes are exacerbated by Akt, S6k, or raptor RNAi. Overall, we demonstrate that increased mTORC1-S6K pathway activation and protein synthesis is an adaptive mechanism by which neuronal cells respond to copper deficiency.

A Gene-Expression Based Comparison of Murine and Human Inhibitory Interneurons in the Cerebellar Cortex and Nuclei

Cerebellar information processing is critically shaped by several types of inhibitory interneurons forming various intra-cerebellar feed-forward and feed-back loops. Evidence gathered over the past decades has focused interest on a non-uniform set of cortical inhibitory interneurons distinct from "classical" Golgi, basket or stellate cells, summarily referred to as PLIs (for Purkinje cell layer interneurons). Similarly, cerebellar nuclear inhibitory interneurons have gained increasing attention. Our understanding of the functions of these cells is still fragmentary. For humans, we lack functional data, and even any dependable morphological classification for these cells. Here, I used publicly available single cell based gene expression data to compare inhibitory interneurons from the cerebellar cortex and inhibitory nuclear neurons of humans and mice. Integration of nuclear and cortical cells revealed transcriptomic similarities between subsets of these cells and suggest known characteristics of cortical cell types may be helpful to devise strategies for the further characterization of nuclear inhibitory interneurons. Comparison of human and murine PLIs indicate that these strongly differ by the expression of genes used to characterize these cells in mice. This limits their utility to identify and classify human PLIs, and leaves the question open as to the number and characteristics of non-Golgi inhibitory interneurons resident in the cerebellar granule cell and Purkinje cell layers in humans.

The interplay between oxidative stress and inflammation supports autistic-related behaviors in Cntnap2 knockout mice

Autism Spectrum Disorder (ASD) is a highly prevalent neurodevelopmental condition characterized by social communication deficits and repetitive/restricted behaviors. Several studies showed that oxidative stress and inflammation may contribute to ASD. Indeed, increased levels of oxygen radicals and pro-inflammatory molecules were described in the brain and peripheral blood of persons with ASD and mouse models. Despite this, a potential direct connection between oxidative stress and inflammation within specific brain areas and ASD-related behaviors has not been investigated in detail yet. Here, we used RT-qPCR, RNA sequencing, metabolomics, immunohistochemistry, and flow cytometry to show that pro-inflammatory molecules were increased in the cerebellum and periphery of mice lacking Cntnap2, a robust model of ASD. In parallel, oxidative stress was present in the cerebellum of mutant animals. Systemic treatment with N-acetyl-cysteine (NAC) rescued cerebellar oxidative stress, inflammation, as well as motor and social impairments in Cntnap2-/- mice, concomitant with enhanced function of microglia cells in NAC-treated mutants. Intriguingly, social deficits, cerebellar inflammation, and microglia dysfunction were induced by NAC in Cntnap2+/+ animals. Our findings suggest that the interplay between oxidative stress and inflammation accompanied by genetic vulnerability may underlie ASD-related behaviors in Cntnap2 mutant mice.

The endocytic adaptor AP-2 maintains Purkinje cell function by balancing cerebellar parallel and climbing fiber synapses

The loss of cerebellar Purkinje cells is a hallmark of neurodegenerative movement disorders, but the mechanisms remain enigmatic. We show that endocytic adaptor protein complex 2 (AP-2) is crucial for Purkinje cell survival. Using mouse genetics, viral tracing, calcium imaging, and kinematic analysis, we demonstrate that loss of the AP-2 μ subunit in Purkinje cells leads to early-onset ataxia and progressive degeneration. Synaptic dysfunction, marked by an overrepresentation of parallel fibers (PFs) over climbing fibers (CFs), precedes Purkinje cell loss. Mechanistically, AP-2 interacts with the PF-enriched protein GRID2IP, and its loss triggers GRID2IP degradation and glutamate δ2 receptor (GLURδ2) accumulation, leading to an excess of PFs while CFs are reduced. The overrepresentation of PFs increases Purkinje cell network activity, which is mitigated by enhancing glutamate clearance with ceftriaxone. These findings highlight the role of AP-2 in regulating GRID2IP levels in Purkinje cells to maintain PF-CF synaptic balance and prevent motor dysfunction.

Neural circuit mechanisms to transform cerebellar population dynamics for motor control in monkeys

We exploit identification of neuron types during extracellular recording to demonstrate how the cerebellar cortex's well-established architecture transforms inputs into outputs. During smooth pursuit eye movements, the floccular complex performs distinct input-output transformations of temporal dynamics and directional response properties. The responses of different interneuron types localize the circuit mechanisms of each transformation. Mossy fibers and unipolar brush cells emphasize eye position dynamics uniformly across the cardinal axes; Purkinje cells and molecular layer interneurons code eye velocity along directionally biased axes; Golgi cells show unmodulated firing. Differential directional response properties of different neuron types localize the directional input-output transformation to the last-order inputs to Purkinje cells. Differential temporal dynamics pinpoint the site of the temporal input-output transformation to granule cells. Specific granule cell population dynamics allow the temporal transformations required in the area we study and generalize to many temporal transformations, providing a complete framework to understand cerebellar circuit computation.

Impact statement

We dissect the circuit computations performed by the floccular complex of the cerebellum during an exemplar sensory-motor behavior, taking advantage of knowledge of the circuit architecture, existence of discrete neuron types, and a newfound ability to identify neuron types from extracellular recordings. Our results describe the contributions of the major neuron types to the cerebellar input-output computations, identify the population dynamics needed in granule cells to support those computations, and to create a basis set to enable temporally-specific motor behavior and motor learning.

Non-allometric expansion and enhanced compartmentalization of Purkinje cell dendrites in the human cerebellum

Purkinje cell (PC) dendrites are optimized to integrate the vast cerebellar input array and drive the sole cortical output. PCs are classically seen as stereotypical computational units, yet mouse PCs are morphologically diverse and those with multi-branched structure can receive non-canonical climbing fiber (CF) multi-innervation that confers independent compartment-specific signaling. While otherwise uncharacterized, human PCs are universally multi-branched. Do they exceed allometry to achieve enhanced integrative capacities relative to mouse PCs? To answer this, we used several comparative histology techniques in adult human and mouse to analyze cellular morphology, parallel fiber (PF) and CF input arrangement, and regional PC demographics. Human PCs are substantially larger than previously described; they exceed allometric constraint by cortical thickness and are the largest neuron in the brain with 6-7cm total dendritic length. Unlike mouse, human PC dendrites ramify horizontally to form a multi-compartment motif that we show can receive multiple CFs. Human spines are denser (6.9 vs 4.9 spines/μm), larger (~0.36 vs 0.29μm), and include an unreported 'spine cluster' structure-features that may be congruent with enhanced PF association and amplification as human-specific adaptations. By extrapolation, human PCs may receive 500,000 to 1 million synaptic inputs compared with 30-40,000 in mouse. Collectively, human PC morphology and input arrangement is quantitatively and qualitatively distinct from rodent. Multi-branched PCs are more prevalent in posterior and lateral cerebellum, co-varying with functional boundaries, supporting the hypothesis that this morphological motif permits expanded input multiplexing and may subserve task-dependent needs for input association.

Spatial profiling of the interplay between cell type- and vision-dependent transcriptomic programs in the visual cortex

How early sensory experience during "critical periods" of postnatal life affects the organization of the mammalian neocortex at the resolution of neuronal cell types is poorly understood. We previously reported that the functional and molecular profiles of layer 2/3 (L2/3) cell types in the primary visual cortex (V1) are vision-dependent [S. Cheng et al., Cell 185, 311-327.e24 (2022)]. Here, we characterize the spatial organization of L2/3 cell types with and without visual experience. Spatial transcriptomic profiling based on 500 genes recapitulates the zonation of L2/3 cell types along the pial-ventricular axis in V1. By applying multitasking theory, we suggest that the spatial zonation of L2/3 cell types is linked to the continuous nature of their gene expression profiles, which can be represented as a 2D manifold bounded by three archetypal cell types. By comparing normally reared and dark reared L2/3 cells, we show that visual deprivation-induced transcriptomic changes comprise two independent gene programs. The first, induced specifically in the visual cortex, includes immediate-early genes and genes associated with metabolic processes. It manifests as a change in cell state that is orthogonal to cell-type-specific gene expression programs. By contrast, the second program impacts L2/3 cell-type identity, regulating a subset of cell-type-specific genes and shifting the distribution of cells within the L2/3 cell-type manifold. Through an integrated analysis of spatial transcriptomics with single-nucleus RNA-seq data, we describe how vision patterns cortical L2/3 cell types during the critical period.

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

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

Mapping the topography of spatial gene expression with interpretable deep learning

Spatially resolved transcriptomics technologies provide high-throughput measurements of gene expression in a tissue slice, but the sparsity of these data complicates analysis of spatial gene expression patterns. We address this issue by deriving a topographic map of a tissue slice-analogous to a map of elevation in a landscape-using a quantity called the isodepth. Contours of constant isodepths enclose domains with distinct cell type composition, while gradients of the isodepth indicate spatial directions of maximum change in expression. We develop GASTON (gradient analysis of spatial transcriptomics organization with neural networks), an unsupervised and interpretable deep learning algorithm that simultaneously learns the isodepth, spatial gradients and piecewise linear expression functions that model both continuous gradients and discontinuous variation in gene expression. We show that GASTON accurately identifies spatial domains and marker genes across several tissues, gradients of neuronal differentiation and firing in the brain, and gradients of metabolism and immune activity in the tumor microenvironment.

Transcript-specific enrichment enables profiling of rare cell states via single-cell RNA sequencing

Single-cell genomics technologies have accelerated our understanding of cell-state heterogeneity in diverse contexts. Although single-cell RNA sequencing identifies rare populations that express specific marker transcript combinations, traditional flow sorting requires cell surface markers with high-fidelity antibodies, limiting our ability to interrogate these populations. In addition, many single-cell studies require the isolation of nuclei from tissue, eliminating the ability to enrich learned rare cell states based on extranuclear protein markers. In the present report, we addressed these limitations by developing Programmable Enrichment via RNA FlowFISH by sequencing (PERFF-seq), a scalable assay that enables scRNA-seq profiling of subpopulations defined by the abundance of specific RNA transcripts. Across immune populations (n = 184,126 cells) and fresh-frozen and formalin-fixed, paraffin-embedded brain tissue (n = 33,145 nuclei), we demonstrated that programmable sorting logic via RNA-based cytometry can isolate rare cell populations and uncover phenotypic heterogeneity via downstream, high-throughput, single-cell genomics analyses.

Axonal pathology differentially affects human Purkinje cell subtypes in the essential tremor cerebellum

The cerebellar cortex is organized into discrete regions populated by molecularly distinct Purkinje cells (PCs), the sole cortical output neurons. While studies in animal models have shown that PC subtypes differ in their vulnerability to disease, our understanding of human PC subtype and vulnerability remains limited. Here, we demonstrate that human cerebellar regions specialized for motor vs cognitive functions (lobule HV vs Crus I) contain distinct PC populations characterized by specific molecular and anatomical features, which show selective vulnerability in essential tremor (ET), a cerebellar degenerative disorder. Using a known PC subtype marker, neurofilament heavy chain (NEFH), we found that motor lobule HV contains PCs with high NEFH expression, while cognitive lobule Crus I contains PCs with low NEFH expression in post-mortem samples from healthy controls. In the same cerebella, PC axons in lobule HV were 2.2-fold thicker than those in Crus I. Across lobules, axon caliber positively correlated with NEFH expression. In ET cerebella, we identified motor lobule-specific PC axon pathology with a 1.5-fold reduction in caliber and increased axon variability in lobule HV, while Crus I axons were unaffected. Tremor severity and duration in ET correlated with axon diameter variability selectively in lobule HV PCs. Given that axonal caliber is a major determinant of neural signaling capacity, our results (1) suggest that disrupted cerebellar corticonuclear signaling is occurring in ET, (2) provide evidence of region-specific PC subtypes in the human cerebellum and offer insight into how selective PC vulnerability may contribute to the pathophysiology of cerebellar degeneration.

Statistical identification of cell type-specific spatially variable genes in spatial transcriptomics

An essential task in spatial transcriptomics is identifying spatially variable genes (SVGs). Here, we present Celina, a statistical method for systematically detecting cell type-specific SVGs (ct-SVGs)-a subset of SVGs exhibiting distinct spatial expression patterns within specific cell types. Celina utilizes a spatially varying coefficient model to accurately capture each gene's spatial expression pattern in relation to the distribution of cell types across tissue locations, ensuring effective type I error control and high power. Celina proves powerful compared to existing methods in single-cell resolution spatial transcriptomics and stands as the only effective solution for spot-resolution spatial transcriptomics. Applied to five real datasets, Celina uncovers ct-SVGs associated with tumor progression and patient survival in lung cancer, identifies metagenes with unique spatial patterns linked to cell proliferation and immune response in kidney cancer, and detects genes preferentially expressed near amyloid-β plaques in an Alzheimer's model.

Tremor in the Age of Omics: An Overview of the Transcriptomic Landscape of Essential Tremor

Essential Tremor (ET) is the most common movement disorder and has a worldwide prevalence of 1%, including 5% of the population over 65 years old. It is characterized by an active, postural or kinetic tremor, primarily affecting the upper limbs, and is diagnosed based on clinical characteristics. The pathological mechanisms of ET, however, are mostly unknown. Moreover, despite its high heritability, genetic studies of ET genetics have yielded mixed results. Transcriptomics is a field that has the potential to reveal valuable insights about the processes and pathogenesis of ET thus providing an avenue for the development of more effective therapies. With the emergence of techniques such as single-cell and single-nucleus RNA sequencing (scRNA-seq and snRNA-seq), molecular and cellular events can now be more closely examined, providing valuable insights into potential causal mechanisms. In this review, we review the growing literature on transcriptomic studies in ET, aiming to identify biological pathways involved and explore possible avenues for further ET research. We emphasized the convergence on shared of biological pathways across several studies, specifically axonal guidance and calcium signaling. These findings posit multiple hypotheses linking both pathways through the regulation of axonal and synaptic plasticity. We conclude that increasing the sample size is vital to uncover the subtleties of ET clinical and pathological heterogeneity. Additionally, integrating Multiomics approaches should provide a comprehensive understanding of the disease's pathophysiology.

Realistic mossy fiber input patterns to unipolar brush cells evoke a continuum of temporal responses comprised of components mediated by different glutamate receptors

Unipolar brush cells (UBCs) are excitatory interneurons in the cerebellar cortex that receive mossy fiber (MF) inputs and excite granule cells. The UBC population responds to brief burst activation of MFs with a continuum of temporal transformations, but it is not known how UBCs transform the diverse range of MF input patterns that occur in vivo. Here, we use cell-attached recordings from UBCs in acute cerebellar slices to examine responses to MF firing patterns that are based on in vivo recordings. We find that MFs evoke a continuum of responses in the UBC population, mediated by three different types of glutamate receptors that each convey a specialized component. AMPARs transmit timing information for single stimuli at up to 5 spikes/s, and for very brief bursts. A combination of mGluR2/3s (inhibitory) and mGluR1s (excitatory) mediates a continuum of delayed, and broadened responses to longer bursts, and to sustained high frequency activation. Variability in the mGluR2/3 component controls the time course of the onset of firing, and variability in the mGluR1 component controls the duration of prolonged firing. We conclude that the combination of glutamate receptor types allows each UBC to simultaneously convey different aspects of MF firing. These findings establish that UBCs are highly flexible circuit elements that provide diverse temporal transformations that are well suited to contribute to specialized processing in different regions of the cerebellar cortex.

Molecular logic for cellular specializations that initiate the auditory parallel processing pathways

The cochlear nuclear complex (CN), the starting point for all central auditory processing, encompasses a suite of neuronal cell types highly specialized for neural coding of acoustic signals. However, the molecular logic governing these specializations remains unknown. By combining single-nucleus RNA sequencing and Patch-seq analysis, we reveal a set of transcriptionally distinct cell populations encompassing all previously observed types and discover multiple hitherto unknown subtypes with anatomical and physiological identity. The resulting comprehensive cell-type taxonomy reconciles anatomical position, morphological, physiological, and molecular criteria, enabling the determination of the molecular basis of the specialized cellular phenotypes in the CN. In particular, CN cell-type identity is encoded in a transcriptional architecture that orchestrates functionally congruent expression across a small set of gene families to customize projection patterns, input-output synaptic communication, and biophysical features required for encoding distinct aspects of acoustic signals. This high-resolution account of cellular heterogeneity from the molecular to the circuit level reveals the molecular logic driving cellular specializations, thus enabling the genetic dissection of auditory processing and hearing disorders with a high specificity.

CSI-GEP: A GPU-based unsupervised machine learning approach for recovering gene expression programs in atlas-scale single-cell RNA-seq data

Exploratory analysis of single-cell RNA sequencing (scRNA-seq) typically relies on hard clustering over two-dimensional projections like uniform manifold approximation and projection (UMAP). However, such methods can severely distort the data and have many arbitrary parameter choices. Methods that can model scRNA-seq data as non-discrete "gene expression programs" (GEPs) can better preserve the data's structure, but currently, they are often not scalable, not consistent across repeated runs, and lack an established method for choosing key parameters. Here, we developed a GPU-based unsupervised learning approach, "consensus and scalable inference of gene expression programs" (CSI-GEP). We show that CSI-GEP can recover ground truth GEPs in real and simulated atlas-scale scRNA-seq datasets, significantly outperforming cutting-edge methods, including GPT-based neural networks. We applied CSI-GEP to a whole mouse brain atlas of 2.2 million cells, disentangling endothelial cell types missed by other methods, and to an integrated scRNA-seq atlas of human tumors and cell lines, discovering mesenchymal-like GEPs unique to cancer cells growing in culture.

Astrocytes in aging

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

Arrested development: the dysfunctional life history of medulloblastoma

Medulloblastoma is a heterogeneous embryonal tumor of the cerebellum comprised of four distinct molecular subgroups that differ in their developmental origins, genomic landscapes, clinical presentation, and survival. Recent characterization of the human fetal cerebellum at single-cell resolution has propelled unprecedented insights into the cellular origins of medulloblastoma subgroups, including those underlying previously elusive groups 3 and 4. In this review, the molecular pathogenesis of medulloblastoma is examined through the lens of cerebellar development. In addition, we discuss how enhanced understanding of medulloblastoma origins has the potential to refine disease modeling for the advancement of treatment and outcomes.

Ndufs4 knockout mice with isolated complex I deficiency engage a futile adaptive brain response

Paediatric Leigh syndrome (LS) is an early-onset and fatal neurodegenerative disorder lacking treatment options. LS is frequently caused by mutations in the NDUFS4 gene, encoding an accessory subunit of mitochondrial complex I (CI), the first complex of the oxidative phosphorylation (OXPHOS) system. Whole-body Ndufs4 knockout (KO) mice (WB-KO mice) are widely used to study isolated CI deficiency, LS pathology and interventions. These animals develop a brain-specific phenotype via an incompletely understood pathomechanism. Here we performed a quantitative analysis of the sub-brain proteome in six-weeks old WB-KO mice vs. wildtype (WT) mice. Brain regions comprised of a brain slice (BrSl), cerebellum (CB), cerebral cortex (CC), hippocampus (HC), inferior colliculus (IC), and superior colliculus (SC). Proteome analysis demonstrated similarities between CC/HC, and between IC/SC, whereas BrSl and CB differed from these two groups and each other. All brain regions displayed greatly reduced levels of two CI structural subunits (NDUFS4, NDUFA12) and an increased level of the CI assembly factor NDUFAF2. The level of CI-Q module subunits was significantly more reduced in IC/SC than in BrSl/CB/CC/HC, whereas other OXPHOS complex levels were not reduced. Gene ontology and pathway analysis demonstrated specific and common proteome changes between brain regions. Across brain regions, upregulation of cold-shock-associated proteins, mitochondrial fatty acid (FA) oxidation and synthesis (mtFAS) were the most prominent. FA-related pathways were predominantly upregulated in CB and HC. Based upon these results, we argue that stimulation of these pathways is futile and pro-pathological and discuss alternative strategies for therapeutic intervention in LS. SIGNIFICANCE: The Ndufs4 knockout mouse model is currently the most relevant and most widely used animal model to study the brain-linked pathophysiology of human Leigh Syndrome (LS) and intervention strategies. We demonstrate that the Ndufs4 knockout brain engages futile and pro-pathological responses. These responses explain both negative and positive outcomes of intervention studies in Leigh Syndrome mice and patients, thereby guiding novel intervention opportunities.

Astrocytic Neuroligin-3 influences gene expression and social behavior, but is dispensable for synapse number

Neuroligin-3 (Nlgn3) is an autism-associated cell-adhesion molecule that interacts with neurexins and is robustly expressed in both neurons and astrocytes. Neuronal Nlgn3 is an essential regulator of synaptic transmission but the function of astrocytic Nlgn3 is largely unknown. Given the high penetrance of Nlgn3 mutations in autism and the emerging role of astrocytes in neuropsychiatric disorders, we here asked whether astrocytic Nlgn3 might shape neural circuit properties in the cerebellum similar to neuronal Nlgn3. Imaging of tagged Nlgn3 protein produced by CRISPR/Cas9-mediated genome editing showed that Nlgn3 is enriched in the cell body but not the fine processes of cerebellar astrocytes (Bergmann glia). Astrocyte-specific knockout of Nlgn3 did not detectably alter the number of synapses, synaptic transmission, or astrocyte morphology in mouse cerebellum. However, spatial transcriptomic analyses revealed a significant shift in gene expression among multiple cerebellar cell types after the deletion of astrocytic Nlgn3. Hence, in contrast to neuronal Nlgn3, astrocytic Nlgn3 in the cerebellum is not involved in shaping synapses but may modulate gene expression in specific brain areas.

Protocol for the generation of single-nuclei RNA-seq libraries and quantification of heterogeneous cell types activated during social interaction

Quantifying immediate early gene expression as a marker of cellular activity in single-nuclei RNA sequencing data allows for the identification of neurons involved in specific behaviors. Here, we present a protocol for generating single-nuclei libraries from the mouse brain and identifying active cell populations following social interactions. We describe steps for the dissection, preparation, and analysis of the prefrontal cortex, hippocampus, and cerebellum. This protocol has the potential to be modified for any brain region or behavior of interest. For complete details on the use and execution of this protocol, please refer to Walker and Frost.1.

TLX3 regulates CGN progenitor proliferation during cerebellum development and its dysfunction can lead to autism

Tlx3, a master regulator of the fate specification of excitatory neurons, is primarily known to function in post-mitotic cells. Although we have previously identified TLX3 expression in the proliferating granule neuron progenitors (GNPs) of cerebellum, its primary role is unknown. Here, we demonstrate that the dysfunction of Tlx3 from the GNPs significantly reduced its proliferation through regulating anti-proliferative genes. Consequently, the altered generation of GNPs resulted in cerebellar hypoplasia, patterning defects, granule neuron-Purkinje ratio imbalance, and aberrant synaptic connections in the cerebellum. This altered cerebellar homeostasis manifested into a typical autism-like behavior in mice with motor, and social function disabilities. We also show the presence of TLX3 variants with uncharacterized mutations in human cases of autism spectrum disorder (ASD). Altogether, our study establishes Tlx3 as a critical gene involved in developing GNPs and that its deletion from the early developmental stage culminates in autism.

Balance Performance in Aged Mice is Dependent on Unipolar Brush Cells

The vestibular processing regions of the cerebellum integrate vestibular information with other sensory modalities and motor signals to regulate balance, gaze stability, and spatial orientation. A class of excitatory glutamatergic interneurons known as unipolar brush cells (UBCs) are highly concentrated within the granule cell layer of these regions. UBCs receive vestibular signals directly from primary vestibular afferents and indirectly from mossy fibers. Each UBC excites numerous granule cells and could contribute to computations necessary for balance-related motor function. Prior research has implicated UBCs in motor function, but their influence on balance performance remains unclear, especially in aged mice that have age-related impairment. Here we tested whether UBCs contribute to motor coordination and balance by disrupting their activity with chemogenetics in aged and young mice. Age-related balance deficits were apparent in mice > 6 months old. Disrupting the activity of a subpopulation of UBCs caused aged mice to fall off a balance beam more frequently and altered swimming behaviors that are sensitive to vestibular dysfunction. These effects were not seen in young (7-week-old) mice. Thus, disrupting the activity of UBCs impairs mice with age-related balance issues and suggest that UBCs are essential for balance and vestibular function in aged mice.

Spall: accurate and robust unveiling cellular landscapes from spatially resolved transcriptomics data using a decomposition network

Recent developments in spatially resolved transcriptomics (SRT) enable the characterization of spatial structures for different tissues. Many decomposition methods have been proposed to depict the cellular distribution within tissues. However, existing computational methods struggle to balance spatial continuity in cell distribution with the preservation of cell-specific characteristics. To address this, we propose Spall, a novel decomposition network that integrates scRNA-seq data with SRT data to accurately infer cell type proportions. Spall introduced the GATv2 module, featuring a flexible dynamic attention mechanism to capture relationships between spots. This improves the identification of cellular distribution patterns in spatial analysis. Additionally, Spall incorporates skip connections to address the loss of cell-specific information, thereby enhancing the prediction capability for rare cell types. Experimental results show that Spall outperforms the state-of-the-art methods in reconstructing cell distribution patterns on multiple datasets. Notably, Spall reveals tumor heterogeneity in human pancreatic ductal adenocarcinoma samples and delineates complex tissue structures, such as the laminar organization of the mouse cerebral cortex and the mouse cerebellum. These findings highlight the ability of Spall to provide reliable low-dimensional embeddings for downstream analyses, offering new opportunities for deciphering tissue structures.

Pharmacological blocking of neutrophil extracellular traps attenuates immunothrombosis and neuroinflammation in cerebral cavernous malformation

Cerebral cavernous malformation (CCM) is a neurovascular disease with symptoms such as strokes, hemorrhages and neurological deficits. With surgery being the only treatment strategy, understanding the molecular mechanisms of CCM is crucial in finding alternative therapeutic options for CCM. Neutrophil extracellular traps (NETs) were recently reported in CCM, and NETs were shown to have positive or negative effects in different disease contexts. In this study, we investigated the roles of NETs in CCM by pharmacologically inhibiting NET formation using Cl-amidine (a peptidyl arginine deiminase inhibitor). We show here that Cl-amidine treatment reduced lesion burden, coagulation and endothelial-to-mesenchymal transition. Furthermore, NETs promoted the activation of microglia and fibroblasts, leading to increased neuroinflammation and a chronic wound microenvironment in CCM. The inhibition of NET formation caused endothelial quiescence and promoted a healthier microenvironment. Our study suggests the inhibition of NETs as a potential therapeutic strategy in CCM.

Neuronal diversity and stereotypy at multiple scales through whole brain morphometry

We conducted a large-scale whole-brain morphometry study by analyzing 3.7 peta-voxels of mouse brain images at the single-cell resolution, producing one of the largest multi-morphometry databases of mammalian brains to date. We registered 204 mouse brains of three major imaging modalities to the Allen Common Coordinate Framework (CCF) atlas, annotated 182,497 neuronal cell bodies, modeled 15,441 dendritic microenvironments, characterized the full morphology of 1876 neurons along with their axonal motifs, and detected 2.63 million axonal varicosities that indicate potential synaptic sites. Our analyzed six levels of information related to neuronal populations, dendritic microenvironments, single-cell full morphology, dendritic and axonal arborization, axonal varicosities, and sub-neuronal structural motifs, along with a quantification of the diversity and stereotypy of patterns at each level. This integrative study provides key anatomical descriptions of neurons and their types across a multiple scales and features, contributing a substantial resource for understanding neuronal diversity in mammalian brains.

Realistic mossy fiber input patterns to unipolar brush cells evoke a continuum of temporal responses comprised of components mediated by different glutamate receptors

Unipolar brush cells (UBCs) are excitatory interneurons in the cerebellar cortex that receive mossy fiber (MF) inputs and excite granule cells. The UBC population responds to brief burst activation of MFs with a continuum of temporal transformations, but it is not known how UBCs transform the diverse range of MF input patterns that occur in vivo. Here we use cell-attached recordings from UBCs in acute cerebellar slices to examine responses to MF firing patterns that are based on in vivo recordings. We find that MFs evoke a continuum of responses in the UBC population, mediated by three different types of glutamate receptors that each convey a specialized component. AMPARs transmit timing information for single stimuli at up to 5 spikes/s, and for very brief bursts. A combination of mGluR2/3s (inhibitory) and mGluR1s (excitatory) mediates a continuum of delayed, and broadened responses to longer bursts, and to sustained high frequency activation. Variability in the mGluR2/3 component controls the time course of the onset of firing, and variability in the mGluR1 component controls the duration of prolonged firing. We conclude that the combination of glutamate receptor types allows each UBC to simultaneously convey different aspects of MF firing. These findings establish that UBCs are highly flexible circuit elements that provide diverse temporal transformations that are well suited to contribute to specialized processing in different regions of the cerebellar cortex.

Understanding isoform expression by pairing long-read sequencing with single-cell and spatial transcriptomics

RNA isoform diversity, produced via alternative splicing, and alternative usage of transcription start and poly(A) sites, results in varied transcripts being derived from the same gene. Distinct isoforms can play important biological roles, including by changing the sequences or expression levels of protein products. The first single-cell approaches to RNA sequencing-and later, spatial approaches-which are now widely used for the identification of differentially expressed genes, rely on short reads and offer the ability to transcriptomically compare different cell types but are limited in their ability to measure differential isoform expression. More recently, long-read sequencing methods have been combined with single-cell and spatial technologies in order to characterize isoform expression. In this review, we provide an overview of the emergence of single-cell and spatial long-read sequencing and discuss the challenges associated with the implementation of these technologies and interpretation of these data. We discuss the opportunities they offer for understanding the relationships between the distinct variable elements of transcript molecules and highlight some of the ways in which they have been used to characterize isoforms' roles in development and pathology. Single-nucleus long-read sequencing, a special case of the single-cell approach, is also discussed. We attempt to cover both the limitations of these technologies and their significant potential for expanding our still-limited understanding of the biological roles of RNA isoforms.

Impact of Intrauterine Insults on Fetal and Postnatal Cerebellar Development in Humans and Rodents

Many children suffer from neurodevelopmental aberrations that have long-term effects. To understand the consequences of pathological processes during particular periods in neurodevelopment, one has to understand the differences in the developmental timelines of brain regions. The cerebellum is one of the first brain structures to differentiate during development but one of the last to achieve maturity. This relatively long period of development underscores its vulnerability to detrimental environmental exposures throughout gestation. Moreover, as postnatal functionality of the cerebellum is multifaceted, enveloping sensorimotor, cognitive, and emotional domains, prenatal disruptions in cerebellar development can result in a large variety of neurological and mental health disorders. Here, we review major intrauterine insults that affect cerebellar development in both humans and rodents, ranging from abuse of toxic chemical agents, such as alcohol, nicotine, cannabis, and opioids, to stress, malnutrition, and infections. Understanding these pathological mechanisms in the context of the different stages of cerebellar development in humans and rodents can help us to identify critical and vulnerable periods and thereby prevent the risk of associated prenatal and early postnatal damage that can lead to lifelong neurological and cognitive disabilities. The aim of the review is to raise awareness and to provide information for obstetricians and other healthcare professionals to eventually design strategies for preventing or rescuing related neurodevelopmental disorders.

Adaptive protein synthesis in genetic models of copper deficiency and childhood neurodegeneration

Rare inherited diseases caused by mutations in the copper transporters SLC31A1 (CTR1) or ATP7A induce copper deficiency in the brain, causing seizures and neurodegeneration in infancy through poorly understood mechanisms. Here, we used multiple model systems to characterize the molecular mechanisms by which neuronal cells respond to copper deficiency. Targeted deletion of CTR1 in neuroblastoma cells produced copper deficiency that was associated with a metabolic shift favoring glycolysis over oxidative phosphorylation. Proteomic and transcriptomic analysis of CTR1 KO cells revealed simultaneous upregulation of mTORC1 and S6K signaling and reduced PERK signaling. Patterns of gene and protein expression and pharmacogenomics show increased activation of the mTORC1-S6K pathway as a pro-survival mechanism, ultimately resulting in increased protein synthesis. Spatial transcriptomic profiling of Atp7a flx/Y :: Vil1 Cre/+ mice identified upregulated protein synthesis machinery and mTORC1-S6K pathway genes in copper-deficient Purkinje neurons in the cerebellum. Genetic epistasis experiments in Drosophila demonstrated that copper deficiency dendritic phenotypes in class IV neurons are partially rescued by increased S6k expression or 4E-BP1 (Thor) RNAi, while epidermis phenotypes are exacerbated by Akt, S6k, or raptor RNAi. Overall, we demonstrate that increased mTORC1-S6K pathway activation and protein synthesis is an adaptive mechanism by which neuronal cells respond to copper deficiency.

Tissue nonspecific alkaline phosphatase deficiency impairs Purkinje cell development and survival in a mouse model of infantile hypophosphatasia

Loss-of-function mutations in the tissue-nonspecific alkaline phosphatase (TNAP) gene can result in hypophosphatasia (HPP), an inherited multi-systemic metabolic disorder that is well-known for skeletal and dental hypomineralization. However, emerging evidence shows that both adult and pediatric patients with HPP suffer from cognitive deficits, higher measures of depression and anxiety, and impaired sensorimotor skills. The cerebellum plays an important role in sensorimotor coordination, cognition, and emotion. To date, the impact of TNAP mutation on the cerebellar circuitry development and function remains poorly understood. The main objective of this study was to investigate the roles of TNAP in cerebellar development and function, with a particular focus on Purkinje cells, in a mouse model of infantile HPP. Male and female wild type (WT) and TNAP knockout (KO) mice underwent behavioral testing on postnatal day 13-14 and were euthanized after completion of behavioral tests. Cerebellar tissues were harvested for gene expression and immunohistochemistry analyses. We found that TNAP mutation resulted in significantly reduced body weight, shorter body length, and impaired sensorimotor functions in both male and female KO mice. These developmental and behavioral deficits were accompanied by abnormal Purkinje cell morphology and dysregulation of genes that regulates synaptic transmission, cellular growth, proliferation, and death. In conclusion, inactivation of TNAP via gene deletion causes developmental delays, sensorimotor impairment, and Purkinje cell maldevelopment. These results shed light on a new perspective of cerebellar dysfunction in HPP.

FOXP Genes Regulate Purkinje Cell Diversity in Cerebellar Development and Evolution

Mammalian cerebellar development is thought to be influenced by distinct Purkinje cell (PC) subtypes. However, the degree of PC heterogeneity and the molecular drivers of this diversity have remained unclear, hindering efforts to manipulate PC diversification and assess its role in cerebellar development. Here, we demonstrate the critical role of Foxp genes in cerebellar development by regulating PC diversification. We identified 11 PC subtypes in the embryonic mouse cerebellum through single-cell RNA sequencing. Using a novel unsupervised method, we mapped these subtypes in three-dimensional space, revealing discrete PC subtypes predictive of adult cerebellar organization, including longitudinal stripes and lobules. These subtypes exhibit unique combinations of Foxp1 , Foxp2 , and Foxp4 expression. Deletion of Foxp2 and Foxp1 disrupts PC diversification, leading to altered cerebellar patterning, including the loss of a specific Foxp1 -expressing subtype and the cerebellar hemisphere. The Foxp1 -expressing PC subtype is much more abundant in the fetal human cerebellum than in mice, but rare in the chick cerebellum, correlating with cerebellar hemisphere size in these species. This highlights the significance of Foxp1 -expressing PCs in cerebellar hemisphere development and evolution. Therefore, our study identifies Foxp genes as key regulators of PC diversity, providing new insights into the developmental and evolutionary underpinnings of the cerebellum.

Dysregulation of zebrin-II cell subtypes in the cerebellum is a shared feature across polyglutamine ataxia mouse models and patients

Spinocerebellar ataxia type 7 (SCA7) is a genetic neurodegenerative disorder caused by a CAG-polyglutamine repeat expansion. Purkinje cells (PCs) are central to the pathology of ataxias, but their low abundance in the cerebellum underrepresents their transcriptomes in sequencing assays. To address this issue, we developed a PC enrichment protocol and sequenced individual nuclei from mice and patients with SCA7. Single-nucleus RNA sequencing in SCA7-266Q mice revealed dysregulation of cell identity genes affecting glia and PCs. Specifically, genes marking zebrin-II PC subtypes accounted for the highest proportion of DEGs in symptomatic SCA7-266Q mice. These transcriptomic changes in SCA7-266Q mice were associated with increased numbers of inhibitory synapses as quantified by immunohistochemistry and reduced spiking of PCs in acute brain slices. Dysregulation of zebrin-II cell subtypes was the predominant signal in PCs of SCA7-266Q mice and was associated with the loss of zebrin-II striping in the cerebellum at motor symptom onset. We furthermore demonstrated zebrin-II stripe degradation in additional mouse models of polyglutamine ataxia and observed decreased zebrin-II expression in the cerebella of patients with SCA7. Our results suggest that a breakdown of zebrin subtype regulation is a shared pathological feature of polyglutamine ataxias.

Structural deviations of the posterior fossa and the cerebellum and their cognitive links in a neurodevelopmental deletion syndrome

High-impact genetic variants associated with neurodevelopmental disorders provide biologically-defined entry points for mechanistic investigation. The 3q29 deletion (3q29Del) is one such variant, conferring a 40-100-fold increased risk for schizophrenia, as well as high risk for autism and intellectual disability. However, the mechanisms leading to neurodevelopmental disability remain largely unknown. Here, we report the first in vivo quantitative neuroimaging study in individuals with 3q29Del (N = 24) and neurotypical controls (N = 1608) using structural MRI. Given prior radiology reports of posterior fossa abnormalities in 3q29Del, we focused our investigation on the cerebellum and its tissue-types and lobules. Additionally, we compared the prevalence of cystic/cyst-like malformations of the posterior fossa between 3q29Del and controls and examined the association between neuroanatomical findings and quantitative traits to probe gene-brain-behavior relationships. 3q29Del participants had smaller cerebellar cortex volumes than controls, before and after correction for intracranial volume (ICV). An anterior-posterior gradient emerged in finer grained lobule-based and voxel-wise analyses. 3q29Del participants also had larger cerebellar white matter volumes than controls following ICV-correction and displayed elevated rates of posterior fossa arachnoid cysts and mega cisterna magna findings independent of cerebellar volume. Cerebellar white matter and subregional gray matter volumes were associated with visual-perception and visual-motor integration skills as well as IQ, while cystic/cyst-like malformations yielded no behavioral link. In summary, we find that abnormal development of cerebellar structures may represent neuroimaging-based biomarkers of cognitive and sensorimotor function in 3q29Del, adding to the growing evidence identifying cerebellar pathology as an intersection point between syndromic and idiopathic forms of neurodevelopmental disabilities.

From Noise to Knowledge: Diffusion Probabilistic Model-Based Neural Inference of Gene Regulatory Networks

Understanding gene regulatory networks (GRNs) is crucial for elucidating cellular mechanisms and advancing therapeutic interventions. Original methods for GRN inference from bulk expression data often struggled with the high dimensionality and inherent noise in the data. Here we introduce RegDiffusion, a new class of Denoising Diffusion Probabilistic Models focusing on the regulatory effects among feature variables. RegDiffusion introduces Gaussian noise to the input data following a diffusion schedule and uses a neural network with a parameterized adjacency matrix to predict the added noise. We show that using this process, GRNs can be learned effectively with a surprisingly simple model architecture. In our benchmark experiments, RegDiffusion shows superior performance compared to several baseline methods in multiple datasets. We also demonstrate that RegDiffusion can infer biologically meaningful regulatory networks from real-world single-cell data sets with over 15,000 genes in under 5 minutes. This work not only introduces a fresh perspective on GRN inference but also highlights the promising capacity of diffusion-based models in the area of single-cell analysis. The RegDiffusion software package and experiment data are available at https://github.com/TuftsBCB/RegDiffusion.

Increased understanding of complex neuronal circuits in the cerebellar cortex

The prevailing belief has been that the fundamental structures of cerebellar neuronal circuits, consisting of a few major neuron types, are simple and well understood. Given that the cerebellum has long been known to be crucial for motor behaviors, these simple yet organized circuit structures seemed beneficial for theoretical studies proposing neural mechanisms underlying cerebellar motor functions and learning. On the other hand, experimental studies using advanced techniques have revealed numerous structural properties that were not traditionally defined. These include subdivided neuronal types and their circuit structures, feedback pathways from output Purkinje cells, and the multidimensional organization of neuronal interactions. With the recent recognition of the cerebellar involvement in non-motor functions, it is possible that these newly identified structural properties, which are potentially capable of generating greater complexity than previously recognized, are associated with increased information capacity. This, in turn, could contribute to the wide range of cerebellar functions. However, it remains largely unknown how such structural properties contribute to cerebellar neural computations through the regulation of neuronal activity or synaptic transmissions. To promote further research into cerebellar circuit structures and their functional significance, we aim to summarize the newly identified structural properties of the cerebellar cortex and discuss future research directions concerning cerebellar circuit structures and their potential functions.

Spatial profiling of the interplay between cell type- and vision-dependent transcriptomic programs in the visual cortex

How early sensory experience during "critical periods" of postnatal life affects the organization of the mammalian neocortex at the resolution of neuronal cell types is poorly understood. We previously reported that the functional and molecular profiles of layer 2/3 (L2/3) cell types in the primary visual cortex (V1) are vision-dependent (Tan et al., Neuron, 108(4), 2020; Cheng et al., Cell, 185(2), 2022). Here, we characterize the spatial organization of L2/3 cell types with and without visual experience. Spatial transcriptomic profiling based on 500 genes recapitulates the zonation of L2/3 cell types along the pial-ventricular axis in V1. By applying multi-tasking theory (Adler et al., Cell Systems, 8, 2019), we suggest that the spatial zonation of L2/3 cell types is linked to the continuous nature of their gene expression profiles, which can be represented as a 2D manifold bounded by three archetypal cell types ("A", "B", and "C"). By comparing normally reared and dark reared L2/3 cells, we show that visual deprivation-induced transcriptomic changes comprise two independent gene programs. The first, induced specifically in the visual cortex, includes immediate-early genes and genes associated with metabolic processes. It manifests as a change in cell state that is orthogonal to cell type-specific gene expression programs. By contrast, the second program impacts L2/3 cell type identity, regulating a subset of cell type-specific genes and shifting the distribution of cells within the L2/3 manifold, with a depression of the B-type and C-type and a gain of the A-type. Through an integrated analysis of spatial transcriptomic measurements with single-nucleus RNA-seq data from our previous study, we describe how vision patterns L2/3 cortical cell types during the postnatal critical period.

Significance statement

Layer 2/3 (L2/3) glutamatergic neurons are important sites of experience-dependent plasticity and learning in the mammalian cortex. Their properties vary continuously with cortical depth and depend upon experience. Here, by applying spatial transcriptomics and different computational approaches in the mouse primary visual cortex, we show that vision regulates orthogonal gene expression programs underlying cell states and cell types. Visual deprivation not only induces an activity-dependent cell state, but also alters the composition of L2/3 cell types, which are appropriately described as a transcriptomic continuum. Our results provide insights into how experience shapes transcriptomes that may, in turn, sculpt brain wiring, function, and behavior.

Balance performance in aged mice is dependent on unipolar brush cells

The vestibular processing regions of the cerebellum integrate vestibular information with other sensory modalities and motor signals to regulate balance, gaze stability, and spatial orientation. A class of excitatory glutamatergic interneurons known as unipolar brush cells (UBCs) are highly concentrated within the granule cell layer of these regions. UBCs receive vestibular signals directly from primary vestibular afferents and indirectly from mossy fibers. Each UBC excites numerous granule cells and could contribute to computations necessary for balance-related motor function. Prior research has implicated UBCs in motor function, but their influence on balance performance remains unclear, especially in aged mice that have age-related impairment. Here we tested whether UBCs contribute to motor coordination and balance by disrupting their activity with chemogenetics in aged and young mice. Age-related balance deficits were apparent in mice > 6 months old. Disrupting the activity of a subpopulation of UBCs caused aged mice to fall off a balance beam more frequently and altered swimming behaviors that are sensitive to vestibular dysfunction. These effects were not seen in young (7-week-old) mice. Thus, disrupting the activity of UBCs impairs mice with age-related balance issues and suggest that UBCs are essential for balance and vestibular function in aged mice.