Construction and experimental validation of a macrophage cell senescence-related gene signature to evaluate the prognosis, immunotherapeutic sensitivity, and chemotherapy response in bladder cancer

Tumor-associated macrophages (TAMs) are pivotal components of tumor microenvironment (TME), and senescent TAMs contribute to the alternation of the profiles of TME. However, the potential biological mechanisms and the prognosis value of senescent macrophages are largely unknown, especially in bladder cancer (BLCA). Based on the single-cell RNA sequencing of a primary BLCA sample, 23 macrophage-related genes were identified. Genomic difference analysis, LASSO, and Cox regression were used to develop the risk model. TCGA-BLCA cohort (n = 406) was utilized as the training cohort, and then, three independent cohorts (n = 90, n = 221, n = 165) from Gene Expression Omnibus, clinical samples from the local hospital (n = 27), and in vitro cell experiments were used for external validation. Aldo-keto reductase family 1 member B (AKR1B1), inhibitor of DNA binding 1 (ID1), and transforming growth factor beta 1 (TGFB1I1) were determined and included in the predictive model. The model serves as a promising tool to evaluate the prognosis in BLCA (pooled hazard ratio = 2.51, 95% confidence interval = [1.43; 4.39]). The model was also effective for the prediction of immunotherapeutic sensitivity and chemotherapy treatment outcomes, which were further confirmed by IMvigor210 cohort (P < 0.01) and GDSC dataset, respectively. Twenty-seven BLCA samples from the local hospital proved that the risk model was associated with the malignant degree (P < 0.05). At last, the human macrophage THP-1 and U937 cells were treated with H2O2 to mimic the senescent process in macrophage, and the expressions of these molecules in the model were detected (all P < 0.05).