spARC recovers human glioma spatial signaling networks with graph filtering

Abstract

AbstractBiological networks operate within architectural frameworks that influence the state and function of cells through niche-specific factors such as exposure to nutrients and metabolites, soluble signaling molecules, and direct cognate cell-cell communication. Spatial omics technologies incorporate environmental information into the study of biological systems, where the spatial coordinates of cells may directly or indirectly encode these micro-anatomical features. However, they suffer from technical artifacts, such as dropout, that impede biological discovery. Current methods that attempt to correct for this fail to adequately integrate highly informative spatial information when recovering gene expression and modelling cell-cell dynamics in situ. To address this oversight, we developed spatial Affinity-graph Recovery of Counts (spARC), a data diffusion-based filtration method that shares information between neighboring cells in tissue and related cells in expression space, to recover gene dynamics and simulate signalling interactions in spatial transcriptomics data. Following validation, we applied spARC to 10 IDH-mutant surgically resected human gliomas across WHO grades II-IV in order to study signaling networks across disease progression. This analysis revealed co-expressed genes that border the interface between tumor and tumor-infiltrated brain, allowing us to characterize global and local structure of glioma. By simulating paracrine signaling in silico, we identified an Osteopontin-CD44 interaction enriched in grade IV relative to grade II and grade III astrocytomas, and validated the clinical relevance of this signaling axis using TCGA.