An image-driven parameter estimation problem for a reaction–diffusion glioma growth model with mass effects

We present a framework for modeling gliomas growth and their mechanical impact on the surrounding brain tissue (the so-called, mass-effect). We employ an Eulerian continuum approach that results in a strongly coupled system of nonlinear Partial Differential Equations (PDEs): a reaction-diffusion model for the tumor growth and a piecewise linearly elastic material for the background tissue. To estimate unknown model parameters and enable patient-specific simulations we formulate and solve a PDE-constrained optimization problem. Our two main goals are the following: (1) to improve the deformable registration from images of brain tumor patients to a common stereotactic space, thereby assisting in the construction of statistical anatomical atlases; and (2) to develop predictive capabilities for glioma growth, after the model parameters are estimated for a given patient. To our knowledge, this is the first attempt in the literature to introduce an adjoint-based, PDE-constrained optimization formulation in the context of image-driven modeling spatio-temporal tumor evolution. In this paper, we present the formulation, and the solution method and we conduct 1D numerical experiments for preliminary evaluation of the overall formulation/methodology.

[1]  M. Westphal,et al.  Glioma invasion in the central nervous system. , 1996, Neurosurgery.

[2]  E R Laws,et al.  Human malignant astrocytoma xenografts migrate in rat brain: A model for central nervous system cancer research , 1989, Journal of neuroscience research.

[3]  Christos Davatzikos,et al.  A Comparative Study of Biomechanical Simulators in Deformable Registration of Brain Tumor Images , 2008, IEEE Transactions on Biomedical Engineering.

[4]  Max Gunzburger,et al.  Perspectives in flow control and optimization , 1987 .

[5]  Dinggang Shen,et al.  Deformable registration of brain tumor images via a statistical model of tumor-induced deformation , 2006, Medical Image Anal..

[6]  T. Deisboeck,et al.  Complex dynamics of tumors: modeling an emerging brain tumor system with coupled reaction–diffusion equations , 2003 .

[7]  D L S McElwain,et al.  A history of the study of solid tumour growth: The contribution of mathematical modelling , 2004, Bulletin of mathematical biology.

[8]  Helen M. Byrne,et al.  A Multiple Scale Model for Tumor Growth , 2005, Multiscale Model. Simul..

[9]  Guangwei Yuan,et al.  Convergence and stability of explicit/implicit schemes for parabolic equations with discontinuous coefficients , 2004 .

[10]  Claudio H. Sibata,et al.  Patient-specific tumor prognosis prediction via multimodality imaging , 1996, Medical Imaging.

[11]  Olivier Pironneau,et al.  A FICTITIOUS DOMAIN BASED GENERAL PDE SOLVER , 2004 .

[12]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[13]  Tamara G. Kolda,et al.  Revisiting Asynchronous Parallel Pattern Search for Nonlinear Optimization , 2005, SIAM J. Optim..

[14]  Willem Hundsdorfer,et al.  Numerical time integration for air pollution models , 1998 .

[15]  G. Biros,et al.  A Framework for Soft Tissue Simulations with Application to Modeling Brain Tumor Mass-Effect in 3 D Images , 2006 .

[16]  R. LeVeque Numerical methods for conservation laws , 1990 .

[17]  K. Rohr,et al.  Biomechanical modeling of the human head for physically based, nonrigid image registration , 1999, IEEE Transactions on Medical Imaging.

[18]  J. Craggs Applied Mathematical Sciences , 1973 .

[19]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[20]  Tamara G. Kolda,et al.  Algorithm 856: APPSPACK 4.0: asynchronous parallel pattern search for derivative-free optimization , 2006, TOMS.

[21]  G. Biros,et al.  Fast Solvers for Soft Tissue Simulation with Application to Construction of Brain Tumor Atlases , 2007 .

[22]  Serge Duckett,et al.  The Pathology of the aging human nervous system , 1991 .

[23]  Philippe Tracqui,et al.  MODELLING THREE-DIMENSIONAL GROWTH OF BRAIN TUMOURS FROM TIME SERIES OF SCANS , 1999 .

[24]  Carol P. Geer,et al.  Interstitial fluid flow along white matter tracts: A potentially important mechanism for the dissemination of primary brain tumors , 1997, Journal of Neuro-Oncology.

[25]  Christos Davatzikos,et al.  A Biomechanical Model of Soft Tissue Deformation, with Applications to Non-rigid Registration of Brain Images with Tumor Pathology , 1998, MICCAI.

[26]  R. Kress Linear Integral Equations , 1989 .

[27]  Helen M. Byrne,et al.  A two-phase model of solid tumour growth , 2003, Appl. Math. Lett..

[28]  Bruce T. Murray,et al.  Modeling tumor growth: a computational approach in a continuum framework , 2005 .

[29]  R. Kauppinen,et al.  Monitoring thymidine kinase and ganciclovir-induced changes in rat malignant glioma in vivo by nuclear magnetic resonance imaging. , 1998, Cancer gene therapy.

[30]  George Biros,et al.  DYNAMIC DATA-DRIVEN INVERSION FOR TERASCALE SIMULATIONS: REAL-TIME IDENTIFICATION OF AIRBORNE CONTAMINANTS , 2005, ACM/IEEE SC 2005 Conference (SC'05).

[31]  J. Murray,et al.  Virtual and real brain tumors: using mathematical modeling to quantify glioma growth and invasion , 2003, Journal of the Neurological Sciences.

[32]  D. L. Sean McElwain,et al.  A Mixture Theory for the Genesis of Residual Stresses in Growing Tissues I: A General Formulation , 2005, SIAM J. Appl. Math..

[33]  M. Chaplain,et al.  Mathematical Modelling of Tissue Invasion , 2003 .

[34]  Christos Davatzikos,et al.  Modeling Glioma Growth and Mass Effect in 3D MR Images of the Brain , 2007, MICCAI.

[35]  Christos Davatzikos,et al.  Nonlinear elastic registration of brain images with tumor pathology using a biomechanical model [MRI] , 1999, IEEE Transactions on Medical Imaging.

[36]  L. G. Stern,et al.  Fractional step methods applied to a chemotaxis model , 2000, Journal of mathematical biology.

[37]  J. Murray,et al.  A quantitative model for differential motility of gliomas in grey and white matter , 2000, Cell proliferation.

[38]  Christos Davatzikos,et al.  Finite Element Modeling of Brain Tumor Mass-Effect from 3D Medical Images , 2005, MICCAI.

[39]  D. Silbergeld,et al.  Assessment of brain tumor cell motility in vivo and in vitro. , 1995, Journal of neurosurgery.

[40]  Thomas S. Deisboeck,et al.  Simulating ‘structure–function’ patterns of malignant brain tumors , 2004 .

[41]  Laurent Capelle,et al.  Continuous growth of mean tumor diameter in a subset of grade II gliomas , 2003, Annals of neurology.

[42]  S Torquato,et al.  Simulated brain tumor growth dynamics using a three-dimensional cellular automaton. , 2000, Journal of theoretical biology.

[43]  Dinggang Shen,et al.  Registration of brain images with tumors: towards the construction of statistical atlases for therapy planning , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[44]  Hervé Delingette,et al.  Realistic simulation of the 3-D growth of brain tumors in MR images coupling diffusion with biomechanical deformation , 2005, IEEE Transactions on Medical Imaging.