Next Article in Journal
Differential Activation of ERK Signaling in HPV-Related Oropharyngeal Squamous Cell Carcinoma
Next Article in Special Issue
Fibroblast-Mediated Immunoregulation of Macrophage Function Is Maintained after Irradiation
Previous Article in Journal
p53-Dependent Apoptotic Effect of Puromycin via Binding of Ribosomal Protein L5 and L11 to MDM2 and Its Combination Effect with RITA or Doxorubicin
Previous Article in Special Issue
Extracellular Vesicles as Transmitters of Hypoxia Tolerance in Solid Cancers
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 11
Issue 4
10.3390/cancers11040583
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Acute Hypoxia Profile is a Stronger Prognostic Factor than Chronic Hypoxia in Advanced Stage Head and Neck Cancer Patients

by
Martijn van der Heijden
1,2,†,
Monique C. de Jong
3,†,
Caroline V. M. Verhagen
1,2,
Reinout H. de Roest
4,
Sebastian Sanduleanu
5,
Frank Hoebers
5,
C. René Leemans
4,
Ruud H. Brakenhoff
4,
Conchita Vens
1,3,
Marcel Verheij
1,3,6 and
Michiel W. M. van den Brekel
2,7,*
1
Division of Cell Biology, The Netherlands Cancer Institute, Plesmanlaan 121, 1006 CX Amsterdam, The Netherlands
2
Department of Head and Neck Oncology and Surgery, The Netherlands Cancer Institute, Plesmanlaan 121, 1006 CX Amsterdam, The Netherlands
3
Department of Radiation Oncology, The Netherlands Cancer Institute, Plesmanlaan 121, 1006 CX Amsterdam, The Netherlands
4
Amsterdam UMC, Vrije Universiteit Amsterdam, Otolaryngology/Head & Neck Surgery, Cancer Center Amsterdam, De Boelelaan 1117, 1081 HV Amsterdam, The Netherlands
5
Department of Radiation Oncology (MAASTRO), GROW – School for Oncology and Developmental Biology, Maastricht University Medical Centre, Dr. Tanslaan 12, 6229 ET Maastricht, The Netherlands
6
Department of Radiation Oncology, Radboud University Medical Center, Geert Grooteplein Zuid 10, 6525 GA Nijmegen, The Netherlands
7
Department of Oral and Maxillofacial Surgery, Amsterdam UMC, Academic Medical Center, Meibergdreef 9, 1105 AZ Amsterdam, The Netherlands
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cancers 2019, 11(4), 583; https://doi.org/10.3390/cancers11040583
Submission received: 6 March 2019 / Revised: 11 April 2019 / Accepted: 19 April 2019 / Published: 25 April 2019
(This article belongs to the Special Issue Tumor Radioresistance)
Browse Figures
Review Reports Versions Notes

Abstract

:
Hypoxic head and neck tumors respond poorly to radiotherapy and can be identified using gene expression profiles. However, it is unknown whether treatment outcome is driven by acute or chronic hypoxia. Gene expression data of 398 head and neck cancers was collected. Four clinical hypoxia profiles were compared to in vitro acute and chronic hypoxia profiles. Chronic and acute hypoxia profiles were tested for their association to outcome using Cox proportional hazard analyses. In an initial set of 224 patients, scores of the four clinical hypoxia profiles correlated with each other and with chronic hypoxia. However, the acute hypoxia profile showed a stronger association with local recurrence after chemoradiotherapy (p = 0.02; HR = 3.1) than the four clinical (chronic hypoxia) profiles (p = 0.2; HR = 0.9). An independent set of 174 patients confirmed that acute hypoxia is a stronger prognostic factor than chronic hypoxia for overall survival, progression-free survival, local and locoregional control. Multivariable analyses accounting for known prognostic factors substantiate this finding (p = 0.045; p = 0.042; p = 0.018 and p = 0.003, respectively). In conclusion, the four clinical hypoxia profiles are related to chronic hypoxia and not acute hypoxia. The acute hypoxia profile shows a stronger association with patient outcome and should be incorporated into existing prediction models.

1. Introduction

The average overall survival for advanced stage head and neck cancer patients is around 50% [1], but this varies greatly among different groups of patients. Clinical (TNM) staging explains survival variation only partially for these patients [2,3,4]. Human papillomavirus (HPV) positive oropharyngeal tumors represent a distinct subgroup of head and neck squamous cell carcinoma (HNSCC) that is associated with a good prognosis [5]. We have previously shown that the addition of HPV status and a prognostic gene expression profile can improve outcome prediction, suggesting that a substantial part of the survival variation is explained by tumor biology [5].
Hypoxia is one of the most studied biological factors affecting prognosis in HNSCC [6]. Tumor cells can become hypoxic by chronic (diffusion limited) and acute (perfusion limited) mechanisms [7]. Perfusion or diffusion limited, intermittent or cycling hypoxia are alternative terms that have been used to better reflect the mechanisms that result in acute hypoxia as referred to in this study. Both acute and chronic hypoxia can have different effects on tumor cells and their microenvironment. It is unclear whether prognosis is mostly impacted by chronic or acute hypoxia [7]. Because oxygen is essential for DNA-damage upon irradiation, hypoxic tumors respond poorly to radiation [8,9]. Since two third of all HNSCC patients receive radiotherapy as part of their treatment, hypoxia can reduce efficacy of the treatment and impact prognosis [10]. A meta-analysis of clinical trials showed that in vivo modification of the hypoxia status during radiotherapy improves survival of HNSCC patients, demonstrating that hypoxia is an important factor in radioresistance [11]. Likewise, hypoxia-inducible factor is a poor prognostic factor in surgically treated HNSCC patients [12]. Unfortunately, the benefit from hypoxia modification therapy was modest and comes with added toxicity [11]. This led to the hypothesis that only patients with hypoxic tumors profit from such a therapeutic intervention, which was further confirmed in two subsequent studies [13,14]. Initiated in 2014, the NIMRAD study aims to ‘prospectively validate a hypoxia gene signature that can be used in clinical practice to personalize treatment and select appropriate patients for hypoxia modifying treatment’ [15].
A robust method to quantify hypoxia in tumors is required to select patients for hypoxia modifying treatment. Different techniques have been applied to evaluate the level of hypoxia in a tumor and its impact on radiotherapy response [6], including an oxygen-sensitive needle probe inserted into the tumor [16,17,18,19], exogenous immunohistochemical markers (e.g., pimonidazole [20]), endogenous biomarkers (e.g., HIF1-alpha [21,22] or carbonic anhydrase IX [14,21]) and imaging techniques like MRI [23] and PET [24]. None of these techniques are currently used in routine clinical practice because they are too complex or insufficiently accurate. Tumor hypoxia can also be deduced from gene expression profiling and many gene expression profiles for hypoxia have been proposed [13,21,25,26,27,28,29].
In HNSCC, four hypoxia profiles have been validated to be prognostic or predictive [13,25,26,27]. Winter et al. [25] obtained a hypoxia 99-gene expression profile that was validated to be associated with recurrence free survival. Buffa et al. [26] used hypoxia-regulated genes to select co-expressed genes in three HNSCC and five breast cancer studies. The resulting 51-gene profile was validated for its prognostic relevance in four independent datasets. Toustrup et al. [13] generated a profile from in vitro experiments together with gene expression data from 58 head and neck cancer biopsies. The resulting 15-gene expression profile proved to be predictive for the response to nimorazole during radiotherapy in 323 HNSCC patients. Eustace et al. [27] reduced the profile from Buffa et al. to a 26-gene profile that predicted regional control in 157 laryngeal cancer patients treated with radiotherapy. Some authors suggest that the use of multiple hypoxia profiles improves the prognostic value [30]. However, none of these studies made a distinction between acute and chronic hypoxia. Gene expression profiles for acute and chronic hypoxia were generated in vitro. Seigneuric et al. used temporal changes in human epithelial mammary cell lines in response to hypoxia to generate gene expression profiles for acute and chronic hypoxia and showed that acute hypoxia is a prognostic factor in breast cancer [31].
In this study we aimed to improve hypoxia-based HNSCC patient prognostication by evaluating the contribution of chronic and acute hypoxia in gene expression profiles. We selected the four clinical HNSCC hypoxia gene expression profiles which are associated with patient outcome in HNSCC [13,25,26,27]. In an initial set of 224 patients, these gene expression profiles were compared amongst each other and to the in vitro chronic and acute hypoxia gene expression profiles of Seigneuric et al. [31]. The prognostic value of these acute and chronic hypoxia expression profiles was then tested in a new cohort of 91 HNSCC patients and validated in an independent cohort of 174 HPV-negative HNSCC patients, all treated with definitive cisplatin based chemoradiotherapy. Lastly, the prognostic value of acute and chronic hypoxia expression profiles was tested in a multivariable analysis with known prognostic factors in HNSCC.

2. Results

2.1. Few Overlapping Genes in the Four Clinical HNSCC Hypoxia Gene Expression Profiles

We compared the four clinical gene expression profiles of Winter et al., Buffa et al., Toustrup et al. and Eustace et al. to acute and chronic hypoxia profiles of Seigneuric et al. (Figure 1). A more extensive overview of the profiles is presented in Supplementary Materials, Tables S1 and S2. The four clinical profiles consisted of 147 unique genes. Of these, 82% was present in only one of the four profiles. Three genes (2%) were present in all four signatures: Aldolase A(ALDOA), Prolyl 4-Hydroxylase Subunit Alpha 1 (P4HA1) and Solute Carrier Family 2 Member 1 (SLC2A1 a.k.a. GLUT-1). Aldolase A is a glycolytic enzyme, the P4HA1 gene encodes a component of a key enzyme in collagen synthesis and the SLC2A1 (a.k.a. GLUT-1) gene encodes a glucose transporter. None of these three genes were present in the Seigneuric chronic in vitro profile, however 9 genes in this profile were present in at least one of the four clinical profiles. The Seigneuric acute profile had no overlapping genes with any of the clinical profiles or the Seigneuric chronic in vitro profile.

2.2. The Four Clinical HNSCC Hypoxia Gene Expression Profiles are Correlated and Resemble Chronic Hypoxia Response

We tested the conformity of the four clinical hypoxia profiles in 224 HNSCC patients. To this end, we combined gene expression data generated from three HNSCC cohorts (Table 1): Pramana (published), de Jong 1 (expanded expression data set with additional samples that could not be used in the original study) and unpublished gene expression data from a patient cohort as described in de Jong 2 et al. [32,33,34,35,36]. Patient characteristics of the cohorts are available in Supplementary Materials, Tables S3–S5. Scores for the four clinical profiles and the Seigneuric acute and chronic hypoxia profiles were generated for all 224 HNSCC patients. The average Spearman correlation between the scores of the four hypoxia profiles was 0.82 and highly significant (range 0.71–0.90, p-values < 0.0001, Figure 2A), demonstrating that the four gene expression profiles rank patients similarly. All four clinical hypoxia profiles were significantly correlated to the Seigneuric chronic hypoxia profile (correlation 0.60, p < 0.0001, Figure 2A,). The average correlation with Seigneuric acute hypoxia was −0.09 (p = 0.2). A clustering based on all scores showed that the four clinical HNSCC hypoxia profiles cluster together with the Seigneuric chronic profile, whereas no correlation was observed with the Seigneuric acute profile (Figure 2B). Together this strongly suggests that the four clinical profiles represent chronic hypoxia and not acute hypoxia.

2.3. The Acute Hypoxia Profile is Associated with Local Control

As the inclusion criteria of the three cohorts were not comparable, we did not attempt to combine the three cohorts for outcome analyses. The largest cohort with Local Control (LC) data available (Pramana cohort, n = 91) was used to assess the prognostic value of acute and chronic hypoxia [32]. Unfortunately, no clinically validated acute hypoxia profiles are available for HNSCC. We therefore used the Seigneuric acute hypoxia profile as a surrogate marker for acute hypoxia in HNSCC. Since the four clinically validated HNSCC hypoxia profiles were highly correlated to each other and to the Seigneuric Chronic profile (Figure 2A), for each patient the scores of the four clinical profiles were averaged to obtain a joint chronic hypoxia score. The median was used to define “High” and “Low” hypoxia groups.
The acute hypoxia profile was significantly associated with local control (p = 0.02). Patients in the acute hypoxia “High” group showed a higher local recurrence rate (Hazard Ratio (HR): 3.1, 95%CI: 1.1−8.6), compared to patients in the acute hypoxia “Low” group. Patients in chronic hypoxia “High” group tended to have a worse local control rate (p = 0.2, HR = 1.9; Figure 3). Kaplan Meier curves and hazard ratios for the individual profiles are provided in Supplementary Materials, Figure S1. To study the combined effect of the acute hypoxia profile and the chronic hypoxia profile we generated three groups: (1) acute and chronic hypoxia profile both low, (2) acute or chronic hypoxia profile high and (3) acute and chronic hypoxia profile both high (Figure 3C). The comparison in Figure 3C shows that, local control rates are worst when both acute and chronic hypoxia profiles are high (p = 0.04).

2.4. Validation of the Acute Hypoxia Profile as a Prognostic Marker for Outcome in HNSCC

HPV is a major prognostic factor in HNSCC and the Pramana cohort contained HPV-positive and HPV-negative tumors. So, to test whether the prognostic value of acute hypoxia is independent of HPV, we aimed to validate the prognostic relevance of acute hypoxia in a HPV-negative HNSCC cohort. To this end, we collected clinical and gene expression data from a study by Van der Heijden et al. [35,36]. This data derived from pre-treatment patient tumor material of two cohorts of HPV-negative advanced stage HNSCC patients, NKI-CRAD and DESIGN, comprising 174 patients in total. All patients were treated with chemo-radiotherapy and an overview of patient characteristics is available in Supplementary Materials, Table S6. Overall Survival (OS), Progression Free Survival (PFS), Local Control (LC), Locoregional Control (LRC) and Distant Metastasis free survival (DM) data were available for this dataset. As before, scores for all 174 patients were calculated. The Seigneuric acute hypoxia profile was used as a marker for acute hypoxia and the scores of the four clinical HNSCC hypoxia profiles (Winter et al., Buffa et al., Toustrup et al. and Eustace et al.) were averaged to obtain a joint chronic hypoxia score. The distribution of the profiles is shown in Supplementary Materials, Figure S2A. Median splits were used to divide patients in “High” and “Low” groups.
Confirming the findings in the Pramana cohort, the acute hypoxia “High” group had a significantly worse LC (p = 0.006, HR = 3.3; Figure 4A), OS (HR = 1.58, p = 0.023), PFS (HR = 1.66, p = 0.009) and LRC; (HR = 2.4, p = 0.008; Supplementary Materials, Figure S2B). The chronic hypoxia “High” group showed a significant difference in PFS only (HR = 1.52, p = 0.03) and a trend for OS (HR = 1.43, p = 0.075; supplementary Materials, Figure S2C), but no significant difference in LC (p = 0.34; HR = 1.43; Figure 4B) and LRC (HR = 1.59, p = 0.141). We next combined the chronic hypoxia and acute hypoxia groups as described earlier to obtain three groups. When both acute and chronic hypoxia expression profiles are “High”, OS, PFS, LC and LRC rates are significantly worse than when both are low (Figure 4C–E and Supplementary Materials, Figure S3).

2.5. Multivariable Analysis Confirms the Prognostic Value of Acute Hypoxia

To determine whether acute hypoxia has additional prognostic value to known prognostic markers, other than HPV, we performed a multivariable in the Van der Heijden cohort. For 149 patients, tumor volume data were available. Within these 149 patients all available variables were tested in a univariable analysis for their association with OS, PFS, LC, LRC and DM (Supplementary Materials, Table S7). To avoid arbitrary cut-offs, acute hypoxia, chronic hypoxia (the average of the four clinical hypoxia scores) and tumor volume were tested as continuous variables.
Acute hypoxia, gender, tumor site, disease stage, cisplatin dose, and tumor volume showed significant associations with multiple patient outcome measures in univariable analysis (Supplementary Materials, Table S7). Chronic hypoxia only showed associations with locoregional control. Variables with significant associations to patient outcome were combined in a multivariable Cox proportional hazard model (Table 2). Previous research has shown an interaction between hypoxia and tumor volume [37], this interaction was therefore incorporated in the multivariable analysis. Results of the multivariable analysis show that acute hypoxia is significantly associated with OS (HR = 3.5, p = 0.045), PFS (HR = 3.3, p = 0.042), LC (HR = 38.2, p = 0.018) and LRC (HR = 27, p = 0.003), whereas chronic hypoxia is not associated with any of the outcomes. Tumor volume (per cm3 increase) is also associated with OS (HR = 1.01, p = 4.12E-05), PFS (HR = 1.008, p = 0.0007) and LC (HR = 1.014, p = 0.027). The interaction of acute hypoxia with tumor volume was significantly associated with OS (HR = 0.98, p = 0.004), PFS (HR = 0.98, p = 0.02), and LRC (HR = 0.95, p = 0.016). The interaction of chronic hypoxia with tumor volume was not significantly associated with any of the reported outcomes. These analyses confirm the prognostic value of acute hypoxia independent from clinical factors, tumor volume and chronic hypoxia.

3. Discussion

The four clinical gene expression profiles for hypoxia, which have been validated to predict outcome in HNSCC, have few overlapping genes. Nevertheless, they classified patients similarly, indicating that they reflect a similar underlying biological process: chronic and not acute hypoxia. Since the clinical hypoxia profiles were correlated, we combined them in a joint chronic hypoxia score. In contrast to the chronic hypoxia profile, the acute profile was associated with local control in the Pramana cohort. The poor association of the acute hypoxia profile with patient outcome was validated in an independent validation cohort (n = 174) and multivariable analysis showed that acute hypoxia is a significant prognostic factor independent of clinical factors, tumor volume and chronic hypoxia.
The phenomenon that different gene expression profiles, with different genes, can describe the same process, has been reported before [38]. Given the fact that over 4,000 genes are hypoxia-influenced, it seems reasonable to assume that multiple robust, but different, hypoxia gene expression profiles can be assembled [39]. It should be noted that all gene expression profiles have been applied with the same method in our study. Due to lack of access to the reference cohort, it was not always possible to apply gene expression profiles as in the original publication. To be able to compare the gene expression profiles, we decided to apply the same method to all gene expression profiles.

3.1. Acute and Chronic Hypoxia

The terms acute and chronic hypoxia are simplified terms to describe a complex spectrum of hypoxic micro environmental alterations in a tumor [40]. While an absolute distinction between the two cannot be made, many suggestions for the separate origin, measurement and treatment of the two entities have been published [41,42,43,44]. Janssen et al. employed various staining protocols to study acute and chronic hypoxia in head and neck tumors [42]. They showed that tumors contained on average 15% acute hypoxic (proliferating cells around non-perfused vessels) and around 30% chronic hypoxic areas (cells at a large distance from blood vessels). The two different areas did not overlap. Also in vitro studies showed that cells that had been under hypoxia for a short time, showed a different gene expression profile compared to cells that were hypoxic for a prolonged time [31]. As described by Lendahl et al. 4,047 genes were hypoxia-regulated in a colon carcinoma cell line. Only 52 genes were specific to the acute (1 or 2 h) hypoxia response. 144 genes were up- or downregulated by both acute and chronic (24 h) hypoxia, whereas the majority of the genes (4,005) were chronic hypoxia specific [39]. Nonetheless, in the past decades research has focused on generating a gene expression profile for ‘hypoxia’ in general, without distinction of acute and chronic. The fact that all four clinical HNSCC hypoxia profiles correlated with chronic hypoxia could be due to the large excess of genes regulated by chronic hypoxia [39]. Also, the methods used to generate the gene expression profile could be an explanation that the profiles correlated with chronic hypoxia. For example, the Toustrup et al. profile is based on Eppendorf probe measurements to find relevant genes. If indeed, the area of chronic hypoxia is on average twice the area of acute hypoxia, as reported by Janssen et al., this could lead to a bias towards genes that are correlated with chronic hypoxia. Winter, Buffa and Eustace et al. started with 10 hypoxia ‘seed genes’ to develop their signatures. In our data, these 10 genes were not correlated with in vitro Seigneuric acute hypoxia but showed correlations to the in vitro Seigneuric chronic hypoxia profile (Supplementary Materials, Table S8).

3.2. Acute Hypoxia and Prognosis

The importance of acute hypoxia has been recognized for decades [45]. For example, Chan et al. showed that a human lung squamous cell carcinoma cell line (H1299) becomes more radioresistant under acute hypoxia than under chronic hypoxia, with respective oxygen enhancement ratios of 1.96 and 1.37 [46]. Unfortunately, conclusive data on the separate and combined prognostic effects of acute and chronic hypoxia in HNSCC are lacking. This might be due to the fact that it is difficult to distinguish both types of hypoxia with immunohistochemistry. We found that patients with high acute or chronic hypoxia expression, had a 3.1 and 1.9 times higher risk of a local recurrence in the Pramana cohort, respectively. This was confirmed in the validation set of 174 patients with hazard ratios of 3.3 and 1.4 for local failure, for acute and chronic hypoxia, respectively. Although in both sets chronic hypoxia was not significant, the effect size appears comparable to previously reported hazard ratios for chronic hypoxia [13,25,26,27]. For chronic hypoxia gene expression profiles, the general deduction is that more hypoxic tumors are approximately twice as likely to recur than less hypoxic tumors. This effect could be underestimated due to a division into two hypoxia groups according to the median. We therefore used the scores as continuous variables in the multivariable analysis, which shows that acute hypoxia is a stronger prognostic factor than chronic hypoxia. Literature has suggested that cells lacking functional p53 are more susceptible to genomic instability and potential tumorigenesis when experiencing reoxygenation after acute hypoxia compared to chronic hypoxia [47]. This might also explain why the effect is more pronounced in the Van der Heijden cohort, which contains only HPV-negative HNSCC in which p53 mutations are highly prevalent. Since acute and chronic hypoxia have a different etiology, knowledge about which type of hypoxia causes radioresistance in a specific patient, could lead to the use and development of strategies targeting acute or chronic hypoxia.

4. Materials and Methods

4.1. Gene Expression Profiles

We selected four gene expression profiles for hypoxia that are clinically validated for head and neck cancer, namely, Winter et al., Buffa et al., Toustrup et al. and Eustace et al. Gene expression profiles generated by Seigneuric et al. were used for in vitro acute and chronic hypoxia. Since 2% oxygen is relatively close to oxygen concentration in normal tissue we used the <0.02% oxygen derived profiles as the in vitro acute and chronic hypoxia profiles. There is no clear cut-off for acute and chronic hypoxia. To simplify the complex spectrum of hypoxia, we used the Seigneuric 1 h time point profile for acute hypoxia, while the Seigneuric 24 h time point profile was used for chronic hypoxia [6].

4.2. Datasets

We initially collected gene expression data of three different HNSCC patient cohorts (Pramana, De Jong 1 and De Jong 2), comprising 224 patients in total (Table 1). Extensive patient characteristics for each cohort are available in the Supplementary Materials, Tables S3–S6 and the original publications [32,33,34]. Clinical and expression data of Pramana et al. has been previously published. Briefly, the dataset of Pramana et al. consisted of 91 HPV-positive and negative HNSCC patients all treated with definitive cisplatin based chemoradiotherapy [32]. Gene expression data were obtained from fresh-frozen pre-treatment material, analyzed using dual-channel Operon microarray slides. The full dataset is available at (GEO database) and extensive methods are available in the original publication.
Data of 52 matched patients of the De Jong 1 cohort was previously published [33]. This dataset was expanded with an additional 47 samples that could not be used in the original study. Briefly, gene expression of 99 fresh-frozen laryngeal and oropharyngeal (HPV status unknown) carcinomas was determined using the Illumina beads microarray platform. All patients were treated with radiotherapy. The full dataset is available at (GEO database, We are in the process of submitting to GEO and will provide the numbers a.s.a.p.) and has been generated as described previously.
Clinical data from the De Jong 2 cohort are as described previously [34]. The cohort consists of 34 laryngeal carcinomas treated with radiotherapy. Collection of tumor material and clinical data were approved by the Institutional Review Board of the Netherlands Cancer Institute and all patients signed informed consent. RNA was extracted from pre-treatment FFPE biopsies using the Roche High Pure miRNA Isolation Kit (REF: 05080576001). Quality and quantity of the total RNA was assessed with the 2100 Bioanalyzer using a Nano chip. TruSeq cDNA libraries were generated using the TruSeq RNA Library Preparation Kit v2 sample preparation kit (Illumina Cat.No RS-122-2001/2). The reads (51bp) were sequenced on a HiSeq2000 using V3 chemistry (Illumina Inc., San Diego CA, USA), aligned against the human genome, build 37, using Tophat (version 2.0.6). HTseq count was used to count the number of reads per gene. Only uniquely mapped reads were counted. Detailed methods are provided in the Supplementary Methods. Data are available at http://ega-archive.org.
To validate our findings, we used gene expression data collected and generated within a separate study as described by van der Heijden et al. and deposited at http://ega-archive.org [35,36]. The “van der Heijden cohort” combines two independent cohorts, the NKI-CRAD and the DESIGN cohort, and comprises 174 HPV negative advanced stage HNSCC patients. All patients were treated with definitive cisplatin-based chemoradiotherapy. RNA was isolated from fresh-frozen material and gene-expression was determined using RNA-sequencing as described previously.

4.3. Data Analysis

All analyses were performed in R 3.4.3 using RStudio 1.1. The gene expression profiles consisted of genes exclusively upregulated under hypoxia. All expression data were normalized as described in the original publications. In the Pramana and both De Jong cohorts, per sample hypoxia scores were calculated for each gene expression profile using the mean normalized expression of all genes in the profile. In order to combine the three different patient cohorts, scores were rank-normalized per profile between 0 and 1 for each patient cohort. Gene-expression scores were calculated using Gene Set Variation Analysis (GSVA) [48] in the van der Heijden dataset. The median was used to split patients into “High” and “Low” hypoxia groups. Cox proportional hazard analyses were used to test associations with patient outcome (Overall Survival (OS), Progression Free Survival (PFS), Local Control (LC), Locoregional Control (LRC) and Distant Metastasis (DM)). Multivariable Cox proportional hazard analyses were performed in the Van der Heijden cohort and no cut-offs were used to test the continuous variables.

5. Conclusions

All four clinical HNSCC hypoxia profiles correlate with chronic and not with acute hypoxia expression profiles. However, the acute hypoxia profile has a stronger association with patient outcome after chemoradiotherapy in both the test and validation cohort. Multivariable analysis shows the value of acute hypoxia in addition to well-known prognostic factors for HNSCC. Acute hypoxia gene expression should therefore be incorporated into existing hypoxia-based prediction models.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/11/4/583/s1. Figure S1. Univariable analyses of individual hypoxia profiles. Figure S2. Distribution and outcome analyses of hypoxia in the Van der Heijden cohort. Figure S3. Combined acute and chronic hypoxia scores in the Van der Heijden cohort (n = 174), Table S1. Overview of the four clinical hypoxia gene expression profiles associated with patient outcome in HNSCC, Table S2. Genes in the used gene expression profiles, Table S3. Patient characteristics Pramana cohort: 91 HNSCC stage III-IV chemoradiotherapy patients, Table S4. Patient characteristics de Jong 1 cohort: 99 larynx/oropharynx radiotherapy patients, Table S5. Patient characteristics de Jong 2 cohort: 34 larynx radiotherapy patients, Table S6. Patient characteristics of the van der Heijden cohort (NKI-CRAD and DESIGN combined), Table S7. Univariate Cox proportional hazard analysis of parameters with patient outcome in the Van der Heijden cohort, Table S8. 10 hypoxia ‘seed genes’ and correlation with in vitro Seigneuric acute and Seigneuric chronic hypoxia profiles

Author Contributions

Conceptualization, M.v.d.H., M.C.d.J., M.V. and M.W.M.v.d.B.; Data curation, M.v.d.H., M.C.d.J., C.V.M.V., R.H.d.R. and S.S.; Formal analysis, M.v.d.H. and M.C.d.J.; Funding acquisition, F.H., C.R.L., R.H.B., C.V., M.V. and M.W.M.v.d.B.; Resources, F.H., C.R.L., R.H.B., C.V., M.V. and M.W.M.v.d.B.; Supervision, M.V. and M.W.M.v.d.B.; Visualization, M.v.d.H. and M.C.d.J.; Writing—original draft, M.v.d.H. and M.C.d.J.; Writing—review & editing, C.V.M.V., R.H.d.R., S.S., F.H., C.R.L., R.H.B., C.V., M.V. and M.W.M.v.d.B.

Funding

This study was funded by: EU FP7 programme (Artforce #257144), Brunel, Dutch Cancer Society (Design - KWF-A6C7072) and the Verwelius foundation.

Acknowledgments

We would like to thank the NKI Genomics Core Facility for performing RNA sequencing (Wim Brugman, Ron Kerkhoven and Iris de Rink), the NKI RHPC facility for providing computational resources and the Core Facility-Molecular Pathology and Biobank (CFMPB) for collecting and preparing tissue samples (Annegien Broeks, Sten Cornelissen, Linde Braaf and Dennis Peters). We thank Brad Wouters for his helpful suggestions.

Conflicts of Interest

The authors declare no potential conflict of interest.

References

  1. Ferlay, J.; Shin, H.R.; Bray, F.; Forman, D.; Mathers, C.; Parkin, D.M. Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int. J. Cancer 2010. [Google Scholar] [CrossRef]
  2. Groome, P.A.; Schulze, K.; Boysen, M.; Hall, S.F.; Mackillop, W.J.; O’Sullivan, B.; Irish, J.C.; Warde, P.R.; Schneider, K.M.; Mackenzie, R.G.; et al. A comparison of published head and neck stage groupings in laryngeal cancer using data from two countries. J. Clin. Epidemiol. 2002. [Google Scholar] [CrossRef]
  3. Groome, P.A.; Schulze, K.M.; Mackillop, W.J.; Grice, B.; Goh, C.; Cummings, B.J.; Hall, S.F.; Liu, F.F.; Payne, D.; Rothwell, D.M.; et al. A comparison of published head and neck stage groupings in carcinomas of the tonsillar region. Cancer 2001, 92, 1484–1494. [Google Scholar] [CrossRef] [Green Version]
  4. Hall, S.F.; Groome, P.A.; Irish, J.; O’Sullivan, B. TNM-based stage groupings in head and neck cancer: Application in cancer of the hypopharynx. Head Neck 2009, 31, 1–8. [Google Scholar] [CrossRef] [Green Version]
  5. de Jong, M.C.; Pramana, J.; Knegjens, J.L.; Balm, A.J.M.; van den Brekel, M.W.M.; Hauptmann, M.; Begg, A.C.; Rasch, C.R. HPV and high-risk gene expression profiles predict response to chemoradiotherapy in head and neck cancer, independent of clinical factors. Radiother. Oncol. 2010, 95, 365–370. [Google Scholar] [CrossRef] [PubMed]
  6. Janssen, H.L.; Haustermans, K.M.; Balm, A.J.; Begg, A.C. Hypoxia in head and neck cancer: How much, how important? Head Neck 2005, 27, 622–638. [Google Scholar] [CrossRef] [PubMed]
  7. Bristow, R.G.; Hill, R.P. Hypoxia and metabolism: Hypoxia, DNA repair and genetic instability. Nat. Rev. Cancer 2008, 8, 180–892. [Google Scholar] [CrossRef] [PubMed]
  8. Ma, N.Y.; Tinganelli, W.; Maier, A.; Durante, M.; Kraft-Weyrather, W. Influence of chronic hypoxia and radiation quality on cell survival. J. Radiat. Res. 2013, 54 (Suppl. 1), i13–i22. [Google Scholar] [CrossRef] [Green Version]
  9. Gray, L.H.; Conger, A.D.; Ebert, M.; Hornsey, S.; Scott, O.C. The concentration of oxygen dissolved in tissues at the time of irradiation as a factor in radiotherapy. Br. J. Radiol. 1953, 26, 638–648. [Google Scholar] [CrossRef] [PubMed]
  10. De Gonzalez, A.B.; Curtis, R.E.; Kry, S.F.; Gilbert, E.; Lamart, S.; Berg, C.D.; Stovall, M.; Ron, E. Proportion of second cancers attributable to radiotherapy treatment in adults: A cohort study in the US SEER cancer registries. Lancet Oncol. 2011, 12, 353–360. [Google Scholar] [CrossRef]
  11. Overgaard, J. Hypoxic modification of radiotherapy in squamous cell carcinoma of the head and neck—A systematic review and meta-analysis. Radiother. Oncol. 2011, 100, 22–32. [Google Scholar] [CrossRef] [PubMed]
  12. Winter, S.C.; Shah, K.A.; Han, C.; Campo, L.; Turley, H.; Leek, R.; Corbridge, R.J.; Cox, G.J.; Harris, A.L. The Relation Between Hypoxia-Inducible Factor (HIF) -1 a and HIF-2 a Expression With Anemia and Outcome in Surgically Treated Head and Neck Cancer. Cancer. 2006. [Google Scholar] [CrossRef]
  13. Toustrup, K.; Sørensen, B.S.; Nordsmark, M.; Busk, M.; Wiuf, C.; Alsner, J.; Overgaard, J. 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–5931. [Google Scholar] [CrossRef]
  14. Rademakers, S.E.; Hoogsteen, I.J.; Rijken, P.F.; Oosterwijk, E.; Terhaard, C.H.; Doornaert, P.A.; Langendijk, J.A.; van den Ende, P.; Takes, R.; De Bree, R.; et al. Pattern of CAIX expression is prognostic for outcome and predicts response to ARCON in patients with laryngeal cancer treated in a phase III randomized trial. Radiother. Oncol. 2013, 108, 517–522. [Google Scholar] [CrossRef]
  15. Thomson, D.; Yang, H.; Baines, H.; Miles, E.; Bolton, S.; West, C.; Slevin, N. NIMRAD—A phase III trial to investigate the use of nimorazole hypoxia modification with intensity-modulated radiotherapy in head and neck cancer. Clin. Oncol. 2014, 26, 344–347. [Google Scholar] [CrossRef] [PubMed]
  16. Nordsmark, M.; Bentzen, S.M.; Rudat, V.; Brizel, D.; Lartigau, E.; Stadler, P.; Becker, A.; Adam, M.; Molls, M.; Dunst, J.; 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. [Google Scholar] [CrossRef]
  17. Brizel, D.M.; Sibley, G.S.; Prosnitz, L.R.; Scher, R.L.; Dewhirst, M.W. Tumor hypoxia adversely affects the prognosis of carcinoma of the head and neck. Int. J. Radiat. Oncol. Biol. Phys. 1997, 38, 285–289. [Google Scholar] [CrossRef]
  18. Nordsmark, M.; Overgaard, J. A confirmatory prognostic study on oxygenation status and loco-regional control in advanced head and neck squamous cell carcinoma treated by radiation therapy. Radiother. Oncol. 2000, 57, 39–43. [Google Scholar] [CrossRef]
  19. Nordsmark, M.; Overgaard, M.; Overgaard, J. Pretreatment oxygenation predicts radiation response in advanced squamous cell carcinoma of the head and neck. Radiother. Oncol. 1996, 41, 31–39. [Google Scholar] [CrossRef]
  20. Kaanders, J.H.; Wijffels, K.I.; Marres, H.A.; Ljungkvist, A.S.; Pop, L.A.; van den Hoogen, F.J.; de Wilde, P.C.; Bussink, J.; Raleigh, J.A.; van der Kogel, A.J. Pimonidazole binding and tumor vascularity predict for treatment outcome in head and neck cancer. Cancer Res. 2002, 62, 7066–7074. [Google Scholar] [PubMed]
  21. Schrijvers, M.L.; van der Laan, B.F.A.M.; de Bock, G.H.; Pattje, W.J.; Mastik, M.F.; Menkema, L.; Langendijk, J.A.; Kluin, P.M.; Schuuring, E.; van der Wal, J.E. Overexpression of Intrinsic Hypoxia Markers HIF1α and CA-IX Predict for Local Recurrence in Stage T1–T2 Glottic Laryngeal Carcinoma Treated With Radiotherapy. Int. J. Radiat. Oncol. Biol. Phys. 2008, 72, 161–169. [Google Scholar] [CrossRef]
  22. Koukourakis, M.I.; Giatromanolaki, A.; Sivridis, E.; Simopoulos, C.; Turley, H.; Talks, K.; Gatter, K.C.; Harris, A.L. Hypoxia-inducible factor (HIF1A and HIF2A), angiogenesis, and chemoradiotherapy outcome of squamous cell head-and-neck cancer. Int. J. Radiat. Oncol. Biol. Phys. 2002, 53, 1192–1202. [Google Scholar] [CrossRef]
  23. Panek, R.; Welsh, L.; Dunlop, A.; Wong, K.H.; Riddell, A.M.; Koh, D.-M.; Schmidt, M.A.; Doran, S.; Mcquaid, D.; Hopkinson, G.; et al. Repeatability and sensitivity of T2* measurements in patients with head and neck squamous cell carcinoma at 3T. J. Magn. Reson. Imaging 2016, 44, 72–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Fleming, I.N.; Manavaki, R.; Blower, P.J.; West, C.; Williams, K.J.; Harris, A.L.; Domarkas, J.; Lord, S.; Baldry, C.; Gilbert, F.J. Imaging tumour hypoxia with positron emission tomography. Br. J. Cancer 2015, 112, 238–250. [Google Scholar] [CrossRef] [PubMed]
  25. Winter, S.C.; Buffa, F.M.; Silva, P.; Miller, C.; Valentine, H.R.; Turley, H.; Shah, K.A.; Cox, G.J.; Corbridge, R.J.; Homer, J.J.; et al. Relation of a hypoxia metagene derived from head and neck cancer to prognosis of multiple cancers. Cancer Res. 2007, 67, 3441–3449. [Google Scholar] [CrossRef] [PubMed]
  26. Buffa, F.M.; Harris, A.L.; West, C.M.; Miller, C.J. Large meta-analysis of multiple cancers reveals a common, compact and highly prognostic hypoxia metagene. Br. J. Cancer 2010, 102, 428–435. [Google Scholar] [CrossRef] [Green Version]
  27. Eustace, A.; Mani, N.; Span, P.N.; Irlam, J.J.; Taylor, J.; Betts, G.N.; Denley, H.; Miller, C.J.; Homer, J.J.; Rojas, A.M.; 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–4888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Koukourakis, M.I.; Bentzen, S.M.; Giatromanolaki, A.; Wilson, G.D.; Daley, F.M.; Saunders, M.I.; Dische, S.; Sivridis, E.; Harris, A.L. Endogenous markers of two separate hypoxia response pathways (hypoxia inducible factor 2 alpha and carbonic anhydrase 9) are associated with radiotherapy failure in head and neck cancer patients recruited in the CHART randomized trial. J. Clin. Oncol. 2006, 24, 727–735. [Google Scholar] [CrossRef]
  29. Rademakers, S.E.; Lok, J.; van der Kogel, A.J.; Bussink, J.; Kaanders, J.H.A.M. Metabolic markers in relation to hypoxia; staining patterns and colocalization of pimonidazole, HIF-1α, CAIX, LDH-5, GLUT-1, MCT1 and MCT4. BMC Cancer 2011, 11, 167. [Google Scholar] [CrossRef]
  30. Tawk, B.; Schwager, C.; Deffaa, O.; Dyckhoff, G.; Warta, R.; Linge, A.; Krause, M.; Weichert, W.; Baumann, M.; Herold-Mende, C.; et al. Comparative analysis of transcriptomics based hypoxia signatures in head- and neck squamous cell carcinoma. Radiother. Oncol. 2016, 118, 350–358. [Google Scholar] [CrossRef]
  31. Seigneuric, R.; Starmans, M.H.W.; Fung, G.; Krishnapuram, B.; Nuyten, D.S.A.; van Erk, A.; Magagnin, M.G.; Rouschop, K.M.; Krishnan, S.; Rao, R.B.; et al. Impact of supervised gene signatures of early hypoxia on patient survival. Radiother. Oncol. 2007, 83, 374–382. [Google Scholar] [CrossRef]
  32. Pramana, J.; Van den Brekel, M.W.M.; van Velthuysen, M.-L.F.; Wessels, L.F.A.; Nuyten, D.S.; Hofland, I.; Atsma, D.; Pimentel, N.; Hoebers, F.J.; Rasch, C.R.; et al. Gene expression profiling to predict outcome after chemoradiation in head and neck cancer. Int. J. Radiat. Oncol. Biol. Phys. 2007, 69, 1544–1552. [Google Scholar] [CrossRef] [PubMed]
  33. De Jong, M.C.; Pramana, J.; Van Der Wal, J.E.; Lacko, M.; Peutz-Kootstra, C.J.; De Jong, J.M.; Takes, R.P.; Kaanders, J.H.; van der Laan, B.F.; Wachters, J.; et al. CD44 expression predicts local recurrence after radiotherapy in larynx cancer. Clin. Cancer Res. 2010, 16, 5329–5338. [Google Scholar] [CrossRef]
  34. De Jong, M.C.; Ten Hoeve, J.J.; Grénman, R.; Wessels, L.F.; Kerkhoven, R.; Te Riele, H.; van den Brekel, M.W.; Verheij, M.; Begg, A.C. Pretreatment microRNA expression impacting on epithelial-to-mesenchymal transition predicts intrinsic radiosensitivity in head and neck cancer cell lines and patients. Clin. Cancer Res. 2015, 21, 5630–5638. [Google Scholar] [CrossRef] [PubMed]
  35. Van der Heijden, M.; Essers, P.; Verheij, M.; Van den Brekel, M.; Vens, C. OC-0487: EMT signatures as a prognostic marker for metastasis in HPV-negative HNSCC. Radiother. Oncol. 2018, 127, S250–S251. [Google Scholar] [CrossRef]
  36. Van der Heijden, M.; Essers, P.; Verhagen, C.; Willems, S.; Sanders, J.; de Roest, R.; Vossen, D.M.; Leemans, C.R.; Verheij, M.; Brakenhoff, R.H.; et al. Epithelial-to-mesenchymal transistion is a prognostic marker for patient outcome in advanced stage HNSCC patients treated with chemoradiotherapy. Manuscript under review.
  37. Linge, A.; Lock, S.; Gudziol, V.; Nowak, A.; Lohaus, F.; Von Neubeck, C.; Jütz, M.; Abdollahi, A.; Debus, J.; Tinhofer, I.; 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–2649. [Google Scholar] [CrossRef] [PubMed]
  38. Roepman, P.; Kemmeren, P.; Wessels, L.F.A.; Slootweg, P.J.; Holstege, F.C.P. Multiple robust signatures for detecting lymph node metastasis in head and neck cancer. Cancer Res. 2006, 66, 2361–2366. [Google Scholar] [CrossRef] [PubMed]
  39. Lendahl, U.; Lee, K.L.; Yang, H.; Poellinger, L. Generating specificity and diversity in the transcriptional response to hypoxia. Nat. Rev. Genet. 2009, 10, 821–832. [Google Scholar] [CrossRef] [PubMed]
  40. Bayer, C.; Shi, K.; Astner, S.T.; Maftei, C.A.; Vaupel, P. Acute versus chronic hypoxia: Why a simplified classification is simply not enough. Int. J. Radiat. Oncol. Biol. Phys. 2011, 80, 965–968. [Google Scholar] [CrossRef]
  41. Bayer, C.; Vaupel, P. Acute versus chronic hypoxia in tumors: Controversial data concerning time frames and biological consequences. Strahlenther. Onkol. 2012, 188, 616–627. [Google Scholar] [CrossRef]
  42. Janssen, H.L.K.; Haustermans, K.M.G.; Sprong, D.; Blommestijn, G.; Hofland, I.; Hoebers, F.J.; Blijweert, E.; Raleigh, J.A.; Semenza, G.L.; Varia, M.A.; et al. HIF-1α, pimonidazole, and iododeoxyuridine to estimate hypoxia and perfusion in human head-and-neck tumors. Int. J. Radiat. Oncol. Biol. Phys. 2002, 54, 1537–1549. [Google Scholar] [CrossRef]
  43. Maftei, C.A.; Bayer, C.; Shi, K.; Vaupel, P. Intra-and intertumor heterogeneities in total, chronic, and acute hypoxia in xenografted squamous cell carcinomas: Detection and quantification using (immuno-)fluorescence techniques. Strahlenther. Onkol. 2012, 188, 606–615. [Google Scholar] [CrossRef]
  44. Wijffels, K.I.; Kaanders, J.H.; Rijken, P.F.; Bussink, J.; van den Hoogen, F.J.; Marres, H.A.; de Wilde, P.C.; Raleigh, J.A.; van der Kogel, A.J. Vascular architecture and hypoxic profiles in human head and neck squamous cell carcinomas. Br. J. Cancer 2000, 83, 674–683. [Google Scholar] [CrossRef] [Green Version]
  45. Brown, J.M. Evidence for acutely hypoxic cells in mouse tumours, and a possible mechanism of reoxygenation. Br. J. Radiol. 1979. [Google Scholar] [CrossRef] [PubMed]
  46. Chan, N.; Koritzinsky, M.; Zhao, H.; Bindra, R.; Glazer, P.M.; Powell, S.; Belmaaza, A.; Wouters, B.; Bristow, R.G. Chronic hypoxia decreases synthesis of homologous recombination proteins to offset chemoresistance and radioresistance. Cancer Res. 2008, 68, 605–614. [Google Scholar] [CrossRef]
  47. Pires, I.M.; Bencokova, Z.; Milani, M.; Folkes, L.K.; Li, J.A.; Stratford, M.R.; Harris, A.L.; Hammond, E.M. Effects of acute versus chronic hypoxia on DNA damage responses and genomic instability. Cancer Res. 2010, 70, 925–935. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Hänzelmann, S.; Castelo, R.; Guinney, J. GSVA: Gene set variation analysis for microarray and RNA-Seq data. BMC Bioinform. 2013, 14. [Google Scholar] [CrossRef]
Cancers 11 00583 g001 550
Figure 1. Venn diagram showing the number of overlapping genes in the four clinical hypoxia profiles and the Seigneuric acute and chronic in vitro profiles.
Figure 1. Venn diagram showing the number of overlapping genes in the four clinical hypoxia profiles and the Seigneuric acute and chronic in vitro profiles.
Cancers 11 00583 g001
Cancers 11 00583 g002 550
Figure 2. Correlations and clustering of the hypoxia scores from the clinical hypoxia profiles and the Seigneuric acute and chronic in vitro hypoxia profiles in 224 HNSCC patients. (A) Spearman correlations (upper right panels) and scatter plots (lower left panels) of all possible pairs of hypoxia profiles for 224 patients. All Spearman correlations were significant at the p < 0.0001 level. (B) Heatmap showing the scores for the four hypoxia profiles and in vitro profiles as indicated in 224 patients.
Figure 2. Correlations and clustering of the hypoxia scores from the clinical hypoxia profiles and the Seigneuric acute and chronic in vitro hypoxia profiles in 224 HNSCC patients. (A) Spearman correlations (upper right panels) and scatter plots (lower left panels) of all possible pairs of hypoxia profiles for 224 patients. All Spearman correlations were significant at the p < 0.0001 level. (B) Heatmap showing the scores for the four hypoxia profiles and in vitro profiles as indicated in 224 patients.
Cancers 11 00583 g002
Cancers 11 00583 g003 550
Figure 3. Combined acute and chronic hypoxia scores. (A) Kaplan–Meier curve showing local control of the “High” versus “Low” groups of the acute hypoxia profile. (B) Kaplan–Meier curve showing local control of the “High” versus “Low” groups of the chronic hypoxia profile. (C) Kaplan–Meier curve showing local control for 3 groups: (1) acute and chronic hypoxia both low, (2) acute or chronic hypoxia high, or (3) acute and chronic hypoxia both high. The crosstab shows the percent of local recurrences per subgroup in 91 chemoradiotherapy patients. Samples were divided into two groups using the median. Cells are colored in a color corresponding with the line color in the Kaplan–Meier curve. P-value represent the log-rank p-value. (D) Crosstab showing the percentage local recurrences per subgroup.
Figure 3. Combined acute and chronic hypoxia scores. (A) Kaplan–Meier curve showing local control of the “High” versus “Low” groups of the acute hypoxia profile. (B) Kaplan–Meier curve showing local control of the “High” versus “Low” groups of the chronic hypoxia profile. (C) Kaplan–Meier curve showing local control for 3 groups: (1) acute and chronic hypoxia both low, (2) acute or chronic hypoxia high, or (3) acute and chronic hypoxia both high. The crosstab shows the percent of local recurrences per subgroup in 91 chemoradiotherapy patients. Samples were divided into two groups using the median. Cells are colored in a color corresponding with the line color in the Kaplan–Meier curve. P-value represent the log-rank p-value. (D) Crosstab showing the percentage local recurrences per subgroup.
Cancers 11 00583 g003
Cancers 11 00583 g004 550
Figure 4. Validation of prognostic value of combined acute and chronic hypoxia scores in the van der Heijden cohort. (A) Kaplan–Meier curve showing local control of the “High” versus “Low” acute hypoxia. (B) Kaplan–Meier curve showing local control of the “High” versus “Low” chronic hypoxia. (C) Kaplan–Meier curve overall survival (C) and local control (D) for 3 groups: acute and chronic hypoxia both low, acute or chronic hypoxia high or acute and chronic hypoxia both high. (E) Crosstab showing the percentage of local recurrences per subgroup in the validation cohort consisting of 174 patients treated with chemoradiotherapy. Cells are colored in a color corresponding with the line color in the Kaplan–Meier curves in panel C and D.
Figure 4. Validation of prognostic value of combined acute and chronic hypoxia scores in the van der Heijden cohort. (A) Kaplan–Meier curve showing local control of the “High” versus “Low” acute hypoxia. (B) Kaplan–Meier curve showing local control of the “High” versus “Low” chronic hypoxia. (C) Kaplan–Meier curve overall survival (C) and local control (D) for 3 groups: acute and chronic hypoxia both low, acute or chronic hypoxia high or acute and chronic hypoxia both high. (E) Crosstab showing the percentage of local recurrences per subgroup in the validation cohort consisting of 174 patients treated with chemoradiotherapy. Cells are colored in a color corresponding with the line color in the Kaplan–Meier curves in panel C and D.
Cancers 11 00583 g004
Table
Table 1. Summary of characteristics of the five HNSCC patient cohorts.
Table 1. Summary of characteristics of the five HNSCC patient cohorts.
UseStudyCohortsNSitesTreatmentMaterialAssay
Combined for classification analysesPramana [32]Stage III-IV HNSCCs91All head and neckChemo-
radiotherapy		FFDual channel Operon microarray
test cohort
De Jong 1 [33]Larynx / oropharynx99Larynx/oropharynxRadiotherapyFFIllumina beads microarray
De Jong 2 [34]T2-3 larynx 34LarynxRadiotherapyFFPERNAseq
Combined for validation cohortVan der Heijden - NKI-CRAD [35,36]Stage III-IV HNSCCs98Larynx, hypopharynx, HPV-neg oropharynxChemo-
radiotherapy		FFRNAseq
Van der Heijden – DESIGN [35,36]Stage III-IV HNSCCs76Larynx, hypopharynx, HPV-neg oropharynxChemo-
radiotherapy		FFRNAseq
FF: Fresh-Frozen; FFPE: Formalin Fixed Paraffin Embedded.
Table
Table 2. Multivariable Cox proportional hazard analysis of parameters with patient outcome in the Van der Heijden cohort.
Table 2. Multivariable Cox proportional hazard analysis of parameters with patient outcome in the Van der Heijden cohort.
Variable OSPFSLCLRCDM
HRp-ValueHRp-ValueHRp-ValueHRp-ValueHRp-Value
Gender
Female0.830.040.560.040.160.0130.190.0040.280.04
MaleREF REF REF REF REF
Tumor site
OropharynxREF REF REF REF REF
Hypopharynx0.640.090.690.150.960.120.970.940.520.16
Larynx0.580.90.790.451.270.561.080.861.080.87
Disease stage
II-III0.830.570.780.453.740.0371.680.30.380.2
IVA-IVBREF REF REF REF REF
Cumulative cisplatin < 200
NoREF REF REF REF REF
Yes2.030.0031.930.0046.40.00032.70.0120.810.64
Tumor Volume1.010.0000441.0080.00061.0140.0341.0070.161.0060.19
Acute Hypoxia4.140.0214.020.01746.40.014300.00225.50.13
Chronic Hypoxia0.830.6510.990.510.481.230.751.260.73
Acute Hypoxia * Tumor Volume0.980.0070.980.030.960.120.960.0170.980.26
Chronic Hypoxia * Tumor Volume10.990.990.710.980.990.520.990.57
OS: Overall Survival; PFS: Progression Free Survival; LC: Local Control; LRC: Locoregional Control; DM: Distant Metastasis Free Survival; HR: Hazard Ratio. *: Interaction between the two variables.

Share and Cite

MDPI and ACS Style

van der Heijden, M.; de Jong, M.C.; Verhagen, C.V.M.; de Roest, R.H.; Sanduleanu, S.; Hoebers, F.; Leemans, C.R.; Brakenhoff, R.H.; Vens, C.; Verheij, M.; et al. Acute Hypoxia Profile is a Stronger Prognostic Factor than Chronic Hypoxia in Advanced Stage Head and Neck Cancer Patients. Cancers 2019, 11, 583. https://doi.org/10.3390/cancers11040583

AMA Style

van der Heijden M, de Jong MC, Verhagen CVM, de Roest RH, Sanduleanu S, Hoebers F, Leemans CR, Brakenhoff RH, Vens C, Verheij M, et al. Acute Hypoxia Profile is a Stronger Prognostic Factor than Chronic Hypoxia in Advanced Stage Head and Neck Cancer Patients. Cancers. 2019; 11(4):583. https://doi.org/10.3390/cancers11040583

Chicago/Turabian Style

van der Heijden, Martijn, Monique C. de Jong, Caroline V. M. Verhagen, Reinout H. de Roest, Sebastian Sanduleanu, Frank Hoebers, C. René Leemans, Ruud H. Brakenhoff, Conchita Vens, Marcel Verheij, and et al. 2019. "Acute Hypoxia Profile is a Stronger Prognostic Factor than Chronic Hypoxia in Advanced Stage Head and Neck Cancer Patients" Cancers 11, no. 4: 583. https://doi.org/10.3390/cancers11040583

APA Style

van der Heijden, M., de Jong, M. C., Verhagen, C. V. M., de Roest, R. H., Sanduleanu, S., Hoebers, F., Leemans, C. R., Brakenhoff, R. H., Vens, C., Verheij, M., & van den Brekel, M. W. M. (2019). Acute Hypoxia Profile is a Stronger Prognostic Factor than Chronic Hypoxia in Advanced Stage Head and Neck Cancer Patients. Cancers, 11(4), 583. https://doi.org/10.3390/cancers11040583

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop