Identification of prostate cancer using multiparametric MR imaging characteristics of prostate tissues referenced to whole mount histopathology

PSA, an outdated biomarker for prostate cancer: In search of a more specific biomarker, citrate takes the spotlight

The prevailing biomarker employed for prostate cancer (PCa) screening and diagnosis is the prostate-specific antigen (PSA). Despite excellent sensitivity, PSA lacks specificity, leading to false positives, unnecessary biopsies and overdiagnosis. Consequently, PSA is increasingly less used by clinicians, thus underscoring the imperative for the identification of new biomarkers. An emerging biomarker in this context is citrate, a molecule secreted by the normal prostate, which has been shown to be inversely correlated with PCa. Here, we discuss about PSA and its usage for PCa diagnosis, its lack of specificity, and the various conditions that can affect its levels. We then provide our vision about what we think would be a valuable addition to our PCa diagnosis toolkit, citrate. We describe the unique citrate metabolic program in the prostate and how this profile is reprogrammed during carcinogenesis. Finally, we summarize the evidence that supports the usage of citrate as a biomarker for PCa diagnosis, as it can be measured in various patient samples and be analyzed by several methods. The unique relationship between citrate and PCa, combined with the stability of citrate levels in other prostate-related conditions and the simplicity of its detection, further accentuates its potential as a biomarker.

The impact of benign tissue within cancerous regions in the prostate: Characterizing sparse and dense prostate cancers on whole-mount histopathology and on multiparametric MRI

PURPOSE

To establish the incidence, size, zonal location and Gleason Score(GS)/Gleason Grade Group(GG) of sparse versus dense prostate cancer (PCa) lesions and to identify the imaging characteristics of sparse versus dense cancers on multiparametric MRI (mpMRI).

METHODS

Seventy-six men with untreated PCa were scanned prior to prostatectomy with endorectal-coil 3 T MRI including T2-weighted imaging, diffusion-weighted imaging and dynamic contrast-enhanced MRI. Cancerous regions were outlined and graded on the whole-mount, processed specimens, with tissue compositions estimated. Regions with cancer comprising <50 % and ≥ 50 % of the tissue were considered sparse and dense respectively. Regions of interest (ROI) were manually drawn on T2-weighted MRI. Within each patient, area-weighted ROI averages were calculated for each imaging measure for each tissue type, GS/GG, and sparse/dense composition.

RESULTS

A large number of cancer regions were identified on histopathology (n = 1193: 939 (peripheral zone (PZ)) and 254 (transition zone (TZ))). Thirty-seven percent of these lesions were sparse. Sparse lesions were primarily low-grade with the majority of PZ and 100 % of TZ sparse lesions ≤GS3 + 3/GG1. Dense lesions were significantly larger than sparse lesions in both PZ and TZ, p < 0.0001. On imaging, 246/45 PZ and 109/8 TZ dense/sparse 2D cancerous ROIs were drawn. Sparse GS3 + 3 and sparse ≥GS3 + 4 cancers did not have significantly different MRI intensities to dense GS3 + 3 cancers, while sparse GS3 + 3/GG1 cancers differed from benign, p < 0.05.

CONCLUSION

Histopathologically identified prostate cancer lesions were sparse in 37 % of cases. Sparse cancers were entirely low grade in TZ and predominantly low-grade in PZ and generally small, thus likely posing lower risk for spread and progression than dense lesions. Sparse lesions were not distinguishable from dense lesions on mpMRI, but could be distinguished from benign tissues.

Full resolution reconstruction of whole-mount sections from digitized individual tissue fragments

Whole-mount sectioning is a technique in histopathology where a full slice of tissue, such as a transversal cross-section of a prostate specimen, is prepared on a large microscope slide without further sectioning into smaller fragments. Although this technique can offer improved correlation with pre-operative imaging and is paramount for multimodal research, it is not commonly employed due to its technical difficulty, associated cost and cumbersome integration in (digital) pathology workflows. In this work, we present a computational tool named PythoStitcher which reconstructs artificial whole-mount sections from digitized tissue fragments, thereby bringing the benefits of whole-mount sections to pathology labs currently unable to employ this technique. Our proposed algorithm consists of a multi-step approach where it (i) automatically determines how fragments need to be reassembled, (ii) iteratively optimizes the stitch using a genetic algorithm and (iii) efficiently reconstructs the final artificial whole-mount section on full resolution (0.25 µm/pixel). PythoStitcher was validated on a total of 198 cases spanning five datasets with a varying number of tissue fragments originating from different organs from multiple centers. PythoStitcher successfully reconstructed the whole-mount section in 86-100% of cases for a given dataset with a residual registration mismatch of 0.65-2.76 mm on automatically selected landmarks. It is expected that our algorithm can aid pathology labs unable to employ whole-mount sectioning through faster clinical case evaluation and improved radiology-pathology correlation workflows.

Quantification of cross-vendor variation in ADC measurements in vendor-specific prostate MRI-protocols

PURPOSE

The purpose of this study was to quantify the variability of Apparent Diffusion Coefficient (ADC) and test if there were statistically significant differences in ADC between MRI systems and sequences.

METHOD

With a two-chamber cylindrical ADC phantom with fixed ADC values (1,000 and 1,600x10^-6 mm²/s) a single-shot (ss) Echo Planar Imaging (EPI), a multi-shot EPI, a reduced field of view DWI (zoom) and a Turbo Spin Echo DWI sequence were tested in six MRI systems from three vendors at 1.5 T and 3 T. Technical parameters were according to Prostate Imaging Reporting and Data System Version 2.1. ADC maps were calculated by vendor specific algorithms. Absolute and relative differences in ADC from the phantom-ADC were calculated and differences between sequences were tested.

RESULTS

At 3 T absolute differences from phantom given ADC (∼1,000 and ∼ 1,600x10^-6 mm²/s) were -83 - 42x10^-6 mm²/s (-8.3%-4.2%) and -48 - 15x10^-6 mm²/s (-3%-0.9%), respectively and at 1.5 T absolute differences were -81 - 26x10^-6 mm²/s (-2.6%-8.1%) and -74 - 67x10^-6 mm²/s (-4.6%-4.2%), respectively. Significant statistical differences in ADC measurements were identified between vendors in all sequences except for ssEPI and zoom at 3 T in the 1,600x10^-6 mm²/s phantom chamber. Significant differences were also identified between ADC measurements at 1.5 T and 3 T in some of the sequences and vendors, but not all.

CONCLUSION

The variation of ADC between different MRI systems and prostate specific DWI sequences is limited in this phantom study and without apparent clinical relevance. However, prospective multicenter studies of prostate cancer patients are needed for further investigation.

Prostate cancer lesion detection, volume quantification and high-grade cancer differentiation using cancer risk maps derived from multiparametric MRI with histopathology as the reference standard

Multi-parametric MRI (mpMRI) has proven itself a clinically useful tool to assess prostate cancer (PCa). Our objective was to generate PCa risk maps to quantify the volume and location of both all PCa and high grade (Gleason grade group ≥ 3) PCa. Such capabilities would aid physicians and patients in treatment decisions, targeting biopsy, and planning focal therapy. A cohort of men with biopsy proven prostate cancer and pre-prostatectomy mpMRI were studied. PCa and benign ROIs (1524) were identified on mpMRI and histopathology with histopathology serving as the reference standard. Logistic regression models were created to differentiate PCa from benign tissues. The MRI images were registered to ensure correct overlay. The cancer models were applied to each image voxel within prostates to create probability maps of cancer and of high-grade cancer. Use of an optimum probability threshold quantified PCa volume for all lesions >0.1 cc. Accuracies were calculated using area under the curve (AUC) for the receiver operating characteristic (ROC). The PCa models utilized apparent diffusion coefficient (ADC), T2 weighted (T2W), dynamic contrast-enhanced MRI (DCE MRI) enhancement slope, and DCE MRI washout as the statistically significant MRI scans. Application of the PCa maps method provided total PCa volume and individual lesion volumes. The AUCs derived from lesion analysis were 0.91 for all PCa and 0.73 for high-grade PCa. At the optimum threshold, the PCa maps detected 135 / 150 (90%) histopathological lesions >0.1 cc. This study showed the feasibility of cancer risk maps, created from pre-prostatectomy, mpMR images validated with histopathology, to detect PCa lesions >0.1 cc. The method quantified the volume of cancer within the prostate. Method improvements were identified by determining root causes for over and underestimation of cancer volumes. The maps have the potential for improved non-invasive capability in quantitative detection, localization, volume estimation, and MRI characterization of PCa.

Mutually communicated model based on multi-parametric MRI for automated segmentation and classification of prostate cancer

BACKGROUND

Multiparametric magnetic resonance imaging (mp-MRI) is introduced and established as a noninvasive alternative for prostate cancer (PCa) detection and characterization.

PURPOSE

To develop and evaluate a mutually communicated deep learning segmentation and classification network (MC-DSCN) based on mp-MRI for prostate segmentation and PCa diagnosis.

METHODS

The proposed MC-DSCN can transfer mutual information between segmentation and classification components and facilitate each other in a bootstrapping way. For classification task, the MC-DSCN can transfer the masks produced by the coarse segmentation component to the classification component to exclude irrelevant regions and facilitate classification. For segmentation task, this model can transfer the high-quality localization information learned by the classification component to the fine segmentation component to mitigate the impact of inaccurate localization on segmentation results. Consecutive MRI exams of patients were retrospectively collected from two medical centers (referred to as center A and B). Two experienced radiologists segmented the prostate regions, and the ground truth of the classification refers to the prostate biopsy results. MC-DSCN was designed, trained, and validated using different combinations of distinct MRI sequences as input (e.g., T2-weighted and apparent diffusion coefficient) and the effect of different architectures on the network's performance was tested and discussed. Data from center A were used for training, validation, and internal testing, while another center's data were used for external testing. The statistical analysis is performed to evaluate the performance of the MC-DSCN. The DeLong test and paired t-test were used to assess the performance of classification and segmentation, respectively.

RESULTS

In total, 134 patients were included. The proposed MC-DSCN outperforms the networks that were designed solely for segmentation or classification. Regarding the segmentation task, the classification localization information helped to improve the IOU in center A: from 84.5% to 87.8% (p < 0.01) and in center B: from 83.8% to 87.1% (p < 0.01), while the area under curve (AUC) of PCa classification was improved in center A: from 0.946 to 0.991 (p < 0.02) and in center B: from 0.926 to 0.955 (p < 0.01) as a result of the additional information provided by the prostate segmentation.

CONCLUSION

The proposed architecture could effectively transfer mutual information between segmentation and classification components and facilitate each other in a bootstrapping way, thus outperforming the networks designed to perform only one task.

Dictionary learning compressed sensing reconstruction: pilot validation of accelerated echo planar J-resolved spectroscopic imaging in prostate cancer

Objectives

This study aimed at developing dictionary learning (DL) based compressed sensing (CS) reconstruction for randomly undersampled five-dimensional (5D) MR Spectroscopic Imaging (3D spatial + 2D spectral) data acquired in prostate cancer patients and healthy controls, and test its feasibility at 8x and 12x undersampling factors.

Materials and methods

Prospectively undersampled 5D echo-planar J-resolved spectroscopic imaging (EP-JRESI) data were acquired in nine prostate cancer (PCa) patients and three healthy males. The 5D EP-JRESI data were reconstructed using DL and compared with gradient sparsity-based Total Variation (TV) and Perona-Malik (PM) methods. A hybrid reconstruction technique, Dictionary Learning-Total Variation (DLTV), was also designed to further improve the quality of reconstructed spectra.

Results

The CS reconstruction of prospectively undersampled (8x and 12x) 5D EP-JRESI data acquired in prostate cancer and healthy subjects were performed using DL, DLTV, TV and PM. It is evident that the hybrid DLTV method can unambiguously resolve 2D J-resolved peaks including myo-inositol, citrate, creatine, spermine and choline.

Conclusion

Improved reconstruction of the accelerated 5D EP-JRESI data was observed using the hybrid DLTV. Accelerated acquisition of in vivo 5D data with as low as 8.33% samples (12x) corresponds to a total scan time of 14 min as opposed to a fully sampled scan that needs a total duration of 2.4 h (TR = 1.2 s, 32 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_{x}$$\end{document}kx×16 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_{y}$$\end{document}ky×8 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_{z}$$\end{document}kz, 512 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${t}_{2}$$\end{document}t2 and 64 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${t}_{1}$$\end{document}t1).

Supplementary Information

The online version contains supplementary material available at 10.1007/s10334-022-01029-z.