A glycosylation risk score comprehensively assists the treatment of bladder neoplasm in the real-world cohort, including the tumor microenvironment, molecular and clinical prognosis

Background: Bladder cancer is a common urological cancer associated high significant morbidity and mortality rates. Immunotherapy has emerged as a promising treatment option, although response rates vary among patients. Glycosylation has been implicated in tumorigenesis and immune regulation. However, our current comprehensive understanding of the role of glycosylation in bladder cancer and its clinical implications is limited. Methods: We constructed a training cohort based on the downloaded TCGA-BLCA dataset, while additional datasets (Xiangya cohort, GSE32894, GSE48075, GSE31684, GSE69795 and E-MTAB-1803) from Xiangya hospital, GEO and ArrayExpress database were obtained and used as validation cohorts. To identify glycosylation-related genes associated with prognosis, univariate Cox regression and LASSO regression were performed. A Cox proportional hazards regression model was then constructed to develop a risk score model. The performance of the risk score was assessed in the training cohort using Kaplan-Meier survival curves and ROC curves, and further validated in multiple validation cohorts. Results: We classified patients in the training cohort into two groups based on glycosylation-related gene expression patterns: Cluster 1 and Cluster 2. Prognostic analysis revealed that Cluster 2 had poorer survival outcomes. Cluster 2 also showed higher levels of immune cell presence in the tumor microenvironment and increased activation in key steps of the cancer immune response cycle. We developed an independent prognostic risk score (p < 0.001) and used it to construct an accurate prognostic prediction nomogram. The high glycosylation risk score group exhibited higher tumor immune cell infiltration, enrichment scores in immune therapy-related pathways, and a tendency towards a basal subtype. Conversely, the low-risk score group had minimal immune cell infiltration and tended to have a luminal subtype. These findings were consistent in our real-world Xiangya cohort. Conclusion: This multi-omics glycosylation score based on these genes reliably confirmed the heterogeneity of bladder cancer tumors, predicted the efficacy of immunotherapy and molecular subtypes, optimizing individual treatment decisions.