Thyroid hormone remodels cortex to coordinate body-wide metabolism and exploration.

Animals adapt to environmental conditions by modifying the function of their internal organs, including the brain. To be adaptive, alterations in behavior must be coordinated with the functional state of organs throughout the body. Here, we find that thyroid hormone-a regulator of metabolism in many peripheral organs-directly activates cell-type-specific transcriptional programs in the frontal cortex of adult male mice. These programs are enriched for axon-guidance genes in glutamatergic projection neurons, synaptic regulatory genes in both astrocytes and neurons, and pro-myelination factors in oligodendrocytes, suggesting widespread plasticity of cortical circuits. Indeed, whole-cell electrophysiology revealed that thyroid hormone alters excitatory and inhibitory synaptic transmission, an effect that requires thyroid hormone-induced gene regulatory programs in presynaptic neurons. Furthermore, thyroid hormone action in the frontal cortex regulates innate exploratory behaviors and causally promotes exploratory decision-making. Thus, thyroid hormone acts directly on the cerebral cortex in males to coordinate exploratory behaviors with whole-body metabolic state.

Scalable imaging-free spatial genomics through computational reconstruction.

Tissue organization arises from the coordinated molecular programs of cells. Spatial genomics maps cells and their molecular programs within the spatial context of tissues. However, current methods measure spatial information through imaging or direct registration, which often require specialized equipment and are limited in scale. Here, we developed an imaging-free spatial transcriptomics method that uses molecular diffusion patterns to computationally reconstruct spatial data. To do so, we utilize a simple experimental protocol on two dimensional barcode arrays to establish an interaction network between barcodes via molecular diffusion. Sequencing these interactions generates a high dimensional matrix of interactions between different spatial barcodes. Then, we perform dimensionality reduction to regenerate a two-dimensional manifold, which represents the spatial locations of the barcode arrays. Surprisingly, we found that the UMAP algorithm, with minimal modifications can faithfully successfully reconstruct the arrays. We demonstrated that this method is compatible with capture array based spatial transcriptomics/genomics methods, Slide-seq and Slide-tags, with high fidelity. We systematically explore the fidelity of the reconstruction through comparisons with experimentally derived ground truth data, and demonstrate that reconstruction generates high quality spatial genomics data. We also scaled this technique to reconstruct high-resolution spatial information over areas up to 1.2 centimeters. This computational reconstruction method effectively converts spatial genomics measurements to molecular biology, enabling spatial transcriptomics with high accessibility, and scalability.

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

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

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

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

Specialized connectivity of molecular layer interneuron subtypes leads to disinhibition and synchronous inhibition of cerebellar Purkinje cells.

Molecular layer interneurons (MLIs) account for approximately 80% of the inhibitory interneurons in the cerebellar cortex and are vital to cerebellar processing. MLIs are thought to primarily inhibit Purkinje cells (PCs) and suppress the plasticity of synapses onto PCs. MLIs also inhibit, and are electrically coupled to, other MLIs, but the functional significance of these connections is not known. Here, we find that two recently recognized MLI subtypes, MLI1 and MLI2, have a highly specialized connectivity that allows them to serve distinct functional roles. MLI1s primarily inhibit PCs, are electrically coupled to each other, fire synchronously with other MLI1s on the millisecond timescale in vivo, and synchronously pause PC firing. MLI2s are not electrically coupled, primarily inhibit MLI1s and disinhibit PCs, and are well suited to gating cerebellar-dependent behavior and learning. The synchronous firing of electrically coupled MLI1s and disinhibition provided by MLI2s require a major re-evaluation of cerebellar processing.

Combinatorial expression of neurexin genes regulates glomerular targeting by olfactory sensory neurons.

Precise connectivity between specific neurons is essential for the formation of the complex neural circuitry necessary for executing intricate motor behaviors and higher cognitive functions. While trans -interactions between synaptic membrane proteins have emerged as crucial elements in orchestrating the assembly of these neural circuits, the synaptic surface proteins involved in neuronal wiring remain largely unknown. Here, using unbiased single-cell transcriptomic and mouse genetic approaches, we uncover that the neurexin family of genes enables olfactory sensory neuron (OSNs) axons to form appropriate synaptic connections with their mitral and tufted (M/T) cell synaptic partners, within the mammalian olfactory system. Neurexin isoforms are differentially expressed within distinct populations of OSNs, resulting in unique pattern of neurexin expression that is specific to each OSN type, and synergistically cooperate to regulate axonal innervation, guiding OSN axons to their designated glomeruli. This process is facilitated through the interactions of neurexins with their postsynaptic partners, including neuroligins, which have distinct expression patterns in M/T cells. Our findings suggest a novel mechanism underpinning the precise assembly of olfactory neural circuits, driven by the trans -interaction between neurexins and their ligands.

Slide-tags enables single-nucleus barcoding for multimodal spatial genomics.

Recent technological innovations have enabled the high-throughput quantification of gene expression and epigenetic regulation within individual cells, transforming our understanding of how complex tissues are constructed1-6. However, missing from these measurements is the ability to routinely and easily spatially localize these profiled cells. We developed a strategy, Slide-tags, in which single nuclei within an intact tissue section are tagged with spatial barcode oligonucleotides derived from DNA-barcoded beads with known positions. These tagged nuclei can then be used as an input into a wide variety of single-nucleus profiling assays. Application of Slide-tags to the mouse hippocampus positioned nuclei at less than 10 µm spatial resolution and delivered whole-transcriptome data that are indistinguishable in quality from ordinary single-nucleus RNA-sequencing data. To demonstrate that Slide-tags can be applied to a wide variety of human tissues, we performed the assay on brain, tonsil and melanoma. We revealed cell-type-specific spatially varying gene expression across cortical layers and spatially contextualized receptor-ligand interactions driving B cell maturation in lymphoid tissue. A major benefit of Slide-tags is that it is easily adaptable to almost any single-cell measurement technology. As a proof of principle, we performed multiomic measurements of open chromatin, RNA and T cell receptor (TCR) sequences in the same cells from metastatic melanoma, identifying transcription factor motifs driving cancer cell state transitions in spatially distinct microenvironments. Slide-tags offers a universal platform for importing the compendium of established single-cell measurements into the spatial genomics repertoire.

The molecular cytoarchitecture of the adult mouse brain.

The function of the mammalian brain relies upon the specification and spatial positioning of diversely specialized cell types. Yet, the molecular identities of the cell types and their positions within individual anatomical structures remain incompletely known. To construct a comprehensive atlas of cell types in each brain structure, we paired high-throughput single-nucleus RNA sequencing with Slide-seq1,2-a recently developed spatial transcriptomics method with near-cellular resolution-across the entire mouse brain. Integration of these datasets revealed the cell type composition of each neuroanatomical structure. Cell type diversity was found to be remarkably high in the midbrain, hindbrain and hypothalamus, with most clusters requiring a combination of at least three discrete gene expression markers to uniquely define them. Using these data, we developed a framework for genetically accessing each cell type, comprehensively characterized neuropeptide and neurotransmitter signalling, elucidated region-specific specializations in activity-regulated gene expression and ascertained the heritability enrichment of neurological and psychiatric phenotypes. These data, available as an online resource ( www.BrainCellData.org ), should find diverse applications across neuroscience, including the construction of new genetic tools and the prioritization of specific cell types and circuits in the study of brain diseases.

A Purkinje cell to parabrachial nucleus pathway enables broad cerebellar influence over the forebrain.

In addition to its motor functions, the cerebellum is involved in emotional regulation, anxiety and affect. We found that suppressing the firing of cerebellar Purkinje cells (PCs) rapidly excites forebrain areas that contribute to such functions (including the amygdala, basal forebrain and septum), but that the classic cerebellar outputs, the deep cerebellar nuclei, do not directly project there. We show that PCs directly inhibit parabrachial nuclei (PBN) neurons that project to numerous forebrain regions. Suppressing the PC-PBN pathway influences many regions in the forebrain and is aversive. Molecular profiling shows that PCs directly inhibit numerous types of PBN neurons that control diverse behaviors that are not involved in motor control. Therefore, the PC-PBN pathway allows the cerebellum to directly regulate activity in the forebrain, and may be an important substrate for cerebellar disorders arising from damage to the posterior vermis.

Cerebellar circuits for disinhibition and synchronous inhibition.

The cerebellar cortex contributes to diverse behaviors by transforming mossy fiber inputs into predictions in the form of Purkinje cell (PC) outputs, and then refining those predictions1. Molecular layer interneurons (MLIs) account for approximately 80% of the inhibitory interneurons in the cerebellar cortex2, and are vital to cerebellar processing1,3. MLIs are thought to primarily inhibit PCs and suppress the plasticity of excitatory synapses onto PCs. MLIs also inhibit, and are electrically coupled to, other MLIs4-7, but the functional significance of these connections is not known1,3. Behavioral studies suggest that cerebellar-dependent learning is gated by disinhibition of PCs, but the source of such disinhibition has not been identified8. Here we find that two recently recognized MLI subtypes2, MLI1 and MLI2, have highly specialized connectivity that allows them to serve very different functional roles. MLI1s primarily inhibit PCs, are electrically coupled to each other, fire synchronously with other MLI1s on the millisecond time scale in vivo, and synchronously pause PC firing. MLI2s are not electrically coupled, they primarily inhibit MLI1s and disinhibit PCs, and are well suited to gating cerebellar-dependent learning8. These findings require a major reevaluation of processing within the cerebellum in which disinhibition, a powerful circuit motif present in the cerebral cortex and elsewhere9-17, greatly increases the computational power and flexibility of the cerebellum. They also suggest that millisecond time scale synchronous firing of electrically-coupled MLI1s helps regulate the output of the cerebellar cortex by synchronously pausing PC firing, which has been shown to evoke precisely-timed firing in PC targets18.

Brain-region-specific changes in neurons and glia and dysregulation of dopamine signaling in Grin2a mutant mice.

A genetically valid animal model could transform our understanding of schizophrenia (SCZ) disease mechanisms. Rare heterozygous loss-of-function (LoF) mutations in GRIN2A, encoding a subunit of the NMDA receptor, greatly increase the risk of SCZ. By transcriptomic, proteomic, and behavioral analyses, we report that heterozygous Grin2a mutant mice show (1) large-scale gene expression changes across multiple brain regions and in neuronal (excitatory and inhibitory) and non-neuronal cells (astrocytes and oligodendrocytes), (2) evidence of hypoactivity in the prefrontal cortex (PFC) and hyperactivity in the hippocampus and striatum, (3) an elevated dopamine signaling in the striatum and hypersensitivity to amphetamine-induced hyperlocomotion (AIH), (4) altered cholesterol biosynthesis in astrocytes, (5) a reduction in glutamatergic receptor signaling proteins in the synapse, and (6) an aberrant locomotor pattern opposite of that induced by antipsychotic drugs. These findings reveal potential pathophysiologic mechanisms, provide support for both the "hypo-glutamate" and "hyper-dopamine" hypotheses of SCZ, and underscore the utility of Grin2a-deficient mice as a genetic model of SCZ.

Early Alzheimer's disease pathology in human cortex involves transient cell states.

Cellular perturbations underlying Alzheimer's disease (AD) are primarily studied in human postmortem samples and model organisms. Here, we generated a single-nucleus atlas from a rare cohort of cortical biopsies from living individuals with varying degrees of AD pathology. We next performed a systematic cross-disease and cross-species integrative analysis to identify a set of cell states that are specific to early AD pathology. These changes-which we refer to as the early cortical amyloid response-were prominent in neurons, wherein we identified a transitional hyperactive state preceding the loss of excitatory neurons, which we confirmed by acute slice physiology on independent biopsy specimens. Microglia overexpressing neuroinflammatory-related processes also expanded as AD pathology increased. Finally, both oligodendrocytes and pyramidal neurons upregulated genes associated with ß-amyloid production and processing during this early hyperactive phase. Our integrative analysis provides an organizing framework for targeting circuit dysfunction, neuroinflammation, and amyloid production early in AD pathogenesis.

Schizophrenia-associated somatic copy-number variants from 12,834 cases reveal recurrent NRXN1 and ABCB11 disruptions.

While germline copy-number variants (CNVs) contribute to schizophrenia (SCZ) risk, the contribution of somatic CNVs (sCNVs)-present in some but not all cells-remains unknown. We identified sCNVs using blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls, filtering sCNVs at loci recurrently mutated in clonal blood disorders. Likely early-developmental sCNVs were more common in cases (0.91%) than controls (0.51%, p = 2.68e-4), with recurrent somatic deletions of exons 1-5 of the NRXN1 gene in five SCZ cases. Hi-C maps revealed ectopic, allele-specific loops forming between a potential cryptic promoter and non-coding cis-regulatory elements upon 5' deletions in NRXN1. We also observed recurrent intragenic deletions of ABCB11, encoding a transporter implicated in anti-psychotic response, in five treatment-resistant SCZ cases and showed that ABCB11 is specifically enriched in neurons forming mesocortical and mesolimbic dopaminergic projections. Our results indicate potential roles of sCNVs in SCZ risk.

Thyroid hormone rewires cortical circuits to coordinate body-wide metabolism and exploratory drive.

Animals adapt to varying environmental conditions by modifying the function of their internal organs, including the brain. To be adaptive, alterations in behavior must be coordinated with the functional state of organs throughout the body. Here we find that thyroid hormone- a prominent regulator of metabolism in many peripheral organs- activates cell-type specific transcriptional programs in anterior regions of cortex of adult mice via direct activation of thyroid hormone receptors. These programs are enriched for axon-guidance genes in glutamatergic projection neurons, synaptic regulators across both astrocytes and neurons, and pro-myelination factors in oligodendrocytes, suggesting widespread remodeling of cortical circuits. Indeed, whole-cell electrophysiology recordings revealed that thyroid hormone induces local transcriptional programs that rewire cortical neural circuits via pre-synaptic mechanisms, resulting in increased excitatory drive with a concomitant sensitization of recruited inhibition. We find that thyroid hormone bidirectionally regulates innate exploratory behaviors and that the transcriptionally mediated circuit changes in anterior cortex causally promote exploratory decision-making. Thus, thyroid hormone acts directly on adult cerebral cortex to coordinate exploratory behaviors with whole-body metabolic state.

Spatiotemporal transcriptomic maps of whole mouse embryos at the onset of organogenesis.

Spatiotemporal orchestration of gene expression is required for proper embryonic development. The use of single-cell technologies has begun to provide improved resolution of early regulatory dynamics, including detailed molecular definitions of most cell states during mouse embryogenesis. Here we used Slide-seq to build spatial transcriptomic maps of complete embryonic day (E) 8.5 and E9.0, and partial E9.5 embryos. To support their utility, we developed sc3D, a tool for reconstructing and exploring three-dimensional 'virtual embryos', which enables the quantitative investigation of regionalized gene expression patterns. Our measurements along the main embryonic axes of the developing neural tube revealed several previously unannotated genes with distinct spatial patterns. We also characterized the conflicting transcriptional identity of 'ectopic' neural tubes that emerge in Tbx6 mutant embryos. Taken together, we present an experimental and computational framework for the spatiotemporal investigation of whole embryonic structures and mutant phenotypes.

An atlas of healthy and injured cell states and niches in the human kidney.

Understanding kidney disease relies on defining the complexity of cell types and states, their associated molecular profiles and interactions within tissue neighbourhoods1. Here we applied multiple single-cell and single-nucleus assays (>400,000 nuclei or cells) and spatial imaging technologies to a broad spectrum of healthy reference kidneys (45 donors) and diseased kidneys (48 patients). This has provided a high-resolution cellular atlas of 51 main cell types, which include rare and previously undescribed cell populations. The multi-omic approach provides detailed transcriptomic profiles, regulatory factors and spatial localizations spanning the entire kidney. We also define 28 cellular states across nephron segments and interstitium that were altered in kidney injury, encompassing cycling, adaptive (successful or maladaptive repair), transitioning and degenerative states. Molecular signatures permitted the localization of these states within injury neighbourhoods using spatial transcriptomics, while large-scale 3D imaging analysis (around 1.2 million neighbourhoods) provided corresponding linkages to active immune responses. These analyses defined biological pathways that are relevant to injury time-course and niches, including signatures underlying epithelial repair that predicted maladaptive states associated with a decline in kidney function. This integrated multimodal spatial cell atlas of healthy and diseased human kidneys represents a comprehensive benchmark of cellular states, neighbourhoods, outcome-associated signatures and publicly available interactive visualizations.

Exposure of iPSC-derived human microglia to brain substrates enables the generation and manipulation of diverse transcriptional states in vitro.

Microglia, the macrophages of the brain parenchyma, are key players in neurodegenerative diseases such as Alzheimer's disease. These cells adopt distinct transcriptional subtypes known as states. Understanding state function, especially in human microglia, has been elusive owing to a lack of tools to model and manipulate these cells. Here, we developed a platform for modeling human microglia transcriptional states in vitro. We found that exposure of human stem-cell-differentiated microglia to synaptosomes, myelin debris, apoptotic neurons or synthetic amyloid-beta fibrils generated transcriptional diversity that mapped to gene signatures identified in human brain microglia, including disease-associated microglia, a state enriched in neurodegenerative diseases. Using a new lentiviral approach, we demonstrated that the transcription factor MITF drives a disease-associated transcriptional signature and a highly phagocytic state. Together, these tools enable the manipulation and functional interrogation of human microglial states in both homeostatic and disease-relevant contexts.

Cortical somatostatin interneuron subtypes form cell-type-specific circuits.

The cardinal classes are a useful simplification of cortical interneuron diversity, but such broad subgroupings gloss over the molecular, morphological, and circuit specificity of interneuron subtypes, most notably among the somatostatin interneuron class. Although there is evidence that this diversity is functionally relevant, the circuit implications of this diversity are unknown. To address this knowledge gap, we designed a series of genetic strategies to target the breadth of somatostatin interneuron subtypes and found that each subtype possesses a unique laminar organization and stereotyped axonal projection pattern. Using these strategies, we examined the afferent and efferent connectivity of three subtypes (two Martinotti and one non-Martinotti) and demonstrated that they possess selective connectivity with intratelecephalic or pyramidal tract neurons. Even when two subtypes targeted the same pyramidal cell type, their synaptic targeting proved selective for particular dendritic compartments. We thus provide evidence that subtypes of somatostatin interneurons form cell-type-specific cortical circuits.

Early Alzheimer's disease pathology in human cortex is associated with a transient phase of distinct cell states.

Cellular perturbations underlying Alzheimer's disease are primarily studied in human postmortem samples and model organisms. Here we generated a single-nucleus atlas from a rare cohort of cortical biopsies from living individuals with varying degrees of Alzheimer's disease pathology. We next performed a systematic cross-disease and cross-species integrative analysis to identify a set of cell states that are specific to early AD pathology. These changes-which we refer to as the Early Cortical Amyloid Response-were prominent in neurons, wherein we identified a transient state of hyperactivity preceding loss of excitatory neurons, which correlated with the selective loss of layer 1 inhibitory neurons. Microglia overexpressing neuroinflammatory-related processes also expanded as AD pathological burden increased. Lastly, both oligodendrocytes and pyramidal neurons upregulated genes associated with amyloid beta production and processing during this early hyperactive phase. Our integrative analysis provides an organizing framework for targeting circuit dysfunction, neuroinflammation, and amyloid production early in AD pathogenesis.