Shared and distinct transcriptomic cell types across neocortical areas

Expression and function of nicotinic acetylcholine receptors in specific neuronal populations: Focus on striatal and prefrontal circuits

Nicotinic acetylcholine receptors (nAChRs) are widely expressed in the central nervous system and play an important role in the control of neural functions including neuronal activity, transmitter release and synaptic plasticity. Although the common subtypes of nAChRs are abundantly expressed throughout the brain, their expression in different brain regions and by individual neuronal types is not homogeneous or incidental. In recent years, several studies have emerged showing that particular subtypes of nAChRs are expressed by specific neuronal populations in which they have major influence on the activity of local circuits and behavior. It has been demonstrated that even nAChRs expressed by relatively rare neuronal types can induce significant changes in behavior and contribute to pathological processes. Depending on the identity and connectivity of the particular nAChRs-expressing neuronal populations, the activation of nAChRs can have distinct or even opposing effects on local neuronal signaling. In this review, we will summarize the available literature describing the expression of individual nicotinic subunits by different neuronal types in two crucial brain regions, the striatum and the prefrontal cortex. The review will also briefly discuss nicotinic expression in non-neuronal, glial cells, as they cannot be ignored as potential targets of nAChRs-modulating drugs. The final section will discuss options that could allow us to target nAChRs in a neuronal-type-specific manner, not only in the experimental field, but also eventually in clinical practice.

Chronic intermittent hypoxia elicits distinct transcriptomic responses among neurons and oligodendrocytes within the brainstem of mice

Chronic intermittent hypoxia (CIH) is a prevalent condition characterized by recurrent episodes of oxygen deprivation, linked to respiratory and neurological disorders. Prolonged CIH is known to have adverse effects, including endothelial dysfunction, chronic inflammation, oxidative stress, and impaired neuronal function. These factors can contribute to serious comorbidities, including metabolic disorders and cardiovascular diseases. To investigate the molecular impact of CIH, we examined male C57BL/6J mice exposed to CIH for 21 days, comparing with normoxic controls. We used single-nucleus RNA sequencing to comprehensively examine the transcriptomic impact of CIH on key cell classes within the brainstem, specifically excitatory neurons, inhibitory neurons, and oligodendrocytes. These cell classes regulate essential physiological functions, including autonomic tone, cardiovascular control, and respiration. Through analysis of 10,995 nuclei isolated from pontine-medullary tissue, we identified seven major cell classes, further subdivided into 24 clusters. Our findings among these cell classes, revealed significant differential gene expression, underscoring their distinct responses to CIH. Notably, neurons exhibited transcriptional dysregulation of genes associated with synaptic transmission, and structural remodeling. In addition, we found dysregulated genes encoding ion channels and inflammatory response. Concurrently, oligodendrocytes exhibited dysregulated genes associated with oxidative phosphorylation and oxidative stress. Utilizing CellChat network analysis, we uncovered CIH-dependent altered patterns of diffusible intercellular signaling. These insights offer a comprehensive transcriptomic cellular atlas of the pons-medulla and provide a fundamental resource for the analysis of molecular adaptations triggered by CIH.NEW & NOTEWORTHY This study on chronic intermittent hypoxia (CIH) from pons-medulla provides initial insights into the molecular effects on excitatory neurons, inhibitory neurons, and oligodendrocytes, highlighting our unbiased approach, in comparison with earlier studies focusing on single target genes. Our findings reveal that CIH affects cell classes distinctly, and the dysregulated genes in distinct cell classes are associated with synaptic transmission, ion channels, inflammation, oxidative stress, and intercellular signaling, advancing our understanding of CIH-induced molecular responses.

Simple synaptic modulations implement diverse novelty computations

Detecting novelty is ethologically useful for an organism's survival. Recent experiments characterize how different types of novelty over timescales from seconds to weeks are reflected in the activity of excitatory and inhibitory neuron types. Here, we introduce a learning mechanism, familiarity-modulated synapses (FMSs), consisting of multiplicative modulations dependent on presynaptic or pre/postsynaptic neuron activity. With FMSs, network responses that encode novelty emerge under unsupervised continual learning and minimal connectivity constraints. Implementing FMSs within an experimentally constrained model of a visual cortical circuit, we demonstrate the generalizability of FMSs by simultaneously fitting absolute, contextual, and omission novelty effects. Our model also reproduces functional diversity within cell subpopulations, leading to experimentally testable predictions about connectivity and synaptic dynamics that can produce both population-level novelty responses and heterogeneous individual neuron signals. Altogether, our findings demonstrate how simple plasticity mechanisms within a cortical circuit structure can produce qualitatively distinct and complex novelty responses.

Higher-order thalamocortical circuits are specified by embryonic cortical progenitor types in the mouse brain

The sensory cortex receives synaptic inputs from both first-order and higher-order thalamic nuclei. First-order inputs relay simple stimulus properties from the periphery, whereas higher-order inputs relay more complex response properties, provide contextual feedback, and modulate plasticity. Here, we reveal that a cortical neuron's higher-order input is determined by the type of progenitor from which it is derived during embryonic development. Within layer 4 (L4) of the mouse primary somatosensory cortex, neurons derived from intermediate progenitors receive stronger higher-order thalamic input and exhibit greater higher-order sensory responses. These effects result from differences in dendritic morphology and levels of the transcription factor Lhx2, which are specified by the L4 neuron's progenitor type. When this mechanism is disrupted, cortical circuits exhibit altered higher-order responses and sensory-evoked plasticity. Therefore, by following distinct trajectories, progenitor types generate diversity in thalamocortical circuitry and may provide a general mechanism for differentially routing information through the cortex.

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

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

Spatial transcriptomic brain imaging reveals the effects of immunomodulation therapy on specific regional brain cells in a mouse dementia model

Increasing evidence of brain-immune crosstalk raises expectations for the efficacy of novel immunotherapies in Alzheimer's disease (AD), but the lack of methods to examine brain tissues makes it difficult to evaluate therapeutics. Here, we investigated the changes in spatial transcriptomic signatures and brain cell types using the 10x Genomics Visium platform in immune-modulated AD models after various treatments. To proceed with an analysis suitable for barcode-based spatial transcriptomics, we first organized a workflow for segmentation of neuroanatomical regions, establishment of appropriate gene combinations, and comprehensive review of altered brain cell signatures. Ultimately, we investigated spatial transcriptomic changes following administration of immunomodulators, NK cell supplements and an anti-CD4 antibody, which ameliorated behavior impairment, and designated brain cells and regions showing probable associations with behavior changes. We provided the customized analytic pipeline into an application named STquantool. Thus, we anticipate that our approach can help researchers interpret the real action of drug candidates by simultaneously investigating the dynamics of all transcripts for the development of novel AD therapeutics.

Improving replicability in single-cell RNA-Seq cell type discovery with Dune

BACKGROUND

Single-cell transcriptome sequencing (scRNA-Seq) has allowed new types of investigations at unprecedented levels of resolution. Among the primary goals of scRNA-Seq is the classification of cells into distinct types. Many approaches build on existing clustering literature to develop tools specific to single-cell. However, almost all of these methods rely on heuristics or user-supplied parameters to control the number of clusters. This affects both the resolution of the clusters within the original dataset as well as their replicability across datasets. While many recommendations exist, in general, there is little assurance that any given set of parameters will represent an optimal choice in the trade-off between cluster resolution and replicability. For instance, another set of parameters may result in more clusters that are also more replicable.

RESULTS

Here, we propose Dune, a new method for optimizing the trade-off between the resolution of the clusters and their replicability. Our method takes as input a set of clustering results-or partitions-on a single dataset and iteratively merges clusters within each partitions in order to maximize their concordance between partitions. As demonstrated on multiple datasets from different platforms, Dune outperforms existing techniques, that rely on hierarchical merging for reducing the number of clusters, in terms of replicability of the resultant merged clusters as well as concordance with ground truth. Dune is available as an R package on Bioconductor: https://www.bioconductor.org/packages/release/bioc/html/Dune.html .

CONCLUSIONS

Cluster refinement by Dune helps improve the robustness of any clustering analysis and reduces the reliance on tuning parameters. This method provides an objective approach for borrowing information across multiple clusterings to generate replicable clusters most likely to represent common biological features across multiple datasets.

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

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

Integrative multi-omic analysis reveals conserved cell-projection deficits in human Down syndrome brains

Down syndrome (DS) is the most common genetic cause of cognitive disability. However, it is largely unclear how triplication of a small gene subset may impinge on diverse aspects of DS brain physiopathology. Here, we took a multi-omic approach and simultaneously analyzed by RNA-seq and proteomics the expression signatures of two diverse regions of human postmortem DS brains. We found that the overexpression of triplicated genes triggered global expression dysregulation, differentially affecting transcripts, miRNAs, and proteins involved in both known and novel biological candidate pathways. Among the latter, we observed an alteration in RNA splicing, specifically modulating the expression of genes involved in cytoskeleton and axonal dynamics in DS brains. Accordingly, we found an alteration in axonal polarization in neurons from DS human iPSCs and mice. Thus, our study provides an integrated multilayer expression database capable of identifying new potential targets to aid in designing future clinical interventions for DS.

Sex-specific single cell-level transcriptomic signatures of Rett syndrome disease progression

Dominant X-linked diseases are uncommon due to female X chromosome inactivation (XCI). While random XCI usually protects females against X-linked mutations, Rett syndrome (RTT) is a female neurodevelopmental disorder caused by heterozygous MECP2 mutation. After 6-18 months of typical neurodevelopment, RTT girls undergo poorly understood regression. We performed longitudinal snRNA-seq on cerebral cortex in a construct-relevant Mecp2e1 mutant mouse model of RTT, revealing transcriptional effects of cell type, mosaicism, and sex on progressive disease phenotypes. Across cell types, we observed sex differences in the number of differentially expressed genes (DEGs) with 6x more DEGs in mutant females than males. Unlike males, female DEGs emerged prior to symptoms, were enriched for homeostatic gene pathways in distinct cell types over time, and correlated with disease phenotypes and human RTT cortical cell transcriptomes. Non-cell-autonomous effects were prominent and dynamic across disease progression of Mecp2e1 mutant females, indicating wild-type-expressing cells normalizing transcriptional homeostasis. These results improve understanding of RTT progression and treatment.

Cell type-specific gene expression dynamics during human brain maturation

The human brain undergoes protracted post-natal maturation, guided by dynamic changes in gene expression. Most studies exploring these processes have used bulk tissue analyses, which mask cell type-specific gene expression dynamics. Here, using single nucleus (sn)RNA-seq on temporal lobe tissue, including samples of African ancestry, we build a joint paediatric and adult atlas of 75 cell subtypes, which we verify with spatial transcriptomics. We explore the differences between paediatric and adult cell types, revealing the genes and pathways that change during brain maturation. Our results highlight excitatory neuron subtypes, including the LTK and FREM subtypes, that show elevated expression of genes associated with cognition and synaptic plasticity in paediatric tissue. The new resources we present here improve our understanding of the brain during its development and contribute to global efforts to build an inclusive brain cell map.

Optogenetic Determination of Dynamic and Cell-Type-Specific Inhibitory Reversal Potentials

The reversal potential refers to the membrane potential at which the net current flow through a channel reverses direction. The reversal potential is determined by transmembrane ion gradients and, in turn, determines how the channel's activity will affect the membrane potential. Traditional investigation into the reversal potential of inhibitory ligand-gated ion channels (EInh) has relied upon the activation of endogenous receptors, such as the GABA-A receptor (GABAAR). There are, however, challenges associated with activating endogenous receptors, including agonist delivery, isolating channel responses, and the effects of receptor saturation and desensitization. Here, we demonstrate the utility of using a light-gated anion channel, stGtACR2, to probe EInh in the rodent brain. Using mice of both sexes, we demonstrate that the properties of this optically activated channel make it a suitable proxy for studying GABAAR receptor-mediated inhibition. We validate this agonist-independent optogenetic strategy in vitro and in vivo and further show how it can accurately capture differences in EInh dynamics following manipulations of endogenous ion fluxes. This allows us to explore distinct resting EInh differences across genetically defined neuronal subpopulations. Using this approach to challenge ion homeostasis mechanisms in neurons, we uncover cell-specific EInh dynamics that are supported by the differential expression of endogenous ion handling mechanisms. Our findings therefore establish an effective optical strategy for revealing novel aspects of inhibitory reversal potentials and thereby expand the repertoire of optogenetics.

Massively parallel in vivo Perturb-seq reveals cell-type-specific transcriptional networks in cortical development

Leveraging AAVs' versatile tropism and labeling capacity, we expanded the scale of in vivo CRISPR screening with single-cell transcriptomic phenotyping across embryonic to adult brains and peripheral nervous systems. Through extensive tests of 86 vectors across AAV serotypes combined with a transposon system, we substantially amplified labeling efficacy and accelerated in vivo gene delivery from weeks to days. Our proof-of-principle in utero screen identified the pleiotropic effects of Foxg1, highlighting its tight regulation of distinct networks essential for cell fate specification of Layer 6 corticothalamic neurons. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% by lentivirus, to achieve analysis of over 30,000 cells in one experiment and enable massively parallel in vivo Perturb-seq. Compatible with various phenotypic measurements (single-cell or spatial multi-omics), it presents a flexible approach to interrogate gene function across cell types in vivo, translating gene variants to their causal function.

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

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

Substance specific EEG patterns in mice undergoing slow anesthesia induction

The exact mechanisms and the neural circuits involved in anesthesia induced unconsciousness are still not fully understood. To elucidate them valid animal models are necessary. Since the most commonly used species in neuroscience are mice, we established a murine model for commonly used anesthetics/sedatives and evaluated the epidural electroencephalographic (EEG) patterns during slow anesthesia induction and emergence. Forty-four mice underwent surgery in which we inserted a central venous catheter and implanted nine intracranial electrodes above the prefrontal, motor, sensory, and visual cortex. After at least one week of recovery, mice were anesthetized either by inhalational sevoflurane or intravenous propofol, ketamine, or dexmedetomidine. We evaluated the loss and return of righting reflex (LORR/RORR) and recorded the electrocorticogram. For spectral analysis we focused on the prefrontal and visual cortex. In addition to analyzing the power spectral density at specific time points we evaluated the changes in the spectral power distribution longitudinally. The median time to LORR after start anesthesia ranged from 1080 [1st quartile: 960; 3rd quartile: 1080]s under sevoflurane anesthesia to 1541 [1455; 1890]s with ketamine. Around LORR sevoflurane as well as propofol induced a decrease in the theta/alpha band and an increase in the beta/gamma band. Dexmedetomidine infusion resulted in a shift towards lower frequencies with an increase in the delta range. Ketamine induced stronger activity in the higher frequencies. Our results showed substance-specific changes in EEG patterns during slow anesthesia induction. These patterns were partially identical to previous observations in humans, but also included significant differences, especially in the low frequencies. Our study emphasizes strengths and limitations of murine models in neuroscience and provides an important basis for future studies investigating complex neurophysiological mechanisms.

Characterizing cell type specific transcriptional differences between the living and postmortem human brain

Single-nucleus RNA sequencing (snRNA-seq) is often used to define gene expression patterns characteristic of brain cell types as well as to identify cell type specific gene expression signatures of neurological and mental illnesses in postmortem human brains. As methods to obtain brain tissue from living individuals emerge, it is essential to characterize gene expression differences associated with tissue originating from either living or postmortem subjects using snRNA-seq, and to assess whether and how such differences may impact snRNA-seq studies of brain tissue. To address this, human prefrontal cortex single nuclei gene expression was generated and compared between 31 samples from living individuals and 21 postmortem samples. The same cell types were consistently identified in living and postmortem nuclei, though for each cell type, a large proportion of genes were differentially expressed between samples from postmortem and living individuals. Notably, estimation of cell type proportions by cell type deconvolution of pseudo-bulk data was found to be more accurate in samples from living individuals. To allow for future integration of living and postmortem brain gene expression, a model was developed that quantifies from gene expression data the probability a human brain tissue sample was obtained postmortem. These probabilities are established as a means to statistically account for the gene expression differences between samples from living and postmortem individuals. Together, the results presented here provide a deep characterization of both differences between snRNA-seq derived from samples from living and postmortem individuals, as well as qualify and account for their effect on common analyses performed on this type of data.

Nfe2l3 promotes neuroprotection and long-distance axon regeneration after injury in vivo

Nuclear factor erythroid 2 like (Nfe2l) gene family members 1-3 mediate cellular response to oxidative stress, including in the central nervous system (CNS). However, neuronal functions of Nfe2l3 are unknown. Here, we comparatively evaluated expression of Nfe2l1, Nfe2l2, and Nfe2l3 in singe cell RNA-seq (scRNA-seq)-profiled cortical and retinal ganglion cell (RGC) CNS projection neurons, investigated whether Nfe2l3 regulates neuroprotection and axon regeneration after CNS injury in vivo, and characterized a gene network associated with Nfe2l3 in neurons. We showed that, Nfe2l3 expression transiently peaks in developing immature cortical and RGC projection neurons, but is nearly abolished in adult neurons and is not upregulated after injury. Furthermore, within the retina, Nfe2l3 is enriched in RGCs, primarily neonatally, and not upregulated in injured RGCs, whereas Nfe2l1 and Nfe2l2 are expressed robustly in other retinal cell types as well and are upregulated in injured RGCs. We also found that, expressing Nfe2l3 in injured RGCs through localized intralocular viral vector delivery promotes neuroprotection and long-distance axon regeneration after optic nerve injury in vivo. Moreover, Nfe2l3 provided a similar extent of neuroprotection and axon regeneration as viral vector-targeting of Pten and Klf9, which are prominent regulators of neuroprotection and long-distance axon regeneration. Finally, we bioinformatically characterized a gene network associated with Nfe2l3 in neurons, which predicted the association of Nfe2l3 with established mechanisms of neuroprotection and axon regeneration. Thus, Nfe2l3 is a novel neuroprotection and axon regeneration-promoting factor with a therapeutic potential for treating CNS injury and disease.

Cell-type-resolved mosaicism reveals clonal dynamics of the human forebrain

Debate remains around the anatomical origins of specific brain cell subtypes and lineage relationships within the human forebrain1-7. Thus, direct observation in the mature human brain is critical for a complete understanding of its structural organization and cellular origins. Here we utilize brain mosaic variation within specific cell types as distinct indicators for clonal dynamics, denoted as cell-type-specific mosaic variant barcode analysis. From four hemispheres and two different human neurotypical donors, we identified 287 and 780 mosaic variants, respectively, that were used to deconvolve clonal dynamics. Clonal spread and allele fractions within the brain reveal that local hippocampal excitatory neurons are more lineage-restricted than resident neocortical excitatory neurons or resident basal ganglia GABAergic inhibitory neurons. Furthermore, simultaneous genome transcriptome analysis at both a cell-type-specific and a single-cell level suggests a dorsal neocortical origin for a subgroup of DLX1+ inhibitory neurons that disperse radially from an origin shared with excitatory neurons. Finally, the distribution of mosaic variants across 17 locations within one parietal lobe reveals that restriction of clonal spread in the anterior-posterior axis precedes restriction in the dorsal-ventral axis for both excitatory and inhibitory neurons. Thus, cell-type-resolved somatic mosaicism can uncover lineage relationships governing the development of the human forebrain.

Multimodal cortical neuronal cell type classification

Since more than a century, neuroscientists have distinguished excitatory (glutamatergic) neurons with long-distance projections from inhibitory (GABAergic) neurons with local projections and established layer-dependent schemes for the ~ 80% excitatory (principal) cells as well as the ~ 20% inhibitory neurons. Whereas, in the early days, mainly morphological criteria were used to define cell types, later supplemented by electrophysiological and neurochemical properties, nowadays. single-cell transcriptomics is the method of choice for cell type classification. Bringing recent insight together, we conclude that despite all established layer- and area-dependent differences, there is a set of reliably identifiable cortical cell types that were named (among others) intratelencephalic (IT), extratelencephalic (ET), and corticothalamic (CT) for the excitatory cells, which altogether comprise ~ 56 transcriptomic cell types (t-types). By the same means, inhibitory neurons were subdivided into parvalbumin (PV), somatostatin (SST), vasoactive intestinal polypeptide (VIP), and "other (i.e. Lamp5/Sncg)" subpopulations, which altogether comprise ~ 60 t-types. The coming years will show which t-types actually translate into "real" cell types that show a common set of multimodal features, including not only transcriptome but also physiology and morphology as well as connectivity and ultimately function. Only with the better knowledge of clear-cut cell types and experimental access to them, we will be able to reveal their specific functions, a task which turned out to be difficult in a part of the brain being so much specialized for cognition as the cerebral cortex.

TrkB transmembrane domain: bridging structural understanding with therapeutic strategy

TrkB (neuronal receptor tyrosine kinase-2, NTRK2) is the receptor for brain-derived neurotrophic factor (BDNF) and is a critical regulator of activity-dependent neuronal plasticity. The past few years have witnessed an increasing understanding of the structure and function of TrkB, including its transmembrane domain (TMD). TrkB interacts with membrane cholesterol, which bidirectionally regulates TrkB signaling. Additionally, TrkB has recently been recognized as a binding target of antidepressant drugs. A variety of different antidepressants, including typical and rapid-acting antidepressants, as well as psychedelic compounds, act as allosteric potentiators of BDNF signaling through TrkB. This suggests that TrkB is the common target of different antidepressant compounds. Although more research is needed, current knowledge suggests that TrkB is a promising target for further drug development.

Brain endothelial GSDMD activation mediates inflammatory BBB breakdown

The blood-brain barrier (BBB) protects the central nervous system from infections or harmful substances1; its impairment can lead to or exacerbate various diseases of the central nervous system2-4. However, the mechanisms of BBB disruption during infection and inflammatory conditions5,6 remain poorly defined. Here we find that activation of the pore-forming protein GSDMD by the cytosolic lipopolysaccharide (LPS) sensor caspase-11 (refs. 7-9), but not by TLR4-induced cytokines, mediates BBB breakdown in response to circulating LPS or during LPS-induced sepsis. Mice deficient in the LBP-CD14 LPS transfer and internalization pathway10-12 resist BBB disruption. Single-cell RNA-sequencing analysis reveals that brain endothelial cells (bECs), which express high levels of GSDMD, have a prominent response to circulating LPS. LPS acting on bECs primes Casp11 and Cd14 expression and induces GSDMD-mediated plasma membrane permeabilization and pyroptosis in vitro and in mice. Electron microscopy shows that this features ultrastructural changes in the disrupted BBB, including pyroptotic endothelia, abnormal appearance of tight junctions and vasculature detachment from the basement membrane. Comprehensive mouse genetic analyses, combined with a bEC-targeting adeno-associated virus system, establish that GSDMD activation in bECs underlies BBB disruption by LPS. Delivery of active GSDMD into bECs bypasses LPS stimulation and opens the BBB. In CASP4-humanized mice, Gram-negative Klebsiella pneumoniae infection disrupts the BBB; this is blocked by expression of a GSDMD-neutralizing nanobody in bECs. Our findings outline a mechanism for inflammatory BBB breakdown, and suggest potential therapies for diseases of the central nervous system associated with BBB impairment.

The breath shape controls intonation of mouse vocalizations

Intonation in speech is the control of vocal pitch to layer expressive meaning to communication, like increasing pitch to indicate a question. Also, stereotyped patterns of pitch are used to create distinct sounds with different denotations, like in tonal languages and, perhaps, the ten sounds in the murine lexicon. A basic tone is created by exhalation through a constricted laryngeal voice box, and it is thought that more complex utterances are produced solely by dynamic changes in laryngeal tension. But perhaps, the shifting pitch also results from altering the swiftness of exhalation. Consistent with the latter model, we describe that intonation in most vocalization types follows deviations in exhalation that appear to be generated by the re-activation of the cardinal breathing muscle for inspiration. We also show that the brainstem vocalization central pattern generator, the iRO, can create this breath pattern. Consequently, ectopic activation of the iRO not only induces phonation, but also the pitch patterns that compose most of the vocalizations in the murine lexicon. These results reveal a novel brainstem mechanism for intonation.

Projection-TAGs enable multiplex projection tracing and multi-modal profiling of projection neurons

Single-cell multiomic techniques have sparked immense interest in developing a comprehensive multi-modal map of diverse neuronal cell types and their brain wide projections. However, investigating the spatial organization, transcriptional and epigenetic landscapes of brain wide projection neurons is hampered by the lack of efficient and easily adoptable tools. Here we introduce Projection-TAGs, a retrograde AAV platform that allows multiplex tagging of projection neurons using RNA barcodes. By using Projection-TAGs, we performed multiplex projection tracing of the mouse cortex and high-throughput single-cell profiling of the transcriptional and epigenetic landscapes of the cortical projection neurons. Projection-TAGs can be leveraged to obtain a snapshot of activity-dependent recruitment of distinct projection neurons and their molecular features in the context of a specific stimulus. Given its flexibility, usability, and compatibility, we envision that Projection-TAGs can be readily applied to build a comprehensive multi-modal map of brain neuronal cell types and their projections.

Whole-cortex in situ sequencing reveals input-dependent area identity

The cerebral cortex is composed of neuronal types with diverse gene expression that are organized into specialized cortical areas. These areas, each with characteristic cytoarchitecture1,2, connectivity3,4 and neuronal activity5,6, are wired into modular networks3,4,7. However, it remains unclear whether these spatial organizations are reflected in neuronal transcriptomic signatures and how such signatures are established in development. Here we used BARseq, a high-throughput in situ sequencing technique, to interrogate the expression of 104 cell-type marker genes in 10.3 million cells, including 4,194,658 cortical neurons over nine mouse forebrain hemispheres, at cellular resolution. De novo clustering of gene expression in single neurons revealed transcriptomic types consistent with previous single-cell RNA sequencing studies8,9. The composition of transcriptomic types is highly predictive of cortical area identity. Moreover, areas with similar compositions of transcriptomic types, which we defined as cortical modules, overlap with areas that are highly connected, suggesting that the same modular organization is reflected in both transcriptomic signatures and connectivity. To explore how the transcriptomic profiles of cortical neurons depend on development, we assessed cell-type distributions after neonatal binocular enucleation. Notably, binocular enucleation caused the shifting of the cell-type compositional profiles of visual areas towards neighbouring cortical areas within the same module, suggesting that peripheral inputs sharpen the distinct transcriptomic identities of areas within cortical modules. Enabled by the high throughput, low cost and reproducibility of BARseq, our study provides a proof of principle for the use of large-scale in situ sequencing to both reveal brain-wide molecular architecture and understand its development.

Functional specialization of hippocampal somatostatin-expressing interneurons

Hippocampal somatostatin-expressing (Sst) GABAergic interneurons (INs) exhibit considerable anatomical and functional heterogeneity. Recent single-cell transcriptome analyses have provided a comprehensive Sst-IN subpopulations census, a plausible molecular ground truth of neuronal identity whose links to specific functionality remain incomplete. Here, we designed an approach to identify and access subpopulations of Sst-INs based on transcriptomic features. Four mouse models based on single or combinatorial Cre- and Flp- expression differentiated functionally distinct subpopulations of CA1 hippocampal Sst-INs that largely tiled the morpho-functional parameter space of the Sst-INs superfamily. Notably, the Sst;;Tac1 intersection revealed a population of bistratified INs that preferentially synapsed onto fast-spiking interneurons (FS-INs) and were sufficient to interrupt their firing. In contrast, the Ndnf;;Nkx2-1 intersection identified a population of oriens lacunosum-moleculare INs that predominantly targeted CA1 pyramidal neurons, avoiding FS-INs. Overall, our results provide a framework to translate neuronal transcriptomic identity into discrete functional subtypes that capture the diverse specializations of hippocampal Sst-INs.

Directed differentiation of functional corticospinal-like neurons from endogenous SOX6+/NG2+ cortical progenitors

Corticospinal neurons (CSN) centrally degenerate in amyotrophic lateral sclerosis (ALS), along with spinal motor neurons, and loss of voluntary motor function in spinal cord injury (SCI) results from damage to CSN axons. For functional regeneration of specifically affected neuronal circuitry in vivo , or for optimally informative disease modeling and/or therapeutic screening in vitro , it is important to reproduce the type or subtype of neurons involved. No such appropriate in vitro models exist with which to investigate CSN selective vulnerability and degeneration in ALS, or to investigate routes to regeneration of CSN circuitry for ALS or SCI, critically limiting the relevance of much research. Here, we identify that the HMG-domain transcription factor Sox6 is expressed by a subset of NG2+ endogenous cortical progenitors in postnatal and adult cortex, and that Sox6 suppresses a latent neurogenic program by repressing inappropriate proneural Neurog2 expression by progenitors. We FACS-purify these genetically accessible progenitors from postnatal mouse cortex and establish a pure culture system to investigate their potential for directed differentiation into CSN. We then employ a multi-component construct with complementary and differentiation-sharpening transcriptional controls (activating Neurog2, Fezf2 , while antagonizing Olig2 with VP16:Olig2 ). We generate corticospinal-like neurons from SOX6+/NG2+ cortical progenitors, and find that these neurons differentiate with remarkable fidelity compared with corticospinal neurons in vivo . They possess appropriate morphological, molecular, transcriptomic, and electrophysiological characteristics, without characteristics of the alternate intracortical or other neuronal subtypes. We identify that these critical specifics of differentiation are not reproduced by commonly employed Neurog2 -driven differentiation. Neurons induced by Neurog2 instead exhibit aberrant multi-axon morphology and express molecular hallmarks of alternate cortical projection subtypes, often in mixed form. Together, this developmentally-based directed differentiation from genetically accessible cortical progenitors sets a precedent and foundation for in vitro mechanistic and therapeutic disease modeling, and toward regenerative neuronal repopulation and circuit repair.

Association cortical areas in the mouse contain a large population of fast-spiking GABAergic neurons that do not express parvalbumin

GABAergic neurons represent 10-15% of the neuronal population of the cortex but exert a powerful control over information flow in cortical circuits. The largest GABAergic class in the neocortex is represented by the parvalbumin-expressing fast-spiking neurons, which provide powerful somatic inhibition to their postsynaptic targets. Recently, the density of parvalbumin interneurons has been shown to be lower in associative areas of the mouse cortex as compared with sensory and motor areas. Modelling work based on these quantifications linked the low-density of parvalbumin interneurons with specific computations of associative cortices. However, it is still unknown whether the total GABAergic population of association cortices is smaller or whether another GABAergic type can compensate for the low density of parvalbumin interneurons. In the present study, we investigated these hypotheses using a combination of neuroanatomy, mouse genetics and neurophysiology. We found that the GABAergic population of association areas is comparable with that of primary sensory areas, and it is enriched of fast-spiking neurons that do not express parvalbumin and were not accounted for by previous quantifications. We developed an intersectional viral strategy to demonstrate that the population of fast-spiking neurons is comparable across cortical regions. Our results provide quantifications of the density of fast-spiking GABAergic neurons and offers new biological constrains to refine current models of cortical computations.

Functional diversities within neurons and astrocytes in the adult rat auditory cortex revealed by single-nucleus RNA sequencing

The mammalian cerebral cortex is composed of a rich diversity of cell types. Cortical cells are organized into networks that rely on their functional diversity to ultimately carry out a variety of sophisticated cognitive functions. To investigate the breadth of transcriptional diverse cell types in the sensory cortex, we have used single-nucleus RNA sequencing (snRNA-seq) in the auditory cortex of the adult rat. A variety of unique excitatory and inhibitory neuron types were identified. In addition, we report for the first time a diversity of astrocytes in the auditory cortex that may represent functionally unique subtypes. Together, these results pave the way for building models of how neurons in the sensory cortex work in concert with astrocytes at synapses to fulfill high-cognitive functions like learning and memory.

Multidimensional analysis of cortical interneuron synaptic features reveals underlying synaptic heterogeneity

Cortical interneurons represent a diverse set of neuronal subtypes characterized in part by their striking degree of synaptic specificity. However, little is known about the extent of synaptic diversity because of the lack of unbiased methods to extract synaptic features among interneuron subtypes. Here, we develop an approach to aggregate image features from fluorescent confocal images of interneuron synapses and their post-synaptic targets, in order to characterize the heterogeneity of synapses at fine scale. We started by training a model that recognizes pre- and post-synaptic compartments and then determines the target of each genetically-identified interneuron synapse in vitro and in vivo. Our model extracts hundreds of spatial and intensity features from each analyzed synapse, constructing a multidimensional data set, consisting of millions of synapses, which allowed us to perform an unsupervised analysis on this dataset, uncovering novel synaptic subgroups. The subgroups were spatially distributed in a highly structured manner that revealed the local underlying topology of the postsynaptic environment. Dendrite-targeting subgroups were clustered onto subdomains of the dendrite along the proximal to distal axis. Soma-targeting subgroups were enriched onto different postsynaptic cell types. We also find that the two main subclasses of interneurons, basket cells and somatostatin interneurons, utilize distinct strategies to enact inhibitory coverage. Thus, our analysis of multidimensional synaptic features establishes a conceptual framework for studying interneuron synaptic diversity.

Disruption of Cholinergic Retinal Waves Alters Visual Cortex Development and Function

Retinal waves represent an early form of patterned spontaneous neural activity in the visual system. These waves originate in the retina before eye-opening and propagate throughout the visual system, influencing the assembly and maturation of subcortical visual brain regions. However, because it is technically challenging to ablate retina-derived cortical waves without inducing compensatory activity, the role these waves play in the development of the visual cortex remains unclear. To address this question, we used targeted conditional genetics to disrupt cholinergic retinal waves and their propagation to select regions of primary visual cortex, which largely prevented compensatory patterned activity. We find that loss of cholinergic retinal waves without compensation impaired the molecular and synaptic maturation of excitatory neurons located in the input layers of visual cortex, as well as layer 1 interneurons. These perinatal molecular and synaptic deficits also relate to functional changes observed at later ages. We find that the loss of perinatal cholinergic retinal waves causes abnormal visual cortex retinotopy, mirroring changes in the retinotopic organization of gene expression, and additionally impairs the processing of visual information. We further show that retinal waves are necessary for higher order processing of sensory information by impacting the state-dependent activity of layer 1 interneurons, a neuronal type that shapes neocortical state-modulation, as well as for state-dependent gain modulation of visual responses of excitatory neurons. Together, these results demonstrate that a brief targeted perinatal disruption of patterned spontaneous activity alters early cortical gene expression as well as synaptic and physiological development, and compromises both fundamental and, notably, higher-order functions of visual cortex after eye-opening.

Expression characteristics and potential function of non-coding RNA in mouse cortical cells

Non-coding RNAs (ncRNAs) play essential regulatory functions in various physiological and pathological processes in the brain. To systematically characterize the ncRNA profile in cortical cells, we downloaded single-cell SMART-Seq v4 data of mouse cerebral cortex. Our results revealed that the ncRNAs alone are sufficient to define the identity of most cortical cell types. We identified 1,600 ncRNAs that exhibited cell type specificity, even yielding to distinguish microglia from perivascular macrophages with ncRNA. Moreover, we characterized cortical layer and region specific ncRNAs, in line with the results by spatial transcriptome (ST) data. By constructing a co-expression network of ncRNAs and protein-coding genes, we predicted the function of ncRNAs. By integrating with genome-wide association studies data, we established associations between cell type-specific ncRNAs and traits related to neurological disorders. Collectively, our study identified differentially expressed ncRNAs at multiple levels and provided the valuable resource to explore the functions and dysfunctions of ncRNAs in cortical cells.

Transcriptomic atlas of midbrain dopamine neurons uncovers differential vulnerability in a Parkinsonism lesion model

Midbrain dopamine (mDA) neurons comprise diverse cells with unique innervation targets and functions. This is illustrated by the selective sensitivity of mDA neurons of the substantia nigra compacta (SNc) in patients with Parkinson's disease, while those in the ventral tegmental area (VTA) are relatively spared. Here, we used single nuclei RNA sequencing (snRNA-seq) of approximately 70,000 mouse midbrain cells to build a high-resolution atlas of mouse mDA neuron diversity at the molecular level. The results showed that differences between mDA neuron groups could best be understood as a continuum without sharp differences between subtypes. Thus, we assigned mDA neurons to several 'territories' and 'neighborhoods' within a shifting gene expression landscape where boundaries are gradual rather than discrete. Based on the enriched gene expression patterns of these territories and neighborhoods, we were able to localize them in the adult mouse midbrain. Moreover, because the underlying mechanisms for the variable sensitivities of diverse mDA neurons to pathological insults are not well understood, we analyzed surviving neurons after partial 6-hydroxydopamine (6-OHDA) lesions to unravel gene expression patterns that correlate with mDA neuron vulnerability and resilience. Together, this atlas provides a basis for further studies on the neurophysiological role of mDA neurons in health and disease.

Prefrontal cortical dynorphin peptidergic transmission constrains threat-driven behavioral and network states

Prefrontal cortical (PFC) circuits provide top-down control of threat reactivity. This includes ventromedial PFC (vmPFC) circuitry, which plays a role in suppressing fear-related behavioral states. Dynorphin (Dyn) has been implicated in mediating negative affect and maladaptive behaviors induced by severe threats and is expressed in limbic circuits, including the vmPFC. However, there is a critical knowledge gap in our understanding of how vmPFC Dyn-expressing neurons and Dyn transmission detect threats and regulate expression of defensive behaviors. Here, we demonstrate that Dyn cells are broadly activated by threats and release Dyn locally in the vmPFC to limit passive defensive behaviors. We further demonstrate that vmPFC Dyn-mediated signaling promotes a switch of vmPFC networks to a fear-related state. In conclusion, we reveal a previously unknown role of vmPFC Dyn neurons and Dyn neuropeptidergic transmission in suppressing defensive behaviors in response to threats via state-driven changes in vmPFC networks.

Detailed phenotyping of Tbr1-2A-CreER knock-in mice demonstrates significant impacts on TBR1 protein levels and axon development

Spatiotemporal control of Cre-mediated recombination has been an invaluable tool for understanding key developmental processes. For example, knock-in of Cre into cell type marker gene loci drives Cre expression under endogenous promoter and enhancer sequences, greatly facilitating the study of diverse neuronal subtypes in the cerebral cortex. However, insertion of exogenous DNA into the genome can have unintended effects on local gene regulation or protein function that must be carefully considered. Here, we analyze a recently generated Tbr1-2A-CreER knock-in mouse line, where a 2A-CreER cassette was inserted in-frame just before the stop codon of the transcription factor gene Tbr1 . Heterozygous TBR1 mutations in humans and mice are known to cause autism or autism-like behavioral phenotypes accompanied by structural brain malformations, most frequently a reduction of the anterior commissure. Thus, it is critical for modified versions of Tbr1 to exhibit true wild-type-like activity. We evaluated the Tbr1-2A-CreER allele for its potential impact on Tbr1 function and complementation to Tbr1 loss-of-function alleles. In mice with one copy of the Tbr1-2A-CreER allele, we identified reduction of TBR1 protein in early postnatal cortex along with thinning of the anterior commissure, suggesting hypersensitivity of this structure to TBR1 dosage. Comparing Tbr1-2A-CreER and Tbr1 -null heterozygous and homozygous mice to Tbr1 -null complementation crosses showed reductions of TBR1 dosage ranging from 28.4% to 95.9%. Using these combinatorial genotypes, we found that low levels of TBR1 protein (∼16%) are sufficient to establish cortical layer positioning, while greater levels (>50%) are required for normal suppression of layer 5 identity. In total, these results strongly support the conclusion that Tbr1-2A-CreER is a hypomorphic allele. We advise caution when interpreting experiments using this allele, such as transcriptomic studies, considering the sensitivity of various corticogenic processes to TBR1 dosage and the association of heterozygous TBR1 mutations with complex neurodevelopmental disorders.

Targeted blockade of aberrant sodium current in a stem cell-derived neuron model of SCN3A encephalopathy

Missense variants in SCN3A encoding the voltage-gated sodium (Na+) channel α subunit Nav1.3 are associated with SCN3A-related neurodevelopmental disorder (SCN3A-NDD), a spectrum of disease that includes epilepsy and malformation of cortical development. How genetic variation in SCN3A leads to pathology remains unclear, as prior electrophysiological work on disease-associated variants has been performed exclusively in heterologous cell systems. To further investigate the mechanisms of SCN3A-NDD pathogenesis, we used CRISPR/Cas9 gene editing to modify a control human induced pluripotent stem cell (iPSC) line to express the recurrent de novo missense variant SCN3A c.2624T>C (p.Ile875Thr). With the established Ngn2 rapid induction protocol, we generated glutamatergic forebrain-like neurons (iNeurons), which we showed to express SCN3A mRNA and Nav1.3-mediated Na+ currents. We performed detailed whole-cell patch clamp recordings to determine the effect of the SCN3A-p.Ile875Thr variant on endogenous Na+ currents in, and intrinsic excitability of, human neurons. Compared to control iNeurons, variant-expressing iNeurons exhibit markedly increased slowly-inactivating/persistent Na+ current, abnormal firing patterns with paroxysmal bursting and plateau-like potentials with action potential failure, and a hyperpolarized voltage threshold for action potential generation. We then validated these findings using a separate iPSC line generated from a patient harbouring the SCN3A-p.Ile875Thr variant compared to a corresponding CRISPR-corrected isogenic control line. Finally, we found that application of the Nav1.3-selective blocker ICA-121431 normalizes action potential threshold and aberrant firing patterns in SCN3A-p.Ile1875Thr iNeurons; in contrast, consistent with action as a Na+ channel blocker, ICA-121431 decreases excitability of control iNeurons. Our findings demonstrate that iNeurons can model the effects of genetic variation in SCN3A yet reveal a complex relationship between gain-of-function at the level of the ion channel versus impact on neuronal excitability. Given the transient expression of SCN3A in the developing human nervous system, selective blockade or suppression of Nav1.3-containing Na+ channels could represent a therapeutic approach towards SCN3A-NDD.

Pianno: a probabilistic framework automating semantic annotation for spatial transcriptomics

Spatial transcriptomics has revolutionized the study of gene expression within tissues, while preserving spatial context. However, annotating spatial spots' biological identity remains a challenge. To tackle this, we introduce Pianno, a Bayesian framework automating structural semantics annotation based on marker genes. Comprehensive evaluations underscore Pianno's remarkable prowess in precisely annotating a wide array of spatial semantics, ranging from diverse anatomical structures to intricate tumor microenvironments, as well as in estimating cell type distributions, across data generated from various spatial transcriptomics platforms. Furthermore, Pianno, in conjunction with clustering approaches, uncovers a region- and species-specific excitatory neuron subtype in the deep layer 3 of the human neocortex, shedding light on cellular evolution in the human neocortex. Overall, Pianno equips researchers with a robust and efficient tool for annotating diverse biological structures, offering new perspectives on spatial transcriptomics data.

Native-state proteomics of Parvalbumin interneurons identifies unique molecular signatures and vulnerabilities to early Alzheimer's pathology

Dysfunction in fast-spiking parvalbumin interneurons (PV-INs) may represent an early pathophysiological perturbation in Alzheimer's Disease (AD). Defining early proteomic alterations in PV-INs can provide key biological and translationally-relevant insights. We used cell-type-specific in-vivo biotinylation of proteins (CIBOP) coupled with mass spectrometry to obtain native-state PV-IN proteomes. PV-IN proteomic signatures include high metabolic and translational activity, with over-representation of AD-risk and cognitive resilience-related proteins. In bulk proteomes, PV-IN proteins were associated with cognitive decline in humans, and with progressive neuropathology in humans and the 5xFAD mouse model of Aβ pathology. PV-IN CIBOP in early stages of Aβ pathology revealed signatures of increased mitochondria and metabolism, synaptic and cytoskeletal disruption and decreased mTOR signaling, not apparent in whole-brain proteomes. Furthermore, we demonstrated pre-synaptic defects in PV-to-excitatory neurotransmission, validating our proteomic findings. Overall, in this study we present native-state proteomes of PV-INs, revealing molecular insights into their unique roles in cognitive resiliency and AD pathogenesis.

Multimodal Nature of the Single-cell Primate Brain Atlas: Morphology, Transcriptome, Electrophysiology, and Connectivity

Primates exhibit complex brain structures that augment cognitive function. The neocortex fulfills high-cognitive functions through billions of connected neurons. These neurons have distinct transcriptomic, morphological, and electrophysiological properties, and their connectivity principles vary. These features endow the primate brain atlas with a multimodal nature. The recent integration of next-generation sequencing with modified patch-clamp techniques is revolutionizing the way to census the primate neocortex, enabling a multimodal neuronal atlas to be established in great detail: (1) single-cell/single-nucleus RNA-seq technology establishes high-throughput transcriptomic references, covering all major transcriptomic cell types; (2) patch-seq links the morphological and electrophysiological features to the transcriptomic reference; (3) multicell patch-clamp delineates the principles of local connectivity. Here, we review the applications of these technologies in the primate neocortex and discuss the current advances and tentative gaps for a comprehensive understanding of the primate neocortex.

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

Transcranial ultrasound activates mechanosensitive cellular signaling and modulates neural dynamics. Given that intrinsic neuronal activity is limited to a couple hundred hertz and often exhibits frequency preference, we examined whether pulsing ultrasound at physiologic pulse repetition frequencies (PRFs) could selectively influence neuronal activity in the mammalian brain. We performed calcium imaging of individual motor cortex neurons, while delivering 0.35 MHz ultrasound at PRFs of 10, 40, and 140 Hz in awake mice. We found that most neurons were preferentially activated by only one of the three PRFs, highlighting unique cellular effects of physiologic PRFs. Further, ultrasound evoked responses were similar between excitatory neurons and parvalbumin positive interneurons regardless of PRFs, indicating that individual cell sensitivity dominates ultrasound-evoked effects, consistent with the heterogeneous mechanosensitive channel expression we found across single neurons in mice and humans. These results highlight the feasibility of tuning ultrasound neuromodulation effects through varying PRFs.

stDiff: a diffusion model for imputing spatial transcriptomics through single-cell transcriptomics

Spatial transcriptomics (ST) has become a powerful tool for exploring the spatial organization of gene expression in tissues. Imaging-based methods, though offering superior spatial resolutions at the single-cell level, are limited in either the number of imaged genes or the sensitivity of gene detection. Existing approaches for enhancing ST rely on the similarity between ST cells and reference single-cell RNA sequencing (scRNA-seq) cells. In contrast, we introduce stDiff, which leverages relationships between gene expression abundance in scRNA-seq data to enhance ST. stDiff employs a conditional diffusion model, capturing gene expression abundance relationships in scRNA-seq data through two Markov processes: one introducing noise to transcriptomics data and the other denoising to recover them. The missing portion of ST is predicted by incorporating the original ST data into the denoising process. In our comprehensive performance evaluation across 16 datasets, utilizing multiple clustering and similarity metrics, stDiff stands out for its exceptional ability to preserve topological structures among cells, positioning itself as a robust solution for cell population identification. Moreover, stDiff's enhancement outcomes closely mirror the actual ST data within the batch space. Across diverse spatial expression patterns, our model accurately reconstructs them, delineating distinct spatial boundaries. This highlights stDiff's capability to unify the observed and predicted segments of ST data for subsequent analysis. We anticipate that stDiff, with its innovative approach, will contribute to advancing ST imputation methodologies.

Parallel pathways carrying direction-and orientation-selective retinal signals to layer 4 of the mouse visual cortex

Parallel visual pathways from the retina to the primary visual cortex (V1) via the lateral geniculate nucleus are common to many mammalian species, including mice, carnivores, and primates. However, it remains unclear which visual features present in both retina and V1 may be inherited from parallel pathways versus extracted by V1 circuits in the mouse. Here, using calcium imaging and rabies circuit tracing, we explore the relationships between tuning of layer 4 (L4) V1 neurons and their retinal ganglion cell (RGC) inputs. We find that subpopulations of L4 V1 neurons differ in their tuning for direction, orientation, spatial frequency, temporal frequency, and speed. Furthermore, we find that direction-tuned L4 V1 neurons receive input from direction-selective RGCs, whereas orientation-tuned L4 V1 neurons receive input from orientation-selective RGCs. These results suggest that direction and orientation tuning of V1 neurons may be partly inherited from parallel pathways originating in the retina.

Specific and comprehensive genetic targeting reveals brain-wide distribution and synaptic input patterns of GABAergic axo-axonic interneurons

Axo-axonic cells (AACs), also called chandelier cells (ChCs) in the cerebral cortex, are the most distinctive type of GABAergic interneurons described in the neocortex, hippocampus, and basolateral amygdala (BLA). AACs selectively innervate glutamatergic projection neurons (PNs) at their axon initial segment (AIS), thus may exert decisive control over PN spiking and regulate PN functional ensembles. However, the brain-wide distribution, synaptic connectivity, and circuit function of AACs remains poorly understood, largely due to the lack of specific and reliable experimental tools. Here, we have established an intersectional genetic strategy that achieves specific and comprehensive targeting of AACs throughout the mouse brain based on their lineage (Nkx2.1) and molecular (Unc5b, Pthlh) markers. We discovered that AACs are deployed across essentially all the pallium-derived brain structures, including not only the dorsal pallium-derived neocortex and medial pallium-derived hippocampal formation, but also the lateral pallium-derived claustrum-insular complex, and the ventral pallium-derived extended amygdaloid complex and olfactory centers. AACs are also abundant in anterior olfactory nucleus, taenia tecta and lateral septum. AACs show characteristic variations in density across neocortical areas and layers and across subregions of the hippocampal formation. Neocortical AACs comprise multiple laminar subtypes with distinct dendritic and axonal arborization patterns. Retrograde monosynaptic tracing from AACs across neocortical, hippocampal and BLA regions reveal shared as well as distinct patterns of synaptic input. Specific and comprehensive targeting of AACs facilitates the study of their developmental genetic program and circuit function across brain structures, providing a ground truth platform for understanding the conservation and variation of a bona fide cell type across brain regions and species.

Neuronal subtype-specific transcriptomic changes in the cerebral neocortex associated with sleep pressure

Sleep is homeostatically regulated by sleep pressure, which increases during wakefulness and dissipates during sleep. Recent studies have suggested that the cerebral neocortex, a six-layered structure composed of various layer- and projection-specific neuronal subtypes, is involved in the representation of sleep pressure governed by transcriptional regulation. Here, we examined the transcriptomic changes in neuronal subtypes in the neocortex upon increased sleep pressure using single-nucleus RNA sequencing datasets and predicted the putative intracellular and intercellular molecules involved in transcriptome alterations. We revealed that sleep deprivation (SD) had the greatest effect on the transcriptome of layer 2 and 3 intratelencephalic (L2/3 IT) neurons among the neocortical glutamatergic neuronal subtypes. The expression of mutant SIK3 (SLP), which is known to increase sleep pressure, also induced profound changes in the transcriptome of L2/3 IT neurons. We identified Junb as a candidate transcription factor involved in the alteration of the L2/3 IT neuronal transcriptome by SD and SIK3 (SLP) expression. Finally, we inferred putative intercellular ligands, including BDNF, LSAMP, and PRNP, which may be involved in SD-induced alteration of the transcriptome of L2/3 IT neurons. We suggest that the transcriptome of L2/3 IT neurons is most impacted by increased sleep pressure among neocortical glutamatergic neuronal subtypes and identify putative molecules involved in such transcriptional alterations.

Reactive astrocytes secrete the chaperone HSPB1 to mediate neuroprotection

Molecular chaperones are protective in neurodegenerative diseases by preventing protein misfolding and aggregation, such as extracellular amyloid plaques and intracellular tau neurofibrillary tangles in Alzheimer's disease (AD). In addition, AD is characterized by an increase in astrocyte reactivity. The chaperone HSPB1 has been proposed as a marker for reactive astrocytes; however, its astrocytic functions in neurodegeneration remain to be elucidated. Here, we identify that HSPB1 is secreted from astrocytes to exert non-cell-autonomous protective functions. We show that in human AD brain, HSPB1 levels increase in astrocytes that cluster around amyloid plaques, as well as in the adjacent extracellular space. Moreover, in conditions that mimic an inflammatory reactive response, astrocytes increase HSPB1 secretion. Concomitantly, astrocytes and neurons can uptake astrocyte-secreted HSPB1, which is accompanied by an attenuation of the inflammatory response in reactive astrocytes and reduced pathological tau inclusions. Our findings highlight a protective mechanism in disease conditions that encompasses the secretion of a chaperone typically regarded as intracellular.

Molecular architecture of primate specific neural circuit formation

The mammalian cortex is a highly evolved brain region, but we still lack a comprehensive understanding of the molecular mechanisms underlying primate-specific neural circuits formation. In this study, we employed spatial transcriptomics to assess gene expression dynamics in the marmoset cortex during development, focusing on key regions and time points. Spatial transcriptomics identified genes that are sexually, spatially, and temporally differentially expressed in the developing marmoset cortex. Our detailed analysis of the visual cortex unveiled dynamic changes in gene expression across layers with distinct projections and functions. Notably, we discovered numerous axon guidance molecules with spatiotemporal expression patterns unique to the developing marmoset prefrontal cortex (PFC), which control PFC neuronal circuits. Among these molecules, PRSS12 (Protease, Serine, 12 (neurotrypsin, motopsin), when ectopically expressed in the mouse prelimbic cortex, caused similar changes in connectivity as observed in the marmoset A32 area. Furthermore, PRSS12 showed similar expression patterns in both marmoset and human PFC during development, suggesting parallels between marmoset and human brain development. The differential expression of axon guidance molecules in the developing PFC, varying by region, likely contributes to the formation of unique circuits observed in primates.

Selective vulnerability of layer 5a corticostriatal neurons in Huntington's disease

The properties of the cell types that are selectively vulnerable in Huntington's disease (HD) cortex, the nature of somatic CAG expansions of mHTT in these cells, and their importance in CNS circuitry have not been delineated. Here, we employed serial fluorescence-activated nuclear sorting (sFANS), deep molecular profiling, and single-nucleus RNA sequencing (snRNA-seq) of motor-cortex samples from thirteen predominantly early stage, clinically diagnosed HD donors and selected samples from cingulate, visual, insular, and prefrontal cortices to demonstrate loss of layer 5a pyramidal neurons in HD. Extensive mHTT CAG expansions occur in vulnerable layer 5a pyramidal cells, and in Betz cells, layers 6a and 6b neurons that are resilient in HD. Retrograde tracing experiments in macaque brains identify layer 5a neurons as corticostriatal pyramidal cells. We propose that enhanced somatic mHTT CAG expansion and altered synaptic function act together to cause corticostriatal disconnection and selective neuronal vulnerability in HD cerebral cortex.

Highly sensitive spatial transcriptomics using FISHnCHIPs of multiple co-expressed genes

High-dimensional, spatially resolved analysis of intact tissue samples promises to transform biomedical research and diagnostics, but existing spatial omics technologies are costly and labor-intensive. We present Fluorescence In Situ Hybridization of Cellular HeterogeneIty and gene expression Programs (FISHnCHIPs) for highly sensitive in situ profiling of cell types and gene expression programs. FISHnCHIPs achieves this by simultaneously imaging ~2-35 co-expressed genes (clustered into modules) that are spatially co-localized in tissues, resulting in similar spatial information as single-gene Fluorescence In Situ Hybridization (FISH), but with ~2-20-fold higher sensitivity. Using FISHnCHIPs, we image up to 53 modules from the mouse kidney and mouse brain, and demonstrate high-speed, large field-of-view profiling of a whole tissue section. FISHnCHIPs also reveals spatially restricted localizations of cancer-associated fibroblasts in a human colorectal cancer biopsy. Overall, FISHnCHIPs enables fast, robust, and scalable cell typing of tissues with normal physiology or undergoing pathogenesis.

Protocol to study microcircuits in the medial entorhinal cortex in mice using multiple patch-clamp recordings and morphological reconstruction

Multiple patch-clamp recordings and morphological reconstruction are powerful approaches for neuronal microcircuitry dissection and cell type classification but are challenging due to the sophisticated expertise needed. Here, we present a protocol for applying these techniques to neurons in the medial entorhinal cortex (MEC) of mice. We detail steps to prepare brain slices containing MEC and perform simultaneous multiple whole-cell recordings, followed by procedures of histological staining and neuronal reconstruction. We then describe how we analyze morphological and electrophysiological features. For complete details on the use and execution of this protocol, please refer to Shi et al.1.

Most discriminative stimuli for functional cell type clustering

Identifying cell types and understanding their functional properties is crucial for unraveling the mechanisms underlying perception and cognition. In the retina, functional types can be identified by carefully selected stimuli, but this requires expert domain knowledge and biases the procedure towards previously known cell types. In the visual cortex, it is still unknown what functional types exist and how to identify them. Thus, for unbiased identification of the functional cell types in retina and visual cortex, new approaches are needed. Here we propose an optimization-based clustering approach using deep predictive models to obtain functional clusters of neurons using Most Discriminative Stimuli (MDS). Our approach alternates between stimulus optimization with cluster reassignment akin to an expectation-maximization algorithm. The algorithm recovers functional clusters in mouse retina, marmoset retina and macaque visual area V4. This demonstrates that our approach can successfully find discriminative stimuli across species, stages of the visual system and recording techniques. The resulting most discriminative stimuli can be used to assign functional cell types fast and on the fly, without the need to train complex predictive models or show a large natural scene dataset, paving the way for experiments that were previously limited by experimental time. Crucially, MDS are interpretable: they visualize the distinctive stimulus patterns that most unambiguously identify a specific type of neuron.

Heterogeneity of layer 4 in visual areas of rhesus macaque cortex

Recently, single-cell RNA-sequencing (scRNA-seq) has enabled unprecedented insights to the cellular landscape of the brains of many different species, among them the rhesus macaque as a key animal model. Building on previous, broader surveys of the macaque brain, we closely examined five immediately neighboring areas within the visual cortex of the rhesus macaque: V1, V2, V4, MT and TEO. To facilitate this, we first devised a novel pipeline for brain spatial archive - the BrainSPACE - which enabled robust archiving and sampling from the whole unfixed brain. SnRNA-sequencing of ~100,000 nuclei from visual areas V1 and V4 revealed conservation within the GABAergic neuron subtypes, while seven and one distinct principle neuron subtypes were detected in V1 and V4, respectively, all most likely located in layer 4. Moreover, using small molecule fluorescence in situ hybridization, we identified cell type density gradients across V1, V2, V4, MT, and TEO appearing to reflect the visual hierarchy. These findings demonstrate an association between the clear areal specializations among neighboring areas with the hierarchical levels within the visual cortex of the rhesus macaque.