Theranostics 2020; 10(21):9779-9788. doi:10.7150/thno.49671

Research Paper

Prediction of clinically relevant Pancreatico-enteric Anastomotic Fistulas after Pancreatoduodenectomy using deep learning of Preoperative Computed Tomography

Wei Mu1,7, Chang Liu2, Feng Gao3, Yafei Qi4, Hong Lu5, Zaiyi Liu6, Xianyi Zhang2, Xiaoli Cai2, Ruo Yun Ji2, Yang Hou2, Jie Tian1,7, Yu Shi2 Corresponding address

1. Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, School of Medicine, Beihang University, Beijing, China, 100191.
2. Department of Radiology, Shengjing Hospital of China Medical University, Shenyang, China.
3. Department of Pancreato-thyroidic Surgery, Shengjing Hospital of China Medical University, Shenyang, China.
4. Department of Pathology, Shengjing Hospital of China Medical University, Shenyang, China.
5. Department of Radiology, Tianjin Medical University Cancer Institute and Hospital, National Clinical Research Center of Cancer, Key Laboratory of Cancer Prevention and Therapy, Tianjin, PR China.
6. Department of Radiology, Guangdong General Hospital, Guangdong Academy of Medical Sciences, 106 Zhongshan Er Road, Guangzhou 510080, China.
7. Key Laboratory of Molecular Imaging, Chinese Academy of Sciences, Beijing, China, 100190.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Citation:
Mu W, Liu C, Gao F, Qi Y, Lu H, Liu Z, Zhang X, Cai X, Ji RY, Hou Y, Tian J, Shi Y. Prediction of clinically relevant Pancreatico-enteric Anastomotic Fistulas after Pancreatoduodenectomy using deep learning of Preoperative Computed Tomography. Theranostics 2020; 10(21):9779-9788. doi:10.7150/thno.49671. Available from http://www.thno.org/v10p9779.htm

File import instruction

Abstract

Rationale: Clinically relevant postoperative pancreatic fistula (CR-POPF) is among the most formidable complications after pancreatoduodenectomy (PD), heightening morbidity/mortality rates. Fistula Risk Score (FRS) is a well-developed predictor, but it is an intraoperative predictor and quantifies >50% patients as intermediate risk. Therefore, an accurate and easy-to-use preoperative index is desired. Herein, we test the hypothesis that quantitative analysis of contrast-enhanced computed tomography (CE-CT) with deep learning could predict CR-POPFs.

Methods: A group of 513 patients underwent pancreatico-enteric anastomosis after PD at three institutions between 2006 and 2019 was retrospectively collected, and formed a training (70%) and a validation dataset (30%) randomly. A convolutional neural network was trained and generated a deep-learning score (DLS) to identify the patients with higher risk of CR-POPF preoperatively using CE-CT images, which was further externally tested in a prospective cohort collected from August 2018 to June 2019 at the fourth institution. The biological underpinnings of DLS were assessed using histomorphological data by multivariate linear regression analysis.

Results: CR-POPFs developed in 95 patients (16.3%) in total. Compared to FRS, the DLS offered significantly greater predictability in training (AUC:0.85 [95% CI, 0.80-0.90] vs. 0.78 [95% CI, 0.72-0.84]; P = 0.03), validation (0.81 [95% CI, 0.72-0.89] vs. 0.76 [95% CI, 0.66-0.84], P = 0.05) and test (0.89 [95% CI, 0.79-0.96] vs. 0.73 [95% CI, 0.61-0.83], P < 0.001) cohorts. Especially in the challenging patients of intermediate risk (FRS: 3-6), the DLS showed significantly higher accuracy (training: 79.9% vs. 61.5% [P = 0.005]; validation: 70.3% vs. 56.3% [P = 0.04]; test: 92.1% vs. 65.8% [P = 0.013]). Additionally, DLS was independently associated with pancreatic fibrosis (coefficients: -0.167), main pancreatic duct (coefficients: -0.445) and remnant volume (coefficients: 0.138) in multivariate linear regression analysis (r2 = 0.512, P < 0.001). The user satisfaction score in the test cohort was 4 out of 5.

Conclusions: Preoperative CT based deep-learning model provides a promising novel method for predicting CR-POPF occurrences after PD, especially at intermediate FRS risk level. This has a potential to be integrated into radiologic reporting system or incorporated into surgical planning software to accommodate the preferences of surgeons to optimize preoperative strategies, intraoperative decision-making, and even postoperative care.

Keywords: Pancreatic fistula, fistula risk score, pancreatoduodenectomy, computed tomography (CT), deep learning

Introduction

Postoperative pancreatic fistula (POPF) is among the most formidable complications after pancreatoduodenectomy (PD), threatening to prolong hospitalization, increasing medical treatment costs, imposing catastrophic abscess or hemorrhage, and heightening morbidity/mortality rates [1-8]. Surgical and technologic advancements have dramatically reduced procedure-related mortality, but rates of clinically relevant postoperative pancreatic fistula (CR-POPF) have remained unchanged (~11-15%) [3]. Managing CR-POPF is thus a high-priority issue for pancreatic surgeons, who need sensitive predictors of CR-POPF to identify high-risk patients and adjust intra- or postoperative care accordingly.

The reported risk factors for CR-POPFs are broadly classifiable as local factors at pancreatic remnant [9], systemic factors (e.g., high body mass index [BMI]), and operative factors (eg, blood loss) [10-13]. Local risk factors at pancreatic remnant, which strongly linked to underlying local histopathologic changes, such as rich viable gland, absence of fibrosis [14], and fatty pancreas [15, 16], represent the most likely determinants directly related to anatomic failure. The Fistula Risk Score (FRS), which incorporated four of the aforementioned parameters: small-sized main pancreatic duct (MPD), soft glandular texture (by surgeon's palpation), high-risk pathology (chronic pancreatitis [CP] or pancreatic adenocarcinoma [PDAC]), and undue intraoperative blood loss, is a well-developed and validated 10-point scale used to intraoperatively predict CR-POPF development after PD [17-19]. Despite the simplicity and convenience of the FRS, it relies on subjective intraoperative findings of surgeons. Moreover, >50% patients qualify as intermediate risk (FRS scores of 3-6), which is a grey zone warranting more objective and reliable predictors.

Beyond FRS, local risk factors via quantitative imaging from standard-of-care contrast-enhanced computed tomography (CE-CT) images, such as pancreatic thickness [20] or remnant volume [21, 22], have been measured manually in dozens of studies conducted during the last decade, showing promising results in predicting CR-POPF events. Recently, the artificial intelligence (AI) based medical image analysis provides an objective and automatic way to capture all important local properties due to its self-learning characteristics, and has achieved success in different fields [23-26]. However, deep-learning model hasn't been investigated in CR-POPF prediction.

Thus, the purpose of this study was to (1) develop and validate an easy-to-use deep-learning model to predict CR-POPF after PD, and (2) determine its diagnostic performance and compare it with FRS, and finally (3) investigate the histomorphologic changes pertaining to deep-learning score (DLS) generated by the deep-learning model.

Methods

Study population

This was a diagnostic, multicenter, multi-cohort study involving four cohorts from four high-volume academic institutions. This study was approved by the Institutional Review Board of each local institution and adhered to ethical standards of the 1964 Helsinki Declaration, including its later amendments. Informed consent was waived in all cohorts. From institutions A-C, 513 patients were obtained retrospectively from radiologic and pathologic archives, and were randomly divided into a training (n = 359) and validation cohort (n = 154) randomly with a ratio of 70/30 to train and validate a deep-learning model based on preoperative CE-CT images. Another cohort (Institution D) was prospectively collected as the external test (n = 70) dataset, all summarized in Figure 1. Notably, for this prospective cohort, patients who were eligible for pancreaticoduodenectomy were consecutively collected during the period. The radiologists who interpreted the radiological images were blinded to the clinical and laboratory results. The surgeons were also blinded to the estimated POPF risk by the deep-learning model. STROBE guidelines [27] for reporting observational studies were applied during study design, training, validation, and reporting of the prediction model.

Clinical data collection and mitigation strategies

Five lead pancreatic surgeons (>20 years of collective pancreatic surgical experience) performed all PDs in conjunction with either pancreaticojejunostomy (PJ) or pancreaticogastrostomy (PG) [28] for a full array of indications. Medical records provided demographic and clinical data, including age, sex, BMI, diabetes mellitus, reported weight loss in the previous 6 months, jaundice, smoking, or alcohol abuse. POPFs were graded in accordance with ISGPF standards [29, 30] as either biochemical or clinically relevant (see Supplemental Table S1).

The four risk factors required for FRS calculations (see Supplemental Table S2) were obtained from operative notes retrospectively or prospectively recorded during surgery by attending surgeons attuned to this study. These risk factors served to generate FRS scores (0-10) individually and thus categorized fistula risk as low (0-2 points), intermediate (3-6 points), or high (7-10 points). Other details are included as Supplemental material S1.

 Figure 1 

Schematic of study design: Data from Institutes A, B and C pertained to clinical information and corresponding imaging details, serving to train and validate the deep learning signature (DLS). Clinical and imaging data contributed by Institute D was used for external test.

Theranostics Image (Click on the image to enlarge.)

For fistula mitigation, less than 1% of patients received biologic sealants; two laminar intraperitoneal drains were routinely placed (100%), and 52 (10%) patients retained an additional retroperitoneal drain. Postoperative somatostatin analogs (eg, prophylactic octreotide) were administered to 225 (43.9%) patients, and trans-anastomotic stents (largely short, internal silicone elastomer tubes) were deployed in 281 (54.8%) patients.

Histology of pancreatic stump

Specimens of the pancreatic stump were evaluated to quantify fibrous tissue, exocrine glandular atrophy (A) [31] and degrees of lipomatosis (L) [32, 33], as detailed Supplemental material S2.

Preoperative CT imaging and segmentation

CT scanning parameters and detailed descriptions are presented in Supplemental Table S3. All patients underwent preoperative multiphasic scans (64-channel multi-detector CT or better) within 4 weeks of surgery at <3-mm minimum slice thickness and in three standard-of-care phases, adhering to current National Comprehensive Cancer Network (NCCN) guidelines [34, 35]. A nonionic contrast agent containing iodine (300mg/mL) was injected at 2-2.5 mL/kg body weight. Median scan delays from injection of contrast to starts of pancreatic parenchymal and portal venous phases were 40-50sec and 65-70sec, respectively. In most patients, the estimated transection line was at superior mesenteric vein, with modifications as needed for individual tumor locations and projected safety margin restrictions.

Volumetric regions of interest (ROIs) in CT images (Figure 2) were segmented separately by four experienced abdominal radiologists (all with over 10 years of experience in pancreatic imaging), using open-source software (3D Slicer version 4.10; www.slicer.org). Segmentations were undertaken in pancreatic phase of transverse sections. In addition, pancreatic thickness, width, and remnant pancreatic volume were measured as previously published fixed CT classifiers (Supplemental material S1).

Deep-learning (DL) model

The pipeline for the CR-POPF deep-learning model is depicted in Supplemental Figure S1. 2D-ROIs were input to the deep-learning model and DLS values were yielded as probabilities of CR-POPF. To further ensure a robust predictive exercise, the DL model received all pancreatic slices in each patient, conveying the average probability of CR-POPF. Further deep-learning training details can be found in Supplemental material S3. Intermediate activation layers were visualized to assess how the network carries the information from input to output to understand the feature extraction. The Gradient-weighted Class Activation Mapping (Grad-CAM) was used to produce a coarse localization map highlighting the important regions in the image for predicting the target concept (CR-POPF or non-CR-POPF). And the reconstructed localization maps were named as positive and negative filters later, which were also used to evaluate the class discrimination [36]. Both Keras toolkit and Python 3.5 were needed to implement this model.

 Figure 2 

Preoperative CT scan of patients with varying CR-POPF risks. Two patients (A, B) ultimately developed CR-POPFs. FRS was 7 in patient A, and DLS was 0.78, both indicating high risk of CR-POPF; whereas a DLS of 0.60 in patient B suggested high probability of CR-POPF (0.85), despite intermediate FRS risk (FRS=4). The other two patients C and D did not develop CR-POPFs. FRS was 2 in patient C, and DLS was 0.05, both conferring low risk of CR-POPF; whereas a DLS of 0.14 in patient D showed low probability of CR-POPF (0.36) at intermediate FRS risk (FRS = 5). The first panel of each subgroup shows CT images and ROI regions used for DLS acquisition. For the second and third rows in each panel, the first and second columns shows the visualization of the activation layers of the ResCNN model, which assess the feature extraction (pancreatic parenchymal region and stump area), and the third column is the localization map that shows the important hot spots contributing to DLS for predicting non-CR-POPF/CR-POPF. The CT images were overlapped to reveal response locations. Corresponding histologic views of pancreatic stumps (patients B and D) are shown in Supplemental Figure S4.

Theranostics Image (Click on the image to enlarge.)

Statistical analysis

Candidate risk factors included radiological DLS, clinical FRS, demographic parameters (sex, BMI, etc.), anastomotic technique (PD or PJ), surgeons and hospital sites, and fistula mitigation strategies (e.g. prophylactic octreotide or transanastomotic stent usage). Risk factors for CR-POPF demonstrating significance (P < 0.05) in univariate logistic regression analysis were then applied to multivariate logistic regression modeling. Strengths of associations were presented as odds ratios (ORs), with 95% confidence intervals (CIs). We also explored the potential for the radiologic DLS to enhance the clinical FRS in predicting CR-POPFs, especially at intermediate levels of FRS risk, pursuing logistic regression analysis in the predictive models. Their comparative performances were plotted as areas under the receiver operating characteristics (ROC) curve (AUC) by Delong method[37]. Diagnostic indices, such as accuracy, sensitivity, and specificity, were obtained with ROC-derived cut-off and compared using McNemar's test. To investigate a potential relation between DLS and histologic or morphologic changes, stepwise multiple linear regression tests were conducted (detailed in Supplemental materials S4).

All computations relied on standard software applications, including R v3.5.0 (R Project for Statistical Computing, Vienna, Austria), SPSS v25 (IBM Corp, Armonk, NY, USA), and MATLAB R2019a (MathWorks, Natick, MA, USA). Statistical significance was set at P<0.05.

Usability testing

Double-blind usability testing addressed five functional aspects: predictive ability (e.g., AUC, accuracy, etc.), learnability, efficiency, satisfaction, and memorability. Following training, an open-source DLS model was released online (https://github.com/lungproject/Pancreas), offering dual tutor-assisted sample cases. The external test cohort was evaluated independently and the DLS values were calculated as directed by the tutor on their own computers.

Results

Patient characteristics

Comparisons of clinical characteristics, intraoperative data, and postoperative histology in patients with and without CR-POPFs are presented in Supplemental Table S4. CR-POPFs developed in 95 of 583 patients (16.3%), with training, validation and test datasets marked by similar outcomes (P = 0.56). Median FRS values in the presence (vs. absence) CR-POPFs were significantly higher for training (6 [IQR: 4.5-7] vs. 3 [IQR: 2-5]), validation cohort (5.5 [IQR: 4.5-7] vs. 3 [IQR: 2-5]) and test cohort (7 [IQR: 5-7] vs. 4 [IQR: 2-6]). The same was true of median DLS values in training (0.54 [IQR: 0.46-0.65] vs. 0.38 [IQR: 0.29-0.45]), validation (0.54 [IQR: 0.46-0.60] vs. 0.38 [IQR: 0.27-0.48) and test cohort (0.57 [IQR: 0.55-0.67] vs. 0.37 [IQR: 0.28-0.44) subjects (all P < 0.001) (Supplemental Figure S2).

Interrater reproducibility for volumetric segmentation agreement was expressed as dice similarity coefficient and hausdorff distance in Supplemental material S5, and clinical outcomes are detailed in Supplemental material S6.

Performance of DLS in predicting CR-POPF

Figure 2 showed CT images of four representative patients with different DLS and FRS. Through multivariate logistic regression analysis, DLS > 0.5 (OR = 12.23, 95% CI, 5.33-8.10; P < 0.001), incremental FRS increases (per point, OR = 1.43, 95% CI, 1.17-1.75; P < 0.001), octreotide use (OR = 3.70, 95% CI, 1.71-8.04; P = 0.001), and weight loss (OR = 0.34, 95% CI, 0.16-0.74; P = 0.006) proved to be independently associated with occurrences of CR-POPF (Supplemental Table S5), and DLS+FRS were generated consequently. DLS (cutoff: 0.5), FRS (cutoff: 5), and DLS+FRS (cutoff: 0.56) models were compared in differentiating CR-POPF presence or absence (Table 1). In terms of AUC, DLS (training: 0.85, 95% CI, 0.80-0.90, validation: 0.81, 95% CI, 0.72-0.89, test: 0.89, 95% CI, 0.79-0.96) significantly outperformed the FRS (training: 0.78, 95% CI, 0.72-0.84, P = 0.03, validation: 0.76, 95% CI, 0.66-0.84, P = 0.05, test: 0.73, 95% CI, 0.61-0.83, P < 0.001) (Figure 3A, B) in all three cohorts. Additionally, DLS proved higher predictive accuracy (training: 81.3% vs. 71.9%, P < 0.001; validation: 76.6% vs. 69.5%, P = 0.025) and specificity (training: 83.2% vs. 71.3%, P < 0.001; validation: 76.9% vs. 68.5%, P = 0.015) than did FRS, without significant compromise in sensitivity (training: 71.4% vs. 75.0%, P = 0.83; validation: 75.0% vs. 75.0%, P = 1.0) in the training and validation cohorts. In test cohort, the accuracy (P < 0.001), specificity (P = 0.008) and sensitivity (P < 0.001) of DLS were all significantly higher than those of FRS. Though the incorporation of FRS into DLS achieved higher AUCs of 0.87 (95% CI, 0.82-0.91), 0.85 (95% CI, 0.77-0.91) and 0.90 (95% CI, 0.80-0.96) in the training, validation and test cohorts, respectively, these improvements were not significant compared to DLS alone (training: P = 0.11, validation: P = 0.22, test: P = 0.52).

Previously published fixed CT classifiers (remnant pancreatic volume, MPD size, and pancreatic width or thickness) failed to surpass the DLS, and no further incremental gains in discriminatory capacity were achieved by adding these to the DLS (Supplemental Figure S3).

FRS risk stratification analysis

Clinical FRS values served to stratify the train cohort by FRS risk level as low (38.7%), intermediate (47.1%) and high (14.2%). In 139 patients of low FRS risk, actual occurrences of CR-POPF accounted for only 5% of patients (n = 7). The accuracies of both DLS (96.4%; n = 134) and FRS (95.0%; n = 132) (P = 0.90) approached 100%. In 51 patients of high-FRS risk, actual CR-POPFs developed in 41.2% (n = 21). The DLS fared slightly (but not significantly) better than the FRS (60.7% vs. 41.2%; P = 0.053) in diagnostic accuracy. Confusion matrices of low- and high-risk patients in all cohorts are shown in Supplemental Tables S5.

In patients of intermediate FRS risk, the diagnostic accuracy of FRS was low in training cohort (61.5%), validation (56.3%) and test (65.8%) cohorts, whereas the DLS performed significantly better (training, 79.9%, P = 0.005; validation, 70.3%, P = 0.04; test, 92.1%, P = 0.013). In specific, the clinical FRS resulted in greater misclassification of training (57/140), validation (25/50) and test (11/33) patients as false positives, compared with the DLS (training: 26/140; validation: 16/50; test: 3/33), conferring specificities of 59.3% vs. 81.4% in the training cohort (P < 0.001), 50.0% vs. 68.0% in the validation cohort (P = 0.015), and 66.7 % vs. 90.9 % in the test cohort (P = 0.038), as detailed in Figure 4 for each FRS score. A more specific confusion matrix is included as Supplemental Table S6; and details of AUC, specificity, and sensitivity are provided in Supplemental Table S7.

Biology underpinnings of DLS

Histologically, the DLS significantly correlated with less fibrosis (training: ρ = -0.60; validation: ρ = -0.65; P < 0.01) and acinar atrophy (training: ρ = -0.62; validation: ρ = -0.68; P < 0.01), as indicated in Figure 3C-F DLS positively correlated with FRS and pancreatic parameters (volume, width, thickness) but negatively correlated with MPD diameter and gland texture (Supplemental Table S8). DLS also corresponded well with the hot-spotted pancreatic parenchyma and stump areas, representing the most important regions contributing to DLS (Figure 2).

Multivariate linear regression analysis (adjusted r2 = 0.51; F = 47.8; P < 0.001) further revealed that the DLS was independently associated with pancreatic fibrosis (coefficient = -0.17; P = 0.029), MPD (coefficient = -0.45; P < 0.001), and remnant pancreatic volume (coefficient = 0.14; P = 0.012), all contributing to 51.2% of DLS variability (Supplemental Table S9).

Usability testing

In the test cohort, the user satisfaction score was 4 when using a 5-point scale. The external test cohort was prepared within 10 minutes, and the cases were run and DLS values listed within 1~2 minutes (see Supplemental material S7 for usability testing results).

Discussion

Recently, deep learning has become a popular tool for imaging analysis, fueled by its hierarchal automated learning capacities and the optimal parametric sets delivered. In utilizing the hidden layers of deep learning, the DLS described herein may well serve as a new, highly accessible and approachable means of routine patient examination. Within 1~2 minutes, DLS can be automatically generated from preoperative CT scans with excellent performance (AUC = 0.90 in external cohort), the latter being mandatory for nearly all patients in this setting.

 Table 1 

Performances of various predictive models in training, validation and test cohorts

AUC (95% CI)Accuracy (95% CI)Sensitivity (95% CI)Specificity (95% CI)
Deep-learning score (DLS)
Training0.85 (0.80, 0.90)81.3 (77.2, 85.2)71.4 (58.9, 82.1)83.2 (78.9, 87.5)
Validation0.81 (0.72, 0.89)76.6 (70.1, 83.1)75.0 (58.3, 91.7)76.9 (70.0, 83.9)
Test0.89 (0.79, 0.96)87.1 (86.8, 87.5)86.7 (59.5, 98.3)87.3 (75.5, 94.7)
Fistula Risk Score (FRS)
Training0.78 (0.72, 0.84)71.9 (67.4, 76.3)75.0 (64.3, 85.7)71.3 (66.3, 76.2)
Validation0.76 (0.66, 0.84)69.5 (62.3, 76.6)75.0 (58.3, 91.7)68.5 (60.8, 76.2)
Test0.73 (0.61, 0.83)72.9 (72.3, 73.4)73.3 (44.9, 92.2)72.7 (59.0, 83.9)
DLS+FRS
Training0.87 (0.82, 0.91)81.3 (77.4, 85.5)82.1 (71.4, 91.1)81.2 (76.9, 85.8)
Validation0.85 (0.77, 0.91)75.9 (68.8, 83.1)79.2 (62.5, 91.7)75.4 (66.9, 83.1)
Test0.90 (0.80, 0.96)90.0 (89.7, 90.3)86.7 (59.5, 98.3)90.9 (80.0, 97.0)

AUC, area under receiver operating characteristic (ROC) curve; ACC, accuracy; CI, confidence interval.

 Figure 3 

Performance of DLS in various cohorts. Receiver operating characteristic (ROC) curves of DLS, FRS, and DLS+FRS prediction models are plotted for the training cohort (A) and validation cohort (B). DLS is also shown to be significantly correlated with (C, D) fibrosis and (E, F) acinar atrophy in these two cohorts.

Theranostics Image (Click on the image to enlarge.)
 Figure 4 

Plotting of clinically relevant postoperative pancreatic fistula (CR-POPF) relative to Fistula Risk Score (FRS) and deep learning score (DLS) values in patients at intermediate risk. In (A) the training cohort (n = 169) and (B) validation cohort (n = 64), totals above each column reflect occurrences according to risk score. Distributions of actual CR-POPFs are further stratified by FRS and DLS, shown as black (only DLS correct), white (DLS and FRS incorrect), gray (only FRS correct), and dotted (DLS and FRS correct) areas. Both DLS and FRS predicted more cases correctly at FRS values of 3 and 4 than at values of 5 and 6 (larger white areas). The DLS-only correctly classified cases (black areas) distributed mainly in the patients without CR-POPF, which was more than the FRS-only correctly classified cases (gray areas). Note. Intermediate risk: FRS is within the range of 3-6 points. 0: No CR-POPF; 1: CR-POPF.

Theranostics Image (Click on the image to enlarge.)

Other traditional methods of quantitative imaging reported to date are burdened by complicated models and tedious manual measurements, whereas the DLS consolidates all local risk factors into one simplified model. We directly compared the DLS with almost all CT classifiers [5, 20-22, 38-41] cited in literature (remnant pancreatic volume, MPD diameter, and pancreatic width or thickness), discounting the common belief that deep learning always prevails. However, none of the above showed any superiority to DLS or any incremental benefit in combination with DLS, perhaps reflecting their innate correlations with DLS. One recent study [42] has claimed good results in 80 training (AUC = 0.82) and 37 test (AUC = 0.76) subjects by harnessing conventional radiomics analysis. Similarly, Kambakamba et al utilized machine learning-based texture analysis, which could achieve an impressive AUC of 0.95 in predicting POPF on 110 patients from a single institution [43]. However, the present effort is the first to explore a deep-learning based method with significantly larger cohorts from different institutions, which is simpler to use and easier to interpret using one single biomarker (DLS, or the probability). Compared to classical textural analysis, the visualization of the developed deep learning model could highlight the important regions in the image, such as the pancreatic parenchyma and stump areas, for predicting CR-POPF or non-CR-POPF, which is helpful in revealing the biology underpinnings of the DLS. Additionally, our study further emphasized the super prediction ability in the sub-cohort with FRS-based intermediate risk, which is a grey zone warranting more objective and reliable predictors.

The radiologic DLS was designed to add more objective local features to the FRS and fully utilize preoperative CT studies, beyond mere recognition of MPD dilatation. Moreover, pancreatic texture (another determinant of FRS) is indirectly ascertained via DLS, perhaps enabling a strong sense of visually glandular fibrosis or atrophy on CT images. Our analysis also indicated that DLS may outperform the FRS at certain FRS risk levels. In low-risk subjects, in whom CR-POPFs are unlikely, FRS performance was excellent; whereas FRS values >6 clearly equated with high risk. Unfortunately, ~50% of cases gravitate to intermediate FRS risk, according to our data and that of others. CR-COPF rates at this intermediary level (~15%) and overall (~16%) are thus indistinguishable, hindering POPF prevention decisions. Contrastingly, the DLS significantly outperforms the FRS in excluding those patients who lack CR-POPFs, significantly improving specificity in the intermediate-risk group without altering sensitivity. Hence, we see no need for DLS calculation at very low FRS values, except if used as an alternative preoperative biomarker. The DLS is otherwise highly recommended at intermediate FRS levels and may help exclude falsely high risk assessments at FRS values > 6. In short, the DLS is intended for preoperative evaluation of CR-POPF risk in any manner deemed surgically appropriate. As in prior studies, octreotide use for mitigation and nutritional status (weight loss or BMI status) were also associated herein with CR-POPF risk. A comprehensive analysis of all parameters is likely to benefit in predicting postoperative CR-POPF.

Visualization of the convolution filters helped us demystify what has been learned in devising the DLS. Hot-spotted regions, such as remnant pancreatic volume and stump areas (Figure 2) were important contributors to the final concept, consistent with traditional CT [5] or MR [44] indices of these morphologic features. At a macroscopic level, DLS (Table S8) was significantly associated with MPD diameter, remnant pancreatic volume, pancreatic thickness or width, and softer gland texture, while showing a negative relation with pancreatic fibrosis and acinar atrophy in histologic preparations. Higher DLS thus corresponds with soft, non-fibrotic, large-sized pancreatic remnants and a non-distended MPD, signifying full tissue viability and fluid productivity and presenting greater challenge in suturing. Conversely, lower DLS attests to a hardened, fibrotic, and atrophic pancreas, with a dilated MPD and limited secretory capacity for easier anastomosis and lower risk of CR-POPF.

We do acknowledge some limitations of this study. The first issue is that our data accrued from four separate and independent centers attended by various research associates. Although different CT modalities were used, the scanners were standardized to 64-slice capability or better, with slice thickness <3 mm. The Hounsfield scale, an international standardization, was also upheld in all CT systems to control inner variabilities of different scanners. The FRS, however, relies on subjective evaluations of individual surgeons, creating opportunities for data inconsistencies. Furthermore, DLS and FRS values were generated exclusively at high-volume institutions and may not carry their predictive weight in lower-volume facilities equipped with earlier generations of scanners or staffed by fewer surgeons.

Conclusion

In conclusion, quantitative preoperative CT assessment using AI is strongly predictive of CR-POPF in patients with pancreaticoenteric anastomoses following PD. The automated scores could reflect histomorphologic features pertaining to pancreatic duct, remnant pancreatic tissue volume and parenchymal fibrosis. The DLS is particularly helpful for patients with intermediate FRS risks in gauging the potential for CR-POPF. Future efforts would focus on its integration into image archiving and communication systems used for radiologic reporting or incorporation into surgical planning software to accommodate the preferences of surgeons, thus optimizing preoperative strategies, intraoperative decision-making, and even postoperative care.

Abbreviations

CR-POPF, clinically relevant postoperative pancreatic fistula; CE-CT, contrast-enhanced computed tomography; PD, pancreatoduodenectomy; PJ, pancreaticojejunostomy; MPD, main pancreatic duct, FRS, Fistula Risk Score; DLS, deep-learning score; AUC, areas under the receiver operating characteristics curve; CI, confidence intervals; IQR, interquartile range; OR, odds ratio.

Supplementary Material

Attachment

Supplementary figures and tables.

Acknowledgements

This paper is supported by the National Natural Science Foundation of China (Nos: 81771802, 81771893) and the fundamental research funds for the central universities (No.YWF-20-BJ-J-1048).

Author Contributions

Acquisition of data: Yu Shi, Chang Liu, Feng Gao, Yafei Qi, Hong Lu, Zaiyi Liu, Xianyi Zhang, Xiaoli Cai, Chang Liu, Ruo Yun Ji;

Analysis and interpretation of data: Wei Mu, Yu Shi, Jie Tian, Feng Gao, Yafei Qi;

Drafting of manuscript: Wei Mu, Yu Shi, Feng Gao, Yafei Qi, Zaiyi Liu, Hong Lu;

Critical revision: Wei Mu, Yu Shi, Chang Liu, Feng Gao, Yafei Qi, Hong Lu, Zaiyi Liu, Xianyi Zhang, Xiaoli Cai, Chang Liu, Ruo Yun Ji, Jie Tian.

Competing Interests

The authors have declared that no competing interest exists.

References

1. Bames SA, Lillemoe KD, Kaufman HS, Sauter PK, Yeo CJ, Talamini MA. et al. Pancreaticoduodenectomy for benign disease. Am J Surg. 1996;171:131-5

2. Nakeeb A, Lillemoe KD, Cameron JL. The role of pancreaticoduodenectomy for locally recurrent or metastatic carcinoma to the periampullary region. J Am Coll Surg. 1995;180:188-92

3. Nahm CB, Connor SJ, Samra JS, Mittal A. Postoperative pancreatic fistula: a review of traditional and emerging concepts. Clin Exp Gastroenterol. 2018;11:105-18

4. Hill JS, Zhou Z, Simons JP, Ng SC, McDade TP, Whalen GF. et al. A simple risk score to predict in-hospital mortality after pancreatic resection for cancer. Ann Surg Oncol. 2010;17:1802-7

5. Kanda M, Fujii T, Suenaga M, Takami H, Hattori M, Inokawa Y. et al. Estimated pancreatic parenchymal remnant volume accurately predicts clinically relevant pancreatic fistula after pancreatoduodenectomy. Surgery. 2014;156:601-10

6. Nakao A, Kanzaki A, Fujii T, Kodera Y, Yamada S, Sugimoto H. et al. Correlation between radiographic classification and pathological grade of portal vein wall invasion in pancreatic head cancer. Ann Surg. 2012;255:103-8

7. Büchler MW, Wagner M, Schmied BM, Uhl W, Friess H, Z'graggen K. Changes in morbidity after pancreatic resection: toward the end of completion pancreatectomy. Arch Surg. 2003;138:1310-4

8. Poon RTP, Lo SH, Fong D, Fan ST, Wong J. Prevention of pancreatic anastomotic leakage after pancreaticoduodenectomy. Am J Surg. 2002;183:42-52

9. Sato R, Kishiwada M, Kuriyama N, Azumi Y, Mizuno S, Usui M. et al. Paradoxical impact of the remnant pancreatic volume and infectious complications on the development of nonalcoholic fatty liver disease after pancreaticoduodenectomy. J Hepatobiliary Pancreat Sci. 2014;21:562-72

10. Shimoda M, Katoh M, Yukihiro I, Kita J, Sawada T, Kubota K. Body mass index is a risk factor of pancreatic fistula after pancreaticoduodenectomy. Am Surg. 2012;78:190-4

11. Zarzavadjian Le Bian A, Fuks D, Montali F, Cesaretti M, Costi R, Wind P. et al. Predicting the Severity of Pancreatic Fistula after Pancreaticoduodenectomy: Overweight and Blood Loss as Independent Risk Factors: Retrospective Analysis of 277 Patients. Surg Infect (Larchmt). 2019;20:486-91

12. Eshmuminov D, Schneider MA, Tschuor C, Raptis DA, Kambakamba P, Muller X. et al. Systematic review and meta-analysis of postoperative pancreatic fistula rates using the updated 2016 International Study Group Pancreatic Fistula definition in patients undergoing pancreatic resection with soft and hard pancreatic texture. HPB (Oxford). 2018;20:992-1003

13. Smits FJ, Molenaar IQ, Besselink MG, Borel Rinkes IHM, van Eijck CHJ, Busch OR. et al. Early recognition of clinically relevant postoperative pancreatic fistula: a systematic review. HPB (Oxford). 2019

14. Kiyochi H, Matsukage S, Nakamura T, Ishida N, Takada Y, Kajiwara S. Pathologic Assessment of Pancreatic Fibrosis for Objective Prediction of Pancreatic Fistula and Management of Prophylactic Drain Removal After Pancreaticoduodenectomy. World J Surg. 2015;39:2967-74

15. Hong W, Ha HI, Lee JW, Lee SM, Kim MJ. Measurement of Pancreatic Fat Fraction by CT Histogram Analysis to Predict Pancreatic Fistula after Pancreaticoduodenectomy. Korean J Radiol. 2019;20:599-608

16. Mungroop TH, Klompmaker S, Groot Koerkamp B, Besselink MG, Dutch Pancreatic Cancer G. Added Value of Body Fat Distribution in Predicting Clinically Significant Pancreatic Fistula in the a-FRS Following Pancreatoduodenectomy Currently Unclear. Ann Surg. 2019;269:e2-e3

17. Callery MP, Pratt WB, Kent TS, Chaikof EL, Vollmer Jr CM. A prospectively validated clinical risk score accurately predicts pancreatic fistula after pancreatoduodenectomy. J Am Coll Surg. 2013;216:1-14

18. Miller BC, Christein JD, Behrman SW, Drebin JA, Pratt WB, Callery MP. et al. A multi-institutional external validation of the fistula risk score for pancreatoduodenectomy. J Gastrointest Surg. 2014;18:172-80

19. McMillan MT, Soi S, Asbun HJ, Ball CG, Bassi C, Beane JD. et al. Risk-adjusted outcomes of clinically relevant pancreatic fistula following pancreatoduodenectomy: a model for performance evaluation. Ann Surg. 2016;264:344-52

20. Okano K, Oshima M, Kakinoki K, Yamamoto N, Akamoto S, Yachida S. et al. Pancreatic thickness as a predictive factor for postoperative pancreatic fistula after distal pancreatectomy using an endopath stapler. Surgery Today. 2013;43:141-7

21. Frozanpor F, Loizou L, Ansorge C, Segersvärd R, Lundell L, Albiin N. Preoperative Pancreas CT/MRI Characteristics Predict Fistula Rate after Pancreaticoduodenectomy. World J Surg. 2012;36:1858-65

22. Kirihara Y, Takahashi N, Hashimoto Y, Sclabas GM, Khan S, Moriya T. et al. Prediction of Pancreatic Anastomotic Failure After Pancreatoduodenectomy: The Use of Preoperative, Quantitative Computed Tomography to Measure Remnant Pancreatic Volume and Body Composition. Ann Surg. 2013;257:512-9

23. Ariji Y, Fukuda M, Kise Y, Nozawa M, Yanashita Y, Fujita H. et al. Contrast-enhanced computed tomography image assessment of cervical lymph node metastasis in patients with oral cancer by using a deep learning system of artificial intelligence. Oral Surg Oral Med Oral Pathol Oral Radiol. 2019;127:458-63

24. Yang X, Wu L, Ye W, Zhao K, Wang Y, Liu W. et al. Deep Learning Signature Based on Staging CT for Preoperative Prediction of Sentinel Lymph Node Metastasis in Breast Cancer. Academic Radiology.

25. Kann BH, Hicks DF, Payabvash S, Mahajan A, Du J, Gupta V. et al. Multi-institutional validation of deep learning for pretreatment identification of extranodal extension in head and neck squamous cell carcinoma. Journal of Clinical Oncology. 2020;38:1304-11

26. Xu J, Xue K, Zhang K. Current status and future trends of clinical diagnoses via image-based deep learning. Theranostics. 2019;9:7556-65

27. von Elm E, Altman DG, Egger M, Pocock SJ, Gotzsche PC, Vandenbroucke JP. et al. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. PLoS Med. 2007;4:e296

28. Shen Y, Jin W. Reconstruction by Pancreaticogastrostomy versus Pancreaticojejunostomy following Pancreaticoduodenectomy: A Meta-Analysis of Randomized Controlled Trials. Gastroenterol Res Pract. 2012;2012:627095

29. Bassi C, Dervenis C, Butturini G, Fingerhut A, Yeo C, Izbicki J. et al. Postoperative pancreatic fistula: an international study group (ISGPF) definition. Surgery. 2005;138:8-13

30. Bassi C, Marchegiani G, Dervenis C, Sarr M, Abu Hilal M, Adham M. et al. The 2016 update of the International Study Group (ISGPS) definition and grading of postoperative pancreatic fistula: 11 Years After. Surgery. 2017;161:584-91

31. Hatano M, Watanabe J, Kushihata F, Tohyama T, Kuroda T, Koizumi M. et al. Quantification of pancreatic stiffness on intraoperative ultrasound elastography and evaluation of its relationship with postoperative pancreatic fistula. Int Surg. 2015;100:497-502

32. Watanabe H, Kanematsu M, Tanaka K, Osada S, Tomita H, Hara A. et al. Fibrosis and Postoperative Fistula of the Pancreas: Correlation with MR Imaging Findings—Preliminary Results. Radiology. 2013;270:791-9

33. Gaujoux S, Cortes A, Couvelard A, Noullet S, Clavel L, Rebours V. et al. Fatty pancreas and increased body mass index are risk factors of pancreatic fistula after pancreaticoduodenectomy. Surgery. 2010;148:15-23

34. Tempero MA, Malafa MP, Al-Hawary M, Asbun H, Bain A, Behrman SW. et al. Pancreatic Adenocarcinoma, Version 2.2017, NCCN Clinical Practice Guidelines in Oncology. J Natl Compr Canc Netw. 2017;15:1028-61

35. Tempero MA, Malafa MP, Chiorean EG, Czito B, Scaife C, Narang AK. et al. Pancreatic Adenocarcinoma, Version 1.2019. J Natl Compr Canc Netw. 2019;17:202-10

36. Montavon G, Samek W, Müller K-R. Methods for interpreting and understanding deep neural networks. Digit Signal Process. 2018;73:1-15

37. DeLong E, DeLong D, Clarke-Pearson D. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biome trics 1988; 44: 837-45. Am J Nephrol. 2017;45:400-8

38. Tranchart H, Gaujoux S, Rebours V, Vullierme M-P, Dokmak S, Levy P. et al. Preoperative CT Scan Helps to Predict the Occurrence of Severe Pancreatic Fistula After Pancreaticoduodenectomy. Ann Surg. 2012;256:139-45

39. Ridolfi C, Angiolini MR, Gavazzi F, Spaggiari P, Tinti MC, Uccelli F. et al. Morphohistological features of pancreatic stump are the main determinant of pancreatic fistula after pancreatoduodenectomy. Biomed Res Int. 2014;2014:641239

40. Wellner UF, Kayser G, Lapshyn H, Sick O, Makowiec F, Höppner J. et al. A simple scoring system based on clinical factors related to pancreatic texture predicts postoperative pancreatic fistula preoperatively. HPB(Oxford). 2010;12:696-702

41. Kakizawa N, Noda H, Watanabe F, Ichida K, Suzuki K, Rikiyama T. A High Abdominal Aortic Calcification Score on CT is a Risk Factor for Postoperative Pancreatic Fistula in Elderly Patients Undergoing Pancreaticoduodenectomy. World J Surg. 2018;42:1129-37

42. Zhang W, Cai W, He B, Xiang N, Fang C, Jia F. A radiomics-based formula for the preoperative prediction of postoperative pancreatic fistula in patients with pancreaticoduodenectomy. Cancer Manag Res. 2018;10:6469-78

43. Kambakamba P, Mannil M, Herrera PE, Müller PC, Kuemmerli C, Linecker M. et al. The potential of machine learning to predict postoperative pancreatic fistula based on preoperative, non-contrast-enhanced CT: A proof-of-principle study. Surgery. 2020;167:448-54

44. Shi Y, Liu Y, Gao F, Liu Y, Tao S, Li Y. et al. Pancreatic Stiffness Quantified with MR Elastography: Relationship to Postoperative Pancreatic Fistula after Pancreaticoenteric Anastomosis. Radiology. 2018;288:476-84

Author contact

Corresponding address Corresponding author: Yu Shi, MD, Department of Radiology, Shengjing Hospital, China Medical University, No. 36, Sanhao Street, Heping District, Shenyang, 110004, P.R. China. Tel.: +86 24 96615 73211; E-mail address: shiyorg.


Received 2020-6-18
Accepted 2020-7-22
Published 2020-8-1