Accelerated Hyper-Versus Normofractionated Radiochemotherapy with Temozolomide in Patients with Glioblastoma: A Multicenter Retrospective Analysis.


 Background and purpose

The standard treatment of glioblastoma (GB) patients consists of surgery followed by normofractionated radiotherapy (NFRT) with concomitant and adjuvant temozolomide chemotherapy. Whether accelerated hyperfractionated radiotherapy (HFRT) yields comparable results to NFRT in combination with temozolomide has only sparsely been investigated. The objective of this study was to compare NFRT with HFRT in a multicenter analysis.
Materials and methods

A total of 484 GB patients from four centers were retrospectively pooled and analyzed. 310 and 174 patients had been treated with NFRT (30×1.8Gy or 30×2Gy) and HFRT (37×1.6Gy or 30×1.8Gy twice/day), respectively. The primary outcome of interest was overall survival (OS) which was correlated with patient-, tumor- and treatment-related variables via univariable and multivariable Cox frailty models. For multivariable modeling, missing covariates were imputed using multiple imputation by chained equations, and a sensitivity analysis was performed on the complete-cases-only dataset.
Results

After a median follow-up of 15.7 months (range 0.8-88.6 months), median OS was 16.9 months (15.0-18.7 months) in the NFRT group and 14.9 months (13.2-17.3 months) in the HFRT group (p=0.26). In multivariable frailty regression, better performance status, complete versus not complete resection, MGMT hypermethylation, IDH mutation, smaller planning target volume, no steroid administration, and salvage therapy were significantly associated with longer OS (all p<0.01). Treatment differences (HFRT versus NFRT) had no significant effect on OS in either univariable or multivariable analysis.
Conclusions

This analysis suggests that HFRT and temozolomide is a safe option for patients wishing to shorten their treatment time and does not affect OS.


Introduction
Glioblastoma (GB) is the most common malignant tumor of the central nervous system in adults [1], accounting for approximately 69% of all malignant adult brain tumors in Germany [2]. The prognosis remains dismal, with a median overall survival (OS) for the whole patient population of not more than 21 months and 5-year OS rates below 10% after the standard-of-care trimodal therapy consisting of surgery followed by radiotherapy with concomitant and adjuvant temozolomide (TMZ) chemotherapy [3][4][5].
Several prospective glioma trials from the pre-TMZ era have evaluated the potential of HFRT [31][32][33][34]. HFRT with or without dose escalation however failed to show any superiority in OS or progression free survival (PFS) in comparison to normofractionated radiotherapy (NFRT) [35]. A more recently published randomized controlled trial comparing chemo-radiation protocols of dose escalated HFRT with NFRT did not point towards any benefit of HFRT in terms of OS [36].
In addition to the evidence for non-inferiority of HFRT in the abovementioned trials, HFRT schemes significantly reduce treatment time, even though the absolute number of radiotherapy treatments is higher in NFRT schemes. Therefore, the rationale to use it as an alternative scheme for patients is to shorten the time of irradiation from six to less than four weeks. The current treatment standard for primary GB is the combination of radiotherapy and TMZ according to the Stupp protocol [37] plus an eventual adjuvant therapy with tumor treating fields which may add a few months to the PFS and OS [4]. To the best of our knowledge, there are no prospective randomized trials investigating the efficacy of HFRT versus NFRT in combination with TMZ, with only a couple of monocentric retrospective studies reporting comparable outcomes [38,39]. To further explore any differences that might exist between HFRT and NFRT, we performed a retrospective multicenter analysis of GB patients which allowed us to model the effects of both schemes on PFS and OS while accounting for known prognostic factors.

Materials and Methods
This work is based on data of patients with newly diagnosed GB treated from 10/2004 to 7/2018 at four tertiary care institutions. Variables of interest were retrospectively collected by each center and inserted into a Microsoft® Access database, which contained pre-specified selection possibilities for the value assignment of categorical variables. The anonymized center-specific databases were then exported and pooled into a single Microsoft® Excel file, which was further processed with R statistical software.
Inclusion criteria for this analysis were at least one follow-up of OS, age ≥18 years and having received either HFRT (37×1.6Gy or 30×1.8Gy twice/day) or NFRT (30×1.8Gy, 30×2Gy or 30×1.8Gy with 2Gy simultaneously integrated boost). This resulted in 484 patients being eligible for analysis of which 310 had received NFRT and 174 HFRT. Differences between the NFRT and HFRT group were assessed using the Wilcoxon rank sum test and Fisher's exact test for continuous and categorical variables, respectively. PFS was calculated as the time difference between the start of treatment and first clinical or radiological progression. The latter had to be confirmed according to the Response Assessment in Neuro-Oncology (RANO) criteria [40] by a dedicated neuroradiologist. Follow-up included clinical and radiological evaluation every 3 months or depending on the patients' performance status.
The primary outcome of interest for this study was OS, with PFS as a secondary outcome, both calculated according to the Kaplan-Meier method. PFS and OS were correlated with patient-, tumor-and treatment-related variables via univariable and multivariable Cox frailty models. The frailty model accounts for clustering of patients within centers and therefore for any unmeasured center-specific factors that may affect the outcome (e.g. patient selection bias) [41,42]. The Cox frailty model can be written as Where is a vector of covariates for patient , ℎ 0 ( ) the baseline hazard function and is the random effect associated with the -th cluster (clinic) that acts as a multiplier on the baseline hazard. A log-normal distribution was assumed for the frailty terms, equivalent to a normal distribution of the random effects exp( ).
For multivariable modelling, we tried to utilize as many variables as possible [43,44] by imputing the covariates with missing information using the multiple imputation by chained equations R package 'mice' [45]. A "missing at random" mechanism was assumed being responsible for missing variables, with all variables given in Table 1 as well as follow-up time and OS being added into the imputation model. Variables were imputed in the order of their number of missing cases. Logistic regression and predictive mean matching were used for imputing binary and continuous variables, respectively. A total of 100 imputation data sets were created. Each was used to fit a Cox frailty regression model, and the final Cox model was obtained by pooling the coefficients of these 100 Cox models together. For sensitivity analysis, a Cox model using only the complete cases (with no missing values) was built and compared to the Cox model resulting from the imputed datasets.
Statistical significance was defined as p-values <0.01. This threshold was chosen based on the conversion between p-values and minimum Bayes factors [46]. Bayes factors (or likelihood ratios in case of simple hypotheses) measure the strength of evidence between two competing hypotheses [47]. In exploratory analyses, a p-value of 0.01 corresponds to a minimum Bayes factor of 1/6.5, providing moderate-strong evidence against the null hypothesis [46]. All analyses were performed within the statistical programming language R version 4.0.3 [48]. Table 1 shows the baseline characteristics of the 484 patients fulfilling the inclusion criteria for this study. There was a significant difference between the NFRT and HFRT group with respect to the treating institution, because HFRT patients mainly stemmed from two hospitals. Other significant differences between the NFRT and HFRT groups existed for planning target volume (PTV) (larger in the HFRT group), steroid administration (more frequent in the HFRT group), and salvage treatment (less frequent in the HFRT group).

Results
Kaplan-Meier survival curves stratified according to treatment protocol are shown in Figure 1. Median follow-up for the whole cohort was 15.7 months (range 0.8-88.6 months). During individual follow-up, tumor progression occurred in 385 out of 410 patients (it was unknown for 74 patients), and 394 of the 484 patients had died. Median PFS was 10.0 months (95% confidence interval [CI] 9.2-11.0) in the NFRT group and 7.9 months (6.6-9.3 months) in the HFRT group, which was almost significantly different (p=0.012). Median OS was 16.9 months (15.0-18.7 months) in the NFRT group and 14.9 months (13.2-17.3 months) in the HFRT group (p=0.26). OS was 16.0 months in both groups combined (95% CI 15.0-17.6 months).
A difference in PFS and OS between centers was evident from univariable Cox regression analysis ( Table 2). Given that these differences might reflect differences in some unmeasured hospital-specific variables, frailty models were used to fit the other univariable models reported in Table 2. Factors significantly associated with longer PFS were higher KPS, complete resection, MGMT hypermethylation, IDH mutation and simultaneous as well as sequential TMZ versus no TMZ administration. There was also a significant association between institutions and PFS, with patients from center 3 having shorter and patients from center 4 having longer PFS. Factors significantly associated with longer OS were younger age, higher KPS, complete resection, MGMT hypermethylation, IDH mutation, simultaneous as well as sequential TMZ versus no TMZ administration, no steroid administration and salvage treatment.
After creating 100 imputation datasets, we fitted a separate multivariable Cox frailty model to each of these datasets and pooled the regression coefficients and associated p-values. The average hazard ratios and p-values are given in Table 3. After controlling for many possible confounders, higher KPS, complete versus not complete resection, MGMT hypermethylation, IDH mutation, smaller PTV size, no steroid administration, and salvage therapy were significantly associated with longer OS (all p<0.01). RT schedule (HFRT versus NFRT) had no significant effect on OS (p=0.323). The frailty term modelling heterogeneity between hospitals was also statistically significant, showing the importance of its inclusion into the model. A multivariable Cox model fit to the original dataset with complete observations only resulted in similar hazard ratios (Table 4). However, the uncertainties were larger due to the smaller number of observed events (258 patients were not considered due to missing variables), and only MGMT hypermethylation, PTV size, and salvage therapy reached the threshold of statistical significance (p<0.01). The frailty term modelling heterogeneity between hospitals was confirmed as almost statistically significant (p=0.013).

Discussion
The main aim of this study was to compare patients who had received HFRT with those having received NFRT with respect to PFS and OS. The strength of this analysis is that patients were pooled from four different hospitals, resulting in a large cohort of 484 patients for which several putatively prognostic factors were known. The analysis indicates that HFRT yielded comparable outcomes to NFRT.
This confirms previous findings that accelerated HFRT achieves treatment effects comparable to NFRT in a shorter time frame [38,39]. The lower α/β ratio of normal brain tissue for late reactions, for which an α/β value of 2Gy is widely accepted, favors a hyperfractionated acceleration instead of a hypofractionated one. The biologically effective dose (BED) inside the tumor of HFRT protocols may vary considerably because of uncertainties about the kick-off time when repopulation sets in, , and the doubling time of repopulation. Radiobiological estimates for GB resulted in = 37 days and = 15.4 days, indicating a negligible influence of fractionation, but these estimates were obtained based on PFS rather than on the more reliable endpoint OS [49]. Using the formula provided by Lee  Thus, the BEDs of the three different schedules are not too different andgiven an additional clonal heterogeneity and clonal selection by radiotherapymay explain why the particular fractionation scheme had no significant effect on PFS or OS.
To the best of our knowledge, there was no favorable patient selection for HFRT protocols. Quite opposite, in one center (center 3) the HFRT scheme was preferably used in patients with involvement of or a tumor position close to the brain stem; this was also the case in some patients from center 1. Furthermore, some of the prognostic factors indicated a less favorable patient selection into the HFRT group: A significantly higher percentage of patients did receive concomitant steroids (76.1% versus 55.3%, Table 1), which is associated with poorer OS [51]. PTV sizes were also larger in the HFRT group, and a significantly lower percentage of patients received salvage therapy after tumor progression. For these reasons, we performed multivariable regression analysis with frailty terms in order to reduce any possible selection bias that is a general problem of retrospective analyses.
We were able to confirm the effect of several prognostic factors in univariable and multivariable analysis. The significant predictors were gross total versus non-complete resection, MGMT hypermethylation, IDH mutation status, and salvage treatment. Mutations in the IDH1 or IDH2 gene may predict a better prognosis by downregulating several hypoxia-inducible factor 1α (HIF-1α) target genes which are involved in glycolytic energy metabolism [52], among them lactate dehydrogenase (LDH) which catalyzes the conversion of pyruvate to lactate and thereby contributes to the acidification of the tumor microenvironment [53]. Accordingly, the LDH serum concentration was shown to be negatively correlated to OS in GB patients [33]. Both univariable and multivariable analyses also showed that OS was comparable between the HFRT and NFRT groups. However, there was a trend towards different PFS rates between the two groups, eventually pointing towards a systematic difference between centers in the evaluation of PFS within the scope of the RANO guidelines. The differentiation between pseudo-progression and real tumor regrowth is a difficult task in practice and our data reflect the real-world differences that exist between different applied criteria. In this respect, the difference in OS between institutions, which is an unambiguous endpoint, is harder to explain, but could be related to hospitalspecific factors such as patient selection. For example, patients from center 3 had on average significantly larger PTVs than patients from the other institutions (all p<0.0001). The negative selection was obvious in this patient cohort also by admission of HFRT in case of proximity to the brainstem. In addition, there might have been other unmeasured differences, so that the use of frailty models was a reasonable choice. This was confirmed by the significance of the frailty terms in multivariable modelling (Tables 3 and 4).
Another modality to reduce treatment time besides HFRT is hypofractionated accelerated radiation therapy. The use of hypofractionated protocols has been established in several randomized controlled trials [54][55][56][57]. This method is however primarily recommended for elderly patients and its role in the treatment of GB in younger individuals and those with good prognosis requires further research [58]. The other potentially interesting point of hypofractionation is the theoretical possibility of overcoming immanent or acquired radioresistance of tumor cells [44]. An alternative way of overcoming radioresistance could be a local dose escalation in biologically active tumor volumes, albeit the clinical data of dose escalation are needed to be proven in phase 3 trials [59]. Another possibility could be an intensification of chemotherapy regimens in selected patients which may also improve the treatment results with acceptable toxicity [60,61].
Despite the large sample size, this study had several limitations. First, the presented data are retrospective in nature, which implies potential uncertainties in some variables, missing data, as well as the possibility of systematic confounding and biases. We took account of the problem of missing variables by using multiple imputation by chained equations, which is preferred over just using the data with complete cases because it makes better use of the full information within the data [44]. Sensitivity analysis with only the complete cases yielded quantitatively similar results as the imputed dataset analysis and therefore confirmed the latter, although uncertainties in regression coefficients were larger, as would be expected due to the smaller sample size. For example, steroid administration was a statistically significant predictor of OS in multivariable regression on the imputed datasets but failed to reach the significance threshold in univariable or multivariable regression on the original dataset. Nevertheless, the obtained hazard ratio >1 is consistent with the detrimental effects of steroids found by Pitter et al. [51] which Klement and Champ [16] have attributed in part to the well-known effects of corticosteroids to raise blood glucose levels. In general, it is reassuring for the quality of our data that the hazard ratios of all variables are consistent with their expected effects on PFS and OS which are known from previous studies and tumor biology.
Another limitation is that the tolerability of the treatments was not directly evaluated because of the retrospective nature of this work. The need for steroids was assessed as an indicator of both tumor progression and treatmentrelated brain edema. Steroid administration was more frequent in the HFRT group. However, since prophylactic steroid administration was also performed, it cannot be considered an accurate surrogate parameter of toxicity. Previous data from Kaul et al. [38] and Lewitzki et al. [39] reported good tolerance of simultaneous temozolomide and HFRT.
In conclusion, our multicenter analysis suggests that radiochemotherapy with HFRT and TMZ is a safe option for patients wishing to shorten their treatment time and does not affect OS, even in the context of larger tumors or poorer performance status. MGMT hypermethylation, smaller tumor size, and salvage therapy were confirmed to have the most significant positive impact on OS. To the best of our knowledge, both HFRT protocols were well tolerated without excessive acute or late toxicity, albeit no dedicated analysis was performed due to the retrospective nature. Analyses of neurocognition and quality of life in particular should also be performed in the future for a better assessment of tolerance. The rate of treatment interruptions was low in all protocol groups.

Conflicts of Interest Statement
The authors declare that this analysis has been conducted without any financial or other relevant conflicts of interest.

Data Availability Statement
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Animal Research Statement
There was no animal research performed.

Consent to Participate Statement
Ethical approval was waived by the local Ethics Committees of University Clinic of Würzburg, University Clinic Charité Berlin and University Clinic of Freiburg in view of the retrospective nature of the study and all the procedures being performed were part of the routine care.

Consent to Publish Statement
Statistical analysis, procession and publishing of complete anonymized patient data were approved by the local Ethics Committees of University Clinic of Würzburg, University Clinic Charité Berlin and University Clinic of Freiburg.

Plant Reproducibility Statement
Our research group is planning a prospective observational trial to prove the data.

Clinical Trials Registration
This research study was conducted retrospectively from data obtained for clinical purposes. No trial protocol and registration are applicable.

Author contributions
All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Victor Lewitzki, Rainer Klement, Ilinca Popp, David Caul, Felix Ehret and Reinhardt Sweeny. The first draft of the manuscript was written by Rainer Klement and Victor Lewitzki and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Funding
The authors did not receive support from any organization for the submitted work.  HFRT: Hyperfractionated radiotherapy NFRT: Normofractionated radiotherapy.