Estimating the extent of glioblastoma invasion

Christian Engwer,Michael Wenske

DOI: https://doi.org/10.48550/arXiv.2001.05369

2020-01-15

Abstract:Glioblastoma Multiforme is a malignant brain tumor with poor prognosis. There have been numerous attempts to model the invasion of tumorous glioma cells via partial differential equations in the form of advection-diffusion-reaction equations. The patient-wise parametrisation of these models, and their validation via experimental data has been found to be difficult, as time sequence measurements are generally missing. Also the clinical interest lies in the actual (invisible) tumor extent for a particular MRI/DTI scan and not in a predictive estimate. Therefore we propose a stationalised approach to estimate the extent of glioblastoma (GBM) invasion at the time of a given MRI/DTI scan. The underlying dynamics can be derived from an instationary GBM model, falling into the wide class of advection-diffusion-reaction equations. The stationalisation is introduced via an analytical solution of the Fisher-KPP equation, the simplest model in the considered model class. We investigate the applicability in 1D and 2D, in the presence of inhomogeneous diffusion coefficients and on a real 3D DTI-dataset.

Cell Behavior

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to estimate the actual invasive range of glioblastoma multiforme (GBM) at a given MRI/DTI scan time point. Specifically, the authors propose a static method to overcome the difficulties in parameterizing existing models and the lack of time - series data, so as to be able to estimate the invisible invasive range of the tumor at a specific time point. This method is based on the non - steady - state GBM model and is realized by introducing the analytical solution of the Fisher - KPP equation for staticization. The paper also explores the applicability of this method in one - dimensional and two - dimensional cases, as well as its performance on non - uniform diffusion coefficients and real 3D DTI data sets.

### Background and Motivation
Glioblastoma is a malignant brain tumor with a poor prognosis. Many studies have attempted to model the invasive behavior of tumor cells through partial differential equations such as the convection - diffusion - reaction equation. However, patient - specific parameterization of these models and verification through experimental data are very difficult because time - series measurement data are usually lacking. In addition, the clinical focus is on the actual (invisible) tumor range at the time of a specific MRI/DTI scan, rather than on predictive estimation.

### Research Objectives
- **Estimate the Tumor Invasive Range**: Propose a static method to estimate the GBM invasive range at a specific time point.
- **Simplify Parameterization**: Reduce the dependence on absolute quantitative parameterization through the static method, thereby alleviating the parameterization problem.
- **Verify the Effectiveness of the Method**: Verify the effectiveness of the static method under different conditions (such as one - dimensional, two - dimensional, non - uniform diffusion coefficients, and real 3D DTI data sets).

### Methods
1. **Model Description**:
- Use the fully anisotropic convection - diffusion - reaction equation to describe the invasive dynamics of GBM.
- The equation is in the form:
\[
\frac{\partial u}{\partial t}-\nabla\cdot(D_t(x)\nabla u)-\nabla\cdot((\nabla\cdot D_t(x))u)=\rho u(1 - u)\quad\text{in }\Omega\times T
\]
- Where \(u\) represents the tumor density, \(D_t(x)\) is the diffusion tensor, and \(\rho\) is the growth rate.

2. **Static Method**:
- Based on the analytical solution of the Fisher - KPP equation, introduce a static term.
- Transform the equation into a static form through the coordinate transformation \(x = x'-vt\):
\[
0=\Delta u+u(1 - u)-v\cdot\nabla u
\]
- The form of the static term \(p(u)\) is:
\[
p(u)=\vert v\vert\sqrt{\frac{2}{3}}(1-\sqrt{u})u
\]

3. **Numerical Method**:
- Use the finite element method for spatial discretization and the implicit Euler method for time discretization.
- Implement numerical simulation through the DUNE software framework.

### Results
- **Gradient Distribution**: The consistency between the gradient distribution of the static method and the analytical solution was verified under different conditions (such as one - dimensional, two - dimensional, non - uniform diffusion coefficients).
- **Error Analysis**: The error of the static method was evaluated by comparing the tumor front edges of the static method and the full - transient simulation.
- **Practical Application**: The effectiveness of the method was verified on a real 3D DTI data set.

### Conclusion
This paper proposes a static method that can estimate the invasive range of GBM at a specific time point, effectively alleviates the parameterization problem, and verifies its effectiveness under various conditions. This provides a new tool for more accurate estimation of the tumor invasive range in clinical practice.