ST Pipeline: an automated pipeline for spatial mapping of unique transcripts

Spatial transcriptomic revealed intratumor heterogeneity and cancer stem cell enrichment in colorectal cancer metastasis

Metastasis is the main cause of mortality in colorectal cancer (CRC) patients. Exploring the mechanisms of metastasis is of great importance in both clinical and fundamental CRC research. CRC is a highly heterogeneous disease with variable therapeutic outcomes of treatment. In this study, we applied spatial transcriptomics (ST) to generate a tissue-wide transcriptome from two primary colorectal cancer tissues and their matched liver metastatic tissues. Spatial RNA information showed intratumoral heterogeneity (ITH) of both primary and metastatic tissues. The comparison of gene expressions across tissues revealed an apparent enrichment of cancer stem cells (CSCs) in metastatic tissues and identified FOXD1 as a novel metastatic CSC marker. Trajectory and pseudo-time analyses revealed distinct evolutionary trajectories and a dedifferentiation-differentiation process during metastasis. CellphoneDB analysis suggested a dominant interaction of CD74-MIF with tumor cells in metastatic tissues. Further analysis confirmed FOXD1 as a maker of CSCs and the predictor of patient survival, especially in metastatic diseases. Our study found ITH of primary and metastatic tissues and provides novel insights into the cellular mechanisms underlying liver metastasis of CRC and foundations for therapeutic strategies for CRC metastasis.

Spatial host-microbiome sequencing reveals niches in the mouse gut

Mucosal and barrier tissues, such as the gut, lung or skin, are composed of a complex network of cells and microbes forming a tight niche that prevents pathogen colonization and supports host-microbiome symbiosis. Characterizing these networks at high molecular and cellular resolution is crucial for understanding homeostasis and disease. Here we present spatial host-microbiome sequencing (SHM-seq), an all-sequencing-based approach that captures tissue histology, polyadenylated RNAs and bacterial 16S sequences directly from a tissue by modifying spatially barcoded glass surfaces to enable simultaneous capture of host transcripts and hypervariable regions of the 16S bacterial ribosomal RNA. We applied our approach to the mouse gut as a model system, used a deep learning approach for data mapping and detected spatial niches defined by cellular composition and microbial geography. We show that subpopulations of gut cells express specific gene programs in different microenvironments characteristic of regional commensal bacteria and impact host-bacteria interactions. SHM-seq should enhance the study of native host-microbe interactions in health and disease.

Spatial metatranscriptomics resolves host-bacteria-fungi interactomes

The interactions of microorganisms among themselves and with their multicellular host take place at the microscale, forming complex networks and spatial patterns. Existing technology does not allow the simultaneous investigation of spatial interactions between a host and the multitude of its colonizing microorganisms, which limits our understanding of host-microorganism interactions within a plant or animal tissue. Here we present spatial metatranscriptomics (SmT), a sequencing-based approach that leverages 16S/18S/ITS/poly-d(T) multimodal arrays for simultaneous host transcriptome- and microbiome-wide characterization of tissues at 55-µm resolution. We showcase SmT in outdoor-grown Arabidopsis thaliana leaves as a model system, and find tissue-scale bacterial and fungal hotspots. By network analysis, we study inter- and intrakingdom spatial interactions among microorganisms, as well as the host response to microbial hotspots. SmT provides an approach for answering fundamental questions on host-microbiome interplay.

Host-pathogen interactions in the Plasmodium-infected mouse liver at spatial and single-cell resolution

Upon infecting its vertebrate host, the malaria parasite initially invades the liver where it undergoes massive replication, whilst remaining clinically silent. The coordination of host responses across the complex liver tissue during malaria infection remains unexplored. Here, we perform spatial transcriptomics in combination with single-nuclei RNA sequencing over multiple time points to delineate host-pathogen interactions across Plasmodium berghei-infected liver tissues. Our data reveals significant changes in spatial gene expression in the malaria-infected tissues. These include changes related to lipid metabolism in the proximity to sites of Plasmodium infection, distinct inflammation programs between lobular zones, and regions with enrichment of different inflammatory cells, which we term 'inflammatory hotspots'. We also observe significant upregulation of genes involved in inflammation in the control liver tissues of mice injected with mosquito salivary gland components. However, this response is considerably delayed compared to that observed in P. berghei-infected mice. Our study establishes a benchmark for investigating transcriptome changes during host-parasite interactions in tissues, it provides informative insights regarding in vivo study design linked to infection and offers a useful tool for the discovery and validation of de novo intervention strategies aimed at malaria liver stage infection.

Spatial multiomics reveals a subpopulation of fibroblasts associated with cancer stemness in human hepatocellular carcinoma

BACKGROUND

Cancer-associated fibroblasts (CAFs) are the prominent cell type in the tumor microenvironment (TME), and CAF subsets have been identified in various tumors. However, how CAFs spatially coordinate other cell populations within the liver TME to promote cancer progression remains unclear.

METHODS

We combined multi-region proteomics (6 patients, 24 samples), 10X Genomics Visium spatial transcriptomics (11 patients, 25 samples), and multiplexed imaging (92 patients, 264 samples) technologies to decipher the expression heterogeneity, functional diversity, spatial distribution, colocalization, and interaction of fibroblasts. The newly identified CAF subpopulation was validated by cells isolated from 5 liver cancer patients and in vitro functional assays.

RESULTS

We identified a liver CAF subpopulation, marked by the expression of COL1A2, COL4A1, COL4A2, CTGF, and FSTL1, and named F5-CAF. F5-CAF is preferentially located within and around tumor nests and colocalizes with cancer cells with higher stemness in hepatocellular carcinoma (HCC). Multiplexed staining of 92 patients and the bulk transcriptome of 371 patients demonstrated that the abundance of F5-CAFs in HCC was associated with a worse prognosis. Further in vitro experiments showed that F5-CAFs isolated from liver cancer patients can promote the proliferation and stemness of HCC cells.

CONCLUSIONS

We identified a CAF subpopulation F5-CAF in liver cancer, which is associated with cancer stemness and unfavorable prognosis. Our results provide potential mechanisms by which the CAF subset in the TME promotes the development of liver cancer by supporting the survival of cancer stem cells.

A Spatial Transcriptomics Browser for Discovering Gene Expression Landscapes across Microscopic Tissue Sections

A crucial feature of life is its spatial organization and compartmentalization on the molecular, cellular, and tissue levels. Spatial transcriptomics (ST) technology has opened a new chapter of the sequencing revolution, emerging rapidly with transformative effects across biology. This technique produces extensive and complex sequencing data, raising the need for computational methods for their comprehensive analysis and interpretation. We developed the ST browser web tool for the interactive discovery of ST images, focusing on different functional aspects such as single gene expression, the expression of functional gene sets, as well as the inspection of the spatial patterns of cell-cell interactions. As a unique feature, our tool applies self-organizing map (SOM) machine learning to the ST data. Our SOM data portrayal method generates individual gene expression landscapes for each spot in the ST image, enabling its downstream analysis with high resolution. The performance of the spatial browser is demonstrated by disentangling the intra-tumoral heterogeneity of melanoma and the microarchitecture of the mouse brain. The integration of machine-learning-based SOM portrayal into an interactive ST analysis environment opens novel perspectives for the comprehensive knowledge mining of the organization and interactions of cellular ecosystems.

Tissue and cellular spatiotemporal dynamics in colon aging

Tissue structure and molecular circuitry in the colon can be profoundly impacted by systemic age-related effects, but many of the underlying molecular cues remain unclear. Here, we built a cellular and spatial atlas of the colon across three anatomical regions and 11 age groups, encompassing ~1,500 mouse gut tissues profiled by spatial transcriptomics and ~400,000 single nucleus RNA-seq profiles. We developed a new computational framework, cSplotch, which learns a hierarchical Bayesian model of spatially resolved cellular expression associated with age, tissue region, and sex, by leveraging histological features to share information across tissue samples and data modalities. Using this model, we identified cellular and molecular gradients along the adult colonic tract and across the main crypt axis, and multicellular programs associated with aging in the large intestine. Our multi-modal framework for the investigation of cell and tissue organization can aid in the understanding of cellular roles in tissue-level pathology.

Spatially Exploring RNA Biology in Archival Formalin-Fixed Paraffin-Embedded Tissues

Spatial transcriptomics has emerged as a powerful tool for dissecting spatial cellular heterogeneity but as of today is largely limited to gene expression analysis. Yet, the life of RNA molecules is multifaceted and dynamic, requiring spatial profiling of different RNA species throughout the life cycle to delve into the intricate RNA biology in complex tissues. Human disease-relevant tissues are commonly preserved as formalin-fixed and paraffin-embedded (FFPE) blocks, representing an important resource for human tissue specimens. The capability to spatially explore RNA biology in FFPE tissues holds transformative potential for human biology research and clinical histopathology. Here, we present Patho-DBiT combining in situ polyadenylation and deterministic barcoding for spatial full coverage transcriptome sequencing, tailored for probing the diverse landscape of RNA species even in clinically archived FFPE samples. It permits spatial co-profiling of gene expression and RNA processing, unveiling region-specific splicing isoforms, and high-sensitivity transcriptomic mapping of clinical tumor FFPE tissues stored for five years. Furthermore, genome-wide single nucleotide RNA variants can be captured to distinguish different malignant clones from non-malignant cells in human lymphomas. Patho-DBiT also maps microRNA-mRNA regulatory networks and RNA splicing dynamics, decoding their roles in spatial tumorigenesis trajectory. High resolution Patho-DBiT at the cellular level reveals a spatial neighborhood and traces the spatiotemporal kinetics driving tumor progression. Patho-DBiT stands poised as a valuable platform to unravel rich RNA biology in FFPE tissues to study human tissue biology and aid in clinical pathology evaluation.

Choroid plexus mast cells drive tumor-associated hydrocephalus

Tumor-associated hydrocephalus (TAH) is a common and lethal complication of brain metastases. Although other factors beyond mechanical obstructions have been suggested, the exact mechanisms are unknown. Using single-nucleus RNA sequencing and spatial transcriptomics, we find that a distinct population of mast cells locate in the choroid plexus and dramatically increase during TAH. Genetic fate tracing and intracranial mast-cell-specific tryptase knockout showed that choroid plexus mast cells (CPMCs) disrupt cilia of choroid plexus epithelia via the tryptase-PAR2-FoxJ1 pathway and consequently increase cerebrospinal fluid production. Mast cells are also found in the human choroid plexus. Levels of tryptase in cerebrospinal fluid are closely associated with clinical severity of TAH. BMS-262084, an inhibitor of tryptase, can cross the blood-brain barrier, inhibit TAH in vivo, and alleviate mast-cell-induced damage of epithelial cilia in a human pluripotent stem-cell-derived choroid plexus organoid model. Collectively, we uncover the function of CPMCs and provide an attractive therapy for TAH.

Spatial transcriptomic analysis of the mouse brain following chronic social defeat stress

Depression is a highly prevalent and disabling mental disorder, involving numerous genetic changes that are associated with abnormal functions in multiple regions of the brain. However, there is little transcriptomic-wide characterization of chronic social defeat stress (CSDS) to comprehensively compare the transcriptional changes in multiple brain regions. Spatial transcriptomics (ST) was used to reveal the spatial difference of gene expression in the control, resilient (RES) and susceptible (SUS) mouse brains, and annotated eight anatomical brain regions and six cell types. The gene expression profiles uncovered that CSDS leads to gene synchrony changes in different brain regions. Then it was identified that inhibitory neurons and synaptic functions in multiple regions were primarily affected by CSDS. The brain regions Hippocampus (HIP), Isocortex, and Amygdala (AMY) present more pronounced transcriptional changes in genes associated with depressive psychiatric disorders than other regions. Signalling communication between these three brain regions may play a critical role in susceptibility to CSDS. Taken together, this study provides important new insights into CSDS susceptibility at the ST level, which offers a new approach for understanding and treating depression.

High-plex protein and whole transcriptome co-mapping at cellular resolution with spatial CITE-seq

In this study, we extended co-indexing of transcriptomes and epitopes (CITE) to the spatial dimension and demonstrated high-plex protein and whole transcriptome co-mapping. We profiled 189 proteins and whole transcriptome in multiple mouse tissue types with spatial CITE sequencing and then further applied the method to measure 273 proteins and transcriptome in human tissues, revealing spatially distinct germinal center reactions in tonsil and early immune activation in skin at the Coronavirus Disease 2019 mRNA vaccine injection site.

A unified pipeline for FISH spatial transcriptomics

High-throughput spatial transcriptomics has emerged as a powerful tool for investigating the spatial distribution of mRNA expression and its effects on cellular function. There is a lack of standardized tools for analyzing spatial transcriptomics data, leading many groups to write their own in-house tools that are often poorly documented and not generalizable. To address this, we have expanded and improved the starfish library and used those tools to create PIPEFISH, a semi-automated and generalizable pipeline that performs transcript annotation for fluorescence in situ hybridization (FISH)-based spatial transcriptomics. We used this pipeline to annotate transcript locations from three real datasets from three different common types of FISH image-based experiments, MERFISH, seqFISH, and targeted in situ sequencing (ISS), and verified that the results were high quality using the internal quality metrics of the pipeline and also a comparison with an orthogonal method of measuring RNA expression. PIPEFISH is a publicly available and open-source tool.

Spatially resolved transcriptomics: a comprehensive review of their technological advances, applications, and challenges

The ability to explore life kingdoms is largely driven by innovations and breakthroughs in technology, from the invention of the microscope 350 years ago to the recent emergence of single-cell sequencing, by which the scientific community has been able to visualize life at an unprecedented resolution. Most recently, the Spatially Resolved Transcriptomics (SRT) technologies have filled the gap in probing the spatial or even three-dimensional organization of the molecular foundation behind the molecular mysteries of life, including the origin of different cellular populations developed from totipotent cells and human diseases. In this review, we introduce recent progresses and challenges on SRT from the perspectives of technologies and bioinformatic tools, as well as the representative SRT applications. With the currently fast-moving progress of the SRT technologies and promising results from early adopted research projects, we can foresee the bright future of such new tools in understanding life at the most profound analytical level.

Temporal dynamics of microglia-astrocyte interaction in neuroprotective glial scar formation after intracerebral hemorrhage

The role of glial scar after intracerebral hemorrhage (ICH) remains unclear. This study aimed to investigate whether microglia-astrocyte interaction affects glial scar formation and explore the specific function of glial scar. We used a pharmacologic approach to induce microglial depletion during different ICH stages and examine how ablating microglia affects astrocytic scar formation. Spatial transcriptomics (ST) analysis was performed to explore the potential ligand-receptor pair in the modulation of microglia-astrocyte interaction and to verify the functional changes of astrocytic scars at different periods. During the early stage, sustained microglial depletion induced disorganized astrocytic scar, enhanced neutrophil infiltration, and impaired tissue repair. ST analysis indicated that microglia-derived insulin like growth factor 1 (IGF1) modulated astrocytic scar formation via mechanistic target of rapamycin (mTOR) signaling activation. Moreover, repopulating microglia (RM) more strongly activated mTOR signaling, facilitating a more protective scar formation. The combination of IGF1 and osteopontin (OPN) was necessary and sufficient for RM function, rather than IGF1 or OPN alone. At the chronic stage of ICH, the overall net effect of astrocytic scar changed from protective to destructive and delayed microglial depletion could partly reverse this. The vital insight gleaned from our data is that sustained microglial depletion may not be a reasonable treatment strategy for early-stage ICH. Inversely, early-stage IGF1/OPN treatment combined with late-stage PLX3397 treatment is a promising therapeutic strategy. This prompts us to consider the complex temporal dynamics and overall net effect of microglia and astrocytes, and develop elaborate treatment strategies at precise time points after ICH.

Celloscope: a probabilistic model for marker-gene-driven cell type deconvolution in spatial transcriptomics data

Spatial transcriptomics maps gene expression across tissues, posing the challenge of determining the spatial arrangement of different cell types. However, spatial transcriptomics spots contain multiple cells. Therefore, the observed signal comes from mixtures of cells of different types. Here, we propose an innovative probabilistic model, Celloscope, that utilizes established prior knowledge on marker genes for cell type deconvolution from spatial transcriptomics data. Celloscope outperforms other methods on simulated data, successfully indicates known brain structures and spatially distinguishes between inhibitory and excitatory neuron types based in mouse brain tissue, and dissects large heterogeneity of immune infiltrate composition in prostate gland tissue.

Spatial epigenome-transcriptome co-profiling of mammalian tissues

Emerging spatial technologies, including spatial transcriptomics and spatial epigenomics, are becoming powerful tools for profiling of cellular states in the tissue context1-5. However, current methods capture only one layer of omics information at a time, precluding the possibility of examining the mechanistic relationship across the central dogma of molecular biology. Here, we present two technologies for spatially resolved, genome-wide, joint profiling of the epigenome and transcriptome by cosequencing chromatin accessibility and gene expression, or histone modifications (H3K27me3, H3K27ac or H3K4me3) and gene expression on the same tissue section at near-single-cell resolution. These were applied to embryonic and juvenile mouse brain, as well as adult human brain, to map how epigenetic mechanisms control transcriptional phenotype and cell dynamics in tissue. Although highly concordant tissue features were identified by either spatial epigenome or spatial transcriptome we also observed distinct patterns, suggesting their differential roles in defining cell states. Linking epigenome to transcriptome pixel by pixel allows the uncovering of new insights in spatial epigenetic priming, differentiation and gene regulation within the tissue architecture. These technologies are of great interest in life science and biomedical research.

Computational solutions for spatial transcriptomics

Transcriptome level expression data connected to the spatial organization of the cells and molecules would allow a comprehensive understanding of how gene expression is connected to the structure and function in the biological systems. The spatial transcriptomics platforms may soon provide such information. However, the current platforms still lack spatial resolution, capture only a fraction of the transcriptome heterogeneity, or lack the throughput for large scale studies. The strengths and weaknesses in current ST platforms and computational solutions need to be taken into account when planning spatial transcriptomics studies. The basis of the computational ST analysis is the solutions developed for single-cell RNA-sequencing data, with advancements taking into account the spatial connectedness of the transcriptomes. The scRNA-seq tools are modified for spatial transcriptomics or new solutions like deep learning-based joint analysis of expression, spatial, and image data are developed to extract biological information in the spatially resolved transcriptomes. The computational ST analysis can reveal remarkable biological insights into spatial patterns of gene expression, cell signaling, and cell type variations in connection with cell type-specific signaling and organization in complex tissues. This review covers the topics that help choosing the platform and computational solutions for spatial transcriptomics research. We focus on the currently available ST methods and platforms and their strengths and limitations. Of the computational solutions, we provide an overview of the analysis steps and tools used in the ST data analysis. The compatibility with the data types and the tools provided by the current ST analysis frameworks are summarized.

Spatially resolved clonal copy number alterations in benign and malignant tissue

Defining the transition from benign to malignant tissue is fundamental to improving early diagnosis of cancer1. Here we use a systematic approach to study spatial genome integrity in situ and describe previously unidentified clonal relationships. We used spatially resolved transcriptomics2 to infer spatial copy number variations in >120,000 regions across multiple organs, in benign and malignant tissues. We demonstrate that genome-wide copy number variation reveals distinct clonal patterns within tumours and in nearby benign tissue using an organ-wide approach focused on the prostate. Our results suggest a model for how genomic instability arises in histologically benign tissue that may represent early events in cancer evolution. We highlight the power of capturing the molecular and spatial continuums in a tissue context and challenge the rationale for treatment paradigms, including focal therapy.

Spacemake: processing and analysis of large-scale spatial transcriptomics data

BACKGROUND

Spatial sequencing methods increasingly gain popularity within RNA biology studies. State-of-the-art techniques quantify messenger RNA expression levels from tissue sections and at the same time register information about the original locations of the molecules in the tissue. The resulting data sets are processed and analyzed by accompanying software that, however, is incompatible across inputs from different technologies.

FINDINGS

Here, we present spacemake, a modular, robust, and scalable spatial transcriptomics pipeline built in Snakemake and Python. Spacemake is designed to handle all major spatial transcriptomics data sets and can be readily configured for other technologies. It can process and analyze several samples in parallel, even if they stem from different experimental methods. Spacemake's unified framework enables reproducible data processing from raw sequencing data to automatically generated downstream analysis reports. Spacemake is built with a modular design and offers additional functionality such as sample merging, saturation analysis, and analysis of long reads as separate modules. Moreover, spacemake employs novoSpaRc to integrate spatial and single-cell transcriptomics data, resulting in increased gene counts for the spatial data set. Spacemake is open source and extendable, and it can be seamlessly integrated with existing computational workflows.

Cell-cell communication inference and analysis in the tumour microenvironments from single-cell transcriptomics: data resources and computational strategies

Carcinomas are complex ecosystems composed of cancer, stromal and immune cells. Communication between these cells and their microenvironments induces cancer progression and causes therapy resistance. In order to improve the treatment of cancers, it is essential to quantify crosstalk between and within various cell types in a tumour microenvironment. Focusing on the coordinated expression patterns of ligands and cognate receptors, cell-cell communication can be inferred through ligand-receptor interactions (LRIs). In this manuscript, we carry out the following work: (i) introduce pipeline for ligand-receptor-mediated intercellular communication estimation from single-cell transcriptomics and list a few available LRI-related databases and visualization tools; (ii) demonstrate seven classical intercellular communication scoring strategies, highlight four types of representative intercellular communication inference methods, including network-based approaches, machine learning-based approaches, spatial information-based approaches and other approaches; (iii) summarize the evaluation and validation avenues for intercellular communication inference and analyze the advantages and limitations for the above four types of cell-cell communication methods; (iv) comment several major challenges while provide further research directions for intercellular communication analysis in the tumour microenvironments. We anticipate that this work helps to better understand intercellular crosstalk and to further develop powerful cell-cell communication estimation tools for tumor-targeted therapy.

spatialGE: quantification and visualization of the tumor microenvironment heterogeneity using spatial transcriptomics

SUMMARY

Spatially resolved transcriptomics promises to increase our understanding of the tumor microenvironment and improve cancer prognosis and therapies. Nonetheless, analytical methods to explore associations between the spatial heterogeneity of the tumor and clinical data are not available. Hence, we have developed spatialGE, a software that provides visualizations and quantification of the tumor microenvironment heterogeneity through gene expression surfaces, spatial heterogeneity statistics that can be compared against clinical information, spot-level cell deconvolution and spatially informed clustering, all using a new data object to store data and resulting analyses simultaneously.

AVAILABILITY AND IMPLEMENTATION

The R package and tutorial/vignette are available at https://github.com/FridleyLab/spatialGE. A script to reproduce the analyses in this manuscript is available in Supplementary information. The Thrane study data included in spatialGE was made available from the public available from the website https://www.spatialresearch.org/resources-published-datasets/doi-10-1158-0008-5472-can-18-0747/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

SM-Omics is an automated platform for high-throughput spatial multi-omics

The spatial organization of cells and molecules plays a key role in tissue function in homeostasis and disease. Spatial transcriptomics has recently emerged as a key technique to capture and positionally barcode RNAs directly in tissues. Here, we advance the application of spatial transcriptomics at scale, by presenting Spatial Multi-Omics (SM-Omics) as a fully automated, high-throughput all-sequencing based platform for combined and spatially resolved transcriptomics and antibody-based protein measurements. SM-Omics uses DNA-barcoded antibodies, immunofluorescence or a combination thereof, to scale and combine spatial transcriptomics and spatial antibody-based multiplex protein detection. SM-Omics allows processing of up to 64 in situ spatial reactions or up to 96 sequencing-ready libraries, of high complexity, in a ~2 days process. We demonstrate SM-Omics in the mouse brain, spleen and colorectal cancer model, showing its broad utility as a high-throughput platform for spatial multi-omics.

Organizing your space: The potential for integrating spatial transcriptomics and 3D imaging data in plants

Plant cells communicate information for the regulation of development and responses to external stresses. A key form of this communication is transcriptional regulation, accomplished via complex gene networks operating both locally and systemically. To fully understand how genes are regulated across plant tissues and organs, high resolution, multi-dimensional spatial transcriptional data must be acquired and placed within a cellular and organismal context. Spatial transcriptomics (ST) typically provides a two-dimensional spatial analysis of gene expression of tissue sections that can be stacked to render three-dimensional data. For example, X-ray and light-sheet microscopy provide sub-micron scale volumetric imaging of cellular morphology of tissues, organs, or potentially entire organisms. Linking these technologies could substantially advance transcriptomics in plant biology and other fields. Here, we review advances in ST and 3D microscopy approaches and describe how these technologies could be combined to provide high resolution, spatially organized plant tissue transcript mapping.

Spatially Resolved Expression of Transposable Elements in Disease and Somatic Tissue with SpatialTE

Spatial transcriptomics (ST) is transforming the way we can study gene expression and its regulation through position-specific resolution within tissues. However, as in bulk RNA-Seq, transposable elements (TEs) are not being studied due to their highly repetitive nature. In recent years, TEs have been recognized as important regulators of gene expression, and thus, TE expression analysis in a spatially resolved manner could further help to understand their role in gene regulation within tissues. We present SpatialTE, a tool to analyze TE expression from ST datasets and show its application in somatic and diseased tissues. The results indicate that TEs have spatially regulated expression patterns and that their expression profiles are spatially altered in ALS disease, indicating that TEs might perform differential regulatory functions within tissue organs. We have made SpatialTE publicly available as open-source software under an MIT license.

Spatial Transcriptomics to define transcriptional patterns of zonation and structural components in the mouse liver

Reconstruction of heterogeneity through single cell transcriptional profiling has greatly advanced our understanding of the spatial liver transcriptome in recent years. However, global transcriptional differences across lobular units remain elusive in physical space. Here, we apply Spatial Transcriptomics to perform transcriptomic analysis across sectioned liver tissue. We confirm that the heterogeneity in this complex tissue is predominantly determined by lobular zonation. By introducing novel computational approaches, we enable transcriptional gradient measurements between tissue structures, including several lobules in a variety of orientations. Further, our data suggests the presence of previously transcriptionally uncharacterized structures within liver tissue, contributing to the overall spatial heterogeneity of the organ. This study demonstrates how comprehensive spatial transcriptomic technologies can be used to delineate extensive spatial gene expression patterns in the liver, indicating its future impact for studies of liver function, development and regeneration as well as its potential in pre-clinical and clinical pathology.

Music of metagenomics-a review of its applications, analysis pipeline, and associated tools

This humble effort highlights the intricate details of metagenomics in a simple, poetic, and rhythmic way. The paper enforces the significance of the research area, provides details about major analytical methods, examines the taxonomy and assembly of genomes, emphasizes some tools, and concludes by celebrating the richness of the ecosystem populated by the "metagenome."

Space: the final frontier — achieving single-cell, spatially resolved transcriptomics in plants

Single-cell RNA-seq is a tool that generates a high resolution of transcriptional data that can be used to understand regulatory networks in biological systems. In plants, several methods have been established for transcriptional analysis in tissue sections, cell types, and/or single cells. These methods typically require cell sorting, transgenic plants, protoplasting, or other damaging or laborious processes. Additionally, the majority of these technologies lose most or all spatial resolution during implementation. Those that offer a high spatial resolution for RNA lack breadth in the number of transcripts characterized. Here, we briefly review the evolution of spatial transcriptomics methods and we highlight recent advances and current challenges in sequencing, imaging, and computational aspects toward achieving 3D spatial transcriptomics of plant tissues with a resolution approaching single cells. We also provide a perspective on the potential opportunities to advance this novel methodology in plants.

Deciphering cell-cell interactions and communication from gene expression

Cell-cell interactions orchestrate organismal development, homeostasis and single-cell functions. When cells do not properly interact or improperly decode molecular messages, disease ensues. Thus, the identification and quantification of intercellular signalling pathways has become a common analysis performed across diverse disciplines. The expansion of protein-protein interaction databases and recent advances in RNA sequencing technologies have enabled routine analyses of intercellular signalling from gene expression measurements of bulk and single-cell data sets. In particular, ligand-receptor pairs can be used to infer intercellular communication from the coordinated expression of their cognate genes. In this Review, we highlight discoveries enabled by analyses of cell-cell interactions from transcriptomic data and review the methods and tools used in this context.

The spatial RNA integrity number assay for in situ evaluation of transcriptome quality

The RNA integrity number (RIN) is a frequently used quality metric to assess the completeness of rRNA, as a proxy for the corresponding mRNA in a tissue. Current methods operate at bulk resolution and provide a single average estimate for the whole sample. Spatial transcriptomics technologies have emerged and shown their value by placing gene expression into a tissue context, resulting in transcriptional information from all tissue regions. Thus, the ability to estimate RNA quality in situ has become of utmost importance to overcome the limitation with a bulk rRNA measurement. Here we show a new tool, the spatial RNA integrity number (sRIN) assay, to assess the rRNA completeness in a tissue wide manner at cellular resolution. We demonstrate the use of sRIN to identify spatial variation in tissue quality prior to more comprehensive spatial transcriptomics workflows.

Applying Machine Learning for Integration of Multi-Modal Genomics Data and Imaging Data to Quantify Heterogeneity in Tumour Tissues

With rapid advances in experimental instruments and protocols, imaging and sequencing data are being generated at an unprecedented rate contributing significantly to the current and coming big biomedical data. Meanwhile, unprecedented advances in computational infrastructure and analysis algorithms are realizing image-based digital diagnosis not only in radiology and cardiology but also oncology and other diseases. Machine learning methods, especially deep learning techniques, are already and broadly implemented in diverse technological and industrial sectors, but their applications in healthcare are just starting. Uniquely in biomedical research, a vast potential exists to integrate genomics data with histopathological imaging data. The integration has the potential to extend the pathologist's limits and boundaries, which may create breakthroughs in diagnosis, treatment, and monitoring at molecular and tissue levels. Moreover, the applications of genomics data are realizing the potential for personalized medicine, making diagnosis, treatment, monitoring, and prognosis more accurate. In this chapter, we discuss machine learning methods readily available for digital pathology applications, new prospects of integrating spatial genomics data on tissues with tissue morphology, and frontier approaches to combining genomics data with pathological imaging data. We present perspectives on how artificial intelligence can be synergized with molecular genomics and imaging to make breakthroughs in biomedical and translational research for computer-aided applications.

TGFβ-blockade uncovers stromal plasticity in tumors by revealing the existence of a subset of interferon-licensed fibroblasts

Despite the increasing interest in targeting stromal elements of the tumor microenvironment, we still face tremendous challenges in developing adequate therapeutics to modify the tumor stromal landscape. A major obstacle to this is our poor understanding of the phenotypic and functional heterogeneity of stromal cells in tumors. Herein, we perform an unbiased interrogation of tumor mesenchymal cells, delineating the co-existence of distinct subsets of cancer-associated fibroblasts (CAFs) in the microenvironment of murine carcinomas, each endowed with unique phenotypic features and functions. Furthermore, our study shows that neutralization of TGFβ in vivo leads to remodeling of CAF dynamics, greatly reducing the frequency and activity of the myofibroblast subset, while promoting the formation of a fibroblast population characterized by strong response to interferon and heightened immunomodulatory properties. These changes correlate with the development of productive anti-tumor immunity and greater efficacy of PD1 immunotherapy. Along with providing the scientific rationale for the evaluation of TGFβ and PD1 co-blockade in the clinical setting, this study also supports the concept of plasticity of the stromal cell landscape in tumors, laying the foundation for future investigations aimed at defining pathways and molecules to program CAF composition for cancer therapy.

High-Spatial-Resolution Multi-Omics Sequencing via Deterministic Barcoding in Tissue

We present deterministic barcoding in tissue for spatial omics sequencing (DBiT-seq) for co-mapping of mRNAs and proteins in a formaldehyde-fixed tissue slide via next-generation sequencing (NGS). Parallel microfluidic channels were used to deliver DNA barcodes to the surface of a tissue slide, and crossflow of two sets of barcodes, A1-50 and B1-50, followed by ligation in situ, yielded a 2D mosaic of tissue pixels, each containing a unique full barcode AB. Application to mouse embryos revealed major tissue types in early organogenesis as well as fine features like microvasculature in a brain and pigmented epithelium in an eye field. Gene expression profiles in 10-μm pixels conformed into the clusters of single-cell transcriptomes, allowing for rapid identification of cell types and spatial distributions. DBiT-seq can be adopted by researchers with no experience in microfluidics and may find applications in a range of fields including developmental biology, cancer biology, neuroscience, and clinical pathology.

SpaCell: integrating tissue morphology and spatial gene expression to predict disease cells

MOTIVATION

Spatial transcriptomics (ST) technology is increasingly being applied because it enables the measurement of spatial gene expression in an intact tissue along with imaging morphology of the same tissue. However, current analysis methods for ST data do not use image pixel information, thus missing the quantitative links between gene expression and tissue morphology.

RESULTS

We developed a user-friendly deep learning software, SpaCell, to integrate millions of pixel intensity values with thousands of gene expression measurements from spatially barcoded spots in a tissue. We show the integration approach outperforms the use of gene-count data alone or imaging data alone to build deep learning models to identify cell types or predict labels of tissue images with high resolution and accuracy.

AVAILABILITY AND IMPLEMENTATION

The SpaCell package is open source under an MIT licence and it is available at https://github.com/BiomedicalMachineLearning/SpaCell.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transcriptional Spatial Profiling of Cancer Tissues in the Era of Immunotherapy: The Potential and Promise

Intratumoral heterogeneity poses a major challenge to making an accurate diagnosis and establishing personalized treatment strategies for cancer patients. Moreover, this heterogeneity might underlie treatment resistance, disease progression, and cancer relapse. For example, while immunotherapies can confer a high success rate, selective pressures coupled with dynamic evolution within a tumour can drive the emergence of drug-resistant clones that allow tumours to persist in certain patients. To improve immunotherapy efficacy, researchers have used transcriptional spatial profiling techniques to identify and subsequently block the source of tumour heterogeneity. In this review, we describe and assess the different technologies available for such profiling within a cancer tissue. We first outline two well-known approaches, in situ hybridization and digital spatial profiling. Then, we highlight the features of an emerging technology known as Visium Spatial Gene Expression Solution. Visium generates quantitative gene expression data and maps them to the tissue architecture. By retaining spatial information, we are well positioned to identify novel biomarkers and perform computational analyses that might inform on novel combinatorial immunotherapies.

Spatial Transcriptomics and In Situ Sequencing to Study Alzheimer's Disease

Although complex inflammatory-like alterations are observed around the amyloid plaques of Alzheimer's disease (AD), little is known about the molecular changes and cellular interactions that characterize this response. We investigate here, in an AD mouse model, the transcriptional changes occurring in tissue domains in a 100-μm diameter around amyloid plaques using spatial transcriptomics. We demonstrate early alterations in a gene co-expression network enriched for myelin and oligodendrocyte genes (OLIGs), whereas a multicellular gene co-expression network of plaque-induced genes (PIGs) involving the complement system, oxidative stress, lysosomes, and inflammation is prominent in the later phase of the disease. We confirm the majority of the observed alterations at the cellular level using in situ sequencing on mouse and human brain sections. Genome-wide spatial transcriptomics analysis provides an unprecedented approach to untangle the dysregulated cellular network in the vicinity of pathogenic hallmarks of AD and other brain diseases.

Multimodal Analysis of Composition and Spatial Architecture in Human Squamous Cell Carcinoma

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

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.

Molecular atlas of the adult mouse brain

Brain maps are essential for integrating information and interpreting the structure-function relationship of circuits and behavior. We aimed to generate a systematic classification of the adult mouse brain based purely on the unbiased identification of spatially defining features by employing whole-brain spatial transcriptomics. We found that the molecular information was sufficient to deduce the complex and detailed neuroanatomical organization of the brain. The unsupervised (non-expert, data-driven) classification revealed new area- and layer-specific subregions, for example in isocortex and hippocampus, and new subdivisions of striatum. The molecular atlas further supports the characterization of the spatial identity of neurons from their single-cell RNA profile, and provides a resource for annotating the brain using a minimal gene set-a brain palette. In summary, we have established a molecular atlas to formally define the spatial organization of brain regions, including the molecular code for mapping and targeting of discrete neuroanatomical domains.

Automation of Spatial Transcriptomics library preparation to enable rapid and robust insights into spatial organization of tissues

BACKGROUND

Interest in studying the spatial distribution of gene expression in tissues is rapidly increasing. Spatial Transcriptomics is a novel sequencing-based technology that generates high-throughput information on the distribution, heterogeneity and co-expression of cells in tissues. Unfortunately, manual preparation of high-quality sequencing libraries is time-consuming and subject to technical variability due to human error during manual pipetting, which results in sample swapping and the accidental introduction of batch effects. All these factors complicate the production and interpretation of biological datasets.

RESULTS

We have integrated an Agilent Bravo Automated Liquid Handling Platform into the Spatial Transcriptomics workflow. Compared to the previously reported Magnatrix 8000+ automated protocol, this approach increases the number of samples processed per run, reduces sample preparation time by 35%, and minimizes batch effects between samples. The new approach is also shown to be highly accurate and almost completely free from technical variability between prepared samples.

CONCLUSIONS

The new automated Spatial Transcriptomics protocol using the Agilent Bravo Automated Liquid Handling Platform rapidly generates high-quality Spatial Transcriptomics libraries. Given the wide use of the Agilent Bravo Automated Liquid Handling Platform in research laboratories and facilities, this will allow many researchers to quickly create robust Spatial Transcriptomics libraries.

Spatial transcriptomics identifies spatially dysregulated expression of GRM3 and USP47 in amyotrophic lateral sclerosis

AIMS

The mechanisms underlying the selective degeneration of motor neurones in amyotrophic lateral sclerosis (ALS) are poorly understood. The aim of this study was to implement spatially resolved RNA sequencing in human post mortem cortical tissue from an ALS patient harbouring the C9orf72 hexanucleotide repeat expansion to identify dysregulated transcripts that may account for differential vulnerabilities of distinct (i) cell types and (ii) brain regions in the pathogenesis of ALS.

METHODS

Using spatial transcriptomics (ST) we analysed the transcriptome of post mortem brain tissue, with spatial resolution down to 100 μm. Validation of these findings was then performed using BaseScope, an adapted, in situ hybridization technique with single-transcript single-cell-resolution, providing extensive regional and cell-type specific confirmation of these dysregulated transcripts. The validation cohort was then extended to include multiple post mortem brain regions and spinal cord tissue from an extended cohort of C9orf72, sporadic ALS (sALS) and SOD1 ALS cases.

RESULTS

We identified sixteen dysregulated transcripts of proteins that have roles within six disease-related pathways. Furthermore, these complementary molecular pathology techniques converged to identify two spatially dysregulated transcripts, GRM3 and USP47, that are commonly dysregulated across sALS, SOD1 and C9orf72 cases alike.

CONCLUSIONS

This study presents the first description of ST in human post mortem cortical tissue from an ALS patient harbouring the C9orf72 hexanucleotide repeat expansion. These data taken together highlight the importance of preserving spatial resolution, facilitating the identification of genes whose dysregulation may in part underlie regional susceptibilities to ALS, crucially highlighting potential therapeutic and diagnostic targets.

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.

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.

Next-generation computational tools for interrogating cancer immunity

The remarkable success of cancer therapies with immune checkpoint blockers is revolutionizing oncology and has sparked intensive basic and translational research into the mechanisms of cancer-immune cell interactions. In parallel, numerous novel cutting-edge technologies for comprehensive molecular and cellular characterization of cancer immunity have been developed, including single-cell sequencing, mass cytometry and multiplexed spatial cellular phenotyping. In order to process, analyse and visualize multidimensional data sets generated by these technologies, computational methods and software tools are required. Here, we review computational tools for interrogating cancer immunity, discuss advantages and limitations of the various methods and provide guidelines to assist in method selection.

Spatiotemporal dynamics of molecular pathology in amyotrophic lateral sclerosis

Paralysis occurring in amyotrophic lateral sclerosis (ALS) results from denervation of skeletal muscle as a consequence of motor neuron degeneration. Interactions between motor neurons and glia contribute to motor neuron loss, but the spatiotemporal ordering of molecular events that drive these processes in intact spinal tissue remains poorly understood. Here, we use spatial transcriptomics to obtain gene expression measurements of mouse spinal cords over the course of disease, as well as of postmortem tissue from ALS patients, to characterize the underlying molecular mechanisms in ALS. We identify pathway dynamics, distinguish regional differences between microglia and astrocyte populations at early time points, and discern perturbations in several transcriptional pathways shared between murine models of ALS and human postmortem spinal cords.

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.

Spatially Resolved Transcriptomics Enables Dissection of Genetic Heterogeneity in Stage III Cutaneous Malignant Melanoma

Cutaneous malignant melanoma (melanoma) is characterized by a high mutational load, extensive intertumoral and intratumoral genetic heterogeneity, and complex tumor microenvironment (TME) interactions. Further insights into the mechanisms underlying melanoma are crucial for understanding tumor progression and responses to treatment. Here we adapted the technology of spatial transcriptomics (ST) to melanoma lymph node biopsies and successfully sequenced the transcriptomes of over 2,200 tissue domains. Deconvolution combined with traditional approaches for dimensional reduction of transcriptome-wide data enabled us to both visualize the transcriptional landscape within the tissue and identify gene expression profiles linked to specific histologic entities. Our unsupervised analysis revealed a complex spatial intratumoral composition of melanoma metastases that was not evident through morphologic annotation. Each biopsy showed distinct gene expression profiles and included examples of the coexistence of multiple melanoma signatures within a single tumor region as well as shared profiles for lymphoid tissue characterized according to their spatial location and gene expression profiles. The lymphoid area in close proximity to the tumor region displayed a specific expression pattern, which may reflect the TME, a key component to fully understanding tumor progression. In conclusion, using the ST technology to generate gene expression profiles reveals a detailed landscape of melanoma metastases. This should inspire researchers to integrate spatial information into analyses aiming to identify the factors underlying tumor progression and therapy outcome.Significance: Applying ST technology to gene expression profiling in melanoma lymph node metastases reveals a complex transcriptional landscape in a spatial context, which is essential for understanding the multiple components of tumor progression and therapy outcome. Cancer Res; 78(20); 5970-9. ©2018 AACR.