Understanding the impact of high-risk human papillomavirus on oropharyngeal squamous cell carcinomas in Taiwan: A retrospective cohort study
============================================================================================================================================

* Guadalupe Lorenzatti Hiles
* Kai-Ping Chang
* Emily L. Bellile
* Chun-I Wang
* Wei-Chen Yen
* Christine M. Goudsmit
* Hannah L. Briggs
* Trey B. Thomas
* Lila Peters
* Macy A. Afsari
* Lisa M. Pinatti
* Anna C. Morris
* Nadine Jawad
* Thomas E. Carey
* Heather M. Walline

## Abstract

**Background and Objectives** Human papillomavirus (HPV)-driven oropharyngeal squamous cell carcinoma (OPSCC) is increasing globally. In Taiwan, HPV-positive OPSCC is obscured by tobacco, alcohol, and betel quid use. We investigated the role of high-risk HPV (hrHPV) in a large retrospective Taiwan OPSCC cohort.

**Methods and Results** The cohort of 541 OPSCCs treated at Chang Gung Memorial Hospital from 1998-2016 consisted of 507 men (94%) and 34 women (6%). Most used tobacco (81%), alcohol (51%), and betel quid (65%). Formalin-fixed, paraffin-embedded tissue was used for p16 staining (a surrogate marker for HPV) and testing for HPV DNA presence and type by Multiplex HPV PCR-MassArray. HPV DNA and/or p16 staining (HPV-positive) was found in 28.4% (150/528) tumors. p16 and HPV DNA were strongly correlated (*F* < 0.0001). HPV16 was present in 82.8%, and HPV58 in 7.5% of HPV-positive tumors. HPV was associated with higher age (55.5 vs. 52.7 years, *p* = 0.004), lower T-stage (*p* = 0.008) better overall survival (OS) (hazard ratio [HR] 0.58 [95% CI 0.42-0.81], *p* = 0.001), and disease-free survival (DFS) (HR 0.54 [95% CI 0.40-0.73], *p* < 0.0001). Alcohol was strongly associated with recurrence and death (OS: HR 2.06 [95% CI 1.54-2.74], *p* < 0.0001; DFS: HR 1.72 [95% CI 1.33-2.24], *p* < 0.0001). OS and DFS in HPV-positive cases decreased for alcohol users (*p* < 0.0001). Obscured by the strong alcohol effect, predictive associations were not found for tobacco or betel quid.

**Conclusions** As with HPV-positive OPSCC globally, HPV is an increasingly important etiological factor in Taiwanese OPSCC. HPV-positive OPSCC has considerable survival benefit, but that is reduced by alcohol, tobacco, and betel quid use. hrHPV is a cancer risk factor in males and females. Vaccinating both sexes with a multivalent vaccine including HPV58, combined with alcohol and tobacco cessation policies will be effective cancer-prevention public health strategies in Taiwan.

## Introduction

The occurrence of oropharyngeal squamous cell carcinoma (OPSCC) is rapidly increasing in North America and Western Europe, accounting for approximately 100,000 new cases worldwide each year [1-3]. In particular, the incidence of OPSCC has been dramatically rising since 1973, at the point of surpassing 5% annual increment in the United States in 2000 [2, 4, 5]. OPSCC has been traditionally associated with tobacco use and excessive alcohol consumption as primary risk factors [6-17]. However, recent behavioral changes in Western countries have promoted a marked drop in the prevalence of these major risk factors [7, 13, 18, 19]. In contrast, high-risk human papillomavirus (hrHPV), HPV genotypes 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, 68, and 73, has become the leading etiologic factor of OPSCC [2, 20-28].

Since the World Health Organization recognized the causative link between 15 hrHPV genotypes and the occurrence of OPSCC in 2007 [29], hrHPV has been accepted as a principal etiological cause of this cancer [2, 20, 22, 24, 30-34]. HPV-driven OPSCC is markedly on the rise [2, 4, 10, 22-24, 35-37]. In the United States, the estimated proportion of positive cases has increased from 20% in 1990 to >70%, where hrHPV now represents the most common cause of OPSCC [2, 7, 27, 36, 38]. Countries in Western Europe have observed similar trends [7, 13, 36, 37, 39-41]. Interestingly, these changes have been accompanied by an increment in the survival rates for OPSCC [36, 42, 43]. HPV-positive patients have a significantly better response to treatment (radiation therapy and chemotherapy as well as surgery) and a more favorable prognosis than those diagnosed with HPV-negative OPSCC [23, 25, 27, 28, 35, 42, 44-50].

Despite these observations, recent studies from Taiwan suggest that its high OPSCC rates continue to increase predominantly due to heavy alcohol drinking, cigarette smoking, and betel quid chewing as etiologic factors [9, 51-54]. The strong influence of these risk habits has limited the search for a viral etiology in this population. There have been a few studies indicating that hrHPV is an emerging risk for head and neck cancer in South-East Asia, with a prevalence of HPV-positive OPSCC reported to be absent or present in 12.6% [55-57] to 34% [57-62] of OPSCCs. In this study, we conducted a retrospective cohort analysis to interrogate the prevalence and significance of HPV-driven OPSCC in tissue samples collected from a single major referral site in Taiwan over a period of 18 years. We evaluated the association between clinical characteristics and traditional risk factors (alcohol, smoking, and betel quid) with HPV-associated OPSCC in Taiwan. HPV results were correlated with risk factor exposure for outcomes and survival analysis.

## Methods

### Case identification and study design

This study was performed on a retrospective cohort of OPSCC cases diagnosed from March 1998 to February 2016 at the Chang Gung Memorial Hospital (CGMH)-Linkou in Taoyuan, Taiwan (Taiwan cohort), as described in Figure 1. CGMH is the largest cancer center and a major referral center in Taiwan. Case selection was not a source of bias as we identified OPSCC tumors with a confirmed primary site in the oropharynx using hospital electronic and pathology records from all patients that underwent curative-intent therapy at this tertiary healthcare center. Patients with unknown primary site (T0 or Tx) were excluded. Primary tumors were biopsied only or resected by surgery and collected for histopathological diagnosis. All patients with available pathology-archived tissue and known clinical records were included in this cohort (Fig 1). In total, 541 OPSCC tumors were retrieved, sectioned, anonymized, and shipped to the University of Michigan for HPV and p16 testing. Five non-squamous cell carcinoma cases were excluded from the study based on their pathological classification and histopathological re-assessment of the submitted sections.

![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/12/2020.10.26.20213231/F1.medium.gif)

[Fig 1.](http://medrxiv.org/content/early/2021/03/12/2020.10.26.20213231/F1)

Fig 1. Flow diagram of cases included in the Taiwan retrospective cohort and study design.
1Results missing due to absent slide or tissue core, or major artifacts that prevented evaluation. 2Results missing due to absent DNA or invalid test. OPSCC, oropharyngeal squamous cell carcinoma; HPV, human papillomavirus; FFPE, formalin-fixed, paraffin-embedded; IHC, immunohistochemical staining; PCR-MA, multiplex PCR-MassArray.

To analyze the contribution of hrHPV on these OPSCC cases, qualitative data on the common risk factors, alcohol consumption, smoking, and betel quid chewing habits were collected. However, because of the retrospective nature of the study, and a change to an electronic record system, data on the quantity of alcohol, smoking, betel quid consumption, or comorbidities could not be retrieved for a large proportion of the patients. Smoking is the only use of tobacco in Taiwan because betel quid preparations do not contain tobacco and tobacco chewing is an extremely uncommon behavior [9, 63]. Demographic information, including patient characteristics (age at diagnosis, and gender), as well as clinical information (stage, tumor site, initial treatment, and outcomes for recurrence, metastasis, and death), were compiled from patient records. Cases were staged at diagnosis according to the seventh edition of the American Joint Cancer Committee (AJCC) [64].

### HPV interrogation

All OPSCC tumors were evaluated for the presence of HPV by two complementary methods: p16 testing and detection of HPV DNA types [65-67] (Fig 1). Results from each determination were blinded to the investigators to avoid bias. Tumors with either p16 and/or HPV DNA positivity were defined as HPV-positive.

#### Detection and genotyping of HPV DNA by Multiplex PCR-MassArray (PCR-MA)

DNA was isolated from tissue curls of formalin-fixed, paraffin-embedded (FFPE) tumor specimens. Two to seven 10-µm FFPE sections were combined in each extraction with AllPrep DNA/RNA FFPE Kit (Catalog No. 80234, QIAGEN, Germantown, MD, US), according to the manufacturer’s recommendations. DNA was eluted in DNase-free water and stored at −20 °C until testing. DNA concentration was measured by a Qubit 2.0 Fluorometer (Catalog No. Q32866, Invitrogen-Thermo Fisher Scientific Inc., Waltham, MA, US) and the Qubit dsDNA HS Assay Kit (Catalog No. Q32851, Invitrogen-Thermo Fisher Scientific Inc., US).

Samples were examined for the presence of HPV DNA and genotyped by PCR-MA analysis, a very sensitive, high-throughput method based on competitive PCR and probe-specific single-base extension coupled with MALDI-TOF mass spectrometry [65-67]. The PCR-MA assay is designed to detect 15 high-risk (HPV 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, 68, and 73), and 2 low-risk (HPV 6 and 11) HPV types, and a possible high-risk subtype (HPV90), as previously described by our laboratory [65-67]. Reactions were prepared with 20 ng of DNA and carried out in quadruplicates in an area physically separated from DNA isolation. Tests were run using a Mass Array 384-format System (Agena Bioscience Inc., San Diego, CA, US). Specimen acceptability was determined using human glyceraldehyde-3-phosphate dehydrogenase (GAPDH) DNA control.

#### p16 testing by immunohistological analysis

FFPE tissue sections (4-µm) were used for p16 immunostaining with a specific antibody against Protein Cyclin-dependent kinase inhibitor 2A (CDKN2A), also known as p16INK4a, as a surrogate for transcriptionally and translationally active HPV [68]. The immunohistochemical (IHC) staining was carried out manually using the clinically validated CINtec-p16 (E6H4) antibody (pre-diluted, Ref. No. 725-4713, Ventana Medical Systems Inc., Tucson, AZ, USA) as stated by the supplier’s protocol. The CINtec-p16 primary antibody was incubated for 1 hour at room temperature followed by washing and appropriate horseradish peroxidase-labeled secondary antibody for 30 minutes at room temperature. All slides were stained with 3,3′-diaminobenzidine for 1-5 minutes, followed by hematoxylin counterstain.

p16 IHC was examined for each slide at 200x and 400x magnification according to the 2018 recommendations of the College of American Pathologists [68]. p16 expression was scored as positive if ≥ 70% of the tumor cells exhibited strong and diffuse nuclear and cytoplasmic p16 immunoreactivity (S1 Fig).

### Statistical analysis

Two-sided Fisher’s exact test (*F*) was used to analyze the relationship between p16 and HPV DNA results. The association between HPV prevalence and study year was evaluated by two-sided simple linear regression and Spearman rank correlation (ρ). Data for HPV status, alcohol consumption, cigarette smoking, betel quid chewing, age, gender, T-stage, N-stage, disease site, initial treatment, and clinical outcomes were collected as covariates. Standard descriptive statistics were performed for each covariate collected. Differences in the distribution of covariates by HPV status were tested by two-sided t-test (continuous measures) or Pearson’s chi-squared test (χ2) (categorical/binary). Time-to-event outcomes were defined beginning from date of pathology diagnosis to death from any cause (Overall Survival), or from date of pathology diagnosis to date of first recurrence or death (Disease-Free Survival); subjects alive with no event were censored at date of last follow-up. The Kaplan-Meier method and log-rank tests were used to estimate survival probabilities and plot survival distributions. Cox proportional hazard models and hazard ratio (HR) estimations for time up to 5 years post-diagnosis were performed to test relative hazards between groups in the whole cohort and in subsets stratified by HPV status or other risk factors (alcohol consumption, cigarette smoking, and/or betel quid chewing), adjusting for age, T- and N-stage. Cases with no HPV status data (p16 or HPV DNA), N = 13, were excluded, leaving 528 cases eligible for Chi-squared and survival calculations. Tests were also performed to examine the variates and survival distributions by gender. Statistical analyses were conducted in SAS v9.4 (SAS Institute Inc., Cary, NC, US) using R v3.6.1 (RStudio, Boston, MA, US) for graph generation, or GraphPad Prism v8.3.0 (GraphPad Software, San Diego, CA, US). Statistical tests were performed using 95% confidence intervals and a 5% significance level.

### Ethics statement

This retrospective study was approved by the Institutional Review Boards of the University of Michigan Medical School and the Chang Gung Memorial Hospital and conducted in compliance with the ethical guidelines of the World Medical Association’s Declaration of Helsinki (1964, amended in 2013) and local regulations. Additional patient consent was not required by the institutional review boards as this OPSCC cohort comprised secondary use of tissue specimens with unidentified chart data. All information stripped of personal identifiers to ensure that the data cannot be linked to individual cases in this cohort, are available in the supplementary S1 Table. The procedures described in this manuscript followed the reporting standards for human subject research of the EQUATOR Network, which are detailed in the STROBE report for this study (S1 Checklist).

## Results

### HPV status and clinical characteristics

A total of 546 OPSCC cases were obtained from an unbiased retrospective chart review of individuals treated with standard of care therapy from March 1998 to February 2016 at the CGMH in Taiwan (Taiwan cohort). Among these cases, five were not OPSCC according to the pathology records and slide review; these were excluded as they did not fulfill our inclusion criteria. Therefore, the final study cohort included 541 OPSCC cases (Fig 1, S1 Table).

The presence of HPV in FFPE tumor sections of oropharyngeal cancer was assessed by IHC staining for p16, a surrogate marker for HPV [68] (Fig 1, S1 Fig, S1 Table). HPV genotypes were identified by PCR-MA [65-67] using tumor genomic DNA (Fig 1, S1 Table). Of the 541 OPSCC tumors tested, p16 was positive in 115 (21.3%), negative in 355 (65.6%), and 71 (13.1%) could not be scored. HPV detection and genotyping showed that 134 (24.8%) tumors were HPV-positive, 379 (68.4%) HPV-negative, and 37 (6.8%) had insufficient DNA (failed to amplify the GAPDH control) (Fig 1, S1 and S2 Tables). Of the 134 positives, HPV16 was found alone in 103 (76.9%) and HPV58 was the second most frequently found in 10 (7.5%) tumors. HPV66, HPV59, HPV45, HPV39, HPV34, HPV31, HPV18, and HPV6 were also found in 13 (9.7%) individual tumors. HPV16 was also present together with other HPV genotypes (HPV6, HPV18, HPV35, HPV58, or HPV59) in 8 (6.0%) tumors (Table 1, S1 Table). These results indicate, that aside from low-risk HPV6, there are diverse oncogenic HPV genotypes in Taiwan OPSCCs.

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/03/12/2020.10.26.20213231/T1)

Table 1. HPV genotypes frequency.

Altogether, p16 overexpression was strongly correlated with HPV status in OPSCC, as the concordance between p16 and HPV DNA testing was 94.9% (423 out of 446, *F* < 0.0001, S2 Table). Therefore, for this study, we defined HPV positivity as either positive by p16 or HPV DNA test. Thus, we had 528 OPSCC tumors with HPV status (positive or negative) called. HPV status (p16 and HPV DNA) data were not obtained for thirteen cases and were not included in the analysis (Fig 1, S1 and S2 Tables). The prevalence of HPV-positive OPSCC in the whole cohort was 28.4% (150 out of 528) (Fig 1). Interestingly, when we examined the yearly occurrence of HPV-positive OPSCC, we found that there was a trend for an increase over the 18 years of study, but it failed to reach statistical significance (Fig 2, S3 Table). However, the same trend is significant when we examined the yearly occurrence of p16 alone (Fig. 2, S3 Table). Nonetheless, this result should be carefully interpreted as 71 cases are missing data for p16 (S1-S3 Tables). We also observed a clear increment in the number of HPV-negative cases that drive the growing incidence of OPSCC rates in Taiwan, thereby obscuring the gradual rise of HPV-positive OPSCCs. Nevertheless, our results demonstrate an increasing role of oncogenic HPV and its causal role as an etiologic factor of OPSCC in this population.

![Fig. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/12/2020.10.26.20213231/F2.medium.gif)

[Fig. 2.](http://medrxiv.org/content/early/2021/03/12/2020.10.26.20213231/F2)

Fig. 2. Yearly HPV occurrence among OPSCC cases by HPV DNA and/or p16 (A) or p16 alone (B).
The graphs show the correlation between the total frequency of HPV-positive (HPV+) and HPV-negative (HPV-) OPSCC cases, and the study years (see S3 Table). HPV status was assessed by (A) HPV DNA and p16 testing (N = 528) or (B) p16 scoring (N = 458). The association was evaluated in the Taiwan cohort from March 1998 to February 2016 by Spearman’s coefficient (ρ) and linear regression (R2). (A) HPV-: ρ = 0.6953, *p* = 0.0014; R2 = 0.5201, *p* = 0.0007. HPV+ ρ = 0.4093, *p* = 0.0917; R2 = 0.1952, *p* = 0.0664. (B) p16-: ρ = 0.6991, *p* = 0.0012; R2 = 0.5555, *p* = 0.0004. p16+ ρ = 0.4741, *p* = 0.0469; R2 = 0.2455, *p* = 0.0365.

Next, we assessed the demographic and clinical features of the whole cohort, as listed in Table 2. In our study, males represented 94% (507 out of 541) of all cases, and the average age at tumor diagnosis was 53.5 years old, with 66% of the tumors, diagnosed in individuals of 41-60 years of age (358 out of 541). Over half of the tumors, 58% (315 out of 541), were biopsied or resected from the tonsils, followed by 24% (132 out of 541) from the soft palate, 17% (90 out of 541) from the base of the tongue, and 1% (4 out of 541) from other non-specified locations in the oropharynx. Unknown primary tumors diagnosed by neck node pathology were not included in the tumor retrieval, which may account for the relatively low incidence of base of tongue tumors in this cohort. The large majority of the cohort had a history of previous exposure to known risk factors: 51% (278 out of 541) drank alcohol, 83% (448 out of 541) smoked, and 65% (349 out of 541) chewed betel quid (Table 2, S4 Table). Most prominently, 87% (468 out of 541) of these individuals used more than one of alcohol, and/or tobacco, and/or betel quid concomitantly. Never-smokers, never-drinkers, and never-betel quid chewers accounted for a small 13% (73 out 541) of the cases (S4 Table). Because of the high density of risk factors, and their combined exposure, we were unable to separate their individual effects. Therefore, only 1% (7 out of 541) of the Taiwan cohort was exposed to alcohol without smoking or betel quid chewing, 11% (60 out of 541) solely smoked, and 1% (5 out of 541) only consumed betel quid (S4 Table).

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/03/12/2020.10.26.20213231/T2)

Table 2. Demographic and clinical characteristics.

We then examined the clinical determinants of HPV status (Table 2). When compared to HPV-negative cases, HPV-positive tumors (HPV DNA-positive and/or p16-positive) had a slightly higher average age at diagnosis (55.5 vs. 52.7 years, *p* = 0.004), with most individuals presenting between 41-60 years old in both groups (67% vs. 63%, *p* = 0.18), and slightly lower T-stage (majority of T2 cases vs. T4, *p* = 0.008). The HPV-positive tumors were also most frequently located in the tonsils (73% vs. 52%, *p* < 0.0001); and presented with lower exposure to all the risk factors including alcohol (37% vs. 56%, *p* < 0.0001), smoking (65% vs. 90%, *p* < 0.0001), and betel quid (35% vs. 75%, *p* < 0.0001) (Table 2, S4 Table). Although the number of females was far lower, representing only 6% (34 out of 541) of the cases (Table 2), they showed pronounced differences with males (S5 Table). The majority, 62% (21 out of 34) of tumors from females, but a minority, 25% (129 out of 507) of tumors from males, were HPV-positive (*p* < 0.0001). The proportion of tonsil tumors was also higher in females than in males (85% vs. 56%, *p* = 0.01), and were more likely to be N3 (*p* = 0.03) (AJCC 7th edition). The females had lower exposure to alcohol (18% vs. 54%, *p* < 0.0001), smoking (21% vs. 87%, *p* < 0.0001), and betel quid (12% vs. 68%, *p* < 0.0001); and were more likely to be never users of these high risk carcinogens. Females also tended to be older at diagnosis (mean: 58.7 years vs. 53.1 years, *p* = 0.002). Among both males and females the majority of cases were of 41-60 years of age (64% v. 66%, *p* = 0.13) (S5 Table). There was also a slightly better prognosis for the female group (OS: log-rank *p* = 0.05; DFS: log-rank *p* = 0.07) (S2 Fig), but given the small group size (N = 34), this observation should be carefully interpreted.

### The role of HPV status and associated risk factors on OPSCC outcome

First, we determined if HPV status had survival benefits on OPSCC by multivariable and Kaplan-Meier analysis (Fig 3, S6 Table). Compared to HPV-negative cases, patients with HPV-positive OPSCC had significantly higher overall survival (HR 0.58, 95% CI 0.42 to 0.81, *p* = 0.001; log-rank *p* < 0.0001) and disease-free survival (HR 0.54, 95% CI 0.40 to 0.73, *p* < 0.0001; log-rank *p* < 0.0001) for up to 5-year post-diagnosis, suggesting that HPV-positivity was an independent predictor for better prognosis.

![Fig 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/12/2020.10.26.20213231/F3.medium.gif)

[Fig 3.](http://medrxiv.org/content/early/2021/03/12/2020.10.26.20213231/F3)

Fig 3. HPV-positive OPSCC is associated with increased survival time.
(A-B) Up to 5-year overall survival (OS) and disease-free survival (DFS) prognostic outcomes of the HPV variable in the whole OPSCC Taiwan cohort. HPV positivity is defined as HPV DNA-positive and/or p16-positive. (A) Table includes the multivariable hazard probabilities analyzed using Cox survival models and hazard ratio (HR) estimations, visualized by forest plots. The complete analysis is found in S6 Table, where estimates were reported for full model with all covariates (HPV status, alcohol, smoking, betel quid, age, N- and T-stage) included as fixed effects. (B) Kaplan-Meier survival analysis. Plots represent the results for up to 5-year OS (left) and DFS (right) comparison between HPV-negative (HPV-) and HPV-positive (HPV+) groups. Log-rank analysis was used to compare the survival distributions (log-rank *p*-values are in the plots). HPV-, HPV-negative; HPV+, HPV-positive.

Similar analyses also revealed that consumption of alcohol was a strong negative prognostic factor for both up to 5-year overall survival (OS) (HR 2.06, 95% CI 1.54 to 2.74, *p* < 0.0001; log-rank *p* < 0.0001) and disease-free survival (DFS) (HR 1.72, 95% CI 1.33 to 2.24, *p* < 0.0001; log-rank *p* < 0.0001) (Fig 4, S6 Table). Surprisingly, smoking, and betel quid chewing had no predictive effects (Fig 4, S6 Table). Smoking did not have statistical significance for worse overall survival (HR 0.76, 95% CI 0.50 to 1.14, *p* = 0.18; log-rank *p* = 0.24) and disease-free survival (HR 0.81, 95% CI 0.56 to 1.17, *p* = 0.26; log-rank *p* = 0.40). Betel quid exposure showed no effect on survival by multivariable analysis (OS: HR 0.92, 95% CI 0.67 to 1.27, *p* = 0.60; log-rank *p* < 0.05 – DFS: HR 0.90, 95% CI 0.66 to 1.21, *p* = 0.46; but by Kaplan-Meier analysis betel quid use, just reached significance log-rank *p* < 0.05).

![Fig 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/12/2020.10.26.20213231/F4.medium.gif)

[Fig 4.](http://medrxiv.org/content/early/2021/03/12/2020.10.26.20213231/F4)

Fig 4. Alcohol is associated with reduced OPSCC survival time.
(A-D) Prognostic outcomes of the alcohol, smoking, and betel quid variables within the whole cohort. We analyzed up to 5-year overall survival (OS) and disease-free survival (DFS) outcomes for high-risk habits. (A) Table includes the multivariable hazard probabilities analyzed using Cox survival models and hazard ratio (HR) estimations, visualized by forest plots, where estimates were reported for full model with all covariates (HPV status, alcohol, smoking, betel quid, age, N- and T-stage) included as fixed effects. The complete analysis is found in S6 Table. (B-D) Kaplan-Meier survival analysis. Plots represent the results for up to 5-year OS (left) and DFS (right) comparison between alcohol (B), smoke (C), and betel quid (D) groups. Log-rank analysis was used to compare the survival distributions (log-rank *p*-values are in the plots).

Additional multivariable survival analysis and Kaplan-Meier estimations for alcohol, smoking, and betel quid groups, controlled for HPV status, were also performed (Fig 5, S3 Fig, S7 Table). Alcohol use had an adverse effect on outcome in both HPV-positive and HPV-negative groups, but HPV was associated with longer OPSCC survival time (survival by alcohol within HPV groups, OS: HPV-positive log-rank *p* = 0.0007, HPV-negative long rank *p* < 0.0001 – DFS: HPV-positive log-rank *p* = 0.02, HPV-negative long-rank *p* = 0.0005). No predictive associations were found for tobacco (survival by smoking within HPV groups, OS: HPV-positive log-rank *p* = 0.44, HPV-negative long rank *p* < 0.28 – DFS: HPV-positive log-rank *p* = 0.49, HPV-negative long-rank *p* = 0.18) or betel quid chewing (survival by betel quid within HPV groups, OS: HPV-positive log-rank *p* = 0.41, HPV-negative log rank *p* < 0.68 – DFS: HPV-positive log-rank *p* = 0.21, HPV-negative long-rank *p* = 0.46) based on HPV status. Most importantly, our data demonstrate that the prognostic benefit of HPV positivity persists in the presence of risk factors, alcohol, and tobacco, and betel quid (DFS by HPV within risk groups: alcohol yes log-rank *p* = 0.008, alcohol no log-rank *p* = 0.007; smoking and/or betel quid yes log-rank *p* = 0.0006, smoking and/or betel quid no log rank *p* = 0.01) (S4 Fig, S8 Table). Furthermore, among non-drinkers, non-smokers, and non-betel quid chewers, HPV-positive OPSCC had the best outcomes. In sum, our study shows that hrHPV has a causal role but also carries a significant prognostic benefit for OPSCC in this Taiwan cohort.

![Fig 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/12/2020.10.26.20213231/F5.medium.gif)

[Fig 5.](http://medrxiv.org/content/early/2021/03/12/2020.10.26.20213231/F5)

Fig 5. In alcohol users, HPV is associated with improved OPSCC survival time.
(A-C) Prognostic outcomes of the alcohol, smoking, and betel quid variables within HPV risk groups. HPV positivity is defined as HPV DNA-positive and/or p16-positive. We analyzed up to 5-year overall survival (OS) and disease-free survival (DFS) outcomes. (A) Table includes the multivariable hazard probabilities analyzed using Cox survival models and hazard ratio (HR) estimations, adjusted for age, T- and N-stage, and visualized by forest plots. (B-C) Kaplan-Meier survival analysis. Plots represent the results for up to 5-year OS (B) and DFS (C) comparison between HPV groups stratified by alcohol groups. Log-rank analysis was used to compare the survival distributions (log-rank p-values are in the plots). The complete analysis is found in S7 Table and S3 Fig. HPV-, HPV-negative; HPV+, HPV-positive.

## Discussion

In Western countries, the prevalence of HPV-driven OPSCC has been rising drastically in the last 4 decades to become one of most common head and neck cancers [2, 4, 10, 22-24, 35-37]. Those affected by OPSCC suffer great losses due to aggressive treatment, morbidity, and death [20, 28, 35, 36, 42, 43, 45, 47-50, 69-70]. Still, the extent of this disease and its public health impact are not well understood outside North America and Western Europe. In this retrospective cohort study, we performed a comprehensive investigation of the impact of HPV-driven OPSCC in a cohort from the largest cancer treatment center in Taiwan from 1998 to 2016. We found that HPV was present in 28.4% of the tumors, with a trend for incremental occurrence over time. HPV16 was the most prevalent genotype (82.8%), followed by HPV58 (7.5%), and other diverse genotypes. HPV-positive OPSCCs occurred in higher proportion in females and presented with different clinical features than their HPV-negative counterparts, including reduced engagement in risk behaviors such as alcohol drinking, cigarette smoking, and betel quid chewing. Additional outcome analysis of the entire cohort showed that HPV-positivity was associated with a notably higher survival rate. Surprisingly, only alcohol, but not smoking or betel quid, were strongly associated with a worse prognosis. The strong prognostic benefit of HPV remained present but reduced in the presence of the associated risk factors alcohol, smoking, and betel quid.

Many studies have demonstrated an increment of HPV-driven OPSCC in numerous countries, showing considerable geographical variability in the proportion of cases over time [36, 71, 72]. Recent estimates calculate the worldwide prevalence of HPV-positive OPSCC between 18% and 35.6% [24, 73]. Reports from the United States, demonstrate a marked incremental variability over the years, from 20% by 1990 to over 70% of OPSCCs being currently caused by hrHPV [2, 7, 27, 36, 38]. Studies at the University of Michigan alone have shown a prevalence of 82.3% [27]. Reports from other developed countries and Western Europe, such as the United Kingdom and Finland, have observed a similar increment in prevalence of HPV-positive OPSCC during the last decades [7, 13, 36, 37, 39-41-72]. Differently, in South-East Asia, recent studies have shown a slower trend for these increments, with the proportion of HPV-positive OPSCC varying from 0% to 34% in various populations [9, 20, 57-62]. In neighboring Hong Kong, hrHPV has been found in 20.8% of tumors [59]. Likewise, a limited number of studies from Taiwan, suggested that hrHPV is a rising etiological factor in head and neck cancer, including OPSCC [55-57]. Particularly, Chien et al. demonstrated that in the early 2000s, 12.6% of the squamous cell tonsillar carcinomas were positive for hrHPV in Taiwan [57]. Here, our study indicates that the proportion of HPV-positive OPSCCs in the cosmopolitan Taiwanese population presented an incremental trend from 1998 to 2016. With an overall prevalence of 28.4%, our results suggest similar proportions in Taiwan to those observed in 1990 in the US [38]. Given this global historical data, we anticipate that Taiwan will also have a significant marked increment going forward in the rates of HPV-related OPSCC. Although the observed prevalence is three-fold lower in our Taiwan cohort compared with Western countries, this is higher than that reported for the US Asian population (12.8%) [46]. Our results also showed high concordance between HPV-DNA and p16 surrogate marker positivity (94.9%, 423 out of 446, *F* < 0.0001, S2 Table). p16 IHC is a robust surrogate marker and predictor for HPV-caused OPSCC, with very low percentage of false negatives (4%) [68]. These, findings indicate that HPV was present in the tested specimens and a likely etiological driver of oropharyngeal carcinogenesis, suggesting that HPV is not likely a passenger virus in most of these tumors. Still, several cases had discordant results: 2.2% (12 out of 541) tumors tested p16 negative but had HPV DNA positive results, while 2.0% (11 out of 541) tumors were p16 positive but had HPV DNA negative results (S2 Table). There are potential explanations for to account for the discordant cases. The p16-negative/HPV DNA positive cases may represent tumors where HPV is an incidental inactive passenger and not a causal driver of disease, or tumors for which p16 has been inactivated by another mechanism. Possible explanations for the p16-positive/HPV DNA negative cases could be that HPV DNA is in fact present in the tumor, but is mutated in the region where amplification/detection occurs for the test, or that p16 is upregulated by a different pathway or mechanism (something other than HPV).

Of important note, our study sheds light on the causal association between HPV genotypes and OPSCC in Taiwan. It is well-recognized that certain high-risk viral genotypes are carcinogenic and highly associated with the development of OPSCC [2, 36], especially HPV16 [3, 11, 29, 42, 73, 74]. In line with the predominant worldwide prevalence of HPV16 in OPSCC, this has also been identified in the majority of cases from Southern China [59], and previous reports in tonsillar squamous cell carcinomas in Taiwan [57]. Coincidentally, in our current study, HPV16 accounted for the vast majority (82.8%) of HPV-positive OPSCCs, followed by HPV58 in 7.5% of cases. Apart from HPV16 and HPV58, other oncogenic hrHPV genotypes were infrequent, and HPV18, which often causes cancer in the Western world, was only present in a very small proportion of the specimens (1.5% alone and 1.5% in combination with HPV16). Interestingly, our previous work has shown similar low proportion of HPV18-positive OPSCCs in the United States [65]. In addition, it has been reported that oral HPV infections in smokers and betel nut chewers, are mainly caused by HPV16. Still, HPV58 has been found at a low percentage in Northern Taiwan, especially in patients with the highest exposure to traditional habits that increase cancer risk [56]. Likewise, a literature search for the genotype prevalence of other HPV-related cancers within East Asia and Taiwan, showed that in Northern China HPV16 is dominant in cervical carcinomas [75, 76]. This prevalence decreases in the Southern regions of the country [20, 77], to become similar to the observed occurrence of other important genotypes in Taiwan, HPV52 and HPV58 [78]. Although HPV16 has been found in over 80% of cervical cancers in Asia, it has been reported in only 50% of the cervical cancers in Taiwan [79]. In this population, HPV58 is the second or the third most common genotype found in cervical malignancies (∼20%, together with HPV18) and cervical HPV infections, but it is rare in other parts of the world [80-83]. HPV52, usually found in cervical cancers in East Asia [59, 84, 85] was included in the test panel, but was not detected in our Taiwanese cohort. This particular distribution of HPV genotypes suggests tissue or geographic specificity for HPV16 and HPV58 in Taiwan, where HPV is becoming an increasingly important etiologic factor of OPSCC.

As part of our study, we also analyzed the clinical determinants of HPV positivity in OPSCC. This malignancy is considered one of the emerging causes of cancer death in Asian and Taiwan males [9, 55]. Curiously, even when most of the cases in our cohort were males (94%), females had most commonly HPV-positive tumors and presented reduced alcohol, smoking, or betel quid habits. There was also a trend for slightly better prognosis within the female group (S2 Fig). While the better prognosis in women may be related to better health utilization of women over men, our results indicate that OPSCC is a disease that also affects Taiwanese women, and more research is necessary to address their specific clinical management. Moreover, HPV-positive tumors were associated with a slightly higher age at diagnosis (mean = 55.5 years) than HPV-negative tumors. This represents another particularity since hrHPV has been historically associated with the onset of OPSCCs in younger, middle-aged individuals [23, 35, 39, 44, 59, 86], even in Taiwan, where the average age at diagnosis has been reported between 40 and 50 years of age [55]. Despite this belief, new studies from Western countries have contradicted these trends by demonstrating that the highest HPV prevalence occurs in patients above 55 and even 70 years of age [19, 87, 88]. This suggests a similar changing epidemiology in Taiwan, with a possible shift of sexual onset at a later age, reduced number of sexual partners in life, or reinfection later in life. Even when our data diverge from previous reports, these agree with studies from Southern China and Western countries indicating a significant correlation between HPV status and earlier primary tumor stage (T-stage) [3, 59]. However, our results demonstrate a higher T-stage within the entire cohort (T3-4 instead of T1-2), although the most common T class among HPV-positive tumors was T2 (46%) but among HPV-negative tumors the most common T-class was T4 (37%). N-stage did not present a significant association with HPV status. These HPV-positive tumors were associated with the tonsils as their primary location, denoting site-specificity, as also indicated in previous observations from Taiwan [55]. Thus, our results indicate that HPV-driven OPSCC is a clinically unique disease with distinctive features to this cohort of Taiwanese patients.

Furthermore, the carcinogenic effect of the risk factors alcohol, smoking, and betel quid has been very well characterized. These are considered the main etiological agents of HPV-negative OPSCC [6-13-15-17]. It has been proposed that the observed increasing rates of OPSCC in East Asia and Taiwan are related to the excessive and extensive use of alcohol, smoking, and betel quid (which does not contain tobacco) [9, 51-54-63]. Previous studies in Taiwan report that exposure to these agents represents a high risk for developing primarily intraoral cancer [9, 57]. For this reason, it was not surprising to observe that above two thirds of the OPSCC cases in our study tested negative for HPV and had a prominent smoking (83%), betel quid chewing (65%), or alcohol drinking (51%) history (Table 2, S4 Table). Earlier findings also showed a similar proportion of alcohol drinkers within the Taiwanese male population [9]. Interestingly, in our cohort, 87% of the individuals with OPSCC used these chemical carcinogens in combination. Only 1% of the cohort consumed alcohol or betel quid alone, and 11% only smoked (S4 Table). Although these historically causative risk factors are still very prominent in Taiwan, we observed that only alcohol consumption is a significant determinant of worse prognosis. This was surprising, as tobacco smoking has been found to be the most important predictor of unfavorable outcome and risk factor in OPSCC, as patients with mutant p53 have a reduced capacity to repair DNA [3, 25, 27, 44, 89]. The reasons for the dwarfed effect of smoking and betel quid as negative drivers remain elusive, but they may reside in the high proportion of smokers in the cohort (83%) or individuals with a complex combination of risk factors, as shown in S4 Table. This makes it impossible to separate the individual effects of the carcinogens, alcohol, tobacco, and betel quid, which could have additive effects on DNA damage and potentiate the negative impact of alcohol. Alternatively, alcohol could also represent a surrogate for lack of social support or socioeconomic disadvantage. Additionally, our findings reveal that HPV-positive tumors were less exposed to these three risk factors, which correlates with previous reports indicating that individuals with HPV-positive OPSCCs are more likely to be never or former smokers, or drinkers [25, 57, 59, 89, 90]. Nonetheless, these are common risk factors for HPV-positive OPSCC, and in Western countries, 10-30% of OPSCC tumors occur in individuals that smoke or drink [27, 91]. In this study, we observed a transformation from cancer that is derived from only smoking, betel quid, and alcohol influence, to an increasing trend of HPV-positive cases. The strong association of HPV with better outcome has been widely reported [23, 25, 27, 28, 35, 42, 44-50]. Our results, showing that HPV-positive OPSCCs have a strikingly better prognosis, agree with these and previous work from South-East Asia and Taiwan [57, 59]. Additional evidence showed that in our Taiwan cohort, even when having poorer outcomes, drinkers, smokers, and betel quid users benefited from the simultaneous presence of HPV, displaying higher OPSCC survival rates than those HPV-negative. Similar interactions have been seen before, especially for smoking or tobacco use and HPV [27, 89]. These results may have a substantial impact on the clinical management of OPSCC patients in Taiwan and their risk stratification. HPV-positive individuals could benefit from the secession of alcohol, smoking, or betel quid habits and therapy de-escalation, as it is currently tested in diverse hospital settings to reduce toxicity and post-treatment morbidity [13, 49, 50, 59, 92, 93].

To our knowledge, this is the first comprehensive study analyzing the impact of HPV-driven OPSCC in Taiwan. Our approach of matching HPV status and prevalence data to clinical features, risk behavior exposure, and clinical outcomes represents a distinctive strength of our study, adding to our understanding of an under-represented ethnic group in Western epidemiological studies. Also, double p16-HPV DNA testing provided clinical relevance to the findings, as the concordance of the results was very high. Because p16 as a surrogate marker for HPV-driven OPSCCs can account for more than 5% false positives [10, 94, 95], we did not rely solely on this test and performed DNA testing for HPV. Our study also suffers from limitations. Cases from a single center, the Chang Gung Memorial Hospital at Linkou, in Taiwan, were included in the study. However, it is the first and largest cancer center in Taiwan, providing cancer care to roughly a quarter of the country’s cancer patients which is a reasonably valid representation of the Taiwanese population. Although, limits to access were not assessed, which could introduce bias. An additional limitation is represented by the lack of both p16 and HPV-DNA results for all the specimens, as several were missing or had insufficient amount of tissue available for analysis (see Fig 1, S1 Table). The retrospective nature of our study also presented challenges to this work. We could not evaluate the impact of therapy because changes to treatment strategies occurred over time, and we were unable to follow changes in sexual behaviors that could perhaps help explain the occurrence of HPV in Taiwan. Data on risk factors were limited, as we could not retrieve the exact amounts and type of alcohol, cigarettes, or betel quid consumed, nor information on previous infections with hrHPV, or comorbidities. For the same reason, disease specific survival times could not be calculated. The AJCC TNM staging system changed to include HPV status in 2017 [96]; however, our cases were classified at diagnosis following previous guidance, and staging changes were not reflected in the reported N-status.

Future research should incorporate efforts to further characterize HPV-driven OPSCC in Taiwan. Since our cohort included a limited number of females, future studies aiming at clarifying the implication of this disease in women, who also suffer from HPV-driven cervical and genital cancers, are needed. Importantly, studies related to the collection of specific epidemiological data will be relevant to guide new public health policies. Since oral hrHPV infections precede the development of OPSCC, research on the natural history of this disease, including the prevalence of HPV genotypes, will help strengthen the current HPV vaccination efforts. The HPV vaccine was recently introduced in Taiwan nationwide, targeting only prepubescent girls. Still, our studies reflect that boys will also benefit from vaccination, as it has the potential to halt the expansion of HPV-positive OPSCC and other cancers [40, 56, 59]. HPV vaccination uptake is presumably low, and initially only the bivalent (HPV 16/18) and tetravalent (HPV 6/11/16/18) vaccines were used. Our data indicate that in this Taiwanese population HPV58 was the second most common hrHPV genotype. The nonavalent vaccine, which protects against HPV58 and other 8 HPV genotypes (HPV 6/11/16/18/31/33/45/52/58), would be the most appropriate choice for this target population based on our study. However, even that vaccine (without assuming cross reactivity) would not cover HPV35, 31, 59, and 66, which we found respectively account for 2.2%, 1.5%, 1.5%, and 1% each. New prospective studies have the potential to shed light on the risk of OPSCC within the vaccinated population and have broad public health implications for control measures. Lastly, multicenter, nationwide studies will provide an understanding of the variables leading to HPV oncogenesis in this population, improving risk definition, outcome prediction, and patient stratification. All these are necessary to provide patients with appropriate care based on their HPV status. Additional investigations into the molecular mechanism of HPV-induced OPSCC are also required to address the risk of recurrence and progression in Taiwan. Based on risk evaluation, HPV-positive OPSCC patients may be candidates for therapy de-intensification. Tailored interventions should also be designed based on longitudinal investigations of the interaction between HPV and risk behaviors, including drinking, smoking, and chewing betel quid.

In conclusion, our retrospective study provides empirical evidence on the impact of hrHPV on OPSCCs in a large cohort from Taiwan. We found that HPV is present and likely an increasing etiological factor in these Taiwanese individuals with OPSCC. These observations may represent continuous behavioral changes in Taiwan. HPV positivity is associated with significantly better outcomes. Involvement in risk habits, alcohol, cigarette, and betel quid use was still widespread, but the substantial prognostic benefit of HPV remained present. Thus, this consistent trend reflects the need for policies and sustained public health interventions aiming to improve the management and prevention of HPV-driven OPSCC in Taiwan.

## Supporting information

S1 Checklist [[supplements/213231_file07.pdf]](pending:yes)

S1 Figure [[supplements/213231_file08.tif]](pending:yes)

S2 Figure [[supplements/213231_file09.tif]](pending:yes)

S3 Figure [[supplements/213231_file10.tif]](pending:yes)

S4 Figure [[supplements/213231_file11.tif]](pending:yes)

S2 Table [[supplements/213231_file12.docx]](pending:yes)

S3 Table [[supplements/213231_file13.docx]](pending:yes)

S4 Table [[supplements/213231_file14.docx]](pending:yes)

S5 Table [[supplements/213231_file15.docx]](pending:yes)

S6 Table [[supplements/213231_file16.docx]](pending:yes)

S7 Table [[supplements/213231_file17.docx]](pending:yes)

S8 Table [[supplements/213231_file18.docx]](pending:yes)

S1 Table [[supplements/213231_file19.xlsx]](pending:yes)

## Data Availability

All relevant data needed to reproduce our findings are included in this manuscript and its Supporting Information, excepting sensitive dates that could allow the identification of the patients in this study.

## Supporting information

**S1 Checklist. Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) report**.

**S1 Fig. Representative p16 immunostaining for OPSCC tissue sections**. Specimens from example p16-negative (P0283) and p16-positive (P0267) tumors are displayed. p16 expression is observed as a brown nuclear and cytoplasmic coloration. Magnification, 200x.

**S2 Fig. Kaplan-Meier plots for survival outcomes by gender**. Up to 5-year overall survival (OS) and disease-free survival (DFS) outcomes were analyzed within the whole cohort by the Kaplan-Meier method and log-rank test (*p*-values).

**S3 Fig. Kaplan-Meier plots for survival outcomes by alcohol, smoking, and betel quid**. (A-C) Comparison of prognostic outcomes of (A) alcohol, (B) smoke, and (C) betel quid between the whole cohort and HPV risk groups. HPV positivity is defined as HPV DNA-positive and/or p16-positive. Up to 5-year overall survival (OS, top) and disease-free survival (DFS, bottom) probabilities were analyzed by the Kaplan-Meier method and log-rank test (*p*-values), as displayed for each risk group. HPV-, HPV-negative; HPV+, HPV-positive.

**S4 Fig. HPV-positive OPSCC is associated with increased disease-free survival time in the presence of other risk factors**. (A-C) Up to 5-year disease-free survival (DFS) prognostic outcome of the HPV variable within alcohol, and smoking and/or betel quid risk groups. HPV positivity is defined as HPV DNA-positive and/or p16-positive. The DFS smoking and betel quid variables were not analyzed individually due to low number of events. (A) Table includes the multivariable hazard probabilities analyzed using Cox survival models and hazard ratio (HR) estimations, adjusted for age, T- and N-stage, and which were visualized by forest plots. The complete analysis is found in S8 Table. (B-C) Kaplan-Meier survival analysis. Plots represent the DFS probabilities of cases stratified by HPV status within the (B) alcohol and (C) smoke and/or betel quid groups. Left, plots showing cases with alcohol or smoke and/or betel quid consumption. Right, plots showing cases without exposition to alcohol or smoke and/or betel quid. Log-rank analysis was used to compare the survival distributions (log-rank *p*-values are in the plots). HPV-, HPV-negative; HPV+, HPV-positive.

**S1 Table. HPV status and patient data**.

**S2 Table. p16 vs. HPV DNA results**.

**S3 Table. Yearly HPV occurrence**.

**S4 Table. Risk factors exposure characteristics, and differences regarding HPV status**.

**S5 Table. Clinical and demographic characteristics by gender**.

**S6 Table. Multivariable survival analysis of the whole cohort**.

**S7 Table. Multivariable survival analysis for alcohol, smoking, and quid within HPV r isk groups, adjusted for age, T- and N-stage**.

**S8 Table. Multivariable disease-free survival analysis by HPV within alcohol, and smoking and/or betel quid risk groups (controlled for alcohol, and smoking and/or betel quid). All models control for age, T- and N-stage**.

## Author Contributions

**Conceptualization:**Thomas E. Carey, Kai-Ping Chang, Heather M. Walline, Guadalupe Lorenzatti Hiles, Emily L. Bellile.

**Data Curation:** Guadalupe Lorenzatti Hiles, Emily L. Bellile, Chun-I Wang, Wei-Chen Yen, Lisa M. Pinatti, Heather M. Walline.

**Formal Analysis**: Guadalupe Lorenzatti Hiles, Emily L. Bellile, Heather M. Walline.

**Funding Acquisition:** Guadalupe Lorenzatti Hiles, Kai-Ping Chang, Thomas E. Carey, Heather M. Walline.

**Investigation:** Guadalupe Lorenzatti Hiles, Chun-I Wang, Wei-Chen Yen, Christine M. Goudsmit, Hannah L. Briggs, Trey B. Thomas, Lila Peters, Macy A. Afsari, Anna C. Morris, Nadine Jawad, Heather M. Walline.

**Project Administration:** Thomas E. Carey, Kai-Ping Chang, Heather M. Walline, Guadalupe Lorenzatti Hiles, Christine M. Goudsmit.

**Resources:** Kai-Ping Chang, Thomas E. Carey, Heather M. Walline.

**Supervision:** Kai-Ping Chang, Thomas E. Carey, Heather M. Walline.

**Validation:** Guadalupe Lorenzatti Hiles, Emily L. Bellile, Hannah L. Briggs, Macy A. Afsari, Heather M. Walline.

**Visualization:** Guadalupe Lorenzatti Hiles, Emily L. Bellile, Macy A. Afsari.

**Writing – Original Draft Preparation:** Guadalupe Lorenzatti Hiles.

**Writing – Review & Editing:** Guadalupe Lorenzatti Hiles, Kai-Ping Chang, Emily L. Bellile, Chun-I Wang, Wei-Chen Yen, Christine M. Goudsmit, Hannah L. Briggs, Trey B. Thomas, Lila Peters, Macy A. Afsari, Lisa M. Pinatti, Anna C. Morris, Nadine Jawad, Thomas E. Carey, Heather M. Walline.

## Data Availability Statement

All relevant data needed to reproduce our findings are included in this manuscript and its Supporting Information, excepting sensitive dates that could allow the identification of the patients in this study.

## Funding

This study was funded by the University of Michigan-Chang Gung Memorial Hospital Pilot Grant ([https://www.rogelcancercenter.org](https://www.rogelcancercenter.org) and [https://www.cgmh.org.tw/en](https://www.cgmh.org.tw/en)) to K-PC and TEC; the National Cancer Institute at the National Institutes of Health ([https://www.cancer.gov](https://www.cancer.gov)), CA194536 to TEC and HMW, and CA194536-S1 to GLH; the Chang Gung Memorial Hospital, CMRPG3H0852, CMRPG3J1251, and CORPG3G0171 to K-PC; the Taiwan Ministry of Science and Technology ([https://www.most.gov.tw](https://www.most.gov.tw)), MOST 108-2314-B-182A-108-MY3 to K-PC; and, funds from the University of Michigan Undergraduate Research Opportunity Program ([https://lsa.umich.edu/urop](https://lsa.umich.edu/urop)) to GLH, TEC and HMW. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Competing interests

The authors have declared that no competing interests exist.

## Acknowledgments

We would like to thank the patients who generously contributed to the study, and the members of the Cancer Center, Chang Gung Memorial Hospital, for their invaluable help.

## Footnotes

*   This manuscript was sent back to us because during the screening process it came as containing too much identifying information. By error, we attached an S1 Table file that still contained the ages of the individuals in our cohort. This has been solved, leaving all minimal data necessary to replicate the study (notice that age and genre are key to the results). A new S1 Table file has been attached. The cases were unidentified in Taiwan and the sample/patient IDs were not known to anyone outside the Taiwanese research team. No other researcher from the Michigan team or outside the research group had access to patient identification. Thank you for your help!

## Abbreviations

CGMH
:   Chang Gung Memorial Hospital
DFS
:   disease-free survival
EQUATOR
:   Enhancing the QUAlity and Transparency Of health Research Network
FFPE
:   formalin-fixed, paraffin-embedded
HPV
:   human papillomavirus
hrHPV
:   high-risk human papillomavirus
IHC
:   immunohistochemical staining
OS
:   overall survival
OPSCC
:   oropharyngeal squamous cell carcinoma
PCR-MA
:   multiplex PCR-MassArray
STROBE
:   Strengthening the Reporting of Observational Studies in Epidemiology

*   Received October 26, 2020.
*   Revision received March 10, 2021.
*   Accepted March 12, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3322/caac.21492&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30207593&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

2.  2.D’Souza G, Kreimer AR, Viscidi R, Pawlita M, Fakhry C, Koch WM, et al. Case-control study of human papillomavirus and oropharyngeal cancer. N Engl J Med. 2007;356(19):1944–56.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa065497&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17494927&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246305500006&link_type=ISI) 

3.  3.Elrefaey S, Massaro MA, Chiocca S, Chiesa F, Ansarin M. HPV in oropharyngeal cancer: the basics to know in clinical practice. Acta Otorhinolaryngol Ital. 2014;34(5):299–309.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25709145&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

4.  4.Shiboski CH, Schmidt BL, Jordan RC. Tongue and tonsil carcinoma: increasing trends in the U.S. population ages 20-44 years. Cancer. 2005;103(9):1843–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cncr.20998&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15772957&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228423100012&link_type=ISI) 

5.  5.Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69(1):7–34.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3322/caac.21551&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30620402&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

6.  6.Blot WJ, McLaughlin JK, Winn DM, Austin DF, Greenberg RS, Preston-Martin S, et al. Smoking and drinking in relation to oral and pharyngeal cancer. Cancer Res. 1988;48(11):3282–7.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI0OC8xMS8zMjgyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

7.  7.Dayyani F, Etzel CJ, Liu M, Ho CH, Lippman SM, Tsao AS. Meta-analysis of the impact of human papillomavirus (HPV) on cancer risk and overall survival in head and neck squamous cell carcinomas (HNSCC). Head Neck Oncol. 2010;2:15.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1758-3284-2-15&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20587061&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

8.  8.Hashibe M, Brennan P, Chuang SC, Boccia S, Castellsague X, Chen C, et al. Interaction between tobacco and alcohol use and the risk of head and neck cancer: pooled analysis in the International Head and Neck Cancer Epidemiology Consortium. Cancer Epidemiol Biomarkers Prev. 2009;18(2):541–50.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY2VicCI7czo1OiJyZXNpZCI7czo4OiIxOC8yLzU0MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzEyLzIwMjAuMTAuMjYuMjAyMTMyMzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

9.  9.Ho PS, Ko YC, Yang YH, Shieh TY, Tsai CC. The incidence of oropharyngeal cancer in Taiwan: an endemic betel quid chewing area. J Oral Pathol Med. 2002;31(4):213–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1034/j.1600-0714.2002.310404.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12076324&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000175549500004&link_type=ISI) 

10. 10.Nauta IH, Rietbergen MM, van Bokhoven AAJD, Bloemena E, Lissenberg-Witte BI, Heideman DAM, et al. Evaluation of the eighth TNM classification on p16-positive oropharyngeal squamous cell carcinomas in the Netherlands and the importance of additional HPV DNA testing. Ann Oncol. 2018;29(5):1273–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/annonc/mdy060&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

11. 11.Sturgis EM, Ang KK. The epidemic of HPV-associated oropharyngeal cancer is here: is it time to change our treatment paradigms? J Natl Compr Canc Netw. 2011;9(6):665–73.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiam5jY24iO3M6NToicmVzaWQiO3M6NzoiOS82LzY2NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzEyLzIwMjAuMTAuMjYuMjAyMTMyMzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

12. 12.Talamini R, Vaccarella S, Barbone F, Tavani A, La Vecchia C, Herrero R, et al. Oral hygiene, dentition, sexual habits and risk of oral cancer. Br J Cancer. 2000;83(9):1238–42.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1054/bjoc.2000.1398&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11027440&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000090107200019&link_type=ISI) 

13. 13.Westra WH. The morphologic profile of HPV-related head and neck squamous carcinoma: implications for diagnosis, prognosis, and clinical management. Head Neck Pathol. 2012;6 Suppl 1:S48–54.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12105-012-0371-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22782223&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

14. 14.Znaor A, Brennan P, Gajalakshmi V, Mathew A, Shanta V, Varghese C, et al. Independent and combined effects of tobacco smoking, chewing and alcohol drinking on the risk of oral, pharyngeal and esophageal cancers in Indian men. Int J Cancer. 2003;105(5):681–6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.11114&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12740918&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183474000016&link_type=ISI) 

15. 15.IARC Working Group on the Evaluation of Carcinogenic Risks to Humans, World Health Organization., International Agency for Research on Cancer. Betel-quid and areca-nut chewing and some areca-nut-derived nitrosamines. Lyon, France Geneva: IARC Press ;Distributed by IARC Press and the World Health Organization Marketing and Dissemination; 2004. x, 334 p. p.
    
    

16. 16.IARC Working Group on the Evaluation of Carcinogenic Risks to Humans., International Agency for Research on Cancer., National Cancer Institute (U.S.). Alcohol drinking. Lyon, France: World Health Organization, International Agency for Research on Cancer; 1988. 416 p. p.
    
    

17. 17.IARC Working Group on the Evaluation of the Carcinogenic Risk of Chemicals to Humans, International Agency for Research on Cancer. Tobacco habits other than smoking, betel-quid and areca-nut chewing, and some related nitrosamines. Lyon, France: International Agency for Research on Cancer; 1985. 291 pages, 4 pages of plates p.
    
    

18. 18.Ng M, Freeman MK, Fleming TD, Robinson M, Dwyer-Lindgren L, Thomson B, et al. Smoking prevalence and cigarette consumption in 187 countries, 1980-2012. JAMA. 2014;311(2):183–92.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2013.284692&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24399557&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000329339400024&link_type=ISI) 

19. 19.Zumsteg ZS, Cook-Wiens G, Yoshida E, Shiao SL, Lee NY, Mita A, et al. Incidence of Oropharyngeal Cancer Among Elderly Patients in the United States. JAMA Oncol. 2016;2(12):1617–23.
    
    

20. 20.Allen CT, Lewis JS, El-Mofty SK, Haughey BH, Nussenbaum B. Human papillomavirus and oropharynx cancer: biology, detection and clinical implications. Laryngoscope. 2010;120(9):1756–72.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/lary.20936&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20669304&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281430600008&link_type=ISI) 

21. 21.Fakhry C, Gillison ML. Clinical implications of human papillomavirus in head and neck cancers. J Clin Oncol. 2006;24(17):2606–11.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjEwOiIyNC8xNy8yNjA2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

22. 22.Gillison ML, Koch WM, Capone RB, Spafford M, Westra WH, Wu L, et al. Evidence for a causal association between human papillomavirus and a subset of head and neck cancers. J Natl Cancer Inst. 2000;92(9):709–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jnci/92.9.709&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10793107&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000086818100013&link_type=ISI) 

23. 23.Gillison ML, D’Souza G, Westra W, Sugar E, Xiao W, Begum S, et al. Distinct risk factor profiles for human papillomavirus type 16-positive and human papillomavirus type 16-negative head and neck cancers. J Natl Cancer Inst. 2008;100(6):407–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jnci/djn025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18334711&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000254296900012&link_type=ISI) 

24. 24.Gillison ML, Chaturvedi AK, Anderson WF, Fakhry C. Epidemiology of Human Papillomavirus-Positive Head and Neck Squamous Cell Carcinoma. J Clin Oncol. 2015;33(29):3235–42.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjEwOiIzMy8yOS8zMjM1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

25. 25.Kumar B, Cordell KG, Lee JS, Worden FP, Prince ME, Tran HH, 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(19):3128–37.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjEwOiIyNi8xOS8zMTI4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

26. 26.Licitra L, Bossi P, Locati LD. A multidisciplinary approach to squamous cell carcinomas of the head and neck: what is new? Curr Opin Oncol. 2006;18(3):253–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/01.cco.0000219254.53091.35&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16552237&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

27. 27.Maxwell JH, Kumar B, Feng FY, Worden FP, Lee JS, Eisbruch A, et al. Tobacco use in human papillomavirus-positive advanced oropharynx cancer patients related to increased risk of distant metastases and tumor recurrence. Clin Cancer Res. 2010;16(4):1226–35.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNsaW5jYW5yZXMiO3M6NToicmVzaWQiO3M6OToiMTYvNC8xMjI2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

28. 28.Worden FP, Kumar B, Lee JS, Wolf GT, Cordell KG, Taylor JM, et al. Chemoselection as a strategy for organ preservation in advanced oropharynx cancer: response and survival positively associated with HPV16 copy number. J Clin Oncol. 2008;26(19):3138–46.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjEwOiIyNi8xOS8zMTM4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

29. 29.IARC Working Group on the Evaluation of Carcinogenic Risks to Humans., International Agency for Research on Cancer. Human papillomaviruses. Lyon, France Geneva, Switzerland: World Health Organization Distributed by WHO Press; 2007. viii, 670 p. p.
    
    

30. 30.Begum S, Cao D, Gillison M, Zahurak M, Westra WH. Tissue distribution of human papillomavirus 16 DNA integration in patients with tonsillar carcinoma. Clin Cancer Res. 2005;11(16):5694–9.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNsaW5jYW5yZXMiO3M6NToicmVzaWQiO3M6MTA6IjExLzE2LzU2OTQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMy8xMi8yMDIwLjEwLjI2LjIwMjEzMjMxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

31. 31.Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El Ghissassi F, et al. A review of human carcinogens--Part B: biological agents. Lancet Oncol. 2009;10(4):321–2.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1470-2045(09)70096-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19350698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264812300014&link_type=ISI) 

32. 32.Mork J, Lie AK, Glattre E, Hallmans G, Jellum E, Koskela P, et al. Human papillomavirus infection as a risk factor for squamous-cell carcinoma of the head and neck. N Engl J Med. 2001;344(15):1125–31.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM200104123441503&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11297703&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167990600003&link_type=ISI) 

33. 33.Rampias T, Sasaki C, Weinberger P, Psyrri A. E6 and e7 gene silencing and transformed phenotype of human papillomavirus 16-positive oropharyngeal cancer cells. J Natl Cancer Inst. 2009;101(6):412–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jnci/djp017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19276448&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264393100009&link_type=ISI) 

34. 34.Stransky N, Egloff AM, Tward AD, Kostic AD, Cibulskis K, Sivachenko A, et al. The mutational landscape of head and neck squamous cell carcinoma. Science. 2011;333(6046):1157–60.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMzMvNjA0Ni8xMTU3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

35. 35.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(4):612–9.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjg6IjI2LzQvNjEyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

36. 36.Chaturvedi AK, Engels EA, Pfeiffer RM, Hernandez BY, Xiao W, Kim E, et al. Human papillomavirus and rising oropharyngeal cancer incidence in the United States. J Clin Oncol. 2011;29(32):4294–301.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjEwOiIyOS8zMi80Mjk0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

37. 37.Ernster JA, Sciotto CG, O’Brien MM, Finch JL, Robinson LJ, Willson T, et al. Rising incidence of oropharyngeal cancer and the role of oncogenic human papilloma virus. Laryngoscope. 2007;117(12):2115–28.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/MLG.0b013e31813e5fbb&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17891052&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000251397200007&link_type=ISI) 

38. 38.Stein AP, Saha S, Yu M, Kimple RJ, Lambert PF. Prevalence of human papillomavirus in oropharyngeal squamous cell carcinoma in the United States across time. Chem Res Toxicol. 2014;27(4):462–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/tx500034c&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24641254&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

39. 39.Chaturvedi AK, Anderson WF, Lortet-Tieulent J, Curado MP, Ferlay J, Franceschi S, et al. Worldwide trends in incidence rates for oral cavity and oropharyngeal cancers. J Clin Oncol. 2013;31(36):4550–9.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjEwOiIzMS8zNi80NTUwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

40. 40.Gooi Z, Chan JY, Fakhry C. The epidemiology of the human papillomavirus related to oropharyngeal head and neck cancer. Laryngoscope. 2016;126(4):894–900.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/lary.25767&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26845348&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

41. 41.Hammarstedt L, Lindquist D, Dahlstrand H, Romanitan M, Dahlgren LO, Joneberg J, et al. Human papillomavirus as a risk factor for the increase in incidence of tonsillar cancer. Int J Cancer. 2006;119(11):2620–3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.22177&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16991119&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241964600018&link_type=ISI) 

42. 42.Boscolo-Rizzo P, Pawlita M, Holzinger D. From HPV-positive towards HPV-driven oropharyngeal squamous cell carcinomas. Cancer Treat Rev. 2016;42:24–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ctrv.2015.10.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26547133&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

43. 43.Gatta G, Botta L, Sánchez MJ, Anderson LA, Pierannunzio D, Licitra L, et al. Prognoses and improvement for head and neck cancers diagnosed in Europe in early 2000s: The EUROCARE-5 population-based study. Eur J Cancer. 2015;51(15):2130–43.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejca.2015.07.043&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26421817&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

44. 44.Ang KK, Harris J, Wheeler R, Weber R, Rosenthal DI, Nguyen-Tân PF, et al. Human papillomavirus and survival of patients with oropharyngeal cancer. N Engl J Med. 2010;363(1):24–35.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa0912217&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20530316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279343000005&link_type=ISI) 

45. 45.Fakhry C, Westra WH, Li S, Cmelak A, Ridge JA, Pinto H, et al. 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(4):261–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jnci/djn011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18270337&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253796600009&link_type=ISI) 

46. 46.Fakhry C, Westra WH, Wang SJ, van Zante A, Zhang Y, Rettig E, et al. The prognostic role of sex, race, and human papillomavirus in oropharyngeal and nonoropharyngeal head and neck squamous cell cancer. Cancer. 2017;123(9):1566–75.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cncr.30353&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28241096&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

47. 47.Haughey BH, Hinni ML, Salassa JR, Hayden RE, Grant DG, Rich JT, et al. Transoral laser microsurgery as primary treatment for advanced-stage oropharyngeal cancer: a United States multicenter study. Head Neck. 2011;33(12):1683–94.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hed.21669&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21284056&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

48. 48.Horne ZD, Glaser SM, Vargo JA, Ferris RL, Balasubramani GK, Clump DA, et al. Confirmation of proposed human papillomavirus risk-adapted staging according to AJCC/UICC TNM criteria for positive oropharyngeal carcinomas. Cancer. 2016;122(13):2021–30.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cncr.30021&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27111669&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

49. 49.Quon H, Cohen MA, Montone KT, Ziober AF, Wang LP, Weinstein GS, et al. Transoral robotic surgery and adjuvant therapy for oropharyngeal carcinomas and the influence of p16 INK4a on treatment outcomes. Laryngoscope. 2013;123(3):635–40.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/lary.22172&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22588642&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

50. 50.Rietbergen MM, Brakenhoff RH, Bloemena E, Witte BI, Snijders PJ, Heideman DA, et al. Human papillomavirus detection and comorbidity: critical issues in selection of patients with oropharyngeal cancer for treatment De-escalation trials. Ann Oncol. 2013;24(11):2740–5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/annonc/mdt319&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23946330&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000326643200006&link_type=ISI) 

51. 51.Chuang SL, Su WW, Chen SL, Yen AM, Wang CP, Fann JC, et al. Population-based screening program for reducing oral cancer mortality in 2,334,299 Taiwanese cigarette smokers and/or betel quid chewers. Cancer. 2017;123(9):1597–609.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cncr.30517&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28055109&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

52. 52.Hsu KY, Tsai YF, Huang CC, Yeh WL, Chang KP, Lin CC, et al. Tobacco-Smoking, Alcohol-Drinking, and Betel-Quid-Chewing Behaviors: Development and Use of a Web-Based Survey System. JMIR Mhealth Uhealth. 2018;6(6):e142.
    
    

53. 53.Lee YA, Li S, Chen Y, Li Q, Chen CJ, Hsu WL, et al. Tobacco smoking, alcohol drinking, betel quid chewing, and the risk of head and neck cancer in an East Asian population. Head Neck. 2019;41(1):92–102.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hed.25383&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30552826&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

54. 54.Shiu MN, Chen TH. Impact of betel quid, tobacco and alcohol on three-stage disease natural history of oral leukoplakia and cancer: implication for prevention of oral cancer. Eur J Cancer Prev. 2004;13(1):39–45.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00008469-200402000-00007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15075787&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000220286000007&link_type=ISI) 

55. 55.Hwang TZ, Hsiao JR, Tsai CR, Chang JS. Incidence trends of human papillomavirus-related head and neck cancer in Taiwan, 1995-2009. Int J Cancer. 2015;137(2):395–408.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.29330&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25395239&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

56. 56.Wang CP, Chen TC, Chen HH, Hsu WL, Chang YL. Prevalence of current oral HPV infection among 100 betel nut chewers or cigarette smokers in Northern Taiwan. J Formos Med Assoc. 2019;118(1 Pt 1):203–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jfma.2018.03.013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29636296&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

57. 57.Chien CY, Su CY, Fang FM, Huang HY, Chuang HC, Chen CM, et al. Lower prevalence but favorable survival for human papillomavirus-related squamous cell carcinoma of tonsil in Taiwan. Oral Oncol. 2008;44(2):174–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.oraloncology.2007.01.018&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17418614&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000254105500011&link_type=ISI) 

58. 58.Huang H, Zhang B, Chen W, Zou SM, Zhang YX, Qiao YL. [Detection of human papillomavirus in oropharyngeal squamous cell carcinoma]. Zhongguo Yi Xue Ke Xue Yuan Xue Bao. 2012;34(6):545–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3881/j.issn.1000-503X.2012.06.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23286396&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

59. 59.Lam EW, Chan JY, Chan AB, Ng CS, Lo ST, Lam VS, et al. Prevalence, Clinicopathological Characteristics, and Outcome of Human Papillomavirus-Associated Oropharyngeal Cancer in Southern Chinese Patients. Cancer Epidemiol Biomarkers Prev. 2016;25(1):165–73.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY2VicCI7czo1OiJyZXNpZCI7czo4OiIyNS8xLzE2NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzEyLzIwMjAuMTAuMjYuMjAyMTMyMzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

60. 60.Li W, Tran N, Lee SC, O’Brien CJ, Tse GM, Scolyer RA, et al. New evidence for geographic variation in the role of human papillomavirus in tonsillar carcinogenesis. Pathology. 2007;39(2):217–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00313020701230823&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17454751&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

61. 61.Maruyama H, Yasui T, Ishikawa-Fujiwara T, Morii E, Yamamoto Y, Yoshii T, et al. Human papillomavirus and p53 mutations in head and neck squamous cell carcinoma among Japanese population. Cancer Sci. 2014;105(4):409–17.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/cas.12369&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24521534&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

62. 62.No JH, Sung MW, Hah JH, Choi SH, Lee MC, Kim HS, et al. Prevalence and prognostic value of human papillomavirus genotypes in tonsillar squamous cell carcinoma: a Korean multicenter study. Cancer. 2015;121(4):535–44.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cncr.29086&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25283642&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

63. 63.Ko YC, Chiang TA, Chang SJ, Hsieh SF. Prevalence of betel quid chewing habit in Taiwan and related sociodemographic factors. J Oral Pathol Med. 1992;21(6):261–4.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1600-0714.1992.tb01007.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1501158&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992JE87100005&link_type=ISI) 

64. 64.Edge SB, American Joint Committee on Cancer., American Cancer Society. AJCC cancer staging handbook: from the AJCC cancer staging manual. 7th ed. New York: Springer; 2010. xix, 718 p. p.
    
    

65. 65.Walline HM, Komarck C, McHugh JB, Byrd SA, Spector ME, Hauff SJ, et al. High-risk human papillomavirus detection in oropharyngeal, nasopharyngeal, and oral cavity cancers: comparison of multiple methods. JAMA Otolaryngol Head Neck Surg. 2013;139(12):1320–7.
    
    

66. 66.Walline HM, Carey TE, Goudsmit CM, Bellile EL, D’Souza G, Peterson LA, et al. High-Risk HPV, Biomarkers, and Outcome in Matched Cohorts of Head and Neck Cancer Patients Positive and Negative for HIV. Mol Cancer Res. 2017;15(2):179–88.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibW9sY2FucmVzIjtzOjU6InJlc2lkIjtzOjg6IjE1LzIvMTc5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

67. 67.Walline HM, Goudsmit CM, McHugh JB, Tang AL, Owen JH, Teh BT, et al. Integration of high-risk human papillomavirus into cellular cancer-related genes in head and neck cancer cell lines. Head Neck. 2017;39(5):840–52.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hed.24729&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28236344&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

68. 68.Lewis JS, Beadle B, Bishop JA, Chernock RD, Colasacco C, Lacchetti C, et al. Human Papillomavirus Testing in Head and Neck Carcinomas: Guideline From the College of American Pathologists. Arch Pathol Lab Med. 2018;142(5):559–97.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.5858/arpa.2017-0286-CP&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

69. 69.Moore EJ, Olsen SM, Laborde RR,García JJ, Walsh FJ, Price DL, et al. Long-term functional and oncologic results of transoral robotic surgery for oropharyngeal squamous cell carcinoma. Mayo Clin Proc. 2012;87(3):219–25.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mayocp.2011.10.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22386176&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

70. 70.Yeh DH, Tam S, Fung K, MacNeil SD, Yoo J, Winquist E, et al. Transoral robotic surgery vs. radiotherapy for management of oropharyngeal squamous cell carcinoma-A systematic review of the literature. Eur J Surg Oncol. 2015;41(12):1603–14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejso.2015.09.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26461255&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

71. 71.Li W, Thompson CH, O’Brien CJ, McNeil EB, Scolyer RA, Cossart YE, et al. Human papillomavirus positivity predicts favourable outcome for squamous carcinoma of the tonsil. Int J Cancer. 2003;106(4):553–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.11261&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12845651&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184672300014&link_type=ISI) 

72. 72.Ndiaye C, Mena M, Alemany L, Arbyn M, Castellsagué X, Laporte L, et al. HPV DNA, E6/E7 mRNA, and p16INK4a detection in head and neck cancers: a systematic review and meta-analysis. Lancet Oncol. 2014;15(12):1319–31.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1470-2045(14)70471-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25439690&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

73. 73.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(2):467–75.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY2VicCI7czo1OiJyZXNpZCI7czo4OiIxNC8yLzQ2NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzEyLzIwMjAuMTAuMjYuMjAyMTMyMzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

74. 74.Amit M, Ilana K, Avraham SP, Binenbaum Y, Bachar G, Billan S, et al. Trends in human papillomavirus-related oropharyngeal cancer in Israel. Head Neck. 2016;38 Suppl 1:E274–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hed.23985&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25545260&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

75. 75.Dai M, Bao YP, Li N, Clifford GM, Vaccarella S, Snijders PJ, et al. Human papillomavirus infection in Shanxi Province, People’s Republic of China: a population-based study. Br J Cancer. 2006;95(1):96–101.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.bjc.6603208&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16773069&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238617900016&link_type=ISI) 

76. 76.Li LK, Dai M, Clifford GM, Yao WQ, Arslan A, Li N, et al. Human papillomavirus infection in Shenyang City, People’s Republic of China: A population-based study. Br J Cancer. 2006;95(11):1593–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.bjc.6603450&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17088908&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000242745800021&link_type=ISI) 

77. 77.Wu RF, Dai M, Qiao YL, Clifford GM, Liu ZH, Arslan A, et al. Human papillomavirus infection in women in Shenzhen City, People’s Republic of China, a population typical of recent Chinese urbanisation. Int J Cancer. 2007;121(6):1306–11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.22726&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17417776&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

78. 78.Lin H, Ma YY, Moh JS, Ou YC, Shen SY, ChangChien CC. High prevalence of genital human papillomavirus type 52 and 58 infection in women attending gynecologic practitioners in South Taiwan. Gynecol Oncol. 2006;101(1):40–5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ygyno.2005.09.028&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16256180&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

79. 79.Lai CH, Huang HJ, Hsueh S, Chao A, Lin CT, Huang SL, et al. Human papillomavirus genotype in cervical cancer: a population-based study. Int J Cancer. 2007;120(9):1999–2006.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.22538&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17266033&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

80. 80.Chen HC, You SL, Hsieh CY, Schiffman M, Lin CY, Pan MH, et al. Prevalence of genotype-specific human papillomavirus infection and cervical neoplasia in Taiwan: a community-based survey of 10,602 women. Int J Cancer. 2011;128(5):1192–203.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.25685&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20853317&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000286960300019&link_type=ISI) 

81. 81.Smith JS, Lindsay L, Hoots B, Keys J, Franceschi S, Winer R, et al. Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical lesions: a meta-analysis update. Int J Cancer. 2007;121(3):621–32.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.22527&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17405118&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000247407400020&link_type=ISI) 

82. 82.Jeng CJ, Phd. l, Ko ML, Ling QD, Shen J, Lin HW, et al. Prevalence of cervical human papillomavirus in Taiwanese women. Clin Invest Med. 2005;28(5):261–6.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16265998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

83. 83.Wang CH, Garvilles RG, Chen CY. Characterization of human papillomavirus infection in north Taiwan. J Med Virol. 2010;82(8):1416–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.21812&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20572075&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

84. 84.Chan PK, Cheung TH, Li WH, Yu MY, Chan MY, Yim SF, et al. Attribution of human papillomavirus types to cervical intraepithelial neoplasia and invasive cancers in Southern China. Int J Cancer. 2012;131(3):692–705.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.26404&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21976212&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000304440100036&link_type=ISI) 

85. 85.Chan PK, Ho WC, Chan MC, Wong MC, Yeung AC, Chor JS, et al. Meta-analysis on prevalence and attribution of human papillomavirus types 52 and 58 in cervical neoplasia worldwide. PLoS One. 2014;9(9):e107573.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0107573&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25229350&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

86. 86.de Sanjosé S, Brotons M, Pavón MA. The natural history of human papillomavirus infection. Best Pract Res Clin Obstet Gynaecol. 2018;47:2–13.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bpobgyn.2017.08.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

87. 87.Gillison ML, Broutian T, Pickard RK, Tong ZY, Xiao W, Kahle L, et al. Prevalence of oral HPV infection in the United States, 2009-2010. JAMA. 2012;307(7):693–703.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2012.101&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22282321&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000300278500026&link_type=ISI) 

88. 88.Steinau M, Saraiya M, Goodman MT, Peters ES, Watson M, Cleveland JL, et al. Human papillomavirus prevalence in oropharyngeal cancer before vaccine introduction, United States. Emerg Infect Dis. 2014;20(5):822–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2005.131311&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24751181&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

89. 89.Vawda N, Banerjee RN, Debenham BJ. Impact of Smoking on Outcomes of HPV-related Oropharyngeal Cancer Treated with Primary Radiation or Surgery. Int J Radiat Oncol Biol Phys. 2019;103(5):1125–31.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijrobp.2018.11.046&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30513378&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

90. 90.Morbini P, Alberizzi P, Ferrario G, Capello G, De Silvestri A, Pedrazzoli P, et al. The evolving landscape of human papillomavirus-related oropharyngeal squamous cell carcinoma at a single institution in Northern Italy. Acta Otorhinolaryngol Ital. 2019;39(1):9–17.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.14639/0392-100X-1905&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30936574&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

91. 91.Smith EM, Rubenstein LM, Haugen TH, Pawlita M, Turek LP. Complex etiology underlies risk and survival in head and neck cancer human papillomavirus, tobacco, and alcohol: a case for multifactor disease. J Oncol. 2012;2012:571862.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2012/571862&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22315596&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

92. 92.Münger K, Baldwin A, Edwards KM, Hayakawa H, Nguyen CL, Owens M, et al. Mechanisms of human papillomavirus-induced oncogenesis. J Virol. 2004;78(21):11451–60.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjExOiI3OC8yMS8xMTQ1MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzEyLzIwMjAuMTAuMjYuMjAyMTMyMzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

93. 93.Parfenov M, Pedamallu CS, Gehlenborg N, Freeman SS, Danilova L, Bristow CA, et al. Characterization of HPV and host genome interactions in primary head and neck cancers. Proc Natl Acad Sci U S A. 2014;111(43):15544–9.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTExLzQzLzE1NTQ0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMC4xMC4yNi4yMDIxMzIzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

94. 94.Bussu F, Sali M, Gallus R, Vellone VG, Zannoni GF, Autorino R, et al. HPV infection in squamous cell carcinomas arising from different mucosal sites of the head and neck region. Is p16 immunohistochemistry a reliable surrogate marker? Br J Cancer. 2013;108(5):1157–62.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/bjc.2013.55&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23403821&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000317038100021&link_type=ISI) 

95. 95.Rietbergen MM, Snijders PJ, Beekzada D, Braakhuis BJ, Brink A, Heideman DA, et al. Molecular characterization of p16-immunopositive but HPV DNA-negative oropharyngeal carcinomas. Int J Cancer. 2014;134(10):2366–72.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.28580&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24249631&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2020.10.26.20213231.atom) 

96. 96.Amin MB, Edge SB, American Joint Committee on Cancer. AJCC cancer staging manual. Eighth edition. ed. Switzerland: Springer; 2017. xvii, 1024 pages p.