A Comprehensive Benchmark of Batch Integration Methods for Spatial Transcriptomics Using a Large-Scale Cancer Atlas

Spatial transcriptomics (ST) enables spatially-resolved gene expression measurement, providing insights into tissue architecture and disease biology. However, batch effects from sequencing protocols, sample processing, and other technical factors can confound biological signals. Although batch correction has been extensively studied in single-cell transcriptomics, spatial integration methods lack rigorous benchmarking on large real-world datasets. This study benchmarks 11 representation-learning methods across three categories—linear, graph-based and probabilistic methods using Owkin's MOSAIC Window dataset, a large-scale spatial transcriptomics atlas of human cancers. Methods are evaluated across three criteria: batch correction, biological conservation, and spatial conservation. We also propose a new integration metric to assess robustness of representations to domain shifts and generalizability to unseen samples. Probabilistic methods (scVIVA, scVI) outperform linear and graph-based approaches in batch correction and biological conservation. On the other hand, graph-based methods excel at spatial conservation but underperform in batch integration. Out-of-Distribution (OOD) evaluation reveals that sophisticated methods show reduced performance on unseen samples while linear methods maintain robust generalization, highlighting trade-offs between integration quality and generalizability that should guide method selection for real-world applications.