Single-cell analysis of chemotherapy-resistant microenvironment identifies a chemo-response biomarker for pancreatic cancer

Ruiling Xu,Jing Shao,Jingyang Liu,Bo Qu,Jing Liu
DOI: https://doi.org/10.21037/jgo-24-93
2024-08-31
Abstract:Background: Pancreatic ductal adenocarcinoma (PDAC) is a highly fatal, aggressive cancer due to limited response chemotherapy. The tumor microenvironment (TME) has emerged as a key player in the development of chemoresistance and in malignant progression. In this study, we hypothesized that chemotherapy response is predictable by chemotherapy-related cell types and their differentially expressed genes (DEGs). Methods: DEGs of chemoresistance cell types were identified via single cell analysis and Wilcoxon test. A chemotherapy response signature was established using a random forest model and validated with public datasets. Bulk cell fraction was analyzed using BayesPrism algorithms. Log-rank test was used to analyze survival of PDAC patients. Results: We found that natural killer (NK) cells, myeloid cells, and erythroid cells were highly infiltrated in chemo-resistant TME. A total of 36 chemoresistance-related DEGs of chemo-resistant cells were identified in chemo-resistant PDAC. Functional enrichment analysis showed that chemoresistance upregulated various inflammation-related pathways, including TGF-β signaling. Based on these features, we constructed a random forest model to predict the response and survival for PDAC patients, which accurately distinguished high-risk and chemoresistant patients with significantly poorer prognosis in both the training and independent validation datasets. Cox regression analysis indicated that predicted labels were an independent prognostic factor in PDAC. Moreover, deconvolution of TME confirmed higher infiltration levels of M2 macrophage and NK cells in predicted chemoresistance. When combined with chemotherapy response related tumor mutations, the model showed excellent ability in predicting chemotherapy response and survival. Conclusions: The TME was closely associated with the chemotherapy response and prognosis of PDAC. Our TME-based random forest model predicted chemotherapy response with complementary knowledge to the response-related genetic mutations.
What problem does this paper attempt to address?