Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Evaluation of different mathematical models and different b-value ranges of diffusion-weighted imaging in peripheral zone prostate cancer detection using b-value up to 4500 s/mm2

  • Zhaoyan Feng,

    Affiliation Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

  • Xiangde Min,

    Affiliation Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

  • Daniel J. A. Margolis,

    Affiliation Department of Radiology, David Geffen School of Medicine at UCLA, Ronald Reagan UCLA Medical Center, Los Angeles, California, United States of America

  • Caohui Duan,

    Affiliation Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, National Center for Magnetic Resonance in Wuhan, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan, Hubei, China

  • Yuping Chen,

    Affiliation Centers for Biomedical Engineering, University of Science and Technology of China, Hefei, Anhui, China

  • Vivek Kumar Sah,

    Affiliation Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

  • Nabin Chaudhary,

    Affiliation Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

  • Basen Li,

    Affiliation Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

  • Zan Ke,

    Affiliation Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

  • Peipei Zhang,

    Affiliation Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

  • Liang Wang

    wang6@tjh.tjmu.edu.cn

    Affiliation Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

Correction

6 Sep 2017: The PLOS ONE Staff (2017) Correction: Evaluation of different mathematical models and different b-value ranges of diffusion-weighted imaging in peripheral zone prostate cancer detection using b-value up to 4500 s/mm2. PLOS ONE 12(9): e0184609. https://doi.org/10.1371/journal.pone.0184609 View correction

Abstract

Objectives

To evaluate the diagnostic performance of different mathematical models and different b-value ranges of diffusion-weighted imaging (DWI) in peripheral zone prostate cancer (PZ PCa) detection.

Methods

Fifty-six patients with histologically proven PZ PCa who underwent DWI-magnetic resonance imaging (MRI) using 21 b-values (0–4500 s/mm2) were included. The mean signal intensities of the regions of interest (ROIs) placed in benign PZs and cancerous tissues on DWI images were fitted using mono-exponential, bi-exponential, stretched-exponential, and kurtosis models. The b-values were divided into four ranges: 0–1000, 0–2000, 0–3200, and 0–4500 s/mm2, grouped as A, B, C, and D, respectively. ADC, <D>, D*, f, DDC, α, Dapp, and Kapp were estimated for each group. The adjusted coefficient of determination (R2) was calculated to measure goodness-of-fit. Receiver operating characteristic curve analysis was performed to evaluate the diagnostic performance of the parameters.

Results

All parameters except D* showed significant differences between cancerous tissues and benign PZs in each group. The area under the curve values (AUCs) of ADC were comparable in groups C and D (p = 0.980) and were significantly higher than those in groups A and B (p< 0.05 for all). The AUCs of ADC and Kapp in groups B and C were similar (p = 0.07 and p = 0.954), and were significantly higher than the other parameters (p< 0.001 for all). The AUCs of ADC in group D was slightly higher than Kapp (p = 0.002), and both were significantly higher than the other parameters (p< 0.001 for all).

Conclusions

ADC derived from conventional mono-exponential high b-value (3200 s/mm2) models is an optimal parameter for PZ PCa detection.

Introduction

Magnetic resonance imaging (MRI) is considered the best imaging technique for the detection, staging and surveillance of prostate cancer (PCa) because of its ability to provide functional and anatomic information [13]. Diffusion-weighted imaging (DWI) is a powerful and dominant component of MRI in peripheral zone (PZ) PCa detection[4].

In traditional diffusion theory, a mono-exponential model, based on the assumption of homogeneous Gaussian diffusion, has been adopted for diffusion analysis in most clinical studies [5, 6]. However, molecular diffusion in biological tissues does not occur freely, but rather faces many forms of hindrance, including that from membranes and intracellular organelles[7]. Therefore, many non-Gaussian diffusion models, including bi-exponential, stretched-exponential, and kurtosis models, have been developed to describe the complicated behavior of water diffusion. Some previous studies have suggested that the bi-exponential model [8, 9], stretched-exponential model [10], and diffusion kurtosis imaging [11] might provide useful information for PCa evaluation. However, the observed benefits from those non-Gaussian diffusion models required validation [12, 13], warranting further studies to explore and compare their roles in the diagnosis of PCa.

The choice of an appropriate b-value is another crucial aspect of DWI that affects PCa detection efficiency. The appropriate b-value is essential for calculating precise diffusion parameters to aid PCa diagnosis [14, 15]. Nevertheless, the maximum b-values used in prostate studies have varied greatly, from 500 to 3000 s/mm2[12, 16, 17]. while the brain delivered values up to 5000 s/mm2[18]. An earlier study investigated the influence of b-value range on DWI parameters and its effect on diagnostic performance for PCa [16]. However, the maximum b-value in that study was 2300 s/mm2; the effect of an ultrahigh b-value (up to 4500 s/mm2) remains unclear.

In this study, we aimed to address the influence of the range of b-values on the discriminatory value of the diffusion parameters of four different mathematical models (mono-exponential, bi-exponential, stretched-exponential and kurtosis) and to determine the optimal parameter and its optimal b-value range for PCa detection.

Materials and methods

Patients

This prospective study was approved by the Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology Institutional Review Board; written informed consent was obtained from each patient prior to examination. Between May 2014 and March 2015, patients with clinically suspected PCa due to either a rising or elevated prostate-specific antigen (PSA) level and/or palpable prostatic nodule, with no prior therapy, were recruited to this study. Initially, 166 consecutive patients were identified. Patients were excluded if a) they had histological confirmation of no cancer (n = 43); b) they had no histological confirmation of cancer (n = 51); c) the images were of poor quality due to susceptibility artifact or movement (n = 4); or d) the cancerous focus was located in the transition zone (n = 12). After these exclusions, 56 patients remained.

MRI protocol

All examinations were performed on a 3.0 T GE 750 discovery system MRI scanner (General Electric, Milwaukee, WI, USA) with a 32–element torso phased-array coil. An endorectal coil was not used. The entire prostate gland and seminal vesicles were imaged in axial, sagittal and coronal slices, using a T2-weighted fast-recovery fast spin echo sequence with the following parameters TR range/TE range, 3322-3849/108-118 ms; slice thickness, 3 mm; inter-slice gap, 0 mm; FOV, 180×180 mm2; and matrix, 320×256. Transversal T1-weighted images were acquired with the following parameters: TR, 446 ms; TE, 8.1 ms; slice thickness, 5 mm; inter-slice gap, 0 mm; FOV, 360×360 mm2; and matrix, 320×224.

Axial DWI was acquired using a single-shot echo-planar imaging pulse sequence with the following parameters: TR, 2500 ms; TE, 84.1 ms; slice thickness, 5 mm; inter-slice gap, 0 mm; FOV, 400×280 mm2; matrix, 128×96, reconstruction matrix size 256 × 256; bandwidth, 250 Hz/pixel; acceleration factor of 2; and 21 b values (number of excitations) 0 (1), 20 (1), 50 (1), 80 (1), 100 (1), 150 (1), 200 (1), 400 (2), 600 (2), 800 (2), 1000 (4), 1200 (4), 1500 (4), 1800 (4), 2000 (6), 2400 (6), 2800 (6), 3200 (8), 3600 (8), 4000 (10), and 4500 (10) s/mm2. The DWI parameter was optimized to ensure a sufficient SNR at the maximum b-values up to 4500mm/s2. The total acquisition time was 10 minutes and 20 seconds. An additional 5mm section thickness axial T2-weighted image was scanned to match the DWI to provide anatomical information.

TRUS-guided biopsy

After the MRI examination, all patients underwent transrectal ultrasound (TRUS)-guided biopsies. All the biopsies were performed by a single urologist with more than 20 years of experience using an ultrasound system (Hawk 2102, BK Medical, Denmark) equipped with a 5.1 MHz endocavitary probe, with a spring-loaded biopsy gun with an 18 gauge core biopsy needle. To match biopsy sextants and MR images, the PZ of the prostate was divided into six regions in MRI, which corresponded to the TRUS-guided biopsy (apex, mid gland, and base on each side of the PZ), and two core biopsies per region were obtained. Cores were individually labeled according to their location with respect to the biopsy scheme. Each specimen was histologically analyzed as cancerous or noncancerous by an experienced pathologist (with > 15 years of experience).

MRI data assessment

MRI assessment was performed according to the methods used by Rosenkrantz et al. and Li et al.[11, 19], and post-processing was performed using the FUNCTOOL (GE Healthcare 4.6) workstation. MR images were interpreted by the consensus of two observers (ZY Feng, and XD Min) with more than 5 years of experience in prostate MRI (each observer had read more than one thousand prostate MRIs). Disagreements between the two observers were resolved by a re-evaluation with the study coordinator, a radiologist with more than 15 years of experience in prostate MRI. On the MR images, the prostate was divided into approximately equal thirds in the cranio-caudal direction, with each third demarcating the left and right basal, mid, and apical sextants, corresponding to the TRUS-guided biopsy. The readers treated each sextant as a separate focus for purposes of the analysis. Regions of interest (ROIs) were placed within the proven cancerous tissues and benign PZs on the T2WI and then automatically copied to DWI. The focus suspicious for PCa with a diameter of not less than 5 mm was delineated on the MRI of the prostate corresponding to the biopsy sextants. If both biopsies obtained in one region were negative or positive for PCa, the region was considered cancer negative or cancer positive respectively. If one biopsy was positive for cancer but the other was negative, the region was still considered cancer positive [11]. The sizes of the ROIs were chosen to be as large as possible with minimal contamination from uninvolved tissues. The SNR for each ROI was calculated for DWI with b-values of 4500 s/mm2. The SNR was defined as the ratio between the average signal intensity of ROIs and the mean of the standard deviation (SD) of the signals in the bilateral internal obturator muscles[20].

Modeling

All analyses in the current study were performed on an ROI level [13, 21]. We performed the non-linear least-square fitting based on the Levenberg-Marquardt algorithm using Matlab (MathWorks, Natick, MA, USA) to fit the following models to the mean signal of each ROI (Fig 1):

  1. Mono-exponential model (1)
  2. Bi-exponential model (2)
  3. Stretched exponential model (3)
  4. Diffusion kurtosis imaging (4)

Where S(b) is the signal intensity at a particular b-value, and S0 is the signal intensity at b = 0 s/mm2. ADC is the diffusion coefficient of the mono-exponential model. <D>, D*, and f are the diffusion parameters of the bi-exponential model: f is the perfusion fraction, <D> is the pure molecular diffusion, and D* is the pseudo-diffusion coefficient. DDC is the diffusion coefficient of the stretched exponential model, and α is the water diffusion heterogeneity index between 0 and 1. Dapp is the diffusion coefficient of the kurtosis model, and Kapp is the kurtosis.

thumbnail
Fig 1. Prostate cancer in a 61-year-old patient with high serum PSA level of 19.43 ng/ml.

(a)Transverse T2-weighted anatomical image shows a hypointense lesion (white arrow) in the right middle peripheral zone of the prostate; (b) A series of b-value images are shown with the corresponding location to the transverse T2-weighted (unites, s/mm2); (c) Measured signal and fitted curves of cancerous tissue and the opposite side benign tissue using maximum b-value of 4500 s/mm2 (group D).

https://doi.org/10.1371/journal.pone.0172127.g001

Statistical analysis

To validate the dependence of the models on b-values, the b-values were divided into four different ranges: 0–1000, 0–2000, 0–3200, and 0–4500 s/mm2, grouped as A, B, C, and D, respectively.

To measure the goodness-of-fit, adjusted coefficients of determination (R2) for the groups were calculated among different models using Eqs (1)–(4). One-way ANOVA statistics and Bonferroni test were used for multiple comparisons of adjusted R2 between models.

ADC, <D>, D*, f, DDC, α, Dapp, and Kapp were estimated in the benign PZs and cancerous tissues. The parameters were calculated using all the b-values in groups A, B, C, and D. Normal distribution was assessed using the Kolmogorov-Smirnov test; if the data presented a normal distribution, an independent t-test was performed to compare the differences in parameters between benign PZs and cancerous tissues. For non-normally distributed variables, the Mann-Whitney U test was used. Next, receiver operating characteristic curve (ROC) analysis was performed, and the area under the curve (AUC) was calculated to evaluate and compare the diagnostic accuracy of the parameters in distinguishing cancerous tissues from benign PZs. The diagnostic performance was expressed as the AUC. Differences in AUCs were assessed using the method of DeLong et al[22].

The adjusted R2 was calculated with Matlab (MathWorks, Natick, MA, USA). The statistical analysis was performed using SPSS software (SPSS for Windows 19.0, Chicago, IL, USA). The ROC analysis was performed using MedCalc version 13.0.0.0 for Windows (MedCalc Software, Mariakerke, Belgium). For all the significance level was set at p< 0.05.

Results

The final study population comprised 56 PZ PCa patients (age range, 52–82 years; median age, 67 years). The median PSA level in these 56 patients was 26.036 ng/mL (range, 2.06–1000 ng/ml). Among the total 336 ROIs, 198 ROIs were evaluated as benign PZs, and 138 ROIs were evaluated as cancerous tissues. The SNR of the ROIs at the highest b-value of 4500 s/mm2was between 8.784 and 64.118. The mean (±SD) ROI SNRs in cancerous tissues and benign PZs were 27.784±10.245 and 17.035±8.107, respectively. The mean (±SD) ROI sizes in cancerous tissues and benign PZs were 169.780±107.914mm2 and 97.57±22.041mm2, respectively.

Goodness-of-fit of the models

The adjusted R2 was calculated to test the goodness-of-fit of the models. The mean values and standard deviations of the adjusted R2 are provided in Table 1. The bi-exponential and the stretched-exponential models provided the highest adjusted R2 among the four models in every group. There were no significant differences between the bi-exponential and stretched-exponential models in any of the four groups (p>0.05 for all). The goodness-of-fit of the kurtosis model was slightly less than that of the bi-exponential and the stretched-exponential models in all groups (p< 0.05 for all). The mono-exponential provided the worst goodness-of-fit in the four groups. Significant differences were observed between kurtosis and the other three models in all groups (p< 0.05 for all). Similarly, the mono-exponential model showed differences from the other three models in all groups (p< 0.05 for all). The mean value of the adjusted R2 of the mono-exponential model decreased with increasing b-value. No clear tendency could be established for the bi-exponential, stretched-exponential or kurtosis models (Fig 2).

thumbnail
Fig 2. Graph showing variations in the adjusted R2 values of the mono-exponential, bi-exponential, stretched-exponential and kurtosis models with an increase in b-value.

The mean value of the adjusted R2 of the mono-exponential model decreased with the increase of b-value, while the mean value of the adjusted R2 of the bi-exponential, stretched-exponential and kurtosis models were stable and excellent with the increase of b-value (It should be noted that, because the adjusted R2 values of the bi-exponential and stretched-exponential models were close, their curves have been superimposed).

https://doi.org/10.1371/journal.pone.0172127.g002

thumbnail
Table 1. Mean values and standard deviation of the adjusted R2.

https://doi.org/10.1371/journal.pone.0172127.t001

Summary of DWI parameters

The mean values and standard deviations of the diffusion parameters measured in groups A, B, C, and D are summarized in Table 2. The mean values of ADC, <D>, f, DDC, α, and Dapp were found to be significantly lower in cancerous tissues than in benign PZs in each group (p< 0.001 for all). Similarly, the Kapp in cancerous tissues was significantly higher than that in benign PZs in each group (p< 0.001 for all). We also found a significant difference in D* between cancerous tissues and benign PZs in groups B and C (p< 0.05), but not in groups A and D (p > 0.05). However, the D* showed a large standard deviation in all groups(Fig 3).

thumbnail
Fig 3.

Boxplot of diffusion parameters calculated using different b-value ranges (groups A, B, C and D). (a-b, d-g), the mean values of ADC, <D>, f, DDC, α, and Dapp were significantly lower in cancerous tissues than in benign PZs in each group (h). The Kapp in cancerous tissues was significantly higher than in benign PZs in each group. The value of D*, which had a large standard deviation, showed no significant difference between cancerous tissues and benign PZs in groups A and D but was significantly different in groups B and C(c).

https://doi.org/10.1371/journal.pone.0172127.g003

thumbnail
Table 2. Diffusion parameters in benign PZs and cancerous tissues.

https://doi.org/10.1371/journal.pone.0172127.t002

Analysis of ROC curves

The AUC, cutoff, sensitivity, and specificity values of the parameters for distinguishing cancerous tissues from benign PZs are reported in Table 3. The ROC curves of each parameter calculated using different b-value ranges are shown in Fig 4.

thumbnail
Fig 4. ROC curve analyses show the diagnostic accuracy of the diffusion parameters in distinguishing between cancerous tissues and benign PZs.

In group A, Kapp had the largest AUC (0.940), but the AUCs of ADC and Kapp were not significantly different (p = 0.070). In groups B and C, the AUCs of ADC and Kapp were comparable and significantly higher than those of the other parameters. In group D, the AUCs of ADC was slightly higher than that of Kapp (0.957 vs 0.953, p = 0.002), and both were significantly higher than those of other parameters (p< 0.001 for all).

https://doi.org/10.1371/journal.pone.0172127.g004

thumbnail
Table 3. Diagnostic performance of parameters calculated with different b-value ranges.

https://doi.org/10.1371/journal.pone.0172127.t003

The parameters ADC, <D>, DDC and Kapp provided the highest AUCs among the four models. The AUCs of ADC in groups C and D were comparable (p = 0.980) and were significantly higher than those of groups A and B (p< 0.05 for all). The AUCs of <D> calculated in groups A, B, and C were significantly higher than in group D (p< 0.001 for all). The AUC of <D> in groups A and B showed no significant difference (p = 0.136), and the value of group B was significantly higher than in groups C and D (p< 0.001 for all). The AUCs of DDC calculated using groups B, C and D showed no significant difference (p > 0.05 for all). The AUCs of DDC calculated using groups B and D were significantly higher than in group A (p = 0.016 and p = 0.014, respectively). The AUCs of DDC calculated using groups A and C showed no significant difference (p = 0.102). The AUCs of Kapp calculated using groups B, C, D were significantly higher than in group A (p< 0.05 for all); however, no significant difference was observed among groups B, C, and D (p > 0.05 for all).

For the parameters calculated using the b-value of group A, the AUCs of DDC and Kapp were not significantly different (p = 0.204); the AUCs of ADC and Kapp were not significantly different (p = 0.070); the AUC of DDC was slightly higher than the AUC of ADC (p = 0.014); and the AUCs of ADC, DDC and Kapp were significantly higher than those of the other parameters (p< 0.05 for all). For the parameters calculated using the b-values of groups B and C, the AUCs of ADC and Kapp were similar (p = 0.07 and p = 0.954, respectively), but were significantly higher than those of the other parameters (p< 0.001 for all). For the parameters calculated using the b-value of group D, the AUC of ADC was slightly higher than that of Kapp (0.957 vs 0.953, p = 0.002); the AUCs of ADC and Kapp were significantly higher than those of the other parameters (p< 0.001 for all).

Discussion

Our study demonstrated that for PZ PCa detection, the parameters derived from the bi-exponential, stretched-exponential, and kurtosis models did not reveal significantly superior diagnostic performance relative to conventional ADC. ADC derived from conventional mono-exponential high b-value (3200 s/mm2) provided excellent diagnostic performance in PZ PCa detection, such that no extra benefit was obtained from cumbersome non-Gaussian models (bi-exponential, stretched-exponential, and kurtosis models) and time-consuming over-high b-value (> 3200 s/mm2).

Appropriate b-values must be carefully selected for the use of DWI in PCa detection. In clinical applications, DWI is typical performed using b-values 800–1500 s/mm2[4, 23]. Nevertheless, the clinical benefits of high-b-value DWI have been explored from head to body[17, 24, 25]. Until now, consensus on the optimal b-value range of DWI for PCa diagnosis has been lacking. As reported in a previous study, b-value distribution mainly influences the repeatability of DWI-derived parameters rather than the diagnostic performance [26]; however, the influence of b-value distribution on diagnostic performance is beyond the scope of this study. Our study focused on the influence of b-value range on the diagnostic performance of the parameters of the various DWI models.

Our study showed that, for PCa detection, the diagnostic performance of DDC and Kapp calculated using b-values up to 2000 s/mm2 outperformed those calculated using b-values up to a maximum of 1000 s/mm2. According to a previous study, DWI with a maximum b-value of 2000 s/mm2 is beneficial for the estimation of the stretched-exponential and kurtosis model parameters[27]. Our results are also in agreement with a study by Mazzoni[16], which found that the Kapp demonstrated higher AUC for PCa detection using a higher b-value (2300 s/mm2) than a low b-value (800 s/mm2). However, the highest b-value they used was only 2300 s/mm2. For brain kurtosis imaging, empirical evidence indicates that maximum b-values of 2000 to 3000 s/mm2 are appropriate[7, 28]. However, no reasonable maximum b-value has been recommended for prostate stretched-exponential and kurtosis imaging. Our results revealed that higher b-values above 2000 s/mm2 provided no extra diagnostic value. Therefore, maximum b-values of approximately 2000 s/mm2 may be appropriate for prostate stretched-exponential and kurtosis model imaging. The suggestion of Rosenkrantz et al., that a maximal b-value of 1500–2000 s/mm2 is suitable for body kurtosis imaging, also supports this conclusion[28]. For the bi-exponential model, <D> provided the highest diagnostic value. The AUC of <D> was comparable in group A and group B and was significantly decreased in group C and D. This phenomenon reminds us to pay attention to the maximum b-value of the bi-exponential model. Some studies showed that the parameters of the bi-exponential model depended heavily on the choice of b-values and that it is critical to select an appropriate range of b-values in prostate imaging to obtain accurate parameters[15, 29]. In contrast to <D>, DDC and Kapp, the diagnostic value of ADC calculated using a maximum b-value of 3200 s/mm2 was superior to those derived using maximum b-values of 2000 s/mm2 and 1000 s/mm2. The increased maximum b-value provided greater contrast between cancerous and non-cancerous tissues, adding to the diagnostic performance. However, when the b-value was higher than 3200 s/mm2, the rapid decrease in the goodness-of-fit (0.871 to 0.800) may outweigh the extra benefit of ADC.

Non-Gaussian diffusion models have been proposed to describe complicated water diffusion behavior and are believed to provide a more accurate description of the DWI signal decay curves obtained using higher b-values [30]. Our results showed that the diagnostic performance of parameters derived from a bi-exponential model was low despite excellent goodness-of-fit. The better fit of the bi-exponential model may result from the model addressing a large number of free parameters. However, the additional parameters may cause over-fitting of the data, resulting in poor repeatability and reliability [12, 13, 30, 31]. A study by Merisaari et al. found that, compared with the mono-exponential, stretched-exponential, and kurtosis models, the parameters of the bi-exponential model demonstrated the worst repeatability and diagnostic performance [26]. The poor repeatability and reliability of the bi-exponential model limits its clinical utility. When calculated with a high b-value (≥2000 mm2/s), the ROC analysis in our study did not reveal significantly superior performance of the parameters derived from the stretched-exponential and kurtosis models over the conventional ADC for PCa detection, although a higher-quality fit to DWI signal decay was provided. These results are consistent with several previous studies [13, 32, 33]. High reliability and repeatability are preconditions for using parameters derived from quantitative diffusion models for disease diagnosis and characterization, and this may be why the mono-exponential model is superior at maintaining robust diagnostic performance. Nevertheless, the additional value for PCa of non-Gaussian diffusion models compared with the mono-exponential model remains controversial and requires further research. It is well known that the mono-exponential, bi-exponential, stretched-exponential, and kurtosis models all belong to phenomenological DWI signal models. However, until now, their physiologic basis remained unclear, necessitating further evaluation of the physiologic basis of DWI[34].

There were several limitations in this study. First, this study took the TRUS-guided biopsy as the standard of reference. The use of whole-mount histopathology would improve the accuracy of the agreement between MR images and histopathology. However, it is unreasonable to expect all of our subjects to undergo prostatectomy. Because TRUS-guided biopsy has good specificity but poor sensitivity, a small tumor focus may be missed, resulting in a false-negative diagnosis. In our study, we focused primarily on large tumors, which could partially obviate this problem. Moreover, TRUS-guided biopsy has been the standard of reference in a large number of studies [11, 19, 35, 36], which also illustrated the feasibility of the method used in our study. Second, the role of DWI parameters in the assessment of the transition zone was not evaluated, because TRUS-guided biopsy may miss some cancers that arise in the transition zone and because 70–75% of PCa are located in the PZ. Additionally, DWI is the dominant sequence for PZ PCa detection, while the dominant sequence for transition zone PCa is T2WI[5]; thus, it is more important to evaluate the diagnostic performance of DWI for PZ PCa. Third, compared with whole mount step-section histopathology, TRUS-guided biopsy may lead to the under-grading of a fraction of cancers, therefore this study did not aim to explore the predictive accuracy of different mathematical models with regard to Gleason score. Fourth, the influence of SNR on diffusion model parameters of different b-value ranges was not evaluated. The issue of low SNR must be considered for diffusion imaging, especially for high-b-value images. Although high b-values up to 4500s/mm2 were used in our study, the SNR of the highest images still reached an extremely high level. The following strategies may have played a positive role in increasing SNR. First, we included imaging at a higher field strength (3.0 Tesla). Second, we increased the number of excitations (10 for the highest b-value image). Third, we used a slice thickness of 5 mm, which is thicker than the 3mm commonly used in clinical settings. Fourth, the resolution we used was 3.125×2.917mm rather than 2.5×2.5mm, which is routinely used in clinical settings. Despite the limitations, we believe our methodical strategies provide sufficient validity for the principal results of our study.

Conclusion

In conclusion, our study suggested that ADC derived from a conventional mono-exponential model using a high b-value (3200s /mm2) DWI is an optimal choice for clinical routine application in PZ PCa detection. For PCa detection, the diagnostic value of parameters derived from non-mono-exponential diffusion models or higher b-values (> 3200 s /mm2) warrant further multicenter exploration. Furthermore, as highly parameterized non-Gaussian diffusion models have a higher information content than the single-parameter mono-exponential model, they may provide useful information for differentiating PCa from other prostatic lesions and might be promising for monitoring PCa progression or response to therapy, which also needs further evaluation.

Supporting information

S1 Data. Computed parameters for benign and cancerous tissues.

https://doi.org/10.1371/journal.pone.0172127.s001

(XLS)

S1 Table. Multiple comparisons for each parameters in Table 2.

https://doi.org/10.1371/journal.pone.0172127.s002

(DOC)

Author Contributions

  1. Conceptualization: ZF XM.
  2. Data curation: ZF XM.
  3. Formal analysis: XM CD YC.
  4. Investigation: ZF XM BL ZK PZ.
  5. Methodology: ZF XM CD YC.
  6. Resources: ZF XM LW.
  7. Software: XM CD YC.
  8. Supervision: LW BL ZK PZ.
  9. Validation: ZF XM.
  10. Visualization: XM.
  11. Writing – original draft: ZF.
  12. Writing – review & editing: ZF DM LW VS NC.

References

  1. 1. Futterer JJ, Briganti A, De Visschere P, Emberton M, Giannarini G, Kirkham A, et al. Can Clinically Significant Prostate Cancer Be Detected with Multiparametric Magnetic Resonance Imaging? A Systematic Review of the Literature. European urology. 2015;68(6):1045–53. pmid:25656808
  2. 2. Vargas HA, Wassberg C, Akin O, Hricak H. MR imaging of treated prostate cancer. Radiology. 2012;262(1):26–42. pmid:22190655
  3. 3. Kobus T, Vos PC, Hambrock T, De Rooij M, Hulsbergen-Van de Kaa CA, Barentsz JO, et al. Prostate cancer aggressiveness: in vivo assessment of MR spectroscopy and diffusion-weighted imaging at 3 T. Radiology. 2012;265(2):457–67. pmid:22843767
  4. 4. Weinreb JC, Barentsz JO, Choyke PL, Cornud F, Haider MA, Macura KJ, et al. PI-RADS Prostate Imaging—Reporting and Data System: 2015, Version 2. European urology. 2015.
  5. 5. Hoeks CM, Barentsz JO, Hambrock T, Yakar D, Somford DM, Heijmink SW, et al. Prostate cancer: multiparametric MR imaging for detection, localization, and staging. Radiology. 2011;261(1):46–66. Epub 2011/09/21. pmid:21931141
  6. 6. Peng Y, Jiang Y, Yang C, Brown JB, Antic T, Sethi I, et al. Quantitative analysis of multiparametric prostate MR images: differentiation between prostate cancer and normal tissue and correlation with Gleason score—a computer-aided diagnosis development study. Radiology. 2013;267(3):787–96. Epub 2013/02/09. pmid:23392430
  7. 7. Jensen JH, Helpern JA. MRI quantification of non-Gaussian water diffusion by kurtosis analysis. NMR Biomed. 2010;23(7):698–710. Epub 2010/07/16. PubMed Central PMCID: PMCPMC2997680. pmid:20632416
  8. 8. Valerio M, Zini C, Fierro D, Giura F, Colarieti A, Giuliani A, et al. 3T multiparametric MRI of the prostate: Does intravoxel incoherent motion diffusion imaging have a role in the detection and stratification of prostate cancer in the peripheral zone? European journal of radiology. 2016;85(4):790–4. pmid:26971425
  9. 9. Shinmoto H, Tamura C, Soga S, Shiomi E, Yoshihara N, Kaji T, et al. An intravoxel incoherent motion diffusion-weighted imaging study of prostate cancer. AJR Am J Roentgenol. 2012;199(4):W496–500. Epub 2012/09/22. pmid:22997399
  10. 10. Liu X, Zhou L, Peng W, Wang H, Zhang Y. Comparison of stretched-Exponential and monoexponential model diffusion-Weighted imaging in prostate cancer and normal tissues. Journal of magnetic resonance imaging: JMRI. 2015;42(4):1078–85. pmid:25727776
  11. 11. Rosenkrantz AB, Sigmund EE, Johnson G, Babb JS, Mussi TC, Melamed J, et al. Prostate cancer: feasibility and preliminary experience of a diffusional kurtosis model for detection and assessment of aggressiveness of peripheral zone cancer. Radiology. 2012;264(1):126–35. pmid:22550312
  12. 12. Merisaari H, Movahedi P, Perez IM, Toivonen J, Pesola M, Taimen P, et al. Fitting methods for intravoxel incoherent motion imaging of prostate cancer on region of interest level: Repeatability and gleason score prediction. Magn Reson Med. 2016. Epub 2016/03/01.
  13. 13. Toivonen J, Merisaari H, Pesola M, Taimen P, Bostrom PJ, Pahikkala T, et al. Mathematical models for diffusion-weighted imaging of prostate cancer using b values up to 2000 s/mm: Correlation with Gleason score and repeatability of region of interest analysis. Magnetic resonance in medicine. 2014.
  14. 14. Thormer G, Otto J, Reiss-Zimmermann M, Seiwerts M, Moche M, Garnov N, et al. Diagnostic value of ADC in patients with prostate cancer: influence of the choice of b values. Eur Radiol. 2012;22(8):1820–8. Epub 2012/04/25. pmid:22527373
  15. 15. Pang Y, Turkbey B, Bernardo M, Kruecker J, Kadoury S, Merino MJ, et al. Intravoxel incoherent motion MR imaging for prostate cancer: an evaluation of perfusion fraction and diffusion coefficient derived from different b-value combinations. Magn Reson Med. 2013;69(2):553–62. Epub 2012/04/11. PubMed Central PMCID: PMCPMC3413736. pmid:22488794
  16. 16. Mazzoni LN, Lucarini S, Chiti S, Busoni S, Gori C, Menchi I. Diffusion-weighted signal models in healthy and cancerous peripheral prostate tissues: comparison of outcomes obtained at different b-values. J Magn Reson Imaging. 2014;39(3):512–8. Epub 2013/06/01. pmid:23723087
  17. 17. Zhang K, Shen Y, Zhang X, Ma L, Wang H, An N, et al. Predicting Prostate Biopsy Outcomes: A Preliminary Investigation on Screening with Ultrahigh B-Value Diffusion-Weighted Imaging as an Innovative Diagnostic Biomarker. PloS one. 2016;11(3):e0151176. Epub 2016/03/11. PubMed Central PMCID: PMCPMC4786278. pmid:26963936
  18. 18. Bai Y, Lin Y, Tian J, Shi D, Cheng J, Haacke EM, et al. Grading of Gliomas by Using Monoexponential, Biexponential, and Stretched Exponential Diffusion-weighted MR Imaging and Diffusion Kurtosis MR Imaging. Radiology. 2016;278(2):496–504. Epub 2015/08/01. PubMed Central PMCID: PMCPMC4734159. pmid:26230975
  19. 19. Li B, Cai W, Lv D, Guo X, Zhang J, Wang X, et al. Comparison of MRS and DWI in the diagnosis of prostate cancer based on sextant analysis. Journal of magnetic resonance imaging: JMRI. 2013;37(1):194–200. pmid:23002033
  20. 20. Kitajima K, Kaji Y, Kuroda K, Sugimura K. High b-value diffusion-weighted imaging in normal and malignant peripheral zone tissue of the prostate: effect of signal-to-noise ratio. Magnetic resonance in medical sciences: MRMS: an official journal of Japan Society of Magnetic Resonance in Medicine. 2008;7(2):93–9. Epub 2008/07/08.
  21. 21. Winfield JM, Orton MR, Collins DJ, Ind TE, Attygalle A, Hazell S, et al. Separation of type and grade in cervical tumours using non-mono-exponential models of diffusion-weighted MRI. European radiology. 2016.
  22. 22. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988;44(3):837–45. pmid:3203132
  23. 23. Barentsz JO, Richenberg J, Clements R, Choyke P, Verma S, Villeirs G, et al. ESUR prostate MR guidelines 2012. Eur Radiol. 2012;22(4):746–57. PubMed Central PMCID: PMC3297750. pmid:22322308
  24. 24. Riva-Amarante E, Jimenez-Huete A, Toledano R, Calero M, Alvarez-Linera J, Escribano J, et al. Usefulness of high b-value diffusion-weighted MRI in the diagnosis of Creutzfeldt-Jakob disease. Neurologia (Barcelona, Spain). 2011;26(6):331–6. Epub 2011/02/25.
  25. 25. Han C, Zhao L, Zhong S, Wu X, Guo J, Zhuang X, et al. A comparison of high b-value vs standard b-value diffusion-weighted magnetic resonance imaging at 3.0 T for medulloblastomas. Br J Radiol. 2015;88(1054):20150220. Epub 2015/08/04. PubMed Central PMCID: PMCPMC4730971. pmid:26235143
  26. 26. Merisaari H, Toivonen J, Pesola M, Taimen P, Bostrom PJ, Pahikkala T, et al. Diffusion-weighted imaging of prostate cancer: effect of b-value distribution on repeatability and cancer characterization. Magn Reson Imaging. 2015;33(10):1212–8. pmid:26220861
  27. 27. Merisaari H, Jambor I. Optimization of b-value distribution for four mathematical models of prostate cancer diffusion-weighted imaging using b values up to 2000 s/mm(2): simulation and repeatability study. Magn Reson Med. 2015;73(5):1954–69. Epub 2014/07/22. pmid:25045885
  28. 28. Rosenkrantz AB, Padhani AR, Chenevert TL, Koh DM, De Keyzer F, Taouli B, et al. Body diffusion kurtosis imaging: Basic principles, applications, and considerations for clinical practice. J Magn Reson Imaging. 2015.
  29. 29. Merisaari H, Jambor I. Optimization of b-value distribution for four mathematical models of prostate cancer diffusion-weighted imaging using b values up to 2000 s/mm(2): Simulation and repeatability study. Magnetic resonance in medicine. 2015;73(5):1954–69. pmid:25045885
  30. 30. Jambor I, Merisaari H, Taimen P, Bostrom P, Minn H, Pesola M, et al. Evaluation of different mathematical models for diffusion-weighted imaging of normal prostate and prostate cancer using high b-values: a repeatability study. Magn Reson Med. 2015;73(5):1988–98. pmid:25046482
  31. 31. Taouli B, Beer AJ, Chenevert T, Collins D, Lehman C, Matos C, et al. Diffusion-weighted imaging outside the brain: Consensus statement from an ISMRM-sponsored workshop. Journal of magnetic resonance imaging: JMRI. 2016.
  32. 32. Roethke MC, Kuder TA, Kuru TH, Fenchel M, Hadaschik BA, Laun FB, et al. Evaluation of Diffusion Kurtosis Imaging Versus Standard Diffusion Imaging for Detection and Grading of Peripheral Zone Prostate Cancer. Invest Radiol. 2015;50(8):483–9. Epub 2015/04/14. pmid:25867657
  33. 33. Quentin M, Blondin D, Klasen J, Lanzman RS, Miese FR, Arsov C, et al. Comparison of different mathematical models of diffusion-weighted prostate MR imaging. Magn Reson Imaging. 2012;30(10):1468–74. Epub 2012/07/24. pmid:22819178
  34. 34. Bourne R, Panagiotaki E. Limitations and Prospects for Diffusion-Weighted MRI of the Prostate. Diagnostics. 2016;6(2). PubMed Central PMCID: PMC4931416.
  35. 35. Fusco R, Sansone M, Petrillo M, Setola SV, Granata V, Botti G, et al. Multiparametric MRI for prostate cancer detection: Preliminary results on quantitative analysis of dynamic contrast enhanced imaging, diffusion-weighted imaging and spectroscopy imaging. Magnetic resonance imaging. 2016;34(7):839–45. pmid:27071309
  36. 36. Tamada T, Sone T, Higashi H, Jo Y, Yamamoto A, Kanki A, et al. Prostate cancer detection in patients with total serum prostate-specific antigen levels of 4–10 ng/mL: diagnostic efficacy of diffusion-weighted imaging, dynamic contrast-enhanced MRI, and T2-weighted imaging. AJR Am J Roentgenol. 2011;197(3):664–70. pmid:21862809