Interpretable and scalable spatial gene set activity analysis with GESSO uncovers functional tissue architecture

Recent advances in spatially resolved transcriptomics (SRT) enabled measurement of sets of pathway genes activity within tissues. However, existing gene set activity scoring methods overlook spatial dependencies among tissue locations, restricting their ability to capture region-specific pathway activities associated with disease pathology or cellular communication. Moreover, these methods lack significance-level inference for activity scores, provide limited interpretability of gene-level contribution to a pathway, and scale poorly to advanced large-size SRT datasets. To address these limitations, we present GESSO (Gene sEt activity Score analysis with Spatial lOcation), a spatially informed gene set scoring method adaptable to diverse SRT platforms. GESSO models gene set activity levels through a graph-regularized matrix decomposition algorithm, jointly inferring spatially coherent gene set activity scores (GASs) and interpretable metagene weights that capture gene-level contributions. It further implements a permutation-based local significance test and a stratified low-resolution approximation that scales to high-resolution SRT datasets such as Visium HD, Stereo-seq, and Xenium Prime. Across 13 datasets from five SRT platforms, GESSO outperformed all existing methods in accuracy, calibration, interpretability, and scalability. Applications revealed novel biological programs, including spatially confined EMT activation within tumor-stroma interfaces, developmental signaling gradients across embryonic tissues, and coordinated B-cell, T-cell, and signaling pathways within germinal centers of human lymph node tissue, revealing the spatial organization of immune function at subregional resolution.