An MRI-informed poromechanical model for organ-scale prediction of glioma growth

Gliomas constitute one of the most aggressive and heterogeneous forms of brain tumors, posing major challenges for understanding their biology and developing effective treatments. Animal models enable the collection of rich longitudinal datasets describing tumor dynamics, which can be integrated within mathematical models to elucidate the biological mechanisms governing tumor growth. While most formulations rely on reaction–diffusion systems with limited insight on tissue deformation and fluid transport, we propose a magnetic resonance imaging (MRI)-informed, poroelastic model to describe C6 glioma growth in rats. We use data from animals ( n = 4 ) that were imaged five times after intracranial injection of cancer cells. Each MRI dataset includes (i) anatomical T 1 -weighted data for brain and tumor segmentation and to assign mechanical properties; (ii) diffusion-weighted MRI, which enables estimation of the fraction of each voxel that is tumor; and (iii) dynamic contrast-enhanced MRI, which informs permeability as well as vascular and liquid fraction maps. Using finite-element simulations, model calibration for each rat uses the Gauss–Newton method informed by the first three MRI datasets. Tumor forecasts are validated by assessing model-data agreement on the remaining two MRI datasets. Our results show relative tumor volume errors between 0.94% and 11.27% at calibration, and prediction errors between 4.73% and 36.03%. Additionally, Dice scores ranged from 0.80 to 0.93 during calibration, and from 0.75 to 0.93 during validation. Thus, our results suggest that our poromechanical model can describe C6 glioma growth. This study provides a first step towards a patient-specific, multiscale model of the spatiotemporal poromechanics underlying glioma progression and therapeutic response.