A multivariate curve resolution analysis of multicenter proton spectroscopic imaging of the prostate for cancer localization and assessment of aggressiveness
Angeliki Stamatelatou,Carlo Giuseppe Bertinetto,Jeroen J Jansen,Geert Postma,Kirsten Margrete Selnaes,Tone F Bathen,Arend Heerschap,Tom W J Scheenen,PCa-MAP consortium,Ulrike I Attenberger,Pascal A T Baltzer,Jurgen J Fütterer,Masoom A Haider,Thomas H Helbich,Berthold Kiefer,Marnix C Maas,Katarzyna J Macura,Daniel J A Margolis,Anwar R Padhani,Stephen H Polanec,Marleen Praet,Tom W Scheenen,Stefan O Schoenberg,Kirsten M Selnaes,Theodorus H van der Kwast,Geert M Villeirs,Trond Viset,Heninrich von Busch
DOI: https://doi.org/10.1002/nbm.5062
Abstract:In this study, we investigated the potential of the multivariate curve resolution alternating least squares (MCR-ALS) algorithm for analyzing three-dimensional (3D) 1 H-MRSI data of the prostate in prostate cancer (PCa) patients. MCR-ALS generates relative intensities of components representing spectral profiles derived from a large training set of patients, providing an interpretable model. Our objectives were to classify magnetic resonance (MR) spectra, differentiating tumor lesions from benign tissue, and to assess PCa aggressiveness. We included multicenter 3D 1 H-MRSI data from 106 PCa patients across eight centers. The patient cohort was divided into a training set (N = 63) and an independent test set (N = 43). Singular value decomposition determined that MR spectra were optimally represented by five components. The profiles of these components were extracted from the training set by MCR-ALS and assigned to specific tissue types. Using these components, MCR-ALS was applied to the test set for a quantitative analysis to discriminate tumor lesions from benign tissue and to assess tumor aggressiveness. Relative intensity maps of the components were reconstructed and compared with histopathology reports. The quantitative analysis demonstrated a significant separation between tumor and benign voxels (t-test, p < 0.001). This result was achieved including voxels with low-quality MR spectra. A receiver operating characteristic analysis of the relative intensity of the tumor component revealed that low- and high-risk tumor lesions could be distinguished with an area under the curve of 0.88. Maps of this component properly identified the extent of tumor lesions. Our study demonstrated that MCR-ALS analysis of 1 H-MRSI of the prostate can reliably identify tumor lesions and assess their aggressiveness. It handled multicenter data with minimal preprocessing and without using prior knowledge or quality control. These findings indicate that MCR-ALS can serve as an automated tool to assess the presence, extent, and aggressiveness of tumor lesions in the prostate, enhancing diagnostic capabilities and treatment planning of PCa patients.