Multiparametric magnetic resonance imaging (mpMRI) is the essential imaging tool in prostate cancer (PCa) diagnosis, with increasing prominence as recent guidelines suggest using mpMRI before systemic biopsy in men with a suspicion of PCa [1]. The primary treatment of PCa is radical prostatectomy, while eligible patients might opt for active surveillance or focal therapy [2], avoiding comorbidities related to radical prostatectomy [3]. The precise estimation of PCa index lesion volume on mpMRI is essential in guiding focal therapy since the overestimation of the lesion size increases the risk of complications, while the underestimation might lead to insufficient disease control.
To date, several studies have explored the accuracy of prostate mpMRI in measuring the PCa index lesion volume with contradictory results [2, 4, 5, 6, 7]. Some of the studies had flaws, including using the ellipsoid formula for determining the ground truth volume [4, 8], using biopsy specimen volume as the ground truth [9], overlooking the individual performance of mpMRI sequences [10], and relying on the correlation without assessing actual volumetric differences [11, 12]. Furthermore, there is little evidence for the impact of an endorectal coil (ERC) in estimating the index lesion volume [8].
In this work, we compared the effectiveness of prostate mpMRI sequences—T2-weighted imaging (T2W), diffusion-weighted imaging (DWI) and apparent diffusion coefficient (ADC) maps, and dynamic contrast-enhanced imaging (DCE)— in assessing the index lesion volume using whole-mount pathology as the ground-truth. Additionally, we evaluated whether using an ERC impacts the performance.
The local ethics committee approved this study and waived the need for informed consent for the retrospective evaluation of anonymized medical data. We reviewed medical records of consecutive PCa patients who underwent prostatectomy in our institution between January 2019 and December 2020. The patients meeting the following inclusion criteria were identified: (1) diagnosis of PCa on whole-mount pathology; (2) having a digitized whole-mount specimen; (3) having a prostate mpMRI scan obtained at 3T mpMRI at our center three months before the radical prostatectomy. Initially, 87 patients were identified based on the inclusion criteria.
We excluded patients based on the following criteria: (1) patients with prominent post-biopsy hemorrhage obscuring the index lesion (n = 2); (3) patients with index lesion volume less than 0.5 mL on whole-mount pathology (n = 6). Furthermore, we excluded 3 PCa patients with a Gleason Score of 3 + 4 without any intraductal pattern or cribriform component as the index lesion was invisible on mpMRI.
All patients underwent prostate mpMRI on a 3.0 Tesla MRI scanner (Skyra, Siemens Medical Systems, Germany); the mpMRI scans were performed with an 18-channel phased-array surface coil and a liquid perfluorocarbon-filled ERC (Medrad, Bayer). All mpMRI scans were acquired at least six weeks after the biopsy to minimize the negative influence of post-biopsy hemorrhage on the diagnostic evaluation. The prostate mpMRI protocol comprised tri-planar T2-weighted imaging, DWI, and DCE imaging from the first to the last. The detailed parameters regarding the MRI sequences are given in Supplementary Table S1.
A radiologist with over 20 years of prostate imaging experience measured the volumes of the index lesion on mpMRI. To ensure that the lesions being measured were the same lesions on whole-mount pathology and mpMRI, we informed the radiologists regarding the location of the index lesion on whole-mount pathology using the sector map proposed by the PI-RADS [13]. The radiologist was blinded to the histopathological volume, yet we allowed the radiologist to evaluate all mpMRI sequences to simulate a clinical workflow following prior studies [5]. The radiologist was free to set the window level at their preference for each patient and sequence.
First, the radiologist manually drew the borders of the index lesion on T2W, DWI (b-value of 2000 s/mm2), ADC maps, and KTrans maps on a dedicated workstation (Syngo Via, Siemens Medical Systems, Germany). Then, the software automatically calculated the index lesion volumes. These volumes were referred to as T2Wvol, DWIvol, ADCvol, and DCEvol. Additionally, we determined the maximum index lesion volume from any of mpMRI sequences (i.e., T2Wvol, DWIvol, ADCvol, and DCEvol) for each patient (maxMRIvol) [7]. For instance, for an individual patient, if the largest lesion size was measured on T2Wvol, then that volume was selected in calculating maxMRIvol.
A genitourinary pathologist evaluated all whole-mount specimens following the guidelines. The pathologist first assessed the macroscopic specimens visually, then 3 mm tissue blocks were created using a purpose-built mega cassette. The tissue blocks were embedded in paraffin, sliced at a thickness of 3 um per block (i.e., a single slice of 3 um was obtained per 3 mm block) perpendicular to the posterior plane, and stained with hematoxylin and eosin. Prepared micro-slices were mounted on a purpose-built mega slide, digitalized, and transferred to in-house software.
The genitourinary pathologist measured the index lesion using planimetry. The software automatically calculated the volume of the index lesion (whole-mountvol). We did not apply a correction factor for tissue shrinkage following prior work [7]. The lesion with the highest Gleason score was defined as the index lesion; if multiple lesions were present, the lesion with the largest volume was accepted as the index lesion [14].
The statistical analyses were performed using the SciPy library of the Python programming language. The normally distributed continuous variables are presented using the mean and standard deviations. The continuous variables without a normal distribution are presented with the median and interquartile range (IQR), and the categorical and ordinal variables are presented with the frequencies and percentages. Next, the Spearman’s Rank-Order correlation was run to investigate the correlations between the histopathological and radiological volumes. Finally, we ran the Friedman and post-hoc Wilcoxon signed-rank tests to compare the index lesion volumes on whole-mount pathology and mpMRI.
The relative volumetric differences amongst mpMRI sequences using whole-mount pathology as the ground truth were assessed using the Friedman and post-hoc Wilcoxon signed-rank tests. Finally, the patients were divided respecting the use of an ERC, and the Man-Whitney-U test was run to compare the relative maxMRIvol differences compared with whole-mount pathology. A P value < 0.05 was accepted as showing a significant result.
In all, 72 men with a mean age of 63.36 ± 11.25 years were enrolled in the study. The median whole-mountvol was 1.97 cm3 (IQR, 3.35), while the median T2Wvol, DWIvol, ADCvol, and DCEvol were 0.68 cm3 (IQR, 1.15), 0.97 cm3 (IQR, 2.21), 0.98 cm3 (IQR, 2.21), 0.82 cm3 (IQR, 1.15) respectively. The Spearman test provided good positive correlations between whole-mountvol and T2Wvol, DWIvol, ADCvol, and DCEvol, with the r values of 0.75, 0.77, 0.78, and 0.75, respectively. The detailed characteristics of the study sample are given in Table 1.
Table 1
The characteristics of the study sample.
VARIABLES | FINDINGS |
---|---|
Age (years) | 63 ± 11.25 |
PI-RADS score | |
PI-RADS 4 | 31 (43%) |
PI-RADS 5 | 41 (57%) |
Lesion zone | |
Peripheral | 60 (83%) |
Transitional | 12 (17%) |
Gleason score | |
3+4 | 44 (61%) |
4+3 | 19 (26%) |
4+4 | 1 (1%) |
4+5 | 8 (11%) |
Prostate volume on MRI (cm3) | 37.26 ± 16.34 |
Prostate volume on pathology (ml) | 45 ± 15 |
PSA level (ng/ml) | 7.08 ± 5.02 |
PSA density (ng/ml2) | 0.14 ± 0.12 |
Endorectal coil | 39 (54%) |
PI-RADS: Prostate Imaging Reporting & Data System; MRI: Magnetic Resonance Imaging; PSA: Prostate Specific Antigen.
The Friedman test demonstrated that mpMRI-derived volumes were significantly lower than the whole-mountvol (P = 0.001) (Table 2). The intra-patient differences of T2Wvol, DWIvol, ADCvol, and DCEvol with respect to whole-mountvol are shown with the Bland-Altman plots in Figure 1. Figure 2 depicts index lesion measurements on mpMRI and digitized whole-mount pathology in a PCa patient.
Table 2
The index lesion volumes on whole-mount pathology and multi-parametric MRI.
VARIABLES | ABSOLUTE INDEX LESION VOLUME (CM3)* | RELATIVE VOLUME DIFFERENCE (%) ** | P VALUE*** |
---|---|---|---|
Whole-mount pathology | 1.97 (3.35) | NA | NA |
T2-weighted imaging | 0.68 (1.15) | 54.32% ± 32.84%(–67.27% – +98.48%) | <0.0001 |
Diffusion-weighted imaging | 0.97 (2.21) | 35.30% ± 41.21%(–105.08% – 95.96%) | <0.0001 |
Apparent diffusion coefficient map | 0.98 (2.21) | 35.96% ± 40.42%(–121.82% – +95.25%) | <0.0001 |
Dynamic contrast-enhanced imaging | 0.82 (1.15) | 41.98 ± 39.83%(–72.73% – 95.76%) | <0.0001 |
Maximum volume on MRI | 1.13 (2.35) | 24.72 ± 45.09%(–121.82% – 95.25%) | <0.0001 |
* Presented with a median and interquartile range; ** Presented with a mean, standard deviation, and range; *** Reflect the comparisons between the index lesion volumes on whole-mount pathology and respective mpMRI sequence.
Bland-Altman plots comparing T2-weighted imaging (T2Wvol), diffusion-weighted imaging (DWIvol), apparent diffusion coefficient (ADCvol), and dynamic contrast-enhanced (DCEvol) imaging in measuring prostate cancer index lesion volume compared with whole-mount pathology. All mpMRI sequences had a significant volumetric underestimation (P < 0.0001).
A 62-year-old man with Gleason Score 3 + 4 prostate cancer. The measurements of the index lesion volume on T2-weighted imaging (a), apparent diffusion coefficient map (b), diffusion-weighted imaging with a b-value of 2000 s/mm2(c), and the grayscale KTrans map of dynamic contrast-enhanced imaging (d) are shown. In addition, the reference digitized whole-mount pathology slice is depicted (e).
The mean relative index lesion volume underestimation for T2Wvol, DWIvol, ADCvol, and DCEvol compared with whole-mountvol were 54.32% ± 32.83% (range, –67.27 – +98.48), 35.3% ± 41.21% (range, –105.08 – +95.06), 35.96% ± 40.41% (range, –121.82% – +95.25%), 41.98% ± 39.82% (range, –72.73 – +95.76).
The Friedman test showed a significant difference in the relative index lesion volume underestimation amongst sequences (P = 0.001). DWIvol and ADCvol had a lower mean relative index lesion volume underestimation than T2Wvol (P = 0.001) and DCEvol (P = 0.001 and 0.005). Additionally, DCEvol demonstrated a lower relative volume underestimation than the T2Wvol (P = 0.001).
The mean relative index lesion volume underestimations of maxMRIvol were 39.16% ± 32.58% (range, –63.38% – +81.53%) and 7.65% ± 51.91% (range, –121.81% – +95.25%) with and without an ERC, respectively (P = 0.002).
The main findings of the study are as follows: (1) T2Wvol, DWIvol, ADCvol, and DCEvol show a good correlation with whole-mount pathology in measuring the index lesion volume; (2) all mpMRI sequences underestimate the index lesion volume compared with whole-mount pathology; (3) T2Wvol leads to the largest index lesion volume underestimation, while DWIvol and ADCvol have the smallest; (4) notably, mpMRI with an ERC leads to a more significant index lesion volume underestimation.
The findings of the present work contrast with the initial studies on the performance of prostate mpMRI in measuring the index lesion volume, which documented that mpMRI overestimated the index lesion volume [11, 15]. Nevertheless, these studies were carried out before the PI-RADS era and had flaws, as mentioned in the introduction. Further, prostate MRI have improved dramatically over the last decade. Indeed, recent studies document that mpMRI underestimates the index lesion volume [16, 17]. Nevertheless, relatively few studies explicitly investigated which mpMRI sequence is the most effective in assessing PCa index lesion volume.
Cornud et al. [17] investigated the performance of the individual mpMRI sequences in assessing PCa index lesion volume in patients who underwent prostate mpMRI at 1.5T with an ERC. Similar to our results, the authors observed that ADC had the best performance in estimating PCa index lesion volume. DCE provided a worse performance than T2-weighted imaging. The authors argue that the poor performance of DCE might relate to the artifacts due to post-biopsy hemorrhage. In addition, we argue that field strength differences between studies might also contribute to these discrepancies.
Bratan et al. [7] and Sun et al. [18] showed that all mpMRI sequences underestimated PCa volumes, with maxMRIvol providing the most accurate estimations. However, Bratan et al. and Sun et al. showed that DCE provided a better index lesion volume estimation than T2-weighted imaging and ADC. Unlike the present work, Sun et al. [18] measured the maximum enhancement area instead of kinetic parametric maps, probably leading to differences between the studies.
Le Nobin et al. [8] also showed that mpMRI significantly underestimated PCa volumes, with T2-weighted imaging providing a better volume estimate than ADC, yet contrasting with our findings. The authors argued that including patients with Gleason Score 3 + 3 PCa in their sample might be the underlying factor. Furthermore, they measured PCa foci other than the index lesion, which we assume might explain the inconsistency.
A few investigators compared the performance of mpMRI and prostate-specific membrane antigen positron emission tomography-computed tomography (PSMA-PET/CT) in estimating PCa volume [19, 21]. For instance, Bettermann et al. [19] found that mpMRI significantly underestimated the index lesion volume, whereas PSMA-PET/CT provided a better estimation compared with whole-mount pathology. Likewise, Zamboglou et al. [20] also found that mpMRI underestimated the index lesion volume on whole-mount pathology. The authors demonstrated that the fusion of mpMRI and PSMA-PET/CT provides a better estimate compared with both modalities alone.
It is recognized that an ERC might deform the prostate gland, which might lead to differences in index lesion volume. Nevertheless, few earlier studies explicitly investigated this disadvantage. In accordance with our results, Bratan et al. [7] found that mpMRI at 3T with an ERC led to reduced performance in assessing PCa volume. Though it might be expected due to potential gland deformation of an ERC, the low performance of mpMRI with an ERC compared without an ERC in assessing PCa index lesion volume was a noteworthy finding.
Several limitations to the present work should be acknowledged. First and foremost, this was a single-center study with all mpMRI scans obtained from a single 3T MRI scanner, which hinders the generalizability. Second, we had a relatively low sample size, preventing us from investigating whether mpMRI-derived index lesion volumes vary across PI-RADS or Gleason scores or lesion location.
Third, we did not apply any correction factor for tissue shrinkage in the present work. Prior studies implemented a broad range of correction factors changing from 1 to 1.5 for compensating tissue shrinkage on whole-mount pathology [10, 11, 21, 22, 23]. On the other hand, in line with our work, some other authors argued that it was not possible to find an optimal correction factor for tissue shrinkage as it varies from patient to patient [7, 17]. Nevertheless, mpMRI-based volumes underestimated the index lesion volume on whole-mount pathology without any correction factor in the present work. Thus, applying a shrinkage factor would only have exacerbated the volume underestimation. Hence, we suggest that this limitation does not hamper the findings of the present work.
In conclusion, the findings of the present work highlight that mpMRI-derived volumes, T2Wvol, DWIvol, ADCvol, DCEvol, and maxMRIvol underestimate PCa index lesion volume compared with whole-mount pathology, with T2Wvol having the largest volume underestimation and maxMRIvol having the smallest. Therefore, using maxMRIvol appears reasonable in assessing PCa index lesion volume. Additionally, using an ERC exacerbates the volume underestimation, and extra care should be exercised while estimating PCa index lesion volume on mpMRI with an ERC.
The additional file for this article can be found as follows:
Table S1The detailed prostate multi-parametric magnetic resonance imaging parameters. DOI: https://doi.org/10.5334/jbsr.2832.s1
All procedures were performed in studies involving human participants following the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration.
The authors have no competing interests to declare.
Dasgupta P, Davis J, Hughes S. NICE guidelines on prostate cancer 2019. BJU Int. 2019; 124: 1–1. DOI: https://doi.org/10.1111/bju.14815
Oto A, Sethi I, Karczmar G, et al. MR Imaging–guided Focal Laser Ablation for Prostate Cancer: Phase I Trial. Radiology. 2013; 267: 932–40. DOI: https://doi.org/10.1148/radiol.13121652
Wilt TJ, MacDonald R, Rutks I, Shamliyan TA, Taylor BC, Kane RL. Systematic Review: Comparative Effectiveness and Harms of Treatments for Clinically Localized Prostate Cancer. Ann Intern Med. 2008; 148: 435. DOI: https://doi.org/10.7326/0003-4819-148-6-200803180-00209
Rud E, Klotz D, Rennesund K, et al. Detection of the index tumour and tumour volume in prostate cancer using T2-weighted and diffusion-weighted magnetic resonance imaging (MRI) alone. BJU Int. 2014; 114: E32–42. DOI: https://doi.org/10.1111/bju.12637
Colvin R, Walker D, Hafron J, et al. Which measurement method should be used for prostate volume for PI-RADS? A comparison of ellipsoid and segmentation methods. Clinical Imaging. 2021; 80: 454–458. DOI: https://doi.org/10.1016/j.clinimag.2021.09.003
Matsugasumi T, Baco E, Palmer S, et al. Prostate Cancer Volume Estimation by Combining Magnetic Resonance Imaging and Targeted Biopsy Proven Cancer Core Length: Correlation with Cancer Volume. J Urol. 2015; 194: 957–65. DOI: https://doi.org/10.1016/j.juro.2015.04.075
Bratan F, Melodelima C, Souchon R, et al. How Accurate Is Multiparametric MR Imaging in Evaluation of Prostate Cancer Volume? Radiology. 2015; 275: 144–54. DOI: https://doi.org/10.1148/radiol.14140524
le Nobin J, Orczyk C, Deng F-M, et al. Prostate tumour volumes: evaluation of the agreement between magnetic resonance imaging and histology using novel co-registration software. BJU Int; 2014. DOI: https://doi.org/10.1111/bju.12750
Liechti MR, Muehlematter UJ, Schneider AF, et al. Manual prostate cancer segmentation in MRI: interreader agreement and volumetric correlation with transperineal template core needle biopsy. Eur Radiol. 2020; 30: 4806–15. DOI: https://doi.org/10.1007/s00330-020-06786-w
Puech P, Potiron E, Lemaitre L, et al. Dynamic contrast-enhanced–magnetic resonance imaging evaluation of intraprostatic prostate cancer: correlation with radical prostatectomy specimens. Urology. 2009; 74(5): 1094–1099. DOI: https://doi.org/10.1016/j.urology.2009.04.102
Mazaheri Y, Hricak H, Fine SW, et al. Prostate tumor volume measurement with combined T2-weighted imaging and diffusion-weighted MR: correlation with pathologic tumor volume. Radiology. 2009; 252: 449. DOI: https://doi.org/10.1148/radiol.2523081423
Isebaert S, Van den Bergh L, Haustermans K, et al. Multiparametric MRI for prostate cancer localization in correlation to whole-mount histopathology. J Magn Reson Imaging. 2013; 37: 1392–1401. DOI: https://doi.org/10.1002/jmri.23938
Weinreb JC, Barentsz JO, Choyke PL, et al. PI-RADS Prostate Imaging – Reporting and Data System: 2015, Version 2. Eur Urol. 2016; 69: 16–40. DOI: https://doi.org/10.1016/j.eururo.2015.08.052
Epstein JI, Egevad L, Amin MB, Delahunt B, Srigley JR, Humphrey PA. The 2014 International Society of Urological Pathology (ISUP) Consensus Conference on Gleason Grading of Prostatic Carcinoma. Am J Surg Pathol. 2016; 40: 244–52. DOI: https://doi.org/10.1097/PAS.0000000000000530
Coakley FV, Kurhanewicz J, Lu Y, et al. Prostate Cancer Tumor Volume: Measurement with Endorectal MR and MR Spectroscopic Imaging. Radiology. 2002; 223: 91–7. DOI: https://doi.org/10.1148/radiol.2231010575
Turkbey B, Mani H, Aras O, et al. Correlation of Magnetic Resonance Imaging Tumor Volume with Histopathology. J Urol. 2012; 188: 1157–63. DOI: https://doi.org/10.1016/j.juro.2012.06.011
Cornud F, Khoury G, Bouazza N, et al. Tumor Target Volume for Focal Therapy of Prostate Cancer—Does Multiparametric Magnetic Resonance Imaging Allow for a Reliable Estimation? J Urol. 2014; 191: 1272–9. DOI: https://doi.org/10.1016/j.juro.2013.12.006
Sun C, Chatterjee A, Yousuf A, et al. Comparison of T2-Weighted Imaging, DWI, and Dynamic Contrast-Enhanced MRI for Calculation of Prostate Cancer Index Lesion Volume: Correlation with Whole-Mount Pathology. AJR Am J Roentgenol. 2019; 212: 351–6. DOI: https://doi.org/10.2214/AJR.18.20147
Bettermann AS, Zamboglou C, Kiefer S, et al. [68Ga-] PSMA-11 PET/CT and multiparametric MRI for gross tumor volume delineation in a slice by slice analysis with whole mount histopathology as a reference standard – Implications for focal radiotherapy planning in primary prostate cancer. Radiother Oncol. 2019; 141: 214–219. DOI: https://doi.org/10.1016/j.radonc.2019.07.005
Zamboglou C, Drendel V, Jilg CA, et al. Comparison of 68Ga-HBED-CC PSMA-PET/CT and multiparametric MRI for gross tumour volume detection in patients with primary prostate cancer based on slice by slice comparison with histopathology. Theranostics. 2017; 7: 228–237. DOI: https://doi.org/10.7150/thno.16638
Jonmarker S, Valdman A, Lindberg A, et al. Tissue shrinkage after fixation with formalin injection of prostatectomy specimens. Virchows Arch. 2006; 449: 297. DOI: https://doi.org/10.1007/s00428-006-0259-5
Schned AR, Wheeler KJ, Hodorowski CA, et al. Tissue-shrinkage correction factor in the calculation of prostate cancer volume. Am J Surg Pathol. 1996; 20: 1501. DOI: https://doi.org/10.1097/00000478-199612000-00009
Stamey TA, McNeal JE, Freiha FS, et al. Morphometric and clinical studies on 68 consecutive radical prostatectomies. J Uro. 1988; 139: 1235. DOI: https://doi.org/10.1016/S0022-5347(17)42876-X