In silico reconstruction of primary and metastatic tumor architecture using GIS-augmented spatial transcriptomics

The tumor microenvironment (TME) comprises different cell populations that interact, contributing to tumor heterogeneity and therapy response. Spatial transcriptomics offers valuable insights into transcriptional complexity and heterogeneity of the TME. We established Geographic Information System (GIS)-augmented In-Silico Reconstruction of Tumor Architecture (GIS-ROTA), a biologically informed analytic framework that integrates pathway or cell type-based enrichment analysis with local Moran's I to uncover functional spatial domains. In our Visium dataset of primary and metastatic estrogen receptor-positive breast tumor samples, GIS-ROTA revealed extensive co-localization of estrogen response with metabolic pathway gene sets and mutual exclusivity with metastasis-related and specific immune-related pathway gene sets. The novelty of our approach lies in considering biological functions prior to identifying any spatial domains, providing direct interpretability and minimizing the subjectivity of interpreting clusters observed from conventional analytic methods. Overall, our GIS-ROTA framework integrates biological knowledge first, yielding spatial patterns with functional relevance and enabling identification of novel targets for development of therapeutic strategies.