Skip to main content

Evaluation of pancreatic tumor development in KPC mice using multi-parametric MRI

Abstract

Background

Pancreatic ductal adenocarcinoma (PDA) is a fatal disease with very poor prognosis. Development of sensitive and noninvasive methods to monitor tumor progression in PDA is a critical and unmet need. Magnetic resonance imaging (MRI) can noninvasively provide information regarding underlying pathophysiological processes such as necrosis, inflammatory changes and fibrotic tissue deposition.

Methods

A genetically engineered KPC mouse model that recapitulates human PDA was used to characterize disease progression. MR measures of T1 and T2 relaxation times, magnetization transfer ratio (MTR), diffusion and chemical exchange saturation transfer were compared in two separate phases i.e. slow and rapid growth phase of tumor. Fibrotic tissue accumulation was assessed histologically using Masson’s trichrome staining. Pearson correlation coefficient (r) was computed to assess the relationship between the fibrotic tissue accumulation and different MR parameters.

Results

There was a negative correlation between amide proton transfer signal intensity and tumor volume (r = − 0.63, p = 0.003) in the slow growth phase of the tumor development. In the terminal stage of rapid growth phase of the tumor development MTR was strongly correlated with tumor volume (r = 0.62, p = 0.008). Finally, MTR was significantly correlated with % fibrosis (r = 0.87; p < 0.01), followed by moderate correlation between tumor volume (r = 0.42); T1 (r = − 0.61), T2 (r = − 0.61) and accumulation of fibrotic tissue.

Conclusions

Here we demonstrated, using multi-parametric MRI (mp-MRI), that MRI parameters changed with tumor progression in a mouse model of PDA. Use of mp-MRI may have the potential to monitor the dynamic changes of tumor microenvironment with increase in tumor size in the transgenic KPC mouse model of pancreatic tumor.

Background

Pancreatic ductal adenocarcinoma (PDA) is the most lethal form of human cancer [1]. Desmoplasia, a hallmark pathologic feature of PDA, is characterized by the presence of a robust stroma containing fibroblasts and inflammatory cells [2]. Diagnosis of pancreatic cancer is usually made at later stages of the disease making it even more difficult to treat. Clear understanding of tumor progression may help in identifying PDA at an early stage. Therefore, in order to characterize the tumor progression, we need 1) to investigate tumor development in an effective preclinical model that closely mimics human disease progression and 2) a sensitive, non-invasive monitoring tool that would provide detailed information regarding the disease progression and can be used in future clinical studies.

As an experimental model for pancreatic cancer, tumors have been implanted subcutaneously and within the pancreas [3,4,5]. However, both of these models do not parallel the human disease progression. A genetically engineered mouse model, having genetic alterations in genes K-rasLSL-G12D/+; Trp53LSL-R172H/+; Cre (KPC) offers an alternative to transplantation models for preclinical therapeutic evaluation as it expresses mutations similar to human pancreatic cells [6] and develops pancreatic tumors in which pathophysiology and molecular features resemble those of human PDA [7]. There is a dire need to develop non-invasive techniques to monitor disease progression as well as study the therapeutic effects of existing and novel chemotherapeutic agents.

Magnetic resonance imaging (MRI) has been proven to be extremely useful in clinical trials for monitoring tumor development [8] and assessing therapeutic effects of novel therapeutic strategies [9,10,11]. Currently, management decisions in almost all phases of diagnosis, treatment, and follow-up rely on gold standard magnetic resonance (MR) measures i.e. spin-lattice relaxation time (T1) and spin-spin relaxation time (T2) [12, 13]. However, T1 and T2 measures alone may not be sensitive enough to investigate the entire spectrum of tumor properties. Recently, diffusion weighted imaging (DWI) [14, 15] and magnetization transfer imaging (MTI) [16,17,18] have been used as complementary imaging techniques to conventional MR measures to characterize adenocarcinoma. Finally, amide proton transfer (APT) imaging based on the chemical exchange saturation transfer (CEST) approach, which provides valuable information regarding the tumor microenvironment, has drawn considerable attention as a novel MRI contrast agent in the field of molecular imaging [19, 20]. To monitor the dynamic microenvironment of the tumor as in PDA, we need to employ mp-MRI to understand the disease progression, and the tumor environment thereby helping to plan the appropriate therapeutic regimen. In fact, several studies have demonstrated that multiparameteric MR based approach enables better quantitation of pathological processes in abdominal solid organs [21,22,23].

Based on the promising results of these studies, we hypothesized the mp-MRI may enable to accurately quantify the disease progression in KPC mice. In this study, we conducted MRI at a high field strength of 14 Tesla (T) utilizing the T1, T2, DWI, MTI and CEST imaging sequences to monitor the progression of pancreatic ductal adenocarcinoma in KPC mice. Additionally, we compared MR parameters with the histologic markers of fibrotic tissue accumulation in pancreatic tumors in the KPC mouse model.

Methods

Animal handling and care

The study was conducted with the approval from our institutional animal care and use committee (IACUC). KPC mice (n = 16) were enrolled in the study when they had a small palpable mass, which was confirmed by ultrasound imaging and MR imaging. All mice were imaged upon enrollment and at the final time point when the tumor reached terminal size (10 mm in any one direction). A subgroup of KPC mice (n = 9) were imaged weekly between the baseline and final time points. All mice were euthanized at the terminal time point and the tumor and associated pancreatic tissue was excised and prepared for histological evaluation.

MR data acquisition

MRI experiments were performed on a 14 T Bruker Avance 600 MHz/89 mm wide-bore vertical MR spectrometer (Bruker Corp., Billerica, MA). A birdcage coil (inner diameter 25 mm) was used to image the animal mounted on a cradle with a respiratory monitoring probe. Animals were anesthetized before being secured to the custom-built cradle. The coil was then inserted vertically into the scanner. For the entire duration of the experiment core body temperature of mice was maintained at 37 °C. The entire MR data acquisition took 50–60 min during which the animals were continuously monitored for respiratory rate. All the images were gated to respiration of the animal. Following the MR acquisition, mice were removed from the coil and allowed to recover. Longitudinal MRI was implemented to measure pancreatic tumor growth in KPC mice (n = 9) over a period of 6 weeks. Mice included in longitudinal MRI had a palpable mass at the baseline measurement, which was measured by ultrasound (US) and later corroborated using MRI.

Anatomical images

The multi-slice MRI protocol, covering the whole tumor, started with fat-suppressed T1 weighted coronal images (repetition time (TR) = 2000 ms, echo time (TE) = 5.49 ms, number of averages (NA) = 1, field of view (FOV) = 30 × 30 mm2, rare factor = 8, matrix size = 128 × 256; yielding spatial resolution of 0.234 × 0.117 mm/pixel) for anatomical reference. Subsequently, 20 axial images were acquired to calculate tumor volume.

T1 mapping

Multiple images using rapid acquisition with refocused echoes (RARE) were acquired using following parameters: TE = 9.66 ms, TR = 5500, 3000, 1500, 1000, 385.8 ms, NA = 1, FOV = 30 × 30 mm2, rare factor = 2, matrix size = 256 × 128 (reconstructed phase encoding steps = 128; acquisition phase encoding steps = 96) yielding spatial resolution of 0.117 × 0.234 mm/pixel. The data acquisition time was approximately 9 min.

T2 mapping

Multiple spin-echo data were acquired in coronal orientation covering the area from liver to kidneys. The quantitative T2 maps were generated using a multi-slice multi echo sequence, with fat signal suppressed, utilizing following parameters: TR = 4000 ms; TE = 12 echoes equally spaced from 6.28 ms to 75.4 ms; NA = 1; FOV = 30 × 30 mm2; matrix size = 256 × 128 (reconstructed phase encoding steps = 128; acquisition phase encoding steps = 91) yielding spatial resolution of 0.117 × 0.234 mm/pixel. To cover the entire abdominal region, 10 contiguous slices were acquired without any inter slice gap. The data acquisition time was approximately 6 min.

Magnetization transfer (MT)

MT ratios (MTR) were acquired using a gradient echo sequence (TR/TE = 625/2 ms, flip angle = 30°) with an off-resonance frequency of 7000 Hz and a saturation pulse block pulse shape, 50 ms width, and 10 μT amplitude. A series of 10 images were acquired with FOV = 30 × 30 mm2, matrix size = 256 × 256 yielding spatial resolution of 0.117 × 0.117 mm/pixel. The acquisition time for data acquisition was approximately 3 min.

Diffusion weighted imaging (DWI)

An echo planar imaging (EPI) diffusion measurement (echo train length = 16, pulse duration = 3.0 ms and diffusion time = 7.46 ms) was performed to acquire series of 10 slices using following parameters: TR = 2500 ms; TE = 17.7 ms; NA = 1; FOV = 30 × 30 mm2; matrix size = 128 × 128 yielding spatial resolution of 0.234 × 0.234 mm/pixel. Diffusion weighted measurements were acquired with 8 different b values (0, 30, 60, 100, 150, 200, 300, 500 s/mm2). The data acquisition time was 2 m 40s.

Chemical exchange transfer saturation (CEST) imaging

On a single 1 mm slice, delineating tumor, amide proton transfer (APT) imaging was performed with respiratory gating using small animal monitoring device (SA instruments, Inc., Stony Brook, NY). CEST imaging data were acquired using RARE sequence (continuous-wave block pulse, B1 = 0.5 μT, duration = 2 s), which was applied at 25 frequency offsets from − 360 Hz to 360 Hz with an interval of 0.5 ppm to estimate a center frequency shift (water saturation shift referencing (WASSR) approach) [24,25,26]. Other imaging parameters were: TR/TE = 2200/7 ms, FOV = 30 × 30 mm2, matrix size = 128 × 128, flip angle = 180o, and number of excitations = 1. For saturation a single slice was acquired with 6 frequency offsets at ±3.0, ±3.5, ±4.0 ppm, with an off-resonance RF pulse applied for 3 s at a power level of 2 μT. Other parameters were: TR/TE = 5000/7 ms, matrix = 128 × 128 (reconstructed phase encoding steps = 128; acquisition phase encoding steps = 96), FOV = 30 × 30 mm2, rare factor = 8. Finally, a control image with the saturation offset at 300 ppm was also acquired. Total acquisition time for each animal was approximately 19 min.

Ultrasound (US) data acquisition

The animal was anaesthetized, placed in the supine position and images of the whole tumor were acquired (550 s probe, Vevo 2100, Fujifilm Visualsonics, Toronto, Canada). Ultrasound images were taken every 0.5 mm, in the transverse plane, throughout the whole tumor. The area of the tumor in each image was determined using Vevo LAB v2.1.0. The tumor volume was calculated by multiplying the area with the inter-slice distance.

Image analyses

All raw MR images were processed using Image-J software (http://imagej.nih.gov/ij/), to measure mean values of the different tumors. Regions of interest (ROI) were drawn to circumscribe the entire tumor (Additional file 1: Figure S1). Anatomical Images were used to measure tumor volume. Tumor volume was measured and reconstructed using Amira (Visualization Sciences Group, Burlington, MA), a 3-D software platform [18]. T1 and T2 maps: Maps were generated using T1 and T2 weighted images. MTR maps: The MTR was measured using the following ratio: (SI0 - SIs/SI0), where SI0 represents the tissue signal intensity without saturation pulse applied while SIs represents the tissue signal intensity with saturation pulse. Diffusion maps: Diffusion weighted MR signal decay was analyzed using mono-exponential model: Sb/S0 = exp.(−b∙ADC). Where Sb is the MRI signal intensity with diffusion weighting b, S0 is the non-diffusion-weighted signal intensity and ADC is the apparent diffusion coefficient. In addition to mono-exponential model, a bi-exponential model was used to estimate intra voxel incoherent motion (IVIM) related parameters of perfusion fraction (or pseudo-diffusion) and diffusion [27]. Three lowest b values of 0, 30 and 60 s/mm2 were used to calculate perfusion component (or pseudo-diffusion) whereas rest of the 5 b values of 100, 150, 200, 300 and 500 s/mm2 were used to calculate tissue diffusivity component. Finally, APT-MR images were quantified using the following equation: [Ssat (− 3.5 ppm) – Ssat (3.5 ppm)]/S0 where Ssat and S0 are the water signal intensities measured with and without saturation pulse.

Histogram analyses

Regions of interest were selected across consecutive three slices having greatest tumor area to increase the coverage and improve the reliability. For each tumor, values from each slice were combined to generate a single histogram. Histograms were then normalized by plotting the percentage of pixels remaining above the specific measurement in the x-axis, generating a cumulative histogram. Finally, the area under curve (AUC) with units of % pixels x parameter on x-axis was used to compare values in two different groups.

Histological analysis

After extracting the tumor, it was immediately embedded in optimum cutting temperature (OCT) medium. Three serial, 8-μm thick sections were cut every 1 mm through the entire tumor (CM1950, Leica Biosystems Inc., Buffalo Grove, Illinois) and the sequential sections at each level were stained with Masson’s trichrome to visualize fibrotic tissue deposition. The three sections, matching with MR images, were used for fibrotic tissue quantification. All sections were examined using a Nikon 80i upright microscope (Nikon, Melville, New York). Images of whole sections stained with Masson’s trichrome were acquired with a 10 x objective lens. Fibrotic tissue, identified by blue staining, was separated by thresholding of the hue (130–190), saturation (20–255), and brightness (10–240) values using ImageJ (ImageJ 1.42 National Institutes of Health, Bethesda, MD). The fibrotic tissue content is presented as the percent area of fibrosis over the whole tumor area. High resolution, 40 μm thick, cryo-images of one of the tumor bearing mice were acquired at Bioinvision® (www.bioinvision.com).

Statistical analysis

Statistical analysis was performed using Graph Pad prism version 6.0 (GraphPad, La Jolla, CA, USA). Tumor volume was compared between baseline and final time points using paired t test. Pearson correlation coefficient (r) with Bonferroni correction was computed to assess the relationship between the tumor volume and different MR parameters. All data were presented as means and standard deviations (SD). Statistical significance was accepted for p < 0.05.

Results

Increase in pancreatic tumor volume with age in KPC mice

Longitudinal MRI demonstrated two distinct phases existed in tumor progression in KPC mice i.e. slow phase and rapid phase. Tumor growth, measured by MRI, in 6 out of 9 KPC mice increased exponentially once it surpassed a threshold value of 250 mm3 (Fig. 1a). Overall, there was a significant increase in tumor mass from baseline measurements (205.5 ± 154.4 mm3, mean ± SD) to the 6-week time point (455.9 ± 137.3 mm3). A moderate correlation was found between volume measurement by US and MRI (r = 0.59) (Fig. 1b). Finally, a high-resolution MR coronal section was compared with a whole body cryo-image (Fig. 1c) and MR axial section was compared with axial US image (Fig. 1d) to compare spatial resolution of different modalities.

Fig. 1
figure 1

a Progressive tumor volumes for individual KPC mice (n = 9) that were imaged longitudinally. There was an exponential increase in tumor volume, in 6 out of 9 animals, once it reached a certain threshold value i.e. 250 mm3 (dotted line) in this study. b Correlation between US and MRI measurements for tumor volume. c Representative whole body image of tumor in KPC mice (Bioinvision, Cleveland, OH) and coronal MR image of the same mouse. d Representative axial MR and US image of tumor in the same KPC mouse. Red arrows point towards pancreatic tumor (c and d)

MR parameters correlate with tumor volume

Multi-parametric data was analyzed based on the size of the tumor i.e. < 250 mm3 and > 250 mm3 as shown in Figs. 2 and 3. We chose a threshold value of 250 mm3 to differentiate between a slow and rapid growth phase in tumor development. Tumor size is an important prognostic factor in various cancers. Indeed, studies have documented a certain cut off size in various tumors enhance the prognosis [28,29,30]. Furthermore, there was a strong correlation between tumor volume and amide proton transfer (APT) signal intensity (r = − 0.63, p = 0.003) in < 250 mm3 whereas no correlation was present between APT signal intensity and tumor volume (r = − 0.03) when the tumor size was > 250 mm3. When the tumor size surpassed the threshold value i.e. 250 mm3, magnetization transfer ratio (MTR) displayed strong positive correlation with the tumor volume (r = 0.62, p = 0.008). Based on intravoxel incoherent motion (IVIM) model, [27, 31, 32] pseudo-diffusion or perfusion component (using low-b value) showed a strong negative correlation with the tumor volume (r = − 0.70, p = 0.003) whereas we did not find any correlation between diffusion (high-b values) and tumor volume (r = 0.12; < 250mm3, r = 0.11; > 250 mm3).

Fig. 2
figure 2

Relationship between tumor volume and different MR parameters. a-d show the relationship between tumor volume and T1, T2, MTR and APT signal intensity respectively, when the tumor volume is less than 250 mm3. e-h show the relationship between tumor volume and T1, T2, MTR and APT signal intensity respectively, when the tumor volume exceeds 250 mm3

Fig. 3
figure 3

Relationship between tumor volume and diffusion measurements. a and b show the relationships between tumor volume and ADC (pseudo-diffusion), and ADC (high-b values), respectively, when the tumor volume is smaller than 250 mm3. c and d show the relationships between tumor volume and ADC (pseudo-diffusion), and ADC (high-b values), respectively, when the tumor volume is larger than 250 mm3

Histogram analyses of tumor progression

To assess pixel-pixel changes in MR parameters regionally, histogram analysis was implemented in the tumor area of KPC mice. Histograms were generated for T1, T2, MTR, and apparent diffusion coefficient (ADC) by combining values from 3 different slices for the pancreatic tumor. There was a general shift toward higher values for MTR and lower values for ADC in a larger tumor group i.e. > 250 mm3 (Fig. 4a). Additionally, the difference between smaller and larger tumor groups was better visualized in a cumulative histogram, which was produced by plotting the percentage of pixels remaining above the x-axis values (Fig. 4b). Finally, area under curve (AUC) was calculated from cumulative histograms for mp-MRI did not show any significant difference between < 250 mm3 and > 250 mm3 (Fig. 5).

Fig. 4
figure 4

a Representative figure of frequency distribution of the pixels of different MR parameters in a KPC mouse at baseline and final time points. b Representative figure of percentage pixels above T1, T2, MTR, ADC and APT value in a KPC mouse at baseline and final time points

Fig. 5
figure 5

Area under curve (AUC) demonstrating differences in multi-parametric MR measures in PDA of < 250 mm3 and > 250 mm3 in 5 KPC mice

Histological analysis

Fibrotic tissue deposition was identified by Masson’s trichrome stained tumors and compared to mp-MR maps (Fig. 6). There was a significant correlation between MTR and % fibrosis (r = 0.87; p < 0.01), followed by moderate correlation between tumor volume (r = 0.42; p = 0.30), T1 (r = 0.59, p = 0.13), T2 (r = − 0.61, p = 0.10) and fibrotic tissue deposition (Fig. 7).

Fig. 6
figure 6

Representative Masson’s Trichrome and H&E stains for a smaller tumor (a) and larger tumor (b) from KPC mice and corresponding anatomic images and colored maps with T1, T2, MTR and ADC measures for pancreatic tumor

Fig. 7
figure 7

Correlation between % fibrosis and tumor volume and other MR parameters. There was moderate correlation between increase in tumor volume, T1, T2 and increase in fibrotic tissue accumulation (a, b, c). MTR % was significantly correlated with increase in fibrotic tissue accumulation (d)

Discussion

The present study evaluated multi-parametric MRI (mp-MRI) to monitor pancreatic tumor progression in the KPC mouse model. MR measurements consisted of T1, T2, ADC, MTR and APT imaging. The results from the present study revealed that 1) the tumor size increases exponentially once it crosses a threshold value i.e. 250 mm3; 2) MTR (%) and ADC (pseudo-diffusion) demonstrated higher correlation with tumor volume compared to other measures especially during the later stages of tumor development; and 3) MTR was significantly correlated with increase in fibrotic tissue accumulation.

We have demonstrated  a moderate correlation between US and MRI measures of tumor volume. Although US has been used extensively as an imaging tool in small animals, poor signal-to-noise ratio (SNR) images is one of the major limitations to ultrasound use. On the other hand, MRI has the ability to provide images in great anatomic detail. Decreased SNR could be one of the reasons that explains a moderate correlation between tumor volume measurements using US and MRI. Additionally, tumor burden or size has been demonstrated to be an important prognostic factor among patients with different types of carcinomas [28,29,30]. Furthermore, it has been suggested that tumors grow in a non-linear fashion [33]. Similarly, in the present study we have demonstrated that the growth of tumors follows a non-linear trend i.e. in two separate phases i.e. slow and rapid. Finally, we found that a cut off size of 250 mm3 in the KPC mouse model may provide enhanced prognosis.

The measurement of any MR parameter in isolation may undermine the dynamic nature of tumor microenvironment. However, the use of mp-MRI enables the evaluation of multiple parameters to give a more representative picture of the tumor microenvironment. Methods such as quantitative T1, T2, ADC and MTR are well established methods used to characterize tumor progression [8, 34,35,36]. T1 and T2 relaxation times of water molecules in tissue have been demonstrated as sensitive indicators of tumor progression as well as responses to different therapeutic agents [35, 37]. T1 is sensitive to factors such as the amount of 1) water in the extracellular space and 2) protein in the water [38]. Similarly, the results of our study demonstrate a moderate positive correlation between T1 and tumor volume in the smaller tumor group (< 250 mm3) suggesting increased water and protein content especially during the earlier stages of tumor development. Conversely, when the tumor size increased beyond the threshold value (250 mm3), T1 failed to provide detailed information regarding the tumor development. In advanced tumor stages there is more fibrotic tissue deposition when compared to earlier stages [39], which can be demonstrated by a decrease in T2 relaxation time. Similarly, in the present study, we revealed a moderate negative correlation between T2 and tumor volume in the larger tumor group and a decrease in T2 relaxation time was also correlated to an increase in fibrotic tissue content. These findings suggested that fibrotic tissue accumulated during later stages of disease progression and that T2 is more sensitive than T1 to tumor progression, especially in terms of fibrotic tissue accumulation.

DWI is sensitive to the thermally driven motion of water molecules along the orientation of additional diffusion sensitizing gradients applied during MR sequence. The signal attenuation coefficient, known as apparent diffusion coefficient (ADC) can be derived from diffusion weighted images. ADC has been shown to be sensitive to tissue structure at the cellular level [40, 41]. Results from the present study revealed only moderate correlation between tumor volume and ADC in larger tumor group suggesting that tumor environment becomes more restrictive as tumor size increases. Moreover, a study by Muraoka et al. demonstrated significant differences in ADC between areas of sparse and dense fibrotic area [8]. Furthermore, ADC is a combined measure of thermally driven molecular movement of water i.e. diffusion and microcirculation of blood in capillaries i.e. perfusion [42] as demonstrated by intravoxel incoherent motion (IVIM) model. Indeed, the IVIM model has been used to study various cancer types [31, 32]. The signal from blood flow is rapidly attenuated at low b values (b < 100–150 s/mm2), whereas higher b values are required to suppress the perfusion contribution [43, 44].

One of the most sensitive measures to monitor the tumor progression, in our animal model, was MTR (%). Previous studies have demonstrated that the transgenic KPC mouse model has the highest degree of fibrotic tissue accumulation as compared to other mouse models evaluated [45]. Furthermore, studies have demonstrated that MTR values are sensitive to fibrotic tissue deposition and have been used to study liver fibrosis and PDA [17, 18, 34]. Additionally, a study by Farr et al. has demonstrated a significant correlation between MTR and fibrotic tissue deposition in the KPC mouse model [18]. Furthermore, using a xenograft mouse model, Li et al. have demonstrated a significant correlation between MTR and fibrotic tissue deposition in different cell lines when the tumor size was close to 10 mm [17]. Similarly, in our study, we have demonstrated a significant correlation between MTR and tumor volume especially when the tumor volume exceeded 250 mm3. Whereas when the tumor volume was less than 250 mm3, we did not find any correlation between MTR and tumor volume. Additionally, we found that there was a moderate correlation between an increase in tumor volume and fibrotic tissue deposition and a strong correlation between MTR and accumulation of fibrotic tissue suggesting that with an increase in tumor volume there is increase in fibrotic tissue deposition and that MTR changes are sensitive to accumulation of fibrotic tissue in the tumor. Therefore, MTR may be a valuable measure in evaluation of novel treatments and more importantly, in determining the stage of cancer and planning the treatment regimen accordingly.

Finally, APT imaging has been used to measure the concentration of endogenous mobile proteins and peptides which are increased in high-grade brain tumors compared with low-grade tumors [46]. We, however, did not find any significant increase in APT signal intensity when comparing tumor in small and large tumor groups. Surprisingly, in the smaller tumor group, we found a significant decrease in the amide proton group with an increase in tumor volume suggesting that these protons have an important role during early stages of tumor development. Once a threshold value is reached, the concentration of amide protons remained unchanged, suggesting a role of other chemical constituents in higher tumor volumes. Indeed, studies of PDA have suggested that there is increased accumulation of hyaluronic acid (HA) rather than an increase in tumor cellularity, which may be one of the major factors leading to an increase in tumor volume [47, 48].

The multi-parametric approach has been utilized before to characterize various tumors and therapeutic effects [18, 49]. To the best of our knowledge this is the first study to characterize the KPC mouse model using a host of MR measurements. Our study has a few limitations that need to be acknowledged. First, due to the dynamic nature of the tumor microenvironment, a host of pathological processes such as necrosis and fibrotic tissue accumulation occur simultaneously with tumor development and future studies need to be conducted in KPC mice at various stages of tumor development. Additionally, MR results need to be correlated with histological measurements in detail. In this study, we did not validate the MR measures with the pathological events at different stages of tumor development. Second, we did not quantify HA accumulation, which correlates with high tumor interstitial fluid pressure (IFP) which in turn collapses the surrounding capillaries. Future studies looking at these features and corroboration with histological measurements are warranted.

Conclusions

We believe that quantitative mp-MRI has a potential role in the monitoring of disease progression and therapeutic evaluation of tumors. In our study, we have demonstrated that MR parameters such as T1 and APT are sensitive to changes in the tumor microenvironment during the early stages of tumor development whereas parameters such as T2, MTR and ADC are sensitive to pathology during the later stages of tumor development. Additionally, our multi-parametric data suggests that changes in advanced MR techniques such as MTR, ADC and APT imaging have the potential to be used in both preclinical and eventually in clinical models to document the underlying pathophysiological processes and thereby initiating tumor targeting therapy.

Abbreviations

APT:

amide proton transfer

AUC:

area under curve

CEST:

chemical exchange saturation transfer

DWI:

diffusion weighted imaging

MTR:

magnetization transfer ratio

PDA:

pancreatic ductal adenocarcinoma

References

  1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67(1):7–30.

    Article  Google Scholar 

  2. Yang S, Wang X, Contino G, Liesa M, Sahin E, Ying H, Bause A, Li Y, Stommel JM, Dell'antonio G, Mautner J, Tonon G, Haigis M, Shirihai OS, Doglioni C, Bardeesy N, Kimmelman AC. Pancreatic cancers require autophagy for tumor growth. Genes Dev. 2011;25(7):717–29.

    Article  CAS  Google Scholar 

  3. Grimm J, Potthast A, Wunder A, Moore A. Magnetic resonance imaging of the pancreas and pancreatic tumors in a mouse orthotopic model of human cancer. Int J Cancer. 2003;106(5):806–11.

    Article  CAS  Google Scholar 

  4. Cui JH, Kruger U, Vogel I, Luttges J, Henne-Bruns D, Kremer B, Kalthoff H. Intact tissue of gastrointestinal cancer specimen orthotopically transplanted into nude mice. Hepatogastroenterology. 1998;45(24):2087–96.

    CAS  PubMed  Google Scholar 

  5. Tan MH, Chu TM. Characterization of the tumorigenic and metastatic properties of a human pancreatic tumor cell line (AsPC-1) implanted orthotopically into nude mice. Tumour Biol. 1985;6(1):89–98.

    CAS  PubMed  Google Scholar 

  6. Hingorani SR, Wang L, Multani AS, Combs C, Deramaudt TB, Hruban RH, Rustgi AK, Chang S, Tuveson DA. Trp53R172H and KrasG12D cooperate to promote chromosomal instability and widely metastatic pancreatic ductal adenocarcinoma in mice. Cancer Cell. 2005;7(5):469–83.

    Article  CAS  Google Scholar 

  7. Hruban RH, Adsay NV, Albores-Saavedra J, Anver MR, Biankin AV, Boivin GP, Furth EE, Furukawa T, Klein A, Klimstra DS, Kloppel G, Lauwers GY, Longnecker DS, Luttges J, Maitra A, Offerhaus GJ, Perez-Gallego L, Redston M, Tuveson DA. Pathology of genetically engineered mouse models of pancreatic exocrine cancer: consensus report and recommendations. Cancer Res. 2006;66(1):95–106.

    Article  CAS  Google Scholar 

  8. Muraoka N, Uematsu H, Kimura H, Imamura Y, Fujiwara Y, Murakami M, Yamaguchi A, Itoh H. Apparent diffusion coefficient in pancreatic cancer: characterization and histopathological correlations. J Magn Reson Imaging. 2008;27(6):1302–8.

    Article  Google Scholar 

  9. Ansari C, Tikhomirov GA, Hong SH, Falconer RA, Loadman PM, Gill JH, Castaneda R, Hazard FK, Tong L, Lenkov OD, Felsher DW, Rao J, Daldrup-Link HE. Development of novel tumor-targeted theranostic nanoparticles activated by membrane-type matrix metalloproteinases for combined cancer magnetic resonance imaging and therapy. Small. 2014;10(3):566–75 417.

    Article  CAS  Google Scholar 

  10. Baron P, Deckers R, Knuttel FM, Bartels LW. T1 and T2 temperature dependence of female human breast adipose tissue at 1.5 T: groundwork for monitoring thermal therapies in the breast. NMR Biomed. 2015;28(11):1463–70.

    Article  Google Scholar 

  11. Hectors SJ, Jacobs I, Heijman E, Keupp J, Berben M, Strijkers GJ, Grull H, Nicolay K. Multiparametric MRI analysis for the evaluation of MR-guided high intensity focused ultrasound tumor treatment. NMR Biomed. 2015;28(9):1125–40.

    Article  Google Scholar 

  12. Shami VM, Mahajan A, Loch MM, Stella AC, Northup PG, White GE, Brock AS, Srinivasan I, de Lange EE, Kahaleh M. Comparison between endoscopic ultrasound and magnetic resonance imaging for the staging of pancreatic cancer. Pancreas. 2011;40(4):567–70.

    Article  Google Scholar 

  13. Kinney T. Evidence-based imaging of pancreatic malignancies. Surg Clin North Am. 2010;90(2):235–49.

    Article  Google Scholar 

  14. Ma C, Li Y, Wang L, Wang Y, Zhang Y, Wang H, Chen S, Lu J. Intravoxel incoherent motion DWI of the pancreatic adenocarcinomas: monoexponential and biexponential apparent diffusion parameters and histopathological correlations. Cancer Imaging. 2017;17(1):12.

    Article  Google Scholar 

  15. De Robertis R, Tinazzi Martini P, Demozzi E, Dal Corso F, Bassi C, Pederzoli P, D'Onofrio M. Diffusion-weighted imaging of pancreatic cancer. World J Radiol. 2015;7(10):319–28.

    Article  Google Scholar 

  16. Li W, Zhang Z, Nicolai J, Yang GY, Omary RA, Larson AC. Quantitative magnetization transfer MRI of desmoplasia in pancreatic ductal adenocarcinoma xenografts. NMR Biomed. 2013;26(12):1688–95.

    Article  Google Scholar 

  17. Li W, Zhang Z, Nicolai J, Yang GY, Omary RA, Larson AC. Magnetization transfer MRI in pancreatic cancer xenograft models. Magn Reson Med. 2012;68(4):1291–7.

    Article  Google Scholar 

  18. Farr N, Wang YN, D'Andrea S, Gravelle KM, Hwang JH, Lee D. Noninvasive characterization of pancreatic tumor mouse models using magnetic resonance imaging. Cancer Med. 2017;6(5):1082–90.

    Article  Google Scholar 

  19. Sagiyama K, Mashimo T, Togao O, Vemireddy V, Hatanpaa KJ, Maher EA, Mickey BE, Pan E, Sherry AD, Bachoo RM, Takahashi M. In vivo chemical exchange saturation transfer imaging allows early detection of a therapeutic response in glioblastoma. Proc Natl Acad Sci U S A. 2014;111(12):4542–7.

    Article  CAS  Google Scholar 

  20. Zhou JY, Payen JF, Wilson DA, Traystman RJ, van Zijl PCM. Using the amide proton signals of intracellular proteins and peptides to detect pH effects in MRI. Nat Med. 2003;9(8):1085–90.

    Article  CAS  Google Scholar 

  21. Watanabe H, Kanematsu M, Tanaka K, Osada S, Tomita H, Hara A, Goshima S, Kondo H, Kawada H, Noda Y, Tanahashi Y, Kawai N, Yoshida K, Moriyama N. Fibrosis and postoperative fistula of the pancreas: correlation with MR imaging findings--preliminary results. Radiology. 2014;270(3):791–9.

    Article  Google Scholar 

  22. Hwang I, Lee JM, Lee KB, Yoon JH, Kiefer B, Han JK, Choi BI. Hepatic steatosis in living liver donor candidates: preoperative assessment by using breath-hold triple-echo MR imaging and 1H MR spectroscopy. Radiology. 2014;271(3):730–8.

    Article  Google Scholar 

  23. Taouli B, Tolia AJ, Losada M, Babb JS, Chan ES, Bannan MA, Tobias H. Diffusion-weighted MRI for quantification of liver fibrosis: preliminary experience. AJR Am J Roentgenol. 2007;189(4):799–806.

    Article  Google Scholar 

  24. Kim M, Gillen J, Landman BA, Zhou J, van Zijl PC. Water saturation shift referencing (WASSR) for chemical exchange saturation transfer (CEST) experiments. Magn Reson Med. 2009;61(6):1441–50.

    Article  Google Scholar 

  25. Jin T, Wang P, Zong X. Kim SG. MR imaging of the amide-proton transfer effect and the pH-insensitive nuclear overhauser effect at 9.4 T. Magn Reson Med. 2013;69(3):760–70.

    Article  CAS  Google Scholar 

  26. Desmond KL, Moosvi F, Stanisz GJ. Mapping of amide, amine, and aliphatic peaks in the CEST spectra of murine xenografts at 7 T. Magn Reson Med. 2014;71(5):1841–53.

    Article  Google Scholar 

  27. Kim S, Decarlo L, Cho GY, Jensen JH, Sodickson DK, Moy L, Formenti S, Schneider RJ, Goldberg JD, Sigmund EE. Interstitial fluid pressure correlates with intravoxel incoherent motion imaging metrics in a mouse mammary carcinoma model. NMR Biomed. 2012;25(5):787–94.

    Article  Google Scholar 

  28. Chen L, Ma X, Li H, Gu L, Li X, Gao Y, Xie Y, Zhang X. Influence of tumor size on oncological outcomes of pathological T3aN0M0 renal cell carcinoma treated by radical nephrectomy. PLoS One. 2017;12(3):e0173953.

    Article  Google Scholar 

  29. Gao Z, Wang C, Xue Q, Wang J, Shen Z, Jiang K, Shen K, Liang B, Yang X, Xie Q, Wang S, Ye Y. The cut-off value of tumor size and appropriate timing of follow-up for management of minimal EUS-suspected gastric gastrointestinal stromal tumors. BMC Gastroenterol. 2017;17(1):8.

    Article  Google Scholar 

  30. Su X, Fang D, Li X, Xiong G, Zhang L, Hao H, Gong Y, Zhang Z, Zhou L. The influence of tumor size on oncologic outcomes for patients with upper tract urothelial carcinoma after radical Nephroureterectomy. Biomed Res Int. 2016;2016:4368943.

    PubMed  PubMed Central  Google Scholar 

  31. Lemke A, Laun FB, Klauss M, Re TJ, Simon D, Delorme S, Schad LR, Stieltjes B. Differentiation of pancreas carcinoma from healthy pancreatic tissue using multiple b-values: comparison of apparent diffusion coefficient and intravoxel incoherent motion derived parameters. Investig Radiol. 2009;44(12):769–75.

    Article  Google Scholar 

  32. Le Bihan D, Breton E, Lallemand D, Aubin ML, Vignaud J, Laval-Jeantet M. Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology. 1988;168(2):497–505.

    Article  CAS  Google Scholar 

  33. Liu X, Johnson S, Liu S, Kanojia D, Yue W, Singh UP, Wang Q, Wang Q, Nie Q, Chen H. Nonlinear growth kinetics of breast cancer stem cells: implications for cancer stem cell targeted therapy. Sci Rep. 2013;3:2473.

    Article  Google Scholar 

  34. Aisen AM, Doi K, Swanson SD. Detection of liver fibrosis with magnetic cross-relaxation. Magn Reson Med. 1994;31(5):551–6.

    Article  CAS  Google Scholar 

  35. Beall PT, Asch BB, Chang DC, Medina D, Hazlewood CF. Distinction of normal, preneoplastic, and neoplastic mouse mammary primary cell cultures by water nuclear magnetic resonance relaxation times. J Natl Cancer Inst. 1980;64(2):335–8.

    Article  CAS  Google Scholar 

  36. Yin T, Peeters R, Feng Y, Liu Y, Yu J, Dymarkowski S, Himmelreich U, Oyen R, Ni Y. Characterization of a rat orthotopic pancreatic head tumor model using three-dimensional and quantitative multi-parametric MRI. NMR Biomed. 2017;30(2).

    Article  Google Scholar 

  37. Gullino PM. Considerations on blood supply and fluid exchange in tumors. Prog Clin Biol Res. 1982;107:1–20.

    CAS  PubMed  Google Scholar 

  38. Rofstad EK, Steinsland E, Kaalhus O, Chang YB, Hovik B, Lyng H. Magnetic resonance imaging of human melanoma xenografts in vivo: proton spin-lattice and spin-spin relaxation times versus fractional tumour water content and fraction of necrotic tumour tissue. Int J Radiat Biol. 1994;65(3):387–401.

    Article  CAS  Google Scholar 

  39. Li W, Hong L, Hu L, Magin RL. Magnetization transfer imaging provides a quantitative measure of chondrogenic differentiation and tissue development. Tissue Eng Part C Methods. 2010;16(6):1407–15.

    Article  CAS  Google Scholar 

  40. Le Bihan D, Delannoy J, Levin RL. Temperature mapping with MR imaging of molecular diffusion: application to hyperthermia. Radiology. 1989;171(3):853–7.

    Article  CAS  Google Scholar 

  41. Latour LL, Svoboda K, Mitra PP, Sotak CH. Time-dependent diffusion of water in a biological model system. Proc Natl Acad Sci U S A. 1994;91(4):1229–33.

    Article  CAS  Google Scholar 

  42. Le Bihan D, Breton E, Lallemand D, Grenier P, Cabanis E. Laval-Jeantet M. MR imaging of intravoxel incoherent motions: application to diffusion and perfusion in neurologic disorders. Radiology. 1986;161(2):401–7.

    Article  CAS  Google Scholar 

  43. Le Bihan D, Turner R. The capillary network: a link between IVIM and classical perfusion. Magn Reson Med. 1992;27(1):171–8.

    Article  CAS  Google Scholar 

  44. Le Bihan D, Breton E, Lallemand D, Desbleds MT, Aubin ML, Vignaud J, Roger B. Contribution of intravoxel incoherent motion (IVIM) imaging to neuroradiology. J Neuroradiol. 1987;14(4):295–312.

    CAS  PubMed  Google Scholar 

  45. Olive KP, Jacobetz MA, Davidson CJ, Gopinathan A, McIntyre D, Honess D, Madhu B, Goldgraben MA, Caldwell ME, Allard D, Frese KK, Denicola G, Feig C, Combs C, Winter SP, Ireland-Zecchini H, Reichelt S, Howat WJ, Chang A, Dhara M, Wang L, Ruckert F, Grutzmann R, Pilarsky C, Izeradjene K, Hingorani SR, Huang P, Davies SE, Plunkett W, Egorin M, Hruban RH, Whitebread N, McGovern K, Adams J, Iacobuzio-Donahue C, Griffiths J, Tuveson DA. Inhibition of hedgehog signaling enhances delivery of chemotherapy in a mouse model of pancreatic cancer. Science. 2009;324(5933):1457–61.

    Article  CAS  Google Scholar 

  46. Wen Z, Hu S, Huang F, Wang X, Guo L, Quan X, Wang S, Zhou J. MR imaging of high-grade brain tumors using endogenous protein and peptide-based contrast. Neuroimage. 2010;51(2):616–22.

    Article  Google Scholar 

  47. DuFort CC, DelGiorno KE, Carlson MA, Osgood RJ, Zhao C, Huang Z, Thompson CB, Connor RJ, Thanos CD, Scott Brockenbrough J, Provenzano PP, Frost GI, Michael Shepard H, Hingorani SR. Interstitial pressure in pancreatic ductal adenocarcinoma is dominated by a gel-fluid phase. Biophys J. 2016;110(9):2106–19.

    Article  CAS  Google Scholar 

  48. Provenzano PP, Cuevas C, Chang AE, Goel VK, Von Hoff DD, Hingorani SR. Enzymatic targeting of the stroma ablates physical barriers to treatment of pancreatic ductal adenocarcinoma. Cancer Cell. 2012;21(3):418–29.

    Article  CAS  Google Scholar 

  49. Hectors SJ, Jacobs I, Strijkers GJ, Nicolay K. Multiparametric MRI analysis for the identification of high intensity focused ultrasound-treated tumor tissue. PLoS One. 2014;9(6):e99936.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work is supported by NIH R01CA188654 and NIH R01CA154451.

Availability of data and materials

All data generated during this study are included in this published article.

Author information

Authors and Affiliations

Authors

Contributions

DL contributed the study concept, and RV and DL designed the study. Data acquisition was performed by RV, JP, YNW, KG, SW. RV performed the statistical analysis. All authors contributed to data analysis, interpretation and approval of final manuscript.

Corresponding author

Correspondence to Donghoon Lee.

Ethics declarations

Ethics approval and consent to participate

All experiments were performed in accordance with the guidelines for the care and use of laboratory animal of the national institutes of health and with approval from our institutional animal care and use committee.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Publisher’s Note

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

Additional file

Additional file 1:

Figure S1. Time course multiparametric MR images with colorized MR maps overlaid on T2 weighted anatomic images obtained at weeks 1, 2, 3 and 4. (TIFF 1522 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Vohra, R., Park, J., Wang, YN. et al. Evaluation of pancreatic tumor development in KPC mice using multi-parametric MRI. Cancer Imaging 18, 41 (2018). https://doi.org/10.1186/s40644-018-0172-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40644-018-0172-6

Keywords