Abstract
Background/Aim: High risk Human papillomavirus (hr-HPV) and smoking are independant risk factors for head and neck squamous cell carcinomas (HNSCC). While hr-HPV+ HNSCC has a better prognosis than smoking-associated HNSCC no systematic data are yet available about the combined risk. Patients and Methods: We performed a meta-analysis to assess the overall survival of HNSCC patients relative to the hr-HPV and smoking status. A literature review up to November 2019 was conducted in PubMed and Cochrane Library using the search terms ‘HPV, Smoking and HNSCC’. Results: Nine out of 748 articles were included, 1,436 out of 2,080 patients were hr-HPV+. The prevalence of hr-HPV+ smokers was 36%. The meta-analysis showed a significantly better 5-year overall survival for HPV+ non-smokers compared to smokers with risk ratio of 1.94 (95% confidence intervaI=1.46-2.58). Conclusion: Smoking is a negative prognostic factor for overall survival in patients with hr-HPV+ HNSCC and should thus be an important part of staging and treatment.
Head and neck cancer is the sixth most common cancer worldwide. Approximately 630,000 new cases are diagnosed annually resulting in more than 350,000 deaths per year. Head and neck cancer accounts for approximately 3-4% of all malignancies in Europe and the United States (1, 2). More than 90% of these cases are head and neck squamous cell carcinomas (HNSCC) (3). Despite their anatomic proximity, HNSCC are a heterogeneous group of tumor entities with differences in etiology, diagnostic and therapeutic approaches as well as prognosis (3). The risk of developing HNSCC is associated with several factors. It is well accepted that HNSCC is etiologically associated with excessive tobacco and alcohol consumption. Both smoking and drinking are important risk factors and have a synergistic effect (4). Furthermore, infection with high-risk human papilloma viruses (hr-HPV) have been identified as a the second major causative agent (5). Approximately 40% of all HNSCCs are reported to be hr-HPV+, as well as approximately 70% of oropharyngeal squamous cell carcinomas (OPSCCs) (6). HPV-16, HPV-18, HPV-31 and HPV-33 are classified as hr- HPV. HPV-16 is the most common HPV subtype and is associated with 86-92% of all HPV+ HNSCCs (5, 6). From an epidemiological perspective, the incidence of smoking-associated HNSCC is declining, while that of hr-HPV-induced cancer is increasing, especially in younger patients (7).Moreover, the underlying cause of HNSCC, i.e. toxins such as tobacco and alcohol, or hr-HPV+, leads to significant differences not only in etiology, but also in prognosis; Hr-HPV+ status being associated with an improved prognosis. Accordingly, several studies suggested a survival difference of up to 26.9% between hr-HPV+ and HPV− HNSCC (8). Because hr-HPV+ tumors differ from and have a more favorable prognosis than HPV− HNSCC, the recently released Eighth Edition of the American Joint Committee on Cancer created a separate staging algorithm for hr-HPV+ OPSCC, distinguishing it from HPV− OPSCC. Ηr-HPV+ OPSCC that used to be designated as ‘advanced stage’ is now categorized as ‘lower stage’, which gives a much more accurate survival prediction (9, 10). Smoking, on the other hand, is correlated with poor prognosis. The negative impact of smoking on survival has been well documented (11-13). The 5-year overall survival (OS) rate for never-smokers compared with current-/former smokers is 77.7% vs. 57.5/58.1% respectively (14). Consequently, our goal was to assess how these two independent prognostic factors interact with each other. To our knowledge, no meta-analysis of survival according to combined HPV/smoking status has been performed. The aim of the current study was to perform a meta-analysis on the survival of patients in these distinct groups and to compare smoking and non-smoking hr-HPV+ groups.
Patients and Methods
Literature search. We conducted a comprehensive search looking for published literature that evaluated survival of HNSCC patients according to both hr-HPV and smoking status in the PubMed (http://www.ncbi.nlm.nih.gov/pubmed) and Cochrane Library (https://www.cochranelibrary.com) databases. This search included articles published up to 12/26/2019. The following search terms were used: ‘HPV, smoking, head and neck cancer’. Additionally, we checked references quoted in original or review articles that may not have been found in the database during the initial literature search. Our search strategy was performed according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses and logged in the Prospective Register of Systematic Reviews (15). We screened all search results and included studies of patients with HNSCC, investigating survival rates according to HPV and smoking status. When abstracts met our inclusion criteria, the articles were read in detail. Exclusion criteria were: Patient survival information not reported, evaluation of survival based on only one marker (HPV or smoking) or the survival of HPV and smoking status was not combined (e.g. both parameters were assessed separately), non-HNSCC primary cancer (e.g. pre-cancerous lesions, esophageal cancer, lung cancer, skin cancer), cell culture or animal models, systematic reviews or case reports and articles written in a language other than English or German. Finally, this meta-analysis included studies with the following criteria: (i) Proportion of the HNSCC subgroups smoking versus nonsmoking; (ii) numerical survival data for the subgroups [5-year overall survival (OS) or Kaplan-Meier curves of OS by subgroup].
Study selection and data extraction. Two Authors (A.M. and V.V) independently selected the studies using the inclusion and exclusion criteria. Relevant data from all eligible publications were extracted: First author’s name, study design, study time, sample number, gender, tumor location, HPV detection method, smoking definition, number of patients in the smoking and non-smoking hr-HPV+ groups. The 5-year OS rate was extracted from the Kaplan-Meier curve using GraphClick (Version 3.0.3, Arizona Software 2010, www.arizona-software.ch/graphclick). In case of discrepancies, findings were re-analyzed and discussed until consensus was reached.
Statistical analysis. The 5-year OS of the smoking and non-smoking groups were tested using the risk ratio (RR). The summary RR was estimated including the 95% confidence interval (CI). The heterogeneity was tested using a chi-squared-based Q test; a p-value greater than 0.05 means homogeneity. The existence of heterogeneity was evaluated by using the I2 index in the meta-analysis, which was represented as a percentage value between 0 and 100. We initially applied a fixed-effects model (Mantel-Haenszel method and chi-squared test) to the data. When there was significant heterogeneity, we used the random-effects model (DerSimonian-Liard method). We compared the RR of the 5-year OS of the two groups using a forest plot. The funnel plot examined publication bias. For all statistical analysis, we used R Version 3.5.1: A language and environment for statistical computing [R Foundation for Statistical Computing, R Core Team (2018), Vienna, Austria].
Risk of bias. Publication bias was assessed using a funnel plot. After formulating a review question in a hypothetical randomized trial (target trial), the risk of bias was analyzed for each study using the ROBINS-I tool (16). The tool includes seven domains: Confounding, selection of participants for the study (pre-intervention), classification of intervention (at intervention), deviation from intended interventions, missing data, measurements of outcome, selection of the reported results (post-intervention). The risk of bias analysis was performed independently by two Authors (T.H. and V.V.). All eligible publications were analyzed in detail. Relevant information was extracted to assess the risk of bias in each of the aforementioned domains using the classifications low, moderate, serious or critical risk of bias. Conflicting results were discussed to reach a consensus.
Results
Literature search. A total of 748 publications were identified. After multiple assessment phases only nine publications met our inclusion criteria (Figure 1), with total of 2,080 patients (mean=231), 1,663 (80%) male and 417 (20%) female. However, only HPV+ patients were included in our study, representing 1,436 patients (range per study=46-484 patients). All nine included studies were cohort studies. The main characteristics of the identified studies are summarized in Table I (11, 17-24). A total of 367 (36%) were HPV+ smokers and 656 (64%) were HPV+ non-smokers. Five studies investigated only OPSCC (17, 19, 21, 23, 24), three studies included patients with HNSCC (11, 20, 22), and one study reported patients with carcinoma of unknown primary (18). Three studies were performed in Europe (11, 17, 18), five studies in North America (20-24), and one study in Asia (19). For determination of HPV states, six studies used the hr-HPV surrogate marker p16 followed by HPV detection using HPV polymerase chain reaction (PCR) or in-situ hybridization (ISH) for DNA detection (11, 17-19, 23, 24). In three studies, HPV was detected using HPV PCR or ISH without studying the p16 status (20-22).
Meta-analysis. All nine studies investigated the 5-year OS of the two hr-HPV+ HNSCC smoking and nonsmoking subgroups (Figure 2). Ηr-HPV+ non-smoking patients had a better prognosis compared with hr-HPV+ smoking patients (no significant heterogeneity, p=0.83, fixed-effect model: RR=1.94, 95% CI=1.46-2.58).
Subgroup analysis. The sub-meta-analyses according to tumor locations included OPSCC and HNSCC other than OPSCC (OPSCC: no significant heterogeneity, p=0.83; fixed effect model, RR=1.86, 95% CI=1.36-2.54, and HNSCC without OPSCC: no significant heterogeneity, p=0.38; fixed effects model: RR=2.23; 95% CI=1.07-4.65). In both sub analyses, RR and CI were essentially unaltered compared with the major meta-analyses of the complete data set.
The sub-meta-analyses for HPV proof showed comparable results to the major meta-analyses (ISH HPV DNA PCR/ISH and p16: no significant heterogeneity, p=0.91; fixed effect model, RR=1.68; 95% CI=1.19-2.38, and HPV DNA PCR or ISH without p16: no significant heterogeneity, p=0.68, fixed effects model, RR=2.58, 95% CI=1.59-4.20).
Sensitivity analysis. The parameters of 50% of the most recent studies (no significant heterogeneity, p=0.59; fixed effects model, RR=1.97, 95% CI=1.45-2.68) (11, 21-24) confirmed the results of the whole meta-analysis with the whole data set.
Publication bias and risk of bias assessment. The funnel plot shape did not reveal any obvious evidence of asymmetry (Figure 3). Egger’s test showed no significant publication bias (p=0.68). The risk of bias was assessed using the Risk of Bias in Non-randomized Studies - of Interventions (ROBINS)-I tool, a tool designed specifically to assess the risk of bias in systematic reviews (25). The mapping of domain-level judgements is shown in Table II. Confounding domains included age, gender, co-morbidities and TNM classification. Alcohol consumption, surgery, radiotherapy and chemotherapy were defined as hypothetical co-interventions in the pre-defined target trial. Usually, no controlling method was applied for all criteria which led to the judgement of moderate to high risk of bias in any included study. The wide variety of criteria for patient inclusion in all studies regarding HPV detection method, HPV subtype and tumor sublocalization of HNSCC resulted in moderate to serious risk of bias. p16 testing followed by DNA ISH or PCR was considered the gold standard for HPV detection and all high-risk HPV subtypes. The HNSCC sublocalizations would have been included in the analysis in an ideal trial. Furthermore, there was a lack of consistent definition of the intervention group ‘smoker’. Some authors defined smoking using standardized measures such as pack-years, while others focused on the timely aspect and defined a smoker as a patient who actively/currently smokes. Serious risk of bias was determined in the ‘deviation of intended intervention’ domain if there was evidence of imbalance of co-interventions among the sub-groups, e.g. higher alcohol consumption in connection with smoking. Moreover, missing data can introduce bias if the same HPV detection methods are not applied to all patients. The measurement of the OS outcome was judged as low-risk of bias in all studies and the selection of reported results carried generally low to moderate risk of bias.
Discussion
Many studies have investigated HNSCC patient outcomes according to HPV or smoking status. It is known that hr-HPV+ cancers tend to have better outcomes than HPV-negative ones and that smokers usually have a worse prognosis than non-smokers (8-12). Recently Chen et al. published a systematic review showing that all included studies reported that smoking was associated with worse outcomes in hr-HPV+ cancer (26). Variability of the reported studies regarding reported outcomes and smoking metrics prohibited further analyses. Yet to our knowledge, there has never been a meta-analysis of data and results regarding prognosis of patients with hr-HPV+ HNSCC based on smoking status. Our meta-analysis compares the outcome of smokers and non-smokers with hrHPV+ HNSCC.
The analysis showed that smoking plays an important prognostic role in HNSCC, even in the presence of HPV. Our results suggest that non-smokers with hr-HPV+ HNSCC have a better 5-year OS compared to their smoking counterparts. Others showed an increased hazard ratio for hr-HPV+ smokers (27-29). Further studies examined the disease-specific and disease-free survival (30-34). Going back to the pathogenesis of the HNSCC, we can find two different pathways for its carcinogenesis. In HPV-associated HNSCC, infection with hr-HPV-types leads to an increased expression of E6 and E7, which degrade and inactivate p53 and enhance expression of retinoblastoma protein, which causes overexpression of p16 (35, 36). On the other hand, HPV− HNSCCs are associated with smoking and alcohol, and exposure to carcinogenic substances such as tobacco-specific nitrosamines that lead to p53 mutation (37, 38). However, in hr-HPV+ carcinomas, smoking is an independent risk factor for reduced treatment efficiency and OS (39-41). The significant negative impact of smoking on prognosis might be associated with increased smoking-related comorbidities, such as pulmonary and cardiovascular diseases, therefore contributing to a worse survival chance after treatment and increased mortality. Furthermore, studies have shown that smoking increases tumor aggressiveness in cancer cells by stimulation proliferation, angiogenesis, migration, invasion and reducing the response to cytotoxic cancer agents such as chemo- or radiotherapy (42, 43). House et al. recently showed that miR-133a-3p is a target of smoking-induced changes in patients with HPV+ OPSCC and alters the expression (44). Therefore, recording the smoking status should be an important part of staging and treatment even in patients with hr-HPV+ HNSCC to predict individual patient risk (45). This also highlights the importance of smoking cessation programs for patients who are actively undergoing treatment as a smoking habit is a potentially modifiable adverse risk factor that might lead to improved outcomes possibly even after diagnosis. The Eighth Edition of the American Joint Committee on Cancer TNM classification distinguishes p16-associated OPSCC from p16− OPSCC (9). The negative impact of smoking is not taken into account in the present classification (39).
HPV detection methods also play an important role. p16 is used as a surrogate marker for the detection of hr-HPV in HNSCC, and p16 immunohistochemistry (IHC) has a sensitivity of 94% and specificity of 82-83%. HPV DNA PCR has a sensitivity of 97-98% and specificity of 84-87%, whilst HPV DNA ISH has a sensitivity of 85-88% and specificity of 83-88% (46, 47). Therefore, the combination of p16 IHC followed by HPV PCR or ISH remains a gold standard for HPV detection in HNSCC to confirm the HPV-related etiology of HNSCC because, as previously demonstrated, p16+ HNSCC with HPV-independent carcinogenesis may exist at a low percentage (10). In the present sub-meta-analysis, six out of nine studies used HPV DNA PCR or ISH to confirm the etiology of HPV-associated HNSCC in patients with p16+ IHC (Table I). Even the sub-meta-analysis including the three studies using only HPV DNA PCR or ISH confirmed the results shown using the whole data set. However, this detection strategy does not confirm p16-associated carcinogenesis and there is the possibility that HPV DNA may also have been an innocent bystander (in cases with p16− HNSCC) (10).
There are some limitations to our study: (a) This meta-analysis is based on data extracted from published literature and not from individual patient data. (b) All included studies were found to carry moderate to serious risk of bias. The vast majority of publications were retrospective cohort studies/non-randomized, controlled trials. These types of studies have greater susceptibility to bias. An important factor to be considered is the role of confounding factors and co-interventions. One problem was the different definitions of smoking (risk of bias analysis). Apart from smoking, there are other factors which also play an important role in prognoses and outcomes of patients with HNSCC such as comorbidities, tumor stage, therapy regimes (e.g. surgery, chemotherapy or radiation). Some studies included different subsites of HNSCC (11, 20, 22), and others focused on OPSCC only (17, 19, 21, 23, 24).
This heterogeneity in patient selection should be noted because in patients with OPSCC, the impact of HPV in tumorigenesis is accepted. In other sub-localizations, there is no difference for HNSCC with confirmed p16/HPV detection, and HPV detection is therefore not required.
Conclusion
Smoking plays a significant role in the prognosis of patients with hr-HPV+ HNSCC. Smokers with hr-HPV+ HNSCC have worse survival outcomes compared to hr-HPV+ nonsmokers. Smoking status should be an important part of staging and be considered when tailoring the treatment strategy even for patients with HPV-associated HNSCC, since these patients are at possible risk of undertreatment.
Acknowledgements
The Authors would like to thank Ulrich Gauger from Private Statistical Office, Berlin, Germany.
Footnotes
Authors’ Contributions
Moonef Alotaibi conducted the systematic review, collected and analyzed the data, wrote the main article. Valeria Valova collected and analyzed the data, conducted the risk of bias analysis, wrote the main article. Toni Hänsel conducted and supervised the risk of bias analysis, revised the main article. Carmen Stromberger revised the main article. Grzegorz Kofla revised the main article. Heidi Olze revised the main paper. Iris Piwonski revised the main article. Andreas Albers revised the main article. Sebastian Ochsenreither revised the main article. Annekatrin Coordes designed the review, conducted the systematic review and supervised the work. All Authors commented on the article at all stages.
This article is freely accessible online.
Conflicts of Interest
The Authors declare that they have no conflicts of interest in regard to this study.
- Received October 4, 2020.
- Revision received November 27, 2020.
- Accepted December 4, 2020.
- Copyright© 2021, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved