Multimodal data analysis reveals that pancreatobiliary-type ampullary adenocarcinoma resembles pancreatic adenocarcinoma and differs from cholangiocarcinoma
Journal of Translational Medicine volume 20, Article number: 272 (2022)
Ampullary adenocarcinoma (AAC) arises from the ampulla of Vater where the pancreatic duct and bile duct join and empty into the duodenum. It can be classified into intestinal and pancreatobiliary types based on histopathology or immunohistochemistry. However, there are no biomarkers for further classification of pancreatobiliary-type AAC which has important implications for its treatment. We aimed to identify the tumor origin of pancreatobiliary-type AAC by systematically analyzing whole-slide images (WSIs), survival data, and genome sequencing data collected from multiple centers.
This study involved three experiments. First, we extracted quantitative and highly interpretable features from the tumor region in WSIs and constructed a histologic classifier to differentiate between pancreatic adenocarcinoma (PAC) and cholangiocarcinoma. The histologic classifier was then applied to patients with pancreatobiliary-type AAC to infer the tumor origin. Secondly, we compared the overall survival of patients with pancreatobiliary-type AAC stratified by the adjuvant chemotherapy regimens designed for PAC or cholangiocarcinoma. Finally, we compared the mutation landscape of pancreatobiliary-type AAC with those of PAC and cholangiocarcinoma.
The histologic classifier accurately classified PAC and cholangiocarcinoma in both the internal and external validation sets (AUC > 0.99). All pancreatobiliary-type AACs (n = 45) were classified as PAC. The patients with pancreatobiliary-type AAC receiving regimens designed for PAC showed more favorable overall survival than those receiving regimens designed for cholangiocarcinoma in a multivariable Cox regression (hazard ratio = 7.24, 95% confidence interval: 1.28–40.78, P = 0.025). The results of mutation analysis showed that the mutation landscape of AAC was very similar to that of PAC but distinct from that of cholangiocarcinoma.
This multi-center study provides compelling evidence that pancreatobiliary-type AAC resembles PAC instead of cholangiocarcinoma in different aspects, which can guide the treatment selection and clinical trials planning for pancreatobiliary-type AAC.
Ampullary adenocarcinoma (AAC) is a rare malignant neoplasm that arises within the ampullary complex , which could originate from three types of epithelial cells: biliary, pancreatic ductal, or duodenal. To date, the suitable chemotherapy regimens for AAC remain in the early exploration stage. The rarity of the disease and complexity of histology were the two main barriers to the exploration of effective chemotherapy regimens for AAC.
AAC can be histologically dichotomized into intestinal and pancreatobiliary types . In most studies, pancreatobiliary-type AAC is found to have worse prognosis than intestinal-type AAC [3, 4]. Cancers that arise from different cellular origins often exhibit different sensitivities to therapeutics . Thus, the chemotherapy regimens for AAC should be different for intestinal and pancreatobiliary types [6, 7]. The recommended regimens for intestinal-type AAC tend to be similar to those for colorectal cancer . However, it is still unclear whether pancreatobiliary-type AAC should be treated like cholangiocarcinoma [8, 9] or pancreatic adenocarcinoma (PAC) [4, 6], as there are no sensitive and specific immunohistochemical markers to determine the tumor origin [10, 11]. Effective techniques to identify the tumor origin of pancreatobiliary-type AAC are in great demand and can greatly promote the development of treatments for this disease in the future.
With the rapid development of sequencing technologies in recent years, researchers have proposed using genomics, transcriptomics, and proteomics to infer tumor origin [10, 12,13,14,15,16]. However, these technologies are not routinely used in clinical practice. Advances in computational pathology have demonstrated great success in classifying cancer types [17,18,19,20], predicting cancer prognosis [21,22,23], and detecting genetic alterations [24,25,26]. Computational analysis of histopathological images can identify informative and quantitative features that could be too subtle for pathologists to notice.
In this study, we hypothesized that the site of origin of pancreatobiliary-type AAC can be directly inferred from hematoxylin and eosin (H&E) whole-slide images (WSIs). To solve this challenging task, we first developed a WSI-based classification model using the patients with established diagnosis of cholangiocarcinoma and PAC and then used this model to classify patients with pancreatobiliary-type AAC. The classification model was validated using cohorts from multiple sites. In addition, to support the findings of the histopathological analysis, we compared the overall survival of patients with pancreatobiliary-type AAC treated with adjuvant chemotherapies designed for either cholangiocarcinoma or PAC and compared the mutation landscape of pancreatobiliary-type AAC with those of cholangiocarcinoma and PAC.
We performed histopathological analysis, survival analysis, and mutation analysis to investigate whether pancreatobiliary-type AAC is separable. The overview of our study design is shown in Fig. 1. For the histopathological analysis, we constructed a classification model using the H&E WSIs of cholangiocarcinoma and PAC from The Cancer Genome Atlas (TCGA). The model was then validated in the held-out TCGA dataset (internal validation set) and the SYSUCC dataset (external validation set) to evaluate classification performance and applied to the pancreatobiliary-type AACs in the SYSUCC and Zhejiang datasets to infer the tumor origin. For the survival analysis, we investigated whether pancreatobiliary-type AAC patients receiving adjuvant chemotherapy regimens designed for either cholangiocarcinoma or PAC had significantly different overall survival. For the mutation analysis, we explored and compared the mutation landscape of AAC with those of PAC and cholangiocarcinoma, using the genome sequencing data from cBioPortal.
The flowcharts of collecting the WSI datasets and genome sequencing datasets are shown in Additional file 1: Figs. S1, S2, respectively. For the histopathological analysis, three datasets were used. The TCGA dataset consisted of 38 cholangiocarcinoma patients and 38 PAC patients. The SYSUCC dataset consisted of 71 cholangiocarcinoma patients, 70 PAC patients, and 35 pancreatobiliary-type AAC patients from SYSUCC in Guangdong Province, China. The Zhejiang dataset consisted of 10 pancreatobiliary-type AAC patients from Beilun People’s Hospital and Ningbo Yinzhou People’s Hospital in Zhejiang Province, China. The clinicopathological characteristics of the three datasets are summarized in Additional file 1: Tables S1–S3. In these datasets, each patient had one WSI. All WSIs were in SVS format, acquired from formalin-fixed and paraffin-embedded tumor samples after surgical resection, and centrally reviewed by three pathologists (WH, JY, and MJ).
In the TCGA dataset, tissue slides were scanned using an Aperio ScanScope scanner at 40 × magnification (0.25 μm per pixel). Since there are much more PAC cases than cholangiocarcinoma cases in TCGA, we randomly selected 38 PAC cases to match the number of cholangiocarcinoma cases while keeping a similar distribution of tumor stage, tumor grade, and sex. The identifiers of the cholangiocarcinoma and PAC cases used in our TCGA dataset are listed in Additional file 1: Tables S4, S5.
In the SYSUCC and Zhejiang datasets, tissue slides were scanned using a TEKSSQRAY SQS-1000 scanner at 20 × magnification (0.20 μm per pixel). The histopathologic phenotype of AAC was determined based on H&E staining and immunohistochemical staining if necessary. Pancreatobiliary-type AACs are characterized by small solid nest of cells with rounded nuclei surrounded by desmoplastic stroma and forming simple or branching rounded glands. On immunohistochemical staining, they mainly express CK7, CK19, MUC1, MUC5AC, and MUC6. Conversely, Intestinal-type AACs are characterized by tall often pseudostratified columnar epithelium with oval nuclei forming elongated glands. They mainly express CK20, CDX2, SATB2, and MUC2.
For the survival analysis, 26 pancreatobiliary-type AAC patients in the SYSUCC dataset received adjuvant chemotherapy regimens designed for cholangiocarcinoma (n = 7) or PAC (n = 19). The regimens designed for PAC included S-1 based chemotherapy, FOLFIRINOX, gemcitabine plus nab-paclitaxel, or gemcitabine plus capecitabine, while the regimens designed for cholangiocarcinoma included gemcitabine plus oxaliplatin, capecitabine, or gemcitabine plus cisplatin. The clinicopathological characteristics of the patients used for survival analysis are summarized in Additional file 1: Table S6. The overall survival data for these patients were retrieved from clinical records. Overall survival was measured as the time interval between the date of surgery and the date of death or last follow-up. The median and range of follow-up time were 23.3 and 9.8–39.6 months.
For the mutation analysis, we used the genome sequencing datasets for PAC , cholangiocarcinoma , and AAC  from cBioPortal (http://www.cbioportal.org/datasets). After excluding some cases (see Additional file 1: Fig. S2), the resulting 88 pancreatobiliary-type AACs, 290 cholangiocarcinomas, and 176 PACs were used for the mutation analysis. The URLs for these datasets are provided below:
Cholangiocarcinoma at https://cbioportal-datahub.s3.amazonaws.com/chol_icgc_2017.tar.gz
Computational pathology workflow
Figure 2 shows the computational pathology workflow. We constructed the classification model through manual annotation of tumor region, color normalization, and feature extraction. We then processed the WSIs of pancreatobiliary-type AAC using the same workflow and input the extracted features into the classification model to divide cancers into either cholangiocarcinoma or PAC. The TCGA dataset was divided into the training and validation sets according to a rough ratio of 4:1. We trained a linear SVM classifier based on the top K features selected by the ANOVA F-value using the training set. The two hyperparameters involved in feature selection and model selection, i.e., the number of features K and the regularization parameter C in the linear SVM, were determined by a grid search scheme with five-fold cross validation in the training set. Then, the best hyperparameters were used to train the final classification model using the whole training set, and the model performance was assessed using the untouched validation sets. The linear SVM classifier and feature selection were implemented with a popular machine learning package in python, Scikit-learn v0.19.1.
Extraction of histopathological image features
The feature extraction process consists of three steps: manual tumor region annotation, color normalization, and feature extraction. The histopathological image features were extracted from the tumor region manually annotated by a pathologist. It is necessary to limit the feature extraction to the tumor region. Otherwise, the classification model built on the WSIs of cholangiocarcinoma and PAC will learn from the features specific to normal pancreas and liver tissue; the normal tissue does not help to identify the origin of tissue of pancreatobiliary-type AAC. The annotated tumor region was cropped into tiles with a size of 2000 × 2000 pixels without overlap to facilitate subsequent analysis. To overcome the undesirable color variations due to, for example, different slide scanners and staining protocols in pathology labs, we used a structure-preserving color normalization algorithm  to transform the color appearance of image tiles into that of a target image preferred selected by a pathologist.
We extracted a total of 150 highly interpretable image features for each WSI using the feature extraction algorithm  we previously developed. These features quantitatively described the size, staining, shape, and density of cell nuclei. Briefly, for each WSI we first segmented all nuclei and then computed 10 nucleus-level features for each segmented nucleus, including nuclear area (denoted by area); the major axis length, minor axis length, and their ratio (major, minor, and ratio); the mean intensity in R, G, and B channels (rMean, gMean, and bMean); the mean, maximal, and minimal distance to its neighbors (distMean, distMax, and distMin). Finally, we aggregated each type of nucleus-level features into 15 image-level features, including a 10-bin histogram and five distribution statistics (mean, standard deviation, skewness, kurtosis, and entropy). As a result, each image was described by 150 features (10 × 15). Using the nucleus-level feature ratio as an example, the corresponding 15 image-level features were denoted by ratio_bin1 to ratio_bin10, ratio_mean, ratio_std, ratio_skewness, ratio_kurtosis, and ratio_entropy. The histogram features from ratio_bin1 to ratio_bin10 represent the proportions of nuclei with shape varying from round to elongated in the tumor region. More details about the feature extraction pipeline are provided in our previous work .
We used a two-sided Mann–Whitney U-test to compare each of the 150 image features between cholangiocarcinoma and PAC. For multiple testing correction, p values were adjusted by the false discovery rate (FDR) procedure according to Benjamini & Hochberg adjustment . We used area under receiver operating characteristic curve (AUC) to assess whether the linear SVM classifier could separate cholangiocarcinoma and PAC. The AUC was computed with the R package pROC v1.15.3. For the pancreatobiliary-type AAC patients who received adjuvant chemotherapy in the SYSUCC dataset, we stratified them into two groups according to the regimens designed for either cholangiocarcinoma or PAC. The Kaplan–Meier method was used to estimate the overall survival, and the log-rank test was performed to test whether the overall survival is significantly different between two groups. A multivariable Cox regression model was used to evaluate the independent prognostic value of the treatment group over sex and histologic grade. For the survival analysis, we used the R package survival v2.43–3. OncoPlot was used to delineate the mutation landscapes of PAC, AAC, and cholangiocarcinoma using the R package maftools v2.6.05. Fisher’s exact test on each gene was performed to detect the differentially mutated genes between two datasets (AAC vs. PAC and AAC vs. cholangiocarcinoma). FDR adjustment was also performed for multiple testing correction. A p value < 0.05 was considered statistically significant, or when FDR adjustment was necessary, an adjusted p value (i.e., q value) < 0.05 was used.
Quantitative image analysis identifies significantly different histopathological features between cholangiocarcinoma and PAC
In most cases, pathologists can easily differentiate cholangiocarcinoma and PAC by looking for the normal tissue components in a tissue sample. However, it would be very challenging for pathologists to perform this differentiation by just looking at the tumor region, which is the case for the diagnosis of the tumor origin of pancreatobiliary-type ACC because tumor surrounding tissue now does not provide any useful information. To find the true morphological differences in the two kinds of tumors, we compared each of the 150 quantitative image features extracted from the tumor region using the Mann–Whitney U-test. To enhance the statistical power, we combined the cholangiocarcinoma and PAC samples in the TCGA dataset and the SYSUCC dataset to increase the sample size (n = 109 for cholangiocarcinoma and n = 108 for PAC). Of the 150 image features, 108 features were found significantly associated with the cancer type (q value < 0.05). Figure 3 illustrates the fold change and adjusted p value for each of the 108 features. The fold change for each feature was defined as the ratio of the median feature value between the two groups (cholangiocarcinoma/PAC). A fold change less than 1 means that the median for a specific feature in the cholangiocarcinoma group is less than that in the PAC group.
For the nuclear size related histogram features which describe the proportion of nuclei in a tissue sample from small size to large size, area_bin1 and area_bin2 had fold changes less than one, whereas area_bin4 and area_bin7-10 had fold changes greater than 1. This indicates that cholangiocarcinoma tends to have larger nuclei than PAC. For nuclear shape related features, ratio_bin1 and ratio_bin2 had fold changes greater than one, while the features ratio_bin6-10 had fold changes less than one. The first few bins represent the proportions of relatively round nuclei, and the latter bins represent the proportions of relatively elongated nuclei. Therefore, we can infer that PAC tends to have more elongated nuclei, which essentially are stromal cells, compared with cholangiocarcinoma. Similar analysis of the nucleus density related features such as distMean_bin1-10 and distMean_mean suggests that nuclei are more densely distributed in PAC than cholangiocarcinoma. To the best our knowledge, this is the first study comparing quantitative nuclear features in cholangiocarcinoma and PAC.
Model predictions reveal that pancreatobiliary-type AAC resembles PAC
We used the patients with an established diagnosis of cholangiocarcinoma or PAC to build a linear SVM classification model to differentiate between the two cancer types. The model was trained with 80% samples of the TCGA dataset and validated using the held-out samples as the internal validation set and the SYSUCC dataset as the external validation set. Additional file 1: Fig. S3 shows the receiver operating characteristic curves for the model in the two validation sets. We can see that the classification model perfectly distinguishes cholangiocarcinoma and PAC (AUC > 0.99 in both validation sets). The model’s outputs for the cholangiocarcinoma and PAC patients in the TCGA and SYSUCC datasets are shown separately in Fig. 4A (model’s outputs for the training samples in the TCGA dataset were also included). We then applied this powerful model to the pancreatobiliary-type AACs in the SYSUCC and Zhejiang datasets. As shown in Fig. 4A, all pancreatobiliary-type AACs were classified as PAC, with the predicted probability of being cholangiocarcinoma less than 0.5. These results indicate that pancreatobiliary-type AAC histologically resembles PAC and differs from cholangiocarcinoma. Figure 4B shows the H&E pathological images of six patients which correspond to the black squares in Fig. 4A.
Patients with pancreatobiliary-type AAC benefit more from the adjuvant chemotherapy designed for PAC
Now that we have found that pancreatobiliary-type AAC was histologically similar to PAC, we further investigated whether patients with pancreatobiliary-type AAC receiving adjuvant chemotherapy regimens designed for PAC had significantly better prognosis than those treated with regimens designed for cholangiocarcinoma. In the SYSUCC dataset, 19 pancreatobiliary-type AAC patients received the adjuvant chemotherapy designed for PAC while seven received the adjuvant chemotherapy designed for cholangiocarcinoma. Kaplan–Meier survival curves for the two groups are presented in Fig. 5. The one-year and two-year overall survival rates were 94.7% and 72% for the PAC treatment group vs. 80% and 53.3% for the cholangiocarcinoma treatment group. Using Kaplan–Meier estimates, a more favorable overall survival was seen in patients treated with the regimens designed for PAC (log-rank test P = 0.0162, Fig. 5). Moreover, the type of adjuvant chemotherapy was significantly associated with survival in a multivariable Cox regression model adjusted for sex and histologic grade (hazard ratio = 7.24, 95% confidence interval: 1.28–40.78, P = 0.025). The results of multivariable of Cox regression are shown in Table 1. Although the survival difference between treatment groups is significant, the statistical power may be limited due to the relatively small sample size.
Somatic mutation analysis was performed to compare the mutation landscape of pancreatobiliary-type AAC with those of PAC and cholangiocarcinoma. Figure 6A shows the top 50 most frequently mutated genes in pancreatobiliary-type AAC. The mutation rates of this same set of genes in PAC and cholangiocarcinoma are shown in Fig. 6B, C (rows are also rearranged according to mutation prevalence). For a more intuitive comparison between pancreatobiliary-type AAC and the other two cancer types, the paired bar charts of the mutation rate for pancreatobiliary-type AAC vs. PAC and pancreatobiliary-type AAC vs. cholangiocarcinoma are shown in Fig. 6D, E. As we can see, both pancreatobiliary-type AAC and PAC had high frequency of mutations in TP53 and KRAS (more than 60%). We also observed that even though TP53 and KRAS were the most frequently mutated genes in cholangiocarcinoma (Fig. 6C), their mutation rates were only 26% and 20%, respectively, which were much lower than those in pancreatobiliary-type AAC and PAC.
As shown in Fig. 6F, the number of differentially mutated genes between pancreatobiliary-type AAC and PAC was 2, while the number between pancreatobiliary-type AAC and cholangiocarcinoma was 34 (Fisher’s exact test q < 0.05). In addition, we also observed that the characteristically mutated genes in PAC, such as KRAS and TP53, were the top 2 most differentially mutated genes between pancreatobiliary-type AAC and cholangiocarcinoma. In contrast, the mutation frequency of the two genes were not significantly different between pancreatobiliary-type AAC and PAC. Moreover, previous studies showed that the CDKN2A and TTN were also the most commonly mutated genes in PAC [32, 33]. The mutation rates of the two genes were not significantly different between pancreatobiliary-type AAC and PAC, whereas they were significantly different between pancreatobiliary-type AAC and cholangiocarcinoma. Together, these data showed the mutation landscape of pancreatobiliary-type AAC was very similar to that of PAC but distinct from that of cholangiocarcinoma.
The difficulty for exploring the regimens for pancreatobiliary-type AAC lies in the lack of reliable techniques to infer the tumor origin. Due to the similar immune-histologic manifestations between pancreatic cancer and biliary cancer, it is difficult to further classify pancreatobiliary-type AAC. Therefore, whether pancreatobiliary-type AAC should be treated like PAC or cholangiocarcinoma is an outstanding issue. To our knowledge, this is the first study that uses computational pathology methods to classify pancreatobiliary-type AAC based on routinely available H&E tissue slides. The experimental results showed that PAC and cholangiocarcinoma are highly distinguishable and that pancreatobiliary-type AAC resembles PAC and differs from cholangiocarcinoma.
Differential diagnosis between PAC and cholangiocarcinoma has significant implications in patient management such as chemotherapy regimens and prognosis but is quite challenging histologically. Both carcinomas present similar histomorphology with infiltrating ductal architecture, mild to moderate nuclear atypia, and dense desmoplastic reaction [34, 35]. Many immunohistochemical markers have been tested to aid pathologists in distinguishing PAC from cholangiocarcinoma [36, 37]. However, because of the overlapping immunohistochemical profiles, most of them are not sensitive and specific enough to be used in clinical practice. In this study, based on ubiquitously available H&E slides, our computational analysis of WSIs identified distinctly different image features between PAC and cholangiocarcinoma. Most importantly, these features are highly interpretable. For instance, we found that cholangiocarcinoma tends to have larger nuclei than PAC and that PAC tends to have denser fibrous tissue than cholangiocarcinoma. These subtle differences cannot be captured by human eyes and have not been reported before.
Previous studies have showed similar clinical outcomes and genomic profiles between pancreatobiliary-type AAC and PAC. Williams et al.  reported that pancreatobiliary-type AAC had very similar overall survival to PAC but was significantly more aggressive than intestinal-type AAC. The median overall survival was 33.3, 31.4, and 71.7 months for pancreatobiliary-type AAC, PAC, and intestinal-type AAC, respectively. By genomic sequencing, Yachida et al.  showed that the prevalence of driver gene mutations was distinct between pancreatobiliary-type AAC and intestinal-type AAC. Different form previous studies, we, for the first time, systematically compared the mutation landscape of pancreatobiliary-type AAC with those of PAC and cholangiocarcinoma. We found the mutation landscape of pancreatobiliary-type AAC resembled that of PAC but significantly differed from that of cholangiocarcinoma. This indicates that pancreatobiliary-type AAC is likely to originate from the pancreatic duct epithelium instead of the biliary epithelium.
The role of adjuvant therapy for AAC remains controversial [39,40,41]. For example, Bonet et al. showed that adjuvant therapy after curative-intent resection of AAC was not associated with improved long-term survival . On the other hand, a collaborative study by the Johns Hopkins Hospital and Mayo Clinic reported that adjuvant chemoradiation therapy would improve the outcomes of AAC . There are several reasons for this discordance. First, many studies did not consider the pancreatobiliary-type and intestinal-type separately, which have distinct genomic characteristics and thus should be considered separately. In addition, the regimens of adjuvant therapy used for AAC were varied. The most commonly used chemotherapy regimens for AAC were single gemcitabine; however, gemcitabine alone may be too weak to improve patient outcomes. Ecker et al. reported that gemcitabine-based chemotherapy for pancreatobiliary-type AAC tended to be associated with better survival, though the association was not statistically significant . Since our results showed that all pancreatobiliary-type AACs were classified as PAC, it is worthy to explore whether more effective adjuvant chemotherapy regimens for PAC such as S-1  and mFOLFIRINOX  can improve the outcomes of pancreatobiliary-type AAC. Based on our in-house pancreatobiliary-type AAC dataset, we indeed observed that the adjuvant chemotherapy regimens designed for PAC (mainly S-1) was significantly associated with a survival benefit compared with those designed for cholangiocarcinoma (gemcitabine plus oxaliplatin, capecitabine, or gemcitabine plus cisplatin). This supports our argument that pancreatobiliary-type AAC should be considered as and managed like PAC.
This study has several limitations. First, intratumor heterogeneity is unavoidable in studies that take tumor samples and may affect the performance of the histologic classifier. To alleviate the impact of intratumor heterogeneity, we only included the tissue slides obtained from surgical resection specimens which contain a much larger tumor area compared with biopsy specimens. Secondly, due to the rarity of pancreatobiliary-type AAC and only a subset of patients receiving adjuvant therapy, the statistical power may be limited in the survival analysis comparing two treatment groups. Further clinical trials are warranted to validate the effectiveness of applying the adjuvant or first-line chemotherapy regimens for PAC to pancreatobiliary-type AAC.
This multicenter study provides a promising histologic model for the classification of cholangiocarcinoma vs. PAC and for inferring the tumor origin of pancreatobiliary-type AAC, using routinely available H&E stained diagnostic slides without the extra expense of, for example, genome sequencing. The analyses of different types of data suggest that pancreatic ductal epithelial cells are likely to be the site of origin for pancreatobiliary-type AAC, which could guide clinicians and researchers to select treatments and plan clinical trials for this disease in the near future.
Availability of data and materials
Quantitative features extracted from all H&E stained whole-slide images, survival data, processed genomic sequencing data, and source code are available on GitHub (https://github.com/chengjun583/PB-type-ampullary-adenocarcinoma). The sample identifiers of the whole-slide images used in the TCGA dataset and the download links to the original genomic sequencing data are available in patient cohorts section. Whole-slide images in the SYSUCC and Zhejiang datasets are available from the authors upon reasonable request.
Area under the receiver operating characteristic curve
False discovery rate
Hematoxylin and eosin
Receiver operating characteristic
Sun yat-sen university cancer center
The cancer genome atlas
Haeberle L, Riemer J, Esposito I. Molecular Pathology of Carcinomas of the Ampullary/Periampullary Region. In: Pancreatic cancer. Springer, New York, NY; 2016;pp. 1–17.
Ang DC, Shia J, Tang LH, Katabi N, Klimstra DS. The utility of immunohistochemistry in subtyping adenocarcinoma of the ampulla of vater. Am J Surg Pathol. 2014;38:1371–9.
Kim WS, Choi DW, Choi SH, Heo JS, Do You D, Lee HG. Clinical significance of pathologic subtype in curatively resected ampulla of vater cancer. J Surg Oncol. 2012;105:266–72.
Williams JL, Chan CK, Toste PA, Elliott IA, Vasquez CR, Sunjaya DB, et al. Association of histopathologic phenotype of periampullary adenocarcinomas with survival. JAMA Surg. 2017. https://doi.org/10.1001/jamasurg.2016.3466.
Gingras MC, Covington KR, Chang DK, Donehower LA, Gill AJ, Ittmann MM, et al. Ampullary cancers harbor ELF3 tumor suppressor gene mutations and exhibit frequent wnt dysregulation. Cell Rep. 2016;14:907–19 (NIH Public Access).
Schiergens TS, Reu S, Neumann J, Renz BW, Niess H, Boeck S, et al. Histomorphologic and molecular phenotypes predict gemcitabine response and overall survival in adenocarcinoma of the ampulla of Vater. Surg (United States). Surgery. 2015;158:151–61.
Okano K, Oshima M, Suto H, Ando Y, Asano E, Kamada H, et al. Ampullary carcinoma of the duodenum: current clinical issues and genomic overview. Surg Today. 2021;52:189–97.
Malka D, Cervera P, Foulon S, Trarbach T, de la Fouchardière C, Boucher E, et al. Gemcitabine and oxaliplatin with or without cetuximab in advanced biliary-tract cancer (BINGO): A randomised, open-label, non-comparative phase 2 trial. Lancet Oncol. 2014;15:819–28.
Lamarca A, Palmer DH, Wasan HS, Ross PJ, Ma YT, Arora A, et al. Second-line FOLFOX chemotherapy versus active symptom control for advanced biliary tract cancer (ABC-06): a phase 3, open-label, randomised, controlled trial. Lancet Oncol. 2021;22:690–701.
Chakraborty S, Ecker BL, Seier K, Aveson VG, Balachandran VP, Drebin JA, et al. Genome-derived classification signature for ampullary adenocarcinoma to improve clinical cancer care. Clin Cancer Res. 2021;27:5891–9.
Hayata Y, Nakagawa H, Kurosaki S, Kawamura S, Matsushita Y, Hayakawa Y, et al. Axin2+ peribiliary glands in the periampullary region generate biliary epithelial stem cells that give rise to ampullary carcinoma. Gastroenterology. 2021;160:2133-2148.e6.
Hayashi H, Takiguchi Y, Minami H, Akiyoshi K, Segawa Y, Ueda H, et al. Site-specific and targeted therapy based on molecular profiling by next-generation sequencing for cancer of unknown primary site: a nonrandomized phase 2 clinical trial. JAMA Oncol. 2020;6:1931–8.
Zhao Y, Pan Z, Namburi S, Pattison A, Posner A, Balachander S, et al. CUP-AI-Dx: a tool for inferring cancer tissue of origin and molecular subtype using rna gene-expression data and artificial intelligence. EBioMedicine. 2020. https://doi.org/10.1016/j.ebiom.2020.103030.
Grewal JK, Tessier-Cloutier B, Jones M, Gakkhar S, Ma Y, Moore R, et al. Application of a neural network whole transcriptome-based pan-cancer method for diagnosis of primary and metastatic cancers. JAMA Netw Open. 2019;2:e192597.
Jiao W, Atwal G, Polak P, Karlic R, Cuppen E, Al-Shahrour F, et al. A deep learning system accurately classifies primary and metastatic cancers using passenger mutation patterns. Nat Commun. 2020. https://doi.org/10.1038/s41467-019-13825-8.
Takenami T, Maeda S, Karasawa H, Suzuki T, Furukawa T, Morikawa T, et al. Novel biomarkers distinguishing pancreatic head cancer from distal cholangiocarcinoma based on proteomic analysis. BMC Cancer. 2019. https://doi.org/10.1186/s12885-019-5548-x.
Cheng J, Han Z, Mehra R, Shao W, Cheng M, Feng Q, et al. Computational analysis of pathological images enables a better diagnosis of TFE3 Xp11.2 translocation renal cell carcinoma. Nat Commun. 2020;11:1–9.
Coudray N, Ocampo PS, Sakellaropoulos T, Narula N, Snuderl M, Fenyö D, et al. Classification and mutation prediction from non–small cell lung cancer histopathology images using deep learning. Nat Med. 2018;24:1.
Campanella G, Hanna MG, Geneslaw L, Miraflor A, Werneck Krauss Silva V, Busam KJ, et al. Clinical-grade computational pathology using weakly supervised deep learning on whole slide images. Nat Med. 2019. https://doi.org/10.1038/s41591-019-0508-1.
Lu MY, Chen TY, Williamson DFK, Zhao M, Shady M, Lipkova J, et al. AI-based pathology predicts origins for cancers of unknown primary. Nature. 2021;594:106–10.
Yu K-H, Zhang C, Berry GJ, Altman RB, Ré C, Rubin DL, et al. Predicting non-small cell lung cancer prognosis by fully automated microscopic pathology image features. Nat Commun. 2016;7:12474.
Cheng J, Zhang J, Han Y, Wang X, Ye X, Meng Y, et al. Integrative analysis of histopathological images and genomic data predicts clear cell renal cell carcinoma prognosis. Cancer Res Inc. 2017;77:91–100.
Cheng J, Mo X, Wang X, Parwani A, Feng Q, Huang K. Identification of topological features in renal tumor microenvironment associated with patient survival. Bioinformatics. 2018;34:1024–30.
Kather JN, Pearson AT, Halama N, Jäger D, Krause J, Loosen SH, et al. Deep learning can predict microsatellite instability directly from histology in gastrointestinal cancer. Nat Med. 2019. https://doi.org/10.1038/s41591-019-0462-y.
Kather JN, Heij LR, Grabsch HI, Loeffler C, Echle A, Muti HS, et al. Pan-cancer image-based detection of clinically actionable genetic alterations. Nat Cancer. 2020;1:789–99.
Cheng J, Liu Y, Huang W, Hong W, Wang L, Zhan X, et al. Computational image analysis identifies histopathological image features associated with somatic mutations and patient survival in gastric adenocarcinoma. Front Oncol. 2021. https://doi.org/10.3389/fonc.2021.623382.
Sanchez-Vega F, Mina M, Armenia J, Chatila WK, Luna A, La KC, et al. Oncogenic signaling pathways in the cancer genome atlas. Cell. 2018;173:321 (NIH Public Access).
Jusakul A, Cutcutache I, Yong CH, Lim JQ, Huang MN, Padmanabhan N, et al. Whole-genome and epigenomic landscapes of etiologically distinct subtypes of cholangiocarcinoma. Cancer Discov. 2017;7:1116 (NIH Public Access).
Gingras MC, Covington KR, Chang DK, Donehower LA, Gill AJ, Ittmann MM, et al. Ampullary cancers harbor ELF3 tumor suppressor gene mutations and exhibit frequent wnt dysregulation. Cell Rep Cell Rep. 2016;14:907–19.
Vahadane A, Peng T, Sethi A, Albarqouni S, Wang L, Baust M, et al. Structure-preserving color normalization and sparse stain separation for histological images. IEEE Trans Med Imaging. 2016;0062:1–1.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57(1):289–300.
Waddell N, Pajic M, Patch AM, Chang DK, Kassahn KS, Bailey P, et al. Whole genomes redefine the mutational landscape of pancreatic cancer. Nature Nature. 2015;518:495–501.
Wang Z, Zou W, Wang F, Zhang G, Chen K, Hu M, et al. Identification of the immune cell infiltration landscape in pancreatic cancer to assist immunotherapy. Future Oncol. 2021;17:4131–43.
Lim JH, Park CK. Pathology of cholangiocarcinoma. Abdom Imaging. 2004. https://doi.org/10.1007/s00261-004-0187-2.
Farrow B, Berger DH, Rowley D. Tumor-derived pancreatic stellate cells promote pancreatic cancer cell invasion through release of thrombospondin-2. J Surg Res. 2009;156:155–60.
Hooper JE, Morgan TK, Grompe M, Sheppard BC, Troxell ML, Corless CL, et al. The novel monoclonal antibody HPC2 and N-cadherin distinguish pancreatic ductal adenocarcinoma from cholangiocarcinoma. Hum Pathol. 2012;43:1583–9 (WB Saunders).
Lok T, Chen L, Lin F, Wang HL. Immunohistochemical distinction between intrahepatic cholangiocarcinoma and pancreatic ductal adenocarcinoma. Hum Pathol. 2014;45:394–400 (WB Saunders).
Yachida S, Wood LD, Suzuki M, Takai E, Totoki Y, Kato M, et al. Genomic sequencing identifies ELF3 as a driver of ampullary carcinoma. Cancer Cell. 2016;29:229–40.
Ecker BL, Vollmer CM, Behrman SW, Allegrini V, Aversa J, Ball CG, et al. Role of adjuvant multimodality therapy after curative-intent resection of ampullary carcinoma. JAMA Surg. 2019;154:706–14.
Bonet M, Rodrigo A, Vázquez S, Carrizo V, Vilardell F, Mira M. Adjuvant therapy for true ampullary cancer: a systematic review. Clin Transl Oncol Clin Transl Oncol. 2020;22:1407–13.
Jin Z, Hartgers ML, Sanhueza CT, Shubert CR, Alberts SR, Truty MJ, et al. Prognostic factors and benefits of adjuvant therapy after pancreatoduodenectomy for ampullary adenocarcinoma: mayo clinic experience. Eur J Surg Oncol. 2018;44:677–83.
Miura JT, Jayakrishnan TT, Amini A, Johnston FM, Tsai S, Erickson B, et al. Defining the role of adjuvant external beam radiotherapy on resected adenocarcinoma of the ampulla of vater. J Gastrointest Surg. 2014;18:2003–8.
Narang AK, Miller RC, Hsu CC, Bhatia S, Pawlik TM, Laheru D, et al. Evaluation of adjuvant chemoradiation therapy for ampullary adenocarcinoma: the johns hopkins hospital—mayo clinic collaborative study. Radiat Oncol. 2011. https://doi.org/10.1186/1748-717X-6-126.
Uesaka K, Boku N, Fukutomi A, Okamura Y, Konishi M, Matsumoto I, et al. Adjuvant chemotherapy of S-1 versus gemcitabine for resected pancreatic cancer: a phase 3, open-label, randomised, non-inferiority trial (JASPAC 01). Lancet. 2016;388:248–57.
Conroy T, Hammel P, Hebbar M, Ben Abdelghani M, Wei AC, Raoul J-L, et al. FOLFIRINOX or gemcitabine as adjuvant therapy for pancreatic cancer. N Engl J Med. 2018;379:2395–406.
This study was supported by National Natural Science Foundation of China (61901275 and 62171290), Shenzhen University Startup Fund (2019131), Shenzhen-Hong Kong Joint Research Program (SGDX20201103095613036), and SZU Top Ranking Project (86000000210).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of SYSUCC (B2021-276-01), and informed consent was waived.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Flowchart of collecting the WSI datasets. Figure S2 Flowchart of collecting the genome sequencing datasets. Figure S3 Performance of the classification model in classifying cholangiocarcinoma and PAC. Table S1 Clinicopathological characteristics of the TCGA dataset. Table S2 Clinicopathological characteristics of the SYSUCC dataset. Table S3 Clinicopathological characteristics of the Zhejiang dataset. Table S4 Identifiers of the cholangiocarcinoma cases used in the TCGA dataset. Table S5 Identifiers of the PAC cases used in the TCGA dataset. Table S6 Clinicopathological characteristics of the SYSUCC AAC dataset for survival analysis.
About this article
Cite this article
Cheng, J., Mao, Y., Hong, W. et al. Multimodal data analysis reveals that pancreatobiliary-type ampullary adenocarcinoma resembles pancreatic adenocarcinoma and differs from cholangiocarcinoma. J Transl Med 20, 272 (2022). https://doi.org/10.1186/s12967-022-03473-w