A novel glycogene-related signature for prognostic prediction and immune microenvironment assessment in kidney renal clear cell carcinoma

Abstract

Background: Kidney Renal Clear Cell Carcinoma (KIRC) is a prevalent urinary malignancies worldwide. Glycosylation is a key post-translational modification that is essential in cancer progression. However, its relationship with prognosis, tumour microenvironment (TME), and treatment response in KIRC remains unclear. Method: Expression profiles and clinical data were retrieved from The Cancer Genome Atlas and Gene Expression Omnibus databases. Consensus clustering, Cox regression, and LASSO regression analyses were conducted to develop an optimal glycogene-related signature. The prognostic relevance of this molecular signature was rigorously analyzed, along with its connections to tumour microenvironment (TME), tumour mutation burden, immune checkpoint activity, cancer-immunity cycle regulation, immunomodulatory gene expression patterns, and therapeutic response profiles. Validation was performed using real-world clinical specimens, quantitative PCR (qPCR), and immunohistochemistry (IHC), supported by cohort analyses from the Human Protein Atlas (HPA) database. Results: A glycogene-associated prognostic scoring system was established to categorize patients into risk-stratified subgroups. Patients in the high-risk cohort exhibited significantly poorer survival outcomes (p < 0.001). By incorporating clinicopathological variables into this framework, we established a predictive nomogram demonstrating strong calibration and a concordance index (C-index) of 0.78. The high-risk subgroup displayed elevated immune infiltration scores (p < 0.001), upregulated expression of immune checkpoint-related genes (p < 0.05), and an increased frequency of somatic mutations (p = 0.043). The risk score positively correlated with cancer-immunity cycle activation and immunotherapy-related signals. The high-risk groups also showed associations with T cell exhaustion, immune-activating genes, chemokines, and receptors. Drug sensitivity analysis revealed that low-risk patients were more sensitive to sorafenib, pazopanib, and erlotinib, whereas high-risk individuals responded better to temsirolimus (p < 0.01). qPCR and IHC analyses consistently revealed distinct expression patterns of MX2 and other key genes across the risk groups, further corroborated by the HPA findings. Conclusion: This glycogene-based signature provides a robust tool for predicting prognosis, TME characteristics, and therapeutic responses in KIRC, offering potential clinical utility in patient management.