Next Article in Journal
Dynamic Contrast Enhanced-MR CEST Urography: An Emerging Tool in the Diagnosis and Management of Upper Urinary Tract Obstruction
Previous Article in Journal
Publisher’s Note: Continued Publication of Tomography by MDPI
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Respiratory Motion Mitigation and Repeatability of Two Diffusion-Weighted MRI Methods Applied to a Murine Model of Spontaneous Pancreatic Cancer

1
Department of Radiology, University of Pennsylvania, Philadelphia, PA 19104, USA
2
Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania, Philadelphia, PA 19104, USA
3
Pancreatic Cancer Research Center, University of Pennsylvania, Philadelphia, PA 19104, USA
4
Abramson Cancer Center, University of Pennsylvania, Philadelphia, PA 19104, USA
*
Author to whom correspondence should be addressed.
Submission received: 12 December 2020 / Accepted: 2 February 2021 / Published: 20 February 2021

Abstract

:
Respiratory motion and increased susceptibility effects at high magnetic fields pose challenges for quantitative diffusion-weighted MRI (DWI) of a mouse abdomen on preclinical MRI systems. We demonstrate the first application of radial k-space-sampled (RAD) DWI of a mouse abdomen using a genetically engineered mouse model of pancreatic ductal adenocarcinoma (PDAC) on a 4.7 T preclinical scanner equipped with moderate gradient capability. RAD DWI was compared with the echo-planar imaging (EPI)-based DWI method with similar voxel volumes and acquisition times over a wide range of b-values (0.64, 535, 1071, 1478, and 2141 mm2/s). The repeatability metrics are assessed in a rigorous test–retest study (n = 10 for each DWI protocol). The four-shot EPI DWI protocol leads to higher signal-to-noise ratio (SNR) in diffusion-weighted images with persisting ghosting artifacts, whereas the RAD DWI protocol produces relatively artifact-free images over all b-values examined. Despite different degrees of motion mitigation, both RAD DWI and EPI DWI allow parametric maps of apparent diffusion coefficients (ADC) to be produced, and the ADC of the PDAC tumor estimated by the two methods are 1.3 ± 0.24 and 1.5 ± 0.28 × 10−3 mm2/s, respectively (p = 0.075, n = 10), and those of a water phantom are 3.2 ± 0.29 and 2.8 ± 0.15 × 10−3 mm2/s, respectively (p = 0.001, n = 10). Bland-Altman plots and probability density function reveal good repeatability for both protocols, whose repeatability metrics do not differ significantly. In conclusion, RAD DWI enables a more effective respiratory motion mitigation but lower SNR, while the performance of EPI DWI is expected to improve with more advanced gradient hardware.

1. Introduction

Genetically engineered mouse (GEM) models of pancreatic ductal adenocarcinoma (PDAC), in which the cancer arises spontaneously in pancreas, [1,2] recapitulate key features of the human disease, including a dense extracellular matrix (stroma), which represents a unique tumor microenvironment in PDAC. These models have become a major tool for the development of biomarkers and treatment pipelines [3,4,5,6,7,8]. Diffusion-weighted MR imaging (DWI) provides quantitative metrics related to the translational mobility of water hindered by microstructures present in biological tissues [9,10]. Due to its sensitivity to changes in the tumor microenvironment and cell density, DWI has been frequently employed in oncology for the diagnosis and detection of treatment responses [11,12,13,14]. Since DWI pulse sequences are sensitized to motion on the micrometer scale, macroscopic (millimeter) scale movement such as respiratory motion can introduce artefacts, e.g., ghosting, resulting in errors in the quantitative measurement of water diffusion in tissues. Located in the retroperitoneum, the pancreas is presumably less susceptible to respiratory motion than other upper abdominal organs, such as the liver. However, 4D CT and 4D MRI studies have shown that the amplitude of pancreas movement during breathing is comparable to that of the liver and lungs [15,16,17]. For example, Ferris et al. demonstrated 5–7-mm root-mean-square displacement of the human pancreas [18]. To minimize the respiratory motion artifacts, clinical DWI employs rapid acquisition protocols, such as single-shot echo-planar imaging (EPI) following diffusion weighting preparation [19]. However, unique challenges arise in the reverse translation of clinical EPI DWI techniques to the preclinical instrument: (1) a higher magnetic field strength leads to proportionally greater magnetic susceptibility effects [20], (2) faster respiration rates of mice exacerbate motion interference, (3) less advanced gradient hardware on preclinical scanners severely limits the echo train length that can be employed in EPI acquisitions, and (4) the small dimensions of the mouse limit the number of elements in array coils [21]; hence, the acceleration factor is usually less than two. These factors combined lead to the requirement of using the multi-shot EPI method in DWI studies or suboptimal image quality in single-shot EPI DWI of mouse abdomens on preclinical scanners [22,23,24,25].
Physical restrictions of respiratory motion, such as casting live mice in alginate-molds [26] or applying a rigid structure [27], appear to be effective for suppressing motion artifacts. However, concerns of respiratory distress in compromised mice may limit their application. New acquisition methods, including spatiotemporal encoding (SPEN) [28], motion-tolerant overlapping-echo detachment (DM-OLED) [29], and the use of intermolecular dipole–dipole interactions [30] have been proposed to minimize the sensitivity to motion, but in each case, motion robustness comes at the expense of other image quality factors, such as reduced spatial resolution, limited b-values, or low SNR. In comparison, radial k-space sampling (RAD) is inherently insensitive to motion due to the over-sampling of the k-space center [31] and does not trade image quality for motion robustness. On clinical scanners, dynamic contrast-enhanced (DCE)-MRI and DWI studies of patients with liver, lung, breast, and prostate cancers benefit from radial k-space sampling for the effective mitigation of respiratory motion interference [32,33,34,35].
RAD DWI has not been investigated in murine models of orthotopically implanted or spontaneous pancreatic cancer or in mouse abdominal organs in general. We report here the first such study in a GEM model of pancreatic cancer in comparison with a multi-shot EPI DWI protocol. We compared the SNR and image quality of the two protocols, i.e., for producing artifact-free diffusion-weighted images of pancreatic cancer in situ and their accuracy in estimating the diffusion coefficients of phantoms. In addition, we performed a rigorous test–retest study as recommended by the quantitative imaging biomarkers alliance (QIBA) [36] and reported the repeatability metrics of the apparent diffusion coefficient (ADC). These data provided reference ADC values of the PDAC in a GEM model and the repeatability metrics of two translational protocols suitable for the DWI of abdominal organs/cancers in freely breathing mice.

2. Methodology

2.1. GEM Model of Pancreatic Cancer

All animal handling procedures were reviewed and approved by the local Institutional Animal Care and Use Committee. We employed a GEM model of PDAC, referred to as KPC mice, in which Kras (“K”) and p53 (“P”) mutations were introduced to the pancreatic epithelium via Cre recombinase (“C”) [1]. The model was bred at the Pancreatic Cancer Mouse Hospital of the Abramson Cancer Center of our institute. KPC mice spontaneously developed premalignant Pancreatic Intraepithelial Neoplastic (“PanIN”) lesions at 7–10 weeks of age, leading to invasive PDAC at 17–19 weeks with high penetrance and reproducible kinetics. Screening for tumors was done via weekly abdominal palpations starting at 11 weeks of age, followed by ultrasound examination (Vevo 2100, VisualSonics, Toronto, ON, Canada) to estimate the tumor sizes. Mice with confirmed tumor masses (both sexes, 18–25 weeks old) were enrolled in the MRI study.

2.2. In Vivo MRI and Test-Retest Study

MRI was performed on a 4.7 T horizontal bore system (Agilent, Palo Alto, CA, USA) equipped with a 12-cm, 25 gauss/cm gradient coil with a maximum slew rate of 11.5 T/cm/s. Animals were anesthetized via free breathing of 1.5–3% isoflurane in oxygen administered through a nose cone. Probes for monitoring the respiration rate and rectal temperature (SAI Inc., Stonybrook, NY, USA) were placed, and each mouse was positioned on a bed that was mounted into a 35 × 40 mm (ID × length) quadrature birdcage RF coil (M2M Imaging, Cleveland, OH, USA). Core body temperature was maintained at 37 ± 0.2 °C during MRI exams by directing a regulated air source into the bore. Respiration rate was maintained in the range of 80–100 breaths/min during imaging by adjustment of the isoflurane level. Two sealed 5-mm NMR tubes (Wilmad Glass, Vineland, NJ, USA) containing water and pure 1-butanol, respectively, were placed underneath the mice in the slotted mouse beds.
Following scanner calibrations and the generation of scout images, contiguous, axial, T2 weighted images spanning the tumor were acquired using a fast spin-echo protocol: effective echo time (TE) = 45 ms, repetition time (TR) = 1.5 s, slice thickness = 0.5 mm, field of view (FOV) = 32 × 32 mm2, slices = 12–20, matrix = 128 × 128, averages = 2, bandwidth = 50 kHz, echo train length = 8, and echo spacing = 15 ms. The radially sampled diffusion-weighted spin-echo protocol (DW-SE-RAD) and diffusion-weighted spin-echo 4-shot EPI (DW-SE-EPI) protocol, both with prospective respiratory gating, were then applied. For DW-SE-RAD: TE = 30 ms, effective TR = 2 to 3 respiratory periods (0.6–0.75 s per respiratory period) based on the number of slices expanding the entire tumor, one radial view (spanning the diameter of the k-space) acquired per TR, FOV = 32 × 32 mm2, slice thickness = 0.5 mm, slices = 12–18, matrix = 64 read × 101 views, averages = 1, bandwidth = 50 kHz, and total acquisition time = 20–38 min, depending on the respiratory rate and effective TR. For DW-SE-EPI: 4-shot EPI with an interleaved phase encoding the trajectory for each shot, TE = 26 ms, TR = 4 respiration periods, same FOV, slice thickness = 1.5 mm, slices = 4–6, matrix = 128 × 128, averages = 16, bandwidth = 250 kHz, and total acquisition time = 20–35 min. The slice thickness and TR employed for EPI represent the minimal values allowed by the gradient hardware. The remaining acquisition parameters were designed to keep the voxel volumes and total acquisition times of the two protocols similar to facilitate a fair comparison. Both DWI protocols employed identical diffusion times (Δ = 14.4 ms); diffusion gradient durations (δ = 9 ms); and five b-values (0.64, 535, 1071, 1478, and 2141 mm2/s). The minimum b-value was selected to minimize the contributions from the intravoxel incoherent motion (micro perfusion) [37]. All diffusion gradients were oriented along the diagonal, gdiff = (1:1:1) to minimize the contributions from intermolecular dipole–dipole interactions [38] to the observed signals. Images were acquired with both positive and negative (-gdiff) diffusion gradients, resulting in a total of 10 images per slice location. Fat saturation was not applied in either protocol to avoid the suppression of 1-butanol 1H signals arising in the aliphatic region (0.9–1.5 ppm) of its proton spectrum.
In clinical studies, the RAD and EPI acquisition methods have been shown to be robust to respiration motion, even in the absence of respiratory gating [32]. The need for such gating was assessed in two tumor-bearing mice in which respiratory-gated and non-gated data were acquired sequentially from each mouse. In the non-gated scans, the TR was set to mimic the effective TR used in the gated scans (e.g., 2 to 3 respiratory periods for RAD).
In the test–retest study, eight mice received DW-SE-RAD and DW-SE-EPI applied sequentially to each mouse during the “test” and “retest” scans. After the test scan, the mice were recovered from anesthesia and returned to the cage to have free access to water and food for 2–4 h. All procedures were repeated during the retest scan. Two additional mice received only RAD DWI for the test and retest scans, and another two received only EPI DWI. Hence, for each DWI protocol, test and retest data were obtained from ten mice (n = 10).

2.3. Data Processing

EPI data were phase-corrected and transformed to images using the vendor software. Radial k-space data were zero- and first-order phase-corrected and regridded to 64 × 64 images offline via a filtered back-projection procedure. Images acquired with positive and negative diffusion gradients were combined by taking the geometric mean of the two in order to minimize the contributions of background and imaging gradients to the signal attenuation [39]. The resulting images were subject to pixel-wise three-parameter least squares fit to Equation (1) using all five b-values:
S ( b ) = S 0 e x p { b D } + A
where S0, D, and A are the adjustable parameters. S0 is the signal intensity in the absence of diffusion weighting, D is the ADC, and A is a constant that accounts for the noise floor and slowly diffusing water contributions to the observed signal. Maps of these parameters were generated for further analysis. All offline data processing was performed using custom codes developed in Interactive Data Language (IDL; Harris Geospatial Solutions, Broomfield, CO, USA). Regions of interest (ROI) representing the tumor, spinal muscle, and two phantoms were manually drawn by trained staff (JC and SP) using ImageJ (National Institute of Health, Bethesda, MD, USA) [40]. Mean ADC values were obtained over the entire tumor volume and for the phantom and spinal muscles over all slices containing the tumor.

2.4. Statistical Analysis

To determine the correlation of ADC obtained from the two DWI protocols, we first applied the generalized linear model (GLM) to the test and the retest data separately to compute the ADC correlation coefficient (R2) between the two protocols. In eight mice, each received both DWI protocols in both the test and retest scans; the combined test and retest data were not suitable for a linear model, because they were not independent. To assess the correlation in the combined data, we employed the generalized estimating equations (GEE) approach to assess the linear regression using the mean ADC of each mouse estimated from the RAD protocol as the response variable and that from the EPI protocol as the independent variable in statistics software R version 3.6.3 (www.r-project.org) (28 May 2020) (R Foundation for Statistical Computing, Vienna, Austria).
For each protocol, Bland-Altman plots were employed to assess the bias and precision of ADC from the test–retest study. Within-subject standard deviation (SDws) and coefficient of variation (CVWS) were calculated as follows:
S D w s   =   [ (   Δ D 2 ) / 2 n ]
C V W S   =   [   ( Δ D m ) 2 / 2 n ]
where ΔD is the difference between the test–retest ADC estimates for a given animal, m is the mean of the measurements for a given pair, n is the number of paired test–retest studies (n = 10) for each DWI protocol, and the sum is taken over n. The repeatability coefficient (RC) was calculated at a confidence level of 95% [36]:
R C =   1.96   2 S D w s 2   =   2.77   ×   S D w s
To illustrate the ADC heterogeneity within a tumor and to compare the pixel-wise ADC distributions between the test and retest scans, plots of the probability density function (PDF) were made using R software.

3. Results

Images from the DW-SE-EPI protocol without prospective gating were motion-corrupted, and meaningful ADC maps could not be generated (data not shown). In comparison, although non-respiratory-gated DW-SE-RAD images exhibit increased noise and streaking artifacts compared to their gated counterparts (Figure 1), both allow the construction of good quality ADC maps, confirming the motion-insensitive benefit of radial k-space sampling. Higher levels of artifacts and potential signal loss from motion during the diffusion time seem to lead to faster apparent signal decay and increased ADC (1.1 ± 0.097 vs. 1.3 ± 0.17 × 10−3 mm2/s in tumors and 3.2 ± 0.13 vs. 4.0 ± 0.30 in the water phantom for gated and non-gated, respectively). These results demonstrate that both RAD and EPI protocols would benefit from prospective gating, which was applied in the subsequent studies presented below.
Representative images generated by the two DWI protocols with prospective gating are depicted in Figure 2. In rows A (RAD) and B (EPI), frames 1–5 are windowed to provide the optimal visualization of the attenuation due to the different levels of diffusion weighting. The signal attenuation with increasing b-values is comparable between the two protocols. The anatomy and tumor are clearly visualized in both datasets. At similar total acquisition times, the EPI protocol has a higher SNR (in the lowest b-value image) of 150 for the water phantom compared to 15 with the RAD protocol. This is due, in large part, to the higher efficiency and longer TR of the EPI acquisition protocol. The average SNRs in tumors are 119 and 20 for EPI and RAD, respectively, a smaller relative difference than that of water, likely due to the shorter T1 in tumors. Geometric distortion in the EPI images was minimal due to optimal shimming and the relatively short echo train employed. Panels Figure 2(A5’,B5’) are the same images as Figure 2(A5,B5), respectively but are windowed to highlight the image artifacts. In the EPI images, ghosting is evident in the phase-encoding direction (left to right) at all b-values and has a comparable magnitude to the desired signal in the high b-value images. Such ghosting is likely due to shot-to-shot phase errors [41] induced by respiration motion. In contrast, no coherent artifacts are evident in the DW-SE-RAD images. The corresponding ADC maps are shown in Figure 2(A6,B6) in color for the tumor and grayscale for the other tissues and phantoms.
The signal from the water phantom (blue arrows) drops to noise levels for the b-values over 1000 mm2/s for both protocols, consistent with the fact that the diffusion coefficient of free water is much higher than those of most biological tissues; therefore, a 1-butanol phantom was included to better match the tissue ADC values [42]. While, in DW-SE-RAD images, the signal from 1-butanol (yellow arrow) remains above the noise at the largest b-value, it is corrupted at all b-values in the DW-SE-EPI images due to the off-resonance effects. The 1H spectrum of 1-butanol consists of four resonances at 3.51, 1.56, 1,25, and 0.8 ppm, with most of the signal centered around 1.56 ppm. The phase correction scheme employed for the generation of EPI images only compensates for the dominant water peak (4.7 ppm), leading to the observed signal corruption.
Diffusion coefficient estimates for the water phantom had good reproducibility for both protocols with mean ± standard deviations of 3.2 ± 0.29 and 2.8 ± 0.15 × 10−3 mm2/s for the DW-SE-RAD and DW-SE-EPI protocols, respectively, compared to the literature value of 3.20 × 10−3 mm2/s at 37 °C [43]. However, the mean difference was found to be significant (p = 0.001, Table 1). This underestimation with EPI was also observed in phantom-only scans in which no mouse was loaded, excluding the possibility that respiratory motion artifacts may be the cause. Since we did not measure the temperature of the water phantom directly, we cannot be certain that they were at 37 °C, even if the core temperatures of the mice were maintained at that temperature. This could lead to uncertainties in the actual ADC value of the water phantom. The diffusion coefficient of the 1-butanol phantom estimated by the RAD DWI protocol was 0.44 ± 0.22 × 10−3 mm2/s, which compares favorably with the reported values (0.3–0.628 × 10−3 mm2/s at 25–30 °C) [42]; however, such an estimation was not feasible for EPI DWI due to the off-resonance effects.
For the ADC estimates of the PDAC tumor, a strong correlation between the two protocols was obtained from the test (R2 = 0.67, p < 0.05) and the retest (R2 = 0.96, p = 0.0001) datasets, separately. When the test and retest data were combined, the GEE metrics also revealed a statistically significant correlation of the tumor ADC estimates between the DW-SE-RAD vs. DW-SE-EPI methods (p = 1.9 × 10−15) (Figure 3).
The statistical metrics generated from the test–retest study are presented in Table 1 for the water phantom, PDAC tumor, and spinal muscle. The RAD DWI protocol yields smaller SDws, CVWS, and RC of the ADC estimates for both tumor and muscle, suggesting a trend of higher repeatability, although the differences are not statistically significant. For the water phantom, the trend is reversed in favor of EPI DWI. These observations suggest that a more effective motion mitigation may underlie the observed trend in the repeatability metrics.
Bland-Altman plots of ADC values generated from the test–retest study are shown in Figure 4 for the tumor and muscle. For both protocols, the biases for both tissues are very small relative to RC, indicating a minimal variation in the ADC estimates in the test–retest setting. The intra-subject difference of the ADC (ΔADC, y-axis values) does not appear to be correlated with the magnitude of their mean (x-axis values) for the tissues investigated here. The distances between the lines representing the mean ± RC in Figure 4 confirm the trend of higher repeatability provided by the RAD compared to the EPI protocol.
Among the ten pairs of test–retest ADC values from each DWI protocol, the best and worst cases, corresponding to the smallest and largest ΔD (or ΔADC) values, respectively, were selected. The two worst cases are identified by red arrows in Figure 4A,C. The PDFs of the tumor ADC values for these cases are shown in Figure 5. The PDFs (analogous to histograms) reveal the heterogeneity and distribution of the ADC values within the tumor. The PDFs from test and retest scans are nearly superimposed in the best-case examples (Figure 5A,B) but are notably different in the worst-case examples, with larger differences from the EPI DWI protocol (Figure 5C,D).

4. Discussion

Many cancers arise from the abdominal organs, such as the liver, and colorectal cancer. While respiratory motion is known to interfere with the quantitative DWI of the mouse abdomen, this challenge can be mitigated by employing subcutaneous (subQ) models, which are produced by inoculating cancer cells under the skin of the hind limb or flank, since both locations are relatively immune to respiratory motion. Recently, there has been an increased use of GEM models, in which tumors arise from and interact with the native tissue environment and, therefore, overcome some key limitations associated with subQ models [44].
To meet the need for mitigating respiratory motions in preclinical DWI, we demonstrated the first application of radially sampled DWI of the mouse abdomen and made comparisons between the radial and EPI-based methods. While both protocols are implemented on a preclinical scanner equipped with moderate gradient hardware, DW-SE-EPI was more affected by the limited gradient capability, leading to requirements for a longer TR, thicker image slices, and multi-shot EPI acquisition. To compare RAD and EPI DWI performances under such constraints, we closely matched the voxel volume and total acquisition time. Our data showed that the radial method is superior in motion mitigation and accuracy in estimating the diffusion coefficient of the water phantom. In comparison, EPI DWI has a more favorable SNR but leads to suboptimal mitigation of motion artifacts. It is possible, though, that the relative merits of the two DWI methods could be different with improvements in the gradient and RF hardware on preclinical MRI systems. For example, a greater maximum gradient strength may enable a shorter diffusion time (Δ) and bipolar diffusion gradients, both of which appear effective in the suppression of motion artifacts in a recent EPI DWI study of murine liver cancer [22].
The lower SNR observed in the DW-SE-RAD images may potentially have been affected by the dispersion of motion-related artifacts into the background regions, since radial methods avoid coherent artifacts (e.g., ghosting) by spreading out the effects throughout the FOV, leading to apparently elevated background noise. Relatively long T1 values of tumors (and water) could also contribute to reduce the SNR in the radial acquisition due to shorter TRs than the EPI protocol. To improve image SNR, it may be possible to utilize trajectory correction schemes [45] and iterative image reconstruction methods [46] in the DW-SE-RAD protocol.
While 1-butanol may have been sufficient for radial-based methods, its off-resonance characteristics led to undesired artifacts in EPI. Our study could have been improved with the utilization of a water-based phantom with the ADC closely matching that of the tissue of interest, such as a polyvinylpyrrolidone (PVP) phantom [47].
Despite different degrees of motion mitigation, both DW-SE-RAD and DW-SE-EPI allow quantitative DWI of the mouse abdomen by producing parametric maps of the ADC. Tumor ADC values estimated from RAD vs. EPI DWI are highly correlated (Figure 3) and do not differ significantly (Table 1). The ADC values from KPC mice provide a reference for future studies and are well within the range reported from human PDAC (1.11–1.50 × 10−3 mm2/s) [48,49], suggesting similarities in the tumor architecture and microenvironment between the GEM model and human disease.
The repeatability metrics provide variability assessments for each DWI protocol. Such assessments would enable the estimation of the sample size required by each DWI method [36] and are motivated by increased consensus within the cancer imaging community to report these metrics for quantitative imaging markers [14,27,50,51,52]. Considering the highly motion-susceptible location of the pancreatic tumor in freely breathing mice, our reproducibility metrics of the DW-SE-RAD protocol (SDws = 0.12 × 10−3 mm2/s and CVWS = 9%) compare favorably to those reported from tumors located in less motion-susceptible locations—for example, the SDws of 0.1 × 10−3 mm2/s from breast cancer xenografts grown in the hind flank of mice [50], SDws of 0.06 × 10−3 mm2/s, and CVWS of 7% from orthotopically implanted breast tumors in mice with restricted respiration [27].
In this work, we endeavored to keep the voxel volumes and scan times similar between the two protocols in order to facilitate performance comparisons. In this regard, it is possible that alternative choices of imaging parameters may have further improved the comparisons, such as increasing the number of segments or the number of b-values for the EPI protocol rather than simply increasing the number of averages. A greater number of EPI segments would also help reduce the gradient duty cycle. However, it is likely that much of the ghosting artifacts observed in the EPI images are due to shot-to-shot phase variations generated by motion. As such, increases in the number of segments would increase the amplitude and/or spatial frequency of phase-encoding ghosts, subsequently increasing the amount of signal overlap. Further, one of the benefits of averaging is the suppression of ghosting artifacts caused by motion. Our observation with a four-shot EPI with no averaging yielded poor-quality images (data not shown).
While exactly matching the scan parameters between the two sequences may make the comparisons fairer in some sense, one of our goals in designing these protocols was to take advantage of the best features of each method. The main motivation for employing EPI methods is the very high acquisition efficiency associated with the technique. Further reductions in the acquisition efficiency due to an increased number of segments may detract from the advantage of this technique.
One disadvantage of the radial DWI protocol is the relatively long scan time, which may make it impractical for applications that require many b-values and/or multiple diffusion gradient directions. Therefore, minimizing the total acquisition time of a DWI study is important for increasing the throughput and permitting multiple imaging markers to be obtained in the same session. To accelerate the acquisition of RAD DWI, several approaches are worth investigating, including diffusion-weighted radial fast spin-echo [53] and diffusion-weighted PROPELLER [54], which combines the desirable features of both the radial and multi-echo methods. Of note, radial sampling’s ability to tolerate under-sampling can be exploited in combination with new image reconstruction models for under-sampled radial data [55]. Consequently, such an approach may allow reducing the acquisition time while maintaining the good quality of under-sampled images.

Author Contributions

Conceptualization, J.C. (Jianbo Cao), S.P., and R.Z.; methodology development, J.C. (Jianbo Cao), S.P., and R.Z.; software, S.P., and H.K.S.; data analysis, J.C. (Jianbo Cao), J.C. (Jinbo Chen), S.P., and R.Z.; resources, C.C.; data curation, J.C. (Jianbo Cao), H.Y., V.C., and S.P.; writing—original draft preparation, S.P., and R.Z.; writing—review and editing, H.K.S., S.P., M.R., and R.Z.; visualization, J.C. (Jianbo Cao), S.P., and H.K.S.; supervision, R.Z.; funding acquisition, M.R. and R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by from NIH U24-CA-231858 (Penn Quantitative Imaging Resource for Pancreatic Cancer).

Institutional Review Board Statement

All animal handling procedures were reviewed and approved by the Animal Care and Use Committee of University of Pennsylvania (Protocol # 806083 approved since 03/12/2017).

Informed Consent Statement

Not applicable (no humans were involved).

Data Availability Statement

Image data will be shared at a data repository being built for Penn Quantitative Imaging Resource for Pancreatic Cancer (https://pennpancreaticcancerimagingresource.github.io/data.html/).

Acknowledgments

We acknowledge the excellent technical support from the Small Animal Imaging Facility (SAIF) of the Radiology Department at the University of Pennsylvania.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ADC or Dapparent diffusion coefficient
CVWSwithin-subject coefficient of variation
DWIdiffusion-weighted MRI
DW-SE-EPIdiffusion-weighted 4-shot spin-echo echo planar imaging protocol
DW-SE-RADdiffusion-weighted radially sampled spin-echo protocol
GEMgenetically engineered mouse
PDACpancreatic ductal adenocarcinoma
PDFprobability density function
RCrepeatability coefficient
SDwswithin-subject standard deviation

References

  1. Hingorani, S.R.; Petricoin, E.F., III; Maitra, A.; Rajapakse, V.; King, C.; Jacobetz, M.A.; Ross, S.; Conrads, T.P.; Veenstra, T.D.; Hitt, B.A.; et al. Preinvasive and invasive ductal pancreatic cancer and its early detection in the mouse. Cancer Cell 2003, 4, 437–450. [Google Scholar] [CrossRef] [Green Version]
  2. Bardeesy, N.; Cheng, K.H.; Berger, J.H.; Chu, G.C.; Pahler, J.; Olson, P.; Hezel, A.F.; Horner, J.; Lauwers, G.Y.; Hanahan, D.; et al. Smad4 is dispensable for normal pancreas development yet critical in progression and tumor biology of pancreas cancer. Genes Dev. 2006, 20, 3130–3146. [Google Scholar] [CrossRef] [Green Version]
  3. Frese, K.K.; Neesse, A.; Cook, N.; Bapiro, T.E.; Lolkema, M.P.; Jodrell, D.I.; Tuveson, D.A. Nab-paclitaxel potentiates gemcitabine activity by reducing cytidine deaminase levels in a mouse model of pancreatic cancer. Cancer Discov. 2012, 2, 260–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Jacobetz, M.A.; Chan, D.S.; Neesse, A.; Bapiro, T.E.; Cook, N.; Frese, K.K.; Feig, C.; Nakagawa, T.; Caldwell, M.E.; Zecchini, H.I.; et al. Hyaluronan impairs vascular function and drug delivery in a mouse model of pancreatic cancer. Gut 2013, 62, 112–120. [Google Scholar] [CrossRef] [PubMed]
  5. Sherman, M.H.; Ruth, T.Y.; Engle, D.D.; Ding, N.; Atkins, A.R.; Tiriac, H.; Collisson, E.A.; Connor, F.; Van Dyke, T.; Kozlov, S.; et al. Vitamin D Receptor-Mediated Stromal Reprogramming Suppresses Pancreatitis and Enhances Pancreatic Cancer Therapy. Cell 2014, 159, 80–93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Lo, A.; Wang, L.C.S.; Scholler, J.; Monslow, J.; Avery, D.; Newick, K.; O’Brien, S.; Evans, R.A.; Bajor, D.J.; Clendenin, C.; et al. Tumor-Promoting Desmoplasia Is Disrupted by Depleting FAP-Expressing Stromal Cells. Cancer Res. 2015, 75, 2800–2810. [Google Scholar] [CrossRef] [Green Version]
  7. Beatty, G.L.; Chiorean, E.G.; Fishman, M.P.; Saboury, B.; Teitelbaum, U.R.; Sun, W.; Huhn, R.D.; Song, W.; Li, D.; Sharp, L.L.; et al. CD40 agonists alter tumor stroma and show efficacy against pancreatic carcinoma in mice and humans. Science 2011, 331, 1612–1616. [Google Scholar] [CrossRef] [Green Version]
  8. Cao, J.; Pickup, S.; Clendenin, C.; Blouw, B.; Choi, H.; Kang, D.; Rosen, M.; O’Dwyer, P.J.; Zhou, R. Dynamic Contrast-enhanced MRI Detects Responses to Stroma-directed Therapy in Mouse Models of Pancreatic Ductal Adenocarcinoma. Clin. Cancer Res. 2019, 25, 2314–2322. [Google Scholar] [CrossRef] [PubMed]
  9. White, N.S.; McDonald, C.R.; Farid, N.; Kuperman, J.; Karow, D.; Schenker-Ahmed, N.M.; Bartsch, H.; Rakow-Penner, R.; Holland, D.; Shabaik, A.; et al. Diffusion-weighted imaging in cancer: Physical foundations and applications of restriction spectrum imaging. Cancer Res. 2014, 74, 4638–4652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Le Bihan, D. Apparent diffusion coefficient and beyond: What diffusion MR imaging can tell us about tissue structure. Radiology 2013, 268, 318–322. [Google Scholar] [CrossRef]
  11. Taouli, B.; Beer, A.J.; Chenevert, T.; Collins, D.; Lehman, C.; Matos, C.; Padhani, A.R.; Rosenkrantz, A.B.; Shukla-Dave, A.; Sigmund, E.; et al. Diffusion-weighted imaging outside the brain: Consensus statement from an ISMRM-sponsored workshop. J. Magn. Reson. Imaging 2016, 44, 521–540. [Google Scholar] [CrossRef]
  12. Roth, Y.; Tichler, T.; Kostenich, G.; Ruiz-Cabello, J.; Maier, S.E.; Cohen, J.S.; Orenstein, A.; Mardor, Y. High-b-value diffusion-weighted MR imaging for pretreatment prediction and early monitoring of tumor response to therapy in mice. Radiology 2004, 232, 685–692. [Google Scholar] [CrossRef]
  13. Koh, D.M.; Scurr, E.; Collins, D.; Kanber, B.; Norman, A.; Leach, M.O.; Husband, J.E. Predicting response of colorectal hepatic metastasis: Value of pretreatment apparent diffusion coefficients. Am. J. Roentgenol. 2007, 188, 1001–1008. [Google Scholar] [CrossRef]
  14. Galbán, C.J.; Ma, B.; Malyarenko, D.; Pickles, M.D.; Heist, K.; Henry, N.L.; Schott, A.F.; Neal, C.H.; Hylton, N.M.; Rehemtulla, A.; et al. Multi-site clinical evaluation of DW-MRI as a treatment response metric for breast cancer patients undergoing neoadjuvant chemotherapy. PLoS ONE 2015, 10, e0122151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Dolde, K.; Dávid, C.; Echner, G.; Floca, R.; Hentschke, C.; Maier, F.; Niebuhr, N.; Ohmstedt, K.; Saito, N.; Alimusaj, M.; et al. 4DMRI-based analysis of inter—And intrafractional pancreas motion and deformation with different immobilization devices. Biomed. Phys. Eng. Express 2019, 5, 025012. [Google Scholar] [CrossRef]
  16. Jones, B.L.; Schefter, T.; Miften, M. Adaptive motion mapping in pancreatic SBRT patients using Fourier transforms. Radiol. Oncol. 2015, 115, 217–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Batista, V.; Richter, D.; Chaudhri, N.; Naumann, P.; Herfarth, K.; Jäkel, O. Significance of intra-fractional motion for pancreatic patients treated with charged particles. Radiat. Oncol. 2018, 13, 120. [Google Scholar] [CrossRef] [Green Version]
  18. Ferris, W.S.; Kissick, M.W.; Bayouth, J.E.; Culberson, W.S.; Smilowitz, J.B. Evaluation of radixact motion synchrony for 3D respiratory motion: Modeling accuracy and dosimetric fidelity. J. Appl. Clin. Med. Phys. 2020, 21, 96–106. [Google Scholar] [CrossRef]
  19. Bammer, R.; Keeling, S.L.; Augustin, M.; Pruessmann, K.P.; Wolf, R.; Stollberger, R.; Hartung, H.P.; Fazekas, F. Improved diffusion-weighted single-shot echo-planar imaging (EPI) in stroke using sensitivity encoding (SENSE). Magn. Reson. Med. 2001, 46, 548–554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Porter, D.A.; Heidemann, R.M. High resolution diffusion-weighted imaging using readout-segmented echo-planar imaging, parallel imaging and a two-dimensional navigator-based reacquisition. Magn. Reson. Med. 2009, 62, 468–475. [Google Scholar] [CrossRef]
  21. Gareis, D.; Wichmann, T.; Lanz, T.; Melkus, G.; Horn, M.; Jakob, P.M. Mouse MRI using phased-array coils. NMR Biomed. 2007, 20, 326–334. [Google Scholar] [CrossRef] [PubMed]
  22. Ribot, E.J.; Trotier, A.J.; Castets, C.R.; Dallaudière, B.; Thiaudière, E.; Franconi, J.M.; Miraux, S. Free-breathing 3D diffusion MRI for high-resolution hepatic metastasis characterization in small animals. Clin. Exp. Metastasis 2016, 33, 167–178. [Google Scholar] [CrossRef]
  23. Lee, H.J.; Luci, J.J.; Tantawy, M.N.; Lee, H.; Nam, K.T.; Peterson, T.E.; Price, R.R. Detecting peritoneal dissemination of ovarian cancer in mice by DWIBS. Magn. Reson. Imaging 2013, 31, 227–234. [Google Scholar] [CrossRef] [Green Version]
  24. Chow, A.M.; Gao, D.S.; Fan, S.J.; Qiao, Z.; Lee, F.Y.; Yang, J.; Man, K.; Wu, E.X. Liver fibrosis: An intravoxel incoherent motion (IVIM) study. J. Magn. Reson. Imaging 2012, 36, 159–167. [Google Scholar] [CrossRef] [PubMed]
  25. Lee, Y.; Kim, H. Assessment of diffusion tensor MR imaging (DTI) in liver fibrosis with minimal confounding effect of hepatic steatosis. Magn. Reson. Med. 2015, 73, 1602–1608. [Google Scholar] [CrossRef] [PubMed]
  26. Sun, J.; Zhang, X.P.; Li, X.T.; Tang, L.; Cui, Y.; Zhang, X.Y.; Sun, Y.S. Applicable apparent diffusion coefficient of an orthotopic mouse model of gastric cancer by improved clinical MRI diffusion weighted imaging. Sci. Rep. 2014, 4, 6072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Ge, X.; Quirk, J.D.; Engelbach, J.A.; Bretthorst, G.L.; Li, S.; Shoghi, K.I.; Garbow, J.R.; Ackerman, J.J. Test-Retest Performance of a 1-Hour Multiparametric MR Image Acquisition Pipeline with Orthotopic Triple-Negative Breast Cancer Patient-Derived Tumor Xenografts. Tomography 2019, 5, 320. [Google Scholar] [PubMed]
  28. Bao, Q.; Liberman, G.; Solomon, E.; Frydman, L. High-resolution diffusion MRI studies of development in pregnant mice visualized by novel spatiotemporal encoding schemes. NMR Biomed. 2020, 33, e4208. [Google Scholar] [CrossRef] [PubMed]
  29. Ma, L.; Cai, C.; Yang, H.; Cai, S.; Qian, J.; Xiao, L.; Zhong, K.; Chen, Z. Motion-tolerant diffusion mapping based on single-shot overlapping-echo detachment (OLED) planar imaging. Magn. Reson. Med. 2018, 80, 200–210. [Google Scholar] [CrossRef]
  30. Kennedy, S.D.; Zhong, J. Diffusion measurements free of motion artifacts using intermolecular dipole-dipole interactions. Magn. Reson. Med. 2004, 52, 1–6. [Google Scholar] [CrossRef]
  31. Zaitsev, M.; Maclaren, J.; Herbst, M. Motion artifacts in MRI: A complex problem with many partial solutions. J. Magn. Reson. Imaging 2015, 42, 887–901. [Google Scholar] [CrossRef]
  32. Lin, W.; Guo, J.; Rosen, M.A.; Song, H.K. Respiratory motion-compensated radial dynamic contrast-enhanced (DCE)-MRI of chest and abdominal lesions. Magn. Reson. Med. 2008, 60, 1135–1146. [Google Scholar] [CrossRef] [Green Version]
  33. Dougherty, L.; Isaac, G.; Rosen, M.A.; Nunes, L.W.; Moate, P.J.; Boston, R.C.; Schnall, M.D.; Song, H.K. High frame-rate simultaneous bilateral breast DCE-MRI. Magn. Reson. Med. 2007, 57, 220–225. [Google Scholar] [CrossRef] [PubMed]
  34. Rosenkrantz, A.B.; Khasgiwala, A.; Doshi, A.M.; Ream, J.M.; Taneja, S.S.; Lepor, H. Detection of prostate cancer local recurrence following radical prostatectomy: Assessment using a continuously acquired radial golden-angle compressed sensing acquisition. Abdom. Radiol. (N. Y.) 2017, 42, 290–297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Heacock, L.; Gao, Y.; Heller, S.L.; Melsaether, A.N.; Babb, J.S.; Block, T.K.; Otazo, R.; Kim, S.G.; Moy, L. Comparison of conventional DCE-MRI and a novel golden-angle radial multicoil compressed sensing method for the evaluation of breast lesion conspicuity. J. Magn. Reson. Imaging 2017, 45, 1746–1752. [Google Scholar] [CrossRef] [PubMed]
  36. Shukla-Dave, A.; Obuchowski, N.A.; Chenevert, T.L.; Jambawalikar, S.; Schwartz, L.H.; Malyarenko, D.; Huang, W.; Noworolski, S.M.; Young, R.J.; Shiroishi, M.S.; et al. Quantitative imaging biomarkers alliance (QIBA) recommendations for improved precision of DWI and DCE-MRI derived biomarkers in multicenter oncology trials. J. Magn. Reson. Imaging 2019, 49, e101–e121. [Google Scholar] [CrossRef] [PubMed]
  37. Kjølby, B.F.; Khan, A.R.; Chuhutin, A.; Pedersen, L.; Jensen, J.B.; Jakobsen, S.; Zeidler, D.; Sangill, R.; Nyengaard, J.R.; Jespersen, S.N.; et al. Fast diffusion kurtosis imaging of fibrotic mouse kidneys. NMR Biomed. 2016, 29, 1709–1719. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Garrett-Roe, S.; Warren, W.S. Numerical studies of intermolecular multiple quantum coherences: High-resolution NMR in inhomogeneous fields and contrast enhancement in MRI. J. Magn. Reson. 2000, 146, 1–13. [Google Scholar] [CrossRef]
  39. Neeman, M.; Freyer, J.P.; Sillerud, L.O. A simple method for obtaining cross-term-free images for diffusion anisotropy studies in NMR microimaging. Magn. Reson. Med. 1991, 21, 138–143. [Google Scholar] [CrossRef]
  40. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef] [PubMed]
  41. Feinberg, D.A.; Oshio, K. Phase errors in multi-shot echo planar imaging. Magn. Reson. Med. 1994, 32, 535–539. [Google Scholar] [CrossRef]
  42. Tofts, P.S.; Lloyd, D.; Clark, C.A.; Barker, G.J.; Parker, G.J.M.; McConville, P.; Baldock, C.; Pope, J.M. Test liquids for quantitative MRI measurements of self-diffusion coefficient in vivo. Magn. Reson. Med. 2000, 43, 368–374. [Google Scholar] [CrossRef] [Green Version]
  43. Patlak, C.S.; Hospod, F.E.; Trowbridge, S.D.; Newman, G.C. Diffusion of radiotracers in normal and ischemic brain slices. J. Cereb. Blood Flow Metab. 1998, 18, 776–802. [Google Scholar] [CrossRef] [Green Version]
  44. Sharpless, N.E.; Depinho, R.A. The mighty mouse: Genetically engineered mouse models in cancer drug development. Nat. Rev. Drug Discov. 2006, 5, 741–754. [Google Scholar] [CrossRef]
  45. Rosenzweig, S.; Holme, H.C.M.; Uecker, M. Simple auto-calibrated gradient delay estimation from few spokes using Radial Intersections (RING). Magn. Reson. Med. 2019, 81, 1898–1906. [Google Scholar] [CrossRef] [PubMed]
  46. Block, K.T.; Uecker, M.; Frahm, J. Undersampled radial MRI with multiple coils. Iterative image reconstruction using a total variation constraint. Magn. Reson. Med. 2007, 57, 1086–1098. [Google Scholar] [CrossRef] [PubMed]
  47. Wagner, F.; Laun, F.B.; Kuder, T.A.; Mlynarska, A.; Maier, F.; Faust, J.; Demberg, K.; Lindemann, L.; Rivkin, B.; Nagel, A.M. Temperature and concentration calibration of aqueous polyvinylpyrrolidone (PVP) solutions for isotropic diffusion MRI phantoms. PLoS ONE 2017, 12, e0179276. [Google Scholar] [CrossRef] [Green Version]
  48. Niwa, T.; Ueno, M.; Ohkawa, S.; Yoshida, T.; Doiuchi, T.; Ito, K.; Inoue, T. Advanced pancreatic cancer: The use of the apparent diffusion coefficient to predict response to chemotherapy. Br. J. Radiol. 2009, 82, 28–34. [Google Scholar] [CrossRef]
  49. 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, 1302–1308. [Google Scholar] [CrossRef] [PubMed]
  50. Whisenant, J.G.; Ayers, G.D.; Loveless, M.E.; Barnes, S.L.; Colvin, D.C.; Yankeelov, T.E. Assessing reproducibility of diffusion-weighted magnetic resonance imaging studies in a murine model of HER2+ breast cancer. Magn. Reson. Imaging 2014, 32, 245–249. [Google Scholar] [CrossRef] [Green Version]
  51. Shoghi, K.I.; Badea, C.T.; Blocker, S.J.; Chenevert, T.L.; Laforest, R.; Lewis, M.T.; Luker, G.D.; Manning, H.C.; Marcus, D.S.; Mowery, Y.M.; et al. Co-Clinical Imaging Resource Program (CIRP): Bridging the Translational Divide to Advance Precision Medicine. Tomography 2020, 6, 273–287. [Google Scholar] [PubMed]
  52. Barnes, A.; Alonzi, R.; Blackledge, M.; Charles-Edwards, G.; Collins, D.J.; Cook, G.; Coutts, G.; Goh, V.; Graves, M.; Kelly, C.; et al. UK quantitative WB-DWI technical workgroup: Consensus meeting recommendations on optimisation, quality control, processing and analysis of quantitative whole-body diffusion-weighted imaging for cancer. Br. J. Radiol. 2018, 91, 20170577. [Google Scholar] [CrossRef] [PubMed]
  53. Trouard, T.P.; Theilmann, R.J.; Altbach, M.I.; Gmitro, A.F. High-resolution diffusion imaging with DIFRAD-FSE (diffusion-weighted radial acquisition of data with fast spin-echo) MRI. Magn. Reson. Med. 1999, 42, 11–18. [Google Scholar] [CrossRef]
  54. Wang, F.N.; Huang, T.Y.; Lin, F.H.; Chuang, T.C.; Chen, N.K.; Chung, H.W.; Chen, C.Y.; Kwong, K.K. PROPELLER EPI: An MRI technique suitable for diffusion tensor imaging at high field strength with reduced geometric distortions. Magn. Reson. Med. 2005, 54, 1232–1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Knoll, F.; Raya, J.G.; Halloran, R.O.; Baete, S.; Sigmund, E.; Bammer, R.; Block, T.; Otazo, R.; Sodickson, D.K. A model-based reconstruction for undersampled radial spin-echo DTI with variational penalties on the diffusion tensor. NMR Biomed. 2015, 28, 353–366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Diffusion-weighted images and the resulting apparent diffusion coefficient (ADC) maps from the radially sampled diffusion-weighted spin-echo (DW-SE-RAD) protocol with and without prospective respiratory gating. The ADC estimates of the tumor (mean ± SD from multiple slices spanning the tumor in units of × 10−3 mm2/s) are 1.1 ± 0.097 vs. 1.3 ± 0.17 (gated vs. non-gated), while estimates of the diffusion coefficient of the water phantom (in units of × 10−3 mm2/s) are 3.2 ± 0.13 vs. 4.0 ± 0.30 (gated vs. non-gated).
Figure 1. Diffusion-weighted images and the resulting apparent diffusion coefficient (ADC) maps from the radially sampled diffusion-weighted spin-echo (DW-SE-RAD) protocol with and without prospective respiratory gating. The ADC estimates of the tumor (mean ± SD from multiple slices spanning the tumor in units of × 10−3 mm2/s) are 1.1 ± 0.097 vs. 1.3 ± 0.17 (gated vs. non-gated), while estimates of the diffusion coefficient of the water phantom (in units of × 10−3 mm2/s) are 3.2 ± 0.13 vs. 4.0 ± 0.30 (gated vs. non-gated).
Tomography 07 00007 g001
Figure 2. Representative diffusion-weighted images and ADC maps obtained from the two diffusion-weighted MRI (DWI) protocols. Rows A and B are diffusion-weighted images at similar slice positions generated by the DW-SE-RAD and diffusion-weighted spin-echo 4-shot echo-planar imaging (DW-SE-EPI) protocol, respectively. For each row, the first five images (1–5) are windowed to highlight the decay due to diffusion weighting, while frames (A5’,B5’) are windowed to optimally display (A5,B5), respectively. (A6,B6) are ADC maps (mm2/s) displayed in color overlay for the tumor and in grayscale for other tissues and phantoms. Blue arrows point to the water phantom, while yellow arrows to the 1-butanol phantom. Due to its low diffusion coefficient (~0.44 × 10−3 mm2/s), 1-butanol appears dark in the gray-scaled ADC map (A6) produced from the DW-SE-RAD protocol.
Figure 2. Representative diffusion-weighted images and ADC maps obtained from the two diffusion-weighted MRI (DWI) protocols. Rows A and B are diffusion-weighted images at similar slice positions generated by the DW-SE-RAD and diffusion-weighted spin-echo 4-shot echo-planar imaging (DW-SE-EPI) protocol, respectively. For each row, the first five images (1–5) are windowed to highlight the decay due to diffusion weighting, while frames (A5’,B5’) are windowed to optimally display (A5,B5), respectively. (A6,B6) are ADC maps (mm2/s) displayed in color overlay for the tumor and in grayscale for other tissues and phantoms. Blue arrows point to the water phantom, while yellow arrows to the 1-butanol phantom. Due to its low diffusion coefficient (~0.44 × 10−3 mm2/s), 1-butanol appears dark in the gray-scaled ADC map (A6) produced from the DW-SE-RAD protocol.
Tomography 07 00007 g002
Figure 3. Correlation analysis of the tumor ADC values estimated by the DW-SE-RAD vs. DW-SE-EPI protocols using the generalized estimating equations approach. The analysis was applied to mice (n = 8) that were subjected to both DWI protocols in both the test and retest scans (detailed in Methods).
Figure 3. Correlation analysis of the tumor ADC values estimated by the DW-SE-RAD vs. DW-SE-EPI protocols using the generalized estimating equations approach. The analysis was applied to mice (n = 8) that were subjected to both DWI protocols in both the test and retest scans (detailed in Methods).
Tomography 07 00007 g003
Figure 4. Bland-Altman plots of the tumor (A,C) and muscle (B,D) ADC measurements obtained from the test–retest experiment for each DWI protocol. The mean of ΔADC and mean ± RC are specified in the figure, while the RC values are reported in Table 1. Red arrows point to the worst cases in each DWI protocol. ΔADC = ADCtest − ADCretest, RC is repeatability coefficient defined in Equation [4].
Figure 4. Bland-Altman plots of the tumor (A,C) and muscle (B,D) ADC measurements obtained from the test–retest experiment for each DWI protocol. The mean of ΔADC and mean ± RC are specified in the figure, while the RC values are reported in Table 1. Red arrows point to the worst cases in each DWI protocol. ΔADC = ADCtest − ADCretest, RC is repeatability coefficient defined in Equation [4].
Tomography 07 00007 g004
Figure 5. Probability density function (PDF) of the pixel-wise tumor ADC values from the best (A,B) and worst cases (C,D) in the test–retest study for the DW-SE-RAD and DW-SE-EPI protocols, respectively (solid line: test and dotted line: retest). The best and worst cases correspond to the smallest and largest ΔD values, respectively, from ten mice subjected to the test–retest study for each protocol.
Figure 5. Probability density function (PDF) of the pixel-wise tumor ADC values from the best (A,B) and worst cases (C,D) in the test–retest study for the DW-SE-RAD and DW-SE-EPI protocols, respectively (solid line: test and dotted line: retest). The best and worst cases correspond to the smallest and largest ΔD values, respectively, from ten mice subjected to the test–retest study for each protocol.
Tomography 07 00007 g005
Table 1. Statistical parameters derived from the apparent diffusion coefficient (ADC) estimates in the test–retest study a.
Table 1. Statistical parameters derived from the apparent diffusion coefficient (ADC) estimates in the test–retest study a.
ProtocolADC
(Mean ± SD)
ΔD
(Mean ± SD)
SDwsCVWSRC
Test bRetest b
Water cDW-SE-RAD3.2 ± 0.293.3 ± 0.27−0.048 ± 0.280.190.0600.53
DW-SE-EPI2.8 ± 0.152.8 ± 0.100.069 ± 0.150.110.0390.31
MuscleDW-SE-RAD1.8 ± 0.291.7 ± 0.250.045 ± 0.200.140.0730.38
DW-SE-EPI1.8 ± 0.221.7 ± 0.250.060 ± 0.260.180.0960.50
TumorDW-SE-RAD1.3 ± 0.191.3 ± 0.29−0.017 ± 0.180.120.0900.34
DW-SE-EPI1.5 ± 0.321.5 ± 0.44−0.082 ± 0.340.240.130.66
a: All values are reported in the unit of × 10−3 mm2/s, except for CVWS and RC, which are unitless. ΔD, SDws, CVWS, and RC are defined in Equations (2)–(4). b: The differences of the ADC between test and retest are not significant (p > 0.05). c: The difference of ADC derived from radially sampled diffusion-weighted spin-echo (DW-SE-RAD) vs. diffusion-weighted spin-echo 4-shot echo-planar imaging (DW-SE-EPI) is significant for the water phantom (p < 0.001) but not significant for the muscle or tumor (p > 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cao, J.; Song, H.K.; Yang, H.; Castillo, V.; Chen, J.; Clendenin, C.; Rosen, M.; Zhou, R.; Pickup, S. Respiratory Motion Mitigation and Repeatability of Two Diffusion-Weighted MRI Methods Applied to a Murine Model of Spontaneous Pancreatic Cancer. Tomography 2021, 7, 66-79. https://0-doi-org.brum.beds.ac.uk/10.3390/tomography7010007

AMA Style

Cao J, Song HK, Yang H, Castillo V, Chen J, Clendenin C, Rosen M, Zhou R, Pickup S. Respiratory Motion Mitigation and Repeatability of Two Diffusion-Weighted MRI Methods Applied to a Murine Model of Spontaneous Pancreatic Cancer. Tomography. 2021; 7(1):66-79. https://0-doi-org.brum.beds.ac.uk/10.3390/tomography7010007

Chicago/Turabian Style

Cao, Jianbo, Hee Kwon Song, Hanwen Yang, Victor Castillo, Jinbo Chen, Cynthia Clendenin, Mark Rosen, Rong Zhou, and Stephen Pickup. 2021. "Respiratory Motion Mitigation and Repeatability of Two Diffusion-Weighted MRI Methods Applied to a Murine Model of Spontaneous Pancreatic Cancer" Tomography 7, no. 1: 66-79. https://0-doi-org.brum.beds.ac.uk/10.3390/tomography7010007

Article Metrics

Back to TopTop