Bioinformatics analysis of prognosis and immune microenvironment of immunological cell death-related gemcitabine-resistance genes in bladder cancer

BACKGROUND: Bladder cancer (BC) is the most common malignant tumor of the urinary system. Gemcitabine resistance partly accounts for treatment failure and recurrence in BC. Immunological cell death (ICD) is correlated with chemoresistance. The prognosis of patients with similar tumor stage still varies in response to chemotherapy, recurrence, and disease progression. Therefore, our study aimed to provide a prognostic model based on ICD-related and gemcitabine-resistance genes for BC. METHODS: The data of BC patients were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs), and differentially expressed gemcitabine resistance-related genes (DEGRRGs) were identified using the edgeR package. The survival-associated DEGRRGs were identified by univariate Cox analysis. A prognostic model was established by univariate Cox regression analysis and validated by GEO dataset. The outcome of low-risk group and high-risk group was analyzed by the Kaplan-Meier curve. The relationship between risk score and immune cell infiltration was investigated using the TIMER online database. RESULTS: The prognosis of patients in the ICD-high group was significantly better than ICD-low group. A prognostic model containing 5 gemcitabine resistance-related ICD-associated genes, including PTPRR, HOXB3, SIGLEC15, UNC5CL, and CASQ1, was established. In both TCGA prognostic model and GEO validation model, patients in the low-risk group had better outcomes than high-risk group. According to the receiver operating characteristic (ROC) curves, the risk score area under ROC curve (AUC) of the TCGA prognostic model were calculated to be 0.705, while the risk score of the GEO validation model were calculated to be 0.716. Patients in the high-risk group had a significantly higher immune score, stromal score, and infiltration of M0 macrophages, M1 macrophages, M2 macrophages, and activated CD4(+) T cells. Patients in the high-risk group had significantly lower infiltration of the regulatory T cells, resting dendritic cell (DCs), and activated DCs. CONCLUSIONS: The present study highlighted the functional role of gemcitabine resistance-related ICD-associated genes, constructed a prognostic score for the outcome evaluation and searched for potential targets to overcome gemcitabine chemoresistance in BC.