Multiphase computed tomography radiomics of pancreatic intraductal papillary mucinous neoplasms to predict malignancy

Stuart L Polk, Jung W Choi, Melissa J McGettigan, Trevor Rose, Abraham Ahmed, Jongphil Kim, Kun Jiang,Yoganand Balagurunathan, Jin Qi, Paola T Farah, Alisha Rathi, Jennifer B Permuth, Daniel Jeong

Abstract

Key words: Radiomics; Intraductal papillary mucinous neoplasm; Multiphase computed tomography; Pancreas; Oncology; Pancreatic cancer

INTRODUCTION

Intraductal papillary mucinous neoplasms (IPMNs) are mucin producing neoplasms that can arise from the pancreatic main duct or its side branches. These are considered non-invasive precursor lesions that can give rise to invasive pancreatic ductal adenocarcinoma. IPMNs are a subset of pancreatic cystic lesions that have been shown in histopathological studies to account for approximately half of resected pancreatic cysts[1-3]. IPMNs are typically detected incidentally on imaging, as pancreatic cystic lesions are reported to be present on up to 2.6% of all computed tomography (CT) studies[4,5]. Despite increased detection of IPMNs on CT, it remains challenging to accurately predict whether IPMNs harbor high grade or invasive disease which would warrant surgical resection. Conversely, there exists a significant risk of major or minor post-operative morbidity (40%) and post-operative mortality(2.3%) with surgical resection of pancreatic cystic lesions[6]. This underscores the implications at stake for accurately identifying an IPMN as malignant or benign - it can lead to either the removal of a potentially life-threatening lesion or the avoidance of unnecessary surgical intervention and its associated consequences and economic costs. While the American College of Radiology (ACR), American Gastroenterology Association (AGA) and International Consensus Guidelines (ICG) (also called the Fukuoka Consensus Guidelines, or Fukuoka Criteria) each provide dedicated guidelines for management of cystic pancreatic lesions, no single guideline provides high sensitivity in diagnosing malignant lesions without also causing an increase in false positives[7-9]. At many institutions including ours, the ICG serve as the basis for evaluating suspected IPMNs on imaging[10,11]. ICG includes standard, qualitative radiologic and clinical features categorized as either “worrisome features” or “highrisk stigmata”, which are used to classify an IPMN as benign (characterized by low or moderate grade pathology) or malignant (characterized by high grade or invasive pathology). Based on this classification using the ICG criteria, a patient will be designated for surveillance with follow-up in the case of a benign IPMN or further work up and/or surgical resection in the case of a malignant IPMN. CT imaging is important as it provides the basis for standard IPMN evaluation per ICG and allows triaging of patients into intervention or surveillance groups.

Several recent studies have evaluated the accuracy of the ICG 2012 criteria in characterizing IPMNs. Tsukagoshiet al[12]reported a diagnostic accuracy of 77% for predicting highvslow-grade dysplasia in IPMNs using the ICG criteria. Sugimotoet al[13]and Yamadaet al[14]attained an area under the curve (AUC) of 67% and 63%respectively for predicting malignantvsbenign IPMNs using the ICG criteria[13,14].Other studies have illustrated similar diagnostic accuracies in the specificity and positive predictive value (PPV) of the ICG criteria, with false positive rates of up to 47% and 66% in predicting malignancy in branch-duct IPMNs (BD-IPMNs)[7,15,16].Furthermore, there is a discrepancy between pre- and post-op diagnosis for nearly 1 out of 3 patients who have undergone pancreatic cyst surgery[17]. While the ICG serve as a valuable diagnostic method for characterizing IPMNs, these studies suggest there is still room for improvement. ICG criteria rely heavily on qualitative metrics to predict IPMN pathology. A quantitative radiological assessment, radiomics, may potentially improve diagnostic accuracy of radiological evaluation of IPMNs.

Radiomics involves extraction and analysis of quantitative features from diagnostic images to produce mineable data[18]. Radiomics features encompass size, shape,texture, heterogeneity, and signal intensity. Using dedicated software and predictive models, these imaging features (or radiophenotypes) are quantified into specific numeric values which can be correlated to clinicopathologic diagnoses and outcomes[18]. Radiomic features possess several advantages over standard qualitative radiologic features. They can represent quantitative and objective measures and reflect tumor heterogeneity and sub-regional environment[19,20]. In addition to reproducibility and stability, they also exhibit substantial prognostic value[21-25].Previous authors have shown that radiomics features can be more strongly linked to outcomes than prevailing clinical data[26-30]. CT is the more common imaging modality for performing radiomics analysis of pancreatic cysts and IPMNs for several reasons.First, CT and magnetic resonance imaging (MRI) display similar accuracy in characterizing pancreatic cysts as either benign or malignant[31,32]. CT also possesses widespread availability, high spatial and temporal resolution, and short scanning durations. Additionally, CT voxels utilize standard Hounsfield Units (HU) and not the relative signal intensity seen with conventional MRI sequences.

Several groups have shown that two dimensional (2D) CT radiomics can help predict malignant histology in IPMNs. Our team showed that a set of textural radiomic features could be used to identify malignant IPMNs with a diagnostic accuracy of 77% using venous phase single mid-lesion images for radiomic analysis[33].In a related study, Hananiaet al[34]identified gray level co-occurrence matrix biomarkers on arterial phase 2D images and achieved a diagnostic accuracy of 82% in predicting malignant or high-grade disease. Chakrabortyet al[35]used 2D radiomic texture features on venous phase images and reported diagnostic accuracy of 74% in identifying high risk IPMNs. The purpose of this proof-of-concept study is to expand on prior 2D work and to evaluate whole lesion radiomic features of IPMNs in predicting malignant histology on multiphase CT. We hypothesized that CT whole lesion radiomics features in conjunction with the ICG could improve the diagnostic accuracy of predicting IPMN histology.

MATERIALS AND METHODS

Fifty-one subjects with pancreatic lesions who had surgical resection at our tertiary cancer center with histology demonstrating intraductal papillary mucinous neoplasm were included in this IRB approved, Health Insurance Portability and Accountability Act (HIPAA) compliant retrospective cohort. The initial cohort was derived from a radiological records database search with matching pathology result availability.Subjects included had surgical resection between 2001 and 2016 and preoperative multiphase CT available. Subjects were excluded who had only noncontrast or single phase post contrast CT imaging or unavailable pathology reports. Subjects with surgical pathology demonstrating high grade dysplasia or invasive tumor were included in the "malignant" cohort while subjects with surgical pathology demonstrating moderate or low grade dysplasia were included in the "benign" cohort.Benign and malignant cohorts were composed of similar aged subjects with a similar proportion of males and females.

Region identification

Whole lesion semi-automated segmentation was performed on preoperative (within 3 mo preceding surgery) axial CT exams by a dedicated oncologic radiologist, using the HealthMyne software platform (Healthmyne, Madison, WI). CT exams were performed utilizing multiphase protocols and segmentations were performed on each slice of the lesion on all available noncontrast, arterial, and venous phase CTs by a board-certified radiologist (Figure 1). Thirty-nine relevant radiomic features listed in Table 1 were extracted from each segmented lesion and analyzed for each contrast phase.

CT protocols

CT exams were performed on different scanners: Siemens Sensation 16 (n= 29),Sensation 40 (n= 8), Sensation 64 (n= 7), SOMATOM Emotion (n= 1), SOMATOM Definition AS (n= 3), (Siemens Healthcare, Erlangen, Germany). CT exams from GE LightSpeed Pro 32 (n= 1) (GE Healthcare, Waukesha, WI) Philips Brilliance 16P (n=1) (Philips Medical Systems, Best, The Netherlands), and Toshiba Aquilion (n= 1)(Toshiba Medical Systems Co, Otawara, Japan) were also used for exams in this study.CT exams were 3 phase studies with noncontrast, arterial, and venous phases obtained. Arterial phase CTs were generally obtained with bolus triggering and an injection rate of 4 mL/s of Iopamidol 76% weight based dosing to achieve late arterial phase imaging approximately 35 s post injection (Bracco Diagnostics Inc., Monroe Township, NJ). Venous phase CTs achieved imaging approximately 60 s post injection. Weight based Iopamidol dosing generally ranged from 75 mL for patients below 55 kg, to 150 mL for patients above 110 kg with gradient increases every 5 kg.Field of view (FOV) ranged from (293-494) mm × (293-494) mm for each phase based on patient size. Matrix was 512 × 512 for each CT. Slice thickness was 4.2 ± 1.2 mm for noncontrast exams, and 3.0 ± 0.4 mm for arterial and venous exams. Noncontrast,arterial, and venous phase voxel volumes were 2.5 ± 0.8 mm3, 1.8 ± 0.5 mm3, and 1.8 ±0.5 mm3, respectively. Region of Interest (ROI) measurements were performed to measure Hounsfield units over the abdominal aorta and suprarenal inferior vena cava(IVC) on each contrast phase to ensure no differences between groups were specifically attributable to contrast phase variability. For contrast phase control purposes, manual regions of interest were drawn over the aorta and IVC and Hounsfield units were recorded to ensure there were no significant differences in contrast phase related variability between groups (Table 2).

Statistical analysis

Univariate analysis of the 39 radiomic features was performed for the pre-contrast,arterial, and venous phase radiomics data. Conventional quantitative and qualitative CT measurements were collected, which included ICG criteria. Categorical variables were compared with the chi-square test and continuous variables compared with the Mann WhitneyUtest. Logistic regression analysis was performed to discriminatepathologically-confirmed malignant and benign IPMN case groups using SAS/STAT version 14.3 (SAS Institute Inc., Cary, NC, United States). Multivariate analysis was performed on data normalized by standard deviation to generate a model including radiomic features and clinical variables. Log 2 transformation was performed on 13 radiomic features prior to model building to approximately conform to normality.Features withPvalues ＜ 0.1 were considered to build multivariable models.Backward elimination was performed with stay level of 0.05. The multivariate model's ability to discriminate was assessed by 5-fold cross validation. The statistical methods of this study were reviewed by our biostatician Jongphil Kim PhD. from H.Lee Moffitt Cancer Center.

Table 1 Radiomic features evaluated in the intraductal papillary mucinous neoplasm study cohort (n = 39)

RESULTS

The IPMN study cohort consisted of 51 patient subjects. Twenty-nine subjects (15 males, age 71.4 ± 8.6 years) comprised the malignant cohort, and 22 subjects (11 males,age 70.0 ± 7.1 years) comprised the benign cohort. Radiomic analysis showed 18/39 noncontrast (Table 3), 19/39 arterial phase (Table 4), and 21/39 venous phase (Table 5) features differentiated malignant from benign IPMNs (P＜ 0.05). The most statistically predictive features in univariate analysis were seen on the venous phase and these included log2 (GLCM_energy), log2 (energy_HU), log2 (GLCM_angular_second_moment), and GLCM_entropy.

Of the conventional CT measurements including ICG criteria (Table 6), main pancreatic duct diameter (P＜ 0.001), enhancing solid component (P= 0.005), and maximum cyst wall thickness (P＜ 0.001) were the specific features with the most significant difference between malignant and benign pathology groups. Multivariate analysisviabackward elimination, including the individual ICG 2017 criteria yielded two significant variables: thickened and enhancing cyst wall and enhancing mural nodule ＜ 5 mm with a combined ICG model AUC (95%CI) of 0.817 (0.709-0.926).

Multivariable analysisviabackward elimination for post contrast radiomics achieved an AUC (95%CI) of 0.87 (0.767-0.974) for a model including arterial phase radiomics features and 0.834 (0.716-0.953) for a model including venous phase radiomics features. Combined multivariable analysis including ICG criteria and radiomic results are shown in Table 7. Of note, the multivariable model including conventional variables and arterial phase radiomics features achieved an AUC(95%CI) of 0.93 (0.85-1.0) with a 5-fold cross validation AUC of 0.90.

No significant differences were noted between malignant and benign cohorts with voxel volumes on noncontrast, arterial, and venous phases (P= 0.5, 0.3, and 0.3)respectively. No significant differences were noted between aortic and IVC ROI Hounsfield unit measurements between groups on each phase (Table 2).

Figure 1 Low and high grade intraductal papillary mucinous neoplasms on computed tomography. A: Axial post contrast arterial phase computed tomography(CT) image magnified at the pancreatic head shows a hypodense well defined pathology proven intraductal papillary mucinous neoplasm (IPMN) with low grade dysplasia (blue arrow). Arterial gray level co-occurrence matrix (GLCM) entropy was 9.31; B: Axial post contrast arterial phase CT image in a different patient with pathology proven IPMN with high grade dysplasia (yellow arrow). Arterial GLCM entropy was 11.52. These lesions could be described similarly using qualitative terminology, however in this study, quantitative radiomics features improved prediction of malignant pathology compared to conventional measures alone.

DISCUSSION

In this retrospective analysis of CT data from patients who underwent resection for IPMN, we demonstrated that a model combining radiomic features and ICG criteria improved diagnostic accuracy compared to ICG criteria alone. Our model yielded an AUC of 0.90 in predicting pathologic status of IPMNs compared to an AUC of 0.82 using the ICG criteria.

Our work is concordant with previous studies by our team and others, attaining similar values of AUC for predicting IPMN malignancy using CT radiomic features[33,34]. Permuthet al[33]integrated worrisome radiologic features from the ICG,radiomic features, and miRNA expression levels into a logistic regression model that demonstrated an AUC = 0.93 for predicting IPMN malignant pathology. Hananiaet al[34]used a cross-validated panel consisting of 10 top-performing markers to attain a linear regression model yielding an AUC = 0.96 for differentiating low-grade and high-grade IPMNs. In line with our study, each of these studies independently identified that the strongest radiomic features for differentiating high-grade and lowgrade IPMNs were textural. In a more recent study, Chakrabortyet al[35]combined radiomic features with five clinical variables associated with IPMN risk to design a predictive model for discriminating between low-risk and high-risk BD-IPMNs,which achieved an AUC = 0.81. The current study adds to these previous works by employing 3D whole lesion analysis as opposed to 2D mid-lesion analysis, and we also evaluated multiphase CT for our radiomic analysis as opposed to solely venous or arterial phase CT. We showed that unique data is contained within radiomics features of each contrast phase. While several individual venous phase features produced the highest univariate AUCs, multivariable analysis showed the arterial phase model had higher multivariable AUC. Our results suggest there is diagnostic benefit in multiphase radiomics evaluation of IPMNs compared solely to venous phase evaluation.

Radiomics allowed us to quantify the appearance of a lesion, which can sometimes be challenging to appreciate qualitatively. We quantified this appearance using entropy applied to the gray level co-occurrence matrix, which is an imaging biomarker that can be used for characterizing tumor phenotype and heterogeneity.Tumors have extensive genetic and phenotypic variation, which manifests as heterogeneity within the same tumor, across metastases from the same tumor, and between patients[36-41]. Entropy reflects tumor heterogeneity and is associated with the extent of differential gene expression within a tumor as well as tumor metabolism,tumor stage, patient prognosis, and treatment response[36-41]. Higher tumor heterogeneity has been associated with poor prognosis, high-grade dysplasia, and resistance to anticancer treatment[42-46]. Tumor heterogeneity can be accurately estimated through genomic analysis, which typically requires invasive and risky tumor biopsies limited to the sampling site, and may not be representative of the entire tumor[42-46]. Radiomics can provide an alternative noninvasive and effective approach to measuring tumor heterogeneity.

In our study, IPMNs with pathology of high-grade dysplasia or invasive diseasehad significantly higher gray level co-occurrence matrix (GLCM) features on CT including energy and entropy. GLCM evaluates the spatial distribution of voxels based on their gray levels or in the case of CT, Hounsfield units. In GLCM, a matrix is constructed quantifying how often voxels with each HU value neighbor other voxels with similar HU. Various calculations can be performed on this matrix which ultimately describes the texture of the voxels within a region. Increased entropy as seen in our malignant cohort, suggests a higher degree of underlying tumor heterogeneity and is therefore associated with higher grade tumors. Qualitatively,higher entropy at a macroscopic level correlates with a “salt and pepper” appearance on CT. There are also underlying histologic implications as seen in (Figure 2), where a higher grade IPMN is shown to have more histologic irregularities compared to a lower grade IPMN. The measurement of entropy with quantitative imaging provides a way to appraise the degree of tumor heterogeneity in pancreatic cancer precursors,and elucidates information in the image which would be difficult to describe using only ICG and conventional imaging descriptors.

Table 2 Computed tomography Hounsfield units measured in the intraductal papillary mucinous neoplasm study cohort by contrast phase

Despite the benefit that radiomics presents to diagnostic medicine, there are challenges concerning methodology and its implementation into clinical practice that must be addressed. First, CT image fields of view are obtained using different settings based on patient body habitus. This can create voxel sizes that vary between patients,which could affect texture features. While digital imaging and communications in medicine (DICOM) datasets can be normalized for pixel size homogeneity, this was not performed in our series. However, we calculated voxel volumes for each CT based on FOV, matrix, and slice thickness data for our cohort and saw no significant differences between groups related to voxel size inhomogeneity. Additionally,multiphase post-contrast abdominal CT involves imaging at standardized times after contrast injection; however, imaging presentation can be sensitive to a patient's cardiac output[47]. The varied times of arrival of contrast to a lesion may slightly alter the post contrast Hounsfield Units, and therefore the quantification of radiomic features in that lesion. We performed standardized ROI Hounsfield measurements of the aorta and IVC on each of the 3 contrast phases and showed there was no significant contrast timing difference affecting our cohorts in each phase.

It is also important to normalize radiomic variables to a uniform scale when working with models. We performed this for multivariate analysis by normalizing feature values by standard deviations so that all variables are on a constant scale similar to previous authors[48]. Additionally, logarithmic standardization may be helpful to reduce the effect of statistical outliers. Degrees of change within each radiomic feature are relative to the group's composition and there are no wellestablished normal ranges for radiomic values within pancreatic IPMNs.Furthermore, radiomic feature nomenclature can vary depending on techniques and available software, highlighting the need for standard operating procedures and universal language in cases of collaboration between different practices, providers,and institutions[49].

Limitations of this proof-of-concept study include its retrospective study design and small sample size. We acknowledge that the main criteria for inclusion were a pathologic diagnosis obtained through surgical resection, which makes the study vulnerable to selection bias. In radiomic studies involving many features and/or biomarkers, it is more ideal to have a large sample size for logistic regression analysis.Also, performing radiomics on post-contrast CT is sensitive to slight changes intiming, causing variability in the quality of the images generated. To a similar point,since it is critical that CT scanner settings be homogenous for radiomics analysis,another limitation to this study is heterogeneity in CT scanner type and settings, as there is known radiomics variability with different scanner settings[49,50]. Our study included CT scans across a relatively long period involving different available CT technologies, whereas a future evaluation could be limited to a single scanner or provide inter-scanner validation. Furthermore, our study included both main-duct and branch-duct IPMNs. Main-duct IPMNs possess a much higher risk of malignant progression than branch-duct IPMNs, and are typically designated directly for surgery regardless of diagnostic evaluation of malignant status[51]. Therefore, future larger prospective radiomics studies evaluating particularly branch-duct IPMNs may have greater implications on management.

Table 3 Noncontrast phase radiomics features (P ＜ 0.05) for the intraductal papillary mucinous neoplasm study cohort

In summary, this study demonstrates the potential utility of whole lesion 3D radiomics as a supplementary method for predicting IPMN malignancy. Our findings suggest that utilizing radiomic features in combination with ICG criteria to evaluate IPMN pathology on CT can improve accuracy in prediction of high- or low-grade tumor histology. Additionally, multi-phase post contrast CT radiomic analysis may have diagnostic benefits over single venous phase CT analysis. We have just started to conduct a larger multi-center radiomic study that will have a training and test cohort and an independent prospective validation cohort and will integrate blood-based biomarker data to allow earlier and more accurate diagnosis of IPMNs of concern. As a result, patients may receive earlier surgical intervention when warranted and avoid unnecessary interventions in the cases of low risk disease.

Table 4 Arterial phase radiomics features (P ＜ 0.05) for the intraductal papillary mucinous neoplasm study cohort

Table 5 Venous phase radiomics features (P ＜ 0.05) for the intraductal papillary mucinous neoplasm study cohort

Table 6 Distribution of lnternational Consensus Guidelines criteria in benign vs malignant intraductal papillary mucinous neoplasm cases

Table 7 Multivariable model analysis

Figure 2 Histology of low grade and high grade intraductal papillary mucinous neoplasm lesions are shown. A: Low grade intraductal papillary mucinous neoplasm (IPMN) hematoxylin and eosin stain (HE) at 20 × magnification; B: Low grade IPMN HE at 100 × magnification. Low grade IPMN is represented by a single layer of mucinous epithelium containing columnar, palisading nuclei and abundant mucinous cytoplasm, with low nuclei/cytoplasmic ratio, minimal cytologic atypia and an absence of obvious mitosis; C: High grade IPMN HE at 20 × magnification; D: High grade IPMN HE at 100 × magnification. In contrast to low grade, high grade IPMN frequently demonstrates decreased cytoplasmic mucin contents, loss of nuclear polarity with nuclear overlapping, higher nuclei/cytoplasmic ratio, significant cytologic atypia and easily identifiable mitosis.

ARTICLE HIGHLIGHTS

This study demonstrates the potential utility of whole lesion radiomics on multiphase CT as a supplementary method for predicting malignancy in IPMN. A larger prospective multi-center trial will be beneficial to further define clinical utility and to improve risk stratification in IPMN.