Skip to main content

Molecular analyses of unselected head and neck cancer cases demonstrates that human papillomavirus transcriptional activity is positively associated with survival and prognosis

Abstract

Background

Human papillomavirus DNA detection in head and neck squamous cell carcinoma has been linked to improved patient prognosis. The main aims of the study was to test the hypotheses that HPV16 E6/E7 oncogene and p53 function within tumours were associated with the widely reported improved patient survival and prognosis in head and neck cancer.

Methods

HPV16 DNA, mRNA and p53 mRNA presence were analysed in a prospective study of 42 unselected HNSCC patients; correlating the data with patient age, tumour staging/grade, treatment response, disease recurrence and survival.

Results

HPV16 DNA and HPV16 mRNA were present in 45.2 % and 21.4 % of patients, respectively. There was a significant positive association between the detection of HPV16 E6/E7 mRNA and p53 mRNA (p = 0.032), but this was not replicated for HPV16 DNA. Five-year disease free survival for the whole cohort was 63 % (CI 52.5–73.5 %). Multivariable analysis revealed only HPV16 E6/E7 mRNA expression to have significant prognostic influence (p = 0.04).

Conclusions

Our study suggests that HPV16 oncogenic transcriptional activity within HNSCC tumours is associated with improved patient survival and better prognosis in a German population. Simple HPV DNA detection alone did not demonstrate this association. The significant association of full-length (wild-type) p53 with HPV16 E6/E7 mRNA is further evidence for a functional relationship, which could contribute to the widely reported improved survival and prognosis. Larger studies are required to validate the frequency of HPV16 mRNA expression in HNSCC.

Peer Review reports

Background

Head and neck squamous cell carcinoma (HNSCC), with a worldwide incidence rate of approximately 500,000 new diagnoses per annum, is the 6th most common form of cancer [1]. Mortality has not improved substantially over the last few decades [2], due to late diagnosis (75 % of cases) and/or recurrent primary malignancies, and remains at 40–50 % at 5 years [3, 4]. Risk factors for the development of HNSCCs include tobacco and alcohol use as well as oral human papillomavirus (HPV) infection [5].

Human papillomavirus (HPV) has an important pathogenic role for a subset of HNSCC [58], leading to the grouping of oropharyngeal tumours into HPV-negative and HPV-positive cancers. Meta-analyses have been used to assess the potential worldwide prevalence of HPV in HNSCC, resulting in an estimate of 26–35.4 % [7, 8]. A recent systematic review of oropharyngeal subsite cancers and HPV prevalence found a significant increase when comparing studies recruiting patients before 2000 and after 2005 respectively (40.5 % [CI 35.1–46.1 %] versus 72.2 % [CI 52.9–85.7]) [9].

HPV status of HNSCC has been shown to affect overall patient prognosis, HPV-positive tumours have an improved response to therapy as reflected in a significant reduction in the risk of progression and improvement in overall survival rate compared to HPV-negative malignancies, irrespective of the treatment regime used [1013]. However, a certain number of patients with HPV-positive tumours have similar clinical outcomes to those seen in HPV-negative disease, the latter being linked to other prognostic factors such as increased EGFR expression [14] or smoking status [15].

As the fraction of HPV-positive head and neck cancers increases, this may highlight the importance of understanding the underlying mechanisms responsible for oncogenesis and the differences between HPV-related and HPV-unrelated disease that determine clinical outcome [16, 17].

The detection of HPV DNA alone, in the absence of evidence for viral gene expression, is not unequivocal molecular evidence that HPV infection either causes or promotes malignant progression in a lesion. In order to address this point, we have used a highly sensitive PCR detection method, together with a commercially available HPV genotyping kit (Linear Array HPV genotyping test, Roche Diagnostics Ltd., UK), to determine the presence of high- and low-risk HPV subtype L1 DNA in HNSCC tumours. Furthermore, we employed a separate PCR assay to detect HPV16 E6/E7 DNA and mRNA expression [18]. The ability of HPV 16 E6 to bind and promote the degradation of p53 has been suggested as one reason for an absence of disruptive p53 mutations in HPV 16 positive HNSCCs [19]. Full-length p53 RNA expression was therefore examined in age-matched groups of HPV16 E6/E7 expressing and non-expressing tumours. Patient data was collected, including patient age, sex, tumour staging, treatment response, recurrence and mortality. The relationship between the detection of HPV DNA, RNA, p53 expression and clinical data was then analysed. To further investigate discordant results, subgroup analysis was undertaken utilizing p16INK4a immunohistochemistry.

Although p16INK4A is the most commonly used biomarker for HPV positive head and neck cancer, specific limitations need to be observed. There is a subgroup of HPV- HNSCCs where expression of p16INK4A could lead to erroneous classification as HPV positive. Immunohistochemical staining itself can produce false positives i.e. intra-observer/inter-observer variation of background stain. In addition, a small but significant cohort of the false positive HPV- samples have p16INK4A mutations that may account for accumulation of inactive p16INK4A [15].

Methods

Patients and specimens

Clinical samples were obtained from unselected patients attending the Bielefeld Academic Teaching Hospital, Department of Otorhinolaryngology, Bielefeld, Germany between 2008 and 2009. All patients gave written informed consent for the study, which was approved by the local research ethics committee (Ethik-Kommission der Ärztekammer Westfalen-Lippe). A sample of the tumour was taken for histological analysis and the remainder was snap-frozen in liquid nitrogen and transported to the UK on dry ice prior to DNA and RNA extraction. A consultant histopathologist with expertise in head and neck pathology reviewed each tissue block to ensure adequate tumor sampling. All experiments were performed in the Department of Pathology at the University of Cambridge.

The study comprised 42 HNSCC tissue samples (patient mean age = 62.9 years, 95 % CI 59.11–66.01, range 46–87). Tumour staging was classified using the TNM classification of malignant tumours [20]. Post-operative follow-up was conducted at regular intervals for a period of five years, assessing response to treatment, disease recurrence and patient mortality. Clinical and histopathological features of the 42 HNSCC patients are shown in Table 1.

Table 1 Clinical, histopathological and follow-up data for 42 HNSCC patients

DNA extraction

Samples (max. 25 mg) were disrupted in a Bullet Blender™ (Next Advance, Averill Park, USA) in 300 μl digestion mix (10 mM Tris, pH 7.5; 10 mM EDTA; 0.5 % SDS; 200 μg/ml Proteinase K) for 5 min and then incubated overnight at 37 °C. Following Proteinase K inactivation at 56 °C for 10 min, the lysate was subjected to a phenol:chloroform extraction (1:1:1 volume) and the DNA precipitated from the supernatant with 1 ml 100 % ethanol. The DNA was then centrifuged (13,000 rpm, 4 °C, 20 min), the pellet washed with 70 % ethanol, air dried and re-suspended in 200 μl PBS. RNase digestion and total genomic DNA isolation was then performed using the DNeasy Blood and Tissue Kit (QIAGEN Ltd, UK), according to the manufacturer’s instructions. DNA was eluted with purified (deionised, double-distilled) H2O, and after purity and concentration were ascertained using a Nanodrop™ 1000 spectrophotometer, stored at −20 °C until use.

RNA extraction and cDNA synthesis

Samples (max. 25 mg) were disrupted in a Bullet Blender™ (Next Advance, Averill Park, USA) in 400 μl TRIzol reagent (Invitrogen Ltd, UK) for 5 min. A further 600 μl of TRIzol was added to each sample and the lysate subjected to a chloroform extraction (200 μl). The supernatant was precipitated with 0.5 ml isopropanol o/n at −20 °C. The RNA was then centrifuged (13,000 rpm, 4 °C, 20 min), the pellet washed with 70 % ethanol, air dried and re-suspended in 50 μl H2O. DNase I digestion (Invitrogen Ltd, UK) was then performed prior to RNA clean-up using a PureLink™ RNA Mini Kit (Invitrogen Ltd, UK), according to the manufacturer’s instructions. RNA was eluted with purified (deionised, double-distilled) H2O, purity and concentration was ascertained using a Nanodrop™ 1000 spectrophotometer.

Up to 5 μg of total RNA was reverse transcribed using BioScript Reverse Transcriptase (Bioline Ltd, London, UK) according to the manufacturer’s instructions. Briefly, RNA was incubated with random primers (Invitrogen Ltd, UK) at 70 °C for 5 min before being placed on ice for 5 min. Following the addition of reaction buffer, dNTPs and reverse transcriptase, cDNA synthesis was performed under the following conditions: 25 °C for 5 min, 42 °C for 30 min (synthesis step) and 70 °C for 10 min (stops reaction). Samples were stored at −20 °C prior to use. Reverse transcriptase negative samples acted as negative controls for subsequent analyses.

PCR methods

L1 single step DNA PCR analysis

A PCR assay using the PGMY09/11 L1 consensus primer set was performed as described previously [21]. The PCR cycling conditions were as follows; denaturing step of 95 °C for 5 min, followed by 40 cycles of 95 °C for 1 min, 55 °C for 1 min and 72 °C for 1 min. This was followed by a final extension period of 10 min at 72 °C. Agarose gel electrophoresis was used to confirm the presence/absence of bands specific for both HPV and human β-globin. Human placental DNA (Sigma-Aldrich, UK), containing 1000 copies of the plasmid pSP64-HPV16 [22], was used as a positive control for this and subsequent HPV16 PCRs of genomic DNA.

L1 nested DNA PCR and direct cycle sequencing

PCR reactions that were negative following amplification with the PGMY09/11 L1 consensus primers were subjected to further PCR amplification using the GP5+/GP6+ primer pair as described previously [23]. The PCR cycling conditions were as follows; denaturing step of 95 °C for 5 min, followed by 30 cycles of 95 °C for 1 min, 40 °C for 2 min and 72 °C for 1.5 min. This was followed by a final extension period of 10 min at 72 °C.

Positive bands identified following agarose gel electrophoresis were excised, the DNA purified using QiaQuick Gel Extraction columns according to the manufacturer’s instructions (QIAGEN Ltd, UK) and sequenced directly (Geneservice Ltd, UK) using the GP5+/GP6+ primers. The sequences were then aligned with known HPV types (NCBI Basic Local Alignment Search Tool).

PGMY-line blot assay/Linear Array HPV genotyping test

The procedure was carried out according to the manufacturer’s instructions and as previously described [24]. Briefly, PCR amplification was carried out with LA HPV GT primers as provided: Each 100 μl reaction consisted of 50 μl working master mix containing MgCl2, KCl, Amplitaq Gold DNA polymerase, uracil-N-glycosilase, deoxynucleotides (dNTPs) and biotinylated PGMY and β-globin primers together with 50 μl of DNA sample. DNA templates were titrated to a concentration of 2–4 ng/μl, i.e. 100–200 ng template DNA per reaction. The Applied Biosystems gold-plated 96-Well GeneAmp PCR System 9700 was programmed as follows: 50 °C for 2 min, 95 °C for 9 min and 40 cycles of 95 °C for 30 s, 55 °C for 1 min, 72 °C for 1 min and finally, at 72 °C for 5 min before holding indefinitely at 72 °C. The 40 cycles had a ramp rate set at 50 %.

The amplicons were denatured and hybridised to a strip containing specific probes for 37 HPV genotypes and β-globin reference lines before undergoing stringent washes.

Colorimetric determination with a Linear Array Detection Kit: the colour change reaction was from streptavidin-horseradish peroxidase mediated precipitation of working substrate. Positive reactions appeared as blue lines on the strip. The strips were interpreted using the HPV reference guide provided.

E6/E7 DNA PCR analysis

A PCR assay using primers specific for HPV16 E6/E7 was performed as described previously [25]. The PCR cycling conditions were as follows; denaturing step of 94 °C for 5 min, followed by 40 cycles of 94 °C for 1 min, 60 °C for 45 s and 72 °C for 1 min. This was followed by a final extension period of 10 min at 72 °C.

E6/E7 cDNA PCR analysis

cDNA generated from extracted HNSCC mRNA was subjected to E6/E7 PCR analysis using the same primers and conditions as used for DNA samples. Expected amplicon sizes were as follows; unspliced RNA, 406 bp: E6*I, 224 bp: E6*II, 107 bp. A PCR for GAPDH was performed, to confirm that the RNA was suitable for amplification, using primers that have been previously described [26]. The PCR cycling conditions were as follows; denaturing step of 95 °C for 5 min, followed by 40 cycles of 95 °C for 1 min, 60 °C for 1 min and 72 °C for 1 min. This was followed by a final extension period of 5 min at 72 °C.

p53 cDNA PCR analysis

p53 amplification was performed on 9 HNSCC samples in which E6/E7 transcripts had been successfully amplified (mean age 63.7 years, 95 % CI 57.41–69.93, range 50–74), together with 12 HNSCC where no E6/E7 transcripts were detectable (mean age 63.2 years, 95 % CI 57.83–68.87, range 46–78). Following unsuccessful amplification of full-length p53 using established primers [27], possibly as a result of RNA fragmentation or degradation, five primer pairs were designed to span the entire open reading frame. PCR was carried out on an Applied Biosystems Veriti using Platinum Taq polymerase (Invitrogen Ltd, UK) and comprised 1x Platinum Taq buffer, 1.5 mM MgCl2, 200 μM each dNTP, 200 nM each primer, and 1 U of Platinum Taq polymerase in a 50 μl reaction. The PCR cycling conditions for all primer pairs were as follows; denaturing step of 94 °C for 2 min, followed by 40 cycles of 94 °C for 1 min, 55 °C for 1 min and 72 °C for 1 min. This was followed by a final extension period of 10 min at 72 °C. cDNA generated from primary human foreskin keratinocytes was used as a positive control in this assay.

Immunohistochemistry analysis

p16INK4a protein is an inhibitor of cyclin dependent kinase and has increased expression with elevated levels of HPV E7. Selected FFPE sections (5 μm) were deparaffinized and antigen target retrieval was performed with citrate-buffer boil (0.1 M Sodium Citrate and 0.1 M Citric acid). Immunohistochemistry was performed as previously described using a mouse monoclonal antibody (BD Biosciences, USA [18].

Statistical analysis

Statistical calculations were performed using SPSS Version 17 (Chicago, IL, USA). Significance of the association between the detection of HPV16 L1 DNA, HPV16 E6/E7 DNA, HPV16 E6/E7 mRNA, p53 mRNA and P16INK4a was calculated using a two-sided Fisher’s exact test. A Student’s t-test was used to compare mean ages. Rates of disease-specific and disease free survival were estimated by means of the Kaplan–Meier method and were compared by the log-rank test. Univariable and multivariable models were developed using Cox regression to investigate size of tumor, nodal status, smoking (never, current or former if smoking cessation >10 years), adjuvant chemoradiotherapy, age, histological differentiation, primary subsite, HPV and p53 status as potential predictors of outcome. Time dependent co-variants were investigated to identify concordance with the proportional hazards assumption.

Results

Detection and identification of HPV L1 sequences

Sensitive detection of HPV L1 sequences was performed using a combination of PGMY09/11, Linear Array HPV genotyping test and nested PCR approaches, as previously described [18]. PGMY09/11 PCR detected HPV in 2/42 samples (4.8 %), Linear Array in 10/42 samples (23.8 %) and the nested PCR approach in 29/42 samples (69.0 %). Direct sequencing confirmed the HPV types detected from PGMY09/11 and GP5+/GP6+ nested PCR approaches. Combined results from L1 detection methods showed the presence of HPV 6 in 9/42 (21.4 %) and HPV 16 in 19/42 (45.2 %) samples. Other HPV sequences were detected in 7 samples, but remained unidentifiable due to mixed sequencing traces. HPV 11, HPV 26 and HPV 40 were each found on one occasion in the external meatus, tonsil and oral cavity regions respectively. The distribution of the HPV type 16 can be seen in Table 2. There was no evident association between the detection of HPV 6 or 16 DNA and tumour site.

Table 2 Distribution of HPV16 in HNSCCs at different anatomical locations

Detection of HPV16 E6/E7 DNA and mRNA

HPV16 E6/E7 DNA was detected in 21/42 samples (50 %). There was a significant association between the detection of HPV16 by E6/E7- and L1-based PCR detection methods (p = 0.0016, two-sided Fisher’s exact test). HPV16 E6/E7 transcripts were detected in 9/42 samples (21.4 %). No significant association was seen between the presence of HPV16 L1 DNA and E6/E7 transcript detection (p = 1) or the presence of HPV16 E6/E7 DNA and E6/E7 transcript detection (p = 1). Of the samples expressing HPV 16 E6/E7, the primary lesion was located in the larynx (x6), hypopharynx (x1), oropharynx (x1) and oral cavity (x1). There was no association between anatomical location and HPV16 E6/E7 expression.

Clinical evaluation and follow-up

Forty-two patients with primary HNSCC were recruited to this study with a 5-year disease specific survival of 72 % (95 % CI 66–76 %). At the time of surgery, the mean patient age was 62.9 years (range 46–87). The median follow-up time was 53 months (range 45–61 months).

Expression of HPV16 mRNA E6/E7 was found to be a significant predictor of disease free survival on univariable analysis but this was not the case when stratifying by presence of HPV16 L1 or E6/7 DNA (Figs. 1 and 2). Disease specific survival shows the same clear trend but does not quite achieve significance. Multivariable analysis investigated the influence of tumour size (T1-2 versus T3-4), nodal status, smoking, primary surgery, primary radiotherapy, primary surgery + adjuvant radiotherapy, age, histological differentiation (well/moderate/poor SCC) and subsite. Forward model selection confirmed mRNA expression as an independent prognostic factor (Table 3).

Fig. 1
figure 1

Disease free survival (DFS) stratified by HPV mRNA/DNA (Log- Rank multiple regression analysis p = 0.04 and p = 0.68, respectively)

Fig. 2
figure 2

Disease specific survival (DSS) stratified by HPV mRNA/DNA (Log- Rank multiple regression analysis p = 0.17 and p = 0.66 respectively)

Table 3 Log-Rank analysis for disease free survival (DFS) and disease specific survival (DSS)

There was no significant difference in the mean age, irrespective of the method of HPV detection used, between HPV-positive and HPV-negative patients.

Detection of p53 mRNA

Amplicons spanning the entire open reading frame of p53 were detected in 6/9 (66.7 %) of HPV16 E6/E7 transcript positive HNSCC samples but only in 2/12 (16.7 %) of the age-matched control group, where E6/E7 transcripts were undetectable (Table 4). There was a significant association between E6/E7 expression and the presence of full-length p53 (p = 0.032, Fisher’s exact test; OR 10). Full-length p53 expression was not found to correlate with anatomical site, HPV16 L1/E6/E7 DNA or p16INK4a.

Table 4 p53 detection analysis

Discussion

In this study, we report the prevalence of HPV16 DNA and mRNA in 42 unselected HNSCC patients, significantly correlating the presence of HPV16 mRNA expression with the successful amplification of p53 (despite a small sample size). Overall, HPV16 DNA was present in 45.2 % of samples, with detectable E6/E7 mRNA in 21.4 %. In two-thirds of the samples expressing E6/E7, p53 expression was also seen. This study has enabled the classification of HNSCCs into subsets based on both their HPV 16 DNA and mRNA status.

HPV16 and HPV6 were the predominant subtypes detected, being present in 47.6 % and 21.4 % of the samples, respectively. Meta-analysis has revealed that HPV16 is present in the majority of HPV-positive HNSCCs, with HPVs 18, 6 and 11 involved to a lesser extent and an overall prevalence of HPV in HNSCC of 25.9 % 27. Large epidemiological studies also suggest that the oropharyngeal site is most associated with high-risk oncogenic HPV subtypes [28, 29].

The higher HPV prevalence observed in this series may be due to the combined use of both sensitive PCR approaches and a commercial assay (Roche Linear Array HPV genotyping test) for the detection of HPV sequences. However, several studies within the meta-analysis found a comparable prevalence, indicating that our findings are within the expected range. There were no cases of HPV 18 positive HNSCC in our study. The predominance of HPV 18 in adenocarcinomas of the cervix indicates a tropism towards glandular tissue [30, 31]. In the context of HNSCC, adenocarcinomas are rare and occur mainly in the salivary glands [32]. No cases of adenocarcinoma were present in this study. The DNA of HPV 26 and HPV 40 were each found on one occasion. HPV 26 has been classified as a probable high-risk type, whereas HPV 40 is a low-risk type [33]. A clear correlation existed between the detection of HPV 16 DNA by L1- and E6/E7-specific PCRs. However, despite 32 samples being concordant, it is noted that 10 samples were positive in one of the PCRs only. This is important for both the estimation of the overall contribution of HPV to HNSCC carcinogenesis and to suggest ‘indirect’ mechanisms by which lesions may arise. The detection of HPV L1 sequences, by Linear Array or PCR amplification, is frequently performed and used as an indicator of HPV-related disease. This has allowed classification of HNSCC tumours into HPV-positive and HPV-negative groups. Several studies have documented significantly better overall outcome in HNSCC patients with tumours positive for HPV DNA compared to their HPV-negative counterparts [911, 34].

One of the most significant findings from this study is that p16 positivity does not correlate with HPV E6/E7 mRNA expression. p16INK4A immunohistochemistry outcomes are classified as positive or negative depending on a variable spectrum of detection (ranging from true negative to background staining to false positive to true positive). A systematic review by Grønhøj Larsen et al. [35] identified thirty-nine different studies published between January 1980 to October 2012 which reported sensitivity and specificity of p16 detection (with respect to HPV detection of DNA by PCR). The results ranged from 73 % to 100 % for sensitivity, and 46 % to 100 % for specificity. The definition of “positive” overexpression of p16 is also determined by different % staining cut-off levels such as  ≥ 5 %, ≥5 % but ≤ 70 %, ≥70 %, ≥75 % or even less specific verbal definitions such as “diffuse and strong nuclear cytoplasmic staining”. There are also no standardisations of the anti-p16 monoclonal antibodies used and different clones have different staining avidities and consequently sensitivities and specificities. Our study used a definition of staining ≥70 % of cell cytoplasm and nuclei, seen per field, which stained strongly with the antibody.

This study demonstrates a significant correlation between HPV status and prognosis, but this was restricted to mRNA expression alone. Although this information is important, recent studies suggest further biomarkers may also be required as a subset of the HPV-positive cohort can respond in a manner normally associated with HPV-negative tumors [16, 17]. It is conceivable that HNSCCs positive for L1 but negative for E6/E7 DNA have lost the selective pressure to retain oncogene expression following other genomic events. Conversely, HNSCCs that are E6/E7 positive may have lost L1 sequences, for example following integration, making L1 detection unreliable in determining the role of HPV infection at the time of diagnosis.

The detection of viral oncogene transcripts remains the gold standard by which tumours are classified as HPV positive [36]. In our study, transcript analysis revealed 21.4 % of samples expressed HPV16 E6/E7 mRNA. Although the patients with HPV16 DNA positive tumours tended to be younger (mean age for HPV L1 DNA positive v negative: 60.4 v 64.9 years; mean age for HPV 16 E6/E7 DNA positive v negative: 60.3 v 65.4 years), there was no difference in age between patients with E6/E7 mRNA positive tumours and their mRNA negative counterparts (mean age 63.7 v. 62.6 years). The detection of L1 and E6/E7 DNA in the absence of E6/E7 mRNA expression may reflect virus present as either a passenger or as latent infection whereas the presence of E6/E7 mRNA provides strong molecular evidence for ongoing HPV transcription (and by implication, translation into proteins), which caused or promoted malignant progression in the lesion. HPV-associated cervical carcinomas commonly have integrated HPV genomes in a single copy [3739]. Therefore, a tumour containing transcriptionally active HPV may well have multiple copies of mature mRNA transcribed from a single copy of HPV DNA. An actively transcribed gene can have mRNA copy numbers several orders of magnitude higher than the DNA from which it originates. This may explain our findings that 4 of the samples were mRNA positive but DNA negative i.e. reflecting insufficiently sensitive DNA detection assays. Sotlar et al. have previously quantified the sensitivity of RT-PCR against DNA-PCR for HPV E6/E7; their data suggest a benefit in favour of the former test with a ten-fold increased detection rate for HPV16 subtype [25].

HPV status, as determined solely by DNA-based detection methods, may be insufficient to predict patient response to treatment and overall outcome. Our study shows that a proportion of HPV DNA positive tumours contain transcriptionally active genomes and this may have more important clinical significance. Other potential prognostic markers, in addition to extensive tobacco and alcohol use, include mutation/over-expression of p53, elevated EGFR levels and over-expression of cyclin D (CCND1) [40, 41]. p53 expression analysis revealed a significant association between HPV 16 E6/E7 expression and the presence of full-length p53. Many studies have revealed inverse correlations between disruptive p53 mutations and HPV in HNSCC [6, 27, 42], mutation of p53 has been shown to result in mis-splicing of transcripts [27, 35, 43] and this may explain the partial amplification seen in the HPV negative samples analysed. However, as a significant number of p53 mutations do not affect full-length expression, this observation may not readily explain all the findings in this study [44]. An alternative method of analysis that may elucidate this in future studies would involve direct sequence analysis of theTP53 gene (exons 5–9) [45]. HPV positive tumours characteristically have wild-type p53; mutations in p53 are not required for tumour development because high-risk HPV E6 targets p53 for degradation [46]. Defects in apoptosis associated with mutated p53 can potentially lead to treatment resistant tumours and may partly explain the improved response seen in HPV associated disease [47].

The five-year clinical observation period for this study, in combination with other studies [4850], suggests the need to place HPV16 mRNA status as the preferred stratification modality when designing future treatment regimes or clinical trials. This has particular importance when we consider the trend towards “de-intensification” of existing chemotherapeutic regimes in order to decrease treatment related co-morbidity [6, 46, 51].

Conclusions

In conclusion, the evidence from this small study suggests that HPV status, as determined solely by DNA-based detection methods, is insufficient to predict patient response to treatment and overall outcome in a German population. Our data supports the hypothesis that HPV positive tumours containing transcriptionally active genomes are a more reliable indicator of disease free survival than HPV DNA detection alone. Therefore larger studies are required to establish the frequency of HPV16 mRNA expression in HNSCC. This will determine whether this group, with or without p53 expression, represents a subset of tumours that have improved treatment response.

Abbreviations

HNSCC, Head & neck squamous cell carcinoma; HPV, human papillomavirus; PCR, polymerase chain reaction

References

  1. de Martel C, Ferlay J, Franceschi S, Vignat J, Bray F, Forman D, Plummer M. Global burden of cancers attributable to infections in 2008: a review andsynthetic analysis. Lancet Oncol. 2012;13:607–15.

    Article  PubMed  Google Scholar 

  2. Forastiere A, Koch W, Trotti A, Sidransky D. Head and neck cancer. N Engl J Med. 2001;345:1890–900.

    Article  CAS  PubMed  Google Scholar 

  3. Hitt R, López-Pousa A, Martínez-Trufero J, Escrig V, Carles J, Rizo A, Isla D, Vega ME, Martí JL, Lobo F, Pastor P, Valentí V, et al. Phase III study comparing cisplatin plus fluorouracil to paclitaxel, cisplatin, and fluorouracil induction chemotherapy followed by chemoradiotherapy in locally advanced head and neck cancer. J Clin Oncol. 2005;23:8636–45.

    Article  PubMed  Google Scholar 

  4. Brockstein B, Haraf DJ, Rademaker AW, Kies MS, Stenson KM, Rosen F, Mittal BB, Pelzer H, Fung BB, Witt ME, Wenig B, Portugal L, et al. Patterns of failure, prognostic factors and survival in locoregionally advanced head and neck cancer treated with concomitant chemoradiotherapy: a 9-year, 337-patient, multi-institutional experience. Ann Oncol. 2004;15:1179–86.

    Article  CAS  PubMed  Google Scholar 

  5. D'Souza G, Kreimer AR, Viscidi R, Pawlita M, Fakhry C, Koch WM, Westra WH, Gillison ML. Case–control study of human papillomavirus and oropharyngeal cancer. N Engl J Med. 2007;356:1944–56.

    Article  PubMed  Google Scholar 

  6. Masterson L, Moualed D, Liu ZW, Howard J, Dwivedi R, Benson R, Tysome JR, Sterling JC, Jani P, Sudhoff H, Goon PKC. De-escalation treatment protocols for HPV associated oropharyngeal squamous cell carcinoma: a systematic review and meta-analysis of current clinical trials. Eur J Cancer. 2014;50(15):2636–48.

    Article  PubMed  Google Scholar 

  7. Kreimer AR, Clifford GM, Boyle P, Franceschi S. Human papillomavirus types in head and neck squamous cell carcinomas worldwide: a systematic review. Cancer Epidemiol Biomarkers Prev. 2005;14:467–75.

    Article  CAS  PubMed  Google Scholar 

  8. Termine N, Panzarella V, Falaschini S, Russo A, Matranga D, Lo Muzio L, Campisi G. HPV in oral squamous cell carcinoma vs head and neck squamous cell carcinoma biopsies: a meta-analysis (1988–2007). Ann Oncol. 2008;19:1681–90.

    Article  CAS  PubMed  Google Scholar 

  9. Mehanna H, Beech T, Nicholson T, El-Hariry I, McConkey C, Paleri V, Roberts S. The prevalence of human papillomavirus in oropharyngeal and nonoropharyngeal head and neck cancer: a systematic review and meta-analysis of trends by time and region. Head Neck. 2013;35:747–55.

    Article  PubMed  Google Scholar 

  10. Fakhry C, Westra WH, Li S, Cmelak A, Ridge JA, Pinto H, Forastiere A, Gillison ML. Improved survival of patients with human papillomavirus-positive head and neck squamous cell carcinoma in a prospective clinical trial. J Natl Cancer Inst. 2008;100:261–9.

    Article  CAS  PubMed  Google Scholar 

  11. Licitra L, Perrone F, Bossi P, Suardi S, Mariani L, Artusi R, Oggionni M, Rossini C, Cantù G, Squadrelli M, Quattrone P, Locati LD, et al. High-risk human papillomavirus affects prognosis in patients with surgically treated oropharyngeal squamous cell carcinoma. J Clin Oncol. 2006;24:5630–6.

    Article  CAS  PubMed  Google Scholar 

  12. Lindquist D, Romanitan M, Hammarstedt L, Näsman A, Dahlstrand H, Lindholm J, Onelöv L, Ramqvist T, Ye W, Munck-Wikland E, Dalianis T. Human papillomavirus is a favourable prognostic factor in tonsillar cancer and its oncogenic role is supported by the expression of E6 and E7. Mol Oncol. 2007;1:350–5.

    Article  PubMed  Google Scholar 

  13. Masterson L, Moualed D, Masood A, Dwivedi RC, Benson R, Sterling JC, Rhodes KM, Sudhoff H, Jani P, Goon P. De-escalation treatment protocols for human papillomavirus-associated oropharyngeal squamous cell carcinoma. Cochrane Database Syst Rev 2014, Issue 2. Art. No.: CD010271. doi: 10.1002/14651858.CD010271.

  14. Kong CS, Narasimhan B, Cao H, Kwok S, Erickson JP, Koong A, Pourmand N, Le QT. The relationship between human papillomavirus status and other molecular prognostic markers in head and neck squamous cell carcinomas. Int J Radiat Oncol Biol Phys. 2009;74:553–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Romagosa C, Simonetti S, López-Vicente L, Mazo A, Lleonart ME, Castellvi J. p16(Ink4a) overexpression in cancer: a tumor suppressor gene associated with senescence and high-grade tumors. Oncogene. 2011;30(18):2087–97.

    Article  CAS  PubMed  Google Scholar 

  16. Chaturvedi AK, Engels EA, Anderson WF, Gillison ML. Incidence trends for human papillomavirus-related and -unrelated oral squamous cell carcinomas in the United States. J Clin Oncol. 2008;26:612–9.

    Article  PubMed  Google Scholar 

  17. Masterson L, Sorgeloos F, Winder D, Lechner M, Marker A, Malhotra S, Sudhoff H, Jani P, Goon PKC, Sterling JC. Deregulation of SYCP2 predicts early stage HPV+ oropharyngeal carcinoma – a prospective whole transcriptome analysis. Cancer Sci. 2015;106(11):1568–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Winder DM, Ball SL, Vaughan K, Hanna N, Woo YL, Fränzer JT, Sterling JC, Stanley MA, Sudhoff H, Goon PK. Sensitive HPV detection in oropharyngeal cancers. BMC Cancer. 2009;9:440.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Westra WH, Taube JM, Poeta ML, Begum S, Sidransky D, Koch WM. Inverse relationship between human papillomavirus-16 infection and disruptive p53 gene mutations in squamous cell carcinoma of the head and neck. Clin Cancer Res. 2008;14:366–9.

    Article  CAS  PubMed  Google Scholar 

  20. Sobin LH, Gospodarowicz MK, Wittekind CH. TNM Classification of Malignant Tumors. 7th ed. Oxford: Wiley-Blackwell; 2009.

    Google Scholar 

  21. Gravitt PE, Peyton CL, Alessi TQ, Wheeler CM, Coutlée F, Hildesheim A, Schiffman MH, Scott DR, Apple RJ. Improved amplification of genital human papillomaviruses. J Clin Microbiol. 2000;38:357–61.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Stanley MA, Browne HM, Appleby M, Minson AC. Properties of a non-tumorigenic human cervical keratinocyte cell line. Int J Cancer. 1989;43:672–76.

    Article  CAS  PubMed  Google Scholar 

  23. de Roda Husman AM, Walboomers JM, van den Brule AJ, Meijer CJ, Snijders PJ. The use of general primers GP5 and GP6 elongated at their 3’ ends with adjacent highly conserved sequences improves human papillomavirus detection by PCR. J Gen Virol. 1995;76:1057–62.

    Article  PubMed  Google Scholar 

  24. Woo YL, Damay I, Stanley M, Crawford R, Sterling J. The use of HPV Linear Array Assay for multiple HPV typing on archival frozen tissue and DNA specimens. J Virol Methods. 2007;142:226–30.

    Article  CAS  PubMed  Google Scholar 

  25. Sotlar K, Stubner A, Diemer D, Sotlar K, Stubner A, Diemer D, Menton S, Menton M, Dietz K, Wallwiener D, Kandolf R, Bültmann B. Detection of high-risk human papillomavirus E6 and E7 oncogene transcripts in cervical scrapes by nested RT-polymerase chain reaction. J Med Virol. 2004;74:107–16.

    Article  CAS  PubMed  Google Scholar 

  26. Eagle RA, Traherne JA, Ashiru O, Wills MR, Trowsdale J. Regulation of NKG2D ligand gene expression. Hum Immunol. 2006;67:159–69.

    Article  CAS  PubMed  Google Scholar 

  27. Hauser U, Balz V, Carey TE, Grénman R, Van Lierop A, Scheckenbach K, Bier H. Reliable detection of p53 aberrations in squamous cell carcinomas of the head and neck requires transcript analysis of the entire coding region. Head Neck. 2002;24:868–73.

    Article  PubMed  Google Scholar 

  28. Licitra L, Zigon G, Gatta G, Sánchez MJ, Berrino F, Working Group EUROCARE. Human papillomavirus in HNSCC: a European epidemiologic perspective. Hematol Oncol Clin North Am. 2008;22:1143–53.

    Article  PubMed  Google Scholar 

  29. Yu GP, Mehta V, Branovan D, Huang Q, Hashibe M, Zhang ZF, Schantz SP. Improved survival among patients with base of tongue and tonsil cancer in the United States. Cancer Causes Control. 2012;23:153–64.

    Article  PubMed  Google Scholar 

  30. Bosch FX, Manos MM, Muñoz N, Sherman M, Jansen AM, Peto J, Schiffman MH, Moreno V, Kurman R, Shah KV. Prevalence of human papillomavirus in cervical cancer: a worldwide perspective. International biological study on cervical cancer (IBSCC) Study Group. J Natl Cancer Inst. 1995;87:796–802.

    Article  CAS  PubMed  Google Scholar 

  31. Clifford GM, Smith JS, Plummer M, Munoz N, Franceschi S. Human papillomavirus types in invasive cervical cancer worldwide: a meta-analysis. Br J Cancer. 2003;88:63–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Canto MT, Devesa SS. Oral cavity and pharynx cancer incidence rates in the United States, 1975–1998. Oral Oncol. 2002;38:610–17.

    Article  PubMed  Google Scholar 

  33. Muñoz N, Bosch FX, de Sanjosé S, Herrero R, Castellsagué X, Shah KV, Snijders PJ, Meijer CJ. International Agency for Research on Cancer Multicenter Cervical Cancer Study Group. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003;348:518–27.

    Article  PubMed  Google Scholar 

  34. Kumar B, Cordell KG, Lee JS, Worden FP, Prince ME, Tran HH, Wolf GT, Urba SG, Chepeha DB, Teknos TN, Eisbruch A, Tsien CI, et al. EGFR, p16, HPV Titer, Bcl-xL and p53, sex, and smoking as indicators of response to therapy and survival in oropharyngeal cancer. J Clin Oncol. 2008;26:3128–37.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Grønhøj Larsen C, Gyldenløve M, Jensen DH, Therkildsen MH, Kiss K, Norrild B. Correlation between human papillomavirus and p16 overexpression in oropharyngeal tumours: a systematic review. Br J Cancer. 2014;110(6):1587–94.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Chung CH, Gillison ML. Human papillomavirus in head and neck cancer: its role in pathogenesis and clinical implications. Clin Cancer Res. 2009;15:6758–62.

    Article  CAS  PubMed  Google Scholar 

  37. Choo KB, Pan CC, Han SH. Integration of human papillomavirus type 16 into cellular DNA of cervical carcinoma: preferential deletion of the E2 gene and invariable retention of the long control region and the E6/E7 open reading frames. Virology. 1987;161:259–61.

    Article  CAS  PubMed  Google Scholar 

  38. Kim CJ, Um SJ, Hwang ES, Park SN, Kim SJ, Namkoong SE, Park JS. The antibody response to HPV proteins and the genomic state of HPVs in patients with cervical cancer. Int J Gynecol Cancer. 1999;9:1–11.

    Article  CAS  PubMed  Google Scholar 

  39. Park JS, Hwang ES, Park SN, Ahn HK, Um SJ, Kim CJ, Kim SJ, Namkoong SE. Physical status and expression of HPV genes in cervical cancers. Gynecol Oncol. 1997;65:121–29.

    Article  CAS  PubMed  Google Scholar 

  40. Castellsagué X, Alemany L, Quer M, Halec G, Quirós B, Tous S. HPV Involvement in Head and Neck Cancers: Comprehensive Assessment of Biomarkers in 3680 Patients. J Natl Cancer Inst. 2016; 108(6). doi: 10.1093/jnci/djv403.

  41. Quon H, Liu FF, Cummings BJ. Potential molecular prognostic markers in head and neck squamous cell carcinomas. Head Neck. 2001;23:147–59.

    Article  CAS  PubMed  Google Scholar 

  42. Nodzenski E, Brachman D, Mick R, Montag A, Graves D, Vokes EE, Weichselbaum RR. Human papilloma virus and p53 in head and neck cancer: clinical correlates and survival. Clin Cancer Res. 1996;2:755–62.

    PubMed  Google Scholar 

  43. Koch WM, Brennan JA, Zahurak M. p53 mutation and locoregional treatment failure in head and neck squamous cell carcinoma. J Natl Cancer Inst. 1996;88:1580–86.

    Article  CAS  PubMed  Google Scholar 

  44. Garritano S, Inga A, Gemignani F, Landi S. More targets, more pathways and more clues for mutant p53. Oncogenesis. 2013;2:e54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Masterson L, Winder DM, Marker A, Sterling JC, Sudhoff H, Moffat D, Goon PKC. Investigating the role of human papillomavirus in squamous cell carcinoma of the temporal bone. Head Neck Oncol. 2013;5:22.

    Google Scholar 

  46. Fakhry C, Gillison ML. Clinical implications of human papillomavirus in head and neck cancers. J Clin Oncol. 2006;24:2606–11.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Lowe SW, Bodis S, McClatchey A, Remington L, Ruley HE, Fisher DE, Housman DE, Jacks T. p53 status and the efficacy of cancer therapy in vivo. Science. 1994;266:807–10.

    Article  CAS  PubMed  Google Scholar 

  48. Shi W, Kato H, Perez-Ordonez B, Pintilie M, Huang S, Shi W, Kato H, Perez-Ordonez B, Pintilie M, Huang S, Hui A, O'Sullivan B, Waldron J, Cummings B, Kim J, Ringash J, Dawson LA, Gullane P, Siu L, Gillison M, Liu FF. Comparative prognostic value of HPV16 E6 mRNA compared with in situ hybridization for human oropharyngeal squamous carcinoma. J Clin Oncol. 2009;27:6213–21.

    Article  PubMed  Google Scholar 

  49. Fallai C, Perrone F, Licitra L, Pilotti S, Locati L, Bossi P, Orlandi E, Palazzi M, Olmi P. Oropharyngeal squamous cell carcinoma treated with radiotherapy or radiochemotherapy: prognostic role of TP53 and HPV status. Int J Radiat Oncol Biol Phys. 2009;75:1053–9.

    Article  PubMed  Google Scholar 

  50. Rotnáglová E, Tachezy R, Saláková M, Procházka B, Košl'abová E, Veselá E, Ludvíková V, Hamšíková E, Klozar J. HPV involvement in tonsillar cancer: prognostic significance and clinically relevant markers. Int J Cancer. 2011;129:101–10.

    Article  PubMed  Google Scholar 

  51. Tribius S, Ihloff AS, Rieckmann T, Petersen C, Hoffmann M. Impact of HPV status on treatment of squamous cell cancer of the oropharynx: what we know and what we need to know. Cancer Lett. 2011;304:71–9.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the patients from Bielefeld Academic Teaching Hospital (Bielefield, Germany) for volunteering clinical samples.

Funding

LM received a research grant from Cancer Research UK (C45051/A14962).

Authors’ contributions

DW, LM, SB, KV carried out all laboratory experiments and drafted the manuscript. ML, LS and HS supplied all samples, co-ordinated patient recruitment and helped to draft the manuscript. JS and PG conceived the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

All patients gave written informed consent for the study, which was approved by the local research ethics committee (LREC) (Ethik-Kommission der Ärztekammer Westfalen-Lippe).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Peter K. C. Goon.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Masterson, L., Winder, D.M., Ball, S.L.R. et al. Molecular analyses of unselected head and neck cancer cases demonstrates that human papillomavirus transcriptional activity is positively associated with survival and prognosis. BMC Cancer 16, 367 (2016). https://doi.org/10.1186/s12885-016-2398-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-016-2398-7

Keywords