Abstract
BACKGROUND: The response to immunotherapy is limited in advanced biliary tract cancer (BTC). Response prediction is a serious challenge in the clinic. METHODS: This study included 60 patients with advanced BTC who received anti-PD-1 treatment. Among these patients, 30 were subjected to 520 gene panel sequencing, and 50 were subjected to multiplex circulating cytokine testing. The entropy and mutation features were analysed via the optimized pipeline based on our previous work. The repeated LASSO algorithm was used to identify the optimal features. The associations between sequence features and cell communications were explored by analysing single-cell transcriptome data from BTC (GSE125449). Cox regression was used to develop the integrated model. Time-dependent C-index, Kaplan‒Meier, and receiver operating characteristic (ROC) curves were used to assess the prediction performance. RESULTS: TP53, NRAS, FBXW7, and APC were identified as prognosis-related genes. The average C-indices of sequence entropy (0.819) and mutation (0.817) for overall survival (OS) were significantly greater than those of tumour mutation burden (TMB, 0.392) and mutation score (0.638). Single-cell transcriptome data revealed that TP53, KRAS, and NRAS were enriched in plasmacytoid dendritic cells (pDCs) and that the communication between pDCs and macrophages was mediated through the CXCL signalling pathway. The integrated model (EM-CXCL10) showed powerful predictive performance for survival status (AUC: 0.863, 95% CI: 0.643-0.972) and objective response rate (AUC: 0.990, 95% CI: 0.822-1.000). CONCLUSIONS: This study constructed a multidimensional strategy that might indicate the prognosis of BTC immunotherapy, enabling the recognition of BTC patients who would benefit from immunotherapy, thereby guiding personalized clinical decision-making.