Skip to main content
  • Research article
  • Open access
  • Published:

Radiomic nomogram based on MRI to predict grade of branching type intraductal papillary mucinous neoplasms of the pancreas: a multicenter study

Abstract

Background

Accurate diagnosis of high-grade branching type intraductal papillary mucinous neoplasms (BD-IPMNs) is challenging in clinical setting. We aimed to construct and validate a nomogram combining clinical characteristics and radiomic features for the preoperative prediction of low and high-grade in BD-IPMNs.

Methods

Two hundred and two patients from three medical centers were enrolled. The high-grade BD-IPMN group comprised patients with high-grade dysplasia and invasive carcinoma in BD-IPMN (n = 50). The training cohort comprised patients from the first medical center (n = 103), and the external independent validation cohorts comprised patients from the second and third medical centers (n = 48 and 51). Within 3 months prior to surgery, all patients were subjected to magnetic resonance examination. The volume of interest was delineated on T1-weighted (T1-w) imaging, T2-weighted (T2-w) imaging, and contrast-enhanced T1-weighted (CET1-w) imaging, respectively, on each tumor slice. Quantitative image features were extracted using MITK software (G.E.). The Mann-Whitney U test or independent-sample t-test, and LASSO regression, were applied for data dimension reduction, after which a radiomic signature was constructed for grade assessment. Based on the training cohort, we developed a combined nomogram model incorporating clinical variables and the radiomic signature. Decision curve analysis (DCA), a receiver operating characteristic curve (ROC), a calibration curve, and the area under the ROC curve (AUC) were used to evaluate the utility of the constructed model based on the external independent validation cohorts.

Results

To predict tumor grade, we developed a nine-feature-combined radiomic signature. For the radiomic signature, the AUC values of high-grade disease were 0.836 in the training cohort, 0.811 in external validation cohort 1, and 0.822 in external validation cohort 2. The CA19–9 level and main pancreatic duct size were identified as independent parameters of high-grade of BD-IPMNs using multivariate logistic regression analysis. The CA19–9 level and main pancreatic duct size were then used to construct the radiomic nomogram. Using the radiomic nomogram, the high-grade disease-associated AUC values were 0.903 (training cohort), 0.884 (external validation cohort 1), and 0.876 (external validation cohort 2). The clinical utility of the developed nomogram was verified using the calibration curve and DCA.

Conclusions

The developed radiomic nomogram model could effectively distinguish high-grade patients with BD-IPMNs preoperatively. This preoperative identification might improve treatment methods and promote personalized therapy in patients with BD-IPMNs.

Introduction

The pancreatic ductal system mucinous epithelium can develop mucin-producing tumors, such as intraductal papillary mucinous neoplasms (IPMNs) of the pancreas. IPMNs represent approximately 21–33% of cystic neoplasms and is one of the precursors of pancreatic cancer that is identifiable radiographically [1,2,3]. In the past two decades, we found that the detection rate and incidence of IPMNs have increased significantly as a result of the use of advanced diagnostic imaging technology [1, 4, 5]. When a patients is suspected of having an IPMN, determining the grade of malignancy for an individual patient permits decisions regarding whether surgery or surveillance is appropriate to be made. However, although clinicians are experienced at diagnosing and treating IPMNs, it is still challenging to distinguish high-grade IPMNs (i.e., high-grade dysplasia (HGD) to invasive carcinoma) from low-grade IPMNs (i.e., low-grade dysplasia (LGD) to intermediate-grade dysplasia (IGD)) before surgery.

Currently, laboratory tests, endoscopy, cytology, and imaging technologies play the main roles in differentiating between high-grade and low-grade IPMNs. The International Association of pancreatic diseases (IAP) recommends active surgical treatment for IPMNs of the main duct type and mixed type, while for IPMNs of the branch duct type, according to the revised 2017 international consensus guidelines, surgery is recommended for tumors with indicative features, such as mural nodules, cyst size > 3 cm, main pancreatic duct (MPD) size > 5 mm, rapid cyst growth (≥ 5 mm over 2 years), and increased serum carbohydrate antigen (CA19–9 > 37 ng/ml) levels [6, 7]. The IAP guidelines noted that the mean prevalence of invasive malignancy in BD-IPMNs was 17.7% (1–37%) [8]. In this setting, a considerable number of patients with benign lesions received unnecessary invasive surgery, and the existing grade assessment system, with unsatisfactory specificity and positive predictive value, remains unreliable [9, 10]. Therefore, there is an urgent need for a highly sensitive and specific preoperative prediction system to help establish individualized treatment decisions. Medical imaging produces essential information for the preoperative assessment of BD-IPMNs. It has been suggested that computed tomography (CT)-derived radiological features could assess the grade of BD-IPMNs objectively [11,12,13]. Radiomics is the quantitative analysis of images that comprise a large number of features combined with machine learning [14]. Radiomics has demonstrated potential utility in oncological imaging in prognosis, detection, and differential diagnosis assessments, such as in the lung [15], breast [16], prostate [17, 18] and liver [19]. However, to the best of our knowledge, there has been no previous research using magnetic resonance (MRI)-derived radiomics in the grade assessment in patients with BD-IPMNs. Some scoring systems or nomograms to predict malignancy from clinical variables have been suggested [20,21,22]; however, these methods had limitations, e.g., the lack of external or internal validation of their clinical efficacy.

Thus, the present study aimed to construct a predictive model that integrated clinical variables and a radiomic signature for preoperative grade assessment in patients with BD-IPMNs. This improved preoperative evaluation model could spare low-grade patients from potentially morbid surgery and permit high-grade patients to undergo resection before transformation into an invasive phenotype [23].

Methods

Workflow

The research workflow is shown in Fig. 1, and comprises four parts: Acquiring the images, segmenting the region of interest (ROI), extracting features, and prediction model construction. We acquired T1-w, T2-w MRI images and CET1-w images, and radiologists outlined the tumor area manually on all image slices. Then, the quantitative radiomic features were extracted from the ROI, after which a machine learning model was established to assess the grade of BD-IPMNs.

Fig. 1
figure 1

Workflow of the development of the radiomic nomogram

The machine learning model-based tumor pathological grade assessment was developed and validated using three separate datasets. Medical center A (n = 103) data were used as the training dataset to construct the tumor pathological grade assessment model. Medical center B and C (n = 48 and 51) data were used as independent validation datasets to test the developed model.

After extraction of the quantitative radiomics features from the ROI, least absolute shrinkage and selection operator (LASSO) regression, Spearman correlation analysis, analysis of variance (ANOVA) tests, and Mann–Whitney U test, were applied to select the best radiomics features from the radiomic signature. Then, to construct the tumor pathological grade assessment nomogram, multivariate logistic regression was used to integrate the clinical variables and the radiomic signature.

Patients

The Institutional Review Boards of the three centers approved this retrospective study. The provision of signed informed consent was waived. This study was conducted following the tenets of the Declaration of Helsinki. Between Mar. 2012 and Feb. 2020, patients who were reported to have BD-IPMNs on pathological assessment and who underwent an MRI scan of their pancreas before surgery were included in this study. If the branch duct-type information was not included on the pathology report, patients with BD-IPMNs were selected according to their pre-operative MRI report. Among them, 91 patients were excluded as follows: Patients lacking complete clinical data (n = 23), suboptimal MRI image quality caused by severe motion artifacts (n = 9,) and patients who did not receive an MRI scan within 3 months before surgery (n = 59)(Fig. S1).

A multidisciplinary team comprising, surgeons, oncologists, and radiologists assessed the patients. Clinical characteristics (e.g., sex, symptoms (abdominal pain or fatigue), age, mural nodule, cyst size, main pancreatic duct (MPD) size, CA19–9 and CEA level) were obtained through review of the clinical data by a surgeon and a radiologist with over 10 years of clinical experience. The size of the MPD was determined at the maximum dilation point of the pancreatic duct. Mural nodules were defined as any solid papillary protuberances in the cyst. The serum level of CEA and CA19–9 were measured using ELISA. The serum CA19–9 level was judged to be elevated if it was higher than the upper limit of normal (37 U/mL). The serum CEA level was judged to be elevated if it was higher than the upper limit of normal (5 U/mL). The pathological grade (LGD to invasive carcinoma) was assessed by two pathologists with over 10 years of experience in abdominal tumor diagnosis according to the 2010 World Health Organization classification [24].

Acquisition of MRI images

Four different scanning images: preoperative T1-w, T2-w, CET1-w arterial phase, and CET1-w portal venous phase were acquired for all images. T1-w and CET1-w images were collected during breath holding and T2 images were collected by breath triggering. The details of MRI image acquisition are shown in Supplementary material I.

Tumor segmentation and Radiomics feature extraction

To extract the features, portal venous phase imaging, arterial phase imaging, and pre-contrast T1-w, T2-w were used (Fig. S2). Radiologists manually contoured the ROIs on the MRI images. If possible, reference images were acquired using contrast-enhanced computed tomography (CT). Image standardization before ROI segmentation was conducted, which makes the thickness of the layer consistent, thus ensuring that the size and bedding of each ROI segmentation are consistent. Before ROI segmentation, the two radiologists who performed the analysis (Y.J.W. and W.Y) were blinded to the clinical outcome. The whole tumor was segmented manually using ITK-SNAP (www.itk-snap.org) on all slices of the tumor [25]. According to previously published studies, when there were multiple tumors, the tumors with the largest diameter were selected for analysis [26]. One of the radiologists drew the tumor boundary, which was verified by the other radiologist.

The MITK software (Medical Imaging Interaction Tookit 3.1.0.A, GE Healthcare) was used to extract the radiomic features from the three-dimensional ROIs, resulting in 328 features. The extracted features included: Histograms, texture parameters, RLM (run length matrix), GLCM (gray level co-occurrence matrix), and form factor parameters. Each factor’s average value was subtracted from all extracted radiometric features, and divided by the standard deviation value (Z score normalization), which eliminates the limitation imposed by each feature’s units. To standardize the different scales used to process the variables, each feature’s average value was subtracted from all the radiomic features in the training data set, and then divided by their standard deviation values, respectively. Then, using the mean and standard deviation values derived from the training dataset, the same normalization method was applied.

To construct a realistic radiomic signature that combines the most appropriate radiomic features we used the LASSO regression method, as described by Monica and Kumamaru, to select the most nonredundant and robust radiomic features [27, 28]. Supplementary material II describes the details of the LASSO method. To predict the grade assessment for each patient, the best selected radiomic features from the training cohort were recalculated using a linear combination of selected features weighted by their respective coefficients. Receiver operating characteristic (ROC) curve and area under the ROC curve (AUC) were used to evaluate the predictive accuracy of the developed radiomic signature. The reproducibility of radiomics feature extraction was assessed using intra- and inter-class correlation coefficients (ICCs). Initially, 30 ROIs were chosen randomly from each sequence. To calculate the intra-observer ICC, reader A segmentation was repeated at a 7-day interval and compared with the original segmentation. Comparing the segment extraction of reader B with that of reader A (first time) was used to calculate the inter-observer ICC. An ICC value greater than 0.8 was considered to show good consistency of feature extraction [29].

Radiomic Nomogram construction

After univariate logistic regression, independent predictors of high-grade BD-IPMNs were selected using multivariate logistic regression analysis, and then the new radiomic nomogram was established using these independent predictors. Next, to test the calibration and recognition performance of the radiomic nomogram, the training and validation sets were used. The calibration curve displayed the performance characteristics of the multimodal radiomic nomogram models graphically. The predictive accuracy of the combined nomogram model was indicated by the degree of overlap between the diagonal in the graph and the calibration curve. The independent validation cohorts of 48 and 51 patients from the second and third medical centers were used to validate the radiomic nomogram. In the three cohorts (training and two independent validation cohorts), the clinical utility of the radiomic nomogram model was assessed using decision curve analysis (DCA). Supplementary material III shows the details of the DCA method.

Statistical analysis

SPSS 18.0 (IBM, Armonk, NY, USA), R software (v. 3.5.1, Vienna, Austria), and MedCalc software v. 15.2.2 (https://www.medcalc.org/) were used to perform all the statistical analyses. To compare continuous variables, Mann–Whitney U-tests and the independent-sample t-test were performed when appropriate. To compare categorical variables between groups, Fisher’s exact test or a chi-squared test was used. To assess the impact of variations between intra- and inter-readers in the extracted radiomics features, intra- and inter-class correlation coefficients (ICC) were determined. In the ROC analysis, the best threshold was determined using the Youden Index. For the developed logistic regression models, the goodness of fit was examined using the Hosmer-Lemeshow test. The “glmnet” package was used to carry out the LASSO regression analysis. The “pROC” package was used to produce the ROC plots. The “rms” package was used to plot the nomogram and calibration curve. For DCA, the “rmda” package was used. Statistical significance was accepted at P <  0.05.

Results

Patients characteristics

Among the three medical centers, we recruited 202 patients diagnosed with BD-IPMNs. The training cohort comprised patients from the first medical center (n = 103). The external independent validation cohorts comprised the patients from the second and third medical centers (n = 48 and 51). The characteristics of the patients in the training and validation datasets were evenly distributed (Table S1). There were no statistically significant differences in BD-IPMN pathological grade assessment and clinical characteristics (sex, age, symptoms, mural nodule, cyst size, MPD size, CA19–9, and CEA) between the training and validation datasets. Pathologically, high-grade BD-IPMNs were detected in 24.8% patients. The detailed distribution of clinical characteristics in the low-grade and high-grade groups is summarized in Table 1. In the training and validation datasets, we noted significant differences for the cyst size, mural nodule, MPD size, and CA19–9 between the low grade and high grade groups.

Table 1 Characteristics of patients with IPMNs in the low-grade group and high-grade group

Selection of Radiomics features and construction of the Radiomics signature

In the training dataset, 334 statistically significant features (P <  0.05) between the low-grade and high-grade groups were identified from 1312 texture features, from which 51 features were identified by Spearman correlation analysis. To construct the radiomics signature, LASSO was used to select the nine most valuable texture features (Fig. S3), including five GLCM features, two histogram features, one texture parameter feature, and one form factor feature. The LASSO regression method derived the coefficient for each selected feature. Then the radiomic features classification and the calculation of texture features after dimension reduction are carried out (Table S2 and Fig. S4). Supplementary material II shows the details of the selected features and the formula used to calculate the radiomic signature. The developed radiomic signature model produced a good result when predicting the histological grade (LGD/IGD vs. HGD/associated invasive carcinoma), resulting in an AUC of 0.836 in the training set (95% confidence interval (CI), 0.750–0.901), 0.811 in validation set 1 (95% CI, 0.671–0.909), and 0.822 in validation set 2 (95% CI, 0.690–0.915). Figure 2a–c shows the ROC curves of the radiomic signature based on the three datasets. Next, we determined the quantitative scores of the radiomic signature for each patient with respect to the classification of grade assessment of BD-IPMNs, to show the effectiveness of the radiomic signature model at the individual level (Fig. 3a–c), the percentage of patients in the low grade category whose rad scores overlap with high grade category were 58.4% in the training set, 60.5% in validation set 1 and 51.4% in validation set 2. Table 2 shows the details of the performance evaluation of the radiomic signature.

Fig. 2
figure 2

Accuracy evaluation of the radiomic signature. The radiomic signature predicted high-grade of BD-IPMNs in the training cohort (a) (AUC = 0.836), external validation cohort 1 (b) (AUC = 0.811) and external validation cohort 2(c) (AUC = 0822)

Fig. 3
figure 3

Rad-score prediction of grade of BD-IPMNs. The rad-score depicted using score dot diagrams in a the training cohort, b external validation cohort 1, and c external validation cohort 2. Red indicates a high grade of BD-IPMNs, blue represents a low grade of BD-IPMNs. A high score indicates a high likelihood of high-grade BD-IPMNs

Table 2 Univariate and multivariate logistic regression analysis of the radiomic signature and preoperative clinical parameters

Construction of the combined Nomogram

In the univariate analysis of the training cohort, the low- and high-grade groups showed significant differences for the CA19–9 level, largest cyst size, size of the MPD, mural nodule, and the radiomic signature. The AUC value for all the significant variables in the univariate analysis showed in Table 3. Multivariate logistic analysis showed that the size of the MPD (odds ratio (OR): 4.263, 95% CI: 1.762–10.3113, P = 0.001), the CA19–9 level (OR: 8.402, 95% CI: 1.622–43.524, P = 0.011), and the radiomic signature (OR: 3.434, 95% CI: 1.638–7.197, P = 0.001) were independent parameters of high-grade of BD-IPMN. Therefore, a radiomics nomogram model incorporating the developed radiomics signature with the size of the MPD and the CA19–9 level was constructed. A weighted number of points was assigned to each factor. The nomogram was then used to calculate the total point score of each patient, which was analyzed for its correlation with the estimated probability of high-grade of BD-IPMNs (Fig. 4a).

Table 3 AUC values for all the significant variables in univariate analysis of three groups
Fig. 4
figure 4

Radiomics nomogram a to predict high-grade BD-IPMNs. The radiomic nomogram was constructed using the rad-score, the CA19–9 level, and the size of the MPD in the data from the training cohort. Calibration for high-grade BD-IPMNs in b the training cohort, c external validation cohort 1, and d external validation cohort 2. The optimal nomogram is indicated by the dotted line reference line. The performance of the radiomics nomogram in high-grade prediction is indicated by the dashed line, and bias in the nomogram is corrected using the solid line

Radiomic Nomogram performance evaluation

ROC analysis was used to confirm the utility of the combined nomogram, resulting in an AUC of 0.903 (95% CI, 0.828–0.952), a sensitivity of 0.734, and a specificity of 0.948 for the training set; an AUC of 0.884 (95% CI, 0.759–0.958), a sensitivity of 0.900, and a specificity of 0.790 for validation set 1; and an AUC of 0.876 (95% CI, 0.754–0.952), a sensitivity of 0.857, and a specificity of 0.837 for validation set 2 (Fig. 5a–c). The AUC values showed that the combined nomogram performed well in the assessment of tumor grade. The calibration curve showed that there was sufficient consistency between the nomogram-estimated grade and the actual observed in the three cohorts (Fig. 4b-d). Table 4 shows the details of the radiomic nomogram’s performance.

Fig. 5
figure 5

Accuracy evaluation of the radiomic nomogram. The accuracy of the radiomic nomogram to predict high-grade of BD-IPMNs was evaluated in the training cohort (a) (AUC = 0.903), external validation cohort 1 (b) (AUC = 0.884), and external validation cohort 2 (c) (AUC = 0.876)

Table 4 Predictive performance of the radiomic signature and radiomic Nomogram

DCA was used to reveal the utility of the combined nomogram for clinical decision making. The clinical utility of the corresponding strategies was demonstrated by the area under the decision curve (Fig. 6a–c). The area occupied by the combined nomogram (red) was larger than that of the radiomic signature (blue) alone, and was larger than those of the “all “(gray) or “none” (black) strategies in the training and the validation sets.

Fig. 6
figure 6

Clinical utility evaluation of the radiomic signature and the radiomic nomogram using DCA curves. Evaluation in the training cohort (a), external validation cohort 1 (b), and external validation cohort 2 (c). The y-axis shows the net benefit. The x-axis shows the threshold probability. The highest net benefit was gained using the radiomic nomogram (red line) compared with the radiomic signature (blue line), the treat-all strategy (gray line), and the treat-none strategy (horizontal black line)

Discussion

In the present study, we constructed a combined nomogram and investigated its ability to predict tumor pathological grade preoperatively in patients with BD-IPMNs. The radiomic signature was based on nine features, and was considered an effective method of preoperative tumor grade assessment, representing a non-invasive imaging biomarker. Combining the radiomics signature with clinical variables (CA19–9 and the MPD size), as a combined nomogram model, significantly improved the predictive performance. The repeatability and reliability of the developed prediction model was confirmed using independent datasets from other institutions.

Preoperative grade assessment of patients with BD-IPMNs is clinically important. According to the current consensus, most of the excised BD-IPMNs are low-grade diseases. The average incidence of invasive malignant tumors in BD-IPMNs is 17.7% (1–37%) [8]; therefore, although many patients might have dangerous signs, such as mural nodules and dilation of the main pancreatic duct, they can still be followed up for a long time without selective surgery. However, currently, the performance of the existing discrimination system is unsatisfactory, and most studies recommend immediate resection for patients with obstructive jaundice, enhanced mural nodules, and main duct dilation > 10 mm [30]. Similarly, there are no biomarkers to predict high-grade BD-IPMNs. Therefore, an advanced discriminative method with high sensitivity and specificity would provide valuable information to determine clinical strategies, resulting in patients with low-grade IPMNs avoiding unnecessary surgery for a benign disease. Of note, the goal in stratifying patients is to identify those that are at risk of developing or harboring invasive malignancy, and so even if patients with “benign” disease undergo resection this may be of treatment benefit so that they do not develop future carcinoma.

In the present study, we analyzed the predictive power of quantitative MRI features to assess the grade of BD-IPMNs. The results showed the potential value of radiomic features. We analyzed both high and low-order radiomic features identified in previous studies. The histogram parameter (low-order feature) is associated to single pixel characteristics, describing the distribution of voxel intensity via common and basic measures [31]. Hoffman et al. performed preoperative MRI texture analysis [32] and showed that intensity histogram-based statistical features and entropy from MRI images could predict the malignancy of BD-IPMNs. High order radiomic features, including GLCM features, mainly assess the spatial relationships among local neighboring pixels [33, 34]. Hanania et al. [13] identified 14 imaging biomarkers within GLCM features that predicted the histopathological grade within cyst contours. Tobaly et al. [35] developed a radiomic model mostly based on high order CT radiomic features, which showed high diagnostic performance in differentiating benign from malignant IPMNs. Our results showed that most of the selected features were high-order features, which was consistent with the results presented by previous studies. Therefore, we hypothesized that the tumor biology and heterogeneity were better represented by high-order features. Previous studies recognized the value of assessing pathological features among radiomic features, such as in nonfunctional pancreatic neuroendocrine tumors, soft-tissue masses, and rectal cancer [36,37,38]. However, it remains challenging to associate a single radiomic feature with the complex biological processes of tumors. Therefore, it is common to construct a multi-factor panel to estimate the results in a radiomic background. Multi-factor based radiomic methods are usually more suitable to describe the complex heterogeneity of BD-IPMNs. The results presented here indicated that our radiomic signature could satisfactorily discriminate low-grade and high-grade of BD-IPMNs in the patients in the training and validation cohorts.

According to previous studies of IPMNs, the symptoms, obstructive jaundice, presence of mural nodules, cyst size, age, and sex were also associated with pathology grade assessment [21, 39]. However, we failed to confirm these results in the present study. In addition, parameters like pain could be subjectively reported according to the different sensitivities of patients to tumor-related abdominal pain, which might also be associated with tumor infiltration into the viscera, blood vessels, or peripheral nerves. We observed that age and sex were not related to grade, which might be related to populations from different ethnicities. Obstructive jaundice was rare in the patients in our study, which could be related to the limited sample size.

The radiomic nomogram constructed in our study incorporated clinical factors (the size of MPD and CA19–9 levels), which might be useful to identify those patients at risk of malignancy and to select the best treatment.

The CA19–9 level is useful for the differential diagnosis of pancreatic carcinoma and benign pancreatic diseases, and increased CA19–9 levels are an independent predictor of malignant BD-IPMNs [20, 40,41,42]. However, in the present study, CA19–9 alone achieved an AUC of 0.762 to discriminate high-grade in patients training set. Thus, CA19–9 alone is not sufficient to evaluate the high-grade of BD-IPMNs accurately. Previous studies showed that main duct dilatation was a significant predictor of malignancy and suggested that the patients with main duct dilatation over 10 mm should undergo resection without further testing or calculation [39]. In the present study, we divided the MPD dilation into four degrees according to previous studies [10, 39], and found that patients with a larger MPD dilation were more likely to have a high-grade of BD-IPMNs.

The DCA curve showed that the radiomic nomogram was superior to the radiomic signature over a large threshold probability range, which indicated that the clinical parameters increased the incremental value for BD-IPMN grade assessment in the training and validation sets.

The established grade assessment model has obvious advantages. For clinicians seeking low-cost techniques to improve patient management, quantitative image analysis is attractive, because it comprises a noninvasive disease assessment at multiple time points and tumor sites compared with conventional diagnostic imaging. In the current study, the results showed that the AUC of the model was further improved to 0.903 by combining the nomogram with clinical characteristics and the radiomic signature, suggesting that advanced imaging technology should be considered together with well-known clinical factors. In addition, the pre-operative data used to establish the grade assessment model is easy to obtain, with less additional costs and provides good results.

The present study had certain limitations. First, because of the retrospective design of our study, there may be selection bias. Second, this study involved a relatively small sample, and the low-grade and high-grade groups were not further separated into patients with low and intermediate-grade dysplasia or high-grade dysplasia and patients with associated invasive carcinoma. However, our cohort was sufficient to build a reliable model. Third, the area of potential peritumoral infiltration were not included in this study, which may result in underestimation of tumor area. In future, studies are still needed for detailed analysis of features extracted from peritumoral area. However, as we focus on the visible tumor area, the results are still reliable. Our study had another limitation: Tumor area segmentation has to be performed manually by radiologists. Automatic segmentation will be more convenient and should become the standard. In the future, we will construct a fully automatic prediction model incorporating automatic segmentation of the pancreas and cyst areas.

Conclusions

The preoperative pathological grade of BD-IPMNs could be predicted effectively using the developed nomogram model combining the radiomic signature and tumor clinical characteristics. The predictive nomogram model represents an accurate and noninvasive assessment method for patients with BD-IPMNs before surgery, which will help clinicians to alter the treatment protocol for each patient and thus obtain improved clinical outcomes in the future.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

BD-IPMNs:

Branching type intraductal papillary mucinous neoplasms

MPD:

Main pancreatic duct

LGD:

Low-grade dysplasia

IGD:

Intermediate-grade dysplasia

HGD:

High-grade dysplasia

References

  1. Gaujoux S, Brennan MF, Gonen M, D'Angelica MI, DeMatteo R, Fong Y, et al. Cystic lesions of the pancreas: changes in the presentation and management of 1,424 patients at a single institution over a 15-year time period. J Am Coll Surg. 2011;212(4):590–600. https://doi.org/10.1016/j.jamcollsurg.

    Article  PubMed  Google Scholar 

  2. Sakorafas GH, Smyrniotis V, Reid-Lombardo KM, Sarr MG. Primary pancreatic cystic neoplasms revisited. Part III. Intraductal papillary mucinous neoplasms. Surg Oncol. 2011;20:e109–18. https://doi.org/10.1016/j.suronc.2011.01.004.

    Article  PubMed  Google Scholar 

  3. de Jong K, Nio CY, Hermans JJ, Dijkgraaf MG, Gouma DJ, van Eijck CHJ, et al. High prevalence of pancreatic cysts detected by screening magnetic resonance imaging examinations. Clin Gastroenterol Hepatol. 2010;8:806–11. https://doi.org/10.1016/j.cgh.2010.05.017.

    Article  PubMed  Google Scholar 

  4. Laffan TA, Horton KM, Klein AP, Berlanstein B, Siegelman SS, Kawamoto S, et al. Prevalence of unsuspected pancreatic cysts on MDCT. AJR Am J Roentgenol. 2008;191:802–7. https://doi.org/10.2214/AJR.07.3340.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Fernandez-del Castillo C, Warshaw AL. Current management of cystic neoplasms of the pancreas. Adv Surg. 2000;34:237–48.

    CAS  PubMed  Google Scholar 

  6. Tanaka M, Castillo CF-D, Kamisawa T, Jang JY, Levy P, Ohtsuka T, et al. Revisions of international consensus Fukuoka guidelines for the management of IPMN of the pancreas. Pancreatology. 2017;17:738–53. https://doi.org/10.1016/j.pan.2017.07.007.

    Article  PubMed  Google Scholar 

  7. Do RKG, Katz SS, Gollub MJ, Li J, LaFemina J, Zabor EC, et al. Interobserver agreement for detection of malignant features of intraductal papillary mucinous neoplasms of the pancreas on MDCT. AJR Am J Roentgenol. 2014;203:973–9. https://doi.org/10.2214/AJR.13.11490.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Tanaka M, Castillo C F-d, Adsay V, Chari S, Falconi M, Jang J-Y, et al. International consensus guidelines 2012 for the management of IPMN and MCN of the pancreas. Pancreatology. 2012;12:183–97. https://doi.org/10.1016/j.pan.2012.04.004.

    Article  PubMed  Google Scholar 

  9. Watanabe Y, Nishihara K, Niina Y, Abe Y, Amaike T, Kibe S, et al. Validity of the management strategy for intraductal papillary mucinous neoplasm advocated by the international consensus guidelines 2012: a retrospective review. Surg Today. 2016;46:1045–52. https://doi.org/10.1007/s00595-015-1292-2.

    Article  PubMed  Google Scholar 

  10. Yamada S, Fujii T, Murotani K, Kanda M, Sugimoto H, Nakayama G, et al. Comparison of the international consensus guidelines for predicting malignancy in intraductal papillary mucinous neoplasms. Surgery. 2016;159:878–84. https://doi.org/10.1016/j.surg.2015.08.042.

    Article  PubMed  Google Scholar 

  11. Chakraborty J, Midya A, Gazit L, Attiyeh M, Langdon-Embry L, Allen PJ, et al. CT radiomics to predict high-risk intraductal papillary mucinous neoplasms of the pancreas. Med Phys. 2018;45:5019–29. https://doi.org/10.1002/mp.13159.

    Article  PubMed  Google Scholar 

  12. Permuth JB, Choi J, Balarunathan Y, Kim J, Chen DT, Chen L, et al. Combining radiomic features with a miRNA classifier may improve prediction of malignant pathology for pancreatic intraductal papillary mucinous neoplasms. Oncotarget. 2016;7:85785–97. https://doi.org/10.18632/oncotarget.11768.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Hanania AN, Bantis LE, Feng Z, Wang H, Tamm EP, Katz MH, et al. Quantitative imaging to evaluate malignant potential of IPMNs. Oncotarget. 2016;7:85776–84. https://doi.org/10.18632/oncotarget.11769.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Avanzo M, Stancanello J, El Naqa I. Beyond imaging: the promise of radiomics. Phys Med. 2017;38:122–39. https://doi.org/10.1016/j.ejmp.2017.05.071.

    Article  PubMed  Google Scholar 

  15. Emaminejad N, Qian W, Guan Y, Tan M, Qiu Y, Liu H, et al. Fusion of Quantitative Image and Genomic Biomarkers to Improve Prognosis Assessment of Early Stage Lung Cancer Patients. IEEE Trans Biomed Eng. 2016;63:1034–43. https://doi.org/10.1109/TBME.2015.2477688.

    Article  PubMed  Google Scholar 

  16. Lloyd MC, Allam-Nandyala P, Purohit CN, Burke N, Coppola D, Bui MM, et al. Using image analysis as a tool for assessment of prognostic and predictive biomarkers for breast cancer: how reliable is it? J Pathol Inform. 2010;1:29. https://doi.org/10.4103/2153-3539.74186.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Cameron A, Khalvati F, Haider MA, Wong A. MAPS: a quantitative Radiomics approach for prostate Cancer detection. IEEE Trans Biomed Eng. 2016;63:1145–56. https://doi.org/10.1109/TBME.2015.2485779.

    Article  PubMed  Google Scholar 

  18. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, They Are Data. Radiology. 2016;278:563–77. https://doi.org/10.1148/radiol.2015151169.

    Article  PubMed  Google Scholar 

  19. He L, Li H, Dudley JA, Maloney TC, Brady SL, Somasundaram E, et al. Machine Learning Prediction of Liver Stiffness Using Clinical and T2-Weighted MRI Radiomic Data. AJR Am J Roentgenol. 2019;213:592–601. https://doi.org/10.2214/AJR.19.21082.

    Article  PubMed  Google Scholar 

  20. Hwang DW, Jang J-Y, Lim C-S, Lee SE, Yoon Y-S, Ahn YJ, et al. Determination of malignant and invasive predictors in branch duct type intraductal papillary mucinous neoplasms of the pancreas: a suggested scoring formula. J Korean Med Sci. 2011;26:740–6. https://doi.org/10.3346/jkms.2011.26.6.740.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Shimizu Y, Yamaue H, Maguchi H, Yamao K, Hirono S, Osanai M, et al. Validation of a nomogram for predicting the probability of carcinoma in patients with intraductal papillary mucinous neoplasm in 180 pancreatic resection patients at 3 high-volume centers. Pancreas. 2015;44:459–64. https://doi.org/10.1097/MPA.0000000000000269.

    Article  PubMed  Google Scholar 

  22. Correa-Gallego C, Do R, Lafemina J, Gonen M, D'Angelica MI, DeMatteo RP, et al. Predicting dysplasia and invasive carcinoma in intraductal papillary mucinous neoplasms of the pancreas: development of a preoperative nomogram. Ann Surg Oncol. 2013;20:4348–55. https://doi.org/10.1245/s10434-013-3207-z.

    Article  PubMed  Google Scholar 

  23. John D. Birkmeyer, Andrea E Siewers, Emily V a Finlayson, Therese a Stukel, F Lee Lucas, Ida Batista. et al. hospital volume and surgical mortality in the United States. N Engl J Med. 2002;346:1128–37. https://doi.org/10.1056/NEJMsa012337.

    Article  Google Scholar 

  24. Lee JE, Choi S-Y, Min JH, Yi BH, Lee MH, Kim SS, et al. Determining Malignant Potential of Intraductal Papillary Mucinous Neoplasm of the Pancreas: CT versus MRI by Using Revised 2017 International Consensus Guidelines. Radiology. 2019;293:134–43. https://doi.org/10.1148/radiol.2019190144.

    Article  PubMed  Google Scholar 

  25. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage. 2006;31:1116–28. https://doi.org/10.1016/j.neuroimage.2006.01.015.

    Article  PubMed  Google Scholar 

  26. Xia W, Chen Y, Zhang R, Yan Z, Zhou X, Zhang B, et al. Radiogenomics of hepatocellular carcinoma: multiregion analysis-based identification of prognostic imaging biomarkers by integrating gene data-a preliminary study. Phys Med Biol. 2018;63:035044. https://doi.org/10.1088/1361-6560/aaa609.

    Article  CAS  PubMed  Google Scholar 

  27. Vasquez MM, Hu C, Roe DJ, Chen Z, Halonen M, Guerra S, et al. Least absolute shrinkage and selection operator type methods for the identification of serum biomarkers of overweight and obesity: simulation and application. BMC Med Res Methodol. 2016;16:154. https://doi.org/10.1186/s12874-016-0254-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Kumamaru KK, Saboo SS, Aghayev A, Cai P, Quesada CG, George E, et al. CT pulmonary angiography-based scoring system to predict the prognosis of acute pulmonary embolism. J Cardiovasc Comput Tomogr. 2016;10:473–9. https://doi.org/10.1016/j.jcct.2016.08.007.

    Article  PubMed  Google Scholar 

  29. Mehta S, Bastero-Caballero RF, Sun Y, Zhu R, Murphy DK, Hardas B, et al. Performance of intraclass correlation coefficient (ICC) as a reliability index under various distributions in scale reliability studies. Stat Med. 2018;37:2734–52. https://doi.org/10.1002/sim.7679.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Tanaka M. Clinical management and surgical decision-making of IPMN of the pancreas. Methods Mol Biol. 2019;1882:9–22. https://doi.org/10.1007/978-1-4939-8879-2_2.

    Article  CAS  PubMed  Google Scholar 

  31. Vidić I, Egnell L, Jerome NP, Teruel JR, Sjøbakk TE, Østlie A, et al. Support vector machine for breast cancer classification using diffusion-weighted MRI histogram features: Preliminary study. J Magn Reson Imaging. 2018;47:1205–16. https://doi.org/10.1002/jmri.25873.

    Article  PubMed  Google Scholar 

  32. Hoffman DH, Ream JM, Hajdu CH, Rosenkrantz AB. Utility of whole-lesion ADC histogram metrics for assessing the malignant potential of pancreatic intraductal papillary mucinous neoplasms (IPMNs). Abdom Radiol (NY). 2017;42:1222–8. https://doi.org/10.1007/s00261-016-1001-7.

    Article  Google Scholar 

  33. Kazerooni AF, Nabil M, Khah HH, Alviri M, Heidari-Sooreshjaani M, Gity M, et al. ADC-derived spatial features can accurately classify adnexal lesions. J Magn Reson Imaging. 2018;47:1061–71. https://doi.org/10.1002/jmri.25854.

    Article  Google Scholar 

  34. Attiyeh MA, Chakraborty J, Gazit L, Langdon-Embry L, Gonen M, Balachandran VP, et al. Preoperative risk prediction for intraductal papillary mucinous neoplasms by quantitative CT image analysis. HPB (Oxford). 2019;21:212–8. https://doi.org/10.1016/j.hpb.2018.07.016.

    Article  Google Scholar 

  35. Tobaly D, Santinha J, Sartoris R, Dioguardi Burgio M, Matos C, Cros J, et al. CT-based Radiomics analysis to predict malignancy in patients with Intraductal papillary mucinous neoplasm (IPMN) of the pancreas. Cancers (Basel). 2020;12(11):3089. https://doi.org/10.3390/cancers12113089.

    Article  CAS  Google Scholar 

  36. Bian Y, Zhao Z, Jiang H, Fang X, Li J, Cao K, et al. Noncontrast Radiomics Approach for Predicting Grades of Nonfunctional Pancreatic Neuroendocrine Tumors. J Magn Reson Imaging. 2020. https://doi.org/10.1002/jmri.27176.

  37. Wang H, Zhang J, Bao S, Liu J, Hou F, Huang Y, et al. Preoperative MRI-based Radiomic machine-learning Nomogram may accurately distinguish between benign and malignant soft-tissue lesions: a two-center study. J Magn Reson Imaging. 2020. https://doi.org/10.1002/jmri.27111.

  38. Ma X, Shen F, Jia Y, Xia Y, Li Q, Jianping L, et al. MRI-based radiomics of rectal cancer: preoperative assessment of the pathological features. BMC Med Imaging. 2019;19:86. https://doi.org/10.1186/s12880-019-0392-7.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Jang J-Y, Park T, Lee S, Kim Y, Lee SY, Kim S-W, et al. Proposed Nomogram Predicting the Individual Risk of Malignancy in the Patients With Branch Duct Type Intraductal Papillary Mucinous Neoplasms of the Pancreas. Ann Surg. 2017;266:1062–8. https://doi.org/10.1097/SLA.0000000000001985.

    Article  PubMed  Google Scholar 

  40. Fritz S, Hackert T, Hinz U, Hartwig W, Büchler MW, Werner J, et al. Role of serum carbohydrate antigen 19-9 and carcinoembryonic antigen in distinguishing between benign and invasive intraductal papillary mucinous neoplasm of the pancreas. Br J Surg. 2011;98:104–10. https://doi.org/10.1002/bjs.7280.

    Article  CAS  PubMed  Google Scholar 

  41. Shin SH, Han DJ, Park KT, Kim YH, Park JB, Kim SC, et al. Validating a simple scoring system to predict malignancy and invasiveness of intraductal papillary mucinous neoplasms of the pancreas. World J Surg. 2010;34:776–83. https://doi.org/10.1007/s00268-010-0416-5.

    Article  PubMed  Google Scholar 

  42. Ingkakul T, Sadakari Y, Ienaga J, Satoh N, Takahata S, Tanaka M, et al. Predictors of the presence of concomitant invasive ductal carcinoma in intraductal papillary mucinous neoplasm of the pancreas. Ann Surg. 2010;251:70–5. https://doi.org/10.1097/SLA.0b013e3181c5ddc3.

    Article  PubMed  Google Scholar 

Download references

Acknowledgments

The authors thank the editors of ELIXIGEN for language editing.

Funding

This work was supported by the Key Research and Development Project of Zhejiang Province (2020C01058).

Author information

Authors and Affiliations

Authors

Contributions

SC: manuscript draft, conception and design, study selection, data extraction, review of final draft. TT: conception and design, study selection, data extraction, provide technical support, review of final draft. QS: collection and assembly of data, study selection, statistical analysis. YW: collection and assembly of data, data analysis and interpretation. ZS: conception and design of the study, review of final draft. WY: collection and assembly of data, study selection. XG: conception and design, review of final draft. Final approval of manuscript: All authors.

Corresponding author

Correspondence to Xiangyang Gong.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the ethics review board of the First Affiliated Hospital of Zhejiang University, the second Affiliated Hospital of Zhejiang University and the Zhejiang Provincial People’s Hospital, the need for informed patient consent for inclusion was waived.

Consent for publication

Publication is approved by all authors and by the responsible authorities where the work was carried out.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1.

The Supplementary Material for this article was available in additional file 1.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cui, S., Tang, T., Su, Q. et al. Radiomic nomogram based on MRI to predict grade of branching type intraductal papillary mucinous neoplasms of the pancreas: a multicenter study. Cancer Imaging 21, 26 (2021). https://doi.org/10.1186/s40644-021-00395-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40644-021-00395-6

Keywords