Mathematical models for diffusion-weighted imaging of prostate cancer using b values up to 2000 s/mm2 : Correlation with Gleason score and repeatability of region of interest analysis

A Pilot Study on Patient-specific Computational Forecasting of Prostate Cancer Growth during Active Surveillance Using an Imaging-informed Biomechanistic Model

Active surveillance (AS) is a suitable management option for newly diagnosed prostate cancer, which usually presents low to intermediate clinical risk. Patients enrolled in AS have their tumor monitored via longitudinal multiparametric MRI (mpMRI), PSA tests, and biopsies. Hence, treatment is prescribed when these tests identify progression to higher-risk prostate cancer. However, current AS protocols rely on detecting tumor progression through direct observation according to population-based monitoring strategies. This approach limits the design of patient-specific AS plans and may delay the detection of tumor progression. Here, we present a pilot study to address these issues by leveraging personalized computational predictions of prostate cancer growth. Our forecasts are obtained with a spatiotemporal biomechanistic model informed by patient-specific longitudinal mpMRI data (T2-weighted MRI and apparent diffusion coefficient maps from diffusion-weighted MRI). Our results show that our technology can represent and forecast the global tumor burden for individual patients, achieving concordance correlation coefficients from 0.93 to 0.99 across our cohort (n = 7). In addition, we identify a model-based biomarker of higher-risk prostate cancer: the mean proliferation activity of the tumor (P = 0.041). Using logistic regression, we construct a prostate cancer risk classifier based on this biomarker that achieves an area under the ROC curve of 0.83. We further show that coupling our tumor forecasts with this prostate cancer risk classifier enables the early identification of prostate cancer progression to higher-risk disease by more than 1 year. Thus, we posit that our predictive technology constitutes a promising clinical decision-making tool to design personalized AS plans for patients with prostate cancer.

SIGNIFICANCE

Personalization of a biomechanistic model of prostate cancer with mpMRI data enables the prediction of tumor progression, thereby showing promise to guide clinical decision-making during AS for each individual patient.

Preoperative prediction of VETC in hepatocellular carcinoma using non-Gaussian diffusion-weighted imaging at high b values: a pilot study

Background

Vessels encapsulating tumor clusters (VETC) have been considered an important cause of hepatocellular carcinoma (HCC) metastasis.

Purpose

To compare the potential of various diffusion parameters derived from the monoexponential model and four non-Gaussian models (DKI, SEM, FROC, and CTRW) in preoperatively predicting the VETC of HCC.

Methods

86 HCC patients (40 VETC-positive and 46 VETC-negative) were prospectively enrolled. Diffusion-weighted images were acquired using six b-values (range from 0 to 3000 s/mm2). Various diffusion parameters derived from diffusion kurtosis (DK), stretched-exponential (SE), fractional-order calculus (FROC), and continuous-time random walk (CTRW) models, together with the conventional apparent diffusion coefficient (ADC) derived from the monoexponential model were calculated. All parameters were compared between VETC-positive and VETC-negative groups using an independent sample t-test or Mann-Whitney U test, and then the parameters with significant differences between the two groups were combined to establish a predictive model by binary logistic regression. Receiver operating characteristic (ROC) analyses were used to assess diagnostic performance.

Results

Among all studied diffusion parameters, only DKI_K and CTRW_α significantly differed between groups (P=0.002 and 0.004, respectively). For predicting the presence of VETC in HCC patients, the combination of DKI_K and CTRW_α had the larger area under the ROC curve (AUC) than the two parameters individually (AUC=0.747 vs. 0.678 and 0.672, respectively).

Conclusion

DKI_K and CTRW_α outperformed traditional ADC for predicting the VETC of HCC.

Truly reproducible uniform estimation of the ADC with multi-b diffusion data- Application in prostate diffusion imaging

PURPOSE

The ADC is a well-established parameter for clinical diagnostic applications, but lacks reproducibility because it is also influenced by the choice diffusion weighting level. A framework is evaluated that is based on multi-b measurement over a wider range of diffusion-weighting levels and higher order tissue diffusion modeling with retrospective, fully reproducible ADC calculation.

METHODS

Averaging effect from curve fitting for various model functions at 20 linearly spaced b-values was determined by means of simulations and theoretical calculations. Simulation and patient multi-b image data were used to compare the new approach for diffusion-weighted image and ADC map reconstruction with and without Rician bias correction to an active clinical trial protocol probing three non-zero b-values.

RESULTS

Averaging effect at a certain b-value varies for model function and maximum b-value used. Images and ADC maps from the novel procedure are on-par with the clinical protocol. Higher order modeling and Rician bias correction is feasible, but comes at the cost of longer computation times.

CONCLUSIONS

Application of the new framework makes higher order modeling more feasible in a clinical setting while still providing patient images and reproducible ADC maps of adequate quality.

Comparison Between Amide Proton Transfer Magnetic Resonance Imaging Using 3-Dimensional Acquisition and Diffusion-Weighted Imaging for Characterization of Prostate Cancer: A Preliminary Study

OBJECTIVE

This study aimed to compare diagnostic performance for tumor detection and for assessment of tumor aggressiveness in prostate cancer (PC) between amide proton transfer magnetic resonance imaging (MRI) with 3-dimensional acquisition (3DAPT) and diffusion-weighted imaging.

METHODS

The subjects were 23 patients with 27 pathologically proven PCs who underwent 3T multiparametric MRI. With reference to the pathology findings, 2 readers in consensus identified the location of PC on multiparametric MRI and measured APT signal intensity (APT SI [%]) and mean apparent diffusion coefficient (ADC) of the benign region and each PC lesion.

RESULTS

The mean ADC showed a significant difference between benign regions and PC lesions (0.74 ± 0.15 vs 1.37 ± 0.21, P < 0.001), whereas APT SI did not ( P = 0.091). Lesion APT SI was significantly higher and lesion ADC was significantly lower in PCs with Gleason group (GG) ≥3 than in PCs with GG ≤2 (3.37 ± 1.30 vs 1.78 ± 0.67, P < 0.001, and 0.71 ± 0.18 vs 0.79 ± 0.10, P = 0.038, respectively). The APT SI was significantly higher in GG3 than in GG1, in GG3 than in GG2, and in GG4 than in GG2 ( P = 0.009, P = 0.001, and P = 0.006, respectively). The area under the curve for separating tumor lesions and benign regions was 0.601 for 3DAPT and 0.983 for ADC ( P < 0.001). The area under the curve for separating tumors with GG ≤2 from tumors with GG ≥3 was 0.912 for 3DAPT and 0.734 for ADC ( P = 0.172).

CONCLUSIONS

In patients with PC, it might be preferable to use ADC to discriminate benign from malignant tissue and use APT SI for assessment of tumor aggressiveness.

Comparison of Diffusion Kurtosis Imaging and Standard Mono-Exponential Apparent Diffusion Coefficient in Diagnosis of Significant Prostate Cancer-A Correlation with Gleason Score Assessed on Whole-Mount Histopathology Specimens

BACKGROUND

The study was undertaken to compare the diagnostic performance of diffusion kurtosis imaging (DKI) with the standard monoexponential (ME) apparent diffusion coefficient (ADC) model in the detection of significant prostate cancer (PCa), using whole-mount histopathology of radical prostatectomy specimens as a reference standard.

METHODS

155 patients with prostate cancer had undergone multiparametric magnetic resonance imaging (mpMRI) at 3T before prostatectomy. Quantitative diffusion parameters-the apparent diffusion coefficient corrected for non-Gaussian behavior (D_app), kurtosis (K), ADC₁₂₀₀, and ADC₂₀₀₀ were correlated with Gleason score and compared between cancerous and benign tissue and between GS ≤ 3 + 3 and GS ≥ 3 + 4 tumors.

RESULTS

The mean values of all diffusion parameters (D_app, K, ADC₁₂₀₀, ADC₂₀₀₀) were significantly different both between malignant and benign tissue and between GS ≤ 3 + 3 and GS ≥ 3 + 4 tumors. Although the kurtosis model was better fitted to DWI data, the diagnostic performance in receiver operating characteristic (ROC) analysis of DKI and the standard ADC model in the detection of significant PCa was similar in the peripheral zone (PZ) and in peripheral and transitional zones (TZ) together. In conclusion, our study was not able to demonstrate a clear superiority of the kurtosis model over standard ADC in the diagnosis of significant PCa in PZ and in both zones combined.

Introducing new plan evaluation indices for prostate dose painting IMRT plans based on apparent diffusion coefficient images

BACKGROUND

Dose painting planning would be more complicated due to different levels of prescribed doses and more complex evaluation with conventional plan quality indices considering uniform dose prescription. Therefore, we tried to introduce new indices for evaluating the dose distribution conformity and homogeneity of treatment volumes based on the tumoral cell density and relative volumes of each lesion in prostate IMRT.

METHODS

CT and MRI scans of 20 male patients having local prostate cancer were used for IMRT DP planning. Apparent diffusion coefficient (ADC) images were imported to a MATLAB program to identify lesion regions based on ADC values automatically. Regions with ADC values lower than 750 mm²/s and regions with ADC values higher than 750 and less than 1500 mm²/s were considered CTV_70Gy (clinical tumor volume with 70 Gy prescribed dose), and CTV_60Gy, respectively. Other regions of the prostate were considered as CTV_53Gy. New plan evaluation indices based on evaluating the homogeneity (IOE(H)), and conformity (IOE(C)) were introduced, considering the relative volume of each lesion and cellular density obtained from ADC images. These indices were compared with conventional homogeneity and conformity indices and IOEs without considering cellular density. Furthermore, tumor control probability (TCP) was calculated for each patient, and the relationship of the assessed indices were evaluated with TCP values.

RESULTS

IOE (H) and IOE (C) with considering cellular density had significantly lower values compared to conventional indices and IOEs without considering cellular density. (P < 0.05). TCP values had a stronger relationship with IOE(H) considering cell density (R² = -0.415), and IOE(C) without considering cell density (R² = 0.624).

CONCLUSION

IOE plan evaluation indices proposed in this study can be used for evaluating prostate IMRT dose painting plans. We suggested to consider cell densities in the IOE(H) calculation formula and it's appropriate to calculate IOE(C) without considering cell density values.

Quicksort leave-pair-out cross-validation for ROC curve analysis

Receiver Operating Characteristic (ROC) curve analysis and area under the ROC curve (AUC) are commonly used performance measures in diagnostic systems. In this work, we assume a setting, where a classifier is inferred from multivariate data to predict the diagnostic outcome for new cases. Cross-validation is a resampling method for estimating the prediction performance of a classifier on data not used for inferring it. Tournament leave-pair-out (TLPO) cross-validation has been shown to be better than other resampling methods at producing a ranking of data that can be used for estimating the ROC curves and areas under them. However, the time complexity of TLPOCV, $$O\left( n^2\right)$$ O n 2 , means that it is impractical in many applications. In this article, a method called quicksort leave-pair-out cross-validation (QLPOCV) is presented in order to decrease the time complexity of obtaining a reliable ranking of data to $$O\left( n\log n\right)$$ O n log n . The proposed method is compared with existing ones in an experimental study, demonstrating that in terms of ROC curves and AUC values QLPOCV produces as accurate performance estimation as TLPOCV, outperforming both k-fold and leave-one-out cross-validation.

Diffusion weighted imaging and diffusion kurtosis imaging in abdominal oncological setting: why and when

This article provides an overview of diffusion kurtosis (DKI) imaging in abdominal oncology. DKI allows for more data on tissue structures than the conventional diffusion model (DWI). However, DKI requires high quality images at b-values greater than 1000 s/mm² and high signal-to-noise ratio (SNR) that traditionally MRI systems are not able to acquire and therefore there are generally amplified anatomical distortions on the images due to less homogeneity of the field. Advances in both hardware and software on modern MRI scanners have currently enabled ultra-high b-value imaging and offered the ability to apply DKI to multiple extracranial sites. Previous studies have evaluated the ability of DKI to characterize and discriminate tumor grade compared to conventional DWI. Additionally, in several studies the DKI sequences used were based on planar echo (EPI) acquisition, which is susceptible to motion, metal and air artefacts and prone to low SNRs and distortions, leading to low quality images for some small lesions, which may affect the accuracy of the results. Another problem is the optimal b-value of DKI, which remains to be explored and not yet standardized, as well as the manual selection of the ROI, which could affect the accuracy of some parameters.

Diffusion and quantification of diffusion of prostate cancer

For assessing a cancer treatment, and for detecting and characterizing cancer, Diffusion-weighted imaging (DWI) is commonly used. The key in DWI's use extracranially has been due to the emergence of of high-gradient amplitude and multichannel coils, parallelimaging, and echo-planar imaging. The benefit has been fewer motion artefacts and high-quality prostate images.Recently, new techniques have been developed to improve the signal-to-noise ratio of DWI with fewer artefacts, allowing an increase in spatial resolution. For apparent diffusion coefficient quantification, non-Gaussian diffusion models have been proposed as additional tools for prostate cancer detection and evaluation of its aggressiveness. More recently, radiomics and machine learning for prostate magnetic resonance imaging have emerged as novel techniques for the non-invasive characterisation of prostate cancer. This review presents recent developments in prostate DWI and discusses its potential use in clinical practice.

How to read biparametric MRI in men with a clinical suspicious of prostate cancer: Pictorial review for beginners with public access to imaging, clinical and histopathological database

Prostate Magnetic Resonance Imaging (MRI) is increasingly being used in men with a clinical suspicion of prostate cancer (PCa). Performing prostate MRI without the use of an intravenous contrast (IV) agent in men with a clinical suspicion of PCa can lead to reduced MRI scan time. Enabling a large array of different medical providers (from mid-level to specialized radiologists) to evaluate and potentially report prostate MRI in men with a clinical suspicion of PCa with a high accuracy could be one way to enable wide adoption of prostate MRI in men with a clinical suspicion of PCa. The aim of this pictorial review is to provide an insight into acquisition, quality control and reporting of prostate MRI performed without IV contrast agent in men with a clinical suspicion of PCa, aimed specifically at radiologists starting reporting prostate MRI, urologists, urology/radiology residents and mid-level medical providers without experience in reporting prostate MRI. Free public access (http://petiv.utu.fi/improd/and http://petiv.utu.fi/multiimprod/) to complete datasets of 161 and 338 men is provided. The imaging datasets are accompanied by clinical, laboratory and histopathological findings. Several topics are simplified in order to provide a solid base for the development of skills needed for an unsupervised review and potential reporting of prostate MRI in men with a clinical suspicion of PCa. The current review represents the first step towards enabling a large array of different medical providers to review and report accurately prostate MRI performed without IV contrast agent in men with a clinical suspicion of PCa.

Diffusion-weighted imaging in prostate cancer

Diffusion-weighted imaging (DWI), a key component in multiparametric MRI (mpMRI), is useful for tumor detection and localization in clinically significant prostate cancer (csPCa). The Prostate Imaging Reporting and Data System versions 2 and 2.1 (PI-RADS v2 and PI-RADS v2.1) emphasize the role of DWI in determining PIRADS Assessment Category in each of the transition and peripheral zones. In addition, several recent studies have demonstrated comparable performance of abbreviated biparametric MRI (bpMRI), which incorporates only T2-weighted imaging and DWI, compared with mpMRI with dynamic contrast-enhanced MRI. Therefore, further optimization of DWI is essential to achieve clinical application of bpMRI for efficient detection of csPC in patients with elevated PSA levels. Although DWI acquisition is routinely performed using single-shot echo-planar imaging, this method suffers from such as susceptibility artifact and anatomic distortion, which remain to be solved. In this review article, we will outline existing problems in standard DWI using the single-shot echo-planar imaging sequence; discuss solutions that employ newly developed imaging techniques, state-of-the-art technologies, and sequences in DWI; and evaluate the current status of quantitative DWI for assessment of tumor aggressiveness in PC.

Quantitative diffusion-weighted imaging and dynamic contrast-enhanced MR imaging for assessment of tumor aggressiveness in prostate cancer at 3T

PURPOSE

To compare diffusion-weighted imaging (DWI) and dynamic contrast-enhanced MR imaging (DCE-MRI) for characterization of prostate cancer (PC).

METHODS

104 PC patients who underwent prostate multiparametric MRI at 3T including DWI and DCE-MRI before MRI-guided biopsy or radical prostatectomy. Apparent diffusion coefficient (ADC) with histogram analysis (mean, 0-25th percentile, skewness, and kurtosis), intravoxel incoherent motion model including D and f; stretched exponential model including distributed diffusion coefficient (DDC) and a; and permeability parameters including K^trans, K_ep, and V_e were obtained from a region of interest placed on the dominant tumor of each patient.

RESULTS

ADC_mean, ADC_0-25, D, DDC, and V_e were significantly lower and K_ep was significantly higher in GS ≥ 3 + 4 tumors (n = 89) than in GS = 3 + 3 tumors (n = 15), and also in GS ≥ 4 + 3 tumors (n = 57) than in GS ≤ 3 + 4 tumors (n = 47) (P < 0.001 to P = 0.040). f was significantly lower in GS ≥ 4 + 3 tumors than in GS ≤ 3 + 4 tumors (P = 0.022), but there was no significant difference between GS = 3 + 3 tumors and GS ≥ 3 + 4 tumors, or between the remaining metrics in both comparisons. In metrics with area under the curve (AUC) >0.80, there was a significant difference in AUC between ADC_0-25 and D, and DDC for separating GS ≤ 3 + 4 tumors from GS ≥ 4 + 3 tumors (P = 0.040 and P = 0.022, respectively). There were no significant differences between metrics with AUC > 0.80 for separating GS = 3 + 3 tumors from GS ≥ 3 + 4 tumors. ADC_0-25 had the highest correlation with Gleason grade (ρ = -0.625, P < 0.001).

CONCLUSIONS

DWI and DCE-MRI showed no apparent clinical superiority of non-Gaussian models or permeability MRI over the mono-exponential model for assessment of tumor aggressiveness in PC.

Artificial intelligence and imaging biomarkers for prostate radiation therapy during and after treatment

Magnetic resonance imaging (MRI) is increasingly used in the management of prostate cancer (PCa). Quantitative MRI (qMRI) parameters, derived from multi-parametric MRI, provide indirect measures of tumour characteristics such as cellularity, angiogenesis and hypoxia. Using Artificial Intelligence (AI), relevant information and patterns can be efficiently identified in these complex data to develop quantitative imaging biomarkers (QIBs) of tumour function and biology. Such QIBs have already demonstrated potential in the diagnosis and staging of PCa. In this review, we explore the role of these QIBs in monitoring treatment response during and after PCa radiotherapy (RT). Recurrence of PCa after RT is not uncommon, and early detection prior to development of metastases provides an opportunity for salvage treatments with curative intent. However, the current method of monitoring treatment response using prostate-specific antigen levels lacks specificity. QIBs, derived from qMRI and developed using AI techniques, can be used to monitor biological changes post-RT providing the potential for accurate and early diagnosis of recurrent disease.

Comparison of diagnostic performance between diffusion kurtosis imaging parameters and mono-exponential ADC for determination of clinically significant cancer in patients with prostate cancer

PURPOSE

To compare the diagnostic performance between diffusion kurtosis imaging (DKI) parameters and mono-exponential apparent diffusion coefficient (ADC) for determination of clinically significant cancer (CSC, Gleason score (GS) ≥ 7) in patients with histologically proven prostate cancer (PCa).

METHODS

A total of 92 patients (mean age: 71.5 years, range: 47-89 years) who had been diagnosed as PCa and undergone 3 T-MRI including DWI (b values, 0, 100, 1000, 2000s/mm²) were included in this study. The DKI parameters, namely apparent diffusion for non-Gaussian distribution (D_app) and apparent kurtosis coefficient (K_app), were calculated by dedicated software using mono-exponential and diffusion kurtosis models for quantitation. The measurement was performed for a whole tumor after segmentation, and pathologic topographic maps or systemic biopsy results served as the reference standard for segmentation. To compare the diagnostic performance of each parameter for determination of CSC, pair-wise comparison of receiver operating characteristic (ROC) curves was performed.

RESULTS

The study population consisted of GS 6 (n = 18), GS 7 (n = 31), GS 8 (n = 25), GS 9 (n = 15) and GS 10 (n = 3) patients. The area under the ROC curve of K_app (0.707, 95% CI 0.603-0.798) for discriminating CSC from non-CSC was not significantly different from those of mono-exponential ADC (0.725, 0.622-0.813, P = 0.2175) or D_app (0.726, 0.623-0.814, P = 0.9628). Diagnostic predictive values of K_app were estimated to a maximum accuracy of 78%, a sensitivity of 86%, and a specificity of 47%, while those of mono-exponential ADC were 75, 81, and 53%, respectively.

CONCLUSION

The DKI parameters showed a diagnostic performance comparable to mono-exponential ADC for determination of CSC in patients with PCa.

Removing rician bias in diffusional kurtosis of the prostate using real-data reconstruction

PURPOSE

To compare prostate diffusional kurtosis imaging (DKI) metrics generated using phase-corrected real data with those generated using magnitude data with and without noise compensation (NC).

METHODS

Diffusion-weighted images were acquired at 3T in 16 prostate cancer patients, measuring 6 b-values (0-1500 s/mm2 ), each acquired with 6 signal averages along 3 diffusion directions, with noise-only images acquired to allow NC. In addition to conventional magnitude averaging, phase-corrected real data were averaged in an attempt to reduce rician noise-bias, with a range of phase-correction low-pass filter (LPF) sizes (8-128 pixels) tested. Each method was also tested using simulations. Pixelwise maps of apparent diffusion (D) and apparent kurtosis (K) were calculated for magnitude data with and without NC and phase-corrected real data. Average values were compared in tumor, normal transition zone (NTZ), and normal peripheral zone (NPZ).

RESULTS

Simulations indicated LPF size can strongly affect K metrics, where 64-pixel LPFs produced accurate metrics. Relative to metrics estimated from magnitude data without NC, median NC K were lower (P < 0.0001) by 6/11/8% in tumor/NPZ/NTZ, 64-LPF real-data K were lower (P < 0.0001) by 4/10/7%, respectively.

CONCLUSION

Compared with magnitude data with NC, phase-corrected real data can produce similar K, although the choice of phase-correction LPF should be chosen carefully.

Histogram analysis from stretched exponential model on diffusion-weighted imaging: evaluation of clinically significant prostate cancer

OBJECTIVE

To evaluate the usefulness of histogram analysis of stretched exponential model (SEM) on diffusion-weighted imaging in evaluating clinically significant prostate cancer (CSC).

METHODS

A total of 85 patients with prostate cancer underwent 3 T multiparametric MRI, followed by radical prostatectomy. Histogram parameters of the tumor from the SEM [distributed diffusion coefficient (DDC) and α] and the monoexponential model [MEM; apparent diffusion coefficient (ADC)] were evaluated. The associations between parameters and Gleason score or Prostate Imaging Reporting and Data System v. 2 were evaluated. The area under the receiver operating characteristics curve was calculated to evaluate diagnostic performance of parameters in predicting CSC.

RESULTS

The values of histogram parameters of DDC and ADC were significantly lower in patients with CSC than in patients without CSC (p < 0.05), except for skewness and kurtosis. The value of the 25th percentile of α was significantly lower in patients with CSC than in patients without CSC (p = 0.014). Histogram parameters of ADC and DDC had significant weak to moderate negative associations with Gleason score or Prostate Imaging Reporting and Data System v. 2 (p < 0.001), except for skewness and kurtosis. For predicting CSC, the area under the curves of mean ADC (0.856), 50th percentile DDC (0.852), and 25th percentile α (0.707) yielded the highest values compared to other histogram parameters from each group.

CONCLUSION

Histogram analysis of the SEM on diffusion-weighted imaging may be a useful quantitative tool for evaluating CSC. However, the SEM did not outperform the MEM.

ADVANCES IN KNOWLEDGE

Histogram parameters of SEM may be useful for evaluating CSC.

Non-Gaussian models of diffusion weighted imaging for detection and characterization of prostate cancer: a systematic review and meta-analysis

The importance of Diffusion Weighted Imaging (DWI) in prostate cancer (PCa) diagnosis have been widely handled in literature. In the last decade, due to the mono-exponential model limitations, several studies investigated non-Gaussian DWI models and their utility in PCa diagnosis. Since their results were often inconsistent and conflicting, we performed a systematic review of studies from 2012 examining the most commonly used Non-Gaussian DWI models for PCa detection and characterization. A meta-analysis was conducted to assess the ability of each Non-Gaussian model to detect PCa lesions and distinguish between low and intermediate/high grade lesions. Weighted mean differences and 95% confidence intervals were calculated and the heterogeneity was estimated using the I2 statistic. 29 studies were selected for the systematic review, whose results showed inconsistence and an unclear idea about the actual usefulness and the added value of the Non-Gaussian model parameters. 12 studies were considered in the meta-analyses, which showed statistical significance for several non-Gaussian parameters for PCa detection, and to a lesser extent for PCa characterization. Our findings showed that Non-Gaussian model parameters may potentially play a role in the detection and characterization of PCa but further studies are required to identify a standardized DWI acquisition protocol for PCa diagnosis.

Repeatability of radiomics and machine learning for DWI: Short-term repeatability study of 112 patients with prostate cancer

PURPOSE

To evaluate repeatability of prostate DWI-derived radiomics and machine learning methods for prostate cancer (PCa) characterization.

METHODS

A total of 112 patients with diagnosed PCa underwent 2 prostate MRI examinations (Scan1 and Scan2) performed on the same day. DWI was performed using 12 b-values (0-2000 s/mm² ), post-processed using kurtosis function, and PCa areas were annotated using whole mount prostatectomy sections. A total of 1694 radiomic features including Sobel, Kirch, Gradient, Zernike Moments, Gabor, Haralick, CoLIAGe, Haar wavelet coefficients, 3D analogue to Laws features, 2D contours, and corner detectors were calculated. Radiomics and 4 feature pruning methods (area under the receiver operator characteristic curve, maximum relevance minimum redundancy, Spearman's ρ, Wilcoxon rank-sum) were evaluated in terms of Scan1-Scan2 repeatability using intraclass correlation coefficient (ICC)(3,1). Classification performance for clinically significant and insignificant PCa with Gleason grade groups 1 versus >1 was evaluated by area under the receiver operator characteristic curve in unseen random 30% data split.

RESULTS

The ICC(3,1) values for conventional radiomics and feature pruning methods were in the range of 0.28-0.90. The machine learning classifications varied between Scan1 and Scan2 with % of same class labels between Scan1 and Scan2 in the range of 61-81%. Surface-to-volume ratio and corner detector-based features were among the most represented features with high repeatability, ICC(3,1) >0.75, consistently high ranking using all 4 feature pruning methods, and classification performance with area under the receiver operator characteristic curve >0.70.

CONCLUSION

Surface-to-volume ratio and corner detectors for prostate DWI led to good classification of unseen data and performed similarly in Scan1 and Scan2 in contrast to multiple conventional radiomic features.

Prostate Cancer Risk Stratification in Men With a Clinical Suspicion of Prostate Cancer Using a Unique Biparametric MRI and Expression of 11 Genes in Apparently Benign Tissue: Evaluation Using Machine-Learning Techniques

BACKGROUND

Accurate risk stratification of men with a clinical suspicion of prostate cancer (cSPCa) remains challenging despite the increasing use of MRI.

PURPOSE

To evaluate the diagnostic accuracy of a unique biparametric MRI protocol (IMPROD bpMRI) combined with clinical and molecular markers in men with cSPCa.

STUDY TYPE

Prospective single-institutional clinical trial (NCT01864135).

SUBJECTS

Eighty men with cSPCa.

FIELD STRENGTH/SEQUENCE

3T, surface array coils. Two T₂ -weighted and three diffusion-weighted imaging (DWI) acquisitions: 1) b-values 0, 100, 200, 300, 500 s/mm² ; 2) b-values 0,1500 s/mm² ; 3) b-values 0, 2000 s/mm² .

ASSESSMENT

IMPROD bpMRI examinations were qualitatively (IMPROD bpMRI Likert score) and quantitatively (DWI-based Gleason grade score) prospectively reported. Men with IMPROD bpMRI Likert 3-5 had two targeted biopsies followed by 12-core systematic biopsies (SB); those with IMPROD bpMRI Likert 1-2 had only SB. Additionally, 2-core from normal-appearing prostate areas were obtained for the mRNA expression of ACSM1, AMACR, CACNA1D, DLX1, PCA3, PLA2G7, RHOU, SPINK1, SPON2, TMPRSS2-ERG, and TDRD1 measured by quantitative reverse-transcription polymerase chain reaction.

STATISTICAL TESTS

Univariate and multivariate analysis using regularized least-squares, feature selection and tournament leave-pair-out cross-validation (TLPOCV), as well as 10 random splits of the data in training-testing sets, were used to evaluate the mRNA, clinical and IMPROD bpMRI parameters in detecting clinically significant prostate cancer (SPCa) defined as Gleason score ≥ 3 + 4. The evaluation metric was the area under the curve (AUC).

RESULTS

IMPROD bpMRI Likert demonstrated the highest TLPOCV AUC of 0.92. The tested clinical variables had AUC 0.56-0.73, while the mRNA and additional IMPROD bpMRI parameters had AUC 0.50-0.67 and 0.65-0.89 respectively. The combination of clinical and mRNA biomarkers produced TLPOCV AUC of 0.87, the highest TLPOCV performance without including IMPROD bpMRI Likert.

DATA CONCLUSION

The qualitative IMPROD bpMRI Likert score demonstrated the highest accuracy for SPCa detection compared with the tested clinical variables and mRNA biomarkers.

LEVEL OF EVIDENCE

1 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2020;51:1540-1553.

Performance of Diffusion Kurtosis Imaging Versus Diffusion Tensor Imaging in Discriminating Between Benign Tissue, Low and High Gleason Grade Prostate Cancer

RATIONALE AND OBJECTIVES

To investigate the performance of diffusion kurtosis imaging (DKI) and diffusion tensor imaging (DTI) in discriminating benign tissue, low- and high-grade prostate adenocarcinoma (PCa).

MATERIALS AND METHODS

Forty-eight patients with biopsy-proven PCa of different Gleason grade (GG), who provided written informed consent, were enrolled. All subjects underwent 3T DWI examinations by using b values 0, 500, 1000, 1500, 2000, and 2500 s/mm² and six gradient directions. Mean diffusivity, fractional anisotropy (FA), apparent kurtosis (K), apparent kurtosis-derived diffusivity (D), and proxy fractional kurtosis anisotropy (KFA) maps were obtained. Regions of interest were selected in PCa, in the contralateral benign zone, and in the peritumoral area. Histogram analysis was performed by measuring mean, 10th, 25th, and 90th (p90) percentile of the whole-lesion volume. Kruskal-Wallis test with Bonferroni correction was used to assess significant differences between different regions of interest. The correlation between diffusion metrics and GG and between DKI and DTI parameters was evaluated with Pearson's test. ROC curve analysis was carried out to analyze the ability of histogram variables to differentiate low- and high-GG PCa.

RESULTS

All metrics significantly discriminated PCa from benign and from peritumoral tissue (except for K, KFA_p90, and FA). K_p90 showed the highest correlation with GG and the best diagnostic ability (area under the curve = 0.84) in discriminating low- from high-risk PCa.

CONCLUSION

Compared to DTI, DKI provides complementary and additional information about prostate cancer tissue, resulting more sensitive to PCa-derived modifications and more accurate in discriminating low- and high-risk PCa.

Radiomics and machine learning of multisequence multiparametric prostate MRI: Towards improved non-invasive prostate cancer characterization

Purpose

To develop and validate a classifier system for prediction of prostate cancer (PCa) Gleason score (GS) using radiomics and texture features of T₂-weighted imaging (T₂w), diffusion weighted imaging (DWI) acquired using high b values, and T₂-mapping (T₂).

Methods

T₂w, DWI (12 b values, 0–2000 s/mm²), and T₂ data sets of 62 patients with histologically confirmed PCa were acquired at 3T using surface array coils. The DWI data sets were post-processed using monoexponential and kurtosis models, while T₂w was standardized to a common scale. Local statistics and 8 different radiomics/texture descriptors were utilized at different configurations to extract a total of 7105 unique per-tumor features. Regularized logistic regression with implicit feature selection and leave pair out cross validation was used to discriminate tumors with 3+3 vs >3+3 GS.

Results

In total, 100 PCa lesions were analysed, of those 20 and 80 had GS of 3+3 and >3+3, respectively. The best model performance was obtained by selecting the top 1% features of T₂w, ADC_m and K with ROC AUC of 0.88 (95% CI of 0.82–0.95). Features from T₂ mapping provided little added value. The most useful texture features were based on the gray-level co-occurrence matrix, Gabor transform, and Zernike moments.

Conclusion

Texture feature analysis of DWI, post-processed using monoexponential and kurtosis models, and T₂w demonstrated good classification performance for GS of PCa. In multisequence setting, the optimal radiomics based texture extraction methods and parameters differed between different image types.

Validation of IMPROD biparametric MRI in men with clinically suspected prostate cancer: A prospective multi-institutional trial

BACKGROUND

Magnetic resonance imaging (MRI) combined with targeted biopsy (TB) is increasingly used in men with clinically suspected prostate cancer (PCa), but the long acquisition times, high costs, and inter-center/reader variability of routine multiparametric prostate MRI limit its wider adoption.

METHODS AND FINDINGS

The aim was to validate a previously developed unique MRI acquisition and reporting protocol, IMPROD biparametric MRI (bpMRI) (NCT01864135), in men with a clinical suspicion of PCa in a multi-institutional trial (NCT02241122). IMPROD bpMRI has average acquisition time of 15 minutes (no endorectal coil, no intravenous contrast use) and consists of T2-weighted imaging and 3 separate diffusion-weighed imaging acquisitions. Between February 1, 2015, and March 31, 2017, 364 men with a clinical suspicion of PCa were enrolled at 4 institutions in Finland. Men with an equivocal to high suspicion (IMPROD bpMRI Likert score 3-5) of PCa had 2 TBs of up to 2 lesions followed by a systematic biopsy (SB). Men with a low to very low suspicion (IMPROD bpMRI Likert score 1-2) had only SB. All data and protocols are freely available. The primary outcome of the trial was diagnostic accuracy-including overall accuracy, sensitivity, specificity, negative predictive value (NPV), and positive predictive value-of IMPROD bpMRI for clinically significant PCa (SPCa), which was defined as a Gleason score ≥ 3 + 4 (Gleason grade group 2 or higher). In total, 338 (338/364, 93%) prospectively enrolled men completed the trial. The accuracy and NPV of IMPROD bpMRI for SPCa were 70% (113/161) and 95% (71/75) (95% CI 87%-98%), respectively. Restricting the biopsy to men with equivocal to highly suspicious IMPROD bpMRI findings would have resulted in a 22% (75/338) reduction in the number of men undergoing biopsy while missing 4 (3%, 4/146) men with SPCa. The main limitation is uncertainty about the true PCa prevalence in the study cohort, since some of the men may have PCa despite having negative biopsy findings.

CONCLUSIONS

IMPROD bpMRI demonstrated a high NPV for SPCa in men with a clinical suspicion of PCa in this prospective multi-institutional clinical trial.

TRIAL REGISTRATION

ClinicalTrials.gov NCT02241122.

An annotated test-retest collection of prostate multiparametric MRI

Multiparametric Magnetic Resonance Imaging (mpMRI) is widely used for characterizing prostate cancer. Standard of care use of mpMRI in clinic relies on visual interpretation of the images by an expert. mpMRI is also increasingly used as a quantitative imaging biomarker of the disease. Little is known about repeatability of such quantitative measurements, and no test-retest datasets have been available publicly to support investigation of the technical characteristics of the MRI-based quantification in the prostate. Here we present an mpMRI dataset consisting of baseline and repeat prostate MRI exams for 15 subjects, manually annotated to define regions corresponding to lesions and anatomical structures, and accompanied by region-based measurements. This dataset aims to support further investigation of the repeatability of mpMRI-derived quantitative prostate measurements, study of the robustness and reliability of the automated analysis approaches, and to support development and validation of new image analysis techniques. The manuscript can also serve as an example of the use of DICOM for standardized encoding of the image annotation and quantification results.

Evaluation of non-monoexponential diffusion models for hepatocellular carcinoma using b values up to 2000 s/mm2 : A short-term repeatability study

BACKGROUND

Non-monoexponential diffusion models are being used increasingly for the characterization and curative effect evaluation of hepatocellular carcinoma (HCC). But the fitting quality of the models and the repeatability of their parameters have not been assessed for HCC.

PURPOSE

To evaluate kurtosis, stretched exponential, and statistical models for diffusion-weighted imaging (DWI) of HCC, using b-values up to 2000 s/mm² , in terms of fitting quality and repeatability.

STUDY TYPE

Prospective.

POPULATION

Eighteen patients with HCC.

FIELD STRENGTH/SEQUENCE

Conventional and DW images (b = 0, 200, 500, 1000, 1500, 2000 s/mm² ) were acquired at 3.0T.

ASSESSMENT

The parameters of the kurtosis, stretched exponential, and statistical models were calculated on regions of interest (ROIs) of each lesion.

STATISTICAL TESTS

The fitting quality was evaluated through comparing the fitting residuals produced on the average data of ROI between different models using a paired t-test or Wilcoxon rank-sum test. Repeatability of the fitted parameters at the median values on the voxelwise data of ROI was assessed using the within coefficient of variation (WCV), the intraclass correlation coefficient (ICC), and the 95% Bland-Altman limits of agreements (BA-LA). The repeatability was divided into four levels: excellent, good, acceptable, and poor, referring to the values of ICC and WCV.

RESULTS

Among three models, the stretched exponential model provided the best fit to HCC (P < 0.05), whereas the statistical model produced the largest fitting residuals (P < 0.05). The repeatability of K from the kurtosis model was excellent (ICC 0.915; WCV 8.79%), while the distributed diffusion coefficient (DDC) from the stretched model was just acceptable (ICC 0.477; WCV 27.83%). The repeatability was good for other diffusion-related parameters.

DATA CONCLUSION

Considering the model fit and repeatability, the kurtosis and stretched exponential models are the preferred models for the description of the DW signals of HCC with respect to the statistical model.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:297-304.

Tournament leave-pair-out cross-validation for receiver operating characteristic analysis

Receiver operating characteristic analysis is widely used for evaluating diagnostic systems. Recent studies have shown that estimating an area under receiver operating characteristic curve with standard cross-validation methods suffers from a large bias. The leave-pair-out cross-validation has been shown to correct this bias. However, while leave-pair-out produces an almost unbiased estimate of area under receiver operating characteristic curve, it does not provide a ranking of the data needed for plotting and analyzing the receiver operating characteristic curve. In this study, we propose a new method called tournament leave-pair-out cross-validation. This method extends leave-pair-out by creating a tournament from pair comparisons to produce a ranking for the data. Tournament leave-pair-out preserves the advantage of leave-pair-out for estimating area under receiver operating characteristic curve, while it also allows performing receiver operating characteristic analyses. We have shown using both synthetic and real-world data that tournament leave-pair-out is as reliable as leave-pair-out for area under receiver operating characteristic curve estimation and confirmed the bias in leave-one-out cross-validation on low-dimensional data. As a case study on receiver operating characteristic analysis, we also evaluate how reliably sensitivity and specificity can be estimated from tournament leave-pair-out receiver operating characteristic curves.

Multiparametric Magnetic Resonance Imaging of the Prostate: Repeatability of Volume and Apparent Diffusion Coefficient Quantification

OBJECTIVES

The aim of this study was to evaluate the repeatability of a region of interest (ROI) volume and mean apparent diffusion coefficient (ADC) in standard-of-care 3 T multiparametric magnetic resonance imaging (mpMRI) of the prostate obtained with the use of endorectal coil.

MATERIALS AND METHODS

This prospective study was Health Insurance Portability and Accountability Act compliant, with institutional review board approval and written informed consent. Men with confirmed or suspected treatment-naive prostate cancer scheduled for mpMRI were offered a repeat mpMRI within 2 weeks. Regions of interest corresponding to the whole prostate gland, the entire peripheral zone (PZ), normal PZ, and suspected tumor ROI (tROI) on axial T2-weighted, dynamic contrast-enhanced subtract, and ADC images were annotated and assessed using Prostate Imaging Reporting and Data System (PI-RADS) v2. Repeatability of the ROI volume for each of the analyzed image types and mean ROI ADC was summarized with repeatability coefficient (RC) and RC%.

RESULTS

A total of 189 subjects were approached to participate in the study. Of 40 patients that gave initial agreement, 15 men underwent 2 mpMRI examinations and completed the study. Peripheral zone tROIs were identified in 11 subjects. Tumor ROI volume was less than 0.5 mL in 8 of 11 subjects. PI-RADS categories were identical between baseline-repeat studies in 11/15 subjects and differed by 1 point in 4/15. Peripheral zone tROI volume RC (RC%) was 233 mm (71%) on axial T2-weighted, 422 mm (112%) on ADC, and 488 mm (119%) on dynamic contrast-enhanced subtract. Apparent diffusion coefficient ROI mean RC (RC%) were 447 × 10 mm/s (42%) in PZ tROI and 471 × 10 mm/s (30%) in normal PZ. Significant difference in repeatability of the tROI volume across series was observed (P < 0.005). The mean ADC RC% was lower than volume RC% for tROI ADC (P < 0.05).

CONCLUSIONS

PI-RADS v2 overall assessment was highly repeatable. Multiparametric magnetic resonance imaging sequences differ in volume measurement repeatability. The mean tROI ADC is more repeatable compared with tROI volume in ADC. Repeatability of prostate ADC is comparable with that in other abdominal organs.

Histogram analysis of stretched-exponential and monoexponential diffusion-weighted imaging models for distinguishing low and intermediate/high gleason scores in prostate carcinoma

BACKGROUND

Noninvasive measures to evaluate the aggressiveness of prostate carcinoma (PCa) may benefit patients.

PURPOSE

To assess the value of stretched-exponential and monoexponential diffusion-weighted imaging (DWI) for predicting the aggressiveness of PCa.

STUDY TYPE

Retrospective study.

SUBJECTS

Seventy-five patients with PCa.

FIELD STRENGTH

3T DWI examinations were performed using b-values of 0, 500, 1000, and 2000 s/mm² .

ASSESSMENT

The research were based on entire-tumor histogram analysis and the reference standard was radical prostectomy.

STATISTICAL TESTS

The correlation analysis was programmed with Spearman's rank-order analysis between the histogram variables and Gleason grade group (GG). Receiver operating characteristic (ROC) regression was used to analyze the ability of these histogram variables to differentiate low-grade (LG) from intermediate/high-grade (HG) PCa.

RESULTS

The percentiles and mean of apparent diffusion coefficient (ADC) and distributed diffusion coefficient (DDC) were correlated with GG (ρ: 0.414-0.593), while there was no significant relation among α value, skewnesses, and kurtosises with GG (ρ:0.034-0.323). HG tumors (ADC:484 ± 136, 592 ± 139, 670 ± 144, 788 ± 146, 895 ± 141 mm² /s; DDC: 410 ± 142, 532 ± 172, 666 ± 193, 786 ± 196, 914 ± 181 mm² /s) had lower values in the 10^th , 25^th , 50^th , 75^th percentiles and means than LG tumors (ADC: 644 ± 779, 737 ± 84, 836 ± 83, 919 ± 82, 997 ± 107 mm² /s; DDC: 552 ± 82, 680 ± 94, 829 ± 112, 931 ± 106, 1045 ± 100 mm² /s). However, there was no difference between LG and HG tumors in α value (0.671 ± 0.041 vs. 0.633 ± 0.114), kurtosises (ADC 0.09 vs. 0.086; DDC -0.033 vs. -0.317), or skewnesses (ADC -0.036 vs. 0.073; DDC -0.063 vs. 0.136). The above statistics were P < 0.01. ADC10 with AUC = 0.840 and DDC10 with AUC = 0.799 were similar in discriminating between LG and HG PCa at P < 0.05.

DATA CONCLUSION

Histogram variables of DDC and ADC may predict the aggressiveness of PCa, while α value does not. The abilities of ADC10 and DDC10 to discriminate LG from HG tumors were similar, and both better than their respective means.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2018;48:491-498.

Accuracy of high b-value diffusion-weighted MRI for prostate cancer detection: a meta-analysis

Background The diagnostic accuracy of diffusion-weighted imaging (DWI) to detect prostate cancer is well-established. DWI provides visual as well as quantitative means of detecting tumor, the apparent diffusion coefficient (ADC). Recently higher b-values have been used to improve DWI's diagnostic performance. Purpose To determine the diagnostic performance of high b-value DWI at detecting prostate cancer and whether quantifying ADC improves accuracy. Material and Methods A comprehensive literature search of published and unpublished databases was performed. Eligible studies had histopathologically proven prostate cancer, DWI sequences using b-values ≥ 1000 s/mm², less than ten patients, and data for creating a 2 × 2 table. Study quality was assessed with QUADAS-2 (Quality Assessment of diagnostic Accuracy Studies). Sensitivity and specificity were calculated and tests for statistical heterogeneity and threshold effect performed. Results were plotted on a summary receiver operating characteristic curve (sROC) and the area under the curve (AUC) determined the diagnostic performance of high b-value DWI. Results Ten studies met eligibility criteria with 13 subsets of data available for analysis, including 522 patients. Pooled sensitivity and specificity were 0.59 (95% confidence interval [CI], 0.57-0.61) and 0.92 (95% CI, 0.91-0.92), respectively, and the sROC AUC was 0.92. Subgroup analysis showed a statistically significant ( P = 0.03) improvement in accuracy when using tumor visual assessment rather than ADC. Conclusion High b-value DWI gives good diagnostic performance for prostate cancer detection and visual assessment of tumor diffusion is significantly more accurate than ROI measurements of ADC.

Optimization of prostate MRI acquisition and post-processing protocol: a pictorial review with access to acquisition protocols

The aim of this review article is to provide insight into the optimization of 1.5-Testla (T) and 3-T prostate magnetic resonance imaging (MRI). An approach for optimization of data quantification, especially diffusion-weighted imaging (DWI), is provided. Benefits and limitations of various pulse sequences are discussed. Importable MRI protocols and access to imaging datasets is provided. Careful optimization of prostate MR acquisition protocol allows the acquisition of high-quality prostate MRI using clinical 1.5-T/3-T MR scanners with an overall acquisition time < 15 min.

Prospective evaluation of 18F-FACBC PET/CT and PET/MRI versus multiparametric MRI in intermediate- to high-risk prostate cancer patients (FLUCIPRO trial)

PURPOSE

The purpose of this study was to evaluate ¹⁸F-FACBC PET/CT, PET/MRI, and multiparametric MRI (mpMRI) in detection of primary prostate cancer (PCa).

METHODS

Twenty-six men with histologically confirmed PCa underwent PET/CT immediately after injection of 369 ± 10 MBq ¹⁸F-FACBC (fluciclovine) followed by PET/MRI started 55 ± 7 min from injection. Maximum standardized uptake values (SUV_max) were measured for both hybrid PET acquisitions. A separate mpMRI was acquired within a week of the PET scans. Logan plots were used to calculate volume of distribution (V_T). The presence of PCa was estimated in 12 regions with radical prostatectomy findings as ground truth. For each imaging modality, area under the curve (AUC) for detection of PCa was determined to predict diagnostic performance. The clinical trial registration number is NCT02002455.

RESULTS

In the visual analysis, 164/312 (53%) regions contained PCa, and 41 tumor foci were identified. PET/CT demonstrated the highest sensitivity at 87% while its specificity was low at 56%. The AUC of both PET/MRI and mpMRI significantly (p < 0.01) outperformed that of PET/CT while no differences were detected between PET/MRI and mpMRI. SUV_max and V_T of Gleason score (GS) >3 + 4 tumors were significantly (p < 0.05) higher than those for GS 3 + 3 and benign hyperplasia. A total of 442 lymph nodes were evaluable for staging, and PET/CT and PET/MRI demonstrated true-positive findings in only 1/7 patients with metastatic lymph nodes.

CONCLUSIONS

Quantitative ¹⁸F-FACBC imaging significantly correlated with GS but failed to outperform MRI in lesion detection. ¹⁸F-FACBC may assist in targeted biopsies in the setting of hybrid imaging with MRI.

Advanced Diffusion-weighted Imaging Modeling for Prostate Cancer Characterization: Correlation with Quantitative Histopathologic Tumor Tissue Composition-A Hypothesis-generating Study

Purpose To correlate quantitative diffusion-weighted imaging (DWI) parameters derived from conventional monoexponential DWI, stretched exponential DWI, diffusion kurtosis imaging (DKI), and diffusion-tensor imaging (DTI) with quantitative histopathologic tumor tissue composition in prostate cancer in a preliminary hypothesis-generating study. Materials and Methods This retrospective institutional review board-approved study included 24 patients with prostate cancer (mean age, 63 years) who underwent magnetic resonance (MR) imaging, including high-b-value DWI and DTI at 3.0 T, before prostatectomy. The following parameters were calculated in index tumors and nontumoral peripheral zone (PZ): apparent diffusion coefficient (ADC) obtained with monoexponential fit (ADC_ME), ADC obtained with stretched exponential modeling (ADC_SE), anomalous exponent (α) obtained at stretched exponential DWI, ADC obtained with DKI modeling (ADC_DKI), kurtosis with DKI, ADC obtained with DTI (ADC_DTI), and fractional anisotropy (FA) at DTI. Parameters in prostate cancer and PZ were compared by using paired Student t tests. Pearson correlations between tumor DWI and quantitative histologic parameters (nuclear, cytoplasmic, cellular, stromal, luminal fractions) were determined. Results All DWI parameters were significantly different between prostate cancer and PZ (P < .012). ADC_ME, ADC_SE, and ADC_DKI all showed significant negative correlation with cytoplasmic and cellular fractions (r = -0.546 to -0.435; P < .034) and positive correlation with stromal fractions (r = 0.619-0.669; P < .001). ADC_DTI and FA showed correlation only with stromal fraction (r = 0.512 and -0.413, respectively; P < .045). α did not correlate with histologic parameters, whereas kurtosis showed significant correlations with histopathologic parameters (r = 0.487, 0.485, -0.422 for cytoplasmic, cellular, and stromal fractions, respectively; P < .040). Conclusion Advanced DWI methods showed significant correlations with histopathologic tissue composition in prostate cancer. These findings should be validated in a larger study. ^© RSNA, 2017 Online supplemental material is available for this article. An earlier incorrect version of this article appeared online. This article was corrected on November 10, 2017.

Model selection for high b-value diffusion-weighted MRI of the prostate

PURPOSE

To assess the abilities of the standard mono-exponential (ME), bi-exponential (BE), diffusion kurtosis (DK) and stretched exponential (SE) models to characterize diffusion signal in malignant and prostatic tissues and determine which of the four models best characterizes these tissues on a per-voxel basis.

MATERIALS AND METHODS

This institutional-review-board-approved, HIPAA-compliant, retrospective study included 55 patients (median age, 61years; range, 42-77years) with untreated, biopsy-proven PCa who underwent endorectal coil MRI at 3-Tesla, diffusion-weighted MRI acquired at eight b-values from 0 to 2000s/mm². Estimated parameters were apparent diffusion coefficent (ME model); diffusion coefficients for the fast (D_fast) and slow (D_slow) components and fraction of fast component, f_fast (BE model); diffusion coefficient D, and kurtosis K (DK model); distributed diffusion coefficient DDC and α for (SE model). For one region-of-interest (ROI) in PZ and another in PCa in each patient, the corrected Akaike information criterion (AICc) and the Akaike weight (w) were calculated for each voxel.

RESULTS

Based on AICc and w, all non-monoexponential models outperformed the ME model in PZ and PCa. The DK model in PZ and SE model in PCa ROIs best fit the greatest average percentages of voxels (39% and 43%, respectively) and had the highest mean w (35±16×10^-2 and 41±22×10^-2, respectively).

CONCLUSION

DK and SE models best fit DWI data in PZ and PCa, and non-ME models consistently outperformed the ME model. Voxel-wise mapping of the preferential model demonstrated that the vast majority of voxels in either tissue type were best fit with one of the non-monoexponential models. At the given SNR levels, the maximum b-value of 2000s/mm² is not sufficiently high to identify the preferred non-monoexponential model.

Diffusion weighted imaging of the prostate-principles, application, and advances

This review article aims to provide an overview on the principles of diffusion-weighted magnetic resonance imaging (DW-MRI) and its applications in the imaging of the prostate. DW-MRI with regards to different applications for prostate cancer (PCa) detection and characterization, local staging as well as for active surveillance (AS) and tumor recurrence after radical prostatectomy (RP) will be discussed. Furthermore, advances in DW-MRI techniques like diffusion kurtosis imaging (DKI) will be presented.

Prostate Cancer: Diffusion-weighted MR Imaging for Detection and Assessment of Aggressiveness-Comparison between Conventional and Kurtosis Models

Purpose To compare standard diffusion-weighted (DW) imaging and diffusion kurtosis (DK) imaging for prostate cancer (PC) detection and characterization in a large patient cohort, with attention to the potential added value of DK imaging. Materials and Methods This retrospective institutional review board-approved study received a waiver of informed consent. Two hundred eighty-five patients with PC underwent 3.0-T phased-array coil prostate magnetic resonance (MR) imaging, including a DK imaging sequence (b values 0, 500, 1000, 1500, and 2000 sec/mm²) before prostatectomy. Maps of apparent diffusion coefficient (ADC) and diffusional kurtosis (K) were derived by using maximal b values of 1000 and 2000 sec/mm², respectively. Mean ADC and K were obtained from volumes of interest (VOIs) placed on each patient's dominant tumor and benign prostate tissue. Metrics were compared between benign and malignant tissue, between Gleason score (GS) ≤ 3 + 3 and GS ≥ 3 + 4 tumors, and between GS ≤ 3 + 4 and GS ≥ 4 + 3 tumors by using paired t tests, analysis of variance, receiver operating characteristic (ROC) analysis, and exact tests. Results ADC and K showed significant differences for benign versus tumor tissues, GS ≤ 3 + 3 versus GS ≥ 3 + 4 tumors, and GS ≤ 3 + 4 versus GS ≥ 4 + 3 tumors (P < .001 for all). ADC and K were highly correlated (r = -0.82; P < .001). Area under the ROC curve was significantly higher (P = .002) for ADC (0.921) than for K (0.902) for benign versus malignant tissue but was similar for GS ≤ 3 + 3 versus GS ≥ 3 + 4 tumors (0.715-0.744) and GS ≤ 3 + 4 versus GS ≥ 4 + 3 tumors (0.694-0.720) (P > .15). ADC and K were concordant for these various outcomes in 80.0%-88.6% of patients; among patients with discordant results, ADC showed better performance than K for GS ≤ 3 + 4 versus GS ≥ 4 + 3 tumors (P = .016) and was similar to K for other outcomes (P > .136). Conclusion ADC and K were highly correlated, had similar diagnostic performance, and were concordant for the various outcomes in the large majority of cases. These observations did not show a clear added value of DK imaging compared with standard DW imaging for clinical PC evaluation. ^© RSNA, 2017 Online supplemental material is available for this article.

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

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.

Molecular imaging and fusion targeted biopsy of the prostate

PURPOSE

This paper provides a review on molecular imaging with positron emission tomography (PET) and magnetic resonance imaging (MRI) for prostate cancer detection and its applications in fusion targeted biopsy of the prostate.

METHODS

Literature search was performed through the PubMed database using the keywords "prostate cancer", "MRI/ultrasound fusion", "molecular imaging", and "targeted biopsy". Estimates in autopsy studies indicate that 50% of men older than 50 years of age have prostate cancer. Systematic transrectal ultrasound (TRUS) guided prostate biopsy is considered the standard method for prostate cancer detection and has a significant sampling error and a low sensitivity. Molecular imaging technology and new biopsy approaches are emerging to improve the detection of prostate cancer.

RESULTS

Molecular imaging with PET and MRI shows promising results in the early detection of prostate cancer. MRI/TRUS fusion targeted biopsy has become a new clinical standard for the diagnosis of prostate cancer. PET molecular image-directed, three-dimensional ultrasound-guided biopsy is a new technology that has great potential for improving prostate cancer detection rate and for distinguishing aggressive prostate cancer from indolent disease.

CONCLUSION

Molecular imaging and fusion targeted biopsy are active research areas in prostate cancer research.

Computer-aided Detection of Prostate Cancer with MRI: Technology and Applications

One in six men will develop prostate cancer in his lifetime. Early detection and accurate diagnosis of the disease can improve cancer survival and reduce treatment costs. Recently, imaging of prostate cancer has greatly advanced since the introduction of multiparametric magnetic resonance imaging (mp-MRI). Mp-MRI consists of T2-weighted sequences combined with functional sequences including dynamic contrast-enhanced MRI, diffusion-weighted MRI, and magnetic resonance spectroscopy imaging. Because of the big data and variations in imaging sequences, detection can be affected by multiple factors such as observer variability and visibility and complexity of the lesions. To improve quantitative assessment of the disease, various computer-aided detection systems have been designed to help radiologists in their clinical practice. This review paper presents an overview of literatures on computer-aided detection of prostate cancer with mp-MRI, which include the technology and its applications. The aim of the survey is threefold: an introduction for those new to the field, an overview for those working in the field, and a reference for those searching for literature on a specific application.

Evaluating Prostate Cancer Using Fractional Tissue Composition of Radical Prostatectomy Specimens and Pre-Operative Diffusional Kurtosis Magnetic Resonance Imaging

BACKGROUND

Evaluating tissue heterogeneity using non-invasive imaging could potentially improve prostate cancer assessment and treatment.

METHODS

20 patients with intermediate/high-risk prostate cancer underwent diffusion kurtosis imaging, including calculation of apparent diffusion (Dapp) and kurtosis (Kapp), prior to radical prostatectomy. Whole-mount tissue composition was quantified into: cellularity, luminal space, and fibromuscular stroma. Peripheral zone tumors were subdivided according to Gleason score.

RESULTS

Peripheral zone tumors had increased cellularity (p<0.0001), decreased fibromuscular stroma (p<0.05) and decreased luminal space (p<0.0001). Gleason score ≥4+3 tumors had significantly increased cellularity and decreased fibromuscular stroma compared to Gleason score ≤3+4 (p<0.05). In tumors, there was a significant positive correlation between median Kapp and cellularity (ρ = 0.50; p<0.05), and a negative correlation with fibromuscular stroma (ρ = -0.45; p<0.05). In normal tissue, median Dapp had a significant positive correlation with luminal space (ρ = 0.65; p<0.05) and a negative correlation with cellularity (ρ = -0.49; p<0.05). Median Kapp and Dapp varied significantly between tumor and normal tissue (p<0.0001), but only median Kapp was significantly different between Gleason score ≥4+3 and ≤3+4 (p<0.05).

CONCLUSIONS

Peripheral zone tumors have increased cellular heterogeneity which is reflected in mean Kapp, while normal prostate has a more homogeneous luminal space and cellularity better represented by Dapp.

Fitting methods for intravoxel incoherent motion imaging of prostate cancer on region of interest level: Repeatability and gleason score prediction

PURPOSE

To evaluate different fitting methods for intravoxel incoherent motion (IVIM) imaging of prostate cancer in the terms of repeatability and Gleason score prediction.

METHODS

Eighty-one patients with histologically confirmed prostate cancer underwent two repeated 3 Tesla diffusion-weighted imaging (DWI) examinations performed using 14 b-values in the range of 0-500 s/mm² and diffusion time of 19.004 ms. Mean signal intensities of regions-of-interest were fitted using five different fitting methods for IVIM as well as monoexponential, kurtosis, and stretched exponential models. The fitting methods and models were evaluated in the terms of fitting quality [Akaike information criteria (AIC)], repeatability, and Gleason score prediction. Tumors were classified into three groups (3 + 3, 3 + 4, > 3 + 4). Machine learning algorithms were used to evaluate the performance of the combined use of the parameters. Simulation studies were performed to evaluate robustness of the fitting methods against noise.

RESULTS

Monoexponential model was preferred over IVIM based on AIC. The "pseudodiffusion" parameters demonstrated low repeatability and clinical value. Median "pseudodiffusion" fraction values were below 8.00%. Combined use of the parameters did not outperform the monoexponential model.

CONCLUSION

Monoexponential model demonstrated the highest repeatability and clinical values in the regions-of-interest based analysis of prostate cancer DWI, b-values in the range of 0-500 s/mm² . Magn Reson Med 77:1249-1264, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Advanced Noncontrast MR Imaging in Musculoskeletal Radiology

Multiple nonmorphologic magnetic resonance sequences are available in musculoskeletal imaging that can provide additional information to better characterize and diagnose musculoskeletal disorders and diseases. These sequences include blood-oxygen-level-dependent (BOLD), arterial spin labeling (ASL), diffusion-weighted imaging (DWI), and diffusion-tensor imaging (DTI). BOLD and ASL provide different methods to evaluate skeletal muscle microperfusion. The BOLD signal reflects the ratio between oxyhemoglobin and deoxyhemoglobin. ASL uses selective tagging of inflowing blood spins in a specific region for calculating local perfusion. DWI and DTI provide information about the structural integrity of soft tissue including muscles and fibers as well as pathologies.

Relaxation along fictitious field, diffusion-weighted imaging, and T2 mapping of prostate cancer: Prediction of cancer aggressiveness

PURPOSE

To evaluate the performance of relaxation along a fictitious field (RAFF) relaxation time (TRAFF ), diffusion-weighted imaging (DWI)-derived parameters, and T2 relaxation time values for prostate cancer (PCa) detection and characterization.

METHODS

Fifty patients underwent 3T MR examination using surface array coils before prostatectomy. DWI was performed using 14 and 12 b values in the ranges of 0-500 s/mm(2) and 0-2000 s/mm(2) , respectively. Repeated MR examination was performed in 16 patients. TRAFF , DWI-derived parameters (monoexponential, kurtosis, biexponential models), and T2 values were measured and averaged over regions of interest placed in PCa and normal tissue. Repeatability of TRAFF and DWI-derived parameters were assessed by coefficient of repeatability and intraclass correlation coefficient ICC(3,1). Areas under the receiver operating characteristic curve (AUCs) for PCa detection and Gleason score classification were estimated. The parameters were correlated with Gleason score groups using Spearman correlation coefficient (ρ).

RESULTS

ICC(3,1) values for TRAFF were in the range of 0.82-0.92. TRAFF values had higher AUC values for Gleason score classification compared with DWI-derived parameters and T2 . The RAFF method demonstrated the highest ρ value (-0.65).

CONCLUSION

In a quantitative region of interest-based analysis, RAFF outperformed DWI ("low" and "high" b values) and T2 mapping in the characterization of PCa.