Skip to main content

Genomic characterization and immunotherapy for microsatellite instability-high in cholangiocarcinoma

Abstract

Background

Microsatellite instability-high (MSI-H) is a unique genomic status in many cancers. However, its role in the genomic features and immunotherapy in cholangiocarcinoma (CCA) is unclear. This study aimed to systematically investigate the genomic characterization and immunotherapy efficacy of MSI-H patients with CCA.

Methods

We enrolled 887 patients with CCA in this study. Tumor samples were collected for next-generation sequencing. Differences in genomic alterations between the MSI-H and microsatellite stability (MSS) groups were analyzed. We also investigated the survival of PD-1 inhibitor-based immunotherapy between two groups of 139 patients with advanced CCA.

Results

Differential genetic alterations between the MSI-H and MSS groups included mutations in ARID1A, ACVR2A, TGFBR2, KMT2D, RNF43, and PBRM1 which were enriched in MSI-H groups. Patients with an MSI-H status have a significantly higher tumor mutation burden (TMB) (median 41.7 vs. 3.1 muts/Mb, P < 0.001) and more positive programmed death ligand 1 (PD-L1) expression (37.5% vs. 11.9%, P < 0.001) than those with an MSS status. Among patients receiving PD-1 inhibitor-based therapy, those with MSI-H had a longer median overall survival (OS, hazard ratio (HR) = 0.17, P = 0.001) and progression-free survival (PFS, HR = 0.14, P < 0.001) than patients with MSS. Integrating MSI-H and PD-L1 expression status (combined positive score ≥ 5) could distinguish the efficacy of immunotherapy.

Conclusions

MSI-H status was associated with a higher TMB value and more positive PD-L1 expression in CCA tumors. Moreover, in patients with advanced CCA who received PD-1 inhibitor-based immunotherapy, MSI-H and positive PD-L1 expression were associated with improved both OS and PFS.

Trial registration

This study was registered on ClinicalTrials.gov on 07/01/2017 (NCT03892577).

Peer Review reports

Background

The mismatch repair (MMR) pathway is essential for maintaining DNA replication fidelity, mutation avoidance, and genomic stability [1]. When the MMR pathway is not functioning correctly, cells are unable to correct errors during DNA replication, leading to an inconsistent number of microsatellite nucleotide repeats and instability in microsatellite regions [2]. This deficiency, known as deficient DNA mismatch repair (dMMR), is characterized by microsatellite instability-high (MSI-H) in tumors and has a unique genomic status and tumor microenvironments [3, 4]. Patients with MSI-H cancers, including gastrointestinal (GI) cancer, have experienced excellent effects from programmed cell death protein-1 (PD-1)-based immune checkpoint inhibitor (ICI) immunotherapy [5,6,7]. Based on this evidence, the Food and Drug Administration (FDA) approved programmed cell death protein-1 (PD-1) inhibitors monotherapy (and cytotoxic T lymphocyte-associated antigen 4 inhibitors) for patients with MSI-H/dMMR colorectal cancer and multiple other cancers, regardless of the location [3, 8,9,10,11].

Cholangiocarcinoma (CCA) is an aggressive and fatal cancer that originates in the bile duct epithelium [12, 13]. Unfortunately, the majority (> 65%) of patients with biliary tract cancer (BTC) have unresectable disease and patients with advanced CCA have a poor survival rate by classical systemic chemotherapy [12, 13]. However, recent advancements in next-generation sequencing (NGS) technologies have revealed genomic features that may help CCA patients receive precision medicine or immunotherapy [2, 12,13,14,15,16]. In particular, PD-1 inhibitor-based immunotherapy has shown encouraging efficacy in a subset of CCA. The Phase 3 TOPZA-1 and KEYNOTE-966 trials have demonstrated the survival benefits of adding PD-1/ligand 1 (PD-L1) inhibitors for patients with advanced CCA [17, 18]. The prevalence of MSI-H/dMMR in CCA is relatively low (1%-3%) [19,20,21,22], but these patients have shown promising efficacy for immunotherapy. Among advanced CCA patients with MSI-H/dMMR treated with the PD-1 inhibitor pembrolizumab in the Keynote-158 study could achieve a 40.9% objective response rate (ORR) [23], indicating that patients with MSI-H are highly sensitive to PD-1 inhibitors. Some case reports have also shown that MSI-H in patients with CCA predicts response to ICI therapy [24,25,26].

The impact of MSI-H status on the genomic profile and immunotherapy response in CCA is very important but has not been well characterized. In this study, we systematically analyzed the genomic characteristics and immunotherapy efficacy of MSI-H by comparing microsatellite stability (MSS) status in patients with CCA. Integrating histopathological features to identify potentially beneficiary patients for immunotherapy with a MSS status is also very important.

Methods

Patient characteristics

This study involved 887 patients who were diagnosed with histologically confirmed CCA and treated at the Peking Union Medical College Hospital (PUMCH) from 2017 to 2020 and The Affiliated Hospital of Qingdao University (QDUH) from 2017 to 2019. All MSI-H CCA patients were enriched by hand searches of medical records as much as possible to allow better comparability between the two groups. The study collected samples including tumor tissues and paired tumor-free tissues, and clinical data such as gender, age, tumor location, tumor size, tumor differentiation, lymph node metastasis, and tumor stage from electronic medical records. Before sample disposal, two independent pathologists reviewed all tumor tissues to confirm the pathological diagnoses. Only samples with an estimated tumor purity of > 20% on histopathological assessment were subjected to genomic profiling. The study obtained written informed consent from all participants and was approved by the Institutional Ethics Review Committee at PUMCH and QDUH (NCT03892577; RWSHBC). This study followed the guidelines of the Declaration of Helsinki and Good Clinical Practice.

Genomic profiling

Genomic DNA was prepared from formalin-fixed, paraffin-embedded tumor samples and matched white blood cells (84.9%) or paracancerous/normal tissues (15.1%) using a DNA extraction kit (QIAamp DNA FFPE Tissue Kit), according to the manufacturer’s protocols. NGS-based cancer sequencing YuanSu (CSYS) panel (n = 824) or WES (n = 63) was performed as previously described in the OrigiMed [27,28,29,30]. Tumor samples were sequenced by CSYS panel with a mean coverage of 900 × and matched normal blood samples were sequenced to a mean unique coverage of 300 × . For WES, the mean coverage was 500 × for the tumor sample and 100 × for the matched normal blood samples on the Illumina NovaSeq6000 Platform (Illumina).

Genomic alterations, including single-nucleotide variations (SNVs), short and long insertions/deletions (indels), copy-number variations (CNVs), gene rearrangements, and gene fusions, were identified as follows: BWA (v0.6.2) was used for aligning raw reads to the human genome reference sequence (hg19); Picard (version 1.47) was used for removing PCR duplicates; MuTect (v1.7) was used for identifying SNVs; PINDEL (V0.2.5) was used for identifying indels polymorphisms; and SnpEff (v3.0) was used for annotating the functional impact of SNVs and indels. Control-FREEC (v9.7) was used to identify CNV regions. Gene fusions/rearrangements were detected using an in-house pipeline and assessed using Integrative Genomics Viewer (v2.4).

The study determined the functional significance of variants in genes by examining databases and published literature, including ClinVar, Catalogue of Somatic Mutations in Cancer, and PubMed. The study reported known or likely drivers and recurrent variants. Pathogenic mutations were defined as variants that clearly affect the function of a gene [30, 31]. Synonymous variants and variants with uncertain significance (VUS) (named VUS) were excluded because of lacking clinical significance [32]. Only pathogenic, likely pathogenic, and some meaningful mutations that are clear, potentially functionally, clinically significant, or functionally unknown were reported as clear somatic variants (named Not-VUS) in tumors.

The tumor mutation burden (TMB) was calculated by counting the total number of coding mutations per megabase as report [27, 33]. TMB values higher than 10 muts/Mb were considered TMB-high (TMB-H), and those less than 10 muts/Mb were considered TMB-low (TMB-L) [34]. The indel ratio is defined as the number of (frameshift) indel mutations per megabase in the target region [3, 35]. MSI status was determined based on the MANTIS score [36]. MSI-H is defined as more than 15% of selected microsatellite loci showing unstable in tumors compared to matched peripheral blood [29].

Fifteen cancer pathways (Additional File 1: Table S1), including phosphoinositide 3-kinase (PI3K), homologous recombination repair deficiency (HRD), wingless/integrated (WNT), fibroblast growth factor (FGF), cell cycle (CellCycle), switch/sucrose nonfermentable (SWI/SNF), base excision repair (BER), homologous recombination repair (HRR), MMR, nucleotide excision repair (NER), non-homologous end joining (NHEJ), Fanconi anemia (FA), checkpoint factor (CPF), translesion synthesis (TLS), and DNA damage response (DDR), were analyzed [26, 37].

PD-1 inhibitor-based immunotherapy cohort of cholangiocarcinoma and efficacy evaluation

In patients who received PD-1 inhibitor-based therapy, follow-up was conducted to evaluate the efficacy and safety of the drugs until overall survival (OS) was determined. Further analysis was performed on patients who meet the following conditions: (1) patients should have an Eastern Cooperative Oncology Group (ECOG) performance status (PS) of 0–2 and normal baseline organ functions; (2) patients had at least one measurable lesion used to assess the therapeutic response according to the Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1 [38]; (3) patients with radiologically or histologically confirmed advanced CCA who received PD-1 inhibitors based combination therapy with other agents of cancer treatment [39,40,41].

A total of 139 patients with advanced CCA who received PD-1 inhibitor-based combination therapy were enrolled in this study. PD-1 inhibitor-based therapy may have been heterogeneous for different PD-1 inhibitors (including pembrolizumab or nivolumab, toripalimab, camrelizumab, sintilimab, or tislelizumab every three weeks intravenously) and combination therapies in a real-world setting [39,40,41]. Patients were followed every six to eight weeks. The treatment effects were evaluated according to the RECIST version 1.1 guidelines [38] by professional radiologists at our center who were blinded to the therapeutic outcomes and clinicopathological features. Durable clinical benefit (DCB) was defined as complete response (CR), partial response (PR), or stable disease (SD) for ≥ 24 weeks; other was defined as no durable benefit (NDB) [42]. Progression-free survival (PFS) was defined as the time from initiating anti-PD-1 treatment to the first documented disease progression or death due to any causes, while OS was defined as the time between the start of anti-PD-1 treatment and death due to any causes.

Identification and classification of potentially actionable alterations of cholangiocarcinoma

The OncoKB database was used to determine potentially actionable targets (PATs) of genetic alterations in druggable target [43]. All actionable alterations were classified as levels 1, 2, 3A/B, and 4 as follows. Level 1 means FDA-recognized biomarker predictive of response to an FDA-approved drug in this indication; while Level 2 means standard care biomarker recommended by the National Comprehensive Cancer Network (NCCN) or other professional guidelines predictive of response to an FDA-approved drug in this indication. Level 3A means compelling clinical evidence supports the biomarker as being predictive of response to a drug in this indication; while Level 3B means standard care or investigational biomarker predictive of response to an FDA-approved or investigational drug in another indication. For Level 4, this means there is compelling biological evidence to support that the biomarker can predict drug response [43]. Usually, the alterations of levels 1, 2, and 3A were deemed actionable targets [28, 30].

Immunohistochemistry (IHC) analysis of PD-L1 and MMR expression

PD-L1 status was estimated by IHC staining of FFPE tissue sections using anti-PD-L1 antibodies clone 22C3 (DAKO, cat#m3666) or 28–8 (Abcam, Cat#ab205921) [29, 44]. The specimen was considered PD-L1 positive expression when the combined positive score (CPS) was either ≥ 1 [45,46,47]. IHC analysis was conducted to detect MMR-related proteins, including MLH1 (clone ES05, cat#IR079, DAKO), PMS2 (clone EP51, cat#IR087, DAKO), MSH2 (clone FE11, cat#IR08561-2, DAKO), and MSH6 (clone EP49, cat#IR086, DAKO). Tumors were classified as dMMR when the expression of at least one MMR protein was lost and as MMR-proficient when all four MMR proteins had positive nuclear staining in tumor cells [48].

Statistical analysis

The study assessed normal distribution was assessed using the Kolmogorov–Smirnov test. Fisher’s exact test was utilized to compare categorical variables between multiple groups. For continuous variable comparisons between two groups, a two-tailed unpaired t-test was used when the data were normally distributed; otherwise, Wilcoxon’s test was used. The study identified variables associated with MSI-H using logistic regression analysis. For survival analysis, univariable and backward stepwise multivariable Cox proportional hazards regression models (if univariate P value of < 0.1) were used to calculate hazard ratios (HRs). Kaplan–Meier plots (log-rank tests) were used to describe the prognostic factors related to PFS and OS. The results are reported as HRs and 95% confidence intervals (CIs). The R package “survival” for survival analysis, “forestplot” for forest maps, and “ComplexHeatmap” for profiling heat maps were used. The function “somaticInteractions” from R package “maftools” was used to detect mutually exclusive or co-occurring mutations in MSI-H and MSS groups [49]. All reported P values were two-tailed, and P < 0.05 was considered statistically significant. All statistical analyses were performed using R version 4.1.1.

Results

Clinical characteristics of patients with CCA

A total of 887 patients with CCA were included in this study; 584 (65.8%) had intrahepatic cholangiocarcinoma (ICC) and 303 (34.2%) had extrahepatic cholangiocarcinoma (ECC). The median patient age was 60 years old (range: 19–89 years old). Among these patients, 60.2% were male, and 5.4% (48/887) were identified as MSI-H. The median TMB value was 3.7 (interquartile range [IQR]: 1.9–6.2) muts/Mb. The clinicopathological characteristics of the patients and individual details are shown in Additional File 2: Table S2 and Additional File 1: Table S3.

Compared with MSS patients, MSI-H patients had a significantly higher frequency of positive PD-L1 expression (CPS ≥ 1, 37.5% vs. 11.9%, P < 0.001) and frequency of TMB-H (100% vs. 8.7%, P < 0.001 by Fisher’s exact test) (Additional File 2: Table S2, Additional File 2: Figure S1). No statistically significant associations were observed between MSI status and clinical characteristics such as gender, age, tumor location, tumor size, tumor differentiation, and tumor stage (P > 0.05 by Fisher’s exact test).

Genomic features of patients with CCA with different MSI statuses

In total, 7076 somatic SNVs and 2116 small indels, 272 fusions, and 1225 CNVs were identified. After removing VUS gene alternation, 4154 somatic SNVs, 1403 small indels, 169 fusions, and 993 CNV were retained for further analysis. As for 48 MSI-H tumors, 1176 somatic SNVs, 640 small indels, only 2 fusions, and 10 CNV were identified excluding VUS gene alternation (Additional File 1: Table S4).

The differential gene alterations significant between the MSI-H and MSS groups were ARID1A (73% vs. 18%), ACVR2A (69% vs. 1%), TGFBR2 (58% vs. 3%), KMT2D (54% vs. 4%), RNF43 (54% vs. 3%), PBRM1 (52% vs. 7%), MLH1 (46% vs. 2%), KMT2C (35% vs. 7%), PIK3CA (35% vs. 6%), and NF1 (33% vs. 5%) (Fig. 1A, Additional File 1: Table S5). We also found that the mutated genes ARID1A, ACVR2A, TGFBR2, RNF43, KMT2D, and PBRM1 were enriched with the MSI-H status of patients (Additional File 2: Figure S2). Among 15 pathways in cancer, the differential mutation pathways between the MSI-H and MSS groups included the DDR (98% vs. 60%), SWI/SNF (96% vs. 33%), HRD (81% vs. 22%), CellCycle (79% vs. 62%), and PI3K (75% vs. 24%) (Additional File 1: Table S6, Fig. 1A). In addition, tumors with an MSI-H status had a higher TMB (median 41.7 vs. 3.1 Mut/Mb, P < 0.001; Fig. 1B) and a higher insertion indel ratio (median 0.1698 vs. 0.0588, P < 0.001 by Wilcoxon’s test) compared with those with MSS (Additional File 2: Figure S3).

Fig. 1
figure 1

Genomic mutation landscape in cholangiocarcinoma distinguished by MSI status. A Genetic profile of significantly altered genes and associated clinicopathological characteristics of all 887 cholangiocarcinoma patients. Only the top 20 altered genes and top 10 pathways of the MSI-H group are shown. B The TMB value in the MSI-H group was significantly higher than that in the MSS group. MSI-H, microsatellite instability-high; MSS, microsatellite stability; TMB, tumor mutation burden

Association between the efficiency of PD-1 inhibitor-based treatment and MSI status

In this cohort, 139 CCA patients underwent PD-1 inhibitor-based immunotherapy (NCT03892577). Among them, 16 (11.5%) had tumors with MSI-H (N = 12) or dMMR (N = 4). Three of the four dMMR patients had MLH1 and PMS2 loss, and the one had MSH2 and MSH6 loss. The median age of the 139 patients was 61 years (range, 30–84 years), and 58.3% were male. 107 (77.0%) were ICC and 32 (23.0%) were ECCs. The maximum tumor diameter ≥ 5 cm was 55 (39.6%) patients. The ECOG PS scores of 0, 1, and 2 were 31.7%, 51.8%, and 16.5%, respectively. The proportions of TNM stage I, II, III, and IV were 2.9%, 11.5%, 32.4%, and 53.2%, respectively. The percentages of patients with PD-L1 CPS ≥ 1 and PD-L1 CPS ≥ 5 were 44.6% and 29.5%, respectively (Additional File 2: Table S7). Chemotherapy (N = 24, 17.3%), lenvatinib (N = 82, 59.0%), apatinib (N = 9, 6.5%), and other drugs (N = 24, 17.3%) were used in combination with PD-1 inhibitors. The proportions of patients treated with pembrolizumab, nivolumab, toripalimab, camrelizumab, and other PD-1 inhibitors were 31.7%, 13.7%, 28.1%, 15.1%, and 11.4%, respectively. Compared with the MSS group in the immunotherapy cohort, the MSI-H group had better ECOG PS scores, more PD-L1 CPS ≥ 5, more combinations of PD-1 inhibitors combined with chemotherapy, more stage IV, and more TMB-H (P < 0.05 by Fisher’s exact test).

The median follow-up of the 139 patients was 14.1 (IQR: 7.7–27.4, range 1.07–70.3) months. As of the last follow-up, 100 (71.9%) patients had died, 39 (28.1%) patients were censored, 37 of them were still alive, and two patients were lost to follow-up. The median OS was 14.8 (95% CI: 12.6–16.9) months. We observed longer median OS (not evaluable [NE] vs. 13.5 m, HR = 0.17 (95% CI: 0.06–0.46, P = 0.001) and longer median PFS (NE vs. 4.0 m, HR = 0.14 (95% CI: 0.05–0.34, P < 0.001) in MSI-H patients compared with MSS patients by log-rank tests (Fig. 2A and B). There were 65 (46.8%) patients who reached DCB after PD-1 inhibitor combination therapy, including 87.5% MSI-H patients and 41.5% MSS patients (OR = 2.11, P < 0.001 by Fisher’s exact test; Fig. 2C). The ORR (68.8% vs. 17.9%, OR = 3.84, p < 0.001) and disease control rate (DCR) (100.0% vs. 65.9%, OR = 1.52, P = 0.003 by Fisher’s exact test) were also higher in the MSI-H group than the MSS group (Additional File 2: Figure S4).

Fig. 2
figure 2

Kaplan‒Meier estimates of overall survival (A), progression-free survival (B), clinical benefit histogram (C) of 139 advanced cholangiocarcinoma patients receiving PD-1 inhibitor-based therapy stratified by MSI status. MSI-H, microsatellite instability-high; MSS, microsatellite stability

Univariate and multivariate survival analyses of the CCA cohort receiving PD-1 inhibitor-based therapy

Due to the small number of MSI-H patients, potential benefits from histopathological and clinical factors should be identified to confirm the efficacy of immunotherapy for patients with advanced CCA.

Eleven potential prognostic variables for OS were selected based on univariate Cox analysis, including age, largest tumor size, ECOG PS, treatment line, combination agents, PD-L1 (CPS ≥ 1), PD-L1 (CPS ≥ 5), MSI status, TMB (≥ 10 muts/Mb), SWI/SNF pathway, and MMR pathway (Table 1). In the multivariate analysis by backward stepwise multivariable Cox proportional hazards regression models, ECOG PS (1–2 vs. 0: HR, 3.096; 95% CI, 1.832–5.263, P < 0.001) and treatment line (≥ 2 vs. 0–1: HR,2.247; 95% CI, 1.289–3.922, P = 0.004), largest tumor size (≥ 5 cm vs. < 5 cm: HR,1.635; 95% CI, 1.048–2.552, P = 0.030) were independently associated with a significantly shorter OS. Conversely, PD-L1 expression (CPS ≥ 5 vs. < 5: HR, 0.424; 95% CI, 0.259–0.694, P = 0.012), combination therapy (chemotherapy vs. target therapy: HR,0.403; 95% CI, 0.191–0.850, P = 0.017) and MSI status (MSI-H vs. MSS: HR, 0.265; 95% CI, 0.093–0.754, P = 0.013) were associated with a significantly longer OS (Table 1).

Table 1 Univariate and multivariate analyses of prognostic factors for progression-free survival (PFS) and overall survival (OS) in patients with advanced cholangiocarcinoma

The median PFS was 5.3 (95% CI, 3.7–6.9) months. Multivariate Cox analysis showed that only PD-L1 expression (CPS ≥ 5 vs. < 5: HR,0.467; 95% CI, 0.278–0.784, P = 0.004) and MSI status (MSI-H vs. MSS: HR, 0.198; 95% CI, 0.069–0.563, P = 0.002) were independent predictors of longer PFS (Table 1).

Guidance of MSI-H combined with PD-L1 expression in the immunotherapy of CCA patients

Among the 139 CCA patients receiving PD-1 inhibitor-based immunotherapy, 29.5% (41/139) harbored high PD-L1 expression (CPS ≥ 5). Patients in the PD-L1 CPS ≥ 5 expression group had longer OS, PFS, and higher DCB rates than those in the PD-L1 CPS < 5 expression group (Additional File 2: Figure S5, ABC). Integrating PD-L1 expression (CPS ≥ 5) and MSI status, we found that the OS of the MSI-H, MSS-PD-L1CPS ≥ 5, and MSS-PD-L1 CPS < 5 groups were NE, 20.6 months, and 10.6 months (P < 0.001 by log-rank tests, Fig. 3A), respectively; the PFS of these groups were NE, 8.3 months, and 3.8 months (P < 0.001 by log-rank tests, Fig. 3B), respectively. Their DCB rates were 87.5%, 64.5%, and 32.1% (P < 0.001 by Fisher’s exact test, Fig. 3C), respectively. It showed the MSI-H status combined with PD-L1 expression could well distinguish the efficacy of immunotherapy.

Fig. 3
figure 3

Kaplan–Meier curves for overall survival, progression-free survival and clinical benefit histogram of three classifications stratified according to MSI-H combined with PD-L1 CPS ≥ 5 (AC). MSI-H, microsatellite instability-high; MSS, microsatellite stability; Mut, mutation; NDB, no durable benefit; OS, overall survival; PFS, progression-free survival; PD-L1, programmed death ligand 1; TMB, tumor mutation burden; WT, wild type

Factors contributing to MSI-H status

We analyzed the association between mutations of 15 pathways and important clinicopathological parameters (like age, gender, tumor location), and MSI-H status, to analyze what factors contribute to MSI-H status. Results showed that the SWI/SNF (OR = 38.07), BER (OR = 19.14), MMR (OR = 10.23), HRR (OR = 9.68), PI3K (OR = 6.37), and WNT (OR = 4.65) pathways have contributed to MSI-H status by multivariate logistic regression (P < 0.05). Meanwhile, tumor location of ICC (OR = 0.24, P = 0.049) had borderline contributed to MSI-H status (Fig. 4). When seeing the SWI/SNF, BER, MMR, and HRR pathway mutation, we may reflect a positive correlation with MSI-H status.

Fig. 4
figure 4

Factors contributing to MSI-H status. Using multivariate logistic regression, the association between pathways and clinical factors, and MSI-H status were analyzed. MSI-H: microsatellite instability-high

Potentially actionable targets (PATs) and co-altered mutations in CCA with different MSI statuses

We used the OncoKB database to annotate the PATs and found that all the MSI-H groups had level I evidence of pembrolizumab use; TMB-H was also considered to have level I evidence of pembrolizumab use. So, only 8.7% of patients in the MSS group (with TMB-H) fulfilled level I evidence for the use of pembrolizumab (Fig. 5A and B). Except for the targets of MSI-H and TMB-H, only 2.08% of patients had a PAT gene mutation (≤ level 3A) in the MSI-H group, compared with 10.25% in the MSS group (Fig. 5C and D). In the MSS group, IDH1, KRAS, BRAF mutation, and FGFR2 fusion were the main target genes (Fig. 5E and F).

Fig. 5
figure 5

Potentially actionable alterations in MSI-H and MSS cancer based on the OncoKB database. A, B Level I of the proportions of MSI-H (100% vs. 0%) and TMB-H (100% vs. 8.7%) in MSI-H and MSS cancer. C, E Flow diagram in the upper part shows the list of translational targets for each OncoKB recommendation level in MSI-H, and the lower part presents the MSS cholangiocarcinoma. The widths of the belts indicate different frequencies for each target at every level. D, F The upper and lower pie plot shows the distribution of OncoKB levels for translational targets in patients with MSI-H or MSS cholangiocarcinoma. MSI-H, microsatellite instability-high; MSS, microsatellite stability; TMB-H, tumor mutation burden-high

We further explored the underlying differences of co-altered mutations in MSI-H and MSS groups. Co-occurrence mutation statuses in the MSI-H group were observed between ARID1A:NOTCH3, ARID1A:KRAS, KMT2D:LRP1B, MSH3:OBSCN, and ATM:POLE. However, NF1:LRP1B was identified exclusively. On the other hand, in the MSS group, TP53:CDKN2A, TP53:TERT, TP53:FBXW7, KRAS:SF3B1, and KRAS:CDKN2A were identified co-occurrence; however, TP53 mutants with STK11, ATM, IDH1, or BAP1 mutations; and KRAS mutation with ERBB2, BRAF, NF1, IDH1, BAP1, or TERT mutations were identified exclusively (Additional File 2: Figure S6).

Discussion

Our study comprehensively investigated the genomic alterations and PATs for precision medicine according to MSI status in 887 patients with CCA. Tumors with MSI-H status were associated with a higher TMB value and more positive PD-L1 expression but fewer PATs. Furthermore, by investigating immunotherapy in 139 advanced CCA patients, we found that MSI-H patients had longer PFS and OS than MSS patients receiving PD-1 inhibitor-based therapy. By integrating MSI-H with positive PD-L1 expression, we identified certain patients with advanced CCA who may benefit from PD-1 inhibitor-based immunotherapy.

MSI-H cancers have distinct genomic features, and this study chose non-VUS gene alterations for further analysis to emphasize clinical translational precision medicine. The mutation frequencies of ARID1A, ACVR2A, TGFBR2, KMT2D, and RNF43 were significantly higher in the MSI-H group of our cohort. We found that the TMB value of the MSI-H group was higher than those of the MSS group. Abnormal pathways, including the DDR, SWI/SNF, and HRD, were highly frequent (> 80%) in the MSI-H group. Similarly, Eluri et al. reported that the frequency of MSI-H in 7565 ICC cases was 1.8%, and the median TMB of MSI-H patients was 21.7 muts/Mb. Meanwhile, the study also showed that genomic alterations in TP53 (59.9%), PRBM1 (37.2%), ARID1A (13.9%), and APC (13.9%) were enriched in patients with MSI-H [50]. In another study, RNF43 and KMT2D mutations frequently occurred in CCA patients with MSI-H and TMB-H status [51]. Similarly, an analysis of the genomic landscape of MSI-H in 11,395 tumors across 30 cancer types showed that mutations in ACVR2A (73%), KMT2D (68%), KMT2B (66%), and MMR-related genes were enriched in the MSI-H group [52].

Tumors with MSI-H have dysfunctions that correct errors introduced in microsatellites, which could lead to more frameshift mutations, higher TMB value, and increased neoantigens formation and tumor-specific antigens which are thought to be unknown and novel to the individual immune system [53]. Then, tumors with MSI-H are highly infiltrated with cytotoxic T-cell lymphocytes like CD8 + T-cells with a Tc1 phenotype and CD4 + T-cells with a Th1 phenotype, to raise tumor-specific immune responses, which leads to sensitivity to ICI of MSI-H tumors [54]. Some studies have found that MSI-H in patients with CCA has promising efficacy in immunotherapy, however, the patients’ numbers were small [21, 23,24,25,26, 50, 55, 56]. In our cohort with 139 patients, we found that CCA patients with MSI-H had longer median OS (HR = 0.17, P = 0.001) and PFS (HR = 0.14, P < 0.001) than CCA patients with MSS, indicating that patients with MSI-H may have more opportunities to benefit from immunotherapy. In the KEYNOTE-158 study, 22 advanced CCA patients with MSI-H/dMMR received pembrolizumab treatment and ORR could reach up to 40.9%, with median PFS and median OS of 4.2 (2.1–24.9) and 19.4 (6.5-NR), respectively; ≥ 3 years OS rate exceeded 30% [23]. However, for the 104 advanced BTC patients with MSS status that received pembrolizumab treatment, the ORR was only 5.8%, with median PFS and median OS of only 2.0 (1.9–2.1) and 7.4 (5.5–9.6) months, respectively; the survival rate of ≥ 3 years was only about 10% [57]. These results suggested that patients with MSI-H have more opportunities to benefit from immunotherapy. In our study, MSI-H status was a protective factor for both OS and PFS of immunotherapy in multivariate analysis. However, another study found that MSI-H (N = 15) was not an independent protective factor (P = 0.162) in refractory CCA with PD-L1 positivity (CPS≥ 1 or tumor proportion score (TPS) ≥ 1%) and treated with pembrolizumab, indicating that PD-L1 expression is also another important confounding factor [58].

Our study found that positive PD-L1 expression was relatively high in patients with MSI-H status. In many types of GI cancers, including BTC, colorectal cancer, and gastric cancer, PD-L1 expression was significantly associated with MSI-H status (P < 0.001) [22, 59, 60]. In our study, we observed that the combination of MSI status and PD-L1 expression could distinguish the prognosis of immunotherapy for both OS and PFS. We found patients with PD-L1 CPS ≥ 5 had a better both OS and PFS than patients with CPS < 5 in multivariate analysis, which is consistent with our previous study [61,62,63,64]. However, the significance of PD-L1 expression as a biomarker for BTC immunotherapy remains controversial. Yoon et al. reported that a PD-L1 CPS ≥ 1 does not distinguish well from the immunotherapy response in BTC [45]. Another study showed that the median PFS and OS of BTC patients did not differ according to the CPS cutoff values (1, 5, and 10) (P > 0.05 for all) [47]. Zhou et al. found a correlation between PD-L1 expression in patients with CPS ≥ 5 and CPS < 5 and the efficacy of anlotinib along with the PD-1 inhibitor; but they observed a similar median PFS (6.80 vs. 6.24 months, P > 0.05) [62]. In another prospective trial, the high PD-L1 expression (CPS ≥ 5) group did not show any differences in median PFS compared with the low PD-L1 expression group (high vs. low: 5.0 vs. 6.4 months, p = 0.81) [65]. Therefore, although patients with PD-L1 CPS ≥ 5 may benefit from PD-1 inhibitor immunotherapy for both OS and PFS in our cohort, the significance of PD-L1 expression and cut-off values still need to be confirmed in prospective and large retrospective cohorts.

We found that the TMB value of the MSI-H group was significantly higher than that of the MSS group, which is a consensus and consistent with the results of other studies involving many cancers [66]. In addition, we also found more indel alterations in the MSI-H CCA group, similar to those observed in colorectal cancer [67]. Meanwhile, our results also showed that mutated pathways, such as SWI/SNF, BER, MMR, and HRR, mostly contributed to the MSI-H status. Many studies have found a high frequency of MSI-H mutations in the SWI/SNF pathway, including intestinal and pancreatic cancers [68,69,70]. Wang et al. defined co-mutations (combinations of HRR-BER and HRR-MMR) in the DDR pathway and found that co-mutated DDR was significantly associated with MSI-H and TMB-H status [37]. Mutations in the MMR pathway surely contribute to the MSI-H status [53]. Similarly, Zhou et al. found that the percentage of MSI-H patients in the HRR-mut group was higher than that in the HRR-WT group in patients with colon cancer [71]. However, in the multivariate Cox analysis of immunotherapy efficacy in our cohort, the prognostic values of TMB, SWI/SNF pathway, or MMR pathway mutations were not significant (P > 0.05).

The percentage of altered genes that can be used as therapeutic targets (level ≤ 3A) according to the OncoKB database in the MSI-H group is low, while the percentage of altered genes in MSS group for potential therapeutic targets could be up to 10%, mainly including IDH1 (ivosidenib), KRAS (adagrasib), FGFR2 fusion (infigratinib), and BRAF mutation (dabrafenib plus trametinib). It showed the spectrum of precision medicine is different between the two groups. Nearly 8.7% of the patients belonged to the TMB-H group in the MSS group, while all (100%) of the MSI-H patients belonged to the TMB-H group, which can also be recommended to use the PD-1 inhibitor pembrolizumab (level I). Similarly, Eluri et al. also observed that MSI-H cases had lower frequencies of previously identified, actionable ICC drivers, such as FGFR2 fusions (0% vs. 9%, p = NS), IDH1 (3.7% vs. 14.5%, p < 0.0001), and IDH2 (0% vs. 4.1%, p = 0.007) [50]. Different gene alterations and co-altered patterns may partly explain the low prevalence of PAT for other precision medicines except for ICI in the MSI-H CCA group [51]. MSI-H patients with CCA seem only amenable to ICI-based therapy.

Our study provided a large cohort and a great resource to comprehensively investigate genomic alterations and immunotherapy in patients with CCA. However, this study has some limitations. First, from a retrospective perspective, the “enrichment” process of patients with MSI-H tumors may have overestimated the proportion of the natural prevalence of MSI-H. However, this study focused on assessing genomic features and comparing immunotherapy efficacy between the two groups, so we wanted to enroll as many MSI-H patients as possible. Second, we just compared DNA-level characteristics according to MSI status, we hope to analyze the multi-omics feature (including RNA sequencing and immune marker by IHC) between MSI-H and MSS groups further. Third, in our immunotherapy cohort, the heterogeneity of PD-1 inhibitor-based therapy may stem from variations across different PD-1 inhibitors and combination therapies. But we found MSI status and PD-L1 CPS ≥ 5 could stratify the cohort prognosis in log-rank tests and multivariable Cox analysis. We encourage validating our findings in other large cohorts.

Conclusions

The proportion of MSI-H in CCA was low. The mutation spectrum was different between MSI-H and MSS status. The MSI-H status of patients with CCA was associated with TMB-H and positive PD-L1 expression. Moreover, MSI-H and positive PD-L1 expression were associated with improved both OS and PFS in patients with advanced CCA who received PD-1 inhibitor-based immunotherapy.

Availability of data and materials

All data supporting the findings of this study are available in this article and its online supplementary material files. Further inquiries can be directed to the corresponding author.

Abbreviations

BER:

Base excision repair

BTC:

Biliary tract cancer

CCA:

Cholangiocarcinoma

CellCycle:

Cell cycle

CI:

Confidence interval

CNV:

Copy-number variation

CPF:

Checkpoint factor

CPS:

Combined positive score

CR:

Complete response

CSYS:

Cancer sequencing YuanSu

DCB:

Durable clinical benefit

DCR:

Disease control rate

DDR:

DNA damage response

dMMR:

Deficient DNA mismatch repair

ECC:

Extrahepatic cholangiocarcinoma

ECOG-PS:

Eastern Cooperative Oncology Group-Performance Status

FA:

Fanconi anemia

FDA:

Food and Drug Administration

FGF:

Fibroblast growth factor

GI:

Gastrointestinal

HR:

Hazard ratio

HRD:

Homologous recombination repair deficiency

HRR:

Homologous recombination repair

ICC:

Intrahepatic cholangiocarcinoma

ICI:

Immune checkpoint inhibitor

Indels:

Insertions/deletions

IQR:

Interquartile range

MMR:

Mismatch repair

MSI-H:

Microsatellite instability-high

MSS:

Microsatellite stability

Mut:

Mutation

NA:

Not available

NCCN:

National Comprehensive Cancer Network

NDB:

No durable benefit

NE:

Not evaluable

NER:

Nucleotide excision repair

NGS:

Next-generation sequencing

NHEJ:

Non-homologous end joining

NSCLC:

Non-small cell lung cancer

OR:

Odd ratio

ORR:

Objective response rate

OS:

Overall survival

PAT:

Potentially actionable target

PD:

Progressive disease

PD-1:

Programmed death 1

PD-L1:

Programmed death-ligand 1

PFS:

Progression-free survival

PI3K:

Phosphoinositide 3-kinase

PR:

Partial response

PUMCH:

Peking Union Medical College Hospital

QDUH:

Qingdao University

RECIST:

Response Evaluation Criteria in Solid Tumors

SD:

Stable disease

SNV:

Single-nucleotide variation

SWI/SNF:

Switch/sucrose nonfermentable

TLS:

Translesion synthesis

TMB:

Tumor mutation burden

TNM:

Tumor node metastasis

TPS:

Tumor proportion score

UK:

Unknow

VUS:

Variants with uncertain significance

WNT:

Wingless/integrated

WT:

Wild type

References

  1. Liu D, Keijzers G, Rasmussen LJ. DNA mismatch repair and its many roles in eukaryotic cells. Mutat Res Rev Mutat Res. 2017;773:174–87.

    Article  PubMed  CAS  Google Scholar 

  2. Eso Y, Shimizu T, Takeda H, Takai A, Marusawa H. Microsatellite instability and immune checkpoint inhibitors: toward precision medicine against gastrointestinal and hepatobiliary cancers. J Gastroenterol. 2020;55(1):15–26.

    Article  PubMed  CAS  Google Scholar 

  3. Mandal R, Samstein RM, Lee KW, Havel JJ, Wang H, Krishna C, Sabio EY, Makarov V, Kuo F, Blecua P, et al. Genetic diversity of tumors with mismatch repair deficiency influences anti-PD-1 immunotherapy response. Science. 2019;364(6439):485–91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Wang J, Xiu J, Farrell A, Baca Y, Arai H, Battaglin F, Kawanishi N, Soni S, Zhang W, Millstein J, et al. Mutational analysis of microsatellite-stable gastrointestinal cancer with high tumour mutational burden: a retrospective cohort study. Lancet Oncol. 2023;24(2):151–61.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Wu T, Zhang X, Liu X, Cai X, Shen T, Pan D, Liang R, Ding R, Hu R, Dong J, et al. Single-cell sequencing reveals the immune microenvironment landscape related to anti-PD-1 resistance in metastatic colorectal cancer with high microsatellite instability. BMC Med. 2023;21(1):161.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Wang Z, Wang X, Xu Y, Li J, Zhang X, Peng Z, Hu Y, Zhao X, Dong K, Zhang B, et al. Mutations of PI3K-AKT-mTOR pathway as predictors for immune cell infiltration and immunotherapy efficacy in dMMR/MSI-H gastric adenocarcinoma. BMC Med. 2022;20(1):133.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Long J, Wang D, Yang X, Wang A, Lin Y, Zheng M, Zhang H, Sang X, Wang H, Hu K, et al. Identification of NOTCH4 mutation as a response biomarker for immune checkpoint inhibitor therapy. BMC Med. 2021;19(1):154.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Le DT, Durham JN, Smith KN, Wang H, Bartlett BR, Aulakh LK, Lu S, Kemberling H, Wilt C, Luber BS, et al. Mismatch repair deficiency predicts response of solid tumors to PD-1 blockade. Science. 2017;357(6349):409–13.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Le DT, Uram JN, Wang H, Bartlett BR, Kemberling H, Eyring AD, Skora AD, Luber BS, Azad NS, Laheru D, et al. PD-1 Blockade in Tumors with Mismatch-Repair Deficiency. New Engl J Med. 2015;372(26):2509–20.

    Article  PubMed  CAS  Google Scholar 

  10. André T, Shiu K-K, Kim TW, Jensen BV, Jensen LH, Punt C, Smith D, Garcia-Carbonero R, Benavides M, Gibbs P, et al. Pembrolizumab in Microsatellite-Instability–High Advanced Colorectal Cancer. New Engl J Med. 2020;383(23):2207–18.

    Article  PubMed  Google Scholar 

  11. Overman MJ, McDermott R, Leach JL, Lonardi S, Lenz HJ, Morse MA, Desai J, Hill A, Axelson M, Moss RA, et al. Nivolumab in patients with metastatic DNA mismatch repair-deficient or microsatellite instability-high colorectal cancer (CheckMate 142): an open-label, multicentre, phase 2 study. Lancet Oncol. 2017;18(9):1182–91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Valle JW, Lamarca A, Goyal L, Barriuso J, Zhu AX. New Horizons for Precision Medicine in Biliary Tract Cancers. Cancer Discov. 2017;7(9):943–62.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Lin Y, Peng L, Dong L, Liu D, Ma J, Lin J, Chen X, Lin P, Song G, Zhang M, et al. Geospatial Immune Heterogeneity Reflects the Diverse Tumor-Immune Interactions in Intrahepatic Cholangiocarcinoma. Cancer Discov. 2022;12(10):2350–71.

    Article  PubMed  CAS  Google Scholar 

  14. Jusakul A, Cutcutache I, Yong CH, Lim JQ, Huang MN, Padmanabhan N, Nellore V, Kongpetch S, Ng AWT, Ng LM, et al. Whole-Genome and Epigenomic Landscapes of Etiologically Distinct Subtypes of Cholangiocarcinoma. Cancer Discov. 2017;7(10):1116–35.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Zou S, Li J, Zhou H, Frech C, Jiang X, Chu JS, Zhao X, Li Y, Li Q, Wang H, et al. Mutational landscape of intrahepatic cholangiocarcinoma. Nat Commun. 2014;5:5696.

    Article  PubMed  CAS  Google Scholar 

  16. Qiu Z, Ji J, Xu Y, Zhu Y, Gao C, Wang G, Li C, Zhang Y, Zhao J, Wang C, et al. Common DNA methylation changes in biliary tract cancers identify subtypes with different immune characteristics and clinical outcomes. BMC Med. 2022;20(1):64.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Oh D-Y, Ruth He A, Qin S, Chen L-T, Okusaka T, Vogel A, Kim Jin W, Suksombooncharoen T, Ah Lee M, Kitano M et al: Durvalumab plus Gemcitabine and Cisplatin in Advanced Biliary Tract Cancer. NEJM Evidence 2022, 1(8):EVIDoa2200015.

  18. Kelley RK, Ueno M, Yoo C, Finn RS, Furuse J, Ren Z, Yau T, Klümpen H-J, Chan SL, Ozaka M et al: Pembrolizumab in combination with gemcitabine and cisplatin compared with gemcitabine and cisplatin alone for patients with advanced biliary tract cancer (KEYNOTE-966): a randomised, double-blind, placebo-controlled, phase 3 trial. Lancet 2023.

  19. Suda R, Sakai N, Matsushita K, Ishige T, Kawasaki Y, Shiko Y, Furukawa K, Mishima T, Nakadai E, Ohtsuka M. Prediction of mismatch repair deficient biliary tract cancer: Role of morphological features and host immune response detected by routine hematoxylin-eosin staining. J Hepatobiliary Pancreat Sci. 2021;28(8):680–91.

    Article  PubMed  Google Scholar 

  20. Goeppert B, Roessler S, Renner M, Singer S, Mehrabi A, Vogel MN, Pathil A, Czink E, Kohler B, Springfeld C, et al. Mismatch repair deficiency is a rare but putative therapeutically relevant finding in non-liver fluke associated cholangiocarcinoma. Br J Cancer. 2019;120(1):109–14.

    Article  PubMed  CAS  Google Scholar 

  21. Yu J, Zhang X, Huang Q, Tan S, Xiong X, Gou H. Rare DNA Mismatch Repair-Related Protein Loss in Patients with Intrahepatic Cholangiocarcinoma and Combined Hepatocellular-Cholangiocarcinoma and Their Response to Immunotherapy. Cancer Manag Res. 2021;13:4283–90.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Mody K, Starr J, Saul M, Poorman K, Weinberg BA, Salem ME, VanderWalde A, Shields AF. Patterns and genomic correlates of PD-L1 expression in patients with biliary tract cancers. J Gastrointest Oncol. 2019;10(6):1099–109.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Maio M, Ascierto PA, Manzyuk L, Motola-Kuba D, Penel N, Cassier PA, Bariani GM, De Jesus Acosta A, Doi T, Longo F, et al. Pembrolizumab in microsatellite instability high or mismatch repair deficient cancers: updated analysis from the phase II KEYNOTE-158 study. Ann Oncol. 2022;33(9):929–38.

    Article  PubMed  CAS  Google Scholar 

  24. Kai Y, Ikezawa K, Takada R, Daiku K, Maeda S, Abe Y, Yamai T, Fukutake N, Nakabori T, Uehara H, et al. Success rate of microsatellite instability examination and complete response with pembrolizumab in biliary tract cancer. JGH Open. 2021;5(6):712–6.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Niger M, Nichetti F, Dell’Angelo F, Cappelletti V, Pircher C, Vismara M, Cotsoglou C, Bhoori S, Vingiani A, Di Bartolomeo M, et al. Acquired Resistance Mechanisms to PD-L1 Blockade in a Patient With Microsatellite Instability-High Extrahepatic Cholangiocarcinoma. JCO Precis Oncol. 2022;6:e2100472.

    Article  PubMed  Google Scholar 

  26. Shi GM, Huang XY, Wu D, Sun HC, Liang F, Ji Y, Chen Y, Yang GH, Lu JC, Meng XL, et al. Toripalimab combined with lenvatinib and GEMOX is a promising regimen as first-line treatment for advanced intrahepatic cholangiocarcinoma: a single-center, single-arm, phase 2 study. Signal Transduct Target Ther. 2023;8(1):106.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Cao J, Chen L, Li H, Chen H, Yao J, Mu S, Liu W, Zhang P, Cheng Y, Liu B, et al. An Accurate and Comprehensive Clinical Sequencing Assay for Cancer Targeted and Immunotherapies. Oncologist. 2019;24(12):e1294–302.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Lin J, Cao Y, Yang X, Li G, Shi Y, Wang D, Long J, Song Y, Mao J, Xie F, et al. Mutational spectrum and precision oncology for biliary tract carcinoma. Theranostics. 2021;11(10):4585–98.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Wu L, Yao H, Chen H, Wang A, Guo K, Gou W, Yu Y, Li X, Yao M, Yuan S, et al. Landscape of somatic alterations in large-scale solid tumors from an Asian population. Nat Commun. 2022;13(1):4264.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Lin J, Shi J, Guo H, Yang X, Jiang Y, Long J, Bai Y, Wang D, Yang X, Wan X, et al. Alterations in DNA Damage Repair Genes in Primary Liver Cancer. Clin Cancer Res. 2019;25(15):4701–11.

    Article  PubMed  CAS  Google Scholar 

  31. Walsh R, Mazzarotto F, Whiffin N, Buchan R, Midwinter W, Wilk A, Li N, Felkin L, Ingold N, Govind R, et al. Quantitative approaches to variant classification increase the yield and precision of genetic testing in Mendelian diseases: the case of hypertrophic cardiomyopathy. Genome Med. 2019;11(1):5.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Chen N, Nanda R, Howard FM, Vidula N, Yen J, Drusbosky LM, Bucheit L: Abstract P5–03–18: Co-occurring alterations in PALB2 germline carriers identified by liquid biopsy in patients with advanced breast cancer. Cancer Res 2023, 83(5_Supplement):P5–03–18-P05–03–18.

  33. Chalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, Schrock A, Campbell B, Shlien A, Chmielecki J, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med. 2017;9(1):34.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Marabelle A, Fakih M, Lopez J, Shah M, Shapira-Frommer R, Nakagawa K, Chung HC, Kindler HL, Lopez-Martin JA, Miller WH Jr, et al. Association of tumour mutational burden with outcomes in patients with advanced solid tumours treated with pembrolizumab: prospective biomarker analysis of the multicohort, open-label, phase 2 KEYNOTE-158 study. Lancet Oncol. 2020;21(10):1353–65.

    Article  PubMed  CAS  Google Scholar 

  35. Kim HD, Ryu MH, Park YS, Lee SY, Moon M, Kang YK. Insertion-deletion rate is a qualitative aspect of the tumor mutation burden associated with the clinical outcomes of gastric cancer patients treated with nivolumab. Gastric Cancer. 2022;25(1):226–34.

    Article  PubMed  CAS  Google Scholar 

  36. Kautto EA, Bonneville R, Miya J, Yu L, Krook MA, Reeser JW, Roychowdhury S. Performance evaluation for rapid detection of pan-cancer microsatellite instability with MANTIS. Oncotarget. 2017;8(5):7452–63.

    Article  PubMed  Google Scholar 

  37. Wang Z, Zhao J, Wang G, Zhang F, Zhang Z, Zhang F, Zhang Y, Dong H, Zhao X, Duan J, et al. Comutations in DNA Damage Response Pathways Serve as Potential Biomarkers for Immune Checkpoint Blockade. Cancer Res. 2018;78(22):6486–96.

    Article  PubMed  CAS  Google Scholar 

  38. Schwartz LH, Litière S, de Vries E, Ford R, Gwyther S, Mandrekar S, Shankar L, Bogaerts J, Chen A, Dancey J, et al. RECIST 1.1-Update and clarification: From the RECIST committee. Eur J Cancer. 2016;62:132–7.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Deng M, Li S, Wang Q, Zhao R, Zou J, Lin W, Mei J, Wei W, Guo R. Real-world outcomes of patients with advanced intrahepatic cholangiocarcinoma treated with programmed cell death protein-1-targeted immunotherapy. Ann Med. 2022;54(1):803–11.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Yang X, Hu Y, Yang K, Wang D, Lin J, Long J, Xie F, Mao J, Bian J, Guan M et al: Cell-free DNA copy number variations predict efficacy of immune checkpoint inhibitor-based therapy in hepatobiliary cancers. J Immunother Cancer 2021, 9(5).

  41. Yang X, Chen B, Wang Y, Wang Y, Long J, Zhang N, Xue J, Xun Z, Zhang L, Cheng J et al: Real-world efficacy and prognostic factors of lenvatinib plus PD-1 inhibitors in 378 unresectable hepatocellular carcinoma patients. Hepatol Int 2023:1–11.

  42. Quispel-Janssen J, van der Noort V, de Vries JF, Zimmerman M, Lalezari F, Thunnissen E, Monkhorst K, Schouten R, Schunselaar L, Disselhorst M, et al. Programmed Death 1 Blockade With Nivolumab in Patients With Recurrent Malignant Pleural Mesothelioma. J Thorac Oncol. 2018;13(10):1569–76.

    Article  PubMed  Google Scholar 

  43. Chakravarty D, Gao J, Phillips SM, Kundra R, Zhang H, Wang J, Rudolph JE, Yaeger R, Soumerai T, Nissan MH et al: OncoKB: A Precision Oncology Knowledge Base. JCO Precis Oncol 2017, 2017.

  44. Lee SH, Lee HS, Lee SH, Woo SM, Kim DU, Bang S. Efficacy and Safety of Pembrolizumab for Gemcitabine/Cisplatin-Refractory Biliary Tract Cancer: A Multicenter Retrospective Study. J Clin Med. 2020;9(6):1769.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Yoon JG, Kim MH, Jang M, Kim H, Hwang HK, Kang CM, Lee WJ, Kang B, Lee CK, Lee MG, et al. Molecular Characterization of Biliary Tract Cancer Predicts Chemotherapy and Programmed Death 1/Programmed Death-Ligand 1 Blockade Responses. Hepatology. 2021;74(4):1914–31.

    Article  PubMed  CAS  Google Scholar 

  46. Arkenau HT, Martin-Liberal J, Calvo E, Penel N, Krebs MG, Herbst RS, Walgren RA, Widau RC, Mi G, Jin J, et al. Ramucirumab Plus Pembrolizumab in Patients with Previously Treated Advanced or Metastatic Biliary Tract Cancer: Nonrandomized, Open-Label, Phase I Trial (JVDF). Oncologist. 2018;23(12):1407-e1136.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Kim H, Kim R, Jo H, Kim HR, Hong J, Ha SY, Park JO, Kim ST. Expression of PD-L1 as a predictive marker of sensitivity to immune checkpoint inhibitors in patients with advanced biliary tract cancer. Therap Adv Gastroenterol. 2022;15:17562848221117638.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  48. Shi J, Yang X, Wang X, Luo Y, Zhou W, Luo H, Bianba Z, Nima Z, Wang Q, Wang H, et al. Prevalence of Epstein-Barr Virus Infection and Mismatch Repair Protein Deficiency and the Correlation of Immune Markers in Tibetan Patients with Gastric Cancer. Biomed Res Int. 2022;2022:2684065.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Gerstung M, Pellagatti A, Malcovati L, Giagounidis A, Porta MG, Jädersten M, Dolatshad H, Verma A, Cross NC, Vyas P, et al. Combining gene mutation with gene expression data improves outcome prediction in myelodysplastic syndromes. Nat Commun. 2015;6:5901.

    Article  PubMed  CAS  Google Scholar 

  50. Eluri M, Danner De Armas A, Ross JS, Sharaf R, Pavlick DC, Chatterjee D, Javle MM. Clinical and mutational profiles of microsatellite instability (MSI-H) in intrahepatic cholangiocarcinoma (iCCA). J Clin Oncol. 2023;41(4_suppl):607–607.

    Article  Google Scholar 

  51. Kendre G, Murugesan K, Brummer T, Segatto O, Saborowski A, Vogel A. Charting co-mutation patterns associated with actionable drivers in intrahepatic cholangiocarcinoma. J Hepatol. 2023;78(3):614–26.

    Article  PubMed  CAS  Google Scholar 

  52. Li Z, Jia Y, Zhu H, Yuan H, Xing X, Xin Y, Ma T, Pang F, Zhang Y, Hu Y, et al. Genomic landscape of microsatellite instability in Chinese tumors: A comparison of Chinese and TCGA cohorts. Int J Cancer. 2022;151(8):1382–93.

    Article  PubMed  CAS  Google Scholar 

  53. Sinicrope FA, Sargent DJ. Molecular pathways: microsatellite instability in colorectal cancer: prognostic, predictive, and therapeutic implications. Clin Cancer Res. 2012;18(6):1506–12.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  54. Taieb J, Svrcek M, Cohen R, Basile D, Tougeron D, Phelip JM. Deficient mismatch repair/microsatellite unstable colorectal cancer: Diagnosis, prognosis and treatment. Eur J Cancer. 2022;175:136–57.

    Article  PubMed  CAS  Google Scholar 

  55. Ueno M, Ikeda M, Morizane C, Kobayashi S, Ohno I, Kondo S, Okano N, Kimura K, Asada S, Namba Y, et al. Nivolumab alone or in combination with cisplatin plus gemcitabine in Japanese patients with unresectable or recurrent biliary tract cancer: a non-randomised, multicentre, open-label, phase 1 study. Lancet Gastroenterol Hepatol. 2019;4(8):611–21.

    Article  PubMed  Google Scholar 

  56. Cheng Z, Zeng T, Yang G, Liu D, Zheng Z, Yuan Z: Genomic and Immune Features in an Intrahepatic Cholangiocarcinoma Patient with Microsatellite Instability-High Suffered Rapid Acquired Resistance to PD-1 Inhibitor. Liver Cancer 2023.

  57. Piha-Paul SA, Oh DY, Ueno M, Malka D, Chung HC, Nagrial A, Kelley RK, Ros W, Italiano A, Nakagawa K, et al. Efficacy and safety of pembrolizumab for the treatment of advanced biliary cancer: Results from the KEYNOTE-158 and KEYNOTE-028 studies. Int J Cancer. 2020;147(8):2190–8.

    Article  PubMed  CAS  Google Scholar 

  58. Park CS, Sung MJ, Kim SJ, Jo JH, Lee HS, Chung MJ, Bang S, Park SW, Song SY, Park JY. Prognostic Factors in Patients Treated with Pembrolizumab as a Second-Line Treatment for Advanced Biliary Tract Cancer. Cancers (Basel). 2022;14(17):1109292.

    Article  Google Scholar 

  59. Rosenbaum MW, Bledsoe JR, Morales-Oyarvide V, Huynh TG, Mino-Kenudson M. PD-L1 expression in colorectal cancer is associated with microsatellite instability, BRAF mutation, medullary morphology and cytotoxic tumor-infiltrating lymphocytes. Mod Pathol. 2016;29(9):1104–12.

    Article  PubMed  CAS  Google Scholar 

  60. Morihiro T, Kuroda S, Kanaya N, Kakiuchi Y, Kubota T, Aoyama K, Tanaka T, Kikuchi S, Nagasaka T, Nishizaki M, et al. PD-L1 expression combined with microsatellite instability/CD8+ tumor infiltrating lymphocytes as a useful prognostic biomarker in gastric cancer. Sci Rep. 2019;9(1):4633.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Mastracci L, Grillo F, Parente P, Gullo I, Campora M, Angerilli V, Rossi C, Sacramento ML, Pennelli G, Vanoli A, et al. PD-L1 evaluation in the gastrointestinal tract: from biological rationale to its clinical application. Pathologica. 2022;114(5):352–64.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Zhou J, Sun Y, Zhang W, Yuan J, Peng Z, Wang W, Gong J, Yang L, Cao Y, Zhao H, et al. Phase Ib study of anlotinib combined with TQB2450 in pretreated advanced biliary tract cancer and biomarker analysis. Hepatology. 2023;77(1):65–76.

    Article  PubMed  Google Scholar 

  63. Ning C, Zhang X, Yang X, Sang X, Zhao H, Wang H. Conversion therapy of stage IVb unresectable gallbladder carcinoma. Hepatobiliary Surg Nutr. 2022;11(2):335–7.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Lin J, Yang X, Long J, Zhao S, Mao J, Wang D, Bai Y, Bian J, Zhang L, Yang X et al: Pembrolizumab combined with lenvatinib as non-first-line therapy in patients with refractory biliary tract carcinoma. Hepatobiliary Surgery and Nutrition; Publish Ahead of Print 2020.

  65. Zeng TM, Yang G, Lou C, Wei W, Tao CJ, Chen XY, Han Q, Cheng Z, Shang PP, Dong YL, et al. Clinical and biomarker analyses of sintilimab plus gemcitabine and cisplatin as first-line treatment for patients with advanced biliary tract cancer. Nat Commun. 2023;14(1):1340.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. Parikh AR, He Y, Hong TS, Corcoran RB, Clark JW, Ryan DP, Zou L, Ting DT, Catenacci DV, Chao J, et al. Analysis of DNA Damage Response Gene Alterations and Tumor Mutational Burden Across 17,486 Tubular Gastrointestinal Carcinomas: Implications for Therapy. Oncologist. 2019;24(10):1340–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  67. Kim JE, Chun S-M, Hong YS. Kim K-p, Kim SY, Kim J, Sung CO, Cho EJ, Kim TW, Jang SJ: Mutation Burden and I Index for Detection of Microsatellite Instability in Colorectal Cancer by Targeted Next-Generation Sequencing. J Mol Diagn. 2019;21(2):241–50.

    Article  PubMed  CAS  Google Scholar 

  68. Li Y, Yang X, Zhu W, Xu Y, Ma J, He C, Wang F. SWI/SNF complex gene variations are associated with a higher tumor mutational burden and a better response to immune checkpoint inhibitor treatment: a pan-cancer analysis of next-generation sequencing data corresponding to 4591 cases. Cancer Cell Int. 2022;22(1):347.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  69. Tokunaga R, Xiu J, Goldberg RM, Philip PA, Seeber A, Battaglin F, Arai H, Lo JH, Puccini A, Naseem M, et al. Gene mutations of SWI/SNF complex and molecular profile in colorectal cancer. J Clin Oncol. 2019;37(15_suppl):3600–3600.

    Article  Google Scholar 

  70. Tan H, Baca Y, Xiu J, Figueredo J, Florou V, Lou E, Marshall J, El-Deiry WS, Korn WM, Walker P, et al. Impact of somatic SWI/SNF alterations on the genomic landscape of pancreatic ductal adenocarcinoma and response to PARP and immune checkpoint inhibitor therapy. J Clin Oncol. 2022;40(16_suppl):e16289–e16289.

    Article  Google Scholar 

  71. Zhou P, Wu X, Chen H, Hu Y, Zhang H, Wu L, Yang Y, Mao B, Wang H. The mutational pattern of homologous recombination-related (HRR) genes in Chinese colon cancer and its relevance to immunotherapy responses. Aging (Albany NY). 2020;13(2):2365–78.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the patients who volunteered to participate in this study and the staff members at the study site who cared for these patients; the members of the data and safety monitoring committee; and the representatives of the sponsors.

Funding

This work was supported by the National High Level Hospital Clinical Research Funding (2022-PUMCH-B-128), CAMS Innovation Fund for Medical Sciences (CIFMS) (2022-I2M-C&T-A-003, 2021-I2M-1-061, and 2021-I2M-1-003), CAMS Clinical and Translational Medicine Research Funds (2019XK320006), CSCO-hengrui Cancer Research Fund (Y-HR2019-0239, Y-HR2020MS-0415, and Y-HR2020QN-0414), CSCO-MSD Cancer Research Fund (Y-MSDZD2021-0213), the Fundamental Research Funds for the Central Universities (3332023011), the Young Scientists Fund of the National Natural Science Foundation of China (82303720), and National Ten-thousand Talent Program. The funders had no role in the study design, data collection and analysis, decision to publish or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

Conception and design: XY1, BL, NZ, JL, XY2, KW, HZ2; development of methodology: XY1, BL, NZ, JL, ZS, LL, AW, LW, XY2, KW, HZ2; acquisition of data (acquired and managed patients, provided facilities, etc.): XY1, NZ, JL, YL, JX, XC, YW1, YW2, ZX, MP, CZ, SW, HS, JP, XH, MG, LH, JS, JZ, ZL, YM, XS, LW, XY2, KW, HZ2; analysis and interpretation of data: XY1, BL, NZ, JL, ZS, LL, LW, XD, AW, WL, XY2, KW, HZ2; writing, review, and/or revision of the manuscript: XY1, BL, NZ, JL, XY2, KW, HZ; administrative, technical, or material support: XY1, BL, NZ, JL, ZS, LL, HZ1, LW, XY2, KW, HZ2; study supervision: LW, XY2, KW, HZ2.

Corresponding authors

Correspondence to Xiaobo Yang, Kai Wang or Haitao Zhao.

Ethics declarations

Ethics approval and consent to participate

The study protocol (JS-1391) was approved by the Institutional Ethics Review Committee at PUMCH and QDUH (NCT03892577; RWSHBC). The study obtained written informed consent from all participants and was approved by the Institutional Ethics Review Committee at PUMCH and QDUH (NCT03892577; RWSHBC).

Consent for publication

Not applicable.

Competing interests

Kai Wang is an employer of OrigiMed Co., Ltd. Kai Wang, Leilei Lu, Zhijian Song and Baofeng Lian, Xiaowei Dong, Aodi Wang, and Wenjin Liu are employees of OrigiMed Co., Ltd. The other authors report no potential conflicts of interest for this work.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: 

Table S1. Fifteen pathways in cancer. Table S3. Individual Clinicopathological characteristics of 887 patients with cholangiocarcinoma. Table S4. Genomic alterations in 48 patients with MSI-H cholangiocarcinoma. Table S5. Different gene alterations in MSI-H and MSS cancer. Table S6. Different pathway alterations in MSI-H and MSS cancer. 

Additional file 2: 

Table S2. Clinicopathological characteristics of 887 patients with cholangiocarcinoma. Figure S1. Representative cases of hematoxylin-eosin staining and high PD-L1 expression of MSI-H patients with cholangiocarcinoma. (ABEF) Two ICC cases, (CDGH) 2 ECC cases. ECC: extrahepatic cholangiocarcinoma; ICC: intrahepatic cholangiocarcinoma. Figure S2. Proportion of mutated cholangiocarcinoma samples by MSI status: (y-axis) proportion in MSI-H-specific gene; (x-axis) proportion in MSS-specific gene. MSI-H: microsatellite instability-high, MSS: microsatellite stability. Figure S3. The indel ratio of cholangiocarcinoma in MSI-H and MSS cancer. MSI-H: microsatellite instability-high, MSS: microsatellite stability. Figure S4. ORR (A) and DCR (B) of 139 advanced cholangiocarcinoma patients receiving PD-1 inhibitor-based therapy stratified by MSI status. DCR: disease control rate, MSI: microsatellite instability, ORR: Objective Response Rate. Figure S5. Kaplan‒Meier estimates of overall survival and progression-free survival and clinical benefit histogram of patients with advanced cholangiocarcinoma receiving PD-1 inhibitor-based therapy stratified by PD-L1 CPS ≥5 (ABC). PD-L1: programmed death ligand 1; MSI-H: microsatellite instability-high; MSS: microsatellite stability; Mut: mutation; NDB: no durable benefit; OS, overall survival; PFS, progression-free survival; PD-L1: programmed death ligand 1; TMB: tumor mutation burden. Figure S6. The spectrum of top 30 co-occurring or exclusively occurring mutations for MSI-H (A) and MSS (B) cholangiocarcinoma. MSI-H: microsatellite instability-high, MSS: microsatellite stability. Table S7. Baseline clinicopathological characteristics of 139 patients with advanced cholangiocarcinoma who received PD-1 inhibitor-based therapy.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, X., Lian, B., Zhang, N. et al. Genomic characterization and immunotherapy for microsatellite instability-high in cholangiocarcinoma. BMC Med 22, 42 (2024). https://doi.org/10.1186/s12916-024-03257-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-024-03257-7

Keywords