Data-driven separation of MRI signal components for tissue characterization

Prediction of overall survival in patients with locally advanced pancreatic cancer using longitudinal diffusion-weighted MRI

Background and purpose

Biomarkers for prediction of outcome in patients with pancreatic cancer are wanted in order to personalize the treatment. This study investigated the value of longitudinal diffusion-weighted magnetic resonance imaging (DWI) for prediction of overall survival (OS) in patients with locally advanced pancreatic cancer (LAPC) treated with stereotactic body radiotherapy (SBRT).

Materials and methods

The study included 45 patients with LAPC who received 5 fractions of 10 Gy on a 1.5T MRI-Linac. DWI was acquired prior to irradiation at each fraction. The analysis included baseline values and time-trends of the apparent diffusion coefficient (ADC) and DWI parameters obtained using a decomposition method. A multivariable Cox proportional hazards model for OS was made using best-subset selection, using cross-validation based on Bootstrap.

Results

The median OS from the first day of SBRT was 15.5 months (95% CI: 13.2-20.6), and the median potential follow-up time was 19.8 months. The best-performing multivariable model for OS included two decomposition-based DWI parameters: one baseline and one time-trend parameter. The C-Harrell index describing the model's discriminating power was 0.754. High baseline ADC values were associated with reduced OS, whereas no association between the ADC time-trend and OS was observed.

Conclusion

Decomposition-based DWI parameters indicated value in the prediction of OS in LAPC. A DWI time-trend parameter was included in the best-performing model, indicating a potential benefit of acquiring longitudinal DWI during the SBRT course. These findings support both baseline and longitudinal DWI as candidate prognostic biomarkers, which may become tools for personalization of the treatment of patients with LAPC.

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 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₂-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 area under curve (AUC) of 0.999,P< 0.001 and predict the tumor volume change with a correlation coefficient of 0.513,P= 0.034. A classification of the human brain metastases into responders and non-responders resulted in an AUC of 0.74,P= 0.065.Significance.A general data processing framework for analyses of longitudinal MRI data has been developed and applications were demonstrated by classification of tumor type and prediction of radiotherapy response. Further, as part of the assessment, the merits of msNMF for tumor tissue decomposition were demonstrated.

Self-labelled encoder-decoder (SLED) for multi-echo gradient echo-based myelin water imaging

PURPOSE

Reconstruction of high quality myelin water imaging (MWI) maps is challenging, particularly for data acquired using multi-echo gradient echo (mGRE) sequences. A non-linear least squares fitting (NLLS) approach has often been applied for MWI. However, this approach may produce maps with limited detail and, in some cases, sub-optimal signal to noise ratio (SNR), due to the nature of the voxel-wise fitting. In this study, we developed a novel, unsupervised learning method called self-labelled encoder-decoder (SLED) to improve gradient echo-based MWI data fitting.

METHODS

Ultra-high resolution, MWI data was collected from five mouse brains with variable levels of myelination, using a mGRE sequence. Imaging data was acquired using a 7T preclinical MRI system. A self-labelled, encoder-decoder network was implemented in TensorFlow for calculation of myelin water fraction (MWF) based on the mGRE signal decay. A simulated MWI phantom was also created to evaluate the performance of MWF estimation.

RESULTS

Compared to NLLS, SLED demonstrated improved MWF estimation, in terms of both stability and accuracy in phantom tests. In addition, SLED produced less noisy MWF maps from high resolution MR microscopy images of mouse brain tissue. It specifically resulted in lower noise amplification for all mouse genotypes that were imaged and yielded mean MWF values in white matter ROIs that were highly correlated with those derived from standard NLLS fitting. Lastly, SLED also exhibited higher tolerance to low SNR data.

CONCLUSION

Due to its unsupervised and self-labeling nature, SLED offers a unique alternative to analyze gradient echo-based MWI data, providing accurate and stable MWF estimations.