Semla: a versatile toolkit for spatially resolved transcriptomics analysis and visualization.

SUMMARY: Spatially resolved transcriptomics technologies generate gene expression data with retained positional information from a tissue section, often accompanied by a corresponding histological image. Computational tools should make it effortless to incorporate spatial information into data analyses and present analysis results in their histological context. Here, we present semla, an R package for processing, analysis, and visualization of spatially resolved transcriptomics data generated by the Visium platform, that includes interactive web applications for data exploration and tissue annotation. AVAILABILITY AND IMPLEMENTATION: The R package semla is available on GitHub (https://github.com/ludvigla/semla), under the MIT License, and deposited on Zenodo (https://doi.org/10.5281/zenodo.8321645). Documentation and tutorials with detailed descriptions of usage can be found at https://ludvigla.github.io/semla/.

Spatial transcriptomics of human placentas reveal distinct RNA patterns associated with morphology and preeclampsia.

Spatial transcriptomics (ST) maps RNA level patterns within a tissue. This technology has not been previously applied to human placental tissue. We demonstrate analysis of human placental samples with ST. Unsupervised clustering revealed that distinct RNA patterns were found corresponding to different morphological structures. Additionally, when focusing upon terminal villi and hemoglobin associated structures, RNA levels differed between placentas from full term healthy pregnancies and those complicated by preeclampsia. The results from this study can provide a benchmark for future ST studies in placenta.

An integrated single cell and spatial transcriptomic map of human white adipose tissue.

To date, single-cell studies of human white adipose tissue (WAT) have been based on small cohort sizes and no cellular consensus nomenclature exists. Herein, we performed a comprehensive meta-analysis of publicly available and newly generated single-cell, single-nucleus, and spatial transcriptomic results from human subcutaneous, omental, and perivascular WAT. Our high-resolution map is built on data from ten studies and allowed us to robustly identify >60 subpopulations of adipocytes, fibroblast and adipogenic progenitors, vascular, and immune cells. Using these results, we deconvolved spatial and bulk transcriptomic data from nine additional cohorts to provide spatial and clinical dimensions to the map. This identified cell-cell interactions as well as relationships between specific cell subtypes and insulin resistance, dyslipidemia, adipocyte volume, and lipolysis upon long-term weight changes. Altogether, our meta-map provides a rich resource defining the cellular and microarchitectural landscape of human WAT and describes the associations between specific cell types and metabolic states.

Solid-phase capture and profiling of open chromatin by spatial ATAC.

Current methods for epigenomic profiling are limited in their ability to obtain genome-wide information with spatial resolution. We introduce spatial ATAC, a method that integrates transposase-accessible chromatin profiling in tissue sections with barcoded solid-phase capture to perform spatially resolved epigenomics. We show that spatial ATAC enables the discovery of the regulatory programs underlying spatial gene expression during mouse organogenesis, lineage differentiation and in human pathology.

Integrated single cell and spatial transcriptomics reveal autoreactive differentiated B cells in joints of early rheumatoid arthritis.

B cells play a significant role in established Rheumatoid Arthritis (RA). However, it is unclear to what extent differentiated B cells are present in joint tissue already at the onset of disease. Here, we studied synovial biopsies (n = 8) captured from untreated patients at time of diagnosis. 3414 index-sorted B cells underwent RNA sequencing and paired tissue pieces were subjected to spatial transcriptomics (n = 4). We performed extensive bioinformatics analyses to dissect the local B cell composition. Select plasma cell immunoglobulin sequences were expressed as monoclonal antibodies and tested by ELISA. Memory and plasma cells were found irrespective of autoantibody status of the patients. Double negative memory B cells were prominent, but did not display a distinct transcriptional profile. The tissue architecture implicate both local B cell maturation via T cell help and plasma cell survival niches with a strong CXCL12-CXCR4 axis. The immunoglobulin sequence analyses revealed clonality between the memory B and plasma cell pools further supporting local maturation. One of the plasma cell-derived antibodies displayed citrulline autoreactivity, demonstrating local autoreactive plasma cell differentiation in joint biopsies captured from untreated early RA. Hence, plasma cell niches are not a consequence of chronic inflammation, but are already present at the time of diagnosis.

Three-dimensional spatial transcriptomics uncovers cell type localizations in the human rheumatoid arthritis synovium.

The inflamed rheumatic joint is a highly heterogeneous and complex tissue with dynamic recruitment and expansion of multiple cell types that interact in multifaceted ways within a localized area. Rheumatoid arthritis synovium has primarily been studied either by immunostaining or by molecular profiling after tissue homogenization. Here, we use Spatial Transcriptomics, where tissue-resident RNA is spatially labeled in situ with barcodes in a transcriptome-wide fashion, to study local tissue interactions at the site of chronic synovial inflammation. We report comprehensive spatial RNA-Seq data coupled to cell type-specific localization patterns at and around organized structures of infiltrating leukocyte cells in the synovium. Combining morphological features and high-throughput spatially resolved transcriptomics may be able to provide higher statistical power and more insights into monitoring disease severity and treatment-specific responses in seropositive and seronegative rheumatoid arthritis.

Spatial mapping reveals human adipocyte subpopulations with distinct sensitivities to insulin.

The contribution of cellular heterogeneity and architecture to white adipose tissue (WAT) function is poorly understood. Herein, we combined spatially resolved transcriptional profiling with single-cell RNA sequencing and image analyses to map human WAT composition and structure. This identified 18 cell classes with unique propensities to form spatially organized homo- and heterotypic clusters. Of these, three constituted mature adipocytes that were similar in size, but distinct in their spatial arrangements and transcriptional profiles. Based on marker genes, we termed these AdipoLEP, AdipoPLIN, and AdipoSAA. We confirmed, in independent datasets, that their respective gene profiles associated differently with both adipocyte and whole-body insulin sensitivity. Corroborating our observations, insulin stimulation in vivo by hyperinsulinemic-euglycemic clamp showed that only AdipoPLIN displayed a transcriptional response to insulin. Altogether, by mining this multimodal resource we identify that human WAT is composed of three classes of mature adipocytes, only one of which is insulin responsive.

Divergent clonal differentiation trajectories establish CD8+ memory T cell heterogeneity during acute viral infections in humans.

The CD8+ T cell response to an antigen is composed of many T cell clones with unique T cell receptors, together forming a heterogeneous repertoire of effector and memory cells. How individual T cell clones contribute to this heterogeneity throughout immune responses remains largely unknown. In this study, we longitudinally track human CD8+ T cell clones expanding in response to yellow fever virus (YFV) vaccination at the single-cell level. We observed a drop in clonal diversity in blood from the acute to memory phase, suggesting that clonal selection shapes the circulating memory repertoire. Clones in the memory phase display biased differentiation trajectories along a gradient from stem cell to terminally differentiated effector memory fates. In secondary responses, YFV- and influenza-specific CD8+ T cell clones are poised to recapitulate skewed differentiation trajectories. Collectively, we show that the sum of distinct clonal phenotypes results in the multifaceted human T cell response to acute viral infections.

Identification and transfer of spatial transcriptomics signatures for cancer diagnosis.

BACKGROUND: Distinguishing ductal carcinoma in situ (DCIS) from invasive ductal carcinoma (IDC) regions in clinical biopsies constitutes a diagnostic challenge. Spatial transcriptomics (ST) is an in situ capturing method, which allows quantification and visualization of transcriptomes in individual tissue sections. In the past, studies have shown that breast cancer samples can be used to study their transcriptomes with spatial resolution in individual tissue sections. Previously, supervised machine learning methods were used in clinical studies to predict the clinical outcomes for cancer types. METHODS: We used four publicly available ST breast cancer datasets from breast tissue sections annotated by pathologists as non-malignant, DCIS, or IDC. We trained and tested a machine learning method (support vector machine) based on the expert annotation as well as based on automatic selection of cell types by their transcriptome profiles. RESULTS: We identified expression signatures for expert annotated regions (non-malignant, DCIS, and IDC) and build machine learning models. Classification results for 798 expression signature transcripts showed high coincidence with the expert pathologist annotation for DCIS (100%) and IDC (96%). Extending our analysis to include all 25,179 expressed transcripts resulted in an accuracy of 99% for DCIS and 98% for IDC. Further, classification based on an automatically identified expression signature covering all ST spots of tissue sections resulted in prediction accuracy of 95% for DCIS and 91% for IDC. CONCLUSIONS: This concept study suggest that the ST signatures learned from expert selected breast cancer tissue sections can be used to identify breast cancer regions in whole tissue sections including regions not trained on. Furthermore, the identified expression signatures can classify cancer regions in tissue sections not used for training with high accuracy. Expert-generated but even automatically generated cancer signatures from ST data might be able to classify breast cancer regions and provide clinical decision support for pathologists in the future.

Exploring inflammatory signatures in arthritic joint biopsies with Spatial Transcriptomics.

Lately it has become possible to analyze transcriptomic profiles in tissue sections with retained cellular context. We aimed to explore synovial biopsies from rheumatoid arthritis (RA) and spondyloarthritis (SpA) patients, using Spatial Transcriptomics (ST) as a proof of principle approach for unbiased mRNA studies at the site of inflammation in these chronic inflammatory diseases. Synovial tissue biopsies from affected joints were studied with ST. The transcriptome data was subjected to differential gene expression analysis (DEA), pathway analysis, immune cell type identification using Xcell analysis and validation with immunohistochemistry (IHC). The ST technology allows selective analyses on areas of interest, thus we analyzed morphologically distinct areas of mononuclear cell infiltrates. The top differentially expressed genes revealed an adaptive immune response profile and T-B cell interactions in RA, while in SpA, the profiles implicate functions associated with tissue repair. With spatially resolved gene expression data, overlaid on high-resolution histological images, we digitally portrayed pre-selected cell types in silico. The RA displayed an overrepresentation of central memory T cells, while in SpA effector memory T cells were most prominent. Consequently, ST allows for deeper understanding of cellular mechanisms and diversity in tissues from chronic inflammatory diseases.

A Spatiotemporal Organ-Wide Gene Expression and Cell Atlas of the Developing Human Heart.

The process of cardiac morphogenesis in humans is incompletely understood. Its full characterization requires a deep exploration of the organ-wide orchestration of gene expression with a single-cell spatial resolution. Here, we present a molecular approach that reveals the comprehensive transcriptional landscape of cell types populating the embryonic heart at three developmental stages and that maps cell-type-specific gene expression to specific anatomical domains. Spatial transcriptomics identified unique gene profiles that correspond to distinct anatomical regions in each developmental stage. Human embryonic cardiac cell types identified by single-cell RNA sequencing confirmed and enriched the spatial annotation of embryonic cardiac gene expression. In situ sequencing was then used to refine these results and create a spatial subcellular map for the three developmental phases. Finally, we generated a publicly available web resource of the human developing heart to facilitate future studies on human cardiogenesis.

High-definition spatial transcriptomics for in situ tissue profiling.

Spatial and molecular characteristics determine tissue function, yet high-resolution methods to capture both concurrently are lacking. Here, we developed high-definition spatial transcriptomics, which captures RNA from histological tissue sections on a dense, spatially barcoded bead array. Each experiment recovers several hundred thousand transcript-coupled spatial barcodes at 2-µm resolution, as demonstrated in mouse brain and primary breast cancer. This opens the way to high-resolution spatial analysis of cells and tissues.

Conbase: a software for unsupervised discovery of clonal somatic mutations in single cells through read phasing.

Accurate variant calling and genotyping represent major limiting factors for downstream applications of single-cell genomics. Here, we report Conbase for the identification of somatic mutations in single-cell DNA sequencing data. Conbase leverages phased read data from multiple samples in a dataset to achieve increased confidence in somatic variant calls and genotype predictions. Comparing the performance of Conbase to three other methods, we find that Conbase performs best in terms of false discovery rate and specificity and provides superior robustness on simulated data, in vitro expanded fibroblasts and clonal lymphocyte populations isolated directly from a healthy human donor.

ST viewer: a tool for analysis and visualization of spatial transcriptomics datasets.

MOTIVATION: Spatial Transcriptomics (ST) is a technique that combines high-resolution imaging with spatially resolved transcriptome-wide sequencing. This novel type of data opens up many possibilities for analysis and visualization, most of which are either not available with standard tools or too complex for normal users. RESULTS: Here, we present a tool, ST Viewer, which allows real-time interaction, analysis and visualization of Spatial Transcriptomics datasets through a seamless and smooth user interface. AVAILABILITY AND IMPLEMENTATION: The ST Viewer is open source under a MIT license and it is available at https://github.com/SpatialTranscriptomicsResearch/st_viewer. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Barcoded solid-phase RNA capture for Spatial Transcriptomics profiling in mammalian tissue sections.

Spatial resolution of gene expression enables gene expression events to be pinpointed to a specific location in biological tissue. Spatially resolved gene expression in tissue sections is traditionally analyzed using immunohistochemistry (IHC) or in situ hybridization (ISH). These technologies are invaluable tools for pathologists and molecular biologists; however, their throughput is limited to the analysis of only a few genes at a time. Recent advances in RNA sequencing (RNA-seq) have made it possible to obtain unbiased high-throughput gene expression data in bulk. Spatial Transcriptomics combines the benefits of traditional spatially resolved technologies with the massive throughput of RNA-seq. Here, we present a protocol describing how to apply the Spatial Transcriptomics technology to mammalian tissue. This protocol combines histological staining and spatially resolved RNA-seq data from intact tissue sections. Once suitable tissue-specific conditions have been established, library construction and sequencing can be completed in ~5-6 d. Data processing takes a few hours, with the exact timing dependent on the sequencing depth. Our method requires no special instruments and can be performed in any laboratory with access to a cryostat, microscope and next-generation sequencing.

ST Spot Detector: a web-based application for automatic spot and tissue detection for spatial Transcriptomics image datasets.

Motiviation: Spatial Transcriptomics (ST) is a method which combines high resolution tissue imaging with high troughput transcriptome sequencing data. This data must be aligned with the images for correct visualization, a process that involves several manual steps. Results: Here we present ST Spot Detector, a web tool that automates and facilitates this alignment through a user friendly interface. Contact: jose.fernandez.navarro@scilifelab.se. Supplementary information: Supplementary data are available at Bioinformatics online.

Spatial detection of fetal marker genes expressed at low level in adult human heart tissue.

Heart failure is a major health problem linked to poor quality of life and high mortality rates. Hence, novel biomarkers, such as fetal marker genes with low expression levels, could potentially differentiate disease states in order to improve therapy. In many studies on heart failure, cardiac biopsies have been analyzed as uniform pieces of tissue with bulk techniques, but this homogenization approach can mask medically relevant phenotypes occurring only in isolated parts of the tissue. This study examines such spatial variations within and between regions of cardiac biopsies. In contrast to standard RNA sequencing, this approach provides a spatially resolved transcriptome- and tissue-wide perspective of the adult human heart, and enables detection of fetal marker genes expressed by minor subpopulations of cells within the tissue. Analysis of patients with heart failure, with preserved ejection fraction, demonstrated spatially divergent expression of fetal genes in cardiac biopsies.