Journal of Biomedicine and Biotechnology

Volume 2012 (2012), Article ID 715812, 9 pages

doi:10.1155/2012/715812

Research Article

Simulating Radiotherapy Effect in High-Grade Glioma by Using Diffusive Modeling and Brain Atlases

Alexandros Roniotis,1,2 Kostas Marias,¹ Vangelis Sakkalis,¹ Georgios C. Manikis,¹ and Michalis Zervakis²

¹Institute of Computer Science, Foundation For Research and Technology-Hellas, 700 13 Heraklion, Crete, Greece

²Department of Electronic & Computer Engineering, Technical University of Crete, 73100 Chania, Greece

Received 21 February 2012; Revised 18 May 2012; Accepted 21 May 2012

Academic Editor: George E. Plopper

Copyright © 2012 Alexandros Roniotis et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Applying diffusive models for simulating the spatiotemporal change of concentration of tumour cells is a modern application of predictive oncology. Diffusive models are used for modelling glioblastoma, the most aggressive type of glioma. This paper presents the results of applying a linear quadratic model for simulating the effects of radiotherapy on an advanced diffusive glioma model. This diffusive model takes into consideration the heterogeneous velocity of glioma in gray and white matter and the anisotropic migration of tumor cells, which is facilitated along white fibers. This work uses normal brain atlases for extracting the proportions of white and gray matter and the diffusion tensors used for anisotropy. The paper also presents the results of applying this glioma model on real clinical datasets.

1. Introduction

Cancer causes more than 13% of all deaths worldwide with an estimation of 12 million deaths in 2030 [1, 2]. In United States, 2.5% of cancer deaths are caused by brain tumors [3], with grade IV astrocytoma, namely, glioblastoma multiforme (GBM), accounting for 23.1% of them [4–6]. Despite extensive research on GBM during the last decades, mortality has not changed significantly over the last years, with average life expectancy ranging 12 months after diagnosis and only 4% of treated patients being alive after 5 years [7–12].

Unfortunately, the detection rates of the exact boundaries of GBM with common imaging techniques, such as magnetic resonance imaging (MRI), X-ray computed tomography (CT), and positron emission tomography (PET), are still poor [12–14]. Unlike solid tumors, for which simple exponential or geometric expansion could represent expansion of tumor volume, the GBM growth rate cannot be determined as the classical doubling rate [15], because GBM does not form a solid and compact mass with cells invading the surrounding lesion.

Clearly, new mathematical formulations are necessary when studying this specific glioma case. Since early 90s, there has been a vast amount of research towards simulating and formulating the mechanisms of GBM development, both in macroscopic and microscopic levels. Microscopic models [16–20] study the intracellular biological interactions in cell level, while macroscopic models [21–27] study the tumor behavior, velocity, and mass deformation with using anatomical information derived from medical images.

Gompertzian and volumetric macroscopic models [28–30] did not produce realistic clinical cancer growth representations as they lack taking invasion of tumor cells into account. The class of macroscopic models that has been widely used for simulating this diffusive behavior is diffusive modeling [31]. These models make use of equations based on diffusion reaction schemes of glioma cell concentration, in order to couple the migration of GBM cells (diffusion term) and the proliferation of cells (reaction/source term).

The first diffusive GBM model was introduced in 1995 by Tracqui [21] who used the diffusion reaction equation (DRE) for simulating the spatiotemporal change of glioma cell concentration. The DRE writes as follows:

where is the tumor concentration in voxel at time , is the diffusion coefficient, and div are the gradient and divergence operators respectively, and is the net cell proliferation rate. Later, Swanson et al. [22] altered this equation by taking into account the high velocity of GBM in myelin sheaths by setting different diffusion coefficients in white and gray matter (WM, GM). The new equation was changed to

where is the local diffusion coefficient, being either or for being in GM and WM respectively. In that model, is assumed to be five times larger than for high grade glioma, according to the clinical observations of Giese et al. [32, 33].

The next change of the DRE equation was proposed by Price et al. who used T2- and Diffusion Tensor Imaging (DTI) MRIs to identify abnormalities caused by GBM [6, 34]. Therefore, Jbabdi et al. [23] introduced diffusion tensors instead of gradient diffusion coefficients. The new advanced equations for simulating anisotropic growth of GBM in WM writes as follows:

where is the local diffusion tensor which describes the directional tumor cell diffusion, that is, a 3-by-3 symmetric positive definite matrix. Continuing, Clatz et al. [25] used biomechanics to simulate the deformation of brain structures caused by the development of GBM. In 2007, Hogea et al. [27] introduced an advection term in (2), with a two-way coupling with the underlying tissue elastic deformation. Konukoglu et al. [26] adapted the model to successive sessions of MRIs by using Eikonal Equations for coupling the reaction-diffusion dynamics with tumor evolution.

The characteristics of GBM lead to aggressive treatment strategies that most often include surgery, irradiation and chemotherapy. The mathematical approaches that were used for simulating GBM growth inevitably raised the need for incorporating techniques for therapy simulation. Including therapy parameters in GBM models could help the clinicians optimize therapy schemes, as they could predict the response of patients to different therapeutic plans.

Indeed, simulating therapy can be performed by adjusting the proliferation term of the DRE, which reflects the rate of proliferating cells minus the dead cells due to treatment (). In order to simulate radiotherapy and chemotherapy, someone should adjust the treatment term . Swanson et al. [35] predicted survival time after resection by estimating the net rates of proliferation and diffusion, and their ratio from pre-treatment gadolinium-enhanced T1-weighted and T2-weighted MR images. Rockne et al. [36] simulated radiotherapy by introducing the radiobiology parameters and into , which are interpreted biologically as repairable single and lethal double-strand breaks to the cell’s DNA, respectively [37]. Lastly, some important works have been proposed for non diffusive models of glioma, taking into account the biological mechanisms involved in tumor and normal tissue [38, 39].

The most commonly used formalisms of are either the geometric law

the logistic law:

or the Gompertzian law:

where is the geometrical growth parameter and is the maximum tumor cell concentration parameter.

In 2011, Roniotis et al. [40, 41] used the proportions of white and GM for calculating the diffusion coefficient of the DRE, rather than using discrete diffusion coefficients [42]. The proportions of white and GM, as well as the diffusion tensor were extracted by the open SRI-24 brain atlas, thus no DTI processing was needed. This study extends our prior successes with in silico prediction of tumor growth to incorporate the effects of radiotherapy in real clinical datasets. We use the proportional model of Wahl et al. and apply the linear quadratic model [43] for modeling radiotherapy. Compared to previous publications dealing with tumor growth modeling, our approach includes several improvements:

(i)the use of proportional local diffusion coefficients extracted from brain atlases, instead of discrete diffusion coefficients,(ii)the application of isotropic diffusion in GM and anisotropic diffusion in WM on the proportional diffusive model [40],(iii)the application of the linear quadratic model on our previous model [40] for simulating radiotherapy,(iv)the simulation of radiotherapy in multiple fractions by using the linear quadratic model,(v)the evaluation of radiotherapy simulation by using Jaccard, Dice, and Volume similarity metrics.

2. Materials and Methods

2.1. Diffusive Modeling Using Brain Matter Proportions

This study makes use of the proportional diffusion model presented in [40], which is an extended version of (3) with

where is a symmetric matrix that describes the anisotropy of cell diffusion along the 3-dimensional directions for each position and denotes the proportional diffusion coefficient. In this case, writes as

where is the directional diffusion weight, which denotes the contribution of the respective direction to the anisotropic migration of GBM cells in position . Thus, denotes the contribution of each axis to the direction of white fibers, while actually denotes the velocity of cells along this direction. and can be approximated by using the method of Harpold et al. [44].

In order to make a total use of the proportions of WM and of GM that can be provided by brain atlases, such as the STI-24 atlas [45, 46], does not take two discrete values and [32]. Instead, the model uses a continuous pattern of coefficients by computing as follows:

This paper extends this idea by improving (5), so as to turn off anisotropic migration in GM tissue. In order to apply anisotropy only to WM tissues, the identity matrix is used for changing (5) to the following formalism:

2.2. The SRI-24 Atlas

The SRI24 atlas [45, 46] is an MRI-based atlas of normal adult human brain anatomy, generated by registering images of 24 normal brains. This atlas provides the proportions of WM and of GM in each position of the brain, which can be registered to the clinical datasets of the patients. The data is provided in 155 slices with pixels, for both stripped and unstripped skull.

Apart from the matter proportions, the atlas provides the dominant eigenvectors of the diffusion tensor. These eigenvectors have been extracted from DTI, by computing the covariance matrices of the distribution of the 3D Gaussian probability which simulates the diffusion of water. Thus, they can be directly used for simulating the anisotropic migration of glioma cells along white matter fibers, as they represent the directions towards which the water diffusion extends mostly.

By mapping clinical MRI data to a brain matter atlas, it is possible to approximate the required tissue information/composition in the tumor area, since this is not possible to extract from the real MRI. Hence, using local diffusion coefficients on real MRI images becomes possible, even if the lesion “hides” the underlying tissue.

2.3. Simulating Radiotherapy Using the Linear Quadratic Model

This paper investigates an extension of the preceding proportional glioma model to include the effects of the radiotherapy using the Linear Quadratic Model. According to the model, the probability of cells surviving following a single dose of radiation was observed to follow this relationship [47]:

where linear (Gy⁻¹) and quadratic (Gy⁻²) are the radiobiology parameters, which are interpreted biologically as repairable single and lethal double-strand breaks to the cell’s DNA, respectively [37]. In our case, where the clinical datasets have been extracted from patients with radiation in a number of fractions with the same dose, (9) turns to the following equation [48]:

where is the dose per fraction and is the total dose for fractions. In general, fast-proliferating cells, like GBM, have a high tissue response . In our experiment, we use a constant ratio that is widely used in clinical applications for highly developing cancer and has been extracted from empirical dose response data [36]. Thus, if (4b) is used for approximating cell proliferation, the overall proliferation term for (3) turns out to be the following:

2.4. Clinical Datasets

The model uses T1- and T2-MRIs (255 × 255 pixel slices) taken from six patients diagnosed with malignant glioblastoma multiforme and cured with radiotherapy. The datasets were provided by Saarland University, Germany, for the needs of the ContraCancrum project [49]. For each dataset there is information provided for the dose and the fractions. Moreover, there are two or more sessions taken on different dates for tracking glioma development and evaluation simulation result.

2.5. The Workflow of Simulation

The collected MRI datasets have different dimensions from the SRI-atlas. Thus, they have to be firstly registered and interpolated to the atlas, before applying the DRE equation. The Optimized Automatic Registration 3D approach by MIPAV is used for registering the clinical datasets to the atlas [50].

After registration, the datasets are delineated by expert radiologists from the University of Saarland, by using DoctorEye. DoctorEye is a platform for annotating, segmenting and visualizing MRI slices [51]. Figure 1(a) presents the 99th slice of one patient dataset after registration and interpolation to the atlas. Figure 1(b) presents the 99th slice of the SRI24 normal brain atlas. Respectively, Figures 1(d), 1(e), and 1(f) present the 99th image of white matter proportion, gray matter proportion and the one dominant eigenvector. Finally, Figure 1(c) presents the 99th slice of the MRI set with the tumor area annotated by the clinician.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1: (a) the 99th slice of the patient dataset after registration and interpolation, (b) the 99th slice of normal brain atlas, (c) the 99th image after skull stripping and segmentation of tumor, (d) the 99th image of white matter, (e) gray matter proportion, and (f) the one dominant eigenvector images extracted from SRI24 atlas.

In order to solve the DRE, the mathematical framework of Finite Differences that has been presented in [52] is used. More specifically, the Crank Nikolson numerical scheme has been implemented for solving the DRE in three dimensions. Parameters , , and initial concentration of cells , have been extracted from bibliography according to [41, 44]. Thus, we used day⁻¹,  mm²/day,  mm²/day, cells/mm³. For the solver of finite differences, we used a grid with dimensions  mm and a stable time step of  h.

3. Results and Discussion

3.1. Radiobiology Parameters

Before applying the diffusive model, it is appropriate to estimate the values of the radiobiology parameters and . The method used is that of Rockne et al. [36], by projecting the temporal tumor curves for different values of on the real datasets. Estimating radius is ideal for spherical tumors however GBM has very complex shape. Thus, instead of tumor radii, the estimated total number of GBM cells is used for identifying the value of that minimizes the error between simulated and actual target total GBM cells.

Figure 2 presents a set of simulated curves extracted for one specific patient after applying 7 different values of . The real tumor total cells are marked on the curve for three different MRI provided (the three green points on day 1, 15, and 70). The patient was treated with a total dose of 60 Gy at 2 Gy/fraction and 5 fractions/week for 6 weeks. Treatment starts 21 days after the first MRI and finishes on day 56. Each different curve is extracted by simulating radiotherapy for 0.003 increments of , starting from zero (no effect) and reaching . One can see that the optimal value for , that minimizes the error from the total number of cells is . The same process is followed for approximating the radiobiology parameters for the rest five patients. Table 1 presents the patient specific parameters for all the patients.

Table 1: Summary of model parameters for all six patients with glioblastoma multiforme included in the study.

Patient
Age

1
47
0.114
60
2
30

2
64
0.003
48
3
16

3
52
0.006
60
2
30

4
41
0.211
58
2
29

5
73
0.000
42
3
14

6
56
0.013
60
2
30

Figure 2: The curves depict the total number of GBM cells for one patient for eight values of ranging from 0 (no effect of radiotherapy) to 0.021 after applying the diffusive model. The same plot shows the three different points of real cell numbers, estimated after segmenting the tumors on three MRI datasets on the respective days. The optimal value of lies around 0.006 .

3.2. Simulation Result

After defining all simulation parameters, the diffusive model has been applied on the clinical datasets of MRIs. Figure 3 presents the simulated results for the 4th patient 21 days after the end of radiotherapy (100 days after diagnosis). The hot areas depict the resulting simulated tumor projected on the real MRI dataset. The real dataset has been extracted on the same day of simulation.

Figure 3: The simulated cell concentration (hot areas) projected on the real patient data on the 100th day after diagnosis and 21 days after the last fraction of radiotherapy.

For this specific patient, radiotherapy caused shrinkage of the tumor, which cannot be clearly seen on the simulation results of Figure 3. In order to have a better understanding of the shrinkage effect of radiotherapy, we use a line intersecting the point with the maximum concentration. Then, concentration is estimated along this line and presented on a graph. This is applied both on the real initial MRIs (day 1) and the simulated final result. Figure 4(a) shows the initial MRI with the intersection line projected, while Figure 4(b) presents the graph of the concentration of GBM along these line. The red curve shows the initial estimated concentration along this line, while the blue curve shows the simulated concentration of radiotherapy along the same line. In order to have a better visualization of the shrinkage effect, we also provide Figure 6 that depicts the change of the total number of cells in time (for 100 days). It is obvious that the tumor has shrunk considerably.

(a)

(b)

Figure 4: The graph presents the tumor cell concentration for the initial day (diagnosis day) and the 100th day after diagnosis (21th day after the last fraction of radiotherapy). These values have been calculated along the yellow line shown on the left MRI slice, which comes across the point with the maximum concentration value.

3.3. Evaluation of Model

We have made simulations on a series of six datasets in order to compare the simulated tumor with the final actual tumor, as annotated by the radiologists. For evaluation, we used a scheme with solid metrics and objective comparison. Therefore, the annotated final tumor serves as a golden ground truth and the Jaccard (JC), Dice (DS) and Volume Similarity (VS) metrics are adopted [53].

The evaluation metrics are defined as:

where TP (true positive) is the number of tumor voxels belonging to both ground truth and simulated result, FP (false positive) is the number of tumor voxels belonging to simulated result and not belonging to ground truth, and FN (false negative) is the number of tumor voxels belonging to ground truth but not belonging to simulated tumor.

Figure 5 depicts the scatter plots of the three metrics for each patient. The average resulting JC, DS, and VS values for all patients are 95.21%, 98.09% and 99.62%. These values are 2.1%, 4.01%, and 1.89% higher than the respective results for discrete diffusion coefficients and 6.45%, 8.49%, and 4.07% than the uniform tumor growth model.

(a)

(b)

(c)

Figure 5: Scatter plot of JC, DS, and VS (%) for six cases where the proportional, the discrete, and the uniform growth models have been applied.

Figure 6: The temporal change of the number of glioma cells for the 4th patient for 100 days. The last radiotherapy session was carried on day 79.

4. Conclusions

Despite the limited number of datasets, the evaluation results indicate that there might be a slight improvement in using the proposed model on glioblastoma multiforme. This is expected due to the several improvements introduced.

By using the proportional diffusion coefficients, all available information is fully exploited for characterizing tissues, without needing to truncate tissue areas to either white or gray matter.

Moreover, the use of atlases contributes to defining the matter proportions on the areas underlying the tumor area. This could not be true without the use of atlas, as the tumor hides the tissue matter needed for defining the diffusion coefficients at the tumor area. Also, by using anisotropy only on the portions of white matter, instead of all type of tissue, we better approximate the glioma migration along white fibers.

Radiotherapy is usually applied in fractions and the linear quadratic model has been adjusted to this, contrary to previous diffusive models which use the quadratic model for a single dose. By using the formalism of (10), the diffusive model incorporates this idea.

Acknowledgments

The present work was partially supported by the community initiative program INTERREG III, project “ΥΠEPΘEN,” financed by the EC through the European Regional Development Fund (ERDF) and by national funds of Greece and Cyprus. This work was also supported in part by the EC under the project TUMOR (FP7-ICT-2009.5.4-247754).

References

WHO, World Health Organization, 2010, http://www.who.int/mediacentre/factsheets/fs297

P. Boyle and B. Levin, “World cancer report 2008,” IARC, 2010

Brain Tumors, American Brain Tumor Association, http://www.abta.org/

Central Brain Tumor Registry of the United States, “Statistical report: primary brain tumors in the United States, 1992–1997,” Tech. Rep., Chicago, Ill, USA, 2000.

M. Tovi, “MR imaging in cerebral gliomas analysis of tumour tissue components.,” Acta Radiologica, Supplement, vol. 384, pp. 1–24, 1993.

S. J. Price, N. G. Burnet, T. Donovan et al., “Diffusion tensor imaging of brain tumours at 3 T: a potential tool for assessing white matter tract invasion?” Clinical Radiology, vol. 58, no. 6, pp. 455–462, 2003.

P. Lissoni, S. Meregalli, L. Nosetto et al., “Increased survival time in brain glioblastomas by a radioneuroendocrine strategy with radiotherapy plus melatonin compared to radiotherapy alone,” Oncology, vol. 53, no. 1, pp. 43–46, 1996.

E. G. van Meir, C. G. Hadjipanayis, A. D. Norden, H. K. Shu, P. Y. Wen, and J. J. Olson, “Exciting new advances in neuro-oncology: the avenue to a cure for malignant glioma,” CA Cancer Journal for Clinicians, vol. 60, no. 3, pp. 166–193, 2010.

D. Krex, B. Klink, C. Hartmann et al., “Long-term survival with glioblastoma multiforme,” Brain, vol. 130, no. 10, pp. 2596–2606, 2007.

University of California and Los Angeles Neuro-Oncology, “How our patients perform: Glioblastoma Multiforme [GBM],” 2010, http://www.neurooncology.ucla.edu/

K. R. Swanson, E. C. Alvord, and J. D. Murray, “Virtual brain tumours (gliomas) enhance the reality of medical imaging and highlight inadequacies of current therapy,” British Journal of Cancer, vol. 86, no. 1, pp. 14–18, 2002.

D. L. Silbergeld, R. C. Rostomily, and E. C. Alvord, “The cause of death in patients with glioblastoma is multifactorial: clinical factors and autopsy findings in 117 cases of supratentorial glioblastoma in adults,” Journal of Neuro-Oncology, vol. 10, no. 2, pp. 179–185, 1991.

A. Claes, A. J. Idema, and P. Wesseling, “Diffuse glioma growth: a guerilla war,” Acta Neuropathologica, vol. 114, no. 5, pp. 443–458, 2007.

S. A. Toms, W. C. Lin, R. J. Weil, M. D. Johnson, E. D. Jansen, and A. Mahadevan-Jansen, “Intraoperative optical spectroscopy identifies infiltrating glioma margins with high sensitivity,” Neurosurgery, vol. 57, no. 4, pp. 382–391, 2005.

F. G. Blankenberg, R. L. Teplitz, W. Ellis et al., “The influence of volumetric tumor doubling time, DNA ploidy, and histologic grade on the survival of patients with intracranial astrocytomas,” American Journal of Neuroradiology, vol. 16, no. 5, pp. 1001–1012, 1995.

H. B. Frieboes, X. Zheng, C. H. Sun, B. Tromberg, R. Gatenby, and V. Cristini, “An integrated computational/experimental model of tumor invasion,” Cancer Research, vol. 66, no. 3, pp. 1597–1604, 2006.

S. Sjöström, U. Andersson, Y. Liu et al., “Genetic variations in EGF and EGFR and glioblastoma outcome,” Neuro-Oncology, vol. 12, no. 8, pp. 815–821, 2010.

V. Cristini, J. Lowengrub, and Q. Nie, “Nonlinear simulation of tumor growth,” Journal of Mathematical Biology, vol. 46, no. 3, pp. 191–224, 2003.

T. S. Deisboeck, M. E. Berens, A. R. Kansal, S. Torquato, A. O. Stemmer-Rachamimov, and E. A. Chiocca, “Pattern of self-organization in tumour systems: complex growth dynamics in a novel brain tumour spheroid model,” Cell Proliferation, vol. 34, no. 2, pp. 115–134, 2001.

D. Drasdo and S. Höhme, “A single-cell-based model of tumor growth in vitro: monolayers and spheroids,” Physical Biology, vol. 2, no. 3, pp. 133–147, 2005.

P. Tracqui, “From passive diffusion to active cellular migration in mathematical models of tumour invasion,” Acta Biotheoretica, vol. 43, no. 4, pp. 443–464, 1995.

K. R. Swanson, E. C. Alvord, and J. D. Murray, “A quantitative model for differential motility of gliomas in grey and white matter,” Cell Proliferation, vol. 33, no. 5, pp. 317–329, 2000.

S. Jbabdi, E. Mandonnet, H. Duffau et al., “Simulation of anisotropic growth of low-grade gliomas using diffusion tensor imaging,” Magnetic Resonance in Medicine, vol. 54, no. 3, pp. 616–624, 2005.

G. S. Stamatakos, E. A. Kolokotroni, D. D. Dionysiou, E. C. Georgiadi, and C. Desmedt, “An advanced discrete state-discrete event multiscale simulation model of the response of a solid tumor to chemotherapy: mimicking a clinical study,” Journal of Theoretical Biology, vol. 266, no. 1, pp. 124–139, 2010.

O. Clatz, M. Sermesant, P. Y. Bondiau et al., “Realistic simulation of the 3-D growth of brain tumors in MR images coupling diffusion with biomechanical deformation,” IEEE Transactions on Medical Imaging, vol. 24, no. 10, pp. 1334–1346, 2005.

E. Konukoglu, O. Clatz, B. H. Menze et al., “Image guided personalization of reaction-diffusion type tumor growth models using modified anisotropic eikonal equations,” IEEE Transactions on Medical Imaging, vol. 29, no. 1, pp. 77–95, 2010.

C. Hogea, C. Davatzikos, and G. Biros, “Modeling glioma growth and mass effect in 3D MR images of the brain,” in Proceedings of the International Conference on Medical Image Computing and Computer Assisted Intervenion (MICCAI '07), vol. 4791, pp. 642–650, 2007.

A. R. Kansal, S. Torquato, G. R. Harsh, E. A. Chiocca, and T. S. Deisboeck, “Simulated brain tumor growth dynamics using a three-dimensional cellular automaton,” Journal of Theoretical Biology, vol. 203, no. 4, pp. 367–382, 2000.

A. Zizzari, Methods on tumor recognition and planning target prediction for the radiotherapy of cancer [Ph.D. thesis], Electronics Department, University of Magdeburg, 2004.

M. Tabatabai, D. K. Williams, and Z. Bursac, “Hyperbolastic growth models: theory and application,” Theoretical Biology and Medical Modelling, vol. 2, no. 1, p. 14, 2005.

A. Roniotis, K. Marias, V. Sakkalis, and M. Zervakis, “Diffusive modelling of glioma evolution: a review,” Journal of Biomedical Science and Engineering, vol. 3, no. 5, pp. 501–508, 2010.

A. Giese and M. Westphal, “Glioma invasion in the central nervous system,” Neurosurgery, vol. 39, no. 2, pp. 235–252, 1996.

A. Giese, R. Bjerkvig, M. E. Berens, and M. Westphal, “Cost of migration: invasion of malignant gliomas and implications for treatment,” Journal of Clinical Oncology, vol. 21, no. 8, pp. 1624–1636, 2003.

T. E. Conturo, N. F. Lori, T. S. Cull et al., “Tracking neuronal fiber pathways in the living human brain,” Proceedings of the National Academy of Sciences of the United States of America, vol. 96, no. 18, pp. 10422–10427, 1999.

K. R. Swanson, R. C. Rostomily, and E. C. Alvord, “A mathematical modelling tool for predicting survival of individual patients following resection of glioblastoma: a proof of principle,” British Journal of Cancer, vol. 98, no. 1, pp. 113–119, 2008.

R. Rockne, J. K. Rockhill, M. Mrugala et al., “Predicting the efficacy of radiotherapy in individual glioblastoma patients in vivo: a mathematical modeling approach,” Physics in Medicine and Biology, vol. 55, no. 12, pp. 3271–3285, 2010.

E. Hall, Radiobiology for the Radiologist, J. B. Lippincott, Philadelphia, Pa, USA, 1994.

G. S. Stamatakos, V. P. Antipas, N. K. Uzunoglu, and R. G. Dale, “A four-dimensional computer simulation model of the in vivo response to radiotherapy of glioblastoma multiforme: studies on the effect of clonogenic cell density,” British Journal of Radiology, vol. 79, no. 941, pp. 389–400, 2006.

V. P. Antipas, G. S. Stamatakos, and N. K. Uzunoglu, “A patient-specific in vivo tumor and normal tissue model for prediction of the response to radiotherapy: a computer simulation approach,” Methods of Information in Medicine, vol. 46, no. 3, pp. 367–375, 2007.

A. Roniotis, G. Manikis, V. Sakkalis, M. E. Zervakis, I. Karatzanis, and K. Marias, “High grade glioma diffusive modeling using statistical tissue information and diffusion tensors extracted from atlases,” IEEE Transactions on Information Technology in Biomedicine, vol. 16, no. 2, pp. 255–263, 2011.

A. Roniotis, V. Sakkalis, I. Karatzanis, M. E. Zervakis, and K. Marias, “In-depth analysis and evaluation of diffusive glioma models,” IEEE Transactions on Information Technology in Biomedicine, vol. 16, no. 3, pp. 299–307, 2012.

V. Sakkalis, A. Roniotis, C. Farmaki, I. Karatzanis, and K. Marias, in Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, vol. 1, pp. 6809–6812.

S. Wahl, T. D. C. Foletto, and G. Feldmann, “A mathematical model for the estimation of treatment cost in cancer radiotherapy,” in Proceedings of the 3rd Southern Conference on Computational Modeling (MCSUL '09), pp. 77–81, November 2009.

H. L. P. Harpold, E. C. Alvord, and K. R. Swanson, “The evolution of mathematical modeling of glioma proliferation and invasion,” Journal of Neuropathology and Experimental Neurology, vol. 66, no. 1, pp. 1–9, 2007.

T. Rohlfing, N. M. Zahr, E. V. Sullivan, and A. Pfefferbaum, “The SRI24 multichannel atlas of normal adult human brain structure,” Human Brain Mapping, vol. 31, no. 5, pp. 798–819, 2010.

http://www.nitrc.org/projects/sri24

G. Steel, Basic Clinical Radiobiology, Arnold, London, UK, 3rd edition, 2002.

M. Partridge, M. Ramos, A. Sardaro, and M. Brada, “Dose escalation for non-small cell lung cancer: analysis and modelling of published literature,” Radiotherapy and Oncology, vol. 99, no. 1, pp. 6–11, 2011.

K. Marias, D. Dionysiou, V. Sakkalis et al., “Clinically driven design of multi-scale cancer models: the ContraCancrum project paradigm,” Interface Focus, vol. 1, no. 3, pp. 450–461, 2011.

http://mipav.cit.nih.gov/

E. Skounakis, C. Farmaki, V. Sakkalis et al., “DoctorEye: a clinically driven multifunctional platform, for accurate processing of tumors in medical images,” The Open Medical Informatics Journal, vol. 4, pp. 105–115, 2010.

A. Roniotis, K. Marias, V. Sakkalis, G. D. Tsibidis, and M. Zervakis, “A complete mathematical study of a 3D model of heterogeneous and anisotropic glioma evolution,” in Proceedings of the 31st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '09), pp. 2807–2810, September 2009.

J. K. Udupa, V. R. LeBlanc, H. Schmidt et al., “A methodology for evaluating image segmentation algorithms,” in Medical Imaging 2002: Image Processing, pp. 266–277, February 2002.