Luminal water imaging: a novel MRI method for prostate cancer diagnosis
Shirin Sabouri1, Silvia Chang2,3, Richard Savdie4, Jing Zhang5, Edward Jones6, Larry Goldenberg4,7, and Piotr Kozlowski2,4,5,7

1Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada, 2Radiology, University of British Columbia, Vancouver, BC, Canada, 3Vancouver General Hospital, Vancouver, BC, Canada, 4Urologic Sciences, University of British Columbia, Vancouver, BC, Canada, 5UBC MRI Research Center, Vancouver, BC, Canada, 6Pathology, Vancouver General Hospital, Vancouver, BC, Canada, 7Vancouver Prostate Centre, Vancouver, BC, Canada

Synopsis

MR multi-exponential T2 mapping can be used for extracting valuable information about tissue composition in prostate. Using this technique, the fractional volume of the luminal space in the prostatic tissue can be determined. Since tissue composition and the amounts of lumen differ between normal and cancerous tissues, this technique can be applied for detection of prostatic tumors. We have investigated the suitability of using MR multi-exponential T2 mapping for detection and staging of prostate cancer. We have acquired and analyzed MR images of 11 patients, and concluded that this technique is highly sensitive and specific in detection of prostatic tumors.

Target Audience

Clinicians and researchers interested in detection and staging of prostate cancer.

Purpose

To introduce an MRI technique for detection and grading of prostate tumors, and characterization of prostatic tissue.

Introduction

MR multi-exponential T2 mapping is a well-known imaging technique that has been applied for tracking histo-pathological changes in organs such as brain1. Evidence of multi-exponential T2 decay from prostatic tissue has been presented before2, but the suitability of this technique for diagnosis of prostate cancer has not been studied in details. Using multi-exponential T2 mapping, the fractional volume of the luminal space, or so called luminal water fraction (LWF), in the prostatic tissue can be determined. Because of the difference of tissue composition and lumen percentage in normal and cancerous prostatic tissues, we hypothesized that MR multi-exponential T2 mapping can be used for the detection and staging of prostate cancer. We have acquired and analyzed MR images of 11 patients, and investigated the accuracy of this technique in detection and evaluation of prostatic tumors by performing multi-parametric statistical analysis.

Methods

Data acquisition was carried out at the University of British Columbia (UBC) MRI Research Centre, using a 3.0T whole body MR scanner [Achieva 3.0T, Philips Medical Systems, Best, The Netherlands]. Eleven patients with biopsy proven cancer underwent MRI scan with an endorectal coil, prior to undergoing prostatectomy. A 3D multi-echo spin echo sequence (TR/TE=3061/25ms, NE=64, FOV=240x240x40mm3, voxel-size=1x1x4mm3, matrix-size=240x240) was used for scanning of the entire prostate gland. Images were analyzed with Matlab [The MathWorks Inc,Natick, MA, USA]. The analysis involved regularized Non-Negative Least Squares (NNLS)3,4 fitting of multi-exponential decay curves, which generated T2 distributions for every pixel (see Figure 1). The following parameters were defined and used to describe the T2 distributions: number of distinguishable T2 components (N) determined by counting the number of peaks in the distribution; geometric mean of the short (T2-short) and long (T2-long) components, as well as the geometric mean of the entire distribution (gmT2); ratio of area under the long component over the total area under the entire distribution (Luminal Water Fraction – LWF); and areas under the short (A1) and long (A2) components. Average values of these parameters were calculated within 255 ROIs manually outlined on digitized images of the whole-mount histology sections registered to MRI images5. Selection of ROIs included: cancerous peripheral (PZ) and transition (TZ) zones, and normal PZ, TZ, Anterior Fibromuscular Stroma (AFMS), and Periurethral Fibromuscular stroma (PFMS). Statistical analyses were performed with MedCalc [MedCalc Software, Mariakerke, Belgium]. Significant differences between MR parameters of tumor and normal tissues were determined with Kruskal-Wallis test. Correlations between MR parameters and Gleason score were evaluated in PZ and TZ with Spearman’s rank correlation test. Sensitivity and specificity were calculated by Receiver Operating Characteristic (ROC) analysis of individual and combined MR parameters.

Results and Discussions

Representative MR parametric maps are shown in Figure 2. Kruskal-Wallis test indicated that the average values of T2-short, gmT2, A1, A2, and LWF were significantly different between tumor and normal tissue in PZ and TZ, suggesting that these five parameters can be used as tumor indicator in prostate. The average values of N were significantly different between the glandular tissue (i.e. PZ and TZ) and non-glandular tissue (i.e. AFMS and PFMS); hence, N can be used for tissue classification and machine learning purposes. The highest correlation with Gleason score was obtained for LWF (-0.734, p <0.001 in PZ, and -0.712, p <0.001 in TZ). The values of sensitivity, specificity, and Area Under the ROC Curve (AUC) demonstrated high accuracy of tumour detection with the proposed technique (Table 1). The results of this pilot study demonstrate the suitability of MR multi-exponential T2 mapping for prostate cancer diagnosis. Several parameters were defined to characterize T2 decay curve of prostatic tissue, and it was shown that the defined parameters were successfully applied for prostate tissue classification, detection of its cancerous tissue with high sensitivity and specificity, and grading of prostate cancer.

Acknowledgements

This study has been supported by the Canadian Institutes of Health Research. We thank Margaret Luk, Laura Barlow, and Alex Mazur for their kind support and assistance in this research.

References

[1] MacKay A, et al. Magn Reson Med 1994;31:673–677. [2] Storås T. H., et al. J. Magn. Reson. Imaging. 2008;28:1166–1172. [3] Bjarnason TA., Mitchell JR. J Magn Reson 2010;206:200–4. [4] Prasloski T., et al. Magn Reson Med. 2012;67(6):1803-14. [5] Uribe CF et al., Magn. Reson. Imaging, 2015;33:577-583.

Figures

NNLS fitting of a multi-exponential T2 signal decay curve (left-hand side), and its corresponding T2 distribution (right-hand side). This data is from a single pixel in the peripheral zone.

Representative maps of MR parameters and histology wholemount slide of the same slice. Scale bar of gmT2, T2-short , and T2-long images are in ‘ms’. Zero pixels on the T2long map indicate mono-exponential decay.

The values of AUC, sensitivity, and specificity calculated from the variable (s) which provides the max AUC for each case.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2489