Causal machine learning–guided personalized immunochemotherapy strategies in intrahepatic cholangiocarcinoma (CMICC).
Jun-Hao Mei, Jian Lu, Kai Zhang, Ying Zhang, Zhen Li, Naijian Ge, Xue Han, Gaojun Teng, Rong Liu, Jinghe Guo, Sui-Xing Zhong, Xiu-Ping Zhang, Ping-Ping Wu, Wei-Fu Lv, Jun Wu, Jianxiong Wu, Chao Wang, Hui Yan, Tian Huang, Jie Liu1
Background:
Immunochemotherapy (IO-chemo) has become the first-line standard care for patients with advanced intrahepatic cholangiocarcinoma (iCCA), but the incremental benefit of adding IO varies greatly among individuals, making it urgent to develop systems for identifying high-benefit populations based on individualized treatment effect (ITE) estimation.
Methods:
This multi-cohort study included patients with iCCA who underwent either IO-chemo or chemo. The discovery cohort (2018–2022) included patients from three centers, while an independent cohort (2017–2023) from seven other centers served for external validation. Target trial emulation was employed to obtain unbiased average treatment effect estimation. We developed a causal machine-learning (ML) model (CMICC) to estimate heterogeneous treatment effects for IO-chemo. Based on predicted ITE, patients were stratified into high-benefit, no-to-moderate-benefit, and negative-benefit groups for each treatment. Counterfactual analyses compared overall survival between model-guided treatment selection and treatment actually received. The Qini and TOC curves were used to evaluate model performance. Model interpretability was assessed using SHAP.
Results:
A total of 1485 patients, 696 in the IO-chemo group and 789 in the chemo group, were included in the discovery cohort; 562 patients, 246 in the IO-chemo group and 316 in the chemo group, were included in the external validation cohort. Using a recursive feature addition approach, we selected 17 of 55 multidimensional variables to build the CMICC model in the discovery cohort. In the high-benefit group, compared with chemo, the hazard ratio (HR) for IO-chemo was 0.39 (95% CI, 0.30, 0.52),