Spatial molecular profiling: platforms, applications and analysis tools

Gene count normalization in single-cell imaging-based spatially resolved transcriptomics

Recent advances in imaging-based spatially resolved transcriptomics (im-SRT) technologies now enable high-throughput profiling of targeted genes and their locations in fixed tissues. Normalization of gene expression data is often needed to account for technical factors that may confound underlying biological signals. Here, we investigate the potential impact of different gene count normalization methods with different targeted gene panels in the analysis and interpretation of im-SRT data. Using different simulated gene panels that overrepresent genes expressed in specific tissue regions or cell types, we demonstrate how normalization methods based on detected gene counts per cell differentially impact normalized gene expression magnitudes in a region- or cell type-specific manner. We show that these normalization-induced effects may reduce the reliability of downstream analyses including differential gene expression, gene fold change, and spatially variable gene analysis, introducing false positive and false negative results when compared to results obtained from gene panels that are more representative of the gene expression of the tissue's component cell types. These effects are not observed with normalization approaches that do not use detected gene counts for gene expression magnitude adjustment, such as with cell volume or cell area normalization. We recommend using non-gene count-based normalization approaches when feasible and evaluating gene panel representativeness before using gene count-based normalization methods if necessary. Overall, we caution that the choice of normalization method and gene panel may impact the biological interpretation of the im-SRT data.

Multiplex protein imaging in tumour biology

Tissue imaging has become much more colourful in the past decade. Advances in both experimental and analytical methods now make it possible to image protein markers in tissue samples in high multiplex. The ability to routinely image 40-50 markers simultaneously, at single-cell or subcellular resolution, has opened up new vistas in the study of tumour biology. Cellular phenotypes, interaction, communication and spatial organization have become amenable to molecular-level analysis, and application to patient cohorts has identified clinically relevant cellular and tissue features in several cancer types. Here, we review the use of multiplex protein imaging methods to study tumour biology, discuss ongoing attempts to combine these approaches with other forms of spatial omics, and highlight challenges in the field.

stAA: adversarial graph autoencoder for spatial clustering task of spatially resolved transcriptomics

With the development of spatially resolved transcriptomics technologies, it is now possible to explore the gene expression profiles of single cells while preserving their spatial context. Spatial clustering plays a key role in spatial transcriptome data analysis. In the past 2 years, several graph neural network-based methods have emerged, which significantly improved the accuracy of spatial clustering. However, accurately identifying the boundaries of spatial domains remains a challenging task. In this article, we propose stAA, an adversarial variational graph autoencoder, to identify spatial domain. stAA generates cell embedding by leveraging gene expression and spatial information using graph neural networks and enforces the distribution of cell embeddings to a prior distribution through Wasserstein distance. The adversarial training process can make cell embeddings better capture spatial domain information and more robust. Moreover, stAA incorporates global graph information into cell embeddings using labels generated by pre-clustering. Our experimental results show that stAA outperforms the state-of-the-art methods and achieves better clustering results across different profiling platforms and various resolutions. We also conducted numerous biological analyses and found that stAA can identify fine-grained structures in tissues, recognize different functional subtypes within tumors and accurately identify developmental trajectories.

Best practices for single-cell analysis across modalities

Recent advances in single-cell technologies have enabled high-throughput molecular profiling of cells across modalities and locations. Single-cell transcriptomics data can now be complemented by chromatin accessibility, surface protein expression, adaptive immune receptor repertoire profiling and spatial information. The increasing availability of single-cell data across modalities has motivated the development of novel computational methods to help analysts derive biological insights. As the field grows, it becomes increasingly difficult to navigate the vast landscape of tools and analysis steps. Here, we summarize independent benchmarking studies of unimodal and multimodal single-cell analysis across modalities to suggest comprehensive best-practice workflows for the most common analysis steps. Where independent benchmarks are not available, we review and contrast popular methods. Our article serves as an entry point for novices in the field of single-cell (multi-)omic analysis and guides advanced users to the most recent best practices.

ScopeViewer: A Browser-Based Solution for Visualizing Spatial Transcriptomics Data

Motivation

Spatial transcriptomics (ST) enables a high-resolution interrogation of molecular characteristics within specific spatial contexts and tissue morphology. Despite its potential, visualization of ST data is a challenging task due to the complexities in handling, sharing and visualizing large image datasets together with molecular information.

Results

We introduce ScopeViewer, a browser-based software designed to overcome these challenges. ScopeViewer offers the following functionalities: (1) It visualizes large image data and associated annotations at various zoom levels, allowing for intricate exploration of the data; (2) It enables dual interactive viewing of the original images along with their annotations, providing a comprehensive understanding of the context; (3) It displays spatial molecular features with optimized bandwidth, ensuring a smooth user experience; and (4) It bolsters data security by circumventing data transfers.

Availability

ScopeViewer is available at: https://datacommons.swmed.edu/scopeviewer.

Spatial Transcriptomics Arena (STAr): an Integrated Platform for Spatial Transcriptomics Methodology Research

The emerging field of spatially resolved transcriptomics (SRT) has revolutionized biomedical research. SRT quantifies expression levels at different spatial locations, providing a new and powerful tool to interrogate novel biological insights. An essential question in the analysis of SRT data is to identify spatially variable (SV) genes; the expression levels of such genes have spatial variation across different tissues. SV genes usually play an important role in underlying biological mechanisms and tissue heterogeneity. Currently, several computational methods have been developed to detect such genes; however, there is a lack of unbiased assessment of these approaches to guide researchers in selecting the appropriate methods for their specific biomedical applications. In addition, it is difficult for researchers to implement different existing methods for either biological study or methodology development. Furthermore, currently available public SRT datasets are scattered across different websites and preprocessed in different ways, posing additional obstacles for quantitative researchers developing computational methods for SRT data analysis. To address these challenges, we designed Spatial Transcriptomics Arena (STAr), an open platform comprising 193 curated datasets from seven technologies, seven statistical methods, and analysis results. This resource allows users to retrieve high-quality datasets, apply or develop spatial gene detection methods, as well as browse and compare spatial gene analysis results. It also enables researchers to comprehensively evaluate SRT methodology research in both simulated and real datasets. Altogether, STAr is an integrated research resource intended to promote reproducible research and accelerate rigorous methodology development, which can eventually lead to an improved understanding of biological processes and diseases. STAr can be accessed at https://lce.biohpc.swmed.edu/star/ .

A molecular classification of gastric cancer associated with distinct clinical outcomes and validated by an XGBoost-based prediction model

Gastric cancer (GC) is a heterogeneous disease and a leading cause of cancer-related deaths. Discovering robust, clinically relevant molecular classifications is critical for guiding personalized therapies for GC. Here, we propose a refined molecular classification scheme for GC using integrated optimal algorithms and multi-omics data. Based on the important features of mRNA, microRNA, and DNA methylation data selected by the multivariate Cox regression model, three subtypes linked to distinct clinical outcomes were identified by combining similarity network fusion and consensus clustering methods. Three subtypes were validated by an extreme gradient boosting machine learning prediction model with 125 differentially expressed genes in multiple independent cohorts. The molecular characteristics of mutation signatures, characteristic gene sets, driver genes, and chemotherapy sensitivity for each subtype were also identified: subtype 1 was associated with favorable prognosis and characterized by high ARID1A and PIK3CA mutations, subtype 2 was associated with a poor prognosis and harbored high recurrent TP53 mutations, and subtype 3 was associated with high CHD1, APOA1 mutations, and a poor prognosis. The proposed three-subtype scheme achieved a better clinical prediction performance (area under the curve value = 0.71) than The Cancer Genome Atlas classification, which may provide a practical subtyping framework to improve the treatment of GC.

Application of Machine Learning in Spatial Proteomics

Spatial proteomics is an interdisciplinary field that investigates the localization and dynamics of proteins, and it has gained extensive attention in recent years, especially the subcellular proteomics. Numerous evidence indicate that the subcellular localization of proteins is associated with various cellular processes and disease progression. Mass spectrometry (MS)-based and imaging-based experimental approaches have been developed to acquire large-scale spatial proteomic data. To allow the reliable analysis of increasingly complex spatial proteomics data, machine learning (ML) methods have been widely used in both MS-based and imaging-based spatial proteomic data analysis pipelines. Here, we comprehensively survey the applications of ML in spatial proteomics from following aspects: (1) data resources for spatial proteome are comprehensively introduced; (2) the roles of different ML algorithms in data analysis pipelines are elaborated; (3) successful applications of spatial proteomics and several analytical tools integrating ML methods are presented; (4) challenges existing in modern ML-based spatial proteomics studies are discussed. This review provides guidelines for researchers seeking to apply ML methods to analyze spatial proteomic data and can facilitate insightful understanding of cell biology as well as the future research in medical and drug discovery communities.

Cancer Genomics

In the past decade, genomics has fundamentally changed our view of cancer biology, allowing comprehensive analyses of mutations, copy number alterations, structural variants, gene expression and DNA methylation profiles in large-scale studies across different cancer types. Efforts like The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) have fostered international collaborations for cancer genomic analyses and have generated public databases that give scientists around the world access to thoroughly curated data, which have been extensively used as a tool for further hypothesis driven research on several aspects of cancer biology. In parallel, some of these findings are being translated into specific clinical benefits for cancer patients. In this review, we provide a brief historical description of the evolution of international public cancer genome projects and related databases, as well as we discuss about their impact on general cancer research.

Knowledge-graph-based cell-cell communication inference for spatially resolved transcriptomic data with SpaTalk

Spatially resolved transcriptomics provides genetic information in space toward elucidation of the spatial architecture in intact organs and the spatially resolved cell-cell communications mediating tissue homeostasis, development, and disease. To facilitate inference of spatially resolved cell-cell communications, we here present SpaTalk, which relies on a graph network and knowledge graph to model and score the ligand-receptor-target signaling network between spatially proximal cells by dissecting cell-type composition through a non-negative linear model and spatial mapping between single-cell transcriptomic and spatially resolved transcriptomic data. The benchmarked performance of SpaTalk on public single-cell spatial transcriptomic datasets is superior to that of existing inference methods. Then we apply SpaTalk to STARmap, Slide-seq, and 10X Visium data, revealing the in-depth communicative mechanisms underlying normal and disease tissues with spatial structure. SpaTalk can uncover spatially resolved cell-cell communications for single-cell and spot-based spatially resolved transcriptomic data universally, providing valuable insights into spatial inter-cellular tissue dynamics.

Cell-cell communication is a vital feature involving numerous biological processes. Here, the authors develop SpaTalk, a cell-cell communication inference method using knowledge graph for spatially resolved transcriptomic data, providing valuable insights into spatial intercellular tissue dynamics.

Discovering Immune-Mediated Mechanisms of Gastric Carcinogenesis Through Single-Cell RNA Sequencing

Single-cell RNA sequencing (scRNAseq) technology is still relatively new in the field of gastric cancer immunology but gaining significant traction. This technology now provides unprecedented insights into the intratumoral and intertumoral heterogeneities at the immunological, cellular, and molecular levels. Within the last few years, a volume of publications reported the usefulness of scRNAseq technology in identifying thus far elusive immunological mechanisms that may promote and impede gastric cancer development. These studies analyzed datasets generated from primary human gastric cancer tissues, metastatic ascites fluid from gastric cancer patients, and laboratory-generated data from in vitro and in vivo models of gastric diseases. In this review, we overview the exciting findings from scRNAseq datasets that uncovered the role of critical immune cells, including T cells, B cells, myeloid cells, mast cells, ILC2s, and other inflammatory stromal cells, like fibroblasts and endothelial cells. In addition, we also provide a synopsis of the initial scRNAseq findings on the interesting epithelial cell responses to inflammation. In summary, these new studies have implicated roles for T and B cells and subsets like NKT cells in tumor development and progression. The current studies identified diverse subsets of macrophages and mast cells in the tumor microenvironment, however, additional studies to determine their roles in promoting cancer growth are needed. Some groups specifically focus on the less prevalent ILC2 cell type that may contribute to early cancer development. ScRNAseq analysis also reveals that stromal cells, e.g., fibroblasts and endothelial cells, regulate inflammation and promote metastasis, making them key targets for future investigations. While evaluating the outcomes, we also highlight the gaps in the current findings and provide an assessment of what this technology holds for gastric cancer research in the coming years. With scRNAseq technology expanding rapidly, we stress the need for periodic review of the findings and assess the available scRNAseq analytical tools to guide future work on immunological mechanisms of gastric carcinogenesis.  

Pediatric Sarcomas: The Next Generation of Molecular Studies

Simple Summary

There has been an incredible amount of discovery in pediatric sarcomas, but much remains to be accomplished. Clinical challenges include diagnostic heterogeneity and the poor outcome of patients with high risk, metastatic, and relapsed disease. The emergence of single cell sequencing has allowed the ability to document tumor cell heterogeneity in amazing detail, but it does not allow the ability to visualize spatial orientation. This problem has been solved by spatial multi-omics, which can be used to map tumors and visualize the distribution of critical transcripts, mutations, and proteins. However, these tools only offer observational data. High-throughput functional genomics provides a powerful way to highlight oncogenic drivers and potential therapy opportunities. Research has been hamstrung by a need for annotated specimens, particularly in post-therapy, relapsed, and metastatic disease, and initial biopsies offer only limited data opportunities. Data complexity, variability, and inconsistency present problems best approached with AI/machine learning. We stand on the threshold of a revolution in cancer cell biology that has the potential for translation into more effective and more directed therapies, particularly for previously recalcitrant diseases.

Abstract

Pediatric sarcomas constitute one of the largest groups of childhood cancers, following hematopoietic, neural, and renal lesions. Partly because of their diversity, they continue to offer challenges in diagnosis and treatment. In spite of the diagnostic, nosologic, and therapeutic gains made with genetic technology, newer means for investigation are needed. This article reviews emerging technology being used to study human neoplasia and how these methods might be applicable to pediatric sarcomas. Methods reviewed include single cell RNA sequencing (scRNAseq), spatial multi-omics, high-throughput functional genomics, and clustered regularly interspersed short palindromic sequence-Cas9 (CRISPR-Cas9) technology. In spite of these advances, the field continues to be challenged by a dearth of properly annotated materials, particularly from recurrences and metastases and pre- and post-treatment samples.

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.

Spatial components of molecular tissue biology

Methods for profiling RNA and protein expression in a spatially resolved manner are rapidly evolving, making it possible to comprehensively characterize cells and tissues in health and disease. To maximize the biological insights obtained using these techniques, it is critical to both clearly articulate the key biological questions in spatial analysis of tissues and develop the requisite computational tools to address them. Developers of analytical tools need to decide on the intrinsic molecular features of each cell that need to be considered, and how cell shape and morphological features are incorporated into the analysis. Also, optimal ways to compare different tissue samples at various length scales are still being sought. Grouping these biological problems and related computational algorithms into classes across length scales, thus characterizing common issues that need to be addressed, will facilitate further progress in spatial transcriptomics and proteomics.

Mixed Effects Machine Learning Models for Colon Cancer Metastasis Prediction using Spatially Localized Immuno-Oncology Markers

Spatially resolved characterization of the transcriptome and proteome promises to provide further clarity on cancer pathogenesis and etiology, which may inform future clinical practice through classifier development for clinical outcomes. However, batch effects may potentially obscure the ability of machine learning methods to derive complex associations within spatial omics data. Profiling thirty-five stage three colon cancer patients using the GeoMX Digital Spatial Profiler, we found that mixed-effects machine learning (MEML) methods† may provide utility for overcoming significant batch effects to communicate key and complex disease associations from spatial information. These results point to further exploration and application of MEML methods within the spatial omics algorithm development life cycle for clinical deployment.

Mapping Gene Expression in the Spatial Dimension

The main function and biological processes of tissues are determined by the combination of gene expression and spatial organization of their cells. RNA sequencing technologies have primarily interrogated gene expression without preserving the native spatial context of cells. However, the emergence of various spatially-resolved transcriptome analysis methods now makes it possible to map the gene expression to specific coordinates within tissues, enabling transcriptional heterogeneity between different regions, and for the localization of specific transcripts and novel spatial markers to be revealed. Hence, spatially-resolved transcriptome analysis technologies have broad utility in research into human disease and developmental biology. Here, recent advances in spatially-resolved transcriptome analysis methods are summarized, including experimental technologies and computational methods. Strengths, challenges, and potential applications of those methods are highlighted, and perspectives in this field are provided.

Medial Arterial Calcification: JACC State-of-the-Art Review

Medial arterial calcification (MAC) is a chronic systemic vascular disorder distinct from atherosclerosis that is frequently but not always associated with diabetes mellitus, chronic kidney disease, and aging. MAC is also a part of more complex phenotypes in numerous less common diseases. The hallmarks of MAC include disseminated and progressive precipitation of calcium phosphate within the medial layer, a prolonged and clinically silent course, and compromise of hemodynamics associated with chronic limb-threatening ischemia. MAC increases the risk of complications during vascular interventions and mitigates their outcomes. With the exception of rare monogenetic defects affecting adenosine triphosphate metabolism, MAC pathogenesis remains unknown, and causal therapy is not available. Implementation of genetics and omics-based approaches in research recognizing the critical importance of calcium phosphate thermodynamics holds promise to unravel MAC molecular pathogenesis and to provide guidance for therapy. The current state of knowledge concerning MAC is reviewed, and future perspectives are outlined.

Bayesian Modeling of Spatial Molecular Profiling Data via Gaussian Process

MOTIVATION

The location, timing, and abundance of gene expression (both mRNA and proteins) within a tissue define the molecular mechanisms of cell functions. Recent technology breakthroughs in spatial molecular profiling, including imaging-based technologies and sequencing-based technologies, have enabled the comprehensive molecular characterization of single cells while preserving their spatial and morphological contexts. This new bioinformatics scenario calls for effective and robust computational methods to identify genes with spatial patterns.

RESULTS

We represent a novel Bayesian hierarchical model to analyze spatial transcriptomics data, with several unique characteristics. It models the zero-inflated and over-dispersed counts by deploying a zero-inflated negative binomial model that greatly increases model stability and robustness. Besides, the Bayesian inference framework allows us to borrow strength in parameter estimation in a de novo fashion. As a result, the proposed model shows competitive performances in accuracy and robustness over existing methods in both simulation studies and two real data applications.

AVAILABILITY

The related R/C ++ source code is available at https://github.com/Minzhe/BOOST-GP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Recovering Spatially-Varying Cell-Specific Gene Co-expression Networks for Single-Cell Spatial Expression Data

Recent advances in single-cell technologies enable spatial expression profiling at the cell level, making it possible to elucidate spatial changes of cell-specific genomic features. The gene co-expression network is an important feature that encodes the gene-gene marginal dependence structure and allows for the functional annotation of highly connected genes. In this paper, we design a simple and computationally efficient two-step algorithm to recover spatially-varying cell-specific gene co-expression networks for single-cell spatial expression data. The algorithm first estimates the gene expression covariance matrix for each cell type and then leverages the spatial locations of cells to construct cell-specific networks. The second step uses expression covariance matrices estimated in step one and label information from neighboring cells as an empirical prior to obtain thresholded Bayesian posterior estimates. After completing estimates for each cell, this algorithm can further predict or interpolate gene co-expression networks on tissue positions where cells are not captured. In the simulation study, the comparison against the traditional cell-type-specific network algorithms and the cell-specific network method but without incorporating spatial information highlights the advantages of the proposed algorithm in estimation accuracy. We also applied our algorithm to real-world datasets and found some meaningful biological results. The accompanied software is available on https://github.com/jingeyu/CSSN.