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.

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.

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.

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.

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.

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.

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.

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

Antisense oligonucleotides (ASOs) dosed into cerebrospinal fluid (CSF) distribute broadly throughout the brain and hold the promise of treating myriad brain diseases by modulating RNA. CNS tissue is not routinely biopsied in living individuals, leading to reliance on CSF biomarkers to inform on drug target engagement. Animal models can link CSF biomarkers to brain parenchyma, but our understanding of how individual cells contribute to bulk tissue signal is limited. Here we employed single nucleus transcriptomics on tissue from mice treated with RNase H1 ASOs against Prnp and Malat1 and macaques treated with an ASO against PRNP . Activity was observed in every cell type, though sometimes with substantial differences in magnitude. Single cell RNA count distributions implied target 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 disease-relevant cells in a neuronal disorder.

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.

Dissection of artifactual and confounding glial signatures by single-cell sequencing of mouse and human brain.

A key aspect of nearly all single-cell sequencing experiments is dissociation of intact tissues into single-cell suspensions. While many protocols have been optimized for optimal cell yield, they have often overlooked the effects that dissociation can have on ex vivo gene expression. Here, we demonstrate that use of enzymatic dissociation on brain tissue induces an aberrant ex vivo gene expression signature, most prominently in microglia, which is prevalent in published literature and can substantially confound downstream analyses. To address this issue, we present a rigorously validated protocol that preserves both in vivo transcriptional profiles and cell-type diversity and yield across tissue types and species. We also identify a similar signature in postmortem human brain single-nucleus RNA-sequencing datasets, and show that this signature is induced in freshly isolated human tissue by exposure to elevated temperatures ex vivo. Together, our results provide a methodological solution for preventing artifactual gene expression changes during fresh tissue digestion and a reference for future deeper analysis of the potential confounding states present in postmortem human samples.

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 atlas of mouse cerebellar cortex comprehensively defines cell types.

The cerebellar cortex is a well-studied brain structure with diverse roles in motor learning, coordination, cognition and autonomic regulation. However,  a complete inventory of cerebellar cell types is currently lacking. Here, using recent advances in high-throughput transcriptional profiling1-3, we molecularly define cell types across individual lobules of the adult mouse cerebellum. Purkinje neurons showed considerable regional specialization, with the greatest diversity occurring in the posterior lobules. For several types of cerebellar interneuron, the molecular variation within each type was more continuous, rather than discrete. In particular, for the unipolar brush cells-an interneuron population previously subdivided into discrete populations-the continuous variation in gene expression was associated with a graded continuum of electrophysiological properties. Notably, we found that molecular layer interneurons were composed of two molecularly and functionally distinct types. Both types show a continuum of morphological variation through the thickness of the molecular layer, but electrophysiological recordings revealed marked differences between the two types in spontaneous firing, excitability and electrical coupling. Together, these findings provide a comprehensive cellular atlas of the cerebellar cortex, and outline a methodological and conceptual framework for the integration of molecular, morphological and physiological ontologies for defining brain cell types.

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.