Determining the optimum tumor control probability model in radiotherapy of glioblastoma multiforme using magnetic resonance imaging data pre and post radiation therapy
Shabnam Banisharif^{1}, Daryoush ShahbaziGahrouei^{1}, Ali Akhavan^{2}, Naser Rasouli^{1}, Saghar ShahbaziGahrouei^{3}
^{1} Department of Medical Physics, School of Medicine, Isfahan University of Medical Sciences, Isfahan, Iran ^{2} Department of Radiation Oncology, School of Medicine, Isfahan University of Medical Sciences, Isfahan, Iran ^{3} MD Student, School of Medicine, Isfahan University of Medical Sciences, Isfahan, Iran
Date of Submission  10Oct2020 
Date of Decision  14Apr2021 
Date of Acceptance  30Aug2021 
Date of Web Publication  18Feb2022 
Correspondence Address: Prof. Daryoush ShahbaziGahrouei Department of Medical Physics and Medical Engineering, Isfahan University of Medical Sciences, Isfahan Iran
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/jrms.JRMS_1138_20
Background: Glioblastoma multiforme (GBM) is the most common and malignant brain tumor. The current standard of care is surgery followed by radiation therapy (RT). Radiotherapy treatment plan evaluation relies on radiobiological models for accurate estimation of tumor control probability (TCP). This study aimed to assess the impact of obtained magnetic resonance imaging (MRI) data before and 12 weeks after RT to achieve the optimum TCP model to improve dose prescriptions in radiation therapy of GBM. Materials and Methods:: In this quasiexperimental study, MR images and its relevant data from 30 patients consisting of 9 females and 21 males (mean age of 46.3 ± 15.8 years) diagnosed with GBM, whose referred for radiotherapy were selected. The data of age, gender, tumor size, volume, and signal intensity using analysis of MRI data pre and postradiotherapy were used for calculating TCP. TCP was calculated from three common radiobiological models including Poisson, linear quadratic, and equivalent uniform dose. The impact of some radiobiological parameters on final TCP in all patients planned with threedimensional conformal radiation therapy was obtained. Results: A statistically significant difference was found among TCP in Poisson model compared to the other two models (P < 0.001). Changes in tumor volume and size after treatment were statistically significant (P < 0.05). Different combinations of radiobiological parameters (α/β and SF_{2} in all models) observed were meaningful (P < 0.05). Conclusion: The results showed that among TCP radiobiological models, the optimum is the Poisson. The results also identified the importance of TCP radiobiological models in order to improve radiotherapy dose prescriptions.
Keywords: Cancer, glioblastoma multiforme, magnetic resonance imaging, radiation therapy, radiobiological models
How to cite this article: Banisharif S, ShahbaziGahrouei D, Akhavan A, Rasouli N, ShahbaziGahrouei S. Determining the optimum tumor control probability model in radiotherapy of glioblastoma multiforme using magnetic resonance imaging data pre and post radiation therapy. J Res Med Sci 2022;27:10 
How to cite this URL: Banisharif S, ShahbaziGahrouei D, Akhavan A, Rasouli N, ShahbaziGahrouei S. Determining the optimum tumor control probability model in radiotherapy of glioblastoma multiforme using magnetic resonance imaging data pre and post radiation therapy. J Res Med Sci [serial online] 2022 [cited 2022 Sep 24];27:10. Available from: https://www.jmsjournal.net/text.asp?2022/27/1/10/337907 
Introduction   
Glioblastoma multiforme (GBM) can originate from the central nervous system and no metastases have been noticed, but it is wildly expanded to parenchyma around the brain. The most common treatment for GBM was established as surgery with maximum safe resection followed by concurrent chemotherapy and conformal radiation therapy.^{[1]}
Magnetic resonance imaging (MRI) is the best standard method for diagnosing GBM, which is widely utilized in the classification of patients.^{[2]} MRI techniques include T_{1}weighted imaging (T_{1}W), T_{1} contrastenhancing (T_{1}CE), T_{2}weighted imaging (T_{2}W), and fluidattenuated inversion recovery (FLAIR) sequences. Furthermore, advanced imaging techniques can provide more information about the size and position of tumors.^{[2]}
Radiation therapy (RT) has a significant effect on the treatment of brain tumors. The goal of RT is to find the best treatment plan for high local tumor control with minimal side effects. The main goal of RT is to deliver 95% of prescribed dose to maximum volume of the tumor and safety of organ at risk.^{[3]} Mathematical models have been expanded to quantify the biological effects of radiation and predict treatment outcome known as concepts tumor control probability (TCP) and normal tissue complication probability (NTCP).^{[4]} TCP modeling plays an important role in assessing the impact of RT treatment strategies.^{[5]} Three of the most common TCP radiobiological models are Poisson, equivalent uniform dose (EUD), and linear quadratic (LQ). Using mentioned radiobiological models, TCP was calculated for threedimensional conformal radiation therapy (3DCRT), and the relative impact of the radiobiological parameters for patients treated with radiation was also evaluated. Then, treatment planning was modified by using the results of TCP radiobiological models and MRI.
Since there is no published results in conventional MRI to distinguish between true progression and pseudoprogression, the Response Assessment in NeuroOncology (RANO) criteria are used for highgrade glioma to define progressive disease according to the radiographic changes.^{[6]}
Many research studies have examined almost similar data from different perspective studies and reported different conclusions.^{[4],[7],[8]} For instance, Majós et al. have looked at MR images after surgery and shortly pre and 6–12 weeks postRT.^{[7]} They found that performing MRI before RT has positive impacts on the management of patients with GBM by reducing the ratio of pseudoprogression assessments and providing additional predictive information. In another study, Ahmad et al.^{[4]} examined the impact of radiobiological parameters on TCP for prostate cancer. They have observed that (surviving fraction) SF_{2} is the dominant predictor to radiation response. Consequently, other researchers assessed the impact of radiation dose escalation from 59.6 to 90 Gy on TCP and NTCP in 10 patients planned with 3DCRT, but they have calculated only TCP using radiobiological models of EUD and LQ.^{[8]}
All previous relevant studies have evaluated image parameters such as signal intensity, tumor dimensions, and size just pre and posttreatment for highgrade glioma in the region of interest (ROI), but none of them have investigated the association of these changes with TCP in order to introduce the appropriate radiobiological models. For this reason, this study aimed to assess the impact of obtained MRI data before and 12 weeks after RT to achieve the optimum TCP model to improve dose prescriptions in radiation therapy of GBM.
Methods   
Study design and participants
In this quasiexperimental study, 30 patients consisting of 9 females and 21 males (mean age of 46.3 ± 15.8 years) having Grade IV glioma GBM according to the World Health Organization (WHO)^{[6]} were included. Patients with fear of the indoor or intolerance MRI and insufficient information were excluded from the study. Pre and postoperative MRI studies were carried out on all studied patients. The study was approved by the university ethical committee (ethical code no. IR.MUI.MED.REC.1399.167).
Referred patients to SeyedAlShohada Hospital, Isfahan, Iran, were scanned with a 1.5 T MR scanner (Siemens Healthcare, Erlangen, Germany) using commercially quadrature head coils before RT and after surgical resection.^{[9]} Approximately 12 weeks after RT complication, MRI process was performed again with similar protocols. Performing MRI 3 months after RT reduces the chance of pseudoprogression and determines tumor residual by evaluating surgical site.^{[6]} The standard protocol consists of sequences of axial T_{1}W and T_{2}W images, fast spin echoaxial (FLAIR) images, and coronal T_{1}W images. Contrast agent of gadoliniumdiethylenetriaminepentaacetic acid was IV injected at standard dose (0.1 mmol/kg) and T_{1}W images were repeated.^{[10]}
3DCRT was performed with Varian treatment planning system (TPS). All patients underwent computed tomography (CT) with bed similar to treatment bed before RT. CT and MR images of patients were used as input data in TPS. The prescription dose of 60 Gy was delivered during 30 fractions with a range of 1.8–2.0 Gy each session.^{[11]}
Gross tumor volume (GTV) is specified with regard to T_{1} contrast enhancement region with a 20–30 mm extra margin, and clinical target volume is obtained by adding a 1.5–2 cm margin to GTV.^{[12]}
To evaluate the response of radiobiological models, cumulative dose–volume histogram (DVH) of calculated treatment plans was extracted from the TPS. Software package (BIOPLAN) was used to calculate TCP. DVH was given as software input data. The MATLAB (The MathWorks Inc., Natick, USA) codes were utilized in all of the models to calculate the TCP.
Three radiobiological models for tumor control are considered: EUD, Poisson, and LQ.
The equivalent uniform dose model
In EUD, if two different distributions doses have the same effect, the number of their clonogenic cells be equal and TCP is calculated by the following formula:^{[8]}
Where the tumor control dose (TCD_{50}) is the tumor dose for controlling 50% of the tumors when the tumor is homogeneously irradiated and γ50 is the normalized dose–response gradient at the dose of 50 Gy; in the second formula, vi is the fractional organ volume receiving a dose Di and a is a tissuespecific parameter that describes the volume effect.
The Poisson model
To predict the random killing of clonogenic cells (the Poisson standard model) by radiation, TCP is calculated using the following equation:^{[13]}
Where D_{50} is the dose at which 50% of tumors are controlled. The number of clonogenic cells is predicted by exp(eγ) and where γ is called γvalue (normalized dose–response gradient) and may be defined as the percentage point change in response for a 1% change in dose.
The linearquadratic model
The concept of survival cells in the irradiated tissue with a total radiation dose of D = n × d (which n is the number of fractions and d is the dose per fraction) based on LQ model is presented as the following formula:^{[14]}
Radiobiological parameters for GBM included radiosensitivity (a): 0.12 Gy^{1}, repair capability (b): 0.015 Gy^{2}, α/β ratio: 8, repopulation doubling time (T_{d}): 15.4 days, the number of clonogeneses (N): 1 × 10^{4}, kickoff time (T_{K}): 27, D_{ref}: 2 Gy, and TCD_{50} which was chosen from similar Published work^{[15]}: 47.2 and used in the present study to determine TCP. The slope of the dose–response curve is calculated by the following formula:
Magnetic resonance imaging data acquisition
MR image ROI was manually drawn by the Oncologist around the tumor to separate nonenhancing lesion, contrastenhancing lesion, necrotic lesion, and normal white matter. The feasibility of tumor volume, area, and size (two perpendicular diameters and signal intensity measurement were evaluated using the routine clinic software (Digimizer and 3DDoctor) for image processing data of each patient were investigated. The volume was calculated using 3DDoctor software and TPS feasibility. Furthermore, two perpendicular diameters, area and signal intensity, were obtained using Digimizer software. According to the RANO criteria, shrinkage more than 50% in tumor size is considered for partial response, and 25% increase in the sum of enhancing target lesions means progressive disease, but the rest of the patients were not included in these two groups.^{[16]} The images were interpreted by expert oncologists and radiologists.
Statistical analysis
Data were analyzed using IBM SPSS (IBM Crop. 2019. IBM SPSS Statistics for Windows, version 26.0. NY, EUA). For quantitative data, the descriptive statistics were calculated for each category. It was presented as the range of observations (the lowest observation [minimum] minus the largest one [maximum]), the median value (Md), the interquartile range (IQR = the 3^{rd} quartile minus the 1^{st} quartile), the average value (mean), the standard deviation, and 95% confidence interval for the mean parameters. Initially, the Shapiro–Wilk test was used to test normality in each group. Despite the data violated the normality test, according to the central limit theorem, the parametric tests were used to compare the differences in the mean values between groups. Since the sphericity test was violated for all the tests, a repeated measures ANOVA with a Greenhouse–Geisser correction was conducted to compare the mean between more than two dependent groups, followed by the Bonferroni test, as the post hoc tests. For dependent groups (pre and postobservations), the paired samples ttest was utilized. Finally, the Spearman correlation was run for investigating the monotonic relationship between two continuous nonnormal variables. A level of 5% was considered to be significant in all the analyses.
Results   
From 30 studied patients (9 females and 21 males, mean age of 46.28 ± 15.8 years), after RT, 4 patients had progressive diseases, 5 patients showed partial response, and the rest of them showed stable conditions to treatment. These observations were based on the RANO criteria.^{[17]} [Figure 1] indicates that there existed no statistically significant differences in tumor volume, area, size, and the signal intensity of the tumor, before and after RT (P > 0.05). On the other hand, for each patient, the TCP values were calculated from three models (Poisson, EUD, and LQ models). Descriptive statistics of computed TCPs from Poisson model (based on three values of α/β) and EUD and LQ models (based on different values of SF_{2}) are presented in [Table 1]. In addition, in [Figure 2], the repeated measures ANOVA with a Greenhouse–Geisser correction determined that mean TCPs in Poisson model differed statistically significantly between minimum, middle, and maximum groups (F_{(1.1,31.9)} = 134.31, P < 0.001, effect size = 0.82). The same results were observed for the EUD model (F_{(1.3,36.2)} = 108.49, P < 0.001, effect size = 0.79) and LQ model (F_{(1.0,29.0)} = 550.04, P < 0.001, effect size = 0.95). Moreover, post hoc tests using the Bonferroni correction revealed that in Poisson model, in the maximum group (α/β = 10.8), the TCPs (15.07 ± 5.20) were significantly more than the minimum (11.24 ± 3.95) and maximum (11.36 ± 4.28) groups (P < 0.01). Therefore, increase in α parameter, followed by increase in tumor sensitivity and decrease in βcoefficient, increases the α/β ratio, and according to the findings, TCP was almost increasing [Figure 2]a. However, in EUD and LQ models, the computed TCPs in the maximum group (SF_{2} = 0.7) were statistically less than the other groups (P < 0.001). As well, the differences in TCPs between EUD and LQ models were investigated by independent samples ttest, individually for SF_{2} = 0.3, 0.4, 0.5, 0.6, and 0.7 Gy. The results in [Figure 3]a show that the obtained TCPs from EUD and LQ models were not significantly different, except for SF_{2} = 0.7 Gy (P < 0.001). Moreover, the computed TCP values in the middle groups were compared among the three models via ANOVA. The results are presented in [Figure 3]b. As expected, in the middle groups, the mean value of computed TCP from the Poisson model (11.36 ± 4.28) was significantly less than the TCPs from the EUD (96.91 ± 10.64) and LQ models (99.99 ± 0.01) (P < 0.001). Finally, the Spearman rankorder correlation was run to determine the relationship between computed TCPs in three models [Table 2]. From the results, the computed TCP values from the Poisson and EUD models were positively correlated (ρ = 0.37, P < 0.05). It can be concluded that increasing one leads to increasing the other and vice versa. However, no monotonic correlation was found in TCP values between Poisson and LQ models (ρ = 0.20, P > 0.05), and the same results were observed for TCP values between EUD and LQ models (r = 029, P > 0.05). On the other side, [Table 3] shows the descriptive statistics of EUD values in three experimental groups. From the oneway repeated measures ANOVA results, there existed no statistically significant differences in EUD values between low, medium, and high groups (F_{(1.0,29.4)} = 2.65, P > 0.05, effect size = 0.08).  Figure 1: The mean ± standard deviation of (a) tumor volume, (b) tumor area, (c) the signal intensity of tumor, and (d) tumor size (calculated by multiplying two diameters perpendicular to each other). Using paired samples ttest, NS = Not significant and n = 30
Click here to view 
 Figure 2: The mean value of TCP (±standard deviation) among minimum, middle, and maximum groups individually for (a) Poisson model, (b) EUD model, and (c) LQ model. Using oneway repeated measures ANOVA with a Greenhouse–Geisser correction, followed by Bonferroni test. In each type of model, the groups with a similar English letter are not statistically significantly different from each other, ***P < 0.001, and n = 30. EUD = Equivalent uniform dose, TCP = Tumor control probability; LQ = Linear quadratic
Click here to view 
 Figure 3: The mean value of computed TCP (±standard deviation). (a) Comparison of TCPs between EUD and LQ models, individually for different values of SF_{2} (0.3, 0.4, 0.5, 0.6, and 0.7 Gray), using paired samples ttest. (b) Comparison of TCPs in the middle group (Poisson model: α/β =8, and EUD and LQ models: SF_{2} = 0.5) among the three models. Using oneway repeated measures ANOVA with a Greenhouse–Geisser correction, followed by Bonferroni test; The groups with a similar English letter are not statistically significantly different from each other. NS = Not significant, ***P < 0.001, and n = 30. EUD = Equivalent uniform dose; TCP = Tumor control probability; LQ = Linear quadratic
Click here to view 
 Table 1: Descriptive statistics of tumor control probability values in three models (Poisson, equivalent uniform dose, and linear quadratic)
Click here to view 
 Table 2: The Spearman correlation in tumor control probability values between Poisson, equivalent uniform dose, and linearquadratic models in the middle groups
Click here to view 
 Table 3: Descriptive statistics of equivalent uniform dose values in three experimental studied groups
Click here to view 
Discussion   
Despite technical advances for patient treatment with GBM, nearly all patients die due to local progression of the tumor; however, improvement in local tumor control is considered to play an important role in patient management. Clinical outcome data have shown that four patterns including complete response, partial response, stable disease, and progression were defined for tumor growth that satisfactorily correlated with patient survival.^{[9],[15]} In this study, with RANO criteria, tumor progression after radiation therapy was found in 4 out of 30 studied patients, most of whom died in <15 months. Partial response was indicated in five patients, and in the remaining 21 patients, stable condition was seen. Treatment response, measured by changes in GTV observed on T_{1}W MR images, can be predicted and quantified based on relationship between model parameters and various measures of the external RT response.
The results demonstrated that volume, area, and two perpendicular diameters significantly reduced after treatment (P < 0.05), but signal intensity remained unchanged [Figure 1]. From point of comparison view among TCP in three studied models, the mean TCP over the patient in EUD was 88% more than Poisson and 3% more than LQ. In addition, in Poisson, it was 88.5% less than the other two models.
As shown in [Table 1], a statistical significance in TCP values between three investigated models (P < 00.1) was obtained. TCP in Poisson was observed statistically less than EUD and LQ. Moreover, TCP values were compared between EUD and LQ for different values of SF_{2} (0.3–0.7) using Mann–Whitney test. The result showed that TCP values in EUD are statistically less than that of LQ for SF_{2} (0.3–0.6), but in SF_{2} = 0.7, TCP values in LQ are more than EUD (P < 0.01). Furthermore, TCP in LQ for SF_{2} (0.30.5) shows almost the same value.
In this model, TCP value corresponds to α/β of 5–10.8. Assuming tumor clones with varying SF_{2} from 0.3 to 0.7 were used. In EUD, it was found that the mean TCP decreases with increasing SF_{2} (P < 0.001) [Figure 2]. If SF_{2} is constant, TCP in EUD model increases with increasing α/β ratio, but it is not very impressive. As LQ shown in [Figure 2], TCP was evaluated with respect to SF_{2} changes, which decreased with increasing SF_{2}. By comparing TCP in EUD and LQ for different SF_{2} values, it is found that TCP in EUD for SF_{2} (0.3–0.6) was lower than LQ, but in SF_{2} = 0.7, the LQ was lower than that in EUD.
The Spearman correlation analysis revealed a fair positive partial correlation between TCP values in Poisson and EUD (r = 0.37, P < 0.05), but the relation between LQ with Poisson and EUD (α/β =8 and SF_{2} = 0.5) showed no statistically significant difference.
EUD expresses cell survival with SF_{2} at dose of 2 Gy, which indicates the dose required to produce a specific effect on tissue and depends on the type of tumor, but it is independent of SF_{2}, γ50 and TCD_{50}. Thus, by changing (α) and α/β ratio, the EUD reduces whereas these changes are not statistically significant [Table 2]. Hence, the effect of α/β on EUD was much greater when compared with SF_{2.}
In studies performed on mathematical models with the aim of estimating radiobiological parameters,^{[14]} no specific value was observed for TCP, but according to their low survival probability, it is expected to be 30%. [Figure 3] shows that TCP in LQ and EUD was more than 90%, which is in conflict with the probability of low survival. Therefore, these two models cannot be considered suitable indicators for calculating TCP. In Poisson model, the average TCP was 11.36 on average, which is lower than the threshold introduced in the literature.^{[7],[8],[14]} The reasons for TCP reduction here consist of lack of dose coverage in PTV. Therefore, the low TCP value corresponding to the mentioned clinical data and low survival rate of the patients, which is proportional to low TCP value, confirming the Poisson model as an appropriate model.
Previous studies^{[7],[8]} focused on patients who had undergone surgical resection to compare their MRI. However, the present study included all patients for MRI pre and post treatment to adapt with radiobiological models. The differences between previous observations^{[8]} and the present study may arise from calculated TCP with the Poisson model to adapt with MRI outcome data. Treatment response is still a serious problem, and Rowe et al.'s study showed that MR images are not enough alone to recognize true progression in patients with GBM.^{[6]}
The survival of patients with glioma depends on ability to respond to treatment. Although previous studies have indicated that MR images pre and post radiotherapy may be an important factor in predicting survival for gliomas, the present study is the first to collect the data from pre and post treatment and their radiobiological models analyses.
Some limitations of this study are distribution of the set of patients into three possibilities of tumor growth (progressive disease, partial response, and stable) examination, which produced low numbers of patients in some groups and also lack of control study group.
Conclusion   
This study offers a combination of clinical observations and mathematical modeling to test the hypothesis that there is a relationship between image changes and TCP. The radiobiological parameter implications in the final TCP calculation were also illustrated. Therefore, improvement of mathematical models based on clinical results might be useful. The findings of this work showed that among the studied TCP radiobiological models using MRI data before and 12 weeks after RT of GBM, the optimum is the Poisson model. The importance of using radiobiological models to improve dose prescriptions in radiation therapy of GBM was another outcome. In the clinics, this study may be used to delineate radiation field and help to evaluate treatment response. Further studies could provide additional information about the clinical performance by working on diffusion and perfusionweighted MR images and their findings in relation to the TCP radiobiological models.
Financial support and sponsorship
This work was financially supported (grant no. 398936, with ethical code no. IR.MUI.MED.REC.1399.167) by the Isfahan University of Medical Sciences, Isfahan, Iran.
Conflicts of interest
There are no conflicts of interest.
References   
1.  Gao Y, Han F, Zhou Z, Cao M, Kaprealian T, Kamrava M, et al. Distortionfree diffusion MRI using an MRIguided TriCobalt 60 radiotherapy system: Sequence verification and preliminary clinical experience. Med Phys 2017;44:535766. 
2.  Shukla G, Alexander GS, Bakas S, Nikam R, Talekar K, Palmer JD, et al. Advanced magnetic resonance imaging in glioblastoma: A review. Chin Clin Oncol 2017;6:40. 
3.  MacDonald SM, Ahmad S, Kachris S, Vogds BJ, DeRouen M, Gitttleman AE, et al. Intensity modulated radiation therapy versus threedimensional conformal radiation therapy for the treatment of high grade glioma: A dosimetric comparison. J Appl Clin Med Phys 2007;8:4760. 
4.  Ahmad S, Vogds BJ, McKenna F, Vlachaki MT. Tumor control probability (TCP) in prostate cancer: Role of radiobiological parameters and radiation dose escalation. J Xray Sci Technol 2009;17:34754. 
5.  CasaresMagaz O, van der Heide UA, Rørvik J, Steenbergen P, Muren LP. A tumour control probability model for radiotherapy of prostate cancer using magnetic resonance imagingbased apparent diffusion coefficient maps. Radiother Oncol 2016;119:1116. 
6.  Rowe LS, Butman JA, Mackey M, Shih JH, CooleyZgela T, Ning H, et al. Differentiating pseudoprogression from true progression: Analysis of radiographic, biologic, and clinical clues in GBM. J Neurooncol 2018;139:14552. 
7.  Majós C, Cos M, Castañer S, Pons A, Gil M, FernándezCoello A, et al. Preradiotherapy MR imaging: A Prospective pilot study of the usefulness of performing an MR examination shortly before radiation therapy in patients with glioblastoma. AJNR Am J Neuroradiol 2016;37:222430. 
8.  De La Fuente Herman T, Ahmad And S, Vlachaki MT. Intensity modulated radiation therapy versus three dimensional conformal radiation therapy for treatment of high grade glioma: A radiobiological modeling study. J Xray Sci Technol 2010;18:393402. 
9.  Saraswathy S, Crawford FW, Lamborn KR, Pirzkall A, Chang S, Cha S, et al. Evaluation of MR markers that predict survival in patients with newly diagnosed GBM prior to adjuvant therapy. J Neurooncol 2009;91:6981. 
10.  Young RJ, Gupta A, Shah AD, Graber JJ, Schweitzer AD, Prager A, et al. Potential role of preoperative conventional MRI including diffusion measurements in assessing epidermal growth factor receptor gene amplification status in patients with glioblastoma. AJNR Am J Neuroradiol 2013;34:22717. 
11.  Kovács Á Emri M, Opposits G, Pisák T, Vandulek C, Glavák C, et al. Changes in functional MRI signals after 3D based radiotherapy of glioblastoma multiforme. J Neurooncol 2015;125:15766. 
12.  Dhermain F. Radiotherapy of highgrade gliomas: Current standards and new concepts, innovations in imaging and radiotherapy, and new therapeutic approaches. Chin J Cancer 2014;33:1624. 
13.  Allen Li X, Alber M, Deasy JO, Jackson A, Ken Jee KW, Marks LB, et al. The use and QA of biologically related models for treatment planning: Short report of the TG166 of the therapy physics committee of the AAPM. Med Phys 2012;39:1386409. 
14.  Zaider M, Minerbo GN. Tumour control probability: A formulation applicable to any temporal protocol of dose delivery. Phys Med Biol 2000;45:27993. 
15.  Okunieff P, Morgan D, Niemierko A, Suit HD. Radiation doseresponse of human tumors. Int J Radiat Oncol Biol Phys 1995;32:122737. 
16.  Okada H, Weller M, Huang R, Finocchiaro G, Gilbert MR, Wick W, et al. Immunotherapy response assessment in neurooncology: A report of the RANO working group. Lancet Oncol 2015;16:e53442. 
17.  Lin NU, Lee EQ, Aoyama H, Barani IJ, Baumert BG, Brown PD, et al. Challenges relating to solid tumour brain metastases in clinical trials, part 1: Patient population, response, and progression. A report from the RANO group. Lancet Oncol 2013;14:e396406. 
[Figure 1], [Figure 2], [Figure 3]
[Table 1], [Table 2], [Table 3]
