Abstract:High throughput spatial transcriptomics (HST) is a rapidly emerging class of experimental technologies that allow for profiling gene expression in tissue samples at or near single‐cell resolution while retaining the spatial location of each sequencing unit within the tissue sample. Through analyzing HST data, we seek to identify sub‐populations of cells within a tissue sample that may inform biological phenomena. Existing computational methods either ignore the spatial heterogeneity in gene expression profiles, fail to account for important statistical features such as skewness, or are heuristic‐based network clustering methods that lack the inferential benefits of statistical modeling. To address this gap, we develop SPRUCE: a Bayesian spatial multivariate finite mixture model based on multivariate skew‐normal distributions, which is capable of identifying distinct cellular sub‐populations in HST data. We further implement a novel combination of Pólya–Gamma data augmentation and spatial random effects to infer spatially correlated mixture component membership probabilities without relying on approximate inference techniques. Via a simulation study, we demonstrate the detrimental inferential effects of ignoring skewness or spatial correlation in HST data. Using publicly available human brain HST data, SPRUCE outperforms existing methods in recovering expertly annotated brain layers. Finally, our application of SPRUCE to human breast cancer HST data indicates that SPRUCE can distinguish distinct cell populations within the tumor microenvironment. An R package spruce for fitting the proposed models is available through The Comprehensive R Archive Network (CRAN). This article is protected by copyright. All rights reserved

spVC for the detection and interpretation of spatial gene expression variation

Spanve: an Statistical Method to Detect Clustering-friendly Spatially Variable Genes in Large-scale Spatial Transcriptomics Data

Spanve: A Statistical Method for Detecting Downstream-Friendly Spatially Variable Genes in Large-Scale Spatial Transcriptomic Data

A Bayesian modified Ising model for identifying spatially variable genes from spatial transcriptomics data

Evaluating spatially variable gene detection methods for spatial transcriptomics data

VISTA Uncovers Missing Gene Expression and Spatial-induced Information for Spatial Transcriptomic Data Analysis

Spatial pattern and differential expression analysis with spatial transcriptomic data

SpatialDE2: Fast and localized variance component analysis of spatial transcriptomics

Identification of cell-type-specific spatially variable genes accounting for excess zeros

Benchmarking Computational Integration Methods for Spatial Transcriptomics Data

A selective review of recent developments in spatially variable gene detection for spatial transcriptomics

DESpace : spatially variable gene detection via differential expression testing of spatial clusters

Addressing the mean-variance relationship in spatially resolved transcriptomics data with spoon

PROST: quantitative identification of spatially variable genes and domain detection in spatial transcriptomics

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Benchmarking algorithms for spatially variable gene identification in spatial transcriptomics

Adding Highly Variable Genes to Spatially Variable Genes Can Improve Cell Type Clustering Performance in Spatial Transcriptomics Data

Categorization of 33 computational methods to detect spatially variable genes from spatially resolved transcriptomics data

Analysis and Visualization of Spatial Transcriptomic Data

Differential gene expression analysis of spatial transcriptomic experiments using spatial mixed models