Prognostic Implications of Immune-Related Gene Pairs Signatures in Bladder Cancer

Abstract

Compelling evidence indicates that immune function is correlated with the prognosis of bladder cancer (BC). Here, we aimed to develop a clinically translatable immune-related gene pairs (IRGPs) prognostic signature to estimate the overall survival (OS) of bladder cancer. From the 251 prognostic-related IRGPs, 37 prognostic-related IRGPs were identified using LASSO regression. We identified IRGPs with the potential to be prognostic markers. The established risk scores divided BC patients into high and low risk score groups, and the survival analysis showed that risk score was related to OS in the TCGA-training set ( $p<0.001$ ; HR = 7.5 [5.3, 10]). ROC curve analysis showed that the AUC for the 1-year, 3-year, and 5-year follow-up was 0.820, 0.883, and 0.879, respectively. The model was verified in the TCGA-testing set and external dataset GSE13507. Multivariate analysis showed that risk score was an independent prognostic predictor in patients with BC. In addition, significant differences were found in gene mutations, copy number variations, and gene expression levels in patients with BC between the high and low risk score groups. Gene set enrichment analysis showed that, in the high-risk score group, multiple immune-related pathways were inhibited, and multiple mesenchymal phenotype-related pathways were activated. Immune infiltration analysis revealed that immune cells associated with poor prognosis of BC were upregulated in the high-risk score group, whereas immune cells associated with a better prognosis of BC were downregulated in the high-risk score group. Other immunoregulatory genes were also differentially expressed between high and low risk score groups. A 37 IRGPs-based risk score signature is presented in this study. This signature can efficiently classify BC patients into high and low risk score groups. This signature can be exploited to select high-risk BC patients for more targeted treatment.