Abstract

Background. Overall survival for HPV-associated oropharyngeal squamous cell carcinoma (OPSCC) has differed by sex, but little is known regarding cancer-specific outcomes. We assessed the independent association of sex with cancer-specific survival in patients with HPV-associated oropharyngeal squamous cell carcinoma (OPSCC). Methods. We identified 14,183 patients from the Surveillance, Epidemiology, and End Results (SEER) program with OPSCC and tumor HPV status. We used Kaplan–Meier methods to compare overall survival (OS) and OPSCC-specific survival (HNCSS) by patient sex and by tumor HPV status. We then separately fit multivariable survival and competing risk models evaluating the association of sex on these outcomes by tumor HPV status and stratified by the use of guideline-concordant OPSCC treatment. Results. A total of 10,210 persons with HPV-positive tumors (72.0%) and 3,973 with HPV-negative tumors (28.0%) were identified. A larger proportion of women had HPV-negative tumors (24.0%) versus HPV-positive tumors (13.2%; ). Women with HPV-positive tumors were less likely to receive guideline-concordant treatment compared to men. In unadjusted survival analyses, women did not differ in OS or HNCSS compared to men for HPV-positive tumors but had worse OS and HNCSS for HPV-negative tumors. After adjustment, men and women with HPV-positive OPSCC did not differ in OS or HNCSS. However, women with HPV-negative tumors faced worse overall survival (hazard ratio (HR) 1.15, 95% CI 1.02–1.29) that persisted even after stratifying for stage-appropriate treatment (HR 1.28, 95% CI 1.11–1.47). Conclusions. Women with HPV-positive OPSCC had similar survival outcomes compared to men, but those with HPV-negative tumors have worse overall and cancer-specific survival.

1. Introduction

For the past 4 decades, the incidence of oropharyngeal squamous cell carcinoma (OPSCC) has steadily increased, growing by 122% in men and 43% in women from 1999 to 2015. [14] This rise is partially attributable to the increase in human papillomavirus (HPV) infection, the most common sexually transmitted infection among men and women in the United States. [58] OPSCC is now the most prevalent HPV-associated cancer. [1] Although not all OPSCC cases are associated with HPV infection, the proportion of HPV-related tumors among all OPSCC cases is increasing. Recent studies estimate that now over 70% of OPSCC cases are positive for p16, a biomarker frequently used to identify HPV-associated tumors. [9].

Incidence and temporal trends related to HPV-associated OPSCC specifically have also differed. Recent estimates have shown that HPV-associated OPSCC has increased annually 0.8% for women and 2.7% for men. [1] Although OPSCC disproportionately affects men and Whites, it appears to be increasing among all sex and race groups. [5, 10, 11].

OPSCC p16 status has been strongly associated with treatment responsiveness and survival outcomes. [1216] However, outcomes of OPSCC in women, particularly those stratified by HPV status, have not been extensively studied and initial studies have conflicting conclusions regarding differences in survival. Although some studies have found that women have better overall survival with HPV-associated OPSCC, [17] others have found no difference. [1820] Still others have found worse outcomes for women. [18, 19].

Previous OPSCC outcome studies have been limited by sample size or lack of data on cause of death, limiting the analysis to overall survival only, and, therefore, are susceptible to bias related to noncancer causes of death related to comorbid illnesses and behavioral factors. In addition, large cancer registry studies have lacked information on HPV status for OPSCC cases. In this study, we used data from a contemporary cohort of OPSCC cases from a national representative cancer registry database to compare overall and cancer-specific survival by sex, according to tumor HPV status. Our hypothesis is that women with OPSCC will have better survival outcomes regardless of HPV status, in line with trends seen in other cancers overall.

2. Materials and Methods

2.1. Patient Selection

We performed a retrospective cohort study using data from the Surveillance, Epidemiology, and End Results (SEER) Head and Neck with HPV Status Database. The SEER program has collected clinicopathologic data on incident cancer cases from 18 population-based registries, which represent approximately 28% of the US population. The Head and Neck with HPV Status Database includes cases diagnosed between 2010 and 2016, for which available HPV data have been collected. [21] Tumor HPV status in this database was determined through submission per the SEER Collaborative Stage Data Collection System, version 02.02–02.05 schemas. This included any testing performed on surgical specimens such as HPV in situ hybridization (ISH), tissue PCR, ISH for E6/7 RNA, real-time PCR for E6/7 RNA, and p16 immunohistochemistry. Blood or serology testing was excluded. HPV status was determined using the applicable collaborative stage site-specific factor 10 for each disease site schema throughout the data collection period. [22].

We identified incident cases of OPSCC using the International Classification of Diseases for Oncology, Third Edition (ICD-O-3), topography codes C01.9–10.9, which includes the base of the tongue, lingual tonsil, soft palate, uvula, tonsillar fossa, tonsillar pillar, tonsil, vallecula, anterior surface of the epiglottis, lateral and posterior wall of the oropharynx, branchial cleft, and oropharynx, not otherwise specified. Cases that were not of squamous cell histology, as determined by ICD-O-3 (8050–8084), were excluded. Next, we identified and included cases with known American Joint Committee on Cancer, 7th edition, clinical stage, which we restaged to the 8th edition based on available T, N, and M information on the individual level. Cases with unknown HPV status were excluded (n = 11,662). A total of 14,183 patients were identified (Figure 1). We then classified tumors according to SEER HPV data as HPV positive (n = 10,210) or HPV negative (n = 3,973).

2.2. Study Variables

The SEER data were used to identify age, race/ethnicity, year of OPSCC diagnosis, median household income by county attribute, proportional educational attainment by county attribute, and insurance type (which we classified as uninsured, any Medicaid, insured, insured with no specifics, or unknown). Tumor characteristics (stage and T, N, and M values) and treatment type were ascertained from SEER on the individual level. We defined surgery as any definitive resection surgery performed, including excisional surgery of the oropharynx, the base of the tongue resection, palatectomy, uvulectomy, and/or tonsillectomy. Radiotherapy was coded for cases receiving beam radiation or a combination of a beam with implants or isotopes. Chemotherapy was included in the database and encoded as “yes” or “no/unknown.” Each case was then encoded with all treatments received, whether surgery, radiotherapy, or chemotherapy alone or in combination with other treatment methods. To determine whether the guideline-concordant treatment was administered by clinical stage, treatments were compared to National ComprehensiveCancer Network (NCCN) guidelines. [23] NCCN guidelines differ based on p16 status and clinical staging. An overview of NCCN guidelines can be found in Figure 2. To accommodate the categorical nature of the data available, we summarized the NCCN guidelines to whether surgery, chemotherapy, radiotherapy, or a combination of the three was recommended for each patient’s p16 status and clinical stage. Treatment for each case was considered to meet recommended guidelines if the patient received at least the recommended modality of treatment.

The primary outcome in this study was head and neck cancer-specific death. Any cause of death due to cancer in head and neck structures was encoded as a head and neck cancer-specific death. Our secondary outcome was death from any cause. Survival times for primary analyses were calculated from months of survival provided in SEER data, defined as the time from date of diagnosis to either death or end of follow-up.

2.3. Statistical Analysis

We tested for differences in baseline and clinical characteristics between women and men stratified by tumor HPV status using the t-test for normal continuous variables and the chi-square test for categorical variables. We then compared the use of cancer treatments and stage-specific guideline-concordant treatments by sex, also stratified by tumor HPV status, using the chi-square test. To further evaluate differences in cancer treatment by sex, we fit multivariable models assessing sex as a predictor of guideline-concordant treatment, adjusting for age, race, educational attainment and income by county attribute, insurance type, cancer stage, and marital status. We then fit unadjusted Cox regression models to compare overall and cancer-specific survival by sex. Adjusted models of overall and cancer-specific survival were then run to measure the association of sex with these outcomes, separately adjusting for race, sociodemographic factors (median county household income, proportional county education attainment, and insurance type), cancer stage, and treatment by tumor HPV status and stratified by the use of guideline-concordant OPSCC treatment. We separately compared outcomes for all HPV-positive and HPV-negative OPSCC cases, then, stratified by tumor stage, and then, further stratified by use of guideline-concordant treatment, using the aforementioned models. We also tested for interactions between sex and tumor stage. Last, to evaluate the effects of potential differences in competing (non-OPSCC-related) causes of death by sex, we fit similar multivariable Fine-Gray regression models, stratified by HPV tumor status, first for all cases, and then limited to patients receiving guideline-concordant treatment. Complete models are provided in Supplementary Tables. All analyses were performed in STATA Version 13 (STATA Corp., College Station, TX). STATA code for all analyses can be found online at this link: https://github.com/derkao/opscc-by-sex. The Mount Sinai Institutional Review Board determined that this study was exempt from human research.

3. Results

3.1. Baseline Characteristics

In total, 14,183 cases of OPSCC with known staging and HPV status were identified, of which 10,210 (72.0%) were HPV positive and 3,973 (28.0%) were HPV negative (Table 1). Women were more likely to have HPV-negative tumors than HPV-positive tumors (24.0% versus 13.2%; ). For HPV-positive OPSCCs, in comparison to men, women were older (), more likely to be non-Caucasian (), early stage (), single (), and more likely to be insured with Medicaid (). Among HPV-negative OPSCC cases, women were also older compared to men at diagnosis (), earlier stage (), and more likely to be single ().

3.2. Treatment by Stage

In univariate analyses, women with HPV-positive OPSCC were significantly less likely to receive guideline-concordant treatment (Table 2; ) compared to men. For HPV-negative tumors, a disparity of similar magnitude was seen; however, the comparison was not significant (). Differences in the use of guideline-concordant treatment for women with HPV-positive tumors appear to stem from less aggressive treatment for earlier-stage (I-III) cancers. In contrast, for HPV-negative tumors, women were less likely to receive guideline-concordant treatment for more advanced-stage cancers (stage III, , stage IV, ). After adjusting for sociodemographic factors, including race, median county household income, proportional county education attainment, insurance type, cancer stage, and marital status, sex was still a significant predictor of guideline-concordant treatment for HPV-positive tumors (; results not otherwise shown) but not HPV-negative tumors ().

3.3. Survival Outcomes

In unadjusted survival analyses (results not otherwise shown for unadjusted survival analyses), men and women with HPV-positive tumors did not differ in overall survival (OS) (HR 1.05, 95% CI 0.90–1.22) or head and neck cancer-specific survival (HNCSS) (HR 1.17, 95% CI 0.95–1.44), defined by death due to cancer in any head and neck structure. However, in HPV-negative cases, women had worse OS (HR 1.27, 95% CI 1.13–1.43) and HNCSS (HR 1.29, 95% CI 1.10–1.51). Worse cancer-specific survival persisted for women after accounting for possible competing risks (HR 1.25, 95% CI 1.07–1.47). Among patients who received guideline-concordant treatment, we found similar trends; women had no difference in outcomes for HPV-positive tumors but had significantly worse HNCSS for HPV-negative tumors even after accounting for competing risks (p = 0.02).

In adjusted analyses, men and women had similar overall survival with HPV-positive tumors overall and among the subset of cases that received guideline-concordant treatment (Table 3, Figure 3(a)). This was true for head and neck cancer-specific survival as well (Table 3, Figure 3(b)). In multivariable competing risk models, HNCSS was also not different by sex after adjustment for all HPV-positive cases or among those receiving guideline-concordant treatment.

Survival differences between women and men were seen with HPV-negative tumors. Adjusted survival analyses for the overall HPV-negative tumor cohort showed decreased OS for women (HR 1.15, 95% CI 1.02–1.29, Figure 3(c)). This effect was larger in the subset of cases that received guideline-concordant treatment (HR 1.28, 95% CI 1.11–1.47). Statistically significant differences were also seen in stage III with guideline-concordant treatment, stage IVA overall, and stage IVB with guideline-concordant treatment. In HNCSS analysis, survival differences were not statistically significant for the overall cohort of HPV-negative tumors (Figure 3(d)). However, HNCSS was worse in women receiving guideline-concordant treatment in all stages combined (HR 1.30, 95% CI 1.07–1.59); this persisted after accounting for competing risks (HR 1.26, 95% CI 1.03–1.55). These differences were most prominent among patients with stage IVA and IVB tumors. In the analyses testing for statistical interactions between patient sex and tumor stage, there were no statistically significant interactions in groups where there were underlying differences in outcomes (OS or HNCSS) by sex.

4. Discussion

In this population-based cohort study of HPV-associated OPSCC, we found no differences in survival outcomes by sex for HPV-positive tumors but worse survival for women for HPV-negative tumors. This provides evidence that contrasts previous findings that women may have a worse prognosis with HPV-associated OPSCC. For HPV-negative tumors, we found evidence of treatment disparities by sex, and survival differences were still pronounced in women who received guideline-concordant treatment. Treatment-related factors such as tolerability or responsiveness or possibly tumor behavior may, therefore, differ for women with HPV-negative OPSCC. Other factors may include preemptive surgery and the timing of treatment, as early surgery of preneoplastic lesions has been shown to be beneficial. [24].

We found differences in treatment patterns for women and men with OPSCC. For both HPV-positive and HPV-negative tumors, women were less likely to receive guideline-concordant treatment. A recent large cohort study (n = 884) compared to head and neck cancer treatment by sex and found that women received less intensive treatment, despite worse cancer-specific survival. [25] However, this study did not examine differences by HPV status. Additionally, a matched-pair analysis of 572 female and male patients, in which treatment was a matched variable, demonstrated no significant differences in survival outcomes between men and women. [26] There are few studies that have examined treatment differences between men and women stratified by HPV status. In an analysis of the National Cancer Database, patterns of treatment for HPV-positive and HPV-negative OPSCC differed for women versus men; however, proportions of untreated patients were similar. [19] It is, therefore, likely that women experience cancer treatment disparities compared to men for OPSCC, similar to what is observed for other tumor types. [2731] These disparities may be partially explained by a greater burden of comorbidities seen in women versus men with OPSCC, as demonstrated in population-based data. [32].

Prior comparisons of differences in prognosis for OPSCC by sex, particularly for those known to be HPV-associated, have been conflicting. In a two-site retrospective analysis of 239 cases, Yin et al. found that women with HPV-positive OPSCC had improved overall survival compared to men. [17] In contrast, three population-based studies have investigated overall survival by sex and tumor HPV status for OPSCC : Faraji et al. conducted a retrospective analysis using the National Cancer Database with cases from 2010 to 2015, Li et al. used the same database with cases from 2010 to 2014, and Razzaghi et al. used the National Program of Cancer Registries from the CDC. [1820] These larger, population-based studies all found overall survival trends consistent with our study for HPV-positive tumors: women with HPV-positive tumors had no difference compared to men.

Differing from our findings, in HPV-negative OPSCC Yin et al. found no difference in overall survival between men and women using two-site retrospective data. Studies from population-based cancer registry data, however, found that women fared worse than men, in agreement with this study. [1719] No previous study has had cause-of-death information, and however, our study shows that OPSCC-specific survival is no different for men and women for HPV-positive tumors but that differences in overall survival for HPV-negative tumors are partly driven by worse cancer-specific survival.

To our knowledge, this is the first study to examine survival between treatment modalities stratified by HPV status and sex. Feinstein et al. performed a retrospective analysis at a single regional veterans’ health center and found no difference in survival between treatment modalities in all OPSCC cases and HPV-positive OPSCC cases alone. [33] A systematic review from 2015 also demonstrated no statistically significant difference in hazard ratios between patients treated with primary surgery vs. primary radiation. [34] Both studies found better survival in HPV-positive tumors, but neither of these studies looked at sex-based differences. These differences highlight the importance of patient counseling and increased awareness among providers regarding sex-related differences in outcomes of patients undergoing treatment for OPSCC.

4.1. Limitations

Strengths of this study include the population-based national dataset that provided a large sample size and heterogeneity of the patient population, enabling generalization to the U.S. population. This database also provided HPV status data that were determined by robust methods. However, our findings must be viewed in light of several limitations. First, the lack of comorbidity data limits the interpretation of our data. Of note, cancer-specific survival data were available and were used in this study. Second, detailed treatment information on chemotherapy and radiotherapy is not available; thus, we cannot interpret our data with respect to predicting response to treatment. Third, our dataset did not include information on behavioral (e.g., smoking) or environmental exposures, which may have impacted cancer outcomes. Fourth, missing HPV data may be a source of selection bias, as previous studies have demonstrated disparities of HPV testing by demographics such as race and insurance. [35] Lastly, this study utilized a relatively large number of statistical tests; therefore, some values that were below the threshold for significance may represent results that occurred at random and not truly significant associations.

5. Conclusions

In this study of population-based cancer data, we found that outcomes for OPSCC did not differ for women compared to men for HPV-positive tumors. Women, however, did experience worse overall and cancer-specific survival for HPV-negative tumors, which was not fully explained by treatment disparities.

Data Availability

We used data from the Surveillance, Epidemiology, and End Results (SEER) Head and Neck with HPV Status Database, a national cancer registry that collects clinicopathologic data on incident cancer cases from 18 population-based registries. Access to data can be requested here: https://seer.cancer.gov/seerstat/databases/hpv/

Conflicts of Interest

Joshua M. Bauml reports a consulting/advisory role with Clovis Oncology, Bristol-Myers Squibb, Merck, AstraZeneca, Genentech, Celgene, Boehringer Ingelheim, Janssen, Takeda, Ayala, Regeneron, Inivata, Foundation Medicine, Novartis, and Guardant Health; Bauml also reports research support provided to his institution for other projects by Merck, Janssen, Carevive Systems, Novartis, Clovis Oncology, AstraZeneca, Takeda, and Bayer. All other authors declare no conflict of interest regarding the publication of this article.

Authors’ Contributions

DK and KS contributed to conception and design. All authors contributed to the acquisition, analysis, or interpretation of data. DK and KS contributed to the drafting of the manuscript. All authors contributed to the critical revision of the manuscript for important intellectual content. DK contributed to statistical analysis. DK and KS obtained funding. None contributed to administrative, technical, or material support. KS contributed to Supervision. RF, DM, TM, BM, and JB have given no additional contributions. All authors contributed to the final approval of the version to be published. DK had full access to the data and took responsibility for the integrity of the data and the accuracy of the data analysis.

Acknowledgments

This work was supported by the Tisch Cancer Institute Summer Scholars Grant.

Supplementary Materials

A PDF containing three tables that have the complete multivariable models for overall survival and head and neck cancer-specific survival for HPV-positive cases (Supplementary Table 1) and HPV-negative cases (Supplementary Table 2), and treatment breakdown by stage, HPV status, and sex (Supplementary Table 3) is included as a Supplementary File. (Supplementary Materials)