Subventricular zones: new key targets for glioblastoma treatment

Background We aimed to identify subventricular zone (SVZ)-related prognostic factors of survival and patterns of recurrence among patients with glioblastoma. Methods Forty-three patients with primary diagnosed glioblastoma treated in our Cancer Center between 2006 and 2010 were identified. All patients received surgical resection, followed by temozolomide-based chemoradiation. Ipsilateral (iSVZ), contralateral (cSVZ) and bilateral (bSVZ) SVZs were retrospectively segmented and radiation dose-volume histograms were generated. Multivariate analysis using the Cox proportional hazards model was assessed to examine the relationship between prognostic factors and time to progression (TTP) or overall survival (OS). Results Median age was 59 years (range: 25–85). Median follow-up, OS and TTP were 22.7 months (range 7.5–69.7 months), 22.7 months (95% CI 14.5–26.2 months) and 6.4 months (95% CI 4.4–9.3 months), respectively. On univariate analysis, initial contact to SVZ was a poor prognostic factor for OS (18.7 vs 41.7 months, p = 0.014) and TTP (4.6 vs 12.9 months, p = 0.002). Patients whose bSVZ volume receiving at least 20 Gy (V20Gy) was greater than 84% had a significantly improved TTP (17.7 months vs 5.2 months, p = 0.017). This radiation dose coverage was compatible with an hippocampal sparing. On multivariate analysis, initial contact to SVZ and V20 Gy to bSVZ lesser than 84% remained poor prognostic factors for TTP (HR = 3.07, p = 0.012 and HR = 2.67, p = 0.047, respectively). Conclusion Our results suggest that contact to SVZ, as well as insufficient bSVZ radiation dose coverage (V20Gy <84%), might be independent poor prognostic factors for TTP. Therefore, targeting SVZ could be of crucial interest for optimizing glioblastoma treatment.


Background
Glioblastoma is the most common primary brain tumor among adults, with poor outcome following chemoradiotherapy, despite deeper insights into molecular biology and significant advances in therapeutics [1]. Like other cancers, glioblastomas are characterized by a high intratumoral heterogeneity in a wide range of phenotypic and functional features [2][3][4]. The stochastic model has traditionally been accepted as the basis of tumor heterogeneity, wherein a small population of genetically unstable clonal cells randomly acquires the appropriate somatic mutations necessary to confer extensive proliferative capabilities [5]. However, recent evidence suggests that according to the alternative cancer stem cell model [6], glioblastomas are rather organized hierarchically as they contain a subpopulation of cancer cells with stem cell characteristics, including self-renewal capacity and multilineage potency, as well as tumor-initiation/proliferation and migration ability [7][8][9][10][11].
The origin of brain tumor stem cells is still controversial, but neural stem cells are highly likely to be candidate as both populations share many similarities, including molecular pathways (Notch receptor activation [12,13], inactivation of PTEN tumor suppression gene [14]) and markers of gene expression (like Nestin, CD133 or Sox [15]). A growing body of evidence has also suggested that any glioma cell could acquire a stem cell phenotype in response to hypoxia, or even radiation, contributing to GBM radioresistance [16][17][18]. Stem cells of various tissues exist within protective niches composed of a number of differentiated cell types providing direct cell contacts and secreting factors that maintain stem cells primarily in a quiescent state [19]. It is then hypothesized that tumor stem cells might arise from normal stem cells that have acquired mutations which enable them to escape from niche control or after deregulation of extrinsic factors within the niche, leading to uncontrolled proliferation of stem cells and tumorigenesis [20][21][22]. The subventricular zone (SVZ) is the largest niche of neurogenesis in the adult mammalian brain [23,24], so that an emerging hypothesis is that brain tumor stem cells could derive from the SVZ [11]. It should be noted that the subgranular zone abutting the hippocampal dentate gyrus is a secondary nich in adult mammalian brains, but where neurogenesis only occurs in foci closely associated with blood vessels [25].
Many data have already provided interesting relationships between SVZ involvement in glioblastoma and outcome or pattern of relapse, strongly suggesting that tumors contacting SVZ are associated with multifocal presentation at diagnosis, multifocal or distant progression [26][27][28][29] and above all decreased survival [27,[29][30][31]. This stresses the need to better understand the role of the SVZ as a potential source of glioblastoma through initiation and promotion, a potential source of relapses through repopulation of tumor, and as a target for glioblastoma treatment.
Glioblastoma stem cells are known to be inherently resistant to chemotherapy and radiotherapy [32,33]. One strategy to overcome glioblastoma stem cells radioresistance could consist in increasing the radiation dose to the SVZ. The impact of incidental radiation dose to the SVZ during the course of a standard chemoradiotherapy on the outcome of patient has already been assessed, but mostly in inhomogeneous series and with somewhat contradictory results [34][35][36][37][38][39][40][41]. These results are all the more confusing that such irradiation to neural stem cells is source of radiation-induced neurotoxicity [42], and is in contradiction with hippocampal sparing strategies [43].
Therefore, the objective of this retrospective study was to identify SVZ-related prognostic factors of survival and patterns of recurrence, with regard to both surgical and radiotherapeutic management, to better integrate SVZ in the global therapeutic strategy of glioblastoma patients, with as few neurotoxicity as possible.

Patient selection
Patients with primary diagnosed glioblastoma consecutively treated in our Comprehensive Cancer Center between April 2006 and March 2010 were identified through the electronically records database. Inclusion criteria were: a histopathologically proved glioblastoma, treated by surgical total or subtotal resection (thus excluding patients with a "biopsy-only") and full course of radiation, with immediate (<48 h) post-operative magnetic resonance image (MRI) to assess tumor resection and with documented radiological follow-up for at least 6 months after surgery. Patient data, including demographics, imaging data, treatment and clinical outcomes were retrospectively collected. The cut-off date for the analysis was September 2014. This study was approved by our Institutional Review Board and our local ethics committee. Each participant/participant's guardian had provided consent for this study.

Imaging data
Pre-operative imaging data (T1-weighted MRI pre-and post-contrast and FLAIR sequence) were assessed to define the tumor volume and the tumor localization with regard to the SVZ. Given the difficulty to differentiate between non-tumoral vasogenic-edema versus tumorinfiltrative area within the non-enhancing FLAIR hyperintensity lesion in T2-weighted sequences, contacting tumors were defined as tumors having a distance of 0 mm between the contrast enhancing tumor edge and the SVZ, as performed in most previous studies [26][27][28][29][30].
Post-operative post-contrast T1-weighted imaging was used to distinguish : patients with gross total resection (GTR) in case of no residual enhancement or near total resection (NTR) in case of only rim enhancement of the resection cavity, from patients with subtotal resection (STR), namely with residual nodular enhancement [44]. Despite a new RANO classification for evaluation of completeness of resection has been proposed recently which accounts for both enhancing and non-enhancing components of tumor [45], the post-contrast T1-wheighted imaging based classification has been used because it is still widely utilized in studies [26,30,46].

Treatment
All patients received a non-biopsy surgical resection, followed by a standard of care adjuvant temozolomidebased chemoradiation, as previously described [1]. Patients were CT-scanned in treatment position with a slice distance of 2.5 mm. Radiotherapy was delivered by a three dimensional conformal radiotherapy, with a standard dose of 60 Gy in 30 daily fractions, 5 days a week. The gross tumor volume (GTV) was the contrastenhancing lesion on T1-weighted MRI or the surgical cavity with residual contrast-enhancement. The clinical target volume (CTV) was built with a 20-mm margin around the GTV, and then was edited to include the FLAIR signal abnormality and adjust it to anatomic barriers. The planning target volume (PTV) was finally built with a 3-mm isotropic margin expansion.
Multimodal MRI was performed 1 month following the completion of chemoradiotherapy, and every 3 months thereafter until relapse.

SVZ delineation and dosimetry data collection
According to Barani's delineation suggestions [47], ipsilateral (iSVZ), contralateral (cSVZ), and bilateral (bSVZ) SVZs were retrospectively segmented by a single physician as a 5-mm lateral margin of the lateral ventricles based on the patient's original treatment planning CT scan. In a neural stem cells sparing attempt, two delineation methods were assessed: with (TH+) [34,36,38,39] and without (TH-) [37] temporal horn (Fig. 1). Dosevolume histograms were retrospectively generated on the original plans, and dose-volume parameters to each SVZ volume were extracted: mean dose and volume of SVZ (%) receiving more than x Gy (VxGy).

Patterns of recurrence
Progression was defined according to the Response Assessment in Neuro-Oncology Working Group (RANO) criteria [48] and all recurrences were centrally reviewed by an experienced neuroradiologist. For pattern of recurrence analysis, recurrence characterization was defined as: local (recurrence epicenter within 2 cm from the initial edge of T1-gadolinium contrast enhancement lesion) or distant recurrence; within or outside the SVZ.

Statistical analysis
The data were summarized through frequency and percentage for categorical variables and by median and range for continuous variables. Comparisons between groups were performed using the chi-square test or the Fisher's exact test for qualitative data, and the Wilcoxon test for paired quantitative data.
All survival times were calculated from the date of surgery and were estimated by the Kaplan Meier methods with 95% confidence intervals (CI), using the following firstevent definitions: progression according to RANO for time to progression (TTP) and death from any cause for overall survival (OS). Patients that were still progression-free for TTP or alive for OS were censored at the time of their last follow-up. Univariate analysis was used to identify prognostic factors for TTP and OS: the log-rank test was used for categorical variables and the Cox proportional hazards model for continuous variables. For SVZ dosimetric data analysis, several dose-volume parameters were used as variables. Patients were also dichotomized using published mean dose to SVZ volumes cut-off of 40 Gy [38] and 43 Gy [34], as well as using optimal cut-off searched with the minimal p-value approach (1000 bootstrap internal validation). Multivariate analysis was assessed using the Cox proportional hazards model including significant covariates in univariate analysis.
All reported p-values were two-sided. For all statistical tests, differences were considered significant at the 5% level.
All statistical analyses were conducted using STATA 13.0 software.

Patient characteristics
Forty-three patients were consecutively included. Median age at surgery was 59 years (range: 25-85). Twenty six patients (63.4%) presented a tumor contacting the SVZ.
All patients completed radiation course, with interruption of temozolomide for only 5 patients (11.36%), due to toxicity. Patient characteristics are listed in Table 1.  Median iSVZ, cSVZ and bSVZ volumes/doses are summarized in Table 2.
On univariate analysis, age was not a prognostic factor for OS or TTP, as well as sex, ECOG performance status or O6-methylguanine-DNA-methyltransferase (MGMT) methylation status. By contrast, initial contact to SVZ was a poor prognostic factor for OS (18.7 vs 41.7 months, p = 0.014) and TTP (4.6 vs 12.9 months, p = 0.002) ( Table 3).

Surgical parameters
The extent of resection (GTR/NTR vs. STR) was not found to be a prognostic factor for either OS or TTP (Table 3).

SVZ dosimetric parameters
Using TH+ method, no dose-volume parameter as a continuous variable was found to be a prognostic factor for either OS or TTP. With TH-method, no dosevolume parameter was found to be a predictive factor for OS, whereas V20Gy to cSVZ (HR = 0.99 (95% CI 0.98-1), p = 0.083) and V20Gy to bSVZ (HR = 0.99 (95% CI 0.98-1), p = 0.099) trended as prognostic factors for TTP.  Additionally, the patient cohort was divided into high and low radiation dose groups based on two published cut-off values for mean dose to SVZ volumes (40 Gy and 43 Gy) ( Table 4). With TH+ method, patients receiving mean dose to iSVZ >43 Gy had a poorer TTP (4.8 vs 12.1 months, p = 0.072) and a poorer OS (20.5 vs 26.2 months, p = 0.087), but not significantly.
With TH-method, patients receiving mean dose to bSVZ >40 Gy had a significantly improved TTP of 9.4 months compared to 4.6 months in patients who received less than 40 Gy (p = 0.023).
We then studied the pattern of recurrence of tumors according to their initial contact to SVZ, and found that 72% (18/25) of patients with initial SVZ contacting tumor recurred within SVZ, while only 57% (8/14) of patients with non-contacting SVZ tumors recurred within SVZ.
Finally, looking at the pattern of recurrence according to the radiation dose delivered to the SVZ, we found that patients with V20Gy to bSVZ ≤84% preferentially recurred within the SVZ (72.4% in case of exclusive local recurrence (21/29), 71.9% in case of local recurrence with or without distant recurrence (23/32), and 71.9% in case of local and/or distant recurrence (23/32)), while on

Discussion
Owing to complex relationships between brain tumor stem cells and neural stem cells, the optimal management of neurogenic regions of patients with glioblastoma still remains controversial. Herein, we focused on the survival rates and the recurrence patterns of glioblastoma with respect to SVZ-related factors.
In spite of inherent biases due to its retrospective design, this study has several strengths. Indeed, it included a uniform cohort of consecutively-treated glioblastoma patients with the same first-line therapy (surgical resection and adjuvant temozolomide-based chemoradiation with uniform dose/fractionation scheme), and uniform follow-up using most updated criteria, of more than 6 months after surgery, to avoid any pseudo-progression.
Volume delineation, plan evaluation and approval were ensured by a single experienced radiation oncologist (ECJM). MGMT methylation status was available for almost the whole population, as well as all classical glioblastoma prognostic factors. Our salvage therapies were not homogenous (of the 38 patients who benefited a salvage treatment, 17 received a re-resection, most of them received in situ carmustine wafers (CW), and the 21 remaining patients received bevacizumab as second line, temozolomide or fotemustine); however this cofounding factor could influence OS but not TTP. Finally, our multivariate analysis aimed to control any other confounding factors. It should there be mentioned that as compared to intensity modulated radiotherapy which uses inverse planning with predefined constraints to structures, 3DCRT uses a number of fields and a beam arrangement defined by the operator, which can introduce some confounding factors in the assessment of dose to SVZ too difficult to take into account and therefore not included in our multivariable model.
Our results suggest first that initial contact to SVZ is a strong prognostic factor for TTP in multivariate analysis. This result is in accordance with several other similar reports [27,[29][30][31], and could be explained because glioblastoma arising in SVZ might have a higher percentage of more potent cells, making them more invasive and infiltrative. This hypothesis is reinforced by an established correlation between SVZ involvement and invasive tumor phenotype with distant or multifocal progression [26][27][28][29].
From a biological point of view, such an invasive phenotype is supported by experiments showing that SVZ neuronal progenitor cells have a higher migratory potential compared to non-SVZ ones [49].
Contrary to several studies [30,44,50], we failed to correlate the surgical status to survival. This could be explained either by an underpowered sample size, or because initial contact to SVZ, even in case of GTR or NTR, would lead to a microscopic residual disease due to the aggressive phenotype of tumor cells from SVZ.
Interestingly, 83.3% of our patients with distant relapse presented a recurrence within SVZ, and all but one had a GTR, highlighting the role of neurogenic regions on initiation, promotion but also repopulation of tumors. Brain stem tumor cells in neurogenic regions seem to constitute a source for a subsequent relapse.
Overall, survival data and recurrence patterns from both literature and our study suggest that a "prophylactic" irradiation to SVZ could be necessary to perform a therapeutic intensification among contacting tumor patients to encompass therapeutic resistance, but also to eradicate potential "niches" of relapse among other glioblastoma patients. Several studies have established a correlation between a dose level to SVZ and survival outcomes [34][35][36][37][38][39][40]. The first ever study from Evers et al.
(UCLA, USA) found that mean dose >43 Gy to bSVZ was an independent factor for PFS (HR = 0.73, 95%CI = [0.57; 0.95]), but not for OS on multivariate analysis [34]. These results were not confirmed in the series from VUMC (The Netherlands) which could not find an association between this dose to SVZ and PFS [35]. This may be due to the more aggressive population of exclusively grade 4 patients in the VUMC study (compared to a mix grade 3/grade 4 glioma population in the UCLA study), which might require a higher dose to bSVZ. Given these heterogeneous results, Lee et al. updated and pooled these data restricted to only glioblastoma patients, but they failed to find a correlation between bSVZ dose and PFS as in Evers's study. However, they found that increasing mean dose to iSVZ (>59.  [38], and more recently by Adeberg et al. with the same cutoff but only in univariate analysis (median PFS : 8.5 versus 5.2 months; p = 0.013) [40]. Interestingly, contradictory results came from a recent study by Elicin et al. which found a worse PFS among patients with high iSVZ dose (>62.25 Gy) in both subgroups of good performance status and SVZ without tumoral contact. A worse PFS was also found among patients with dose to cSVZ >59.2 Gy (7.1 vs 10.4 months), but it was not confirmed in multivariate analysis [39]. Gupta et al. also found a worse PFS and OS with higher cSVZ doses, but patients in higher dose group had actually larger tumors, more often crossing the midline and with less GTR so that it was not confirmed in multivariate analysis [37]. On the contrary, Adeberg et al. found a better PFS among patients receiving a mean dose to cSVZ ≥30 Gy (HR in multivariate analysis of 0.45, 95%CI = [0.20-0.98]) [40]. A summary table of all previous results and our main results is presented (Table 6). However, several limits should be noted in these retrospective studies, which could explain some of these discordances. Population was not always homogenous : grade 3 gliomas were included in Evers' study [34] and patients could receive other drugs than temozolomide in two studies [34,36]. Methylation status, which is a strong predictive factor of outcome [51], was missing or insufficiently reported except in two studies [37,40].
Variability in the radiotherapy target volume delineation and dose-prescription between studies (one vs. two dose-level prescription depending on European or American guidelines) but also within a series [34,36,40], and variability of the SVZ delineation between studies (with [37] vs. without temporal horn [34,36,38,39]) and within a series (3-5 mm thick), should also be mentioned. Above all, only mean dose to SVZ was considered in all studies, and thresholds were defined using only median mean doses or 75 th percentile.
Instead, our approach consisted of exploring dose distributions to SVZ during radiotherapy through a more complete analysis of the SVZ dose-volume histogram, in order to better understand the role of low/mid doses to SVZ in addition to high dose. This innovative approach had never been assessed before. We also tried to find an optimal cut-off for relevant dose-volume parameters, through a statistical method. It is there worth mentioning that due to the exploratory nature of the analysis, no adjustment for multiple testing was used. This led us to identify a V20Gy >84% to bSVZ (without temporal horn) as an independent prognostic factor for better TTP. The multivariate analysis included contact to SVZ and radiation dose to SVZ as co-variables and could then take into account the fact that patients with contacting tumors inherently received a higher dose to SVZ. Interestingly, in univariate analysis, a mean dose to bSVZ >40Gy was also a factor for better TTP, in the same way as Evers [34]. Last, it should be noted that a poorer TTP and OS has been found among patients with mean dose to iSVZ (with temporal horn) >43 Gy; however this did not reach statistical significance, and additionnally this was probably due to a larger proportion of contacting tumors in the iSVZ high dose group (80%) compared to the low dose group (30%).
Another feature of our approach was to explore two SVZ delineation methods (based on a standard and reproducible 5-mm lateral margin of the lateral ventricle [47]) : with (TH+) or without (TH-) temporal horn. The aim was to assess an hippocampal sparing method (TH-), as most neurogenic regions with potential brain tumor stem cells stand lateral to the body of lateral ventricles and biologically equivalent doses in 2-Gy fractions to 40% of the bilateral hippocampi greater than 7.3 Gy was shown to be associated with long-term neuro-cognitive impairment [52]. Hippocampal avoidance during whole-brain radiotherapy for brain metastases was associated with preservation of memory and quality of life in a phase II trial [43]. Hippocampal avoidance is not assessed in glioblastoma trials because optimal PTV coverage is the main objective. But when it comes to performing "prophylactic" irradiation to SVZ, neurocognitive preservation should be envisaged, as we confirmed that SVZ irradiation is compatible with hippocampal sparing through TH-delineation method. Chen et al. argued that iSVZ dose ≥40 Gy with TH+ method did not correlate with worsened patient KPS score at the end of radiation treatment, however neurocognitive performances are not precisely assessed in KPS [38]. We thus believe that delineation for glioblastoma patients should then include bSVZ delineation without temporal horn, as well as hippocampi delineation, with use of intensity modulated radiotherapy planning to ensure a sufficient coverage to bSVZ while ensuring an hippocampal avoidance according to previous constraints.
Regarding patterns of recurrence, larger proportion of patients with initial SVZ contacting tumor recurred within SVZ compared to patients with non-contacting tumors. This seems to be due to the in-field pattern of recurrence typically seen in glioblastoma following radiotherapy as most relapses occur within 2 cm beyond the contrast enhancement, and almost exclusively in the edema area [53][54][55]. Additionally, patients with V20Gy >84% turned out to relapse away from SVZ contrary to those with V20Gy ≤84%, who mainly relapsed within SVZ. Radiobiological considerations could provide some interesting explanations to this pattern, as there is strong evidence that sublethal irradiation promotes glioblastoma cells migration, via an increase in the αvβ3 integrin expression [56]. Therefore, hypothesizing that SVZ is a potential source of recurrences through repopulation of tumor by remaining brain tumor stem cells following chemoradiation [11], the level of bSVZ irradiation could explane the location of recurrence within vs without SVZ.
Overall, these data could support two hypotheses. The simplest one would be that patients who received a V20Gy >84% to bSVZ really benefited from such a prophylactic irradiation through a brain tumor stem cells eradication. An additional and/or alternative one would be that some of patients who received a higher bSVZ irradiation experienced an "artificial" increase in TTP via a more durable transition through a migratory phenotype, before engaging towards a proliferative one. Then a strategy to encompass such a mechanism of radioresistance would be to associate prophylactic irradiation of bSVZ with intensity modulated radiotherapy (and hippocampal sparing) ensuring a V20Gy >84%, and inhibitors of migration which could target integrinmediated signalling pathway. This strategy has never been assessed but could provide promising results. Meanwhile, some more robust data should soon come from three ongoing or recently concluded clinical trials regarding the impact of SVZ irradiation in GBM (NCT01478854 exploring a deliberate sparing of neurogenesis niches, contrary to NCT02039778 and NCT02177578, exploring a deliberate SVZ irradiation).

Conclusion
To conclude, our data from a homogeneous patient cohort suggest that contact to SVZ could be an independent poor prognostic factor for TTP, as well as bSVZ insufficient dose coverage such as a V20Gy ≤84%. We also provide data suggesting that the pattern of recurrence is influenced by contact to SVZ and by radiation dose to bSVZ. These results highlight the prominent role of the SVZ in the prognostic of glioblastoma patients through mechanisms of therapeutic resistance, and therefore raise the question of their role in the optimal management of glioblastoma. Further prospective clinical trials assessing prophylactic irradiation of bSVZ by intensity modulated irradiation, with or without migration-inhibitors, would be critical in determining whether such strategy could improve devastating outcomes of glioblastoma patients.