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

An interpretable Bayesian clustering approach with feature selection for analyzing spatially resolved transcriptomics data

Bayesian Nonparametric Clustering with Feature Selection for Spatially Resolved Transcriptomics Data

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Generalized Bayesian nonparametric clustering framework for high-dimensional spatial omics data

Bayesian hidden mark interaction model for detecting spatially variable genes in imaging-based spatially resolved transcriptomics data

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

BayeSMART: Bayesian Clustering of Multi-sample Spatially Resolved Transcriptomics Data

Benchmarking Spatial Clustering Methods with Spatially Resolved Transcriptomics Data

Benchmarking Computational Integration Methods for Spatial Transcriptomics Data

Joint dimension reduction and clustering analysis of single-cell RNA-seq and spatial transcriptomics data

BFAST: Joint Dimension Reduction and Spatial Clustering with Bayesian Factor Analysis for Zero-Inflated Spatial Transcriptomics Data.

Accurately deciphering spatial domains for spatially resolved transcriptomics with stCluster

A comprehensive overview of graph neural network-based approaches to clustering for spatial transcriptomics

Probabilistic cell/domain-type assignment of spatial transcriptomics data with SpatialAnno

Bayesian modeling of spatial molecular profiling data via Gaussian process

Benchmarking cell-type clustering methods for spatially resolved transcriptomics data

Bayesian Integrative Region Segmentation in Spatially Resolved Transcriptomic Studies

Benchmarking clustering, alignment, and integration methods for spatial transcriptomics

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

A Robust Statistical Approach for Finding Informative Spatially Associated Pathways