Citations

Abstract Background and Aims Cholangiocarcinoma (CCA) is a rare and aggressive malignancy originating from the bile ducts, characterized by a poor prognosis and limited treatment options. Potentially curative interventions include surgical resection with negative margins (R0) and liver transplantation. However, despite these treatments, the likelihood of post-surgical recurrence remains high, with 50–70% of patients experiencing disease recurrence within five years. This highlights the urgent need for reliable prognostic biomarkers to guide clinical decision-making and personalize postoperative care. Numerous prognostic biomarkers identified in tumor tissue for CCA have been proposed, but many studies suffer from limitations such as single-center designs, small patient cohorts, and semi-quantitative techniques such as immunohistochemistry. These methods often rely on observer interpretation and antibody specificity, which can compromise objectivity and reproducibility. To overcome these limitations, we aimed to identify and validate robust quantitative prognostic biomarkers for CCA by analyzing tumor mRNA expression, integrating them using machine learning (ML), and developing a web-based tool to predict tumor recurrence and overall survival (OS). Method A systematic literature review was conducted to identify candidate prognostic tissue biomarkers for CCA, which were further validated using Cox proportional hazards (PH) OS analyses in five independent transcriptomic cohorts. An international multicenter cohort of 221 patients undergoing tumor resection was assembled, and tumor mRNA expression was analyzed using high-throughput TaqMan OpenArray RT-qPCR. Associations between clinical and molecular features and recurrence-free survival (RFS) and OS were assessed using Cox PH analyses. Six ML survival models were trained and evaluated using nested cross-validation, followed by feature reduction. Shapley Additive exPlanations (SHAP) values were applied to assess feature importance and enhance model interpretability. Results Among 496 candidate CCA prognostic biomarkers identified from the literature, 52 showed prognostic value in Cox PH analyses across transcriptomic cohorts. In the multicenter validation cohort, integrative Cox PH analyses incorporating clinical variables and gene expression identified 11 genes significantly associated with RFS and 25 with OS. The 11 RFS-associated genes, combined with five clinical variables, were used to train six ML survival models, among which the random survival forest (RSF) model demonstrated the highest predictive accuracy. Following training of RSF models across all combinations of four genes and five clinical variables, and subsequent feature importance–based selection, a final model incorporating three genes (FSCN1, LDHA, and SLC2A1) and four clinical variables (CCA subtype, resection margin status, biological sex, and age at surgery) achieved a testing time-dependent AUC of 0.82. This model significantly stratified patients into high- and low recurrence-risk groups, with a median RFS (mRFS) of 8 months and 70 months, respectively. In addition, an OS prediction model using the same feature set achieved a testing AUC of 0.775 and significantly stratified patients according to 5-year OS risk, with a median OS (mOS) of 23 months in the high-risk group and 77 months in the low-risk group. Conclusion We identified and validated a robust panel of prognostic tissue biomarkers in patients with resected CCA. Integration of these biomarkers with clinical variables into ML-based prognostic models enabled accurate prediction of RFS and OS. The resulting web-based tool allows clinicians to input RT-qPCR results for the validated biomarkers and obtain precise, individualized prognostic predictions for postoperative recurrence and survival.

Background and Aims: Cholangiocarcinoma (CCA) is a rare but aggressive malignancy with poorly understood genetic susceptibility. To date, genome-wide association studies (GWAS) investigating germline variants associated with CCA risk remain limited. We aimed to identify genetic risk loci for CCA and its clinical subtypes through comprehensive GWAS and post-GWAS analyses. Approach and Results: We conducted a GWAS of 2366 CCA cases and 11,750 controls of European ancestry. Genome-wide significant loci ( p <5×10 −8 ) were identified and further examined through fine-mapping, functional annotation, and HLA imputation. Subgroup analyses were conducted by CCA subtypes and primary sclerosing cholangitis (PSC) status. Cross-trait linkage disequilibrium score regression and Mendelian randomization were employed to investigate the shared genetic architecture and potential causal relationships with a diverse range of traits. We identified 1 new genome-wide significant variant, rs535777 (OR=1.44), near HLA-DRB1/DQA1 associated with CCA, and 2 variants associated with extrahepatic CCA: rs116224263 (OR=0.17) in LINC02506 at 4p15.1 and rs6914950 (OR=1.63) near HLA-DRB1/DQB1 . Stratified analyses revealed rs2395184 (OR=3.51) near HLA-DRA/DRB5 associated with PSC-related CCA, and rs142674434 (OR=2.98) in THSD7A at 7p21.3 associated with non-PSC-related CCA. HLA imputation uncovered new amino acid residues associated with disease risk. Cross-trait analyses identified shared genetic signals between CCA and anthropometric, lipidemic, lifestyle, and medical traits. Mendelian randomization supported putative causal associations for 12 traits with CCA or its subtypes. Conclusions: Our large-scale GWAS highlights new genetic variants and HLA-linked mechanisms underlying CCA susceptibility. Integrating multi-step post-GWAS approaches enhances understanding of CCA pathogenesis and may facilitate the development of risk biomarkers for early detection and precision prevention strategies.