Epithelial-to-mesenchymal transition is a prognostic marker for patient outcome in advanced stage HNSCC patients treated with chemoradiotherapy.

BACKGROUND: The prognosis of patients with HPV-negative advanced stage head and neck squamous cell carcinoma (HNSCC) remains poor. No prognostic markers other than TNM staging are routinely used in clinic. Epithelial-to-mesenchymal transition (EMT) has been shown to be a strong prognostic factor in other cancer types. The purpose of this study was to determine the role of EMT in HPV-negative HNSCC outcomes. METHODS: Pretreatment tumor material from patients of two cohorts, totalling 174 cisplatin-based chemoradiotherapy treated HPV-negative HNSCC patients, was RNA-sequenced. Seven different EMT gene expression signatures were used for EMT status classification and generation of HNSCC-specific EMT models using Random Forest machine learning. RESULTS: Mesenchymal classification by all EMT signatures consistently enriched for poor prognosis patients in both cohorts of 98 and 76 patients. Uni- and multivariate analyses show important HR of 1.6-5.8, thereby revealing EMT's role in HNSCC outcome. Discordant classification by these signatures prompted the generation of an HNSCC-specific EMT profile based on the concordantly classified samples in the first cohort (cross-validation AUC > 0.98). The independent validation cohort confirmed the association of mesenchymal classification by the HNSCC-EMT model with poor overall survival (HR = 3.39, p < 0.005) and progression free survival (HR = 3.01, p < 0.005) in multivariate analysis with TNM. Analysis of an additional HNSCC cohort from PET-positive patients with metastatic disease prior to treatment further supports this relationship and reveals a strong link of EMT to the propensity to metastasize. CONCLUSIONS: EMT in HPV-negative HNSCC co-defines patient outcome after chemoradiotherapy. The generated HNSCC-EMT prediction models can function as strong prognostic biomarkers.

Genetic Factors Associated with a Poor Outcome in Head and Neck Cancer Patients Receiving Definitive Chemoradiotherapy.

About half of advanced stage head and neck squamous cell carcinoma (HNSCC) patients can be cured by chemoradiotherapy. Patient outcome may be partially determined by the genetic alterations in HNSCC, rendering these alterations promising candidate prognostic factors and/or therapeutic targets. However, their relevance in patient outcome prognosis remains to be assessed in patients that receive standard-of-care chemoradiotherapy. We therefore tested whether frequent genetic alterations were associated with progression free survival (PFS) in advanced stage HNSCC patients who were uniformly treated with definitive platinum-based chemoradiotherapy. To this end, we performed targeted DNA sequencing on frozen pre-treatment tumor biopsy material from 77 patients with advanced stage oro- and hypopharyngeal carcinoma. This provided somatic point mutation and copy number aberration data of 556 genes. The most frequently mutated genes, TP53 (62%), CCND1 (51%), CDKN2A (30%) and PIK3CA (21%), were not associated with PFS. However, co-occurring CCND1 and CDKN2A mutations were associated with short PFS (HR 2.24, p = 0.028) in HPV-negative tumors. Furthermore, tumor mutational burden (sum of somatic point mutations) showed a trend towards decreased PFS (HR 1.9, p = 0.089), and chromosomal instability (CIN) was associated with shorter PFS (HR 2.3, p = 0.023), independent of HPV status. Our results show that tumor mutational burden, CIN markers, and co-occurring CCND1 and CDKN2A mutations are associated with chemoradiotherapy outcomes in advanced stage oro- and hypopharyngeal HNSCC patients, thereby highlighting their prognostic potential. Given their poor prognosis association and link to biological targets, they may also identify patients for novel targeted therapies and immunotherapies.

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

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.

Drug Sensitivity Prediction Models Reveal a Link between DNA Repair Defects and Poor Prognosis in HNSCC.

Head and neck squamous cell carcinoma (HNSCC) is characterized by the frequent manifestation of DNA crosslink repair defects. We established novel expression-based DNA repair defect markers to determine the clinical impact of such repair defects. Using hypersensitivity to the DNA crosslinking agents, mitomycin C and olaparib, as proxies for functional DNA repair defects in a panel of 25 HNSCC cell lines, we applied machine learning to define gene expression models that predict repair defects. The expression profiles established predicted hypersensitivity to DNA-damaging agents and were associated with mutations in crosslink repair genes, as well as downregulation of DNA damage response and repair genes, in two independent datasets. The prognostic value of the repair defect prediction profiles was assessed in two retrospective cohorts with a total of 180 patients with advanced HPV-negative HNSCC, who were treated with cisplatin-based chemoradiotherapy. DNA repair defects, as predicted by the profiles, were associated with poor outcome in both patient cohorts. The poor prognosis association was particularly strong in normoxic tumor samples and was linked to an increased risk of distant metastasis. In vitro, only crosslink repair-defective HNSCC cell lines are highly migratory and invasive. This phenotype could also be induced in cells by inhibiting rad51 in repair competent and reduced by DNA-PK inhibition. In conclusion, DNA crosslink repair prediction expression profiles reveal a poor prognosis association in HNSCC. SIGNIFICANCE: This study uses innovative machine learning-based approaches to derive models that predict the effect of DNA repair defects on treatment outcome in HNSCC.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/79/21/5597/F1.large.jpg.

Biological Determinants of Chemo-Radiotherapy Response in HPV-Negative Head and Neck Cancer: A Multicentric External Validation.

Purpose: Tumor markers that are related to hypoxia, proliferation, DNA damage repair and stem cell-ness, have a prognostic value in advanced stage HNSCC patients when assessed individually. Here we aimed to evaluate and validate this in a multifactorial context and assess interrelation and the combined role of these biological factors in determining chemo-radiotherapy response in HPV-negative advanced HNSCC. Methods: RNA sequencing data of pre-treatment biopsy material from 197 HPV-negative advanced stage HNSCC patients treated with definitive chemoradiotherapy was analyzed. Biological parameter scores were assigned to patient samples using previously generated and described gene expression signatures. Locoregional control rates were used to assess the role of these biological parameters in radiation response and compared to distant metastasis data. Biological factors were ranked according to their clinical impact using bootstrapping methods and multivariate Cox regression analyses that included clinical variables. Multivariate Cox regression analyses comprising all biological variables were used to define their relative role among all factors when combined. Results: Only few biomarker scores correlate with each other, underscoring their independence. The different biological factors do not correlate or cluster, except for the two stem cell markers CD44 and SLC3A2 (r = 0.4, p < 0.001) and acute hypoxia prediction scores which correlated with T-cell infiltration score, CD8+ T cell abundance and proliferation scores (r = 0.52, 0.56, and 0.6, respectively with p < 0.001). Locoregional control association analyses revealed that chronic (Hazard Ratio (HR) = 3.9) and acute hypoxia (HR = 1.9), followed by stem cell-ness (CD44/SLC3A2; HR = 2.2/2.3), were the strongest and most robust determinants of radiation response. Furthermore, multivariable analysis, considering other biological and clinical factors, reveal a significant role for EGFR expression (HR = 2.9, p < 0.05) and T-cell infiltration (CD8+T-cells: HR = 2.2, p < 0.05; CD8+T-cells/Treg: HR = 2.6, p < 0.01) signatures in locoregional control of chemoradiotherapy-treated HNSCC. Conclusion: Tumor acute and chronic hypoxia, stem cell-ness, and CD8+ T-cell parameters are relevant and largely independent biological factors that together contribute to locoregional control. The combined analyses illustrate the additive value of multifactorial analyses and support a role for EGFR expression analysis and immune cell markers in addition to previously validated biomarkers. This external validation underscores the relevance of biological factors in determining chemoradiotherapy outcome in HNSCC.

Role of variant allele fraction and rare SNP filtering to improve cellular DNA repair endpoint association.

BACKGROUND: Large cancer genome studies continue to reveal new players in treatment response and tumorigenesis. The discrimination of functional alterations from the abundance of passenger genetic alterations still poses challenges and determines DNA sequence variant selection procedures. Here we evaluate variant selection strategies that select homozygous variants and rare SNPs and assess its value in detecting tumor cells with DNA repair defects. METHODS: To this end we employed a panel of 29 patient-derived head and neck squamous cell carcinoma (HNSCC) cell lines, of which a subset harbors DNA repair defects. Mitomycin C (MMC) sensitivity was used as functional endpoint of DNA crosslink repair deficiency. 556 genes including the Fanconi anemia (FA) and homologous recombination (HR) genes, whose products strongly determine MMC response, were capture-sequenced. RESULTS: We show a strong association between MMC sensitivity, thus loss of DNA repair function, and the presence of homozygous and rare SNPs in the relevant FA/HR genes. Excluding such selection criteria impedes the discrimination of crosslink repair status by mutation analysis. Applied to all KEGG pathways, we find that the association with MMC sensitivity is strongest in the KEGG FA pathway, therefore also demonstrating the value of such selection strategies for exploratory analyses. Variant analyses in 56 clinical samples demonstrate that homozygous variants occur more frequently in tumor suppressor genes than oncogenes further supporting the role of a homozygosity criterion to improve gene function association or tumor suppressor gene identification studies. CONCLUSION: Together our data show that the detection of relevant genes or of repair pathway defected tumor cells can be improved by the consideration of allele zygosity and SNP allele frequencies.

Comparative genomic analysis of oral versus laryngeal and pharyngeal cancer.

OBJECTIVE: Locally advanced oral squamous cell carcinoma (OSCC) shows lower locoregional control and disease specific survival rates than laryngeal and pharyngeal squamous cell carcinoma (L/P-SCC) after definitive chemoradiotherapy treatment. Despite clinical factors, this can point towards a different tumor biology that could impact chemoradiotherapy response rates. This prompted us to compare the mutational profiles of OSCC with L/P-SCC. METHODS: We performed target capture DNA sequencing on 111 HPV-negative HNSCC samples (NKI dataset), 55 oral and 56 laryngeal/pharyngeal, and identified somatic point mutations and copy number aberrations. We next expanded our analysis with 276 OSCC and 134 L/P-SCC sample data from The Cancer Genome Atlas (TCGA dataset). We focused our analyses on genes that are frequently mutated in HNSCC. RESULTS: The mutational profiles of OSCC and L/P-SCC showed many similarities. However, OSCC was significantly enriched for CASP8 (NKI: 15% vs 0%; TCGA: 17% vs 2%) and HRAS (TCGA: 10% vs 1%) mutations. LAMA2 (TCGA: 5% vs 19%) and NSD1 (TCGA: 7% vs 25%) mutations were enriched in L/P-SCC. Overall, we find that OSCC had fewer somatic point mutations and copy number aberrations than L/P-SCC. Interestingly, L/P-SCC scored higher in mutational and genomic scar signatures associated with homologous recombination DNA repair defects. CONCLUSION: Despite showing a similar mutational profile, our comparative genomic analysis revealed distinctive features in OSCC and L/P-SCC. Some of these genes and cellular processes are likely to affect the cellular response to radiation or cisplatin. Genomic characterizations may guide or enable personalized treatment in the future.

Fanconi anemia and homologous recombination gene variants are associated with functional DNA repair defects in vitro and poor outcome in patients with advanced head and neck squamous cell carcinoma.

Mutations in Fanconi Anemia or Homologous Recombination (FA/HR) genes can cause DNA repair defects and could therefore impact cancer treatment response and patient outcome. Their functional impact and clinical relevance in head and neck squamous cell carcinoma (HNSCC) is unknown. We therefore questioned whether functional FA/HR defects occurred in HNSCC and whether they are associated with FA/HR variants. We assayed a panel of 29 patient-derived HNSCC cell lines and found that a considerable fraction is hypersensitive to the crosslinker Mitomycin C and PARP inhibitors, a functional measure of FA/HR defects. DNA sequencing showed that these hypersensitivities are associated with the presence of bi-allelic rare germline and somatic FA/HR gene variants. We next questioned whether such variants are associated with prognosis and treatment response in HNSCC patients. DNA sequencing of 77 advanced stage HNSCC tumors revealed a 19% incidence of such variants. Importantly, these variants were associated with a poor prognosis (p = 0.027; HR = 2.6, 1.1-6.0) but favorable response to high cumulative cisplatin dose. We show how an integrated in vitro functional repair and genomic analysis can improve the prognostic value of genetic biomarkers. We conclude that repair defects are marked and frequent in HNSCC and are associated with clinical outcome.

Acute isodense intracerebral haematoma due to coagulopathy associated with prostate cancer.

Prostate cancer can be complicated by spontaneous intracranial haematoma. In the acute setting, an intracranial haematoma can appear isodense on CT scanning, especially if a coagulopathy is present. Diffuse intravascular coagulation (DIC) can be the underlying pathophysiological mechanism that leads to spontaneous bleeding. It is important to recognize the possibility of DIC in cancer patients, since it has important therapeutical implications that may be life saving.

Extent of radiosensitization by the PARP inhibitor olaparib depends on its dose, the radiation dose and the integrity of the homologous recombination pathway of tumor cells.

BACKGROUND AND PURPOSE: The PARP inhibitor olaparib is currently tested in clinical phase 1 trials to define safe dose levels in combination with RT. However, certain clinically relevant insights are still lacking. Here we test, while comparing to single agent activity, the olaparib dose and genetic background dependence of olaparib-mediated radiosensitization. MATERIALS AND METHODS: Long-term growth inhibition and clonogenic assays were used to assess radiosensitization in BRCA2-deficient and BRCA2-complemented cells and in a panel of human head and neck squamous cell carcinoma cell lines. RESULTS: The extent of radiosensitization greatly depended on the olaparib dose, the radiation dose and the homologous recombination status of cells. Olaparib concentrations that resulted in radiosensitization prevented PAR induction by irradiation. Seven hours olaparib exposures were sufficient for radiosensitization. Importantly, the radiosensitizing effects can be observed at much lower olaparib doses than the single agent effects. CONCLUSION: Extrapolation of these data to the clinic suggests that low olaparib doses are sufficient to cause radiosensitization, underlining the potential of the treatment. Here we show that drug doses achieving radiosensitization can greatly differ from those achieving single agent activities, an important consideration when developing combined radiotherapy strategies with novel targeted agents.

Incidence Changes of Human Papillomavirus in Oropharyngeal Squamous Cell Carcinoma and Effects on Survival in the Netherlands Cancer Institute, 1980-2009.

AIM: Human papillomavirus (HPV) is a risk factor for oropharyngeal squamous cell carcinoma (OPSCC), with an increasing incidence. The present study aimed to determine the changing incidence of HPV in patients with OPSCC in the period 1980-2009 and its influence on survival. PATIENTS AND METHODS: We randomly sampled 158 patients from a cohort of 828 patients with OPSCC stratified by decade (1980-1989, 1990-1999, 2000-2009). Formalin-fixed paraffin-embedded material was tested for HPV DNA by SPF-10 polymerase chain reaction (PCR) and immunohistochemically stained for p16 and p53. RESULTS: DNA from 146 patients was suitable for HPV detection. HPV DNA was detected in 13/47 (28%), 18/47 (38%), and 20/52 (38%) patients in the cohorts of 1980-1989, 1990-1999, and 2000-2009, respectively (p-value for trend=0.269). Lack of further increase during the most recent decade is inconsistent with the rising incidence and higher prevalence reported in other Western countries. Patients with HPV-positive OPSCC had a better survival in spite of higher tumor stage.

Robust BRCA1-like classification of copy number profiles of samples repeated across different datasets and platforms.

Breast cancers with BRCA1 germline mutation have a characteristic DNA copy number (CN) pattern. We developed a test that assigns CN profiles to be 'BRCA1-like' or 'non-BRCA1-like', which refers to resembling a BRCA1-mutated tumor or resembling a tumor without a BRCA1 mutation, respectively. Approximately one third of the BRCA1-like breast cancers have a BRCA1 mutation, one third has hypermethylation of the BRCA1 promoter and one third has an unknown reason for being BRCA1-like. This classification is indicative of patients' response to high dose alkylating and platinum containing chemotherapy regimens, which targets the inability of BRCA1 deficient cells to repair DNA double strand breaks. We investigated whether this classification can be reliably obtained with next generation sequencing and copy number platforms other than the bacterial artificial chromosome (BAC) array Comparative Genomic Hybridization (aCGH) on which it was originally developed. We investigated samples from 230 breast cancer patients for which a CN profile had been generated on two to five platforms, comprising low coverage CN sequencing, CN extraction from targeted sequencing panels (CopywriteR), Affymetrix SNP6.0, 135K/720K oligonucleotide aCGH, Affymetrix Oncoscan FFPE (MIP) technology, 3K BAC and 32K BAC aCGH. Pairwise comparison of genomic position-mapped profiles from the original aCGH platform and other platforms revealed concordance. For most cases, biological differences between samples exceeded the differences between platforms within one sample. We observed the same classification across different platforms in over 80% of the patients and kappa values of at least 0.36. Differential classification could be attributed to CN profiles that were not strongly associated to one class. In conclusion, we have shown that the genomic regions that define our BRCA1-like classifier are robustly measured by different CN profiling technologies, providing the possibility to retro- and prospectively investigate BRCA1-like classification across a wide range of CN platforms.