SPICEMIX enables integrative single-cell spatial modeling of cell identity

Computational modeling for deciphering tissue microenvironment heterogeneity from spatially resolved transcriptomics

Spatial transcriptomics techniques, while measuring gene expression, retain spatial location information, aiding in situ studies of organismal tissue architecture and the progression of pathological processes. These techniques generate vast amounts of omics data, necessitating the development of computational methods to reveal the underlying tissue microenvironment heterogeneity. The main directions in spatial transcriptomics data analysis are spatial domain detection and spatial deconvolution, which can identify spatial functional regions and parse the distribution of cell types in spatial transcriptomics data by integrating single-cell transcriptomics data. In these two research directions, many computational methods have been successively proposed. This article will categorize them into three types: machine learning-based methods, probabilistic models-based methods, and deep learning-based methods. It will list and discuss the representative algorithms of each type along with their advantages and disadvantages and describe the datasets and evaluation metrics used to assess these computational methods, facilitating researchers in selecting suitable computational methods according to their research needs. Finally, combining the latest technological developments and the advantages and disadvantages of current algorithms, this article will look forward to the future directions of computational method development.

Dissecting Spatiotemporal Structures in Spatial Transcriptomics via Diffusion-Based Adversarial Learning

Recent advancements in spatial transcriptomics (ST) technologies offer unprecedented opportunities to unveil the spatial heterogeneity of gene expression and cell states within tissues. Despite these capabilities of the ST data, accurately dissecting spatiotemporal structures (e.g., spatial domains, temporal trajectories, and functional interactions) remains challenging. Here, we introduce a computational framework, PearlST (partial differential equation [PDE]-enhanced adversarial graph autoencoder of ST), for accurate inference of spatiotemporal structures from the ST data using PDE-enhanced adversarial graph autoencoder. PearlST employs contrastive learning to extract histological image features, integrates a PDE-based diffusion model to enhance characterization of spatial features at domain boundaries, and learns the latent low-dimensional embeddings via Wasserstein adversarial regularized graph autoencoders. Comparative analyses across multiple ST datasets with varying resolutions demonstrate that PearlST outperforms existing methods in spatial clustering, trajectory inference, and pseudotime analysis. Furthermore, PearlST elucidates functional regulations of the latent features by linking intercellular ligand-receptor interactions to most contributing genes of the low-dimensional embeddings, as illustrated in a human breast cancer dataset. Overall, PearlST proves to be a powerful tool for extracting interpretable latent features and dissecting intricate spatiotemporal structures in ST data across various biological contexts.

Cross-domain information fusion for enhanced cell population delineation in single-cell spatial-omics data

Cell population delineation and identification is an essential step in single-cell and spatial-omics studies. Spatial-omics technologies can simultaneously measure information from three complementary domains related to this task: expression levels of a panel of molecular biomarkers at single-cell resolution, relative positions of cells, and images of tissue sections, but existing computational methods for performing this task on single-cell spatial-omics datasets often relinquish information from one or more domains. The additional reliance on the availability of "atlas" training or reference datasets limits cell type discovery to well-defined but limited cell population labels, thus posing major challenges for using these methods in practice. Successful integration of all three domains presents an opportunity for uncovering cell populations that are functionally stratified by their spatial contexts at cellular and tissue levels: the key motivation for employing spatial-omics technologies in the first place. In this work, we introduce Cell Spatio- and Neighborhood-informed Annotation and Patterning (CellSNAP), a self-supervised computational method that learns a representation vector for each cell in tissue samples measured by spatial-omics technologies at the single-cell or finer resolution. The learned representation vector fuses information about the corresponding cell across all three aforementioned domains. By applying CellSNAP to datasets spanning both spatial proteomic and spatial transcriptomic modalities, and across different tissue types and disease settings, we show that CellSNAP markedly enhances de novo discovery of biologically relevant cell populations at fine granularity, beyond current approaches, by fully integrating cells' molecular profiles with cellular neighborhood and tissue image information.

GAGE-seq concurrently profiles multiscale 3D genome organization and gene expression in single cells

The organization of mammalian genomes features a complex, multiscale three-dimensional (3D) architecture, whose functional significance remains elusive because of limited single-cell technologies that can concurrently profile genome organization and transcriptional activities. Here, we introduce genome architecture and gene expression by sequencing (GAGE-seq), a scalable, robust single-cell co-assay measuring 3D genome structure and transcriptome simultaneously within the same cell. Applied to mouse brain cortex and human bone marrow CD34+ cells, GAGE-seq characterized the intricate relationships between 3D genome and gene expression, showing that multiscale 3D genome features inform cell-type-specific gene expression and link regulatory elements to target genes. Integration with spatial transcriptomic data revealed in situ 3D genome variations in mouse cortex. Observations in human hematopoiesis unveiled discordant changes between 3D genome organization and gene expression, underscoring a complex, temporal interplay at the single-cell level. GAGE-seq provides a powerful, cost-effective approach for exploring genome structure and gene expression relationships at the single-cell level across diverse biological contexts.

Emerging Roles of Spatial Transcriptomics in Liver Research

Spatial transcriptomics, leveraging sequencing- and imaging-based techniques, has emerged as a groundbreaking technology for mapping gene expression within the complex architectures of tissues. This approach provides an in-depth understanding of cellular and molecular dynamics across various states of healthy and diseased livers. Through the integration of sophisticated bioinformatics strategies, it enables detailed exploration of cellular heterogeneity, transitions in cell states, and intricate cell-cell interactions with remarkable precision. In liver research, spatial transcriptomics has been particularly revelatory, identifying distinct zonated functions of hepatocytes that are crucial for understanding the metabolic and detoxification processes of the liver. Moreover, this technology has unveiled new insights into the pathogenesis of liver diseases, such as the role of lipid-associated macrophages in steatosis and endothelial cell signals in liver regeneration and repair. In the domain of liver cancer, spatial transcriptomics has proven instrumental in delineating intratumor heterogeneity, identifying supportive microenvironmental niches and revealing the complex interplay between tumor cells and the immune system as well as susceptibility to immune checkpoint inhibitors. In conclusion, spatial transcriptomics represents a significant advance in hepatology, promising to enhance our understanding and treatment of liver diseases.

Challenges and perspectives in computational deconvolution of genomics data

Deciphering cell-type heterogeneity is crucial for systematically understanding tissue homeostasis and its dysregulation in diseases. Computational deconvolution is an efficient approach for estimating cell-type abundances from a variety of omics data. Despite substantial methodological progress in computational deconvolution in recent years, challenges are still outstanding. Here we enlist four important challenges related to computational deconvolution: the quality of the reference data, generation of ground truth data, limitations of computational methodologies, and benchmarking design and implementation. Finally, we make recommendations on reference data generation, new directions of computational methodologies, and strategies to promote rigorous benchmarking.

BANKSY unifies cell typing and tissue domain segmentation for scalable spatial omics data analysis

Spatial omics data are clustered to define both cell types and tissue domains. We present Building Aggregates with a Neighborhood Kernel and Spatial Yardstick (BANKSY), an algorithm that unifies these two spatial clustering problems by embedding cells in a product space of their own and the local neighborhood transcriptome, representing cell state and microenvironment, respectively. BANKSY's spatial feature augmentation strategy improved performance on both tasks when tested on diverse RNA (imaging, sequencing) and protein (imaging) datasets. BANKSY revealed unexpected niche-dependent cell states in the mouse brain and outperformed competing methods on domain segmentation and cell typing benchmarks. BANKSY can also be used for quality control of spatial transcriptomics data and for spatially aware batch effect correction. Importantly, it is substantially faster and more scalable than existing methods, enabling the processing of millions of cell datasets. In summary, BANKSY provides an accurate, biologically motivated, scalable and versatile framework for analyzing spatially resolved omics data.

Representing and extracting knowledge from single-cell data

Single-cell analysis is currently one of the most high-resolution techniques to study biology. The large complex datasets that have been generated have spurred numerous developments in computational biology, in particular the use of advanced statistics and machine learning. This review attempts to explain the deeper theoretical concepts that underpin current state-of-the-art analysis methods. Single-cell analysis is covered from cell, through instruments, to current and upcoming models. The aim of this review is to spread concepts which are not yet in common use, especially from topology and generative processes, and how new statistical models can be developed to capture more of biology. This opens epistemological questions regarding our ontology and models, and some pointers will be given to how natural language processing (NLP) may help overcome our cognitive limitations for understanding single-cell data.

MENDER: fast and scalable tissue structure identification in spatial omics data

Tissue structure identification is a crucial task in spatial omics data analysis, for which increasingly complex models, such as Graph Neural Networks and Bayesian networks, are employed. However, whether increased model complexity can effectively lead to improved performance is a notable question in the field. Inspired by the consistent observation of cellular neighborhood structures across various spatial technologies, we propose Multi-range cEll coNtext DEciphereR (MENDER), for tissue structure identification. Applied on datasets of 3 brain regions and a whole-brain atlas, MENDER, with biology-driven design, offers substantial improvements over modern complex models while automatically aligning labels across slices, despite using much less running time than the second-fastest. MENDER's identification power allows the uncovering of previously overlooked spatial domains that exhibit strong associations with brain aging. MENDER's scalability makes it freely appliable on a million-level brain spatial atlas. MENDER's discriminative power enables the differentiation of breast cancer patient subtypes obscured by single-cell analysis.

Charting the tumor microenvironment with spatial profiling technologies

In recent years technologies that can achieve readouts at cellular resolution such as single-cell RNA sequencing (scRNA-seq) have provided a comprehensive characterization of the cellular proportions and phenotypes that populate the tumor microenvironment (TME). However, because of the sample dissociation steps required by these protocols, they fail to capture information related to the intricate spatial context in which cells operate as well as their dense networks of interactions. Spatial profiling technologies have recently emerged as a valuable way to investigate the physical organization of cells crowding the TME in intact tissues. In this review we first discuss how spatial profiling technologies have propelled TME characterization, and then explore their potential to improve both diagnosis and prognosis for cancer patients in the clinic.

Mapping the topography of spatial gene expression with interpretable deep learning

Spatially resolved transcriptomics technologies provide high-throughput measurements of gene expression in a tissue slice, but the sparsity of this data complicates the analysis of spatial gene expression patterns such as gene expression gradients. We address these issues by deriving a topographic map of a tissue slice-analogous to a map of elevation in a landscape-using a novel quantity called the isodepth. Contours of constant isodepth enclose spatial domains with distinct cell type composition, while gradients of the isodepth indicate spatial directions of maximum change in gene expression. We develop GASTON, an unsupervised and interpretable deep learning algorithm that simultaneously learns the isodepth, spatial gene expression gradients, and piecewise linear functions of the isodepth that model both continuous gradients and discontinuous spatial variation in the expression of individual genes. We validate GASTON by showing that it accurately identifies spatial domains and marker genes across several biological systems. In SRT data from the brain, GASTON reveals gradients of neuronal differentiation and firing, and in SRT data from a tumor sample, GASTON infers gradients of metabolic activity and epithelial-mesenchymal transition (EMT)-related gene expression in the tumor microenvironment.

Spatial transcriptomics in development and disease

The proper functioning of diverse biological systems depends on the spatial organization of their cells, a critical factor for biological processes like shaping intricate tissue functions and precisely determining cell fate. Nonetheless, conventional bulk or single-cell RNA sequencing methods were incapable of simultaneously capturing both gene expression profiles and the spatial locations of cells. Hence, a multitude of spatially resolved technologies have emerged, offering a novel dimension for investigating regional gene expression, spatial domains, and interactions between cells. Spatial transcriptomics (ST) is a method that maps gene expression in tissue while preserving spatial information. It can reveal cellular heterogeneity, spatial organization and functional interactions in complex biological systems. ST can also complement and integrate with other omics methods to provide a more comprehensive and holistic view of biological systems at multiple levels of resolution. Since the advent of ST, new methods offering higher throughput and resolution have become available, holding significant potential to expedite fresh insights into comprehending biological complexity. Consequently, a rapid increase in associated research has occurred, using these technologies to unravel the spatial complexity during developmental processes or disease conditions. In this review, we summarize the recent advancement of ST in historical, technical, and application contexts. We compare different types of ST methods based on their principles and workflows, and present the bioinformatics tools for analyzing and integrating ST data with other modalities. We also highlight the applications of ST in various domains of biomedical research, especially development and diseases. Finally, we discuss the current limitations and challenges in the field, and propose the future directions of ST.

SONAR enables cell type deconvolution with spatially weighted Poisson-Gamma model for spatial transcriptomics

Recent advancements in spatial transcriptomic technologies have enabled the measurement of whole transcriptome profiles with preserved spatial context. However, limited by spatial resolution, the measured expressions at each spot are often from a mixture of multiple cells. Computational deconvolution methods designed for spatial transcriptomic data rarely make use of the valuable spatial information as well as the neighboring similarity information. Here, we propose SONAR, a Spatially weighted pOissoN-gAmma Regression model for cell-type deconvolution with spatial transcriptomic data. SONAR directly models the raw counts of spatial transcriptomic data and applies a geographically weighted regression framework that incorporates neighboring information to enhance local estimation of regional cell type composition. In addition, SONAR applies an additional elastic weighting step to adaptively filter dissimilar neighbors, which effectively prevents the introduction of local estimation bias in transition regions with sharp boundaries. We demonstrate the performance of SONAR over other state-of-the-art methods on synthetic data with various spatial patterns. We find that SONAR can accurately map region-specific cell types in real spatial transcriptomic data including mouse brain, human heart and human pancreatic ductal adenocarcinoma. We further show that SONAR can reveal the detailed distributions and fine-grained co-localization of immune cells within the microenvironment at the tumor-normal tissue margin in human liver cancer.

Concurrent profiling of multiscale 3D genome organization and gene expression in single mammalian cells

The organization of mammalian genomes within the nucleus features a complex, multiscale three-dimensional (3D) architecture. The functional significance of these 3D genome features, however, remains largely elusive due to limited single-cell technologies that can concurrently profile genome organization and transcriptional activities. Here, we report GAGE-seq, a highly scalable, robust single-cell co-assay that simultaneously measures 3D genome structure and transcriptome within the same cell. Employing GAGE-seq on mouse brain cortex and human bone marrow CD34+ cells, we comprehensively characterized the intricate relationships between 3D genome and gene expression. We found that these multiscale 3D genome features collectively inform cell type-specific gene expressions, hence contributing to defining cell identity at the single-cell level. Integration of GAGE-seq data with spatial transcriptomic data revealed in situ variations of the 3D genome in mouse cortex. Moreover, our observations of lineage commitment in normal human hematopoiesis unveiled notable discordant changes between 3D genome organization and gene expression, underscoring a complex, temporal interplay at the single-cell level that is more nuanced than previously appreciated. Together, GAGE-seq provides a powerful, cost-effective approach for interrogating genome structure and gene expression relationships at the single-cell level across diverse biological contexts.

BayesTME: An end-to-end method for multiscale spatial transcriptional profiling of the tissue microenvironment

Spatial variation in cellular phenotypes underlies heterogeneity in immune recognition and response to therapy in cancer and many other diseases. Spatial transcriptomics holds the potential to quantify such variation, but existing analysis methods are limited by their focus on individual tasks such as spot deconvolution. We present BayesTME, an end-to-end Bayesian method for analyzing spatial transcriptomics data. BayesTME unifies several previously distinct analysis goals under a single, holistic generative model. This unified approach enables BayesTME to deconvolve spots into cell phenotypes without any need for paired single-cell RNA-seq. BayesTME then goes beyond spot deconvolution to uncover spatial expression patterns among coordinated subsets of genes within phenotypes, which we term spatial transcriptional programs. BayesTME achieves state-of-the-art performance across myriad benchmarks. On human and zebrafish melanoma tissues, BayesTME identifies spatial transcriptional programs that capture fundamental biological phenomena such as bilateral symmetry and tumor-associated fibroblast and macrophage reprogramming. BayesTME is open source.

RETROFIT: Reference-free deconvolution of cell-type mixtures in spatial transcriptomics

Spatial transcriptomics (ST) profiles gene expression in intact tissues. However, ST data measured at each spatial location may represent gene expression of multiple cell types, making it difficult to identify cell-type-specific transcriptional variation across spatial contexts. Existing cell-type deconvolutions of ST data often require single-cell transcriptomic references, which can be limited by availability, completeness and platform effect of such references. We present RETROFIT, a reference-free Bayesian method that produces sparse and interpretable solutions to deconvolve cell types underlying each location independent of single-cell transcriptomic references. Results from synthetic and real ST datasets acquired by Slide-seq and Visium platforms demonstrate that RETROFIT outperforms existing reference-based and reference-free methods in estimating cell-type composition and reconstructing gene expression. Applying RETROFIT to human intestinal development ST data reveals spatiotemporal patterns of cellular composition and transcriptional specificity. RETROFIT is available at https://bioconductor.org/packages/release/bioc/html/retrofit.html.

A comprehensive benchmarking with practical guidelines for cellular deconvolution of spatial transcriptomics

Spatial transcriptomics technologies are used to profile transcriptomes while preserving spatial information, which enables high-resolution characterization of transcriptional patterns and reconstruction of tissue architecture. Due to the existence of low-resolution spots in recent spatial transcriptomics technologies, uncovering cellular heterogeneity is crucial for disentangling the spatial patterns of cell types, and many related methods have been proposed. Here, we benchmark 18 existing methods resolving a cellular deconvolution task with 50 real-world and simulated datasets by evaluating the accuracy, robustness, and usability of the methods. We compare these methods comprehensively using different metrics, resolutions, spatial transcriptomics technologies, spot numbers, and gene numbers. In terms of performance, CARD, Cell2location, and Tangram are the best methods for conducting the cellular deconvolution task. To refine our comparative results, we provide decision-tree-style guidelines and recommendations for method selection and their additional features, which will help users easily choose the best method for fulfilling their concerns.

Present and future challenges for the investigation of transgenerational epigenetic inheritance

Epigenetic pathways are essential in different biological processes and in phenotype-environment interactions in response to different stressors and they can induce phenotypic plasticity. They encompass several processes that are mitotically and, in some cases, meiotically heritable, so they can be transferred to subsequent generations via the germline. Transgenerational Epigenetic Inheritance (TEI) describes the phenomenon that phenotypic traits, such as changes in fertility, metabolic function, or behavior, induced by environmental factors (e.g., parental care, pathogens, pollutants, climate change), can be transferred to offspring generations via epigenetic mechanisms. Investigations on TEI contribute to deciphering the role of epigenetic mechanisms in adaptation, adversity, and evolution. However, molecular mechanisms underlying the transmission of epigenetic changes between generations, and the downstream chain of events leading to persistent phenotypic changes, remain unclear. Therefore, inter-, (transmission of information between parental and offspring generation via direct exposure) and transgenerational (transmission of information through several generations with disappearance of the triggering factor) consequences of epigenetic modifications remain major issues in the field of modern biology. In this article, we review and describe the major gaps and issues still encountered in the TEI field: the general challenges faced in epigenetic research; deciphering the key epigenetic mechanisms in inheritance processes; identifying the relevant drivers for TEI and implement a collaborative and multi-disciplinary approach to study TEI. Finally, we provide suggestions on how to overcome these challenges and ultimately be able to identify the specific contribution of epigenetics in transgenerational inheritance and use the correct tools for environmental science investigation and biomarkers identification.