Image-based prediction of tumor growth

Participants: Thierry Colin, Olivier Saut, Julie Joie, Patricio Andrés Cumsille Atala, Julien Jouganous, Guillaume Lefèvre.

Clinicians involved: J. Palussière, F. Cornelis, G. Kantor, H. Loiseau and H. Fathallah .


The idea of this project is to predict the time development of a clinical tumor or the response to a treatment using nonlinear PDE models and a series of CT-Scan or MRI of a patient as inputs. We start by extracting from our large model a "simple" PDE model involving few parameters (let's say around 5 independent parameters). Then we try to find the "best" values of the parameters that allow to match with the series of image by solving an optimization problem; then me make a prediction using this set of parameters. We have tested this strategy on metastasis to the lung from distant tumors (kidney, bladder, thyroid). We dispose of time series of CT-scans of patients that are only under monitoring (i.e. without treatment). We use only two CT-scans to parametrize our problem and try to recover the following ones.


Lung metastasis
The test case presented below concerns a metastasis to the lung of a bladder tumor. On the left are the 3 CT-scans that constitute our data set. No treatment was given to the patient during this period. We used only the image of June and September to perform the simulation. In the middle plot, measured volumes of the metastasis (based on the scan) in time are the circle points while the continuous line is the volume predicted by the simulation. On the right images, the tumor given by the simulation for September and December is in red (top row). Other examples for metastasis to the lung are given in the publications.

T. Colin, A. Iollo, D. Lombardi, O. Saut System Identification in Tumor Growth Modeling Using Semi-empirical Eigenfunctions. Math. Models Methods Appl. Sci. 22, 1250003 (2012).
T. Colin, A. Iollo, D. Lombardi and O. Saut, Prediction of the Evolution of Thyroidal Lung Nodules Using a Mathematical Model, ERCIM News, No 82, pp. 37-38, July 2010.



Lung metastasis growth prediction. Left: Data. Right: Model simulations and predictions.


Liver metastasis
We have also some preliminary results concerning metastasis to the lung of a GIST (Gastro-Intestinal Stromal Tumor). When the metastasis is discovered, the patient receives a targeted therapy called Imatinib until he escapes the treatment. He then receives a second therapy with Sunitinib until the next escape time. In the example below we show on the left the clinical data: time curve of the volume of the metastasis measured on the successive CT-scans and the CT-scans corresponding (from left to right and from top to bottom) to the control by the first treatment, then the escape and then the control but the second treatment. Right images is a fit of the model to the volume data.



Liver metastasis. Left: Data provided by the Institut Bergonié. Right: Model fit.


Glioblastoma (GBM)
With H. Fathallah (University of Alabama at Birmingham), we have developed a 3D model of glioblastoma expansion that shows the three layer structure of GBM: a necrotic core, a proliferative rim surrounded by a cloud of invasive cells. These elements can be seen on different sequences of MRI (T1, T1 gado and Flair). It is still an open problem to characterize precisely these elements on these images. Below is shown a simulation of GBM expansion together with an MRI and a biopsy of a glioblastoma.

Th. Colin, A. Iollo, J.-B. Lagaert and O. Saut, An inverse problem for the the recovery of the vascularization of a tumor.
Th. Colin, H. Fathallah, J.-B. Lagaert, O. Saut, A Multilayer Model for GBM: Effects of Invasive Cells and Anti-Angiogenesis on Growth. Submitted.


Glioblastoma. Left: Model simulations. Right: Data provided by the University of Alabama.