Abstract
Tumor mutation burden (TMB) has emerged as an essential biomarker for assessing the efficacy of cancer immunotherapy. However, due to the inherent complexity of tumors, TMB is not always correlated with the responsiveness of immune checkpoint inhibitors (ICIs). Thus, refining the interpretation and contextualization of TMB is a requisite for enhancing clinical outcomes. In this study, we conducted a comprehensive investigation of the relationship between TMB and multi-omics data across 33 human cancer types. Our analysis revealed distinct biological changes associated with varying TMB statuses in STAD, COAD, and UCEC. While multi-omics data offer an opportunity to dissect the intricacies of tumors, extracting meaningful biological insights from such massive information remains a formidable challenge. To address this, we developed and implemented the PGLCN, a biologically informed graph neural network based on pathway interaction information. This model facilitates the stratification of patients into subgroups with distinct TMB statuses and enables the evaluation of driver biological processes through enhanced interpretability. By integrating multi-omics data for TMB prediction, our PGLCN model outperformed previous traditional machine learning methodologies, demonstrating superior TMB status prediction accuracy (STAD AUC: 0.976 ± 0.007; COAD AUC: 0.994 ± 0.007; UCEC AUC: 0.947 ± 0.023) and enhanced interpretability (BA-House: 1.0; BA-Community: 0.999; BA-Grid: 0.994; Tree-Cycles: 0.917; Tree-Grids: 0.867). Furthermore, the biological interpretability inherent to PGLCN identified the Toll-like receptor family and DNA repair pathways as potential combined biomarkers in conjunction with TMB status in gastric cancer. This finding suggests a potential synergistic targeting strategy with immunotherapy for gastric cancer, thus advancing the field of precision oncology.