Decomposition-based framework for tumor classification and prediction of treatment response from longitudinal MRI

Objective: In the field of radiation oncology, the benefit of MRI goes beyond that of providing high soft-tissue contrast images for staging and treatment planning. With the recent clinical introduction of hybrid MRI linear accelerators (MR-Linacs) it has become feasible to map physiological parameters describing diffusion, perfusion, and relaxation during the entire course of radiotherapy, for example. However, advanced data analysis tools are required for extracting qualified prognostic and predictive imaging biomarkers from longitudinal MRI data. In this study, we propose a new prediction framework tailored to exploit temporal dynamics of tissue features from repeated measurements. We demonstrate the framework using a newly developed decomposition method for tumor characterization. Approach: Two previously published MRI datasets with multiple measurements during and after radiotherapy, were used for development and testing: T$_2$-weighted multi-echo images obtained for two mouse models of pancreatic cancer, and diffusion-weighted images for patients with brain metastases. Initially, the data was decomposed using the novel monotonous slope non-negative matrix factorization (msNMF) tailored for MR data. The following processing consisted of a tumor heterogeneity assessment using descriptive statistical measures, robust linear modelling to capture temporal changes of these, and finally logistic regression analysis for stratification of tumors and volumetric outcome. Main Results: The framework was able to classify the two pancreatic tumor types with an \textit{area under curve} AUC of 0.999, $P