Bayesian Hyperspherical Graph Mixture-of-Experts Deciphers Cell–Cell Interaction in Spatial Transcriptomics

Spatial transcriptomics (ST) technologies have transformed our understanding of tissue biology by capturing gene expression with spatial context, enabling systematic analysis of cell-cell interactions (CCIs) and spatial domains in complex tissues. However, existing computational approaches often rely on fixed proximity graphs, curated ligand-receptor databases, or deep graph neural networks that are prone to over-smoothing and lack principled uncertainty quantification. These limitations hinder the discovery of heterogeneous, directional, and long-range CCIs essential for interpreting tissue organization and disease mechanisms. Here, we present B-HGME (Bayesian Hyperspherical Graph Mixture of Experts), a scalable, unsupervised framework that jointly delineates spatial domains and infers CCI networks from ST data with principled uncertainty estimates. B-HGME integrates spatial and gene-regulatory graphs into a dual-scale structure, encodes cell representations on a unit hypersphere via coupled message passing, and decodes edges using a Bayesian mixture-of-experts governed by a Dirichlet-regularized gating network. This design enables the model to capture multi-scale, directional, and biologically coherent interactions while avoiding the over-smoothing and posterior collapse of conventional models. The hyperspherical embedding geometry ensures angular similarity is preserved in high dimensions, and edge-level credibility is derived from the Bayesian posterior, facilitating interpretable and confident CCI inference. Across multiple datasets from six major ST platforms, B-HGME consistently achieves state-of-the-art spatial clustering accuracy and uncovers biologically coherent and diverse CCIs, including novel interactions beyond curated ligand-receptor pairs. B-HGME's hyperspherical embeddings accurately localize canonical astrocytic and laminar markers. Inferred ligand-receptor circuits not only recover known pathways but also reveal previously uncharacterized interactions.