Abstract:Abstract Spatially resolved transcriptomics (SRT) technologies simultaneously measure gene expression and spatial location of cells in a 2D tissue slice, enabling the study of spatial patterns of gene expression. Current approaches to model spatial variation in gene expression assume either that gene expression is determined by discrete cell types or that gene expression varies continuously across a tissue slice. However, neither of these modeling assumptions adequately represent spatial variation in gene expression: the first assumption ignores continuous variation within a spatial region containing cells of the same type, while the second assumption does not allow for discontinuous changes in expression, e.g., due to a sharp change in cell type composition. We propose a model of spatial patterns in gene expression that incorporates both discontinuous and continuous spatial variation in gene expression. Specifically, we model the expression of a gene in a layered tissue slice as a piecewise linear function of a single spatial coordinate with potential discontinuities at tissue layer boundaries. We formulate the problem of inferring tissue layer boundaries and gene expression functions for all genes simultaneously. We derive a dynamic programming algorithm to find the optimal boundaries of tissue layers when these boundaries are lines parallel to a coordinate axis. We generalize this algorithm to arbitrarily curved tissue layers by transforming the tissue geometry to be axis-aligned using the theory of conformal maps from complex analysis. We implement these algorithms in a method called Belayer. Applying Belayer to simulated data and to spatial transcriptomics data from the human dorsolateral prefrontal cortex, we demonstrate that Belayer achieves both higher accuracy in clustering tissue layers compared to state-of-the-art SRT clustering methods, and higher accuracy in identifying cortical layer marker genes compared to commonly-used methods for identification of spatially varying genes.

Partial alignment of multislice spatially resolved transcriptomics data

PASTE2: Partial Alignment of Multi-slice Spatially Resolved Transcriptomics Data

Alignment and Integration of Spatial Transcriptomics Data

Revealing Spatial Multimodal Heterogeneity in Tissues with SpaTrio

DeST-OT: Alignment of Spatiotemporal Transcriptomics Data

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Multi-Slice Spatial Transcriptomics Data Integration Analysis with STG3Net

Automated Integration of Multi-Slice Spatial Transcriptomics Data in 2D and 3D

Statistical batch-aware embedded integration, dimension reduction and alignment for spatial transcriptomics

Mitigating autocorrelation during spatially resolved transcriptomics data analysis

STEP: Spatial Transcriptomics Embedding Procedure for Multi-scale Biological Heterogeneities Revelation in Multiple Samples

Probabilistic embedding, clustering, and alignment for integrating spatial transcriptomics data with PRECAST

Statistical and machine learning methods for spatially resolved transcriptomics with histology

SANTO: a coarse-to-fine alignment and stitching method for spatial omics

Belayer: Modeling discrete and continuous spatial variation in gene expression from spatially resolved transcriptomics

STalign: Alignment of spatial transcriptomics data using diffeomorphic metric mapping

Integrating spatially-resolved transcriptomics data across tissues and individuals: challenges and opportunities

Spatial-linked Alignment Tool (SLAT) for Aligning Heterogenous Slices

Enhanced spatially resolved transcriptomics analysis by matching between expression profiles and spatial topology

SPIRAL: integrating and aligning spatially resolved transcriptomics data across different experiments, conditions, and technologies

Latent feature extraction with a prior-based self-attention framework for spatial transcriptomics