Integrating gene expression and imaging data across Visium capture areas with visiumStitched.

Background: Visium is a widely-used spatially-resolved transcriptomics assay available from 10x Genomics. Standard Visium capture areas (6.5mm by 6.5mm) limit the survey of larger tissue structures, but combining overlapping images and associated gene expression data allow for more complex study designs. Current software can handle nested or partial image overlaps, but is designed for merging up to two capture areas, and cannot account for some technical scenarios related to capture area alignment. Results: We generated Visium data from a postmortem human tissue sample such that two capture areas were partially overlapping and a third one was adjacent. We developed the R/Bioconductor package visiumStitched, which facilitates stitching the images together with Fiji (ImageJ), and constructing SpatialExperiment R objects with the stitched images and gene expression data. visiumStitched constructs an artificial hexagonal array grid which allows seamless downstream analyses such as spatially-aware clustering without discarding data from overlapping spots. Data stitched with visiumStitched can then be interactively visualized with spatialLIBD. Conclusions: visiumStitched provides a simple, but flexible framework to handle various multi-capture area study design scenarios. Specifically, it resolves a data processing step without disrupting analysis workflows and without discarding data from overlapping spots. visiumStiched relies on affine transformations by Fiji, which have limitations and are less accurate when aligning against an atlas or other situations. visiumStiched provides an easy-to-use solution which expands possibilities for designing multi-capture area study designs.

A data-driven single-cell and spatial transcriptomic map of the human prefrontal cortex.

The molecular organization of the human neocortex historically has been studied in the context of its histological layers. However, emerging spatial transcriptomic technologies have enabled unbiased identification of transcriptionally defined spatial domains that move beyond classic cytoarchitecture. We used the Visium spatial gene expression platform to generate a data-driven molecular neuroanatomical atlas across the anterior-posterior axis of the human dorsolateral prefrontal cortex. Integration with paired single-nucleus RNA-sequencing data revealed distinct cell type compositions and cell-cell interactions across spatial domains. Using PsychENCODE and publicly available data, we mapped the enrichment of cell types and genes associated with neuropsychiatric disorders to discrete spatial domains.

Transcriptomic characterization of human lateral septum neurons reveals conserved and divergent marker genes across species.

The lateral septum (LS) is a midline, subcortical structure, which regulates social behaviors that are frequently impaired in neurodevelopmental disorders including schizophrenia and autism spectrum disorder. Mouse studies have identified neuronal populations within the LS that express a variety of molecular markers, including vasopressin receptor, oxytocin receptor, and corticotropin releasing hormone receptor, which control specific facets of social behavior. Despite its critical role in regulating social behavior and notable gene expression patterns, comprehensive molecular profiling of the human LS has not been performed. Here, we conducted single nucleus RNA-sequencing (snRNA-seq) to generate the first transcriptomic profiles of the human LS using postmortem human brain tissue samples from 3 neurotypical donors. Our analysis identified 5 transcriptionally distinct neuronal cell types within the human LS that are enriched for TRPC4, the gene encoding Trp-related protein 4. Differential expression analysis revealed a distinct LS neuronal cell type that is enriched for OPRM1, the gene encoding the µ-opioid receptor. Leveraging recently generated mouse LS snRNA-seq datasets, we conducted a cross-species analysis. Our results demonstrate that TRPC4 enrichment in the LS is highly conserved between human and mouse, while FREM2, which encodes FRAS1 related extracellular matrix protein 2, is enriched only in the human LS. Together, these results highlight transcriptional heterogeneity of the human LS, and identify robust marker genes for the human LS.

An integrated single-nucleus and spatial transcriptomics atlas reveals the molecular landscape of the human hippocampus.

The hippocampus contains many unique cell types, which serve the structure's specialized functions, including learning, memory and cognition. These cells have distinct spatial topography, morphology, physiology, and connectivity, highlighting the need for transcriptome-wide profiling strategies that retain cytoarchitectural organization. Here, we generated spatially-resolved transcriptomics (SRT) and single-nucleus RNA-sequencing (snRNA-seq) data from adjacent tissue sections of the anterior human hippocampus across ten adult neurotypical donors. We defined molecular profiles for hippocampal cell types and spatial domains. Using non-negative matrix factorization and transfer learning, we integrated these data to define gene expression patterns within the snRNA-seq data and infer the expression of these patterns in the SRT data. With this approach, we leveraged existing rodent datasets that feature information on circuit connectivity and neural activity induction to make predictions about axonal projection targets and likelihood of ensemble recruitment in spatially-defined cellular populations of the human hippocampus. Finally, we integrated genome-wide association studies with transcriptomic data to identify enrichment of genetic components for neurodevelopmental, neuropsychiatric, and neurodegenerative disorders across cell types, spatial domains, and gene expression patterns of the human hippocampus. To make this comprehensive molecular atlas accessible to the scientific community, both raw and processed data are freely available, including through interactive web applications.

Spatial domain detection using contrastive self-supervised learning for spatial multi-omics technologies.

Recent advances in spatially-resolved single-omics and multi-omics technologies have led to the emergence of computational tools to detect or predict spatial domains. Additionally, histological images and immunofluorescence (IF) staining of proteins and cell types provide multiple perspectives and a more complete understanding of tissue architecture. Here, we introduce Proust, a scalable tool to predict discrete domains using spatial multi-omics data by combining the low-dimensional representation of biological profiles based on graph-based contrastive self-supervised learning. Our scalable method integrates multiple data modalities, such as RNA, protein, and H&E images, and predicts spatial domains within tissue samples. Through the integration of multiple modalities, Proust consistently demonstrates enhanced accuracy in detecting spatial domains, as evidenced across various benchmark datasets and technological platforms.

The gene expression landscape of the human locus coeruleus revealed by single-nucleus and spatially-resolved transcriptomics.

Norepinephrine (NE) neurons in the locus coeruleus (LC) make long-range projections throughout the central nervous system, playing critical roles in arousal and mood, as well as various components of cognition including attention, learning, and memory. The LC-NE system is also implicated in multiple neurological and neuropsychiatric disorders. Importantly, LC-NE neurons are highly sensitive to degeneration in both Alzheimer's and Parkinson's disease. Despite the clinical importance of the brain region and the prominent role of LC-NE neurons in a variety of brain and behavioral functions, a detailed molecular characterization of the LC is lacking. Here, we used a combination of spatially-resolved transcriptomics and single-nucleus RNA-sequencing to characterize the molecular landscape of the LC region and the transcriptomic profile of LC-NE neurons in the human brain. We provide a freely accessible resource of these data in web-accessible and downloadable formats.

TrkB-dependent regulation of molecular signaling across septal cell types.

The lateral septum (LS), a GABAergic structure located in the basal forebrain, is implicated in social behavior, learning, and memory. We previously demonstrated that expression of tropomyosin kinase receptor B (TrkB) in LS neurons is required for social novelty recognition. To better understand molecular mechanisms by which TrkB signaling controls behavior, we locally knocked down TrkB in LS and used bulk RNA-sequencing to identify changes in gene expression downstream of TrkB. TrkB knockdown induces upregulation of genes associated with inflammation and immune responses, and downregulation of genes associated with synaptic signaling and plasticity. Next, we generated one of the first atlases of molecular profiles for LS cell types using single nucleus RNA-sequencing (snRNA-seq). We identified markers for the septum broadly, and the LS specifically, as well as for all neuronal cell types. We then investigated whether the differentially expressed genes (DEGs) induced by TrkB knockdown map to specific LS cell types. Enrichment testing identified that downregulated DEGs are broadly expressed across neuronal clusters. Enrichment analyses of these DEGs demonstrated that downregulated genes are uniquely expressed in the LS, and associated with either synaptic plasticity or neurodevelopmental disorders. Upregulated genes are enriched in LS microglia, associated with immune response and inflammation, and linked to both neurodegenerative disease and neuropsychiatric disorders. In addition, many of these genes are implicated in regulating social behaviors. In summary, the findings implicate TrkB signaling in the LS as a critical regulator of gene networks associated with psychiatric disorders that display social deficits, including schizophrenia and autism, and with neurodegenerative diseases, including Alzheimer's.

Spatially-resolved transcriptomics of human dentate gyrus across postnatal lifespan reveals heterogeneity in markers for proliferation, extracellular matrix, and neuroinflammation.

The dentate gyrus of the anterior hippocampus is important for many human cognitive functions, including regulation of learning, memory, and mood. However, the postnatal development and aging of the dentate gyrus throughout the human lifespan has yet to be fully characterized in the same molecular and spatial detail as other species. Here, we generated a spatially-resolved molecular atlas of the dentate gyrus in postmortem human tissue using the 10x Genomics Visium platform to retain extranuclear transcripts and identify changes in molecular topography across the postnatal lifespan. We found enriched expression of extracellular matrix markers during infancy and increased expression of GABAergic cell-type markers GAD1, LAMP5, and CCK after infancy. While we identified a conserved gene signature for mouse neuroblasts in the granule cell layer (GCL), many of those genes are not specific to the GCL, and we found no evidence of signatures for other granule cell lineage stages at the GCL post-infancy. We identified a wide-spread hippocampal aging signature and an age-dependent increase in neuroinflammation associated genes. Our findings suggest major changes to the putative neurogenic niche after infancy and identify molecular foci of brain aging in glial and neuropil enriched tissue.

Data-driven identification of total RNA expression genes for estimation of RNA abundance in heterogeneous cell types highlighted in brain tissue.

We define and identify a new class of control genes for next-generation sequencing called total RNA expression genes (TREGs), which correlate with total RNA abundance in cell types of different sizes and transcriptional activity. We provide a data-driven method to identify TREGs from single-cell RNA sequencing data, allowing the estimation of total amount of RNA when restricted to quantifying a limited number of genes. We demonstrate our method in postmortem human brain using multiplex single-molecule fluorescent in situ hybridization and compare candidate TREGs against classic housekeeping genes. We identify AKT3 as a top TREG across five brain regions.

TrkB-dependent regulation of molecular signaling across septal cell types.

The lateral septum (LS), a GABAergic structure located in the basal forebrain, is implicated in social behavior, learning and memory. We previously demonstrated that expression of tropomyosin kinase receptor B (TrkB) in LS neurons is required for social novelty recognition. To better understand molecular mechanisms by which TrkB signaling controls behavior, we locally knocked down TrkB in LS and used bulk RNA-sequencing to identify changes in gene expression downstream of TrkB. TrkB knockdown induces upregulation of genes associated with inflammation and immune responses, and downregulation of genes associated with synaptic signaling and plasticity. Next, we generated one of the first atlases of molecular profiles for LS cell types using single nucleus RNA-sequencing (snRNA-seq). We identified markers for the septum broadly, and the LS specifically, as well as for all neuronal cell types. We then investigated whether the differentially expressed genes (DEGs) induced by TrkB knockdown map to specific LS cell types. Enrichment testing identified that downregulated DEGs are broadly expressed across neuronal clusters. Enrichment analyses of these DEGs demonstrated that downregulated genes are uniquely expressed in the LS, and associated with either synaptic plasticity or neurodevelopmental disorders. Upregulated genes are enriched in LS microglia, associated with immune response and inflammation, and linked to both neurodegenerative disease and neuropsychiatric disorders. In addition, many of these genes are implicated in regulating social behaviors. In summary, the findings implicate TrkB signaling in the LS as a critical regulator of gene networks associated with psychiatric disorders that display social deficits, including schizophrenia and autism, and with neurodegenerative diseases, including Alzheimer's.

Integrated single cell and unsupervised spatial transcriptomic analysis defines molecular anatomy of the human dorsolateral prefrontal cortex.

Generation of a molecular neuroanatomical map of the human prefrontal cortex reveals novel spatial domains and cell-cell interactions relevant for psychiatric disease. The molecular organization of the human neocortex has been historically studied in the context of its histological layers. However, emerging spatial transcriptomic technologies have enabled unbiased identification of transcriptionally-defined spatial domains that move beyond classic cytoarchitecture. Here we used the Visium spatial gene expression platform to generate a data-driven molecular neuroanatomical atlas across the anterior-posterior axis of the human dorsolateral prefrontal cortex (DLPFC). Integration with paired single nucleus RNA-sequencing data revealed distinct cell type compositions and cell-cell interactions across spatial domains. Using PsychENCODE and publicly available data, we map the enrichment of cell types and genes associated with neuropsychiatric disorders to discrete spatial domains. Finally, we provide resources for the scientific community to explore these integrated spatial and single cell datasets at research.libd.org/spatialDLPFC/.

SUFI: an automated approach to spectral unmixing of fluorescent multiplex images captured in mouse and post-mortem human brain tissues.

BACKGROUND: Multispectral fluorescence imaging coupled with linear unmixing is a form of image data collection and analysis that allows for measuring multiple molecular signals in a single biological sample. Multiple fluorescent dyes, each measuring a unique molecule, are simultaneously measured and subsequently "unmixed" to provide a read-out for each molecular signal. This strategy allows for measuring highly multiplexed signals in a single data capture session, such as multiple proteins or RNAs in tissue slices or cultured cells, but can often result in mixed signals and bleed-through problems across dyes. Existing spectral unmixing algorithms are not optimized for challenging biological specimens such as post-mortem human brain tissue, and often require manual intervention to extract spectral signatures. We therefore developed an intuitive, automated, and flexible package called SUFI: spectral unmixing of fluorescent images. RESULTS: This package unmixes multispectral fluorescence images by automating the extraction of spectral signatures using vertex component analysis, and then performs one of three unmixing algorithms derived from remote sensing. We evaluate these remote sensing algorithms' performances on four unique biological datasets and compare the results to unmixing results obtained using ZEN Black software (Zeiss). We lastly integrate our unmixing pipeline into the computational tool dotdotdot, which is used to quantify individual RNA transcripts at single cell resolution in intact tissues and perform differential expression analysis, and thereby provide an end-to-end solution for multispectral fluorescence image analysis and quantification. CONCLUSIONS: In summary, we provide a robust, automated pipeline to assist biologists with improved spectral unmixing of multispectral fluorescence images.

The basolateral amygdala to lateral septum circuit is critical for regulating social novelty in mice.

The lateral septum (LS) is a basal forebrain GABAergic region that is implicated in social novelty. However, the neural circuits and cell signaling pathways that converge on the LS to mediate social behaviors aren't well understood. Multiple lines of evidence suggest that signaling of brain-derived neurotrophic factor (BDNF) through its receptor TrkB plays important roles in social behavior. BDNF is not locally produced in LS, but we demonstrate that nearly all LS GABAergic neurons express TrkB. Local TrkB knock-down in LS neurons decreased social novelty recognition and reduced recruitment of neural activity in LS neurons in response to social novelty. Since BDNF is not synthesized in LS, we investigated which inputs to LS could serve as potential BDNF sources for controlling social novelty recognition. We demonstrate that selectively ablating inputs to LS from the basolateral amygdala (BLA), but not from ventral CA1 (vCA1), impairs social novelty recognition. Moreover, depleting BDNF selectively in BLA-LS projection neurons phenocopied the decrease in social novelty recognition caused by either local LS TrkB knockdown or ablation of BLA-LS inputs. These data support the hypothesis that BLA-LS projection neurons serve as a critical source of BDNF for activating TrkB signaling in LS neurons to control social novelty recognition.

Influence of Alzheimer's disease related neuropathology on local microenvironment gene expression in the human inferior temporal cortex.

Neuropathological lesions in the brains of individuals affected with neurodegenerative disorders are hypothesized to trigger molecular and cellular processes that disturb homeostasis of local microenvironments. Here, we applied the 10x Genomics Visium Spatial Proteogenomics (Visium-SPG) platform, which couples spatial gene expression with immunofluorescence protein co-detection, to evaluate its ability to quantify changes in spatial gene expression with respect to amyloid-ß (Aß) and hyperphosphorylated tau (pTau) pathology in post-mortem human brain tissue from individuals with Alzheimer's disease (AD). We identified transcriptomic signatures associated with proximity to Aß in the human inferior temporal cortex (ITC) during late-stage AD, which we further investigated at cellular resolution with combined immunofluorescence and single molecule fluorescent in situ hybridization (smFISH). The study provides a data analysis workflow for Visium-SPG, and the data represent a proof-of-principal for the power of multi-omic profiling in identifying changes in molecular dynamics that are spatially-associated with pathology in the human brain. We provide the scientific community with web-based, interactive resources to access the datasets of the spatially resolved AD-related transcriptomes at https://research.libd.org/Visium_SPG_AD/.

VistoSeg: Processing utilities for high-resolution images for spatially resolved transcriptomics data.

Spatially resolved transcriptomics (SRT) is a growing field that links gene expression to anatomical context. SRT approaches that use next-generation sequencing (NGS) combine RNA sequencing with histological or fluorescent imaging to generate spatial maps of gene expression in intact tissue sections. These technologies directly couple gene expression measurements with high-resolution histological or immunofluorescent images that contain rich morphological information about the tissue under study. While broad access to NGS-based spatial transcriptomic technology is now commercially available through the Visium platform from the vendor 10× Genomics, computational tools for extracting image-derived metrics for integration with gene expression data remain limited. We developed VistoSeg as a MATLAB pipeline to process, analyze and interactively visualize the high-resolution images generated in the Visium platform. VistoSeg outputs can be easily integrated with accompanying transcriptomic data to facilitate downstream analyses in common programing languages including R and Python. VistoSeg provides user-friendly tools for integrating image-derived metrics from histological and immunofluorescent images with spatially resolved gene expression data. Integration of this data enhances the ability to understand the transcriptional landscape within tissue architecture. VistoSeg is freely available at http://research.libd.org/VistoSeg/.

Decoding Shared Versus Divergent Transcriptomic Signatures Across Cortico-Amygdala Circuitry in PTSD and Depressive Disorders.

OBJECTIVE: Posttraumatic stress disorder (PTSD) is a debilitating neuropsychiatric disease that is highly comorbid with major depressive disorder (MDD) and bipolar disorder. The overlap in symptoms is hypothesized to stem from partially shared genetics and underlying neurobiological mechanisms. To delineate conservation between transcriptional patterns across PTSD and MDD, the authors examined gene expression in the human cortex and amygdala in these disorders. METHODS: RNA sequencing was performed in the postmortem brain of two prefrontal cortex regions and two amygdala regions from donors diagnosed with PTSD (N=107) or MDD (N=109) as well as from neurotypical donors (N=109). RESULTS: The authors identified a limited number of differentially expressed genes (DEGs) specific to PTSD, with nearly all mapping to cortical versus amygdala regions. PTSD-specific DEGs were enriched in gene sets associated with downregulated immune-related pathways and microglia as well as with subpopulations of GABAergic inhibitory neurons. While a greater number of DEGs associated with MDD were identified, most overlapped with PTSD, and only a few were MDD specific. The authors used weighted gene coexpression network analysis as an orthogonal approach to confirm the observed cellular and molecular associations. CONCLUSIONS: These findings provide supporting evidence for involvement of decreased immune signaling and neuroinflammation in MDD and PTSD pathophysiology, and extend evidence that GABAergic neurons have functional significance in PTSD.

spatialLIBD: an R/Bioconductor package to visualize spatially-resolved transcriptomics data.

BACKGROUND: Spatially-resolved transcriptomics has now enabled the quantification of high-throughput and transcriptome-wide gene expression in intact tissue while also retaining the spatial coordinates. Incorporating the precise spatial mapping of gene activity advances our understanding of intact tissue-specific biological processes. In order to interpret these novel spatial data types, interactive visualization tools are necessary. RESULTS: We describe spatialLIBD, an R/Bioconductor package to interactively explore spatially-resolved transcriptomics data generated with the 10x Genomics Visium platform. The package contains functions to interactively access, visualize, and inspect the observed spatial gene expression data and data-driven clusters identified with supervised or unsupervised analyses, either on the user's computer or through a web application. CONCLUSIONS: spatialLIBD is available at https://bioconductor.org/packages/spatialLIBD . It is fully compatible with SpatialExperiment and the Bioconductor ecosystem. Its functionality facilitates analyzing and interactively exploring spatially-resolved data from the Visium platform.

Dissecting transcriptomic signatures of neuronal differentiation and maturation using iPSCs.

Human induced pluripotent stem cells (hiPSCs) are a powerful model of neural differentiation and maturation. We present a hiPSC transcriptomics resource on corticogenesis from 5 iPSC donor and 13 subclonal lines across 9 time points over 5 broad conditions: self-renewal, early neuronal differentiation, neural precursor cells (NPCs), assembled rosettes, and differentiated neuronal cells. We identify widespread changes in the expression of both individual features and global patterns of transcription. We next demonstrate that co-culturing human NPCs with rodent astrocytes results in mutually synergistic maturation, and that cell type-specific expression data can be extracted using only sequencing read alignments without cell sorting. We lastly adapt a previously generated RNA deconvolution approach to single-cell expression data to estimate the relative neuronal maturity of iPSC-derived neuronal cultures and human brain tissue. Using many public datasets, we demonstrate neuronal cultures are maturationally heterogeneous but contain subsets of neurons more mature than previously observed.