Article Text

Original research
Integrative analysis reveals early and distinct genetic and epigenetic changes in intraductal papillary and tubulopapillary cholangiocarcinogenesis
  1. Benjamin Goeppert1,2,
  2. Damian Stichel3,4,5,
  3. Reka Toth6,
  4. Sarah Fritzsche1,
  5. Moritz Anton Loeffler1,
  6. Anna Melissa Schlitter7,
  7. Olaf Neumann1,
  8. Yassen Assenov6,
  9. Monika Nadja Vogel8,
  10. Arianeb Mehrabi2,9,
  11. Katrin Hoffmann2,9,
  12. Bruno Köhler2,10,
  13. Christoph Springfeld2,10,
  14. Dieter Weichenhan6,
  15. Christoph Plass4,6,
  16. Irene Esposito11,
  17. Peter Schirmacher1,2,
  18. Andreas von Deimling3,4,5,
  19. Stephanie Roessler1,2
  1. 1Institute of Pathology, University Hospital Heidelberg, Heidelberg, Germany
  2. 2Liver Cancer Center Heidelberg (LCCH), Heidelberg, Germany
  3. 3Clinical Cooperation Unit Neuropathology, German Cancer Research Center (DKFZ), Heidelberg, Germany
  4. 4German Consortium for Translational Cancer Research (DKTK), Heidelberg, Germany
  5. 5Department of Neuropathology, University Hospital Heidelberg, Heidelberg, Germany
  6. 6Cancer Epigenomics, German Cancer Research Center (DKFZ), Heidelberg, Germany
  7. 7Institute of Pathology, Technische Universitat of Munich, Munich, Germany
  8. 8Diagnostic and Interventional Radiology, Thoraxklinik at University Hospital Heidelberg, Heidelberg, Germany
  9. 9Department of General, Visceral and Transplantation Surgery, University Hospital Heidelberg, Heidelberg, Germany
  10. 10Department of Medical Oncology, National Center of Tumor Diseases, Heidelberg, Germany
  11. 11Institute of Pathology, Heinrich-Heine-Universitat Dusseldorf, Dusseldorf, Germany
  1. Correspondence to Dr Stephanie Roessler, Institute of Pathology, University Hospital Heidelberg, Heidelberg, Germany; stephanie.roessler{at}med.uni-heidelberg.de; Prof. Benjamin Goeppert, Institute of Pathology, University Hospital Heidelberg, Heidelberg, Germany; benjamin.goeppert{at}med.uni-heidelberg.de

Abstract

Objective A detailed understanding of the molecular alterations in different forms of cholangiocarcinogenesis is crucial for a better understanding of cholangiocarcinoma (CCA) and may pave the way to early diagnosis and better treatment options.

Design We analysed a clinicopathologically well-characterised patient cohort (n=54) with high-grade intraductal papillary (IPNB) or tubulopapillary (ITPN) neoplastic precursor lesions of the biliary tract and correlated the results with an independent non-IPNB/ITPN associated CCA cohort (n=294). The triplet sample set of non-neoplastic biliary epithelium, precursor and invasive CCA was analysed by next generation sequencing, DNA copy number and genome-wide methylation profiling.

Results Patients with invasive CCA arising from IPNB/ITPN had better prognosis than patients with CCA not associated with IPNB/ITPN. ITPN was localised mostly intrahepatic, whereas IPNB was mostly of extrahepatic origin. IPNB/ITPN were equally associated with small-duct and large-duct type intrahepatic CCA. IPNB exhibited mutational profiles of extrahepatic CCA, while ITPN had significantly fewer mutations. Most mutations were shared between precursor lesions and corresponding invasive CCA but ROBO2 mutations occurred exclusively in invasive CCA and CTNNB1 mutations were mainly present in precursor lesions. In addition, IPNB and ITPN differed in their DNA methylation profiles and analyses of latent methylation components suggested that IPNB and ITPN may have different cells-of-origin.

Conclusion Integrative analysis revealed that IPNB and ITPN harbour distinct early genetic alterations, IPNB are enriched in mutations typical for extrahepatic CCA, whereas ITPN exhibited few genetic alterations and showed distinct epigenetic profiles. In conclusion, IPNB/ITPN may represent a distinctive, intermediate form of intrahepatic and extrahepatic cholangiocarcinogenesis.

  • cholangiocarcinoma
  • bilary duct carcinoma
  • gene mutation
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Significance of this study

What is already known on this subject?

  • Cholangiocarcinoma (CCA) is a highly aggressive disease with poor prognosis due to late diagnosis. Intraductal papillary (IPNB) and tubulopapillary (ITPN) neoplasms of the bile duct have been described as precursor lesions of CCA.

  • The genetic and epigenetic alterations involved in the stepwise progression from non-neoplastic biliary epithelium via precursor lesion to invasive CCA are poorly understood.

What are the new findings?

  • Patients with invasive CCA arising from IPNB/ITPN had better prognosis than patients with non-IPNB/ITPN-CCA.

  • IPNB and ITPN harbour distinct early genetic alterations.

  • IPNB are enriched in mutations typical for extrahepatic CCA, whereas ITPN exhibited few genetic alterations.

  • IPNB and ITPN showed distinct epigenetic profiles and they may represent a distinctive, intermediate form of intrahepatic and extrahepatic cholangiocarcinogenesis.

  • Certain intrahepatic IPNB and ITPN may represent a to-date not recognised precursor lesion for intrahepatic CCA of the small-duct type.

How might it impact on clinical practice in the foreseeable future?

  • These newly identified alterations corroborate IPNB and ITPN as unique forms of cholangiocarcinogenesis.

  • The presence of IPNB or ITPN may serve as an indicator for better patient survival.

  • The here identified molecular alterations may facilitate early diagnosis and improve treatment options for patients with CCA.

Introduction

Cholangiocarcinoma (CCA) is a highly aggressive disease with poor prognosis due to late diagnosis and chemotherapy resistance.1–3 Based on the anatomic location, CCA is classified as intrahepatic (iCCA), perihilar (pCCA) or distal CCA (dCCA). Irrespective of the localisation, most CCA are unresectable at diagnosis and thus, early detection prior to malignant transformation as non-invasive precursor lesion may dramatically improve patient outcome. Two main subtypes of biliary precursor lesions that precede invasive CCA are nowadays recognised and described in detail in the current, 5th WHO classification of Digestive System Tumours: the microscopic biliary intraepithelial neoplasm (BilIN) and the macroscopically visible intraductal papillary or tubulopapillary neoplasms of the bile duct, IPNB and ITPN, respectively.4 5 While BilIN are frequently encountered in resection specimens of the biliary tract, IPNB and in particular ITPN are less frequent. However, the exact frequencies of IPNB and ITPN are unknown, and misclassification is still an issue. In addition, frequency of intraductal papillary and tubulopapillary biliary precursors seem to depend on geographical and aetiological factors.6–8 Moreover, there might be an underestimation of IPNB/ITPN frequency, because some cases may be misdiagnosed as papillary CCA. A more detailed understanding of cholangiocarcinogenesis may allow standardised screening approaches of high-risk patients and may lead to targeted therapies for unresectable patients. While IPNB has been regarded as the biliary counterpart of intraductal papillary mucinous neoplasm (IPMN) of the pancreas based on their histological similarities,9–11 there is rising evidence that IPMN and IPNB are not identical, in particular regarding morphology, molecular biology and clinical course.12 A recent Taiwanese study of 37 IPNB identified frequent mutations of KRAS (49%), GNAS (32%), RNF43 (24%), APC (24%), TP53 (24%) and CTNNB1 (11%) suggesting that IPNB is a heterogeneous disease involving multiple mutations leading to activation of MAP kinase, Wnt/β-catenin and cAMP signalling.13 However, it appears that mutation frequencies vary between geographic locations as GNAS mutation was a rare event affecting only 2.5% of IPNB in a European study.14 Furthermore, PRKACA and PRKACB gene fusions were exclusively found in intraductal oncocytic papillary neoplasms of the pancreas and bile duct suggesting that oncocytic papillary neoplasms may represent a molecularly distinct subtype.15

DNA methylation signatures have recently gained attention due to their capability of differentiating tumour and tissue types. This resulted in the recent development of promising algorithms to characterise cancers of unknown primary, brain tumours and sarcomas according to their epigenetic signatures.16–18 In addition, integrative genetic and epigenetic analyses may reveal certain subgroups in tumours that are morphologically identical but display distinct molecular patterns including targetable options. Taken together, only few genetic changes have been described in previous studies. Moreover, changes of protein expression of some candidate genes have been described in immunohistochemical studies. However, until now, a systematic and integrative approach in characterising this branch of cholangiocarcinogenesis is lacking. Here, we present a comprehensive and well-characterised single-centre cohort of solely IPNB or ITPN-associated CCA, which we analysed using an integrative histomorphological, genetic and epigenetic approach. We identified genetic differences of IPNB and ITPN, demonstrated that IPNB and ITPN have both similar and distinct mutational profiles compared to the corresponding invasive CCA and correlated these results with clinicopathological data. In addition, we analysed the epigenetic profiles of IPNB and ITPN using DNA methylation profiling indicating different cells-of-origin of IPNB and ITPN.

Materials and methods

Study population and histomorphological subclassification

This study included 54 patients with high-grade intraductal neoplasms of the bile duct and 9 patients with ITPN of the pancreas (ITPN-P; table 1, online supplemental table S1). In detail, 44 patients with IPNB and 10 patients with ITPN were included. IPNB was diagnosed applying the current histomorphological and immunohistochemical criteria provided by WHO5 and ITPN was diagnosed using the criteria of AFIP and keynote publications.19 20 In brief, the definition of ITPN as biliary intraductal neoplasms that are predominantly tubular/trabecular, non-mucinous and negative for mucin 5AC (MUC5AC) was used. Corresponding normal bile duct samples of 41 cases with intraductal neoplasm of the bile duct and of all 9 ITPN-P cases were available. In addition, corresponding invasive CCA tumour tissues of 34 patients with IPNB or ITPN were available. CCA were diagnosed and subgrouped using the current criteria provided by WHO and AFIP.19 For iCCA subgrouping (small-duct versus large-duct) WHO criteria were used (see also online supplemental table S2). All tissue samples were provided by the Tissue Bank of the National Center for Tumor Diseases (NCT, Heidelberg, Germany) in accordance with the regulations of the NCT Tissue Bank.

Supplemental material

Table 1

Clinicopathological characteristics of patients with intraductal neoplasms of the biliary tree and associated invasive CCA (n=54)

For clinical analysis, the study population of patients with intraductal neoplasms was compared with an independent set of patients with CCA without intraductal neoplasms. All patients included underwent bile duct or liver surgery at the University Hospital Heidelberg between 1995 and 2016. In total, this cohort included 402 patients and for 294 of these staging and survival data were available.

Each non-neoplastic, precursor and invasive tumour sample was histologically confirmed by at least two experienced consultant pathologists.

Patients or the public were not involved in the design, or conduct, or reporting, or dissemination plans of our research.

Tissue microarrays and immunohistochemical analysis

A tissue microarray (TMA) was fabricated for all IPNB and ITPN samples. The morphological classification into IPNB and ITPN and histomorphological subtyping was supported by immunohistochemical analyses of MUC1, MUC2, MUC5AC, MUC6 and CDX2. In addition, the differential diagnosis of neuroendocrine neoplasia was excluded by performing immunohistochemistry of synaptophysin and chromogranin A for all cases. Technical details of TMA fabrication and immunohistochemistry are provided in the online supplemental material.

Statistical analysis

Statistical analysis and visualisation were performed using the computing environment R (http://www.R-project.org/) and GraphPad Prism 6. Median survival and corresponding 95% confidence interval (CI) were calculated by Kaplan-Meier survival analysis, and the survival distributions for each category were compared using the log-rank test. All reported p values were two-sided and p<0.05 was considered statistically significant. Enrichment analysis was performed using Fisher’s exact or χ² test as appropriate.

Graphical representations of mutational differences by oncoplots and circos plots were created by the publicly available R-packages ComplexHeatmap, circlize and factoextra (http://www.R-project.org/; https://bioconductor.org/).21 Due to limitations in visualising multidimensional data, only pairwise co-occurrences are represented in the circos plots. Frequencies shown by the oncoplots are on a sample per gene basis and do not consider that some genes may contain more than one mutation in the same tumour sample.

Data availability

The Infinium MethylationEPIC BeadChip data generated in this study are available at the Gene Expression Omnibus (GEO) repository under accession GSE156299. All relevant data are available from the authors on request.

Online supplemental methods including TMA fabrication, diagnostic criteria, immunohistochemistry, genomic profiling and DNA methylation profiling are available online.

Results

Clinicopathological characteristics of the study cohort

This study included 54 patients with high-grade intraductal neoplasms of the bile duct with intrahepatic, perihilar and distal anatomical location (table 1 and figure 1A). DNA sequencing, methylation and copy number alteration (CNA) analyses were performed for these precursor lesions, the corresponding non-neoplastic bile duct and invasive CCA (figure 1A–C). The intraductal neoplasms of the bile duct were further subdivided based on histomorphological and immunohistochemical characterisation into 44 cases of high-grade IPNB and 10 cases of high-grade ITPN (figure 1D). In addition, it was confirmed that neuroendocrine neoplasias have not been included in this study (online supplemental figure S1). Concomitant invasive CCA was overall present in 92.6% (50/54) of patients, that is, in 93.2% (41/44) of IPNB and 90% (9/10) of ITPN, respectively, and follow-up data were available for 42 patients (table 1). Consistent with other European studies, underlying hepatobiliary diseases associated with the development of CCA included unspecific chronic cholecystitis, cholangitis and pancreatitis. In addition, specific underlying hepatobiliary diseases, such as primary sclerosing cholangitis, IgG4-associated cholangitis, viral hepatitis B, hemochromatosis or Caroli syndrome, were detected overall in 16.7% of patients with IPNB/ITPN (online supplemental table S3).

Figure 1

Study design, sample and patient characterisation. (A) A total of 54 patients with IPNB or ITPN neoplasms of the biliary tree were included in this study. The heatmap indicates the available DNA sequencing and DNA methylation/CNA data for IPNB/ITPN precursor lesions, invasive CCA and normal tissue. For 45 precursors lesions all three analyses, DNA sequencing, DNA methylation and CNA data were obtained. (B) For DNA sequencing, precursor lesions of 48 patients and 26 corresponding invasive CCA samples were analysed. (C) Genome-wide DNA methylation and genome-wide CNA profiling was performed for precursor lesions of 51 patients and compared with 41 corresponding non-neoplastic bile duct (normal) and 27 corresponding invasive CCA samples. (D) Immunohistochemical and histomorphological subtyping was performed for the entire cohort; immunohistochemical characterisation of one representative IPNB (top row) and one representative ITPN (bottom row) case is depicted. (E) Kaplan-Meier survival curves of IPNB and ITPN cases. (F) Kaplan-Meier survival curves of IPNB/ITPN-CCA cases with UICC2 or UICC3 compared with non-IPNB/ITPN-CCA cases of an independent cohort with matched numbers of UICC2 or UICC3, respectively. CCA, cholangiocarcinoma; CNA, copy number alteration; IPNB, intraductal papillary neoplasm of the bile duct; ITPN, intraductal tubulopapillary neoplasm of the bile duct; NA, not available.

Survival analysis showed that the patients with IPNB-associated and ITPN-associated CCA of our cohort had high 5-year survival rates of 67.6% and 75%, respectively (figure 1E). To test if the favourable outcome of IPNB/ITPN-associated CCA was due to early detection, we compared the survival of patients with or without IPNB/ITPN-associated CCA and drew on our CCA cohort of 294 patients without intraductal neoplasms of the bile duct (non-IPNB/ITPN-CCA) for which staging and survival data were available. These 294 patients with CCA but without intraductal neoplasms of the bile duct had mainly UICC2 (n=154) or UICC3 (n=81) stage (online supplemental figure S2A). As UICC staging is strongly associated with patient outcome, we compared the UICC stages of patients with IPNB/ITPN and non-IPNB/ITPN-associated CCA. This revealed that UICC staging was associated with patient outcome in CCA without IPNB/ITPN but not in IPNB/ITPN-associated CCA (online supplemental figure S2A,B). In addition, we combined patients into low-stage (UICC 0–2) and high-stage (UICC 3–4) patient groups to increase case numbers in the IPNB/ITPN-associated CCA group. This revealed that UICC staging was significantly associated with overall survival in non-IPNB/ITPN-CCA but not in patients with IPNB/ITPN-CCA (online supplemental figure S2C,D). Last, we compared the overall survival of patients with IPNB/ITPN-CCA with patients without IPNB/ITPN-CCA by combining the largest UICC staging group of patients (UICC 2–3; intermediate group) and by using proportional numbers of UICC 2 and UICC 3 for each cohort. This comparative retrospective analysis showed that patients with IPNB/ITPN-CCA have significantly better outcome than patients without IPNB/ITPN-CCA (log-rank p<0.001; figure 1F).

Next, correlation analyses of clinicopathological data revealed that patients with IPNB were older than patients with ITPN (p=0.038) but no significant differences in sex or tumour staging were detected. Only 18.2% (8/44) of IPNB but 80.0% (8/10) of ITPN had an intrahepatic localisation (Fisher’s exact p<0.001, table 1). For 16 patients with intrahepatic intraductal high-grade precursor lesion, 11 corresponding invasive iCCA were available and subtyping according to the current WHO classification5 showed that 45.5% (5/11) were small-duct and 54.5% (6/11) were large-duct iCCA. Of these 11 cases, 4 cases had IPNB with large-duct iCCA, 2 cases had IPNB with small-duct iCCA, 2 cases had ITPN with large-duct iCCA and 3 cases had ITPN with small-duct iCCA (Fisher’s exact p=0.567, online supplemental figures S3–S5). Histomorphology of all iCCA with corresponding IPNB or ITPN is depicted in online supplemental figure S4 and online supplemental figure S5 for large-duct and small-duct type, respectively. Histologically, IPNB had a high frequency of pancreatobiliary subtype (65.9%, 29/44), followed by intestinal (25%, 11/44) and gastric histology (9.1%, 4/44). Oncocytic differentiation was not observed in our IPNB cohort. Interestingly, only one ITPN was oncocytic (10%, 1/10), while the remaining ITPN (90%, 9/10) were not otherwise specified and intestinal or gastric histology was not seen in our ITPN cohort (Fisher’s exact test, p=0.037; table 1). In addition, IPNB subgrouping according to the Japan-Korea Cooperative Study Group22 showed no significant correlation to other clinicopathological or molecular variables (online supplemental table S4). In conclusion, ITPN was significantly more frequent intrahepatic compared with IPNB and for both precursor types an association with iCCA of the small-duct subtype was seen (online supplemental figure S5).

Mutational profiles of intraductal papillary and tubulopapillary neoplasms of the bile duct

The 54 precursor lesions were analysed by massive parallel next generation sequencing and 48 samples (39 IPNB and 9 ITPN) passed the DNA quality control (figure 1A). In 83.3% (40/48) of precursor lesions, at least one genomic alteration was detected. A total of 137 non-synonymous mutations in 24 genes and 39 CNAs in 13 of the 40 targeted genes were identified. Mutations with the highest frequency were TP53 (45.8%), KRAS (22.9%), SMAD4 (20.8%) and CDKN2A (20.8%; figure 2). The inter-relationship analysis by circos plot revealed that TP53 mutations are linked to KRAS and SMAD4 mutations (figure 2A). Next, the comparison of the 39 IPNB and 9 ITPN revealed that no mutation was detected in 44.4% (4/9) of ITPN but only in 10.3% (4/39) of IPNB (figure 2B). Furthermore, the number of mutations per sample was significantly higher in IPNB with a mean of 2.39 mutations per sample (95% CI 1.91 to 2.86) compared to ITPN with a mean of 0.89 mutations per sample (95% CI 0.08 to 1.70; Mann Whitney U-test, p=0.006).

Figure 2

Massive parallel sequencing reveals frequent genetic alterations in intraductal neoplasms of the bile duct. (A) Circos plot representing co-occurrence of mutations in all 48 non-invasive IPNB and ITPN samples. A band connecting genes represents co-occurring mutations in a given patient. The width of the band represents the frequency of this mutation pair within the dataset. (B) Oncoplot of genomic alterations in IPNB (n=39) and ITPN (n=9) samples. Missense mutations, in-frame mutations, truncations, amplificationsand deletions of 48 patients with intraductal neoplasms are shown. AMP, amplification; DEL, deletion; INFRAME, in-frame mutation; IPNB, intraductal papillary neoplasm of the bile duct; ITPN, intraductal tubulopapillary neoplasm of the bile duct; TRUNC, truncation.

Next, we compared the mutational profiles of our IPNB/ITPN cases with publicly available data of these precursor lesions and found similar mutation frequencies for TP53, SMAD4, FBXW7 and BRAF (online supplemental table S5).13 23–25 Compared to a large cohort of patients with CCA, the here observed mutation frequencies of IPNB and ITPN mirrored the ones of invasive dCCA and iCCA, respectively.26 Consistent with the low frequency of intrahepatic IPNB and high frequency of intrahepatic ITPN, no IDH1 mutation was observed in IPNB but one ITPN case with associated small-duct type iCCA (online supplemental figure S5E) harboured a hotspot IDH1:p.Arg132Leu mutation which typically occurs in iCCA of the small-duct type (figure 2B, online supplemental table S6). In addition, one intrahepatic ITPN and its associated small-duct iCCA (online supplemental figure S5C) both contained a FGFR2:p.Tyr286Cys mutation and one intrahepatic IPNB and the associated small-duct iCCA showed both the identical ERBB2:p.Ala771_Met774dup mutation (online supplemental figure S5A).

One CTNNB1 amplification and six missense mutations were identified in IPNB exclusively (figure 2B, online supplemental table S6). Thus, ITPN exhibit distinct mutation profiles and harbour overall significantly fewer mutations than IPNB suggesting distinct oncogenic mechanisms in ITPN-derived cholangiocarcinogenesis.

Invasive CCA evolve from IPNB or ITPN and distinct mutational changes occur during tumour evolution

The mutational profiles of corresponding invasive CCA were available for 26 cases. Overall, most mutations were shared between the corresponding precursor lesion and invasive CCA samples (figure 3). Mutation frequencies of TP53 (42.3% vs 38.5%), KRAS (30.8% vs 34.6%) and SMAD4 (19.2% vs 15.4%) were almost identical in precursor lesions and invasive CCA, whereas CTNNB1 (19.2% vs 3.8%), CDKN2A (23.1% vs 11.5%), ROBO2 (0% vs 19.2%) and FBXW7 (7.7% vs 19.2%) differed between precursor lesions and corresponding invasive CCA (online supplemental table S7). Of the 60 distinct mutations observed in invasive CCA and the 57 distinct mutations observed in the precursor lesions, 36 mutations were common for both (figure 3C), supporting the hypothesis that invasive CCA evolved from the corresponding precursor lesion. Furthermore, we performed case-by-case comparison of each alteration distinguishing mutations observed in the precursor lesion only, in the corresponding invasive CCA only or shared by the precursor lesion and corresponding invasive CCA. Among the genes with mutations shared between precursor lesion and invasive CCA were TP53 and KRAS (figure 3D). Mutations in several genes decreased during cholangiocarcinogenesis, including CTNNB1, CDKN2A and SMAD4 which were predominantly observed in the precursor lesion only. In contrast, alterations of other genes increased in CCA development via IPNB or ITPN, in particular ROBO2 and FBXW7 (figure 3D). Thus, CTNNB1 and CDKN2A may be involved in the development of IPNB/ITPN, whereas ROBO2, FBXW7 and CDKN2B may drive late cholangiocarcinogenesis and transformation to invasive CCA.

Figure 3

IPNB/ITPN and corresponding invasive CCA tissue samples share genetic molecular profiles. (A) Oncoplot of genomic alterations in IPNB/ITPN (n=26) and (B) in the corresponding invasive CCA samples (n=26). Missense mutations, in-frame mutations, truncations, amplifications and deletions are shown. Genes (rows) of both oncoplots (A,B) were sorted by the total number of mutations in all samples. In addition, the order of samples (columns) in both oncoplots is based on the number of alterations in the precursor samples; thus, invasive CCA are sorted identical to the corresponding precursor lesions. (C) Overlap of specific mutations, that is, mutation of a certain amino acid in the respective gene, present in IPNB/ITPN precursor lesions and corresponding invasive CCA. (D) Shown are the proportions of mutations observed exclusively in the precursor lesions, exclusively in the invasive CCA or shared by the precursor lesion and corresponding invasive CCA in the indicated genes. (E) Diagram depicting the mutational alterations during progression from non-neoplastic normal bile to precursor and invasive CCA. AMP, amplification; DEL, deletion; INFRAME, in-frame mutation; IPNB, intraductal papillary neoplasm of the bile duct; ITPN, intraductal tubulopapillary neoplasm of the bile duct; ITPN, intraductal tubulopapillary neoplasm of the bile duct; TRUNC, truncation.

Genome-wide copy number alterations

Genomic chromosomal imbalances have been observed in CCA by multiple studies.27–30 Using DNA methylation microarrays, we analysed the genomic profiles of normal bile duct (n=41), IPNB (n=41), ITPN (n=10) and the associated invasive CCA (iCCA: n=9; pCCA: n=8; dCCA: n=10) obtained by 850k DNA methylation analyses (figures 1 and 4). As expected, no significant CNA were observed in the normal bile duct controls (figure 4A). IPNB exhibited frequent deletions on chromosome 9p, 17p and 18q, whereas ITPN showed a high frequency of 1q gain and 6q loss (figure 4B,C). Interestingly, these genomic imbalances reflected the clinicopathologically observed differences in the anatomic localisation of IPNB/ITPN. Consistent with the notion that 80.0% of ITPN were associated with iCCA or had intrahepatic localisation, invasive iCCA genomes exhibited a high frequency of 1q amplification and 6q deletion (figure 4D). However, deletion of chromosome 18q, which is frequently detected in IPNB but not in ITPN, was observed in pCCA but not in iCCA (figure 4D,E). In addition, deletions of 9p and 17p, frequently detected in IPNB, were also frequent in the corresponding invasive dCCA tumour samples (figure 4F). Thus, IPNB/ITPN exhibit frequent CNA which appear to be stable during the transition from precursor to invasive carcinoma and specific to the anatomical location of the invasive carcinoma.

Figure 4

Copy number alteration profiles of normal, IPNB or ITPN and invasive CCA samples. (A) Copy number alterations of normal bile duct controls, (B) IPNB, (C) ITPN, (D) iCCA, (E) pCCA and (F) dCCA samples are shown along the chromosomes on the x-axes. The relative frequency of observed gains (green) and losses (red) are depicted above and below the horizontal line, respectively. CCA, cholangiocarcinoma; dCCA, distal CCA; iCCA, intrahepatic CCA; IPNB, intraductal papillary neoplasm of the bile duct; ITPN, intraductal tubulopapillary neoplasm of the bile duct; pCCA, perihilar CCA.

DNA methylation profiles reveal distinct epigenetic profiles in IPNB and ITPN

Epigenetic differences between IPNB and ITPN were detected in terms of divergent DNA methylation profiles. In detail, we compared DNA methylation patterns of non-neoplastic (normal) bile duct epithelia to corresponding biliary precursor lesions and invasive CCA. In addition, we included ITPN of the pancreas (ITPN-P) for comparison in the DNA methylation analysis (online supplemental table S1). Unsupervised clustering of all samples (n=137) resulted in two main clusters which included the non-neoplastic normal bile duct samples and the precursors with corresponding invasive CCA samples, respectively (online supplemental figure S6A). To directly compare the triplet samples of normal bile duct, precursor and invasive CCA, we restricted the clustering analysis to patients for whom all three samples were available (n=24). This resulted in separation into three clusters containing mainly the normal bile duct samples (cluster II), ITPN, iCCA and pCCA samples (cluster I) or IPNB and dCCA samples suggesting that the main difference could be observed between normal and precursor or invasive CCA samples (figure 5A). Furthermore, we observed that 14/27 precursor lesions and corresponding invasive CCA clustered directly next to each other proposing that paired precursor lesions and corresponding invasive CCA exhibit similar DNA methylation profiles (online supplemental figure S6B). Visualisation of DNA methylation patterns by two-dimensional t-SNE plot showed that non-neoplastic bile duct samples separated well from IPNB/ITPN samples (figure 5B). IPNB and ITPN did not appear as distinct groups in the t-SNE plot but subtle differences were detectable (figure 5B). To assess the cell type heterogeneity and to trace the cell‐of‐origin in the precursor samples, we used MeDeCom to decompose methylation data into latent methylation components (LMCs) and analysed the proportions of the resulting LMCs in each sample.31 Following this strategy, we found that IPNB, ITPN and ITPN-P samples separated into two main clusters of which one cluster contained all ITPN and 5/41 IPNB cases, whereas the second cluster contained the remaining 36 IPNB cases and no ITPN case (Fisher’s exact test p<0.001; figure 5C). Interestingly, ITPN-P cases were distributed equally in both clusters suggesting that the cells-of-origin of IPNB, ITPN and ITPN-P differ. Furthermore, it appeared that ITPN had high proportions of LMC1, whereas IPNB samples separated further into three groups with high LMC2, LMC3 or LMC4, respectively. To further dissect the similarities between precursor lesions and corresponding invasive CCA, we performed MeDeCom analysis of the 27 paired precursor and invasive CCA samples and found that all ITPN fell into the same cluster (figure 5D). Thus, all ITPN may have the same cell-of-origin but IPNB may arise from at least three different cell-of-origin types.

Figure 5

DNA methylation profiles of IPNB and ITPN precursor lesions. (A) Unsupervised hierarchical clustering analysis of the triplet normal bile duct, precursor and corresponding invasive CCA of 24 patients. The total number of samples is 72. (B) t-SNE plot of the DNA methylation profiles of bile duct controls of the corresponding IPNB, ITPN and ITPN-P cases (n=50), IPNB (n=41), ITPN (n=10) and ITPN-P (n=9). (C) MeDeCom analysis of IPNB (n=41), ITPN (n=10) and ITPN-P (n=9) samples and (D) of 27 IPNB and ITPN samples with corresponding invasive CCA. CCA, cholangiocarcinoma; IPNB, intraductal papillary neoplasm of the bile duct; ITPN, intraductal tubulopapillary neoplasm of the bile duct; ITPN-P, ITPN of the pancreas.

Discussion

In this study, we characterised a comprehensive cohort of patients with IPNB or ITPN precursor lesions and associated invasive CCA and analysed the genetic and epigenetic alterations. By this integrative analysis and by using intraindividual triplet specimens of non-neoplastic tissue, precursor lesions and invasive CCA, we could delineate distinct molecular patterns in IPNB/ITPN-associated cholangiocarcinogenesis over time and correlate the findings with morphological subtyping and clinicopathological data. Consistent with previous studies of a cohort from Thailand and a mixed American, Asian and European cohort, we demonstrated that the simultaneous presence of IPNB or ITPN precursor lesions is a significant favourable prognostic factor regarding patient survival.25 32

We found that genome-wide DNA methylation profiles of ITPN differ from IPNB profiles, suggesting different epigenetic mechanisms in ITPN and IPNB. In addition, we performed MeDeCom analysis to dissect methylation patterns into LMCs, thereby gaining information on possible cell‐type composition differences. In conclusion, we could reveal (1) that overall DNA methylation patterns between IPNB and ITPN show only modest separation, (2) that IPNB and ITPN may have different cells-of-origin and (3) that IPNB show a higher inter-tumour heterogeneity compared with ITPN.

Panel sequencing data of our cohort confirmed well-known candidate genes that were already discovered in CCA before, such as TP53, KRAS, SMAD4, CDKN2A and ERBB2, to also play a pivotal role in cholangiocarcinogenesis via IPNB.26 33–35 Furthermore, we could deduce a genetic evolution plot presenting distinct mutations that emerge in IPNB/ITPN and vanish during transformation to invasive CCA. In particular, CTNNB1 and CDKN2A mutations occurred in a significant number of the precursor lesions and were then mostly lost in invasive CCA. In contrast, ROBO1, ROBO2 and FBXW7 gene alterations were increasingly detected in late cholangiocarcinogenesis, that is, invasive CCA, and were mostly absent in IPNB/ITPN. Noteworthy, the tumour suppressor gene FBXW7 was previously identified as an epigenetically regulated gene in CCA, and it could be shown that FBXW7 is a direct target for downregulation by miR-200a and miR-429, respectively.36 In addition to epigenetic silencing, genomic loss through chromosome 4q deletion may be an alternative mechanism for shutting down FBXW7.30 The tumour suppressive function of the E3-ubiquitin ligase FBXW7 has been suggested to depend on MYC activation, as a dominant negative form of MYC efficiently inhibited AKT/FBXW7-driven tumourigenesis in a CCA mouse model.37 Furthermore, ROBO1 and ROBO2 mutations were exclusively detected in invasive carcinoma. Recurrent ROBO receptor mutations have been also identified by other studies on CCA but ROBO receptor function has been mainly studied in neuronal development implicating a critical repulsive cue for axon path finding on binding of the Slit ligands.38–40 However, SLIT-ROBO signalling has more recently also been implicated in cancer cell migration, invasion and tumour angiogenesis.41–45 Thereby, ROBO proteins may either inhibit or promote tumourigenesis depending on the tumour entity and subtype. Our data suggest that ROBO1/2 and FBXW7 gene function may play a role in late tumour development and tumour progression of CCA.

Interestingly, all four IPNB cases with CTNNB1 mutation lost the CTNNB1 mutation during cholangiocarcinogenesis, which indicates that either CTNNB1 gene function is not necessary in late cholangiocarcinogenesis or other components of the ß-catenin signalling pathway are functionally crucial at this stage of CCA development. In line with our observation, nuclear ß-catenin accumulation was previously detected in IPNB lesions but not in invasive CCA of the same patients.14 46 Consistently, loss of IPNB-derived CTNNB1 and APC mutations during full malignant transformation has been observed by sequencing approaches.14 46 In contrast to IPNB lesions, neither ITPN nor CCA specimens of our cohort exhibited CTNNB1 mutations. Moreover, we have found in a previous study that epigenetic mechanisms, in particular DNA methylation patterns, play a pivotal role in Wnt pathway regulation of non-IPNB/ITPN associated pCCA and iCCA.47 Earlier, it could be shown that altered expression of CTNNB1 in iCCA may not be caused by genetic mutation.48 Our here presented data confirm these results and expand this concept to IPNB-derived cholangiocarcinogenesis. Similar to CTNNB1, CDKN2A alterations may be involved in early IPNB development but may not be required in invasive CCA. Thus, CTNNB1 may be involved in the development of IPNB and CDKN2A may be involved in the development of IPNB and ITPN, but both are possibly not needed for further progression to invasive CCA, whereas ROBO2, FBXW7 and CDKN2B may exert functional roles in late cholangiocarcinogenesis.

Although our sequencing results are overall consistent with previous studies, we did not find any GNAS mutation in our cohort (online supplemental table S5).13 49 Whether this finding reflects geographic, ethnic or etiological differences cannot be solved here. ITPN harboured significantly fewer mutations than IPNB suggesting distinct oncogenic mechanisms specific to each of these distinct precursor subtypes. This is in line with the observation that ITPN were more often localised intrahepatic in our cohort. Noteworthy, one intrahepatic ITPN was the precursor of a small-duct iCCA with IDH1 mutation (online supplemental table S5 and figure S5E), one intrahepatic ITPN and the associated small-duct iCCA (online supplemental figure S5C) both had the identical FGFR2 mutation and one intrahepatic IPNB and the associated small-duct iCCA showed both the identical ERBB2 mutation (online supplemental figure S5A). Interestingly, in one ITPN of oncocytic histomorphological subtype, no mutation was found in our sequencing approach, suggesting epigenetic instead of genetic alterations or alternative molecular pathways particularly involved in this form of cholangiocarcinogenesis. Indeed, PRKACA and PRKACB gene fusions have been found, for example, in pancreatic and biliary oncocytic papillary neoplasms and may be involved in ITPN development as well.15

According to the current WHO classification, iCCA is subdivided into the small-duct and large-duct type. Since large-duct iCCA are frequently accompanied by high-grade precursor lesions, most often BilIN, we were expecting that IPNB/ITPN would be more frequent in large-duct iCCA.5 Surprisingly, we detected a high proportion of small-duct iCCA in IPNB (33%, 2/6) and ITPN (60%, 3/5) cases. This finding is of particular importance as precursor lesions could not be identified for small-duct iCCA to date.

Molecular data on ITPN of the bile duct are scarce. This is the first comprehensive study identifying distinct genomic and epigenomic alterations in a Western-world cohort including ITPN precursor lesions. In conclusion, our integrative morphomolecular analyses define the molecular landscape of IPNB/ITPN-associated CCA development and propose a putative precursor lesion for small-duct iCCA.

Acknowledgments

We are grateful to Veronika Geißler, Nina Wilhelm and Carolin Kerber (Tissue Bank of the National Center for Tumor Diseases (NCT) Heidelberg, Germany and Institute of Pathology Heidelberg) and Angelika Fraas (Institute of Pathology Heidelberg) for excellent technical assistance. Tissue samples were provided by the tissue bank of the National Center for Tumor Diseases (NCT; Heidelberg, Germany) in accordance with the regulations of the tissue bank and with approval of the Ethics Committee of Heidelberg University (S-206/2005; S-207/2005; S-519/2019).

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Twitter @IEspositoPATH

  • Contributors BG initiated the study. BG and SR designed the research project. BG, AM, KH and AMS provided tissue samples and clinicopathological data. BG, DS, SF, ON and SR performed research experiments. BG, DS, SF, MAL, RT, YA, MNV, ON, DW, CP, BK, CS, IE, AVD and SR analysed data or tissue samples. BG and SR wrote the manuscript. All authors read, commented on and approved the manuscript.

  • Funding This work was supported by Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—Project-ID 314905040—SFB/TR 209 Liver Cancer (B01 to SR and Z01, INF to PS) and the European Union’s Horizon 2020 research and innovation programme under grant agreement No 667273 (HEP-CAR) to PS and SR. This study was in part supported by funds from German Cancer Aid (Deutsche Krebshilfe, project no. 70113922) to SR and from the Helmholtz Foundation to CP. We thank the Genomics and Proteomics Core Facility (GPCF) of DKFZ for performing the DNA methylation profiling.

  • Competing interests PS: Grant, boards and presentations from Novartis, and boards from Incyte. CP: Advisory board honoraria from BioMedX. BG: Advisory board from Novartis.

  • Patient consent for publication Not required.

  • Ethics approval The study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki and was approved by the ethics committee of the University of Heidelberg (S-206/2005, S-207/2005 and S-519/2019).

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data availability statement The Infinium MethylationEPIC BeadChip data generated in this study are available at the Gene Expression Omnibus (GEO) repository under accession GSE156299 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE156299). Data are publicly available without any restrictions.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.