Automatic MRI-Driven Model Calibration for Advanced Brain Tumor Progression Analysis

Our objective is the calibration of mathematical tumor growth models from a single multiparametric scan. The target problem is the analysis of preoperative Glioblastoma (GBM) scans. To this end, we present a fully automatic tumor-growth calibration methodology that integrates a single-species reaction-diffusion partial differential equation (PDE) model for tumor progression with multiparametric Magnetic Resonance Imaging (mpMRI) scans to robustly extract patient specific biomarkers i.e., estimates for (i) the tumor cell proliferation rate, (ii) the tumor cell migration rate, and (iii) the original, localized site(s) of tumor initiation. Our method is based on a sparse reconstruction algorithm for the tumor initial location (TIL). This problem is particularly challenging due to nonlinearity, ill-posedeness, and ill conditioning. We propose a coarse-to-fine multi-resolution continuation scheme with parameter decomposition to stabilize the inversion. We demonstrate robustness and practicality of our method by applying the proposed method to clinical data of 206 GBM patients. We analyze the extracted biomarkers and relate tumor origin with patient overall survival by mapping the former into a common atlas space. We present preliminary results that suggest improved accuracy for prediction of patient overall survival when a set of imaging features is augmented with estimated biophysical parameters. All extracted features, tumor initial positions, and biophysical growth parameters are made publicly available for further analysis. To our knowledge, this is the first fully automatic scheme that can handle multifocal tumors and can localize the TIL to a few millimeters.

[1]  Christos Davatzikos,et al.  An image-driven parameter estimation problem for a reaction–diffusion glioma growth model with mass effects , 2008, Journal of mathematical biology.

[2]  George Biros,et al.  Where did the tumor start? An inverse solver with sparse localization for tumor growth models , 2019, Inverse problems.

[3]  G. Strang On the Construction and Comparison of Difference Schemes , 1968 .

[4]  Kurt Keutzer,et al.  A Novel Domain Adaptation Framework for Medical Image Segmentation , 2018, BrainLes@MICCAI.

[5]  Thomas Welzel,et al.  A comparison of long-term survivors and short-term survivors with glioblastoma, subventricular zone involvement: a predictive factor for survival? , 2014, Radiation oncology.

[6]  Klaudius Scheufele Coupling Schemes and Inexact Newton for Multi-Physics and Coupled Optimization Problems , 2019 .

[7]  Kristin R. Swanson,et al.  The Evolution of Mathematical Modeling of Glioma Proliferation and Invasion , 2007, Journal of neuropathology and experimental neurology.

[8]  G. Biros,et al.  An inverse problem formulation for parameter estimation of a reaction–diffusion model of low grade gliomas , 2014, Journal of mathematical biology.

[9]  Alfredo Quinones-Hinojosa,et al.  Relationship of glioblastoma multiforme to the lateral ventricles predicts survival following tumor resection , 2008, Journal of Neuro-Oncology.

[10]  Brent D. Weinberg,et al.  Location of subventricular zone recurrence and its radiation dose predicts survival in patients with glioblastoma , 2018, Journal of Neuro-Oncology.

[11]  Christos Davatzikos,et al.  Coupling brain-tumor biophysical models and diffeomorphic image registration. , 2017, Computer methods in applied mechanics and engineering.

[12]  George Biros,et al.  Image-Driven Biophysical Tumor Growth Model Calibration , 2019, SIAM J. Sci. Comput..

[13]  Christos Davatzikos,et al.  A Framework for Scalable Biophysics-based Image Analysis , 2017, SC17: International Conference for High Performance Computing, Networking, Storage and Analysis.

[14]  Adolf Pfefferbaum,et al.  The SRI24 multichannel atlas of normal adult human brain structure , 2009, Human brain mapping.