Adult gliomas are aggressive brain tumours associated with low patient survival rates and limited life expectancy. The most important hallmark of this type of tumour is its invasive behaviour, characterized by a markedly phenotypic plasticity, infiltrative tumour morphologies and the ability of malignant progression from low- to high-grade tumour types. Indeed, the widespread infiltration of healthy brain tissue by glioma cells is largely responsible for poor prognosis and the difficulty of finding curative therapies. Meanwhile, mathematical models have been established to analyse potential mechanisms of glioma invasion. In this review, we start with a brief introduction to current biological knowledge about glioma invasion, and then critically review and highlight future challenges for mathematical models of glioma invasion.

Standard radiobiology theory of radiation response assumes a uniform innate radiosensitivity of tumors. However, experimental data show that there is significant intratumoral heterogeneity of radiosensitivity. Therefore, a model with heterogeneity was developed and tested using existing experimental data to show the potential effects from the presence of an intratumoral distribution of radiosensitivity on radiation therapy response over a protracted radiation therapy treatment course. The standard radiation response curve was modified to account for a distribution of radiosensitivity, and for variations in the repopulation rates of the tumor cell subpopulations. Experimental data from the literature were incorporated to determine the boundaries of the model. The proposed model was then used to show the changes in radiosensitivity of the tumor during treatment, and the effects of fraction size, α/β ratio and variation of the repopulation rates of tumor cells. In the presence of an intratumoral distribution of radiosensitivity, there is rapid selection of radiation-resistant cells over a course of fractionated radiation therapy. Standard treatment fractionation regimes result in the near-complete replacement of the initial population of sensitive cells with a population of more resistant cells. Further, as treatment progresses, the tumor becomes more resistant to further radiation treatment, making each fractional dose less efficacious. A wider initial distribution induces increased radiation resistance. Hypofractionation is more efficient in a heterogeneous tumor, with increased cell kill for biologically equivalent doses, while inducing less resistance. The model also shows that a higher growth rate in resistant cells can account for the accelerated repopulation that is seen during the clinical treatment of patients. Modeling of tumor cell survival with radiosensitivity heterogeneity alters the predicted tumor response, and explains the induction of radiation resistance by radiation treatment, the development of accelerated repopulation, and the potential beneficial effects of hypofractionation. Tumor response to treatment may be better predicted by assaying for the distribution of radiosensitivity, or the extreme of the radiosensitivity, rather than measuring the initial, general radiation sensitivity of the untreated tumor.

Tumor heterogeneity is widely considered to be a determinant factor in tumor progression and in particular in its recurrence after therapy. Unfortunately, current medical techniques are unable to deduce clinically relevant information about tumor heterogeneity by means of non-invasive methods. As a consequence, when radiotherapy is used as a treatment of choice, radiation dosimetries are prescribed under the assumption that the malignancy targeted is of a homogeneous nature. In this work we discuss the effects of different radiation dose distributions on heterogeneous tumors by means of an individual cell-based model. To that end, a case is considered where two tumor cell phenotypes are present, which we assume to strongly differ in their respective cell cycle duration and radiosensitivity properties. We show herein that, as a result of such differences, the spatial distribution of the corresponding phenotypes, whence the resulting tumor heterogeneity can be predicted as growth proceeds. In particular, we show that if we start from a situation where a majority of ordinary cancer cells (CCs) and a minority of cancer stem cells (CSCs) are randomly distributed, and we assume that the length of CSC cycle is significantly longer than that of CCs, then CSCs become concentrated at an inner region as tumor grows. As a consequence we obtain that if CSCs are assumed to be more resistant to radiation than CCs, heterogeneous dosimetries can be selected to enhance tumor control by boosting radiation in the region occupied by the more radioresistant tumor cell phenotype. It is also shown that, when compared with homogeneous dose distributions as those being currently delivered in clinical practice, such heterogeneous radiation dosimetries fare always better than their homogeneous counterparts. Finally, limitations to our assumptions and their resulting clinical implications will be discussed.

The choice of any radiotherapy treatment plan is usually made after the evaluation of a few preliminary isodose distributions obtained from different beam configurations. Despite considerable advances in planning techniques, such final decision remains a challenging task that would greatly benefit from efficient and reliable assessment tools. For any dosimetric plan considered, data on dose-volume histograms supplied by treatment planning systems are used to provide estimates on planning target coverage as well as on sparing of organs at risk and the remaining healthy tissue. These partial metrics are then combined into a dose distribution index (DDI), which provides a unified, easy-to-read score for each competing radiotherapy plan. To assess the performance of the proposed scoring system, DDI figures for fifty brain cancer patients were retrospectively evaluated. Patients were divided in three groups depending on tumor location and malignancy. For each patient, three tentative plans were designed and recorded during planning, one of which was eventually selected for treatment. We thus were able to compare the plans with better DDI scores and those actually delivered. When planning target coverage and organs at risk sparing are considered as equally important, the tentative plan with the highest DDI score is shown to coincide with that actually delivered in 32 of the 50 patients considered. In 15 (respectively 3) of the remaining 18 cases, the plan with highest DDI value still coincides with that actually selected, provided that organs at risk sparing is given higher priority (respectively, lower priority) than target coverage. DDI provides a straightforward and non-subjective tool for dosimetric comparison of tentative radiotherapy plans. In particular, DDI readily quantifies differences among competing plans with similar-looking dose-volume histograms and can be easily implemented for any tumor type and localization, irrespective of the planning system and irradiation technique considered. Moreover, DDI permits to estimate the dosimetry impact of different priorities being assigned to sparing of organs at risk or to better target coverage.