Recent advances in spatially resolved transcriptomics have enabled comprehensive measurements of gene expression patterns while retaining the spatial context of the tissue microenvironment. Deciphering the spatial context of spots in a tissue needs to use their spatial information carefully. To this end, we develop a graph attention auto-encoder framework STAGATE to accurately identify spatial domains by learning low-dimensional latent embeddings via integrating spatial information and gene expression profiles. To better characterize the spatial similarity at the boundary of spatial domains, STAGATE adopts an attention mechanism to adaptively learn the similarity of neighboring spots, and an optional cell type-aware module through integrating the pre-clustering of gene expressions. We validate STAGATE on diverse spatial transcriptomics datasets generated by different platforms with different spatial resolutions. STAGATE could substantially improve the identification accuracy of spatial domains, and denoise the data while preserving spatial expression patterns. Importantly, STAGATE could be extended to multiple consecutive sections to reduce batch effects between sections and extracting three-dimensional (3D) expression domains from the reconstructed 3D tissue effectively.

Whole-body regeneration of planarians is a natural wonder but how it occurs remains elusive. It requires coordinated responses from each cell in the remaining tissue with spatial awareness to regenerate new cells and missing body parts. While previous studies identified new genes essential to regeneration, a more efficient screening approach that can identify regeneration-associated genes in the spatial context is needed. Here, we present a comprehensive three-dimensional spatiotemporal transcriptomic landscape of planarian regeneration. We describe a pluripotent neoblast subtype, and show that depletion of its marker gene makes planarians more susceptible to sub-lethal radiation. Furthermore, we identified spatial gene expression modules essential for tissue development. Functional analysis of hub genes in spatial modules, such as plk1, shows their important roles in regeneration. Our three-dimensional transcriptomic atlas provides a powerful tool for deciphering regeneration and identifying homeostasis-related genes, and provides a publicly available online spatiotemporal analysis resource for planarian regeneration research.