ORIGINAL RESEARCH article

Front. Physiol., 21 February 2020

Sec. Gastrointestinal Sciences

Volume 11 - 2020 | https://doi.org/10.3389/fphys.2020.00008

Multiparametric Mapping Magnetic Resonance Imaging of Pancreatic Disease

  • 1. Department of Radiology, Beijing Chaoyang Hospital, Capital Medical University, Beijing, China

  • 2. Division of Digestive and Liver Diseases, Cedars-Sinai Medical Center, Los Angeles, CA, United States

  • 3. Department of Bioengineering, University of California, Los Angeles, Los Angeles, CA, United States

  • 4. Cedars-Sinai Biomedical Imaging Research Institute, Los Angeles, CA, United States

  • 5. Department of Nuclear Science and Engineering, Siemens Healthineers, Princeton, NJ, United States

  • 6. Department of Gastrointestinal Malignancies, Cedars-Sinai Medical Center, Los Angeles, CA, United States

Abstract

Background:

Current magnetic resonance imaging (MRI) of pancreatic disease is qualitative in nature. Quantitative imaging offers several advantages, including increased reproducibility and sensitivity to detect mild or diffuse disease. The role of multiparametric mapping MRI in characterizing various tissue types in pancreatic disease such as chronic pancreatitis (CP) and pancreatic ductal adenocarcinoma (PDAC) has rarely been evaluated.

Purpose:

To evaluate the feasibility of multiparametric mapping [T1, T2, and apparent diffusion coefficient (ADC)] in defining tissue characteristics that occur in CP and PDAC to improve disease diagnosis.

Materials and Methods::

Pancreatic MRI was performed in 17 patients with PDAC undergoing therapy, 7 patients with CP, and 29 healthy volunteers with no pancreatic disease. T1 modified Look-Locker Inversion Recovery (T1 MOLLI), T2-prepared gradient-echo, and multi-slice single-shot echo-planar diffusion weighted imaging (SS-EPI DWI) sequences were used for data acquisition. Regions of interest (ROIs) of pancreas in PDAC, CP, and control subjects were outlined by an experienced radiologist. One-way analysis of variance (ANOVA) was used to compare the difference between groups and regions of the pancreas, and Tukey tests were used for multiple comparison testing within groups. Receiver operator characteristic (ROC) curves were analyzed, and the areas under the curves (AUCs) were calculated using single parameter and combined parameters, respectively.

Results:

T1, T2, and ADC values of the entire pancreas among PDAC, CP, and control subjects; and between upstream and downstream portions of the pancreas in PDAC patients were all significantly different (p < 0.05). The AUC values were 0.90 for T1, 0.55 for T2, and 0.71 for ADC for independent prediction of PDAC. By combining T1, T2, and ADC, the AUC value was 0.94 (sensitivity 91.54%, specificity 85.81%, 95% CI: 0.92–0.96), which yielded higher accuracy than any one parameter only (p < 0.001).

Conclusion:

Multiparametric mapping MRI is feasible for the evaluation of the differences between PDAC, CP, and normal pancreas tissues. The combination of multiple parameters of T1, T2, and ADC provides a higher accuracy than any single parameter alone in tissue characterization of the pancreas.

Introduction

Pancreatic cancer is the most common pancreatic malignant neoplasm and is the third most common cause of cancer-related deaths from NCI data in 2018. In the United States, pancreatic cancer accounts for about 3% of all cancers in the United States and about 7% of all cancer deaths (American Cancer Society, 2018; Siegel et al., 2018). PDAC is the major subtype of exocrine tumor and constitutes more than 90% of all pancreatic malignancies. Because of the tumor’s unique microenvironment and aggressive nature, PDAC has a relatively poor response to the conventional systemic chemotherapy and poor prognosis, with an overall 5-year survival rate of 8.5% (Yamamoto et al., 2015). Studies found that CP has markedly increased risk for pancreatic cancer (Hao et al., 2017). Five years after diagnosis, CP has a nearly eightfold risk for pancreatic cancer (Kirkegard et al., 2017). Inflammation participates in the development of tumor initiation, progression, treatment response, metastasis, and prognosis (Shi and Xue, 2019). The typical histopathologic features of CP contain acinar cell atrophy, pancreatic fibrosis, leukocyte infiltration, fatty replacement, and distorted and blocked ducts. These findings in CP and PDAC suggest that there are similar radiologic appearances. Further, the upstream pancreas (toward the tail end of the tumor) can have changes of CP due to duct obstruction.

Contrast-enhanced CT was used in PDAC detection, staging, and evaluation of prognosis. A previous study also found that CT radiomics could predict PDAC SMAD4 status and tumor stromal content (Attiyeh et al., 2019). X-ray radiation and iodine allergy are the major risks of CT. Endoscopic ultrasonography (EUS) is the most sensitive non-operative imaging method for the detection of pancreatic cancer and showed to be superior to CT (Luz et al., 2014; Singh and Faulx, 2016). EUS-guided fine needle aspiration (FNB) can achieve cytological information. However, it is invasive and highly operator dependent.

Several studies aimed to quantitatively differentiate PDAC from CP or autoimmune pancreatitis using MR techniques (Park et al., 2009; Fukukura et al., 2012; Yin et al., 2017). DWI plays an important role in the identification of PDAC lesions from the background of pancreatic parenchyma (Ichikawa et al., 2007; Wang et al., 2011; Fukukura et al., 2012; Hayano et al., 2016). With high b-value DWI, PDAC lesions can be reliably detected as an increased focal hyper-intensity area (Fukukura et al., 2012).

Hecht et al. (2017) found that the ADC value was significantly lower in tumors with dense fibrosis and may serve as a biomarker of fibrosis architecture. Choi et al. (2016) found that DWI with ADC value was a promising method to differentiate PDAC from mass-forming autoimmune pancreatitis. However, the ADC value only reflects one aspect of the differences between PDAC and CP vs. the normal control pancreas.

T1 and T2 relaxation times are valuable as quantitative parameters to characterize different tissues, especially in myocardial and liver diseases (Apprich et al., 2012; Kali et al., 2015; Blystad et al., 2017; Chen et al., 2017; Vietti Violi et al., 2019). T1 mapping was shown to improve the diagnosis of myocarditis, infarction, iron overload, and amyloidosis (Guo et al., 2009; Karamitsos et al., 2013; Kali et al., 2015). Previous studies found that T1 mapping combined Gd–EOB–DTPA-enhanced MRI can be used to predict the pathologic grading of hepatocellular carcinoma (Cieszanowski et al., 2012; Banerjee et al., 2014; Chen et al., 2017). Recently, Banerjee et al. (2014) and Cassinotto et al. (2015) found that T1 values are strongly correlated with liver fibrosis and liver biopsy in a population of 79 patients. Cieszanowski et al. (2012) demonstrated significantly higher sensitivity and accuracy of T2 relaxation times than ADC values (99.0 and 89.3% vs. 79.0 and 80.9%, respectively) for diagnosing hepatic malignancy. Some studies focused on using the T1 value in the diagnosis and classification of CP and found that T1 could provide quantitative metrics for determining the presence and severity of acinar cell loss and aid in the diagnosis of CP (Tirkes et al., 2017, 2018; Wang et al., 2018). In an animal model of pancreatic cancer, Yin et al. (2017) found that multiparametric MRI was able to characterize pancreatic masses, suggesting that T1, T2, and ADC mapping may have a direct clinical application in patients with PDAC.

In this study, we evaluate the feasibility of non-contrast multiparametric mapping (T1, T2, and ADC) in defining tissue characteristics that occur in CP and PDAC needed for advances in specific diagnosis. The pathological changes in these diseases are complex including changes in cellular density, blood supply, fibrosis, edema, and inflammation. Our pilot study was designed to determine if there are quantitative methods, which can provide specific biomarkers using our novel method non-contract MRI methods based on the hypothesis that a combination of quantitative measures of relaxation time (T1 and T2 values) and ADC values will provide biomarkers that distinguish PDAC, CP, and normal control pancreas.

Materials and Methods

Study Population

During the 18-month study period from October, 2017 to April, 2018, patients with PDAC, CP, and normal control pancreas were recruited into this study. All PDAC patients were confirmed by histopathology using tissue obtained by endoscopic ultrasonography (EUS)-guided FNA. The PDAC patients were undergoing neo-adjuvant chemotherapy at the time of the research MRI.

All the patients with CP had diagnosis established by magnetic resonance cholangiopancreatography (MRCP) or endoscopic retrograde cholangiopancreatography (ERCP) with Cambridge classification for CP (Freeny, 1989; Nattermann et al., 1993; Sahai et al., 1998; Tirkes et al., 2018).

The normal control group had no history of acute pancreatitis, CP, diabetes, pancreatic surgery, and no family history of cancer. Patients with pancreatic cystic lesions, benign tumors, or marked pancreatic atrophy or fat degeneration on MRI images were excluded. This prospective study was approved by the local institutional review board. Written informed consents were obtained from all participants.

Magnetic Resonance Imaging Technique

All subjects were scanned on a 3.0T MR scanner (Biograph mMR, Siemens Healthcare Sector, Erlangen, Germany) with an 18-channel phase array surface coil and were placed head-first, supine position in the magnet. Conventional qualitative sequences and non-contrast quantitative mapping sequences were run for each subject.

Conventional qualitative sequences were transversal T1-weighted three-dimensional volumetric interpolated breath-hold examination (VIBE) with Dixon fat saturation (T1-VIBE-DIXON), and T2 HASTE in transverse, coronal, and sagittal orientations. The parameters of T1-VIBE-DIXON are repetition time (TR) = 4.15 ms; echo time (TE) = 1.39/2.65 ms; flip angle = 9°; field of view (FOV) = 247 × 330 mm; acquisition matrix = 320 × 180; echo train length (ETL) = 2; slice thickness = 3 mm; iPAT acceleration factor = 3. T2 HASTE was performed with the following acquisition parameters: TR = 1,000 ms; TE = 99 ms; flip angle = 105°; FOV = 226 × 330 mm; matrix = 256 × 176; ETL = 109; slice thickness = 5 mm; slice gap = 1 mm; iPAT acceleration factor = 2. MRCP was performed with the following parameters: TR = 8,903 ms; TE = 701 ms; flip angle = 100°; FOV = 300 × 300 mm; acquisition matrix = 384 × 384; ETL = 2; slice thickness = 1 mm; iPAT acceleration factor = 2.

Non-contrast quantitative MRI sequences consisted of 2D MOLLI Trufi with motion correction for T1 mapping, 2D T2-prepared FLASH for T2 mapping, and multi-slice single-shot echo-planar imaging (SS-EPI) for DWI and ADC mapping. T1 MOLLI was acquired using a three-point tool to localize the slice with a single 10-s breath-hold and the following parameters: acquisition window = 280.5 ms; TE = 1.12 ms; echo spacing = 2.7 ms, simulated R–R interval = 1,000 ms; flip angle = 35°; FOV = 390 × 390 mm; matrix = 192 × 144; iPAT acceleration factor = 2. T2-prepared FLASH was acquired, and the slices of T1 MOLLI were copied with a single 10-s breath and the following parameters: acquisition window = 207.4 ms; TR = 3.15 ms, TE = 1.32 ms; echo spacing = 3.1 ms, simulated R–R interval = 1,000 ms, duration of T2 preparations = 0, 30, and 55 ms; flip angle, 12°; FOV, 390 × 390 mm; matrix, 192 × 144; iPAT acceleration factor = 2. For T1 MOLLI and T2-prepared FLASH, three slices located at, above, and below the tumor area were acquired. DWI covered from the dome of the diaphragm to the lower edge of the kidney using the following parameters: b-values = 50, 400, and 800 s/mm2, TR = 4,500 ms; TE = 47 ms; flip angle, 12°; FOV, 390 × 390 mm; matrix, 192 × 144; ETL = 45; slice thickness = 5 mm; slice gap = 1 mm; iPAT acceleration factor = 2, number of averages = 2 (b = 50 s/mm2), 4 (b = 400 s/mm2), and 6 (b = 800 s/mm2). The parameters of all sequences are listed in Table 1.

TABLE 1

ParametersT1-VIBE- DIXONT2 HASTEMRCPT1 mapping (ms)T2 mapping (ms)DWI
Slice thickness (mm)351556
Gap (mm)0105.55.51
# Slices/partitions7246803350
Repetition time (ms)4.151,0008,9032.73.154,500
Echo time (ms)1.39/2.65 (OP/IP)997011.121.3247
Echo train length2109245
Acquisition matrix320 × 180256 × 176384 × 384192 × 126192 × 144172 × 132
Flip angle (degree)9105100351290
FOV (mm)247 × 330226 × 330300 × 300390 × 390390 × 390306 × 399
iPAT acceleration factor322222
Scan time18-s single breath-hold22-s free breathing10-min free breathing10 s × 3 breath-hold10 s × 3 breath-hold3:50 min free breathing
Additional informationRespiratory triggerT2 prepared = 0, 30, and 55 msb = 50 (2 average), 300 (4 average), 800 (6 average)

Parameters of all MRI sequences used in the study.

Imaging Analysis

T1, T2, and ADC maps of a whole pancreas with PDAC, CP, and normal control pancreas were analyzed and measured along the margin of the pancreas by one radiologist with 20 years of experience in abdominal MRI (LW). We used pixel-wise methods to obtain the mean value of the whole pancreas and avoided areas of necrosis and blood vessels. Tumors were localized with reference to other sequences such as T1-VIBE DIXON, T2 HASTE, DWI images, or the previous contrast-enhanced CT, MRI, or PET/CT within 1 week. In addition, the freehand ROIs (regions of interest) were located within the tumor margin avoiding areas of necrosis. The mean values of tumor, upstream and downstream pancreas (Figures 1, 2), CP (Figure 3), normal pancreatic head, body, and tail (Figure 4) were also obtained based on pixel-wise methods of three slices on T1 and T2 mapping and all the slices on the ADC maps. The tumor’s size was defined as the largest diameter in axial images according to the RECIST 1.1 criteria.

FIGURE 1

FIGURE 2

FIGURE 3

FIGURE 4

Statistical Analysis

All the data were analyzed on SPSS v22.0 (Armonk, NY, United States; IBM Corp) and MedCalc (MedCalc for Windows, version 16.2.0.0, Mariakerke, Belgium). Data from the regions of interest were tabulated in a Microsoft Excel worksheet (Microsoft Corporation, Seattle, WA, United States). Descriptive statistics (mean values and the standard deviation) were obtained for the whole pancreas with PDAC, CP, and normal control pancreas, for PDAC mass, upstream and downstream pancreas, and for the head, body, and tail of the normal control pancreas. ANOVA tests were used to compare quantitative parameters between groups, and post hoc Tukey tests were used for multiple comparison testing within groups. The value of p < 0.05 was considered significant. ROC curve analysis and logistic regression was performed to evaluate the sensitivity and specificity of single parameters and each combination of multiple parameters in the prediction of PDAC after chemotherapy, and the AUCs were calculated. The value of multiparametric mapping over single parameter mapping was determined in MedCalc, using the Z-statistic to compare the AUCs of the different ROC curves.

Results

Demographics

A total of 53 imaging studies were included (17 patients with PDAC, 7 patients with CP, and 29 volunteers with normal healthy pancreas). Among 17 patients with PDAC, 9 patients were male, and 8 patients were female. The mean age was 65 years old with a range of 46–80 years.

All patients with PDAC were confirmed by FNA and/or surgical specimen. Six lesions were moderately differentiated adenocarcinomas, three lesions were poorly differentiated adenocarcinomas, and the other seven lesions were read as invasive ductal adenocarcinoma, but the grade was not reported.

Based on expert radiologist review of MR images (LW), nine (52.9%) PDACs were located in the pancreatic head, three (17.6%) in the pancreatic neck, two (11.8%) PDACs in the pancreatic body, and three (17.6%) in the pancreatic tail. The median size for all tumors was 3.2 cm, with a range from 1.3 to 6.7 cm.

Among the seven patients (four males and three females) with CP, the mean age was 53 years old with a range of 30–72 years. The mean width of MPD was 4.0 mm with a range from 1.6 to 6.7 mm.

Among the 29 volunteers with normal healthy pancreas, 13 were male and 16 were female. The mean age was 50 years old with a range of 20–72 years.

Characteristics, Quantification, and Comparison of T1, T2, and ADC Maps

For the whole pancreas with PDAC, CP, and normal control pancreas, T1, T2, and ADC values are listed in Table 2 (Figure 5). Significant differences were found between the three groups (p < 0.001). The whole pancreas with PDAC showed the highest T1 value compared with the CP and normal control pancreas with significant differences found (p < 0.001). T1 values for CP were significantly higher than those for the normal pancreas (p < 0.001). Significant differences were found when comparing T2 and ADC values of the whole pancreas with PDAC to normal pancreas (p < 0.001). In addition, T2 and ADC values of CP had significant differences compared to the normal pancreas (p < 0.001 and p = 0.003). However, no significant differences were found in the whole pancreas in T2 and ADC between PDAC and CP (p = 0.053 and p = 0.171).

TABLE 2

T1 value (mean ± SD) (ms)T2 value (mean ± SD) (ms)ADC value (mean ± SD) (×10–3 mm2/s)
Whole pancreas with PDAC1,675.6 ± 238.163.7 ± 6.11.519 ± 0.189
Chronic pancreatitis1,324.0 ± 222.957.7 ± 10.71.345 ± 0.174
Normal control pancreas860 ± 70.948.0 ± 2.91.127 ± 0.158
p-Value*<0.001<0.001<0.001

T1, T2, and ADC values of the whole pancreas with PDAC, chronic pancreatitis (CP), and normal control pancreas.

*p-Value is the difference between three groups.

FIGURE 5

In patients with PDAC, the T1, T2, and ADC values of mass, upstream and downstream pancreas are listed in Table 3 (Figure 6). The PDAC mass had the highest T1, T2, and ADC values compared with those of the upstream and downstream pancreas. The T1 values showed a significant difference (p < 0.001) when comparing the PDAC tumor with the downstream pancreas. The downstream pancreas showed statistically significant differences with the upstream pancreas (p < 0.001). No significant differences were found between the PDAC and upstream pancreas (p > 0.05). The T2 values were statistically significant (p = 0.029) between the PDAC and downstream pancreas, but no significant differences were observed between the PDAC and upstream pancreas (p = 0.732), and between the upstream and downstream pancreas (p = 0.175). The ADC values were significantly different when comparing PDAC with the upstream (p = 0.048) and downstream pancreas (p = 0.003). However, no significant difference was found between the upstream and downstream pancreas (p = 0.471).

TABLE 3

T1 value (mean ± SD) (ms)T2 value (mean ± SD) (ms)ADC value (mean ± SD) (× 10–3 mm2/s)
PDAC1,816.5 ± 208.564.9 ± 7.61.525 ± 0.243
Downstream pancreas1,133.1 ± 225.754.9 ± 6.91.231 ± 0.168
Upstream pancreas1,598.5 ± 292.061.0 ± 9.21.355 ± 0.194
p-Value*<0.0010.0370.003

T1, T2, and ADC values of the PDAC, downstream, and upstream.

*p-Value is the difference between three groups.

FIGURE 6

In normal pancreas, the head, body, and tail showed relatively homogenous appearances (Figure 7). The T1, T2, and ADC values are listed in Table 4. No significant differences in T1, T2, and ADC values were found between the normal pancreatic head, body, and tail (p > 0.05).

FIGURE 7

TABLE 4

T1 value (mean ± SD) (ms)T2 value (mean ± SD) (ms)ADC value (mean ± SD) (×10–3 mm2/s)
Normal pancreatic head846.3 ± 74.647.5 ± 2.91.146 ± 0.177
Normal pancreatic body854.6 ± 85.648.1 ± 4.21.129 ± 0.172
Normal pancreatic tail870.0 ± 83.247.9 ± 3.51.108 ± 0.184
p-Value*0.5630.8320.802

T1, T2, and ADC values of normal pancreatic head, body, and tail.

*p-Value is the difference between three groups.

Comparisons of the parameters between PDAC mass and non-tumor parts of the pancreas in patients with PDAC, CP, and normal pancreas were performed (Table 5). T1 values were significantly different (p < 0.001) when comparing PDAC with non-tumor parts of the pancreas, CP, and normal pancreas. T2 and ADC value differences were statistically significant (p < 0.05) when comparing PDAC with non-tumor pancreas and normal pancreas. ADC values also shows statistically significant differences (p < 0.05) when comparing non-tumor pancreas and CP with normal pancreas. However, no significant differences between CP and non-tumor pancreas were found in T1, T2, and ADC values (p > 0.05). T2 and ADC values also showed no significant differences between PDAC and CP (p > 0.05).

TABLE 5

Comparison groupsT1-valueT2-valueADC
PDAC vs. non-tumor pancreas<0.0010.0260.003
PDAC vs. chronic pancreatitis<0.0010.6240.117
PDAC vs. normal control pancreas<0.001<0.001<0.001
Non-tumor pancreas vs. chronic pancreatitis0.6360.7200.945
Non-tumor pancreas vs. normal control pancreas<0.0010.2660.015
Chronic pancreatitis vs. normal control pancreas<0.0010.0710.025
Upstream pancreas vs. normal control pancreas<0.001<0.0010.003
Downstream pancreas vs. normal control pancreas<0.0010.0030.435

p-Values of one-way analysis of variance (ANOVA) in multiple comparisons of PDAC, non-tumor pancreas, CP, and normal control pancreas.

Differentiation of PDAC Mass From Non-tumor Pancreas

Evaluation by ROC curves and the AUCs for each parameter and combinations of parameters was performed (Figure 8). Using single parameters, T1 values yielded the greatest AUC (95% CI) [0.901 (0.874–0.928)] compared to T2 [0.552 (0.507–0.598)] and ADC [0.712 (0.672–0.752)]. Based on the ROC analysis, the cutoff values of T1, T2, and ADC for the differentiation of PDAC mass undergoing chemotherapy from non-tumor pancreas were 1,494.0, 59.3 ms, and 1.201 × 10–3 mm2/s, respectively. The sensitivity and specificity were for T1 (96.55 and 78.39%), T2 (87.77 and 31.61%), and ADC values (91.60 and 47.74%) were calculated. When combining two parameters, the highest AUC (95% CI) was obtained when combining T1 and ADC [0.934 (0.914–0.953)] compared to a combination of T1 and T2 [0.913 (0.887–0.938)], or a combination of T2 and ADC [0.725 (0.686–0.763)]. When T1, T2, and ADC were combined, the AUC (95% CI) was 0.937 (0.918–0.956) with a sensitivity of 91.54% and a specificity of 85.81%, and this was significantly higher than that using any single parameter (p < 0.001). There was no statistical difference between a combination of T1, T2, and ADC vs. T1 and ADC (p = 0.158). The sensitivity and specificity of each curve is listed in Table 6.

FIGURE 8

TABLE 6

Sensitivity (%)Specificity (%)
T196.5578.39
T287.7731.61
ADC91.6047.74
T1 + T296.5581.61
T1 + ADC93.7382.90
T2 + ADC34.595.8
T1 + T2 + ADC91.5485.81

Sensitivity and specificity of each single parameter and combinations of multi-parameters.

Discussion

In this study, significant differences were found when defining the tissue characteristics of the whole pancreas with PDAC, CP, and normal control pancreas, and comparing PDAC mass, upstream and downstream pancreas using T1, T2, and ADC values. The differences between the measured T1, T2, and ADC values in whole pancreas with PDAC, CP, and normal pancreas may reflect the development of CP to PDAC in some extent and may provide a new follow-up method of CP. The difference between non-tumor pancreas, especially upstream pancreas and downstream pancreas, and in normal control pancreas was pathologically proven to correlate with the inflammatory pathological change, which plays an important role on the tumor response of chemotherapy, recurrence, and 5-year survival rate (Hayano et al., 2016). Furthermore, we found that a combination of T1, T2, and ADC values was superior in characterizing tissue properties when compared to any single parameter.

Normal healthy pancreas has the lowest T1, T2, and ADC values compared to the PDAC group and CP group. No significant differences were found between the pancreatic head, body, and tail (p > 0.05). Normal healthy pancreas has a short T1 value compared with other abdominal organs due to the presence of a high amount of acinar protein and rough endoplasmic reticulum in the pancreatic cells (Tirkes et al., 2017; Noda et al., 2019). In our study, the mean ADC value showed a slight trend of reduction from head to tail (head 1.146 × 10–3 mm2/s, body 1.129 × 10–3 mm2/s, and tail 1.108 × 10–3 mm2/s) due to the heterogeneity of pancreatic tissue composition, and this is consistent with the prior report (Schoennagel et al., 2011).

Chronic pancreatitis results in irreversible pancreatic structural damage, and the pathological changes include pancreatic calcification, chronic inflammation, and fibrosis, which correlate with an increase in T1 signal intensity in pancreatic parenchyma (Yoon et al., 2016). An increased T1 value was demonstrated with progressive disease (normal controls 865 ± 220 ms vs. mild CP 1,075 ± 221 ms vs. severe CP 1,350 ± 139 ms, p < 0.0001) (Wang et al., 2018) in a previous study. In our study, the T1 value in CP was markedly higher than that of the normal control pancreas (p < 0.001) and lower than that of the whole pancreas with PDAC (p < 0.001). CP also has longer T2 and higher ADC than those of the normal control pancreas, which may be related with edema, fibrosis, gland atrophy, and fat infiltration.

The T1, T2, and ADC values of the whole pancreas with PDAC were highest compared with the other two groups for the malignant tumor, fibrosis, and COP in the upstream pancreas and TACP around the tumor. The differences in T1 values were statistically significant within the whole pancreas with PDAC, CP, and normal control pancreas (p < 0.001). However, T2 and ADC values showed no significant differences between the whole pancreas with PDAC and CP for the pathological changes of CP and upstream pancreas.

Pancreatic ductal adenocarcinoma tends to obstruct MPD, which results in COP, fibrosis, and atrophy of the upstream pancreas (Balkwill and Mantovani, 2001; Coussens and Werb, 2002; Mantovani et al., 2008). COP may be caused by the tumor obstruction of MPD, which leads to increased intraductal pressure and damage to the duct membrane or rupture of the secondary ducts resulting in interstitial extravasation of the activated pancreatic enzymes, the recurrence of inflammation, and interstitial damage results in fibrosis hyperplasia and CP. Another hypothesis is that tumor cells may secrete plasminogen-activating enzymes, which may, in turn, activate trypsinogen-inducing auto-digestion. In some patients, tumor can arise from pre-existing CP (Leal and Liby, 2018). On the other hand, the downstream pancreas may be less influenced than the upstream pancreas.

Significant differences (p < 0.05) when comparing non-tumor pancreas with PDAC mass, upstream pancreas with downstream pancreas, and downstream pancreas with normal control pancreas can be interpreted by COP or TACP (Imamura et al., 1995). TACP was found in the adjacent parenchyma of the PDAC mass. In addition, this can be used to interpret why T1 and T2 values in the downstream pancreas showed significant differences when compared with the normal control pancreas (p < 0.05). The results also proved that the inflammatory changes and fibrosis in the upstream pancreas are markedly more severe than the downstream pancreas (Bali et al., 2011).

Receiver operator characteristic analysis showed that for the single parameters, the T1 value is the most accurate in differentiating the tumor and non-tumor area compared with that of the T2 value and ADC value (p < 0.001). The T2 and ADC values have some limitations when a large amount of fibrosis and the change in tumor cellular density are present in the tumor and upstream pancreas. In our study, we found that if three parameters were combined, the AUC was statistically significantly higher compared with any single parameter and two combined parameters (p < 0.05) with the exception of the combined T1 and ADC values (p = 0.158). Therefore, it appears that the T1 value is critical in tissue characterization, while the ADC value can provide additional helpful information but has lower sensitivity and specificity than the T1 value. Based on our data, the T2 value only provides minimal information in tissue characterization.

Our study has some limitations. First, the relatively small numbers of patients may affect the statistical outcomes. These findings will need to be validated in larger studies. Second, patients were not sub-grouped according to the grade of cancer; therefore, the T1, T2, and ADC values may show wide variability. Subgrouping of patients was not possible due to the small pilot nature of this study. Third, most patients with PDAC underwent chemotherapy at the time of study. This may have altered the tumor characteristics. We found from our research group experience that patients who were treatment naïve are difficult to recruit. Fourth, the ROIs were obtained according to the previous imaging data and other non-contrast MR images, and there was an interval between the MR scan and surgery. The tumor margin could not be validated immediately, so we selected a relatively smaller ROI within the tumor, avoiding the necrosis area. Fifth, T1 and T2 mapping only included three slices because we used the T1 MOLLI sequence. This would not reflect the global features of the tumor. Future studies should recruit patients prior to initiation of treatment and over the course of the treatment so that there are independent predictors of prognosis in the responder and non-responder. Finally, the T1, T2, and ADC maps had different resolutions at the time of image acquisition, which required retrospective calibration to keep the pixel numbers consistent. In future studies, we will keep the same T1, T2, and DWI image resolutions and continuous or three-dimensional T1 and T2 mapping sequences in imaging acquisition to decrease the bias. For the multiparametric mapping, which can demonstrate the differences of PDAC mass undergoing chemotherapy, upstream and downstream pancreas, we believe that it can be used as a biomarker to predict tumor or longitudinal follow-up after chemotherapy.

Conclusion

In conclusion, multiparametric mapping is feasible for the evaluation of the differences between PDAC, CP, and normal pancreas tissue. The combination of multiple parameters of T1, T2, and ADC provides a higher accuracy compared to the result with any single parameter in tissue characterization of the pancreas (the main points of this article are listed in Figure 9).

FIGURE 9

Statements

Data availability statement

All datasets generated for this study are included in the article/supplementary material.

Ethics statement

The studies involving human participants were reviewed and approved by the Cedars-Sinai Institutional Review Board. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

LW, DL, and SP conceived and designed the study. LW and NW analyzed and interpreted the data. NW, ZD, ZZ, and AH collected the data. LW and SG wrote the manuscript. SG, ZF, TJ, AC, FH, SL, and AW critically revised the manuscript. DL and SP approved the final version of the manuscript. LW and YX contributed to statistical analysis.

Funding

This research was supported by the Biomedical Imaging Research Institute and Division of Digestive and Liver Diseases of Cedars-Sinai Medical Center, and National Institutes of Health US U01 DK108314-05.

Acknowledgments

The research group thanks all the patients and normal volunteers who participated in this study. In addition, we thank Xiaoming Bi, Edward Gill, Joselyn Ricafrente, and Laura Smith who assisted with this study and the MR scan.

Conflict of interest

FH was employed by company Siemens Healthineers. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

  • ADC

    apparent diffusion coefficient

  • ANOVA

    one-way analysis of variance

  • AUC

    area under curve

  • COP

    chronic obstructed pancreatitis

  • CP

    chronic pancreatitis

  • DWI

    diffusion weighted imaging

  • MPD

    main pancreatic duct

  • MRI

    magnetic resonance imaging

  • PDAC

    pancreatic ductal adenocarcinoma

  • ROC

    receiver operator characteristic

  • T1 MOLLI

    T1 modified Look-Locker Inversion Recovery

  • T1-VIBE-DIXON

    T1-weighted volumetric interpolated breath-hold examination

  • T2 HASTE

    T2 half-Fourier acquisition single-shot turbo spin–echo sequence

  • TACP

    tumor-associated chronic pancreatitis.

References

  • 1

    American Cancer Society (2018). Cancer Facts & Figures 2018 American Cancer Society.New York, NY: American Cancer Society.

  • 2

    ApprichS.MamischT. C.WelschG. H.StelzenederD.AlbersC.TotzkeU.et al (2012). Quantitative T2 mapping of the patella at 3.0T is sensitive to early cartilage degeneration, but also to loading of the knee.Eur. J. Radiol.81:e00438-43. 10.1016/j.ejrad.2011.03.069

  • 3

    AttiyehM. A.ChakrabortyJ.McintyreC. A.KappagantulaR.ChouY.AskanG.et al (2019). CT radiomics associations with genotype and stromal content in pancreatic ductal adenocarcinoma.Abdom. Radiol.443148–3157. 10.1007/s00261-019-02112-1

  • 4

    BaliM. A.MetensT.DenolinV.DelhayeM.DemetterP.ClossetJ.et al (2011). Tumoral and nontumoral pancreas: correlation between quantitative dynamic contrast-enhanced MR imaging and histopathologic parameters.Radiology261456–466. 10.1148/radiol.11103515

  • 5

    BalkwillF.MantovaniA. (2001). Inflammation and cancer: back to Virchow?Lancet357539–545. 10.1016/s0140-6736(00)04046-0

  • 6

    BanerjeeR.PavlidesM.TunnicliffeE. M.PiechnikS. K.SaraniaN.PhilipsR.et al (2014). Multiparametric magnetic resonance for the non-invasive diagnosis of liver disease.J. Hepatol.6069–77. 10.1016/j.jhep.2013.09.002

  • 7

    BlystadI.WarntjesJ. B. M.SmedbyO.LundbergP.LarssonE. M.TisellA. (2017). Quantitative MRI for analysis of peritumoral edema in malignant gliomas.PLoS One12:e0177135. 10.1371/journal.pone.0177135

  • 8

    CassinottoC.FeldisM.VergniolJ.MouriesA.CochetH.LapuyadeB.et al (2015). MR relaxometry in chronic liver diseases: comparison of T1 mapping, T2 mapping, and diffusion-weighted imaging for assessing cirrhosis diagnosis and severity.Eur. J. Radiol.841459–1465. 10.1016/j.ejrad.2015.05.019

  • 9

    ChenC.ChenJ.XiaC.HuangZ.SongB. (2017). T1 mapping combined with Gd-EOB-DTPA-enhanced magnetic resonance imaging in predicting the pathologic grading of hepatocellular carcinoma.J. Biol. Regul. Homeostat. Agents311029–1036.

  • 10

    ChoiS.-Y.KimS. H.KangT. W.SongK. D.ParkH. J.ChoiY.-H. (2016). Differentiating mass-forming autoimmune pancreatitis from pancreatic ductal adenocarcinoma on the basis of contrast-enhanced MRI and DWI findings.Am. J. Roentgenol.206291–300. 10.2214/AJR.15.14974

  • 11

    CieszanowskiA.Anysz-GrodzickaA.SzeszkowskiW.KaczynskiB.MajE.GornickaB.et al (2012). Characterization of focal liver lesions using quantitative techniques: comparison of apparent diffusion coefficient values and T2 relaxation times.Eur. Radiol.222514–2524. 10.1007/s00330-012-2519-x

  • 12

    CoussensL. M.WerbZ. (2002). Inflammation and cancer.Nature420860–867.

  • 13

    FreenyP. C. (1989). Classification of pancreatitis.Radiol. Clin. North Am.271–3.

  • 14

    FukukuraY.TakumiK.KamimuraK.ShindoT.KumagaeY.TateyamaA.et al (2012). Pancreatic adenocarcinoma: variability of diffusion-weighted MR imaging findings.Radiology263732–740. 10.1148/radiol.12111222

  • 15

    GuoH.AuW. Y.CheungJ. S.KimD.JensenJ. H.KhongP. L.et al (2009). Myocardial T2 quantitation in patients with iron overload at 3 Tesla.J. Magn. Reson. Imaging30394–400. 10.1002/jmri.21851

  • 16

    HaoL.ZengX. P.XinL.WangD.PanJ.BiY. W.et al (2017). Incidence of and risk factors for pancreatic cancer in chronic pancreatitis: a cohort of 1656 patients.Dig. Liver Dis.491249–1256. 10.1016/j.dld.2017.07.001

  • 17

    HayanoK.MiuraF.WadaK.SuzukiK.TakeshitaK.AmanoH.et al (2016). Diffusion-weighted MR imaging of pancreatic cancer and inflammation: prognostic significance of pancreatic inflammation in pancreatic cancer patients.Pancreatology16121–126. 10.1016/j.pan.2015.10.004

  • 18

    HechtE. M.LiuM. Z.PrinceM. R.JambawalikarS.RemottiH. E.WeisbergS. W.et al (2017). Can diffusion-weighted imaging serve as a biomarker of fibrosis in pancreatic adenocarcinoma?J. Magn. Reson. Imaging46393–402. 10.1002/jmri.25581

  • 19

    IchikawaT.ErturkS. M.MotosugiU.SouH.IinoH.ArakiT.et al (2007). High-b value diffusion-weighted MRI for detecting pancreatic adenocarcinoma: preliminary results.Am. J. Roentgenol.188409–414. 10.2214/ajr.05.1918

  • 20

    ImamuraT.IguchiH.ManabeT.OhshioG.YoshimuraT.WangZ.et al (1995). Quantitative analysis of collagen and collagen subtypes I, III, and V in human pancreatic cancer, tumor-associated chronic pancreatitis, and alcoholic chronic pancreatitis.Pancreas11357–364. 10.1097/00006676-199511000-00007

  • 21

    KaliA.ChoiE.-Y.SharifB.KimY. J.BiX.SpottiswoodeB.et al (2015). Native T1 mapping by 3-T CMR imaging for characterization of chronic myocardial infarctions.JACC81019–1030. 10.1016/j.jcmg.2015.04.018

  • 22

    KaramitsosT. D.PiechnikS. K.BanypersadS. M.FontanaM.NtusiN. B.FerreiraV. M.et al (2013). Noncontrast T1 mapping for the diagnosis of cardiac amyloidosis.JACC6488–497. 10.1016/j.jcmg.2012.11.013

  • 23

    KirkegardJ.MortensenF. V.Cronin-FentonD. (2017). Chronic pancreatitis and pancreatic cancer risk: a systematic review and meta-analysis.Am. J. Gastroenterol.1121366–1372. 10.1038/ajg.2017.218

  • 24

    LealA. S.LibyK. T. (2018). Murine models of pancreatitis leading to the development of pancreatic cancer.Curr. Protoc. Pharmacol.83:e48. 10.1002/cpph.48

  • 25

    LuzL. P.Al-HaddadM. A.SeyM. S.DewittJ. M. (2014). Applications of endoscopic ultrasound in pancreatic cancer.World J. Gastroenterol.207808–7818. 10.3748/wjg.v20.i24.7808

  • 26

    MantovaniA.AllavenaP.SicaA.BalkwillF. (2008). Cancer-related inflammation.Nature454:436. 10.1038/nature07205

  • 27

    NattermannC.GoldschmidtA. J.DancygierH. (1993). Endosonography in chronic pancreatitis–a comparison between endoscopic retrograde pancreatography and endoscopic ultrasonography.Endoscopy25565–570. 10.1055/s-2007-1010406

  • 28

    NodaY.GoshimaS.TsujiY.KajitaK.KawadaH.KawaiN.et al (2019). Correlation of quantitative pancreatic T1 value and HbA1c value in subjects with normal and impaired glucose tolerance.J. Magn. Reson. Imaging49711–718. 10.1002/jmri.26242

  • 29

    ParkH. S.LeeJ. M.ChoiH. K.HongS. H.HanJ. K.ChoiB. I. (2009). Preoperative evaluation of pancreatic cancer: comparison of gadolinium-enhanced dynamic MRI with MR cholangiopancreatography versus MDCT.J. Magn. Reson. Imaging30586–595. 10.1002/jmri.21889

  • 30

    SahaiA. V.ZimmermanM.AabakkenL.TarnaskyP. R.CunninghamJ. T.Van VelseA.et al (1998). Prospective assessment of the ability of endoscopic ultrasound to diagnose, exclude, or establish the severity of chronic pancreatitis found by endoscopic retrograde cholangiopancreatography.Gastrointest. Endosc.4818–25. 10.1016/s0016-5107(98)70123-3

  • 31

    SchoennagelB. P.HabermannC. R.RoeschM.HahneJ. D.ArndtC.KleibelerL.et al (2011). Diffusion-weighted imaging of the healthy pancreas: apparent diffusion coefficient values of the normal head, body, and tail calculated from different sets of b-values.J. Magn. Reson. Imaging34861–865. 10.1002/jmri.22743

  • 32

    ShiJ.XueJ. (2019). Inflammation and development of pancreatic ductal adenocarcinoma.Chin. Clin. Oncol.8:19. 10.21037/cco.2019.04.02

  • 33

    SiegelR. L.MillerK. D.JemalA. (2018). Cancer statistics, 2018.CA Cancer J. Clin.687–30. 10.3322/caac.21442

  • 34

    SinghA.FaulxA. L. (2016). Endoscopic evaluation in the workup of pancreatic cancer.Surg. Clin. North Am.961257–1270. 10.1016/j.suc.2016.07.006

  • 35

    TirkesT.LinC.CuiE.DengY.TerritoP. R.SandrasegaranK.et al (2018). Quantitative MR evaluation of chronic pancreatitis: extracellular volume fraction and MR relaxometry.AJR Am. J. Roentgenol.210533–542. 10.2214/AJR.17.18606

  • 36

    TirkesT.LinC.FogelE. L.ShermanS. S.WangQ.SandrasegaranK. (2017). T1 mapping for diagnosis of mild chronic pancreatitis.J. Magn. Reson. Imaging451171–1176. 10.1002/jmri.25428

  • 37

    Vietti VioliN.HilbertT.BastiaansenJ. A. M.KnebelJ. F.LedouxJ. B.StemmerA.et al (2019). Patient respiratory-triggered quantitative T2 mapping in the pancreas.J. Magn. Reson. Imaging50410–416. 10.1002/jmri.26612

  • 38

    WangM.GaoF.WangX.LiuY.JiR.CangL.et al (2018). Magnetic resonance elastography and T1 mapping for early diagnosis and classification of chronic pancreatitis.J. Magn. Reson. Imaging10.1002/jmri.26008[Epub ahead of print],

  • 39

    WangY.ChenZ. E.NikolaidisP.MccarthyR. J.MerrickL.SternickL. A.et al (2011). Diffusion-weighted magnetic resonance imaging of pancreatic adenocarcinomas: association with histopathology and tumor grade.J. Magn. Reson. Imaging33136–142. 10.1002/jmri.22414

  • 40

    YamamotoT.YagiS.KinoshitaH.SakamotoY.OkadaK.UryuharaK.et al (2015). Long-term survival after resection of pancreatic cancer: a single-center retrospective analysis.World J. Gastroenterol.21262–268. 10.3748/wjg.v21.i1.262

  • 41

    YinT.PeetersR.FengY.LiuY.YuJ.DymarkowskiS.et al (2017). Characterization of a rat orthotopic pancreatic head tumor model using three-dimensional and quantitative multi-parametric MRI.NMR Biomed.30:e3676. 10.1002/nbm.3676

  • 42

    YoonJ. H.LeeJ. M.LeeK. B.KimS. W.KangM. J.JangJ. Y.et al (2016). Pancreatic steatosis and fibrosis: quantitative assessment with preoperative multiparametric MR imaging.Radiology279140–150. 10.1148/radiol.2015142254

Summary

Keywords

pancreatic ductal adenocarcinoma, chronic pancreatitis, magnetic resonance imaging, parametric mapping, T1, T2, ADC

Citation

Wang L, Gaddam S, Wang N, Xie Y, Deng Z, Zhou Z, Fan Z, Jiang T, Christodoulou AG, Han F, Lo SK, Wachsman AM, Hendifar AE, Pandol SJ and Li D (2020) Multiparametric Mapping Magnetic Resonance Imaging of Pancreatic Disease. Front. Physiol. 11:8. doi: 10.3389/fphys.2020.00008

Received

12 July 2019

Accepted

09 January 2020

Published

21 February 2020

Volume

11 - 2020

Edited by

Peter Hegyi, University of Szeged, Hungary

Reviewed by

Nandor Faluhelyi, University of Pécs Medical School, Hungary; Shin Hamada, Tohoku University, Japan

Updates

Copyright

*Correspondence: Stephen J. Pandol, Debiao Li,

†These authors have contributed equally to this work

This article was submitted to Gastrointestinal Sciences, a section of the journal Frontiers in Physiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics