Advertisement
Personalized Medicine and Imaging

Development and Validation of a Gene Signature for Patients with Head and Neck Carcinomas Treated by Postoperative Radio(chemo)therapy

Stefan Schmidt, Annett Linge, Alex Zwanenburg, Stefan Leger, Fabian Lohaus, Constanze Krenn, Steffen Appold, Volker Gudziol, Alexander Nowak, Cläre von Neubeck, Inge Tinhofer, Volker Budach, Ali Sak, Martin Stuschke, Panagiotis Balermpas, Claus Rödel, Hatice Bunea, Anca-Ligia Grosu, Amir Abdollahi, Jürgen Debus, Ute Ganswindt, Claus Belka, Steffi Pigorsch, Stephanie E. Combs, David Mönnich, Daniel Zips, Gustavo B. Baretton, Frank Buchholz, Michael Baumann, Mechthild Krause and Steffen Löck
DOI: 10.1158/1078-0432.CCR-17-2345 Published March 2018

Abstract

Purpose: The aim of this study was to identify and independently validate a novel gene signature predicting locoregional tumor control (LRC) for treatment individualization of patients with locally advanced HPV-negative head and neck squamous cell carcinomas (HNSCC) who are treated with postoperative radio(chemo)therapy (PORT-C).

Experimental Design: Gene expression analyses were performed using NanoString technology on a multicenter training cohort of 130 patients and an independent validation cohort of 121 patients. The analyzed gene set was composed of genes with a previously reported association with radio(chemo)sensitivity or resistance to radio(chemo)therapy. Gene selection and model building were performed comparing several machine-learning algorithms.

Results: We identified a 7-gene signature consisting of the three individual genes HILPDA, CD24, TCF3, and one metagene combining the highly correlated genes SERPINE1, INHBA, P4HA2, and ACTN1. The 7-gene signature was used, in combination with clinical parameters, to fit a multivariable Cox model to the training data (concordance index, ci = 0.82), which was successfully validated (ci = 0.71). The signature showed improved performance compared with clinical parameters alone (ci = 0.66) and with a previously published model including hypoxia-associated genes and cancer stem cell markers (ci = 0.65). It was used to stratify patients into groups with low and high risk of recurrence, leading to significant differences in LRC in training and validation (P < 0.001).

Conclusions: We have identified and validated the first hypothesis-based gene signature for HPV-negative HNSCC treated by PORT-C including genes related to several radiobiological aspects. A prospective validation is planned in an ongoing prospective clinical trial before potential application in clinical trials for patient stratification. Clin Cancer Res; 24(6); 1364–74. ©2018 AACR.

Translational Relevance

Patients with HPV-positive, locally advanced head and neck squamous cell carcinomas (HNSCC) show a very good locoregional tumor control (LRC) after postoperative radiochemotherapy (PORT-C) and are therefore candidates for trials on treatment deescalation to reduce toxicity. For patients with HPV-negative HNSCC, additional biomarkers are urgently needed to identify subgroups of patients, (i) who are unlikely to respond to standard PORT-C and may benefit from treatment escalation, or (ii) who will likely not develop locoregional recurrences. We developed and independently validated a 7-gene signature prognostic for LRC of HPV-negative tumors, which is based on several radiobiological parameters or mechanisms. The prognostic performance of this radiobiology-based signature combined with clinical parameters was higher than that of a model containing hypoxia-associated genes and CSC markers only. After additional prospective validation, the 7-gene signature may be applied in clinical trials for patient stratification.

Introduction

Head and neck cancer is the sixth most frequently occurring tumor entity worldwide (1) with an overall 5-year survival rate of about 50% (2). Patients with resectable, locally advanced head and neck squamous cell carcinomas (HNSCC) who are at high risk for tumor recurrence are being routinely treated with postoperative radiochemotherapy (PORT-C). According to the results of three randomized clinical trials, concurrent chemotherapy leads to improved locoregional tumor control (LRC) and prolonged overall survival (OS) compared with postoperative radiotherapy (PORT) alone (3–5). Within the last years, radiotherapy of locally advanced HNSCC has been further improved through the development of new treatment techniques, such as intensity-modulated radiotherapy. Despite the increase in treatment efficacy, patients show a very heterogeneous treatment response. Therefore, the consideration of the individual tumor biology by appropriate biomarkers in addition to well-established clinical parameters may further improve patient stratification for treatment escalation or deescalation strategies.

Besides the consumption of alcohol and tobacco as well-known risk factors for the development of HNSCC, infection with the human papillomavirus (HPV) has been identified as another independent parameter. Also, the incidence of HPV infection in HNSCC has been increasing within the last decade (6). Preclinical and clinical studies have shown that HPV-positive HNSCCs are more radiosensitive than HPV-negative tumors (7, 8). To investigate the impact of HPV in patients who receive PORT-C and to identify additional biomarkers for patient selection, a retrospective, multicenter study of the German Cancer Consortium Radiation Oncology Group (DKTK-ROG) was conducted (9–13). For this cohort, we have shown that patients with HPV16 DNA–positive tumors have superior LRC and OS compared with patients with HPV-negative tumors (9). In particular, 98% of the HPV-positive and only 80% of the HPV-negative oropharyngeal tumors were locoregionally controlled. For patients with HPV DNA–negative tumors, additional biomarkers are urgently needed to identify subgroups of patients, who are unlikely to respond to PORT-C and may benefit from treatment escalation, or who are not anticipated to develop locoregional recurrences.

Tumor hypoxia has been shown to be correlated with increased radioresistance (14). For patients with locally advanced HNSCC, pretreatment hypoxia was significantly associated with low tumor control and OS after primary radio(chemo)therapy compared with patients with highly oxygenated tumors (15, 16). Several hypoxia gene classifiers have been developed in the last decade to assess hypoxia or hypoxia-related changes on the transcriptional level using routinely taken pretreatment biopsies (17, 18). We have recently shown their prognostic validity for patients at high risk of loco-regional failure receiving PORT-C (12). The association of hypoxia and LRC after PORT-C is unexpected as the gross tumor has been removed and subsequently remaining tumor cells are very unlikely to differ in hypoxia (12). This suggests that hypoxia impacts LRC not only by a direct biochemical effect on cellular radioresistance but also by other radiobiological mechanisms (12, 19). Recent studies reported that hypoxia as an external factor also favors increased radioresistance of cancer stem cells (CSC) and invasive tumor growth (reviewed in refs. 20, 21), and CSCs are known to play a major role in radioresistance and tumor recurrence (reviewed in ref. 22). The putative CSC markers CD44, SLC3A2, and MET were shown to be prognostic for LRC in patients who received PORT-C (12). The combined application of hypoxia-associated gene panels and CSC markers further improved patient stratification regarding their risk of locoregional treatment failure (23, 24), which was independently validated (25).

A gene signature including additional radiobiological aspects may predict patient outcome with even higher accuracy. In the literature, the number of gene panels for stratification of patients with HNSCC is steadily growing. However, to the best of our knowledge, they have been developed for patients who received primary radiotherapy (26, 27) or have not been linked to a specific treatment (28–30). Gene signatures prognostic for the response of patients with locally advanced HNSCC to PORT-C covering a broad spectrum of radiobiological aspects are still missing.

Therefore, the major aim of this study was to develop and validate a gene signature and corresponding statistical model for patient stratification beyond HPV infection status to improve the risk assessment for patients with locally advanced HNSCC who receive PORT-C. For the development of this gene signature, a gene set was composed in-house using a hypothesis-driven approach. The gene set incorporated genes that cover many radiobiologically important aspects such as DNA repair, cell cycle, epithelial–mesenchymal transition, CSC markers, hypoxia, proliferation, invasion, metastasis, as well as genes that were reported to be involved in cisplatinum resistance. The signature was developed and independently validated on two large patient cohorts for the primary endpoint LRC and the secondary endpoints OS and freedom from distant metastases (DM). To find the optimal results, internal validation methods were applied, and several statistical methods were compared, including advanced machine learning techniques.

Materials and Methods

Patients

Two different cohorts of patients with locally advanced HNSCC were being considered for this study. The training cohort consisted of 221 patients who were treated with PORT-C between 2004 and 2012 within the 9 partner sites of the DKTK-ROG. Inclusion criteria, data collection, handling and analyses of biomaterial were previously described (9, 12). Briefly, all patients received curatively intended cisplatinum-based PORT-C according to standard protocols with a minimum follow-up of 24 months and presented with a tumor stage pT4 and/or >3 positive lymph nodes and/or positive microscopic resection margins and/or extracapsular spread. The validation cohort consisted of 152 patients who were enrolled by the following criteria: not included in the previous DKTK-ROG training cohort, histologically proven HNSCC, treatment between 1999 and 2006 with PORT or PORT-C according to standard radiotherapy protocols with curative intention (25).

Preparation of biomaterials and biomarker analyses

Formalin-fixed paraffin-embedded (FFPE) blocks of the primary tumor specimens (removed by surgery) were first subjected to hematoxylin and eosin staining to histologically confirm the presence of squamous cell carcinoma. Afterwards, they were processed under standardized procedures for biomarker investigations. DNA extraction and PCR-array based analyses of HPV status have been performed as described previously (9). Briefly, genomic DNA was extracted from 5-μm FFPE sections using the QIAamp DNA FFPE tissue kit (Qiagen). HPV DNA analyses including genotyping were performed using the LCD-Array HPV 3.5 kit (CHIPRON GmbH) according to the manufacturer's instructions.

For both cohorts, gene expression analyses were performed consecutively using nanoString elements technology (nanoString Technologies, Seattle, WA, USA) as described in (12, 25). Briefly, total RNA as well as reporter and capture probes specific to the genes of interest were mixed and incubated at 62°C for 22 hours. Samples were then kept at 4 °C for a maximum of 18 hours and subjected to the nCounter system. Raw counts were logarithmized and then normalized by subtracting the mean of the log-transformed counts of the reference genes ACTR3, B2M, GNB2L1, NDFIP1, POLR2A, RPL11, and RPL37A. Because of insufficient tumor material or too low RNA yield, some of the samples had to be omitted from the analysis. In the training and validation cohort, NanoString and HPV analyses could be performed for 195 and 142 samples, respectively. The expression levels of 178 genes were evaluated by NanoString analyses for both cohorts. The genes were selected by a literature search on a hypothesis-driven basis. Genes were included that have previously been reported to be associated with sensitivity or resistance to radio(chemo)therapy, that is, genes involved in proliferation, invasion, and metastasis; tumor hypoxia–associated genes, genes encoding for putative CSC markers, and DNA repair as well as genes that have been associated with cisplatinum resistance (see Supplementary Table S1).

Study design

The aim of this study was to develop and validate a gene signature for patient stratification beyond HPV infection status to further improve the risk assessment for patients with locally advanced HNSCC who receive PORT-C. Therefore, only patients with HPV16 DNA–negative tumors and available NanoString gene expressions were included (N = 130/221 training, N = 121/152 validation). The study design is presented in Fig. 1. Prognostic models including the following parameters should be compared: the identified gene signature alone, clinical parameters alone, and the combination of clinical parameters and the identified gene signature.

Figure 1.
Figure 1.

Study design.

Statistics and clinical endpoints

The primary endpoint was locoregional tumor control (LRC). Secondary endpoints were freedom from distant metastases (DM) and overall survival (OS). All endpoints were calculated from the first day of radiotherapy to the date of event or censoring. Death was considered a competing risk for locoregional recurrence and DM, while locoregional recurrence and DM did not cause censoring. Survival curves were estimated by the Kaplan–Meier method and compared by log-rank tests. Differences between the training and validation cohort were evaluated by Mann–Whitney U tests for continuous variables and by χ2 tests for categorical variables. Descriptive analyses and the described statistical tests were performed using SPSS 23 (IBM Corporation). A statistical framework was developed to identify gene signatures and corresponding prognostic models to optimally and robustly predict the primary and secondary endpoints. This framework is described in the Supplementary Materials in detail. To evaluate the prognostic performance of the developed models, the concordance index (ci) was calculated (31). While ci = 0.5 is obtained for a noninformative model, ci = 1.0 represents a perfectly predicting model. To compare the performance between nested multivariable Cox models, the likelihood ratio test was applied. The framework to determine gene signatures and corresponding prognostic models was implemented in R Statistics version 3.3.2 (R Foundation for Statistical Computing, Vienna; refs. 32 and 33) and Python (Python Software Foundation. Python Language Reference, version 2.7). An overview of the used programs and packages is given in Supplementary Table S2. For all analyses, two-sided tests were performed and P values below 0.05 were considered as statistically significant.

Statistical framework to identify gene signatures and perform model predictions

The statistical framework to identify gene signatures consists of four main steps, which are described in detail in Supplementary sections S1–S3: (i) the gene expression data are preprocessed. Genes are removed from analysis, if their median expression is below twice the median negative control in the training cohort. The expression of each gene is z-transformed to mean 0 and SD 1 on the training cohort, which is favorable for most machine-learning algorithms. The gene expressions of the validation cohort are transformed on the basis of the means and SDs of the training cohort; (ii) a feature selection method is applied to select the most relevant genes using internal 3-fold cross validation on the training cohort (repeated 333 times). The genes are combined to an ensemble signature based on their frequency of occurrence and their importance (Fig. 2). The resulting signature is then used to build prognostic models on 1,000 bootstrap samples of the training cohort to predict the considered outcome. Several feature selection methods, prognostic models, and different signature sizes (1–10) are compared and the best signature is chosen using the out of the bag data of the bootstrap samples; (iii) to increase the robustness of the signature, genes which are highly correlated to one of the signature genes in the training cohort are combined with this gene to create a new metagene (median expression of the highly correlated genes). The resulting metagene replaces the original gene within the gene signature; (iv) finally, the model is validated using the independent validation cohort. The 95% confidence interval (CI) of the ci is estimated from 1,000 bootstrap samples of the validation cohort. Finally, the validation is declared successful if the 95% CI does not contain 0.5.

Figure 2.
Figure 2.

Cross-validation scheme for identifying the ensemble gene signature. The training cohort was randomly split into 3 equal parts. Each part was used for internal validation and the remaining patients for internal training. This was repeated 333 times. Feature selection was performed on each internal training sample, and a prognostic model was trained using the selected genes. This model was subsequently internally validated. Finally, the occurrence and importance of the genes as well as the validation ci of all cross-validation experiments were used to define the ensemble gene signature.

Data and materials availability

The final models and the raw genomic data used for creating the models are available upon request.

Results

Patient cohorts

In this retrospective study, a multicenter training cohort of 130 patients and an independent, monocenter validation cohort of 121 patients with HPV16 DNA–negative locally advanced HNSCC were available for the development of a gene signature to predict the clinical endpoints LRC, OS, and DM. Patient data, treatment parameters, and tumor characteristics of both patient cohorts were published previously (9, 34) and are summarized in Table 1. Patients in the validation cohort were treated with PORT (n = 90) or PORT-C (n = 31), while all patients of the training cohort received PORT-C as the standard treatment. The training cohort included 44.6% patients with oropharyngeal and 37.7% patients with oral cavity carcinomas. In the validation cohort, 21.5% of the patients have been diagnosed with oropharyngeal and 62.0% with oral cavity carcinomas. Patients in the validation cohort showed lower LRC (statistical trend) and OS, while the incidence of DM was not significantly different. Actuarial rates of LRC, freedom from DM, and OS two years after radiotherapy for the training and validation cohort were 83.8% versus 75.0% (P = 0.096), 79.0% versus 82.8% (P = 0.72), and 76.4% versus 64.9% (P = 0.042), respectively.

View this table:
Table 1.

Patient characteristics for the training and validation cohort

7-Gene signature predicts LRC for HPV16 DNA–negative tumors

To identify a prognostic gene signature for the primary endpoint LRC, the four steps (i)–(iv) of the statistical framework outlined in Materials and Methods were performed.

During the preprocessing step (i), the genes FGFR2, PROM1, and TAF7L were removed from the analysis. The mean validation ci from the 3-fold internal cross validation of step (ii) ranged between 0.57 and 0.68 and was similar between different feature selection methods and statistical models (see Supplementary Fig. S1 for signature size 4). An ensemble gene signature was determined for each combination of feature selection method and prognostic model, as described in Supplementary Section S4. The performance of these signatures was evaluated using 1,000 bootstrap samples of the whole training cohort (see Fig. 3). The highest mean ci of 0.78 was obtained for a signature, which contained the genes SERPINE1, CD24, HILPDA, and TCF3. For the final prediction model, Cox regression was chosen, as it is the most simple of the well performing models. Signature size 4 was chosen based on the mean ci and the signature score of the genes (Supplementary Section S4; Supplementary Fig. S2). The signature score was highest for SERPINE1, followed by CD24, HILPDA, and TCF3, which showed a similar score (see Supplementary Fig. S3). To improve the robustness of the identified 4-gene ensemble signature, it was extended by other highly correlated genes in step (iii), as described in Supplementary Section S4. INHBA, ACTN1, and P4HA2 were found to be highly correlated with SERPINE1, while for CD24, TCF3, and HILPDA no additional correlated genes were found (Supplementary Table S3). Thus, our final 7-gene signature for LRC consisted of the genes SERPINE1, INHBA, ACTN1, P4HA2, CD24, TCF3, and HILPDA. For evaluation, the median of the z-transformed, reference-gene normalized expression of SERPINE1, INHBA, ACTN1, and P4HA2 was considered as a new metagene variable (Supplementary Section S4). The whole training cohort was used to fit the final Cox regression model, leading to a training ci of 0.81 (95% CI, 0.75–0.88). The resulting model parameters are shown in Table 2. In the last step (iv), the validation of the final model was performed on the independent validation cohort. A validation ci of 0.69 (0.60–0.77) was obtained, which represents a successful validation of the gene signature for the endpoint LRC.

Figure 3.
Figure 3.

Performance of ensemble gene signatures for locoregional tumor control on the training cohort. For each combination of feature selection algorithm and statistical model, the mean out-of-the-bag (oob) validation ci of the training cohort and its 95% CI are shown. Performance for the endpoint locoregional tumor control was estimated using 1,000 bootstrap samples of the entire training cohort with signature size 4.

View this table:
Table 2.

Multivariable Cox regression of locoregional tumor control

Patient stratification into groups of low and high risk of recurrence was performed for the final Cox model depending on the risk score, which is given by the linear predictor of the model. The optimal cut-off was chosen based on the training cohort at a risk score of 0.10. This cutoff led to the highest fraction of patient stratifications with a significant difference in LRC, based on 1,000 bootstrap samples of the training cohort (992/1,000). Patients in the high risk group showed significantly lower LRC than patients in the low risk group, both for the training (P < 0.001) and the validation cohort (P = 0.001; see Supplementary Fig. S4).

Inclusion of clinical parameters to the 7-gene LRC signature

For the training cohort, it was shown that the established clinical parameters tumor localization and ECE status were significantly correlated with LRC or the secondary endpoints (9, 12). Using only these two parameters in a multivariable Cox model resulted in a lower performance [training: ci = 0.61 (0.53–0.74), validation: ci = 0.66 (0.57–0.74)] compared with the 7-gene signature. Finally, a multivariable Cox model including both, the clinical parameters ECE status and tumor localization (oral cavity vs. others) as well as the 7-gene signature, increased the training ci to 0.82 (0.77–0.89) and the validation ci to 0.71 (0.62–0.78; see Table 2). While in training the clinical Cox model was significantly improved by adding the 7-gene signature (P < 0.001), adding the clinical parameters to the 7-gene signature resulted in only small improvements (P = 0.53). The difference in validation ci was not statistically significant. An additional validation was performed using only those patients who received concurrent chemotherapy, leading to similar results (validation ci = 0.72).

The extended model was used to stratify the patients into two risk groups (cutoff = 0.37), leading to highly significant differences in LRC for the training (P < 0.001) and the validation cohort (P < 0.001). The corresponding Kaplan–Meier curves are presented in Fig. 4 together with a heatmap of the signature for the training cohort (see also Supplementary Fig. S5).

Figure 4.
Figure 4.

Patient stratification by the 7-gene signature and clinical parameters for locoregional tumor control. Kaplan–Meier estimates of locoregional tumor control (LRC) are shown for the training cohort (A) and the validation cohort (B). Patients were stratified into a low-risk group (LR) or a high-risk group (HR) by the linear predictor of the multivariable Cox model which included the 7-gene signature and the clinical parameters ECE status and tumor localization. The cut-off risk score (0.37) was determined on the training cohort and applied to the validation cohort. C, Heatmap of the 7-gene signature as well as ECE status (0, light; 1, dark), localization (oral cavity, dark; others, light), risk group (LR, light; HR, dark), and LRC during follow-up (yes, light; no, dark) for the training cohort.

Comparison with models based on CSC markers and hypoxia classifiers

In a previous study, it was shown that the expression of CSC markers and hypoxia-related genes were prognostic in patients with locally advanced HPV16 DNA–negative HNSCC, who were treated by PORT-C (12). These results were validated in ref. 25. Here, the performance of these models was compared with the 7-gene signature. While in ref. 25, the best performing model, consisting of ECE status, tumor localization, CD44>0.2, and the 15-gene hypoxia classifier (17), showed a validation ci of 0.65 (0.54–0.74), the identified 7-gene signature combined with the clinical parameters led to a validation ci of 0.71 (0.62–0.78).

7-Gene signature predicts for secondary endpoints

As secondary endpoints, overall survival (OS) and freedom from distant metastases (DM) were considered. The 7-gene signature determined for LRC, combined with the clinical features ECE status and tumor localization, was trained and validated for OS and DM (Supplementary Tables S4 and S5). For OS, training and validation led to a ci of 0.71 (0.65–0.79) and 0.64 (0.57–0.70), respectively. For DM, the ci was 0.69 (0.64–0.80) for training and 0.63 (0.52–0.73) for validation. Cox models including only the clinical features ECE status and tumor localization led to a validation ci of 0.60 (0.54–0.66) for OS and of 0.61 (0.52–0.71) for DM, respectively. Hence, the 7-gene signature could improve the prognostic performance also for the secondary endpoints OS and DM compared with clinical parameters alone. Validation on the subgroup of patients receiving concurrent chemotherapy led to a higher ci for the 7-gene signature [OS: 0.72 (0.56–0.84), DM: 0.74 (0.52–0.90)]. Kaplan–Meier survival analyses are presented in Supplementary Figs. S6 and S7.

Discussion

The overall aim of this study was to identify and validate a gene signature for the stratification of patients with HPV-negative, locally advanced HNSCC who are treated by PORT-C based on the clinical endpoint LRC. We identified a 7-gene signature, which contains genes from an extended in-house gene set compared with previous work, which showed that patients with HPV-negative tumors could be further stratified by the expression of CSC markers and hypoxia-associated genes (12). In addition to the HPV infection status, CSC marker expression levels and tumor hypoxia-associated genes, this gene set included genes related to DNA repair, cell-cycle regulation, epithelial–mesenchymal transition, proliferation, or invasion. A statistical framework was developed with the objective of identifying a gene signature which accurately and robustly predicts the risk of locoregional failure. The framework contains data preprocessing, internal cross validation, signature selection, model building, and independent validation.

The identified 7-gene signature contained the genes SERPINE1, INHBA, ACTN1, and P4HA2 (which were combined into a single metagene due to high mutual correlation) as well as the genes CD24, TCF3, and HILPDA. SERPINE1 (also known as PAI-1), HILPDA, INHBA, and P4HA2 are being induced by the hypoxia-inducible factor HIF1 leading to extracellular matrix remodeling (35–37). SERPINE1 plays a role in enhanced migration and cell proliferation as well as decreased cisplatinum-induced apoptosis (38, 39). In a prospective clinical study including 190 patients, high expression of SERPINE1 has been shown to be associated with poor local recurrence-free, progression-free, and cancer-specific survival (38). In a panel of head and neck xenograft tumors, SERPINE1 expression levels were overexpressed prior to treatment mainly in hypoxic tumors (40). After fractionated irradiation, a correlation between SERPINE1 expression levels and local tumor control was found in vivo (40). In addition, mild hypoxia has been shown to induce SERPINE1 expression via the hypoxia-inducible factor HIF1 (37). SERPINE1 is functionally associated with INHBA (41) and ACTN1 (42). In the hypoxia-associated signatures, HILPDA (also known as HIG2) and P4HA2 have also been included (17, 43). HILPDA has been shown to promote proliferation and invasion (44). In the literature, conflicting data exist for CD24, which is expressed in different tumor entities such as breast cancer and cervical cancer and has shown to be associated with increased tumor growth and progression (45). CD24 has also been shown to be involved in cisplatinum resistance (46) and a shortened progression-free survival was observed for several tumor entities with higher expression (45, 47). In contrast, CD24 overexpression has been shown to be correlated with better survival in patients with oral carcinoma (48). They further showed that CD24−/− mice are able to develop progressive oral cancer. Lack of the surface protein CD24 resulted in the expansion of a highly immunosuppressive CD11b+Gr1+ myeloid cell population leading to oral cancer progression. To date, very little is known about the transcription factor 3 (TCF3) and its potential role in cancer. TCF3 is a critical cell signaling molecule (49) and has been shown to promote cell migration and wound repair (50). In contrast, TCF3 was found to be a cell-intrinsic inhibitor of pluripotent self-renewal through limiting the steady-state levels of self-renewal factors such as Oct-4, Sox2, and Nanog in mouse embryonic stem cells (51). Lack of TCF3 leads to increasing levels of Nanog and other self-renewal genes, minimizing the response to differentiation stimuli (51). According to the factors of the final Cox model, overexpression of SERPINE1 (as well as the highly correlated genes INHBA, ACTN1, and P4HA2) and HILPDA increased the risk for locoregional failure, which is in line with the literature (38, 39). In contrast, a high expression of CD24 led to decreased risk of recurrence. For oral cavity cancer it has been shown that CD24 dampens the functional expansion of myeloid-derived suppressor cells and gives rise to a more favorable prognosis as described above (48), which is in line with our findings. The final Cox model also predicted that a high expression of TCF3 is related to improved LRC, which may be due to its role in the suppression of self-renewal genes (51). However, the functional role of TCF3 in HNSCC needs to be explored in further mechanistic studies.

The identified 7-gene signature showed a good prognostic ability for the endpoint LRC on the validation cohort (ci = 0.69). When combined with the clinical parameters ECE status and tumor localization, its performance could be further improved (ci = 0.71). This indicates that the combination of well-established clinical parameters and prognostic biomarkers may lead to a more accurate prognosis than each of them alone. The model including only the clinical parameters showed the lowest validation performance (ci = 0.66). In the Cox model combining clinical parameters with the 7-gene signature, most signature genes were significantly associated with LRC. While this may be expected on the training cohort, the impact of the 7-gene signature in validation is less clear, as the relevant improvement in ci by 0.05 was not statistically significant. Evaluating the signature combined with HPV status for all patients increased the validation ci to 0.74, which is similar to or even higher than in other studies (30, 52).

The final Cox model showed a better performance on the training cohort (ci=0.82) than on the validation cohort (ci=0.71). This difference is expected, since the final Cox model is adjusted to the training cohort and potential overfitting might occur. In addition, the validation of the proposed 7-gene signature might be impeded by the significant differences between both patient cohorts. Patients in the validation cohort were clinically characterized by a higher percentage of prognostically favorable R0 resections of primary tumors and less lymph nodes with ECE. On the other hand, the validation cohort had a higher percentage of prognostically unfavorable oral cavity tumors, much less concurrent chemotherapy (31/121) than the training cohort and was treated with outdated radiation technologies (53). These negative prognostic factors outbalanced the positive ones resulting in worse outcome in terms of LRC (P = 0.096) and OS (P = 0.042; ref. 25). Lack of concurrent chemotherapy may impede the validation of genes related to cisplatin resistance. For the 7-gene signature, however, only CD24 has been reported to be strongly involved in resistance to cisplatin (46), but also in other mechanisms (48).

On the basis of the final Cox model, a risk score was calculated for each patient, which allowed stratification into groups of low and high risk of recurrence. However, mean gene expressions (Supplementary Table S6) as well as clinical parameters were significantly different between the training and validation cohort. These differences caused a shift in the risk score, such that the stratification cutoff, which was based on the training cohort, led to imbalanced patient risk groups for the validation cohort. While in training approximately 45% of the patients were stratified in the low-risk group and 55% in the high-risk group, for the validation cohort only about 12% of the patients were classified as high risk. Such imbalances may be caused by the differing tumor and treatment characteristics between the cohorts. In addition to clinical reasons, differences in gene expression might also be caused by several biomaterial-related factors such as storage time of FFPE material (3–18 years) or batch effects and stability of reagents and consumables (Supplementary Table S7). Renormalizing the validation data to the training data, as described in refs. 12 and 54, gives the same fraction of patients in the low- and high-risk group and similar LRC rates for both cohorts (Supplementary Fig. S5). However, to apply this renormalization method for individual patient prognosis within clinical trials, the inclusion of reference samples may be required, for which the expected gene expression levels are known. This methodology will be applied to the planned prospective validation of the 7-gene signature. In addition, the application of broadly available and cost-effective PCR-based methods may further improve biomarker stability.

In this study, several algorithms for gene selection and risk prediction were compared. Feature selection algorithms based on mutual information, such as MIFS and MRMR, typically led to a higher ci than simple univariable methods such as Pearson or Spearman correlations (Fig. 3). This behavior can be expected, as the more complex algorithms do not only account for the correlation of the gene expressions to outcome but also consider correlations between the selected genes. Therefore, each gene in the signature represents additional information, which increases the performance of the signature. The performance of prediction models, ranging from the well-known Cox model to complex random forests, was similar on the training cohort. Therefore, the performance of the signature was finally assessed by multivariable Cox regression, which allows easy interpretation. Most of the considered models require additional hyperparameters, such as the regularization parameters Embedded Image and Embedded Image for penalized Cox models or node size and node depth for random forests (see Supplementary Section S3). In an initial experiment, these parameters were chosen based on their default values given in the used software packages and then tuned by a grid search using 2-fold internal cross validation on the training cohort. The resulting parameters were applied in this study and are reported in Supplementary Table S8. While random forests did not outperform simple Cox regression in this study, this may not hold in other situations (55).

The presented 7-gene signature was identified for patients with HPV16 DNA–negative tumors and the primary endpoint LRC. However, it also improved the prognostic value of the clinical parameters for the secondary endpoint OS, while for DM, no significant difference was observed. In particular, for patients receiving concurrent chemotherapy, the validation performance of the 7-gene signature was improved by 10%. This may further enhance the clinical potential of this signature.

A limitation of this study might also be the limited number of genes contained in the initial gene set. Although this has been composed on a hypothesis-driven basis and comprehensive literature search, it may not include all genes of radiobiological relevance. For example CD44, which has been shown to be a prognosticator for LRC in patients with locally advanced HNSCC who received PORT-C (12), had to be omitted from the NanoString analysis due to incorrect probe design. Since the set-up of our gene set, other genes have been shown to be prognostic for outcome in HNSCC. For example, TCGA analyses (56) suggested several genes, related to HPV status. Of these genes CCND1, NOTCH1, YAP1, and SOX2 were found to overlap with our gene set. In the TCGA dataset, patients with CCND1-overexpressing tumors, who received surgery with or without postoperative radiochemotherapy showed worse prognosis. In our study, CCND1 had no impact on the primary endpoint LRC (P = 0.72). Therefore, it was not selected in the gene signature. However, CCND1 showed a significant correlation to the secondary endpoints OS and DM using univariable Cox regression for all 195 patients. For the subgroup of patients with HPV-negative tumors, CCND1 neither correlated with OS nor with DM. This could be explained by the strong correlation of CCND1 with the HPV status in our cohort. In contrast, YAP1 was significantly associated with LRC in our study, but was rated only at rank 14 such that it was not included in the 7-gene signature. NOTCH1 and SOX2 were not related to LRC. Another example is PD-L1, which was strongly associated with local failure in HPV-negative HNSCC (13, 57), but not included in our gene set. To consider these novel developments and identify further biomarkers, whole transcriptome analyses supplemented by whole methylome analyses might be performed and potentially further improve patient stratification.

Currently, an adaptive clinical biomarker matrix trial is set-up within the DKTK-ROG for dose escalation and deescalation in HNSCC. In the first stage, patients with HPV-positive tumors treated by PORT-C will receive a 10% lower radiation dose of the standard concurrent radiochemotherapy schedule. In the second stage, the 7-gene signature is one candidate biomarker for selecting patients with high-risk HPV–negative tumors for dose escalation. To reduce toxicities, especially at higher doses, proton therapy will be considered (58).

In conclusion, this study introduces a novel 7-gene signature predicting LRC for patients with locally advanced HNSCC treated by PORT-C. A prognostic Cox model was trained on a large multicenter patient cohort and independently validated. Although the validation cohort differed in many aspects from the training cohort, a successful validation was achieved, which indicates the robustness of the signature. Prospective validation of the signature is planned within an ongoing prospective clinical trial of the DKTK-ROG before regular application in clinical trials for patient stratification.

Disclosure of Potential Conflicts of Interest

V. Gudziol and V. Budach are consultant/advisory board members for Bristol-Myers Squibb. I. Tinhofer is a consultant/advisory board member for Merck Serono. C. Belka reports receiving commercial research grants from and is a consultant/advisory board member for Merck Darmstadt. S.E. Combs reports receiving speakers bureau honoraria from Acccuray and Brainlab, and is a consultant/advisory board member for Bristol-Myers Squibb and Daiichi Sankyo. No potential conflicts of interest were disclosed by the other authors.

Authors' Contributions

Conception and design: A. Linge, F. Lohaus, A. Abdollahi, M. Baumann, M. Krause, S. Löck

Development of methodology: S. Schmidt, A. Linge, A. Zwanenburg, F. Lohaus, S. Löck

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): A. Linge, F. Lohaus, C. Krenn, S. Appold, V. Gudziol, A. Nowak, C. von Neubeck, I. Tinhofer, V. Budach, A. Sak, M. Stuschke, P. Balermpas, C. Rödel, H. Bunea, A.-L. Grosu, A. Abdollahi, J. Debus, U. Ganswindt, C. Belka, S. Pigorsch, S.E. Combs, D. Mönnich, D. Zips, G.B. Baretton

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): S. Schmidt, A. Linge, A. Zwanenburg, S. Leger, C. Belka, M. Baumann, S. Löck

Writing, review, and/or revision of the manuscript: S. Schmidt, A. Linge, A. Zwanenburg, S. Leger, F. Lohaus, C. Krenn, S. Appold, V. Gudziol, A. Nowak, I. Tinhofer, V. Budach, P. Balermpas, C. Rödel, H. Bunea, A.L. Grosu, J. Debus, U. Ganswindt, C. Belka, S. Pigorsch, S.E. Combs, D. Mönnich, D. Zips, G.B. Baretton, M. Baumann, M. Krause, S. Löck

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): A. Linge, A. Zwanenburg, J. Debus, G.B. Baretton, F. Buchholz, S. Löck

Study supervision: J. Debus, S.U. Pigorsch, M. Baumann, M. Krause, S. Löck

Acknowledgments

The study was financed by a Joint Funding Grant within the German Cancer Consortium (DKTK) awarded to the DKTK-ROG (principal investigator: M. Baumann). The DKTK is funded as one of the National German Health Centers by the Federal German Ministry of Education and Research. The authors gratefully acknowledge the excellent technical assistance by Mrs. Sigrid Balschukat, Mrs. Daniela Friede, and Mrs. Liane Stolz-Kieslich. The authors thank all pathologists, head and neck surgeons, and maxillofacial surgeons at the treatment centers who provided tumor material and data for this study.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Footnotes

  • Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).

  • Received August 11, 2017.
  • Revision received November 4, 2017.
  • Accepted December 29, 2017.
  • Published first January 3, 2018.

References

  1. 1.↵
    1. Duray A,
    2. Descamps G,
    3. Decaestecker C,
    4. Remmelink M,
    5. Sirtaine N,
    6. Lechien J,
    7. et al.
    Human papillomavirus DNA strongly correlates with a poorer prognosis in oral cavity carcinoma. Laryngoscope 2012;122:1558–65.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Leemans CR,
    2. Braakhuis BJM,
    3. Brakenhoff RH
    . The molecular biology of head and neck cancer. Nat Rev Cancer 2011;11:9–22.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Fietkau R,
    2. Lautenschläger C,
    3. Sauer R,
    4. Dunst J,
    5. Becker A,
    6. Baumann M,
    7. et al.
    Postoperative concurrent radiochemotherapy versus radiotherapy in high-risk SCCA of the head and neck: results of the German phase III trial ARO 96–3. J Clin Oncol 2006;24:5507–5507.
    OpenUrl
  4. 4.↵
    1. Cooper JS,
    2. Pajak TF,
    3. Forastiere AA,
    4. Jacobs J,
    5. Campbell BH,
    6. Saxman SB,
    7. et al.
    Postoperative concurrent radiotherapy and chemotherapy for high-risk squamous-cell carcinoma of the Head and Neck. N Engl J Med 2004;350:1937–44.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Bernier J,
    2. Domenge C,
    3. Ozsahin M,
    4. Matuszewska K,
    5. Lefèbvre J-L,
    6. Greiner RH,
    7. et al.
    Postoperative irradiation with or without concomitant chemotherapy for locally advanced Head and Neck Cancer. N Engl J Med 2004;350:1945–52.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Marur S,
    2. D'Souza G,
    3. Westra WH,
    4. Forastiere AA
    . HPV-associated head and neck cancer: a virus-related cancer epidemic. Lancet Oncol 2010;11:781–9.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Kimple RJ,
    2. Smith MA,
    3. Blitzer GC,
    4. Torres AD,
    5. Martin JA,
    6. Yang RZ,
    7. et al.
    Enhanced radiation sensitivity in HPV-positive head and neck cancer. Cancer Res 2013;73:4791–800.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Gupta AK,
    2. Lee JH,
    3. Wilke WW,
    4. Quon H,
    5. Smith G,
    6. Maity A,
    7. et al.
    Radiation response in two HPV-infected head-and-neck cancer cell lines in comparison to a non–HPV-Infected cell line and relationship to signaling through AKT. Int J Radiat Oncol 2009;74:928–33.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Lohaus F,
    2. Linge A,
    3. Tinhofer I,
    4. Budach V,
    5. Gkika E,
    6. Stuschke M,
    7. et al.
    HPV16 DNA status is a strong prognosticator of loco-regional control after postoperative radiochemotherapy of locally advanced oropharyngeal carcinoma: results from a multicentre explorative study of the German cancer consortium radiation oncology group. Radiother Oncol 2014;113:317–23.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Balermpas P,
    2. Rödel F,
    3. Rödel C,
    4. Krause M,
    5. Linge A,
    6. Lohaus F,
    7. et al.
    CD8+ tumour-infiltrating lymphocytes in relation to HPV status and clinical outcome in patients with head and neck cancer after postoperative chemoradiotherapy: A multicentre study of the German cancer consortium radiation oncology group (DKTK-ROG). Int J Cancer 2016;138:171–81.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Tinhofer I,
    2. Budach V,
    3. Saki M,
    4. Konschak R,
    5. Niehr F,
    6. Jöhrens K,
    7. et al.
    Targeted next-generation sequencing of locally advanced squamous cell carcinomas of the head and neck reveals druggable targets for improving adjuvant chemoradiation. Eur J Cancer 2016;57:78–86.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Linge A,
    2. Löck S,
    3. Gudziol V,
    4. Nowak A,
    5. Lohaus F,
    6. von Neubeck C,
    7. et al.
    Low cancer stem cell marker expression and low hypoxia identify good prognosis subgroups in HPV(-) HNSCC after postoperative radiochemotherapy: a multicenter study of the DKTK-ROG. Clin Cancer Res 2016;22:2639–49.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Balermpas P,
    2. Rödel F,
    3. Krause M,
    4. Linge A,
    5. Lohaus F,
    6. Baumann M,
    7. et al.
    The PD-1/PD-L1 axis and human papilloma virus in patients with head and neck cancer after adjuvant chemoradiotherapy: a multicentre study of the German cancer consortium radiation oncology group (DKTK-ROG). Int J Cancer 2017;141:594–603.
    OpenUrl
  14. 14.↵
    1. Overgaard J.
    Hypoxic radiosensitization: adored and ignored. J Clin Oncol 2007;25:4066–74.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Nordsmark M,
    2. Overgaard M,
    3. Overgaard J
    . Pretreatment oxygenation predicts radiation response in advanced squamous cell carcinoma of the head and neck. Radiother Oncol 1996;41:31–9.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Nordsmark M,
    2. Bentzen SM,
    3. Rudat V,
    4. Brizel D,
    5. Lartigau E,
    6. Stadler P,
    7. et al.
    Prognostic value of tumor oxygenation in 397 head and neck tumors after primary radiation therapy. An international multi-center study. Radiother Oncol 2005;77:18–24.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Toustrup K,
    2. Sørensen BS,
    3. Nordsmark M,
    4. Busk M,
    5. Wiuf C,
    6. Alsner J,
    7. et al.
    Development of a hypoxia gene expression classifier with predictive impact for hypoxic modification of radiotherapy in head and neck cancer. Cancer Res 2011;71:5923–31.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Buffa FM,
    2. Harris AL,
    3. West CM,
    4. Miller CJ
    . Large meta-analysis of multiple cancers reveals a common, compact and highly prognostic hypoxia metagene. Br J Cancer 2010;102:428–35.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Yaromina A,
    2. Thames H,
    3. Zhou X,
    4. Hering S,
    5. Eicheler W,
    6. Dörfler A,
    7. et al.
    Radiobiological hypoxia, histological parameters of tumour microenvironment and local tumour control after fractionated irradiation. Radiother Oncol 2010;96:116–22.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Bütof R,
    2. Dubrovska A,
    3. Baumann M
    . Clinical perspectives of cancer stem cell research in radiation oncology. Radiother Oncol 2013;108:388–96.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Peitzsch C,
    2. Perrin R,
    3. Hill RP,
    4. Dubrovska A,
    5. Kurth I
    . Hypoxia as a biomarker for radioresistant cancer stem cells. Int J Radiat Biol 2014;90:636–52.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Baumann M,
    2. Krause M,
    3. Hill R
    . Exploring the role of cancer stem cells in radioresistance. Nat Rev Cancer 2008;8:545–54.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Yaromina A,
    2. Krause M,
    3. Baumann M
    . Individualization of cancer treatment from radiotherapy perspective. Mol Oncol 2012;6:211–21.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Wu Y,
    2. Posner MR,
    3. Schumaker LM,
    4. Nikitakis N,
    5. Goloubeva O,
    6. Tan M,
    7. et al.
    Novel biomarker panel predicts prognosis in human papillomavirus-negative oropharyngeal cancer. Cancer 2012;118:1811–7.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Linge A,
    2. Löck S,
    3. Krenn C,
    4. Appold S,
    5. Lohaus F,
    6. Nowak A,
    7. et al.
    Independent validation of the prognostic value of cancer stem cell marker expression and hypoxia-induced gene expression for patients with locally advanced HNSCC after postoperative radiotherapy. Clin Transl Radiat Oncol 2016;1:19–26.
    OpenUrl
  26. 26.↵
    1. De Cecco L,
    2. Nicolau M,
    3. Giannoccaro M,
    4. Daidone MG,
    5. Bossi P,
    6. Locati L,
    7. et al.
    Head and neck cancer subtypes with biological and clinical relevance: Meta-analysis of gene-expression data. Oncotarget 2015;6:9627–42.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. van den Broek GB,
    2. Wildeman M,
    3. Rasch CRN,
    4. Armstrong N,
    5. Schuuring E,
    6. Begg AC,
    7. et al.
    Molecular markers predict outcome in squamous cell carcinoma of the head and neck after concomitant cisplatin-based chemoradiation. Int J Cancer 2009;124:2643–50.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Agada FO,
    2. Patmore H,
    3. Alhamarneh O,
    4. Stafford ND,
    5. Greenman J
    . Genetic profile of head and neck squamous cell carcinoma: clinical implications. J Laryngol Otol 2009;123:266–72.
    OpenUrlPubMed
  29. 29.↵
    1. Ginos MA,
    2. Page GP,
    3. Michalowicz BS,
    4. Patel KJ,
    5. Volker SE,
    6. Pambuccian SE,
    7. et al.
    Identification of a gene expression signature associated with recurrent disease in squamous cell carcinoma of the head and neck. Cancer Res 2004;64:55–63.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. De Cecco L,
    2. Bossi P,
    3. Locati L,
    4. Canevari S,
    5. Licitra L
    . Comprehensive gene expression meta-analysis of head and neck squamous cell carcinoma microarray data defines a robust survival predictor. Ann Oncol 2014;25:1628–35.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Harrell FE
    . Regression modeling strategies. 1st edition. New York, NY: Springer; 2001.
  32. 32.
    R Development Core Team. R: a language and environment for statistical computing [Internet]. Vienna, Austria; 2008. Available from: http://www.r-project.org.
  33. 33.
    RStudio Team. RStudio: integrated development environment for R [Internet]. Boston, MA; 2015. Available from: http://www.rstudio.com/.
  34. 34.↵
    1. Linge A,
    2. Lohaus F,
    3. Löck S,
    4. Nowak A,
    5. Gudziol V,
    6. Valentini C,
    7. et al.
    HPV status, cancer stem cell marker expression, hypoxia gene signatures and tumour volume identify good prognosis subgroups in patients with HNSCC after primary radiochemotherapy: a multicentre retrospective study of the DKTK-ROG. Radiother Oncol 2016;121:364–73.
    OpenUrl
  35. 35.↵
    1. Gilkes DM,
    2. Bajpai S,
    3. Chaturvedi P,
    4. Wirtz D,
    5. Semenza GL
    . Hypoxia-inducible Factor 1 (HIF-1) promotes extracellular matrix remodeling under hypoxic conditions by inducing P4HA1, P4HA2, and PLOD2 expression in fibroblasts. J Biol Chem 2013;288:10819–29.
    OpenUrlAbstract/FREE Full Text
  36. 36.↵
    1. Depoix CL,
    2. de Selliers I,
    3. Hubinont C,
    4. Debieve F
    . HIF1A and EPAS1 potentiate hypoxia-induced upregulation of inhibin alpha chain expression in human term cytotrophoblasts in vitro. Mol Hum Reprod 2017;23:199–209.
    OpenUrl
  37. 37.↵
    1. Kietzmann T,
    2. Roth U,
    3. Jungermann K
    . Induction of the plasminogen activator inhibitor-1 gene expression by mild hypoxia via a hypoxia response element binding the hypoxia- inducible factor-1 in rat hepatocytes. Am Soc Hematol 1999;94:4177–85.
    OpenUrl
  38. 38.↵
    1. Pavón MA,
    2. Arroyo-Solera I,
    3. Téllez-Gabriel M,
    4. León X,
    5. Virós D,
    6. López M,
    7. et al.
    Enhanced cell migration and apoptosis resistance may underlie the association between high SERPINE1 expression and poor outcome in head and neck carcinoma patients. Oncotarget 2015;6:29016–33.
    OpenUrl
  39. 39.↵
    1. Pavón MA,
    2. Arroyo-Solera I,
    3. Céspedes MV,
    4. Casanova I,
    5. León X,
    6. Mangues R
    . uPA/uPAR and SERPINE1 in head and neck cancer: role in tumor resistance, metastasis, prognosis and therapy. Oncotarget 2016;7:57351–66.
    OpenUrl
  40. 40.↵
    1. Bayer C,
    2. Schilling D,
    3. Hoetzel J,
    4. Egermann HP,
    5. Zips D,
    6. Yaromina A,
    7. et al.
    PAI-1 levels predict response to fractionated irradiation in 10 human squamous cell carcinoma lines of the head and neck. Radiother Oncol 2008;86:361–8.
    OpenUrlPubMed
  41. 41.↵
    1. Wang Y-H,
    2. Tian X,
    3. Liu O-S,
    4. Fang X-D,
    5. Quan H-Z,
    6. Xie S,
    7. et al.
    Gene profiling analysis for patients with oral verrucous carcinoma and oral squamous cell carcinoma. Int J Clin Exp Med 2014;7:1845–52.
    OpenUrl
  42. 42.↵
    1. Menicanin D,
    2. Bartold PM,
    3. Zannettino ACW,
    4. Gronthos S
    . Genomic profiling of mesenchymal stem cells. Stem Cell Rev Reports 2009;5:36–50.
    OpenUrl
  43. 43.↵
    1. Eustace A,
    2. Mani N,
    3. Span PN,
    4. Irlam JJ,
    5. Taylor J,
    6. Betts GNJ,
    7. et al.
    A 26-gene hypoxia signature predicts benefit from hypoxia-modifying therapy in laryngeal cancer but not bladder cancer. Clin Cancer Res 2013;19:4879–88.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Mao X,
    2. Wang C,
    3. Liu D,
    4. Zhang X,
    5. Wang L,
    6. Yan M,
    7. et al.
    Hypoxia upregulates HIG2 expression and contributes to bevacizumab resistance in glioblastoma. Oncotarget 2016;7:47808–20.
    OpenUrl
  45. 45.↵
    1. Lee J-H,
    2. Kim S-H,
    3. Lee E-S,
    4. Kim Y-S
    . CD24 overexpression in cancer development and progression: a meta-analysis. Oncol Rep 2009;22:1149–56.
    OpenUrlPubMed
  46. 46.↵
    1. Modur V,
    2. Joshi P,
    3. Nie D,
    4. Robbins KT,
    5. Khan AU,
    6. Rao K
    . CD24 expression may play a role as a predictive indicator and a modulator of cisplatin treatment response in head and neck squamous cellular carcinoma. PLoS One 2016;11:e0156651.
    OpenUrl
  47. 47.↵
    1. Smith SC
    . The metastasis-associated gene CD24 Is Regulated by Ral GTPase and is a mediator of cell proliferation and survival in human cancer. Cancer Res 2006;66:1917–22.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Fugle CW,
    2. Zhang Y,
    3. Hong F,
    4. Sun S,
    5. Westwater C,
    6. Rachidi S,
    7. et al.
    CD24 blunts oral squamous cancer development and dampens the functional expansion of myeloid-derived suppressor cells. Oncoimmunology 2016;5:e1226719.
    OpenUrl
  49. 49.↵
    1. Zhang Z,
    2. Filho MSS,
    3. Nör JE
    . The biology of head and neck cancer stem cells. Oral Oncol 2012;48:1–9.
    OpenUrlPubMed
  50. 50.↵
    1. Miao Q,
    2. Ku AT,
    3. Nishino Y,
    4. Howard JM,
    5. Rao AS,
    6. Shaver TM,
    7. et al.
    Tcf3 promotes cell migration and wound repair through regulation of lipocalin 2. Nat Commun 2014;5:4088.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Yi F,
    2. Pereira L,
    3. Merrill BJ
    . Tcf3 functions as a steady-state limiter of transcriptional programs of mouse embryonic stem cell self-renewal. Stem Cells 2008;26:1951–60.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Aerts HJWL,
    2. Velazquez ER,
    3. Leijenaar RTH,
    4. Parmar C,
    5. Grossmann P,
    6. Carvalho S,
    7. et al.
    Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun 2014;5:4006.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Baumann M,
    2. Krause M,
    3. Overgaard J,
    4. Debus J,
    5. Bentzen SM,
    6. Daartz J,
    7. et al.
    Radiation oncology in the era of precision medicine. Nat Rev Cancer 2016;16:234–49.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Toustrup K,
    2. Sørensen BS,
    3. Metwally MAH,
    4. Tramm T,
    5. Mortensen LS,
    6. Overgaard J,
    7. et al.
    Validation of a 15-gene hypoxia classifier in head and neck cancer for prospective use in clinical trials. Acta Oncol 2016;55:1091–8.
    OpenUrl
  55. 55.↵
    1. Zheng G.
    Identification of candidate colon cancer biomarkers by applying a random forest approach on microarray data. Oncol Rep 2012;28:1036–42.
    OpenUrlPubMed
  56. 56.↵
    Cancer Genome Atlas Network. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature 2015;517:576–82.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Skinner HD,
    2. Giri U,
    3. Yang LP,
    4. Kumar M,
    5. Liu Y,
    6. Story MD,
    7. et al.
    Integrative analysis identifies a novel AXL-PI3 kinase-PD-L1 signaling axis associated with radiation resistance in head and neck cancer. Clin Cancer Res 2017;23:2713–22.
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    1. Jakobi A,
    2. Bandurska-Luque A,
    3. Stützer K,
    4. Haase R,
    5. Löck S,
    6. Wack L-J,
    7. et al.
    Identification of patient benefit from proton therapy for advanced head and neck cancer patients based on individual and subgroup normal tissue complication probability analysis. Int J Radiat Oncol Biol Phys 2015;92:1165–74.
    OpenUrlCrossRefPubMed
View Abstract
Back to top
View this article with LENS

Open full page PDF
Article Alerts
Email Article
Citation Tools
Development and Validation of a Gene Signature for Patients with Head and Neck Carcinomas Treated by Postoperative Radio(chemo)therapy
Stefan Schmidt, Annett Linge, Alex Zwanenburg, Stefan Leger, Fabian Lohaus, Constanze Krenn, Steffen Appold, Volker Gudziol, Alexander Nowak, Cläre von Neubeck, Inge Tinhofer, Volker Budach, Ali Sak, Martin Stuschke, Panagiotis Balermpas, Claus Rödel, Hatice Bunea, Anca-Ligia Grosu, Amir Abdollahi, Jürgen Debus, Ute Ganswindt, Claus Belka, Steffi Pigorsch, Stephanie E. Combs, David Mönnich, Daniel Zips, Gustavo B. Baretton, Frank Buchholz, Michael Baumann, Mechthild Krause, Steffen Löck and for the DKTK-ROG
Clin Cancer Res March 15 2018 (24) (6) 1364-1374; DOI: 10.1158/1078-0432.CCR-17-2345
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo

Jump to section

Advertisement