bioRxivCancer BiologyPreprint
A Hierarchical Spatial Graph Neural Network Resolves Immunogenic and Tolerogenic Tertiary Lymphoid Structures in Renal Cell Carcinoma
Abstract
Tertiary lymphoid structures (TLS) in the tumour microenvironment span a functional spectrum from immunogenic -- driving germinal centre reactions and anti-tumour immunity -- to tolerogenic, harbouring regulatory T cells and suppressive myeloid populations. Distinguishing these states is clinically critical: immunogenic TLS predict ICI response whereas tolerogenic TLS may promote immune evasion. Bulk transcriptomics conflates productive TLS with exhausted immune infiltrates, masking this distinction. We present a hierarchical graph neural network (GNN) that operates directly on 10x Visium spatial transcriptomics graphs to classify TLS functional state at the cluster level. Using a three-scale architecture combining graph attention (GAT) and differentiable pooling (DiffPool), the model hierarchically aggregates spot-level signals into niche- and region-level representations before predicting immunogenic versus tolerogenic state. Trained on 915 TLS clusters from 24 renal cell carcinoma (RCC) Visium samples (GSE175540), the model achieves a validation AUC-ROC of 0.718 and a clinical AUC of 0.908 on IgG-validated samples from the BIONIKK cohort. Zero-shot transfer to an independent multi-cancer Visium cohort (GSE203612; breast, liver, ovarian, pancreatic, uterine) correctly identifies hepatocellular carcinoma as harbouring the most tolerogenic TLS, consistent with the known immunosuppressive liver tumour microenvironment. Spatial decomposition of CXCL13 across TLS and non-TLS compartments reveals that 85% of tissue CXCL13 signal originates from non-TLS parenchyma, where it co-expresses primarily with exhaustion markers (mean Spearman rho = 0.233) rather than Tfh markers (CXCR5 rho = 0.039) -- a pattern consistent with the paradoxical association of bulk CXCL13 with worse overall survival in TCGA-KIRC (HR = 1.38, p < 0.001). Code and processed data are deposited at GitHub and Zenodo.
- organism: human
- assay: spatial transcriptomics
- research focus: cell type annotation
- research focus: spatial domain detection
- approach: deep learning
- approach: graph neural network
- Citations
- 0
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Spatiotemporal Mapping of Tertiary Lymphoid Structure Heterogeneity Shapes Immune Niches and Clinical Outcomes in Intrahepatic Cholangiocarcinoma
Liaoliao Gao, Jie Mei, … Xuanwen Bao·Science AdvancesPeer-reviewed·Sep 9, 2026
- Spatial Decoding of Tertiary Lymphoid Structure Maturation in Non-Small Cell Lung Cancer Using Deep Neural Networks
Bob Chen, Conrad Foo, … Lisa M. McGinnis·bioRxiv·Jan 20, 2026
- Spatially Organized Tertiary Lymphoid Structures Emerge in Small Cell Lung Cancer and Associate with Improved Survival
Yingying Cao, Michael Nirula, … Anish Thomas·bioRxiv·Aug 11, 2026
- Multi-resolution deep learning characterizestertiary lymphoid structures in solid tumors
M. van Rijthoven, S. Obahor, … K. Silina·bioRxiv·Apr 3, 2023
- Single-cell spatial atlas of tertiary lymphoid structures in ovarian cancer
J. Sarkkinen, A. Junquera Mencia, … A. Farkkila·bioRxiv·May 16, 2023
- The landscape of tertiary lymphoid structures in endometrial cancer revealed through harmonized multi-level transcriptomics
Marta Requesens, Koen Brummel, … Marco de Bruyn·bioRxiv·Aug 20, 2025