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.

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.

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.

Spatial genomics enables multi-modal study of clonal heterogeneity in tissues.

The state and behaviour of a cell can be influenced by both genetic and environmental factors. In particular, tumour progression is determined by underlying genetic aberrations1-4 as well as the makeup of the tumour microenvironment5,6. Quantifying the contributions of these factors requires new technologies that can accurately measure the spatial location of genomic sequence together with phenotypic readouts. Here we developed slide-DNA-seq, a method for capturing spatially resolved DNA sequences from intact tissue sections. We demonstrate that this method accurately preserves local tumour architecture and enables the de novo discovery of distinct tumour clones and their copy number alterations. We then apply slide-DNA-seq to a mouse model of metastasis and a primary human cancer, revealing that clonal populations are confined to distinct spatial regions. Moreover, through integration with spatial transcriptomics, we uncover distinct sets of genes that are associated with clone-specific genetic aberrations, the local tumour microenvironment, or both. Together, this multi-modal spatial genomics approach provides a versatile platform for quantifying how cell-intrinsic and cell-extrinsic factors contribute to gene expression, protein abundance and other cellular phenotypes.

A transcriptomic and epigenomic cell atlas of the mouse primary motor cortex.

Single-cell transcriptomics can provide quantitative molecular signatures for large, unbiased samples of the diverse cell types in the brain1-3. With the proliferation of multi-omics datasets, a major challenge is to validate and integrate results into a biological understanding of cell-type organization. Here we generated transcriptomes and epigenomes from more than 500,000 individual cells in the mouse primary motor cortex, a structure that has an evolutionarily conserved role in locomotion. We developed computational and statistical methods to integrate multimodal data and quantitatively validate cell-type reproducibility. The resulting reference atlas-containing over 56 neuronal cell types that are highly replicable across analysis methods, sequencing technologies and modalities-is a comprehensive molecular and genomic account of the diverse neuronal and non-neuronal cell types in the mouse primary motor cortex. The atlas includes a population of excitatory neurons that resemble pyramidal cells in layer 4 in other cortical regions4. We further discovered thousands of concordant marker genes and gene regulatory elements for these cell types. Our results highlight the complex molecular regulation of cell types in the brain and will directly enable the design of reagents to target specific cell types in the mouse primary motor cortex for functional analysis.

A single-cell map of antisense oligonucleotide activity in the brain.

Antisense oligonucleotides (ASOs) dosed into cerebrospinal fluid (CSF) distribute broadly throughout the central nervous system (CNS). By modulating RNA, they hold the promise of targeting root molecular causes of disease and hold potential to treat myriad CNS disorders. Realization of this potential requires that ASOs must be active in the disease-relevant cells, and ideally, that monitorable biomarkers also reflect ASO activity in these cells. The biodistribution and activity of such centrally delivered ASOs have been deeply characterized in rodent and non-human primate (NHP) models, but usually only in bulk tissue, limiting our understanding of the distribution of ASO activity across individual cells and across diverse CNS cell types. Moreover, in human clinical trials, target engagement is usually monitorable only in a single compartment, CSF. We sought a deeper understanding of how individual cells and cell types contribute to bulk tissue signal in the CNS, and how these are linked to CSF biomarker outcomes. We employed single nucleus transcriptomics on tissue from mice treated with RNase H1 ASOs against Prnp and Malat1 and NHPs treated with an ASO against PRNP. Pharmacologic activity was observed in every cell type, though sometimes with substantial differences in magnitude. Single cell RNA count distributions implied target RNA suppression in every single sequenced cell, rather than intense knockdown in only some cells. Duration of action up to 12 weeks post-dose differed across cell types, being shorter in microglia than in neurons. Suppression in neurons was generally similar to, or more robust than, the bulk tissue. In macaques, PrP in CSF was lowered 40% in conjunction with PRNP knockdown across all cell types including neurons, arguing that a CSF biomarker readout is likely to reflect ASO pharmacodynamic effect in disease-relevant cells in a neuronal disorder. Our results provide a reference dataset for ASO activity distribution in the CNS and establish single nucleus sequencing as a method for evaluating cell type specificity of oligonucleotide therapeutics and other modalities.

Antisense oligonucleotide (ASO) drugs are a type of chemically modified DNA that can be injected into cerebrospinal fluid in order to enter brain cells and reduce the amount of RNA from a specific gene. The brain is a complex mixture of hundreds of billions of cells. When an ASO lowers a target gene's RNA by 50%, is that a 50% reduction in 100% of cells, or a 100% reduction in 50% of cells? Are the many different cell types of the brain affected equally? This new study uses single cell RNA sequencing to answer these questions, finding that ASOs are broadly active across cell types and individual cells, and linking reduction of target protein in cerebrospinal fluid to disease-relevant cells.

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.

The cell type composition of the adult mouse brain revealed by single cell and spatial genomics.

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-seq with Slide-seq-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 signaling, 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 (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.

Slide-tags: scalable, single-nucleus barcoding for multi-modal 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 constructed. Missing from these measurements, however, is the ability to routinely and easily spatially localise 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 input into a wide variety of single-nucleus profiling assays. Application of Slide-tags to the mouse hippocampus positioned nuclei at less than 10 micron spatial resolution, and delivered whole-transcriptome data that was indistinguishable in quality from ordinary snRNA-seq. 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 contextualised receptor-ligand interactions driving B-cell maturation in lymphoid tissue. A major benefit of Slide-tags is that it is easily adaptable to virtually any single-cell measurement technology. As proof of principle, we performed multiomic measurements of open chromatin, RNA, and T-cell receptor sequences in the same cells from metastatic melanoma. We identified spatially distinct tumour subpopulations to be differentially infiltrated by an expanded T-cell clone and undergoing cell state transition driven by spatially clustered accessible transcription factor motifs. Slide-tags offers a universal platform for importing the compendium of established single-cell measurements into the spatial genomics repertoire.

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.

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.

Candelabrum cells are ubiquitous cerebellar cortex interneurons with specialized circuit properties.

To understand how the cerebellar cortex transforms mossy fiber (MF) inputs into Purkinje cell (PC) outputs, it is vital to delineate the elements of this circuit. Candelabrum cells (CCs) are enigmatic interneurons of the cerebellar cortex that have been identified based on their morphology, but their electrophysiological properties, synaptic connections and function remain unknown. Here, we clarify these properties using electrophysiology, single-nucleus RNA sequencing, in situ hybridization and serial electron microscopy in mice. We find that CCs are the most abundant PC layer interneuron. They are GABAergic, molecularly distinct and present in all cerebellar lobules. Their high resistance renders CC firing highly sensitive to synaptic inputs. CCs are excited by MFs and granule cells and are strongly inhibited by PCs. CCs in turn primarily inhibit molecular layer interneurons, which leads to PC disinhibition. Thus, inputs, outputs and local signals converge onto CCs to allow them to assume a unique role in controlling cerebellar output.

Single-cell genomic profiling of human dopamine neurons identifies a population that selectively degenerates in Parkinson's disease.

The loss of dopamine (DA) neurons within the substantia nigra pars compacta (SNpc) is a defining pathological hallmark of Parkinson's disease (PD). Nevertheless, the molecular features associated with DA neuron vulnerability have not yet been fully identified. Here, we developed a protocol to enrich and transcriptionally profile DA neurons from patients with PD and matched controls, sampling a total of 387,483 nuclei, including 22,048 DA neuron profiles. We identified ten populations and spatially localized each within the SNpc using Slide-seq. A single subtype, marked by the expression of the gene AGTR1 and spatially confined to the ventral tier of SNpc, was highly susceptible to loss in PD and showed the strongest upregulation of targets of TP53 and NR2F2, nominating molecular processes associated with degeneration. This same vulnerable population was specifically enriched for the heritable risk associated with PD, highlighting the importance of cell-intrinsic processes in determining the differential vulnerability of DA neurons to PD-associated degeneration.