This paper is in the following e-collection/theme issue:

Surveillance Reports (137) Cancer Epidemiology, Cancer Surveillance and Infodemiology (77) Cancer Survivorship (258) Descriptive Epidemiology and Population Size Estimates (107)

Published on in Vol 9 (2023)

Preprints (earlier versions) of this paper are available at https://preprints.jmir.org/preprint/47579, first published .
Cause of Death in Patients with Oropharyngeal Carcinoma by Human Papillomavirus Status: Comparative Data Analysis

Cause of Death in Patients with Oropharyngeal Carcinoma by Human Papillomavirus Status: Comparative Data Analysis

Cause of Death in Patients with Oropharyngeal Carcinoma by Human Papillomavirus Status: Comparative Data Analysis

Authors of this article:

Dong-Dong Zhang1 Author Orcid Image ;   Min Lei1 Author Orcid Image ;   Yue Wang1 Author Orcid Image ;   Pei-Ji Zeng1 Author Orcid Image ;   Yong-Jun Hong1 Author Orcid Image ;   Cheng-Fu Cai1, 2, 3, 4, 5, 6 Author Orcid Image
Article Authors Cited by (1) Tweetations (2) Metrics

    Original Paper

    1Department of Otolaryngology-Head and Neck Surgery, School of Medicine, Zhongshan Hospital, Xiamen University, Xiamen, China

    2School of Medicine, Xiamen University, Xiamen, China

    3Xiamen Key Laboratory of Otolaryngology Head and Neck Surgery, Xiamen, Xiamen, China

    4Department of Otorhinolaryngology, Head and Neck Surgery, Xiamen University, Xiamen, China

    5Department of Otolaryngology Head and Neck Surgery, Fujian Medical University, Xiamen, China

    6College of Otorhinolaryngology Head and Neck Surgery, Xiamen Haicang Hospital, Xiamen, China

    Corresponding Author:

    Cheng-Fu Cai, MD

    Department of Otolaryngology-Head and Neck Surgery

    School of Medicine, Zhongshan Hospital

    Xiamen University

    Hubin South Road 201-209

    Xiamen, 361003

    China

    Phone: 86 13616011099

    Email: caichengfuxmu@126.com


    Background: The incidence of oropharyngeal squamous cell carcinomas (OPSCC) has increased in recent decades, and human papillomavirus (HPV) infection is the main cause of OPSCC. The data regarding causes of death (CODs) are vitally important in informing follow-up strategies and revising treatment strategies to deal with any possible preventable treatment-related COD. However, limited studies have assessed the competing COD by HPV status in patients with OPSCC.

    Objective: We aimed to analyze the distribution of the competing COD according to HPV status in OPSCC.

    Methods: We retrospectively included stage I-IVB patients with OPSCC from the Surveillance, Epidemiology, and End Results database between 2010 and 2015. The association between HPV status and head and neck cancer–specific mortality (HNCSM), second primary cancer mortality (SPCM), and noncancer-caused mortality (NCCM) were analyzed. The chi-square test, Kaplan-Meier analysis, and Fine and Gray model were used for statistical analysis.

    Results: We included 5852 patients in this study and 73.2% (n=4283) of them had HPV-related tumors. A total of 1537 (26.3%) patients died, including 789 (51.3%), 333 (21.7%), and 415 (27%) patients who died from head and neck cancer, second cancer, and noncancer causes, respectively. The 5-year HNCSM, SPCM, NCCM, and overall mortality were 14.7%, 6.5%, 7.7%, and 26.4%, respectively. Those with HPV-positive disease had a lower cumulative incidence of HNCSM (subdistribution hazard ratio [sHR] 0.362, 95% CI 0.315-0.417; P<.001), SPCM (sHR 0.400, 95% CI 0.321-0.496; P<.001), and NCCM (sHR 0.460, 95% CI 0.378-0.560; P<.001) than those with HPV-negative disease. The 5-year risk of HNCSM was 26.9% and 10.7% in those with HPV-negative and HPV-positive disease, respectively (P<.001). The 5-year risk of SPCM was 12.4% and 4.6% in those with HPV-negative and HPV-positive disease, respectively (P<.001). The 5-year risk of NCCM of death was 13.7% and 5.8% in those with HPV-negative and HPV-positive disease, respectively (P<.001). Using the Fine and Gray competing-risks model, our results show that those with HPV-negative tumors had a significantly higher risk of HNCSM (P<.001), SPCM (P<.001), and NCCM (P<.001) than those with HPV-negative tumors.

    Conclusions: HPV-positive OPSCC has a lower NCSM, SPCM, and NCCM as compared to those with HPV-negative OPSCC. HPV positivity is a favorable prognostic factor in the context of overcoming cancer as well as in terms of reducing the risk of other CODs in OPSCC. Our finding supports the need to tailor patient follow-up based on the HPV status of patients with OPSCC.

    JMIR Public Health Surveill 2023;9:e47579

    doi:10.2196/47579

    Keywords

    oropharyngeal carcinomas; human papillomavirus; causes of death; survival; competing mortality; human papillomavirus; comparative data; analysis; cell death; carcinoma; oropharyngeal carcinoma; surveillance; database; human papillomavirus-related tumor 


    Human papillomavirus (HPV) is the most prevalent sexually transmitted infection globally and is associated with both anogenital and oropharyngeal cancers. It is estimated that over 80% of individuals will be affected by HPV at some point in their lifetime [1]. Like other types of head and neck squamous cell carcinomas (HNSCCs), oropharyngeal squamous cell carcinoma (OPSCC) has been historically associated with alcohol and tobacco consumption. However, infection with carcinogenic HPV infection has emerged as an important risk factor that has driven an increase in the incidence of OPSCC. In recent decades, OPSCC has been one of the malignant tumors with rapidly rising incidences in Western countries [2-5]. In Italy, the incidence of HPV-related OPSCC increased from 16.7% in 2000-2006 to 46.1% in 2013-2018 [6]. Moreover, HPV now accounts for 71% and 51.8% of all OPSCCs in the United States and the United Kingdom, respectively [7-9]. The etiological relation with OPSCC by HPV status was also found in the Asian population [10,11]. The existing evidence shows that HPV status impacts the prognosis and treatment decision-making of OPSCC [12,13]. Patients with HPV-positive disease had a significantly better prognosis than those with HPV-negative disease [12,13]. Moreover, several studies also explored the treatment de-escalation for HPV-positive OPSCC [14,15]. These differences have a direct effect on clinical trials that will determine treatment therapeutic standards in the future clinical practice of OPSCC.

    As survivorship from OPSCC continues to increase in recent decades [16-18], the main participants in the medical field should identify the patients with OPSCC with the highest risk of dying as well as the risk of their specific cause of death (COD). In addition to HPV, smoking and alcohol use are also closely related to the development of OPSCC, which may cause competing mortality in this population. In 2 studies from the Netherlands and the US cancer registry showed that the survival outcomes of head and neck cancer (HNC) were inferior to that of a healthy population of similar gender and age [19,20]. Smoking and alcohol use are the causative factors for cardiopulmonary diseases and second primary tumors. A previous Surveillance, Epidemiology, and End Result (SEER) study showed that life expectancy for patients with HNC was significantly shorter than expected after excluding the COD related to HNC [21]. The other COD may contribute to the differences in the survival of patients with HNC [22,23]. This information is vitally important in informing follow-up strategies and revising treatment strategies to deal with any possible preventable treatment-related COD.

    Because HPV-positive patients with OPSCC are usually younger and consume less alcohol and tobacco, the risk of cardiopulmonary disease and secondary primary tumor in this population may be lower [24,25]. The differences in survival rates for patients with HNSCC can primarily be attributed to an unequal burden of comorbidities associated with tobacco and alcohol use, treatment-related sequelae, and a high risk for the development of second primary cancer. Moreover, the surveillance for HPV-related tumors is similar to HPV-negative tumors in the current guidelines. Given the rapid rise in the incidence of potentially HPV-related OPSCC, it is crucial to provide a more comprehensive description of the COD in OPSCC. Understanding the COD by different HPV status in OPSCC could bring even more positive impact on public health, potentially impacting the follow-up strategies and revising treatment strategies while reducing treatment-related morbidity. However, limited studies have assessed the competing COD by HPV status in patients with OPSCC [26]. The objective of this study was to analyze the differences in the competing COD depending on HPV status in OPSCC.


    Patient Selection

    The National Cancer Institute SEER database collects population-level data on approximately 28% of the US population from 18 registries with the morbidity, mortality, and disease status of millions of patients with cancer [27]. We used the SEER databases to study patients diagnosed with OPSCC between 2010 and 2015 (International Classification of Disease for Oncology ICD-O C019, C090, C091, C098-C104, C108, C109, C142; histology codes 8073-8079 and 8083-8084). Registrars for the SEER database were required to enter HPV status beginning in 2010 and we have obtained HPV-related variables through the use of a data user agreement. The data of the retrospective cohort in the present study were obtained from the SEER HPV database, which includes the demographic data, clinicopathological characteristics, treatment, HPV status, and follow-up of patients with HNSCC [27]. We included patients who met the following criterion: (1) diagnosed between 2010 and 2015; (2) diagnosed with stage I-IVB OPSCC according to the 7th American Joint Committee on Cancer (AJCC) staging system; (3) available data regarding HPV status, race, tumor grade, AJCC staging, tumor stage, nodal stage, and metastatic stage; and (4) available data on surgery, radiotherapy, and chemotherapy. Patients with unknown HPV status, second or more malignant tumors, distant stage, and treated with radioactive implants or radioisotopes were excluded. The patient selection flowchart is listed in Figure 1.

    ‎
    Figure 1. Flow diagram of this study’s cohort. HPV: human papillomavirus; M: metastatic stage; N: nodal stage; NOS: not otherwise specified; NX: unknown nodal stage; RT: radiotherapy; T: tumor stage; TX: unknown tumor stage.

    Ethics Approval

    The SEER database is a public database, and we have obtained permission to access the database (approval number 18275-Nov2021). The present study did not require ethical approval by the institutional review board of the Zhongshan Hospital, Xiamen University, because of the nonidentifying information in the SEER database. The requirement for informed consent was waived because all data used in this study were deidentified.

    Data Collection

    The following variables were included: age, gender, race, tumor grade, HPV status, AJCC 7th staging, and the receipt of treatment. HPV status was determined by p16 immunohistochemistry, tissue polymerase chain reaction, in situ hybridization, real-time polymerase chain reaction for E6/7 RNA, or in situ hybridization for E6/7 RNA [28]. The primary outcomes were HNC-specific mortality (HNCSM), second primary cancer mortality (SPCM), and noncancer-caused mortality (NCCM). Any COD due to HNC was encoded as HNCSM. Events for SPCM were death from the second primary cancer. Events for NCCM were death from noncancer diseases. Our secondary outcome was overall mortality (OM). OM was defined as death from all causes.

    Statistical Analysis

    Categorical variables were represented as absolute values with percentages, and the chi-square test was performed for the comparison between patients who were HPV-negative and HPV-positive. Kaplan-Meier analysis was used to generate the survival curves to compare differences in survival between groups over time, and the equality of these curves was compared by the log-rank test. Due to the competitive relationship among multiple outcomes in clinical survival data, the traditional Kaplan-Meier analysis will misjudge the cumulative mortality rate [29]. Therefore, competing-risks analyses using the Fine and Gray model were used to estimate the effects of covariates on HNCSM, SPCM, and NCCM. We conducted stratum-specific analyses to estimate the effects of HPV status on HNCSM, SPCM, and NCCM in accordance with gender, age, race, AJCC staging, grade, surgery, radiotherapy, and chemotherapy. All statistical analyses were conducted using SPSS (version 22.0; IBM Corp) and Stata/SE (version 14; StataCorp). A P value <.05 was considered statistically significant.


    Patient Baseline Characteristics

    From 2010 to 2015, there were 14,688 cases of HNC recorded with the HPV status in the SEER HPV database; 9496 of these cases were oropharyngeal cancer and 5852 of these cases were stage I-IVB OPSCC and had treatment data available. In this cohort (n=5852; Table 1), 4283 patients (73.2%) had HPV-positive OPSCC and 1569 patients (26.8%) had HPV-negative OPSCC. Of these patients, 84.1% (n=4922), 89.2% (n=5220), and 69.9% (n=4088) were male, White, and had stage T4 disease, respectively. Regarding treatment, 89.1% (n=5216), 71.3% (n=4174), and 55% (n=2626) of patients received radiotherapy, chemotherapy, and surgery, respectively.

    Table 1. Patient baseline characteristics.
    VariablesPatients, nHPVa-negative, n (%)HPV-positive, n (%)P value
    Gender<.001

    		Male49221213 (77.3)3709 (86.6)

    		Female930356 (22.7)574 (13.4)
    Age (years)<.001

    		<50840207 (13.2)633 (14.8)

    		50-643425847 (54)2578 (60.2)

    		≥651587515 (32.8)1072 (25)
    Race<.001

    		White52201288 (82.1)3932 (91.8)

    		Black414212 (13.5)202 (4.7)

    		Other21869 (4.4)149 (3.5)
    Grade<.001

    		Well differentiated233113 (7.2)120 (2.8)

    		Moderately differentiated2240754 (48.1)486 (34.7)

    		Poorly or undifferentiated3379702 (44.7)2677 (62.5)
    AJCCb stage<.001

    		I21893 (5.9)125 (2.9)

    		II384108 (6.9)276 (6.4)

    		III1162342 (21.8)820 (19.1)

    		IVA and IVB40881026 (65.4)3062 (71.5)
    Tc stage<.001

    		T11630360 (22.9)1270 (29.7)

    		T22322527 (33.6)1795 (41.9)

    		T31057342 (21.8)715 (16.7)

    		T4843340 (21.7)503 (11.7)
    Nd stage<.001

    		N0935338 (21.5)597 (13.9)

    		N11096326 (20.8)770 (18)

    		N23562836 (53.3)2726 (63.6)

    		N325969 (4.4)190 (4.4)
    Surgery<.001

    		No3226964 (61.4)2262 (52.8)

    		Yes2626605 (38.6)2021 (47.2)
    Radiotherapy<.001

    		No636228 (14.5)408 (9.5)

    		Yes52161341 (85.5)3875 (90.5)
    Chemotherapy.02

    		No1678487 (31)1191 (27.8)

    		Yes41741082 (69)3092 (72.2)
    Causes of death (N=1537).50

    		Head and neck cancer789378（52.6）411 (50.2)

    		Second cancer333156 (21.7)177 (21.6)

    		Noncancer cause disease415184 (25.6)231 (28.2)

    aHPV: human papillomavirus.

    bAJCC: American Joint Committee on Cancer.

    cT: tumor stage.

    dN: nodal stage.

    Patients with HPV-positive disease were more likely to be male, aged 50-64 years, White, have higher tumor grade, have early tumor stage, have advanced nodal stage, and have stage IV disease (all P<.001). Moreover, those with HPV-positive disease were more likely to receive surgery (P<.001) and radiotherapy (P=.015; Table 1). Similar distributions in HNCSM (P=.36), SPCM (P=.96), and NCCM (P=.26) were found according to HPV status (Figure 2).

    ‎
    Figure 2. The distributions in (A) head and neck cancer–specific mortality, (B) second primary cancer mortality, and (C) noncancer-caused mortality according to human papillomavirus status. HPV: human papillomavirus.

    Survival

    With a median follow-up time of 52 months (range, 0-107 months), a total of 1537 (26.3%) patients died. There were 789 (51.3%), 333 (21.7%), and 415 (27%) patients who died from HNC, second cancer, and noncancer causes, respectively. The survival curves have been listed in Figure 3. The 5-year HNCSM, SPCM, NCCM, and OM were 14.7%, 6.5%, 7.7%, and 26.4%, respectively.

    ‎
    Figure 3. The curves of (A) head and neck cancer–specific mortality, (B) second primary cancer mortality, (C) noncancer-caused mortality, and (D) overall mortality in the entire cohort. HNCSM: head and neck cancer–specific mortality; NCCM: noncancer-caused mortality; OM: overall mortality; SPCM: second primary cancer mortality.

    Competing Mortality

    The cumulative incidence of competing CODs based on HPV status regarding HNCSM, SPCM, and NCCM are listed in Figure 4. Those with HPV-positive disease had a lower cumulative incidence of HNCSM (subdistribution hazard ratio [sHR] 0.362, 95% CI 0.315-0.417; P<.001; Figure 4A), SPCM (sHR 0.400; 95% CI 0.321-0.496; P<.001; Figure 4B), and NCCM (sHR 0.460; 95% CI 0.378-0.560; P<.001; Figure 4C) than those with HPV-negative disease. Using Kaplan-Meier analysis, the results showed that patients who were HPV-negative were more likely to gain experience from second cancer and noncancer causes. The 5-year risk of HNCSM was 26.9% and 10.7% in those with HPV-negative and HPV-positive disease, respectively (P<.001, Figure 5A). The 5-year risk of SPCM was 12.4% and 4.6% in those with HPV-negative and HPV-positive disease, respectively (P<.001, Figure 5B). The 5-year risk of NCCM of death was 13.7% and 5.8% in those with HPV-negative and HPV-positive disease, respectively (P<.001, Figure 5C). The 5-year OM was 26.1% and 10.7% in those with HPV-negative and HPV-positive, respectively (P<.001; Figure 5D).

    ‎
    Figure 4. The cumulative incidence of (A) head and neck cancer–specific mortality, (B) second primary cancer mortality, and (C) noncancer-caused mortality by human papillomavirus status using competing-risks regression. HNCSM: head and neck cancer–specific mortality; HPV: human papillomavirus; NCCM: noncancer-caused mortality; SPCM: second primary cancer mortality.
    ‎
    Figure 5. The effect of human papillomavirus status on (A) head and neck cancer–specific survival, (B) overall survival, (C) second primary cancer survival, and (D) noncancer-caused survival. HNCSM: head and neck cancer–specific mortality; HPV: human papillomavirus; NCCM: noncancer-caused mortality; OM: overall mortality; SPCM: second primary cancer mortality.

    Fine and Gray Competing-Risks Model

    Table 2 shows the sHR and 95% CI for selected covariates using competing-risks regression for different COD. The results of the multivariate analysis showed that the variables related to HNCSM were female, older age, Black, advanced stage, HPV-negative, receipt of surgery, and receipt of radiotherapy (all P<.05). Advanced stage, HPV-negative, receipt of surgery, and receipt of radiotherapy were the independent risk factors associated with SPCM (all P<.05). Older age, Black race, HPV-negative, receipt of surgery, and receipt of chemotherapy were the independent risk factors associated with NCCM (all P<.05). Stratum-specific analyses also showed that patients with HPV-positive tumor had lower HNCSM, SPCM, and NCCM than those with HPV-negative tumors after stratification by different demographic, clinicopathologic, and treatment characteristics (Table S1 in Multimedia Appendix 1).

    Table 2. Multivariate prognostic analyses using the competing-risks model.
    VariablesHNCa-specific mortalitySecond primary cancer mortalityNoncancer-caused mortality
    sHRb95% CIP valuesHR95% CIP valuesHR95% CIP value
    Gender

    		Male (reference)1N/AcN/A1N/AN/A1N/AN/A

    		Female1.2181.012-1.467.040.7290.528-1.005.050.9930.771-1.279.96
    Age (years)

    		<50 (reference)1N/AN/A1N/AN/A1N/AN/A

    		50-641.4971.160-1.932.0021.2310.872-1.737.241.2510.861-1.819.24

    		≥652.1761.665-2.844<.0011.260.868-1.830.223.1392.164-4.554<.001
    Race

    		White (reference)1N/AN/A1N/AN/A1N/AN/A

    		Black1.5251.215-1.913<.0011.1670.811-1.677.411.5481.131-2.120.006

    		Other1.0490.721-1.525.080.9310.523-1.658.811.2720.787-2.055.33
    AJCCd stage

    		I (reference)1N/AN/A1N/AN/A1N/AN/A

    		II2.1471.117-4.126.021.1240.435-2.901.810.9860.575-1.691.96

    		III2.9691.637-5.384<.0010.4030.615-3.201.420.8950.545-1.469.66

    		IVA and IVB4.5152.534-8.043<.0012.5831.168-5.713.020.8750.542-1.412.59
    Grade

    		Well differentiated (reference)1N/AN/A1N/AN/A1N/AN/A

    		Moderately differentiated0.7970.581-1.094.161.0830.638-1.839.771.0590.669-1.676.81

    		Poorly or undifferentiated0.5900.430-0.810.0010.8140.480-1.381.450.8470.532-1.351.49
    HPVe status

    		Negative (reference)1N/AN/A1N/AN/A1N/AN/A

    		Positive0.4390.378-0.508<.0010.4370.349-0.547<.0010.5840.472-0.721<.001
    Surgery

    		No (reference)1N/AN/A1N/AN/A1N/AN/A

    		Yes0.4630.389-0.552<.0010.5240.400-0.689<.0010.5720.456-0.716<.001
    Radiotherapy

    		No (reference)1N/AN/A1N/AN/A1N/AN/A

    		Yes0.4070.318-0.524<.0010.3870.267-0.561<.0010.9320.666-1.304.68
    Chemotherapy

    		No (reference)1N/AN/A1N/AN/A1N/AN/A

    		Yes0.9870.801-1.215.901.1520.834-1.592.390.6880.525-0.877.003

    aHNC: head and neck cancer.

    bsHR: subdistribution hazard ratio.

    cN/A: not applicable.

    dAJCC: American Joint Committee on Cancer.

    eHPV: human papillomavirus.


    Principal Findings

    In this study, we investigated the competing mortality in OPSCC by HPV status using a large cohort. Our results indicated that patients with HPV-negative tumors had a significantly higher risk of HNCSM, SPCM, and NCCM than those with HPV-positive tumors. The main finding of our study was that HPV-positive disease was an independent factor for better survival outcomes across all COD.

    HPV status has been confirmed to be an important prognostic and predictive biomarker in OPSCC [12,13]. In recent years, several studies also have explored the relationship between HPV status and the second tumor in OPSCC. Gorphe et al [30] analyzed 888 patients with OPSCC between 2011 and 2020, they found that the 5-year competing risks of a second cancer were 16.1% and 49.9% in those patients who were HPV-positive and HPV-negative, respectively (P<.001). In addition, a study from Bosshart et al [31] that included 91 patients also found that HPV-positive had a lower risk of second cancer than those with HPV-negative tumors (14% vs 38.2%). Since the time of COD was only recorded in the SEER database, we could not obtain the exact time of diagnosis of the second primary tumor. However, Culié et al [32] included 340 patients with OPSCC from 14 French hospitals, and they found that patients who were HPV-positive developed low second tumors than patients who were HPV-negative (3.3% vs 11.7%, P=.02). However, the COD was similar in the 2 groups.

    In our study, we found HNC-related death was the main COD in OPSCC (14.7%), followed by noncancer causes (7.7%), and second cancer (6.5%). Alcohol and tobacco use have been confirmed as risk factors for HNSCC. HPV-related OPSCC cannot be considered in the same field of cancerization as HPV-negative OPSCC and other HNSCCs [2-5]. However, we observed a quite similar distribution of CODs between the 2 groups, with primary cancer as the major COD. Our study suggested that the overall survival pattern of patients with OPSCC depends not only on their clinicopathological characteristics, but also significantly varies based on their age and the presence of comorbidities.

    Several previous studies have assessed CODs in patients with HNSCC. They found that age, gender, race or ethnicity, presence of comorbidities, BMI, or type of treatment were the risk factors in the competing COD in patients with HNSCC [22,33-37]. In our study, we analyzed whether the HPV status would influence the COD of patients with OPSCC. A population-based cohort from Denmark was the first to investigate COD by HPV status in 1521 patients with OPSCC [26]. There were 54.3% of patients were HPV-positive and they found that HPV-positive was associated with better survival across all CODs [26]. Lop et al [38] included 423 patients with OPSCC and 12.5% were HPV-positive. They found that patients with HPV-positive disease had better OPSCC-specific survival (P<.001) and second primary cancer survival (P<.001) than those with HPV-negative tumors. However, noncancer-related causes survival was similar between the 2 groups (P=.13). The difference in patient characteristics may affect the distribution of CODs in the previous studies. Those with HPV-positive disease had a significantly lower cumulative incidence of death from all CODs than those with HPV-negative disease. Therefore, the inferior survival in patients with HPV-negative tumors was reflected in the higher risk of HNCSM, SPCM, and NCCM. Patients with HPV-negative OPSCC may benefit from screening for the risk to die of noncancer death or early detection of second cancers. Our study suggests that OPSCC should undergo more different effective therapies, more prevention methods to avoid noncancer deaths, and more screening tools for early detection of the second tumor according to different HPV status. They may benefit from increased research funding.

    There were several studies indicating that those with HPV-positive OPSCC have a favorable survival than those with HPV-negative disease regardless of the treatment received [39,40]. In our study, 89.8% of patients were locally advanced-stage OPSCC, we also found that the receipt of treatment was associated with all CODs. The probability of treatment compliance was not recorded in the SEER database. Hess et al [41] found that treatment compliance in patients who were HPV-negative was significantly lower than in those with patients who were HPV-positive. Therefore, treatment compliance by HPV status might also have a negative effect on cancer-specific survival in OPSCC.

    Screening for secondary cancer has an important effect in reducing mortality, especially for the HPV-negative population with a higher level of tobacco use [42]. Our data indicate that there were 51.3% and 21.7% patients who died from HNC and second cancer. Therefore, more frequent follow-ups should be performed to increase the detection of recurrence and secondary cancers. Moreover, our study shows an increased risk of SPCM, and NCCM in patients who were HPV-positive than in those with HPV-negative tumors. Our finding is important because it adds to the current knowledge regarding the differences in CODs among HPV status and supports the need to individualize follow-ups according to the HPV status.

    In our study, 415 patients died of noncancer-related diseases, and 127 (30.6%) died from heart disease. In our study, patients with HPV-positive OPSCC tend to be younger, White, and male. Previous studies also showed that women of the Black race had a higher risk of cardiovascular disease in the general population and were cancer survivors [43,44]. According to the differences regarding epidemiological characteristics by HPV status may affect the burden of comorbidity, thus affecting the NCCM.

    HPV-related tumors have a significant burden on health insurance [45]. While the HPV vaccination has been shown to prevent cervical cancer, HPV vaccination also has the potential to prevent a significant proportion of HPV-related OPSCC. We should note that 84.1% of patients in our study were male, and 85.6% were aged 50 years or older. The vaccine only prevents new HPV infections and does not treat existing infections [46]. HPV vaccination is generally not recommended for those older than age 26 years who are unvaccinated, recent guidelines have a permissive recommendation for HPV vaccination in unvaccinated adults aged 27-45 years [46]. Additionally, 47.7% of patients continued to smoke after their OPSCC diagnosis [47]. Therefore, current smoking status should be monitored throughout the treatment course of OPSCC, and smoking-cessation assistance should be provided as necessary, especially for patients with negative HPV status and a history of smoking. Therefore, further studies are needed to verify the efficacy of HPV vaccination in reducing the incidence of OPSCC and to investigate the effect of smoking cessation on various COD in patients with OPSCC.

    We acknowledge certain limitations in our study. First, retrospective studies have inherent weaknesses of selection bias. Second, the lack of data regarding other sexually transmitted infections limits the interpretation of our data because behaviors that put people at risk for HPV also increase their risk for other sexually transmitted infections. Third, the lack of comorbidity data also limits the interpretation of our data. Of note, HNCSM, SPCM, and NCCM data were available and were used in this study. Fourth, the consumption of smoking and alcohol was not recorded in the SEER database. The high rates of HPV-positive OPSCC in our study could potentially be attributed to the low prevalence of smoking habits and differences in sexual behavior when compared with other geographic regions with a low incidence of HPV-positive OPSCC [37,48]. This difference in the incidence of HPV-positive tumors could limit the generalization of our findings to regions with a low incidence of HPV-positive tumors and high levels of tobacco and alcohol use. Moreover, the median follow-up time was 52 months in our study and second cancer or noncancer death is likely to develop after this delay. Finally, the SEER does not have specific records on HPV status detection methods. However, a previous study showed that the overall HPV status (positive or negative) had high degrees of accuracy in the SEER database [49]. Further research addressing these limitations could provide a more comprehensive understanding of the relationship between HPV status and competing mortality in patients with OPSCC.

    Conclusions

    In conclusion, our study suggests that HPV-positive OPSCC has a lower NCSM, SPCM, and NCCM as compared to those with HPV-negative OPSCC. HPV positivity is a favorable prognostic factor in the context of overcoming cancer as well as in terms of reducing the risk of other CODs in OPSCC. These disparities have direct implications on clinical trials aiming to define future therapeutic standards in clinical practice, including treatment intensification or de-escalation according to response rates, investigation of adjuvant treatment after complete response, or follow-up strategies involving imaging modalities and duration. Moreover, medical interventions tailored to the intercurrent CODs may independently provide survival benefits in addition to cancer treatments.

    Acknowledgments

    This project was supported by the Innovation Project in Fujian Province (2019-CXB-29), the Natural Science Foundation of Fujian Science and Technology Department (2020J02060), the Key Medical and Health Project of Xiamen Science and Technology Bureau (3502Z20204009), the Natural Science Foundation of Fujian Province of China (2022J05006), and the National Natural Science Foundation of China (32070761).

    Data Availability

    The data set used in the data analysis is publicly available.

    Authors' Contributions

    DDZ and ML drafted the paper. CFC acquired the data sets. CFC conceived this study. YW, PJZ, and YJH conducted the statistical analyses. DDZ, ML, and CFC participated in this study’s design. All authors read and approved the final paper.

    Conflicts of Interest

    None declared

    Multimedia Appendix 1

    Table S1. Stratum-specific analyses to determine the effect of HPV status according to different demography and clinicopathological characteristics using the Fine and Gray competing-risks model.

    DOCX File , 27 KB
    1. Scott-Wittenborn N, Fakhry C. Epidemiology of HPV related malignancies. Semin Radiat Oncol. 2021;31(4):286-296. [FREE Full text] [CrossRef] [Medline]
    2. Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM. Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer. 2010;127(12):2893-2917. [FREE Full text] [CrossRef] [Medline]
    3. Blomberg M, Nielsen A, Munk C, Kjaer SK. Trends in head and neck cancer incidence in Denmark, 1978-2007: focus on human papillomavirus associated sites. Int J Cancer. 2011;129(3):733-741. [FREE Full text] [CrossRef] [Medline]
    4. Carlander ALF, Larsen CG, Jensen DH, Garnæs E, Kiss K, Andersen L, et al. Continuing rise in oropharyngeal cancer in a high HPV prevalence area: a Danish population-based study from 2011 to 2014. Eur J Cancer. 2017;70:75-82. [CrossRef] [Medline]
    5. Tardini E, Zhang X, Canahuate G, Wentzel A, Mohamed ASR, Van Dijk L, et al. Optimal treatment selection in sequential systemic and locoregional therapy of oropharyngeal squamous carcinomas: deep Q-learning with a patient-physician digital twin dyad. J Med Internet Res. 2022;24(4):e29455. [FREE Full text] [CrossRef] [Medline]
    6. Del Mistro A, Frayle H, Menegaldo A, Favaretto N, Gori S, Nicolai P, et al. Age-independent increasing prevalence of human papillomavirus-driven oropharyngeal carcinomas in North-East Italy. Sci Rep. 2020;10(1):9320. [FREE Full text] [CrossRef] [Medline]
    7. 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-4301. [FREE Full text] [CrossRef] [Medline]
    8. Senkomago V, Henley SJ, Thomas CC, Mix JM, Markowitz LE, Saraiya M. Human papillomavirus-attributable cancers—United States, 2012-2016. MMWR Morb Mortal Wkly Rep. 2019;68(33):724-728. [FREE Full text] [CrossRef] [Medline]
    9. Schache AG, Powell NG, Cuschieri KS, Robinson M, Leary S, Mehanna H, et al. HPV-related oropharynx cancer in the United Kingdom: an evolution in the understanding of disease etiology. Cancer Res. 2016;76(22):6598-6606. [FREE Full text] [CrossRef] [Medline]
    10. 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-417. [FREE Full text] [CrossRef] [Medline]
    11. Chor JSY, Vlantis AC, Chow TL, Fung SC, Ng FYH, Lau CH, et al. The role of human papillomavirus in head and neck squamous cell carcinoma: a case control study on a Southern Chinese population. J Med Virol. 2016;88(5):877-887. [CrossRef] [Medline]
    12. Wittekindt C, Klussmann JP. Tumor staging and HPV-related oropharyngeal cancer. Recent Results Cancer Res. 2017;206:123-133. [CrossRef] [Medline]
    13. Wagner S, Reder H, Sharma SJ, Würdemann N, Wittekindt C, Klußmann JP. HPV-driven oropharyngeal carcinoma - incidence, trends, diagnosis and therapy. Urologe. 2018;57(12):1457-1463. Das HPV-getriebene Oropharynxkarzinom – Inzidenz, Trends, Diagnose und Therapie. [CrossRef] [Medline]
    14. Petrelli F, Luciani A, Ghidini A, Cherri S, Gamba P, Maddalo M, et al. Treatment de-escalation for HPV+ oropharyngeal cancer: a systematic review and meta-analysis. Head Neck. 2022;44(5):1255-1266. [CrossRef] [Medline]
    15. Patel RR, Ludmir EB, Augustyn A, Zaorsky NG, Lehrer EJ, Ryali R, et al. De-intensification of therapy in human papillomavirus associated oropharyngeal cancer: a systematic review of prospective trials. Oral Oncol. 2020;103:104608. [CrossRef] [Medline]
    16. Fakhry C, Andersen KK, Eisele DW, Gillison ML. Oropharyngeal cancer survivorship in Denmark, 1977-2012. Oral Oncol. 2015;51(11):982-984. [CrossRef] [Medline]
    17. Duman-Lubberding S, van Uden-Kraan CF, Peek N, Cuijpers P, Leemans CR, de Leeuw IMV. An eHealth application in head and neck cancer survivorship care: health care professionals' perspectives. J Med Internet Res. 2015;17(10):e235. [FREE Full text] [CrossRef] [Medline]
    18. Mujcic A, Blankers M, Boon B, de Leeuw IMV, Smit F, van Laar M, et al. Effectiveness, cost-effectiveness, and cost-utility of a digital smoking cessation intervention for cancer survivors: health economic evaluation and outcomes of a pragmatic randomized controlled trial. J Med Internet Res. 2022;24(3):e27588. [FREE Full text] [CrossRef] [Medline]
    19. van der Schroeff MP, van de Schans SAM, Piccirillo JF, Langeveld TPM, de Jong RJB, Janssen-Heijnen MLG. Conditional relative survival in head and neck squamous cell carcinoma: permanent excess mortality risk for long-term survivors. Head Neck. 2010;32(12):1613-1618. [CrossRef] [Medline]
    20. Massa ST, Osazuwa-Peters N, Christopher KM, Arnold LD, Schootman M, Walker RJ, et al. Competing causes of death in the head and neck cancer population. Oral Oncol. 2017;65:8-15. [CrossRef] [Medline]
    21. Massa ST, Cass LM, Osazuwa-Peters N, Christopher KM, Walker RJ, Varvares MA. Decreased cancer-independent life expectancy in the head and neck cancer population. Head Neck. 2017;39(9):1845-1853. [CrossRef] [Medline]
    22. Mell LK, Dignam JJ, Salama JK, Cohen EEW, Polite BN, Dandekar V, et al. Predictors of competing mortality in advanced head and neck cancer. J Clin Oncol. 2010;28(1):15-20. [FREE Full text] [CrossRef] [Medline]
    23. Baxi SS, Pinheiro LC, Patil SM, Pfister DG, Oeffinger KC, Elkin EB. Causes of death in long-term survivors of head and neck cancer. Cancer. 2014;120(10):1507-1513. [FREE Full text] [CrossRef] [Medline]
    24. Larsen CG, Jensen DH, Carlander ALF, Kiss K, Andersen L, Olsen CH, et al. Novel nomograms for survival and progression in HPV+ and HPV- oropharyngeal cancer: a population-based study of 1,542 consecutive patients. Oncotarget. 2016;7(44):71761-71772. [FREE Full text] [CrossRef] [Medline]
    25. 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-269. [FREE Full text] [CrossRef] [Medline]
    26. Nørregaard C, Grønhøj C, Jensen D, Friborg J, Andersen E, von Buchwald C. Cause-specific mortality in HPV+ and HPV- oropharyngeal cancer patients: insights from a population-based cohort. Cancer Med. 2018;7(1):87-94. [FREE Full text] [CrossRef] [Medline]
    27. SEER*Stat database: incidence—SEER research plus data (specialized head and neck fields), 18 registries, nov 2020 sub (2000-2018)—linked to county attributes—total U.S., 1969-2019 counties, National Cancer Institute, DCCPS, Surveillance Research Program, released april 2021, based on the november 2020 submission. Surveillance, Epidemiology, and End Results (SEER) Program (www.seer.cancer.gov). URL: https://seer.cancer.gov/data-software/documentation/seerstat/nov2020/ [accessed 2023-08-08]
    28. Fullerton ZH, Butler SS, Mahal BA, Muralidhar V, Schoenfeld JD, Tishler RB, et al. Short-term mortality risks among patients with oropharynx cancer by human papillomavirus status. Cancer. 2020;126(7):1424-1433. [FREE Full text] [CrossRef] [Medline]
    29. Putter H, Fiocco M, Geskus RB. Tutorial in biostatistics: competing risks and multi-state models. Stat Med. 2007;26(11):2389-2430. [CrossRef] [Medline]
    30. Gorphe P, Classe M, Ammari S, Garcia G, Even C, Casiraghi O, et al. Patterns of disease events and causes of death in patients with HPV-positive versus HPV-negative oropharyngeal carcinoma. Radiother Oncol. 2022;168:40-45. [FREE Full text] [CrossRef] [Medline]
    31. Bosshart SL, Morand GB, Broglie MA. Frequency and localization of second primary tumors in patients with oropharyngeal carcinoma-the influence of the human papilloma virus. Cancers (Basel). 2021;13(8):1755. [FREE Full text] [CrossRef] [Medline]
    32. Culié D, Rousseau A, Pretet JL, Guily JLS. HPV status and therapeutic initial strategy impact on survival and oncologic outcomes: 5-year results from the multicentric prospective cohort of oropharyngeal cancers papillophar. Eur Arch Otorhinolaryngol. 2022;279(6):3071-3078. [CrossRef] [Medline]
    33. Väisänen JA, Alho OP, Koivunen PT, Läärä E. Cause-specific mortality in patients with head and neck cancer: long-term follow-up of a population-based cohort from 1986 to 2012 accounting for competing risks. Oral Oncol. 2018;79:20-26. [CrossRef] [Medline]
    34. Läärä E, Korpi JT, Pitkänen H, Alho OP, Kantola S. Competing risks analysis of cause-specific mortality in patients with oral squamous cell carcinoma. Head Neck. 2017;39(1):56-62. [CrossRef] [Medline]
    35. Rose BS, Jeong JH, Nath SK, Lu SM, Mell LK. Population-based study of competing mortality in head and neck cancer. J Clin Oncol. 2011;29(26):3503-3509. [FREE Full text] [CrossRef] [Medline]
    36. Shen W, Sakamoto N, Yang L. Cancer-specific mortality and competing mortality in patients with head and neck squamous cell carcinoma: a competing risk analysis. Ann Surg Oncol. 2015;22(1):264-271. [CrossRef] [Medline]
    37. Kwon M, Roh JL, Song J, Lee SW, Kim SB, Choi SH, et al. Noncancer health events as a leading cause of competing mortality in advanced head and neck cancer. Ann Oncol. 2014;25(6):1208-1214. [FREE Full text] [CrossRef] [Medline]
    38. Lop J, García J, López M, Taberna M, Mena M, Alemany L, et al. Competing mortality in oropharyngeal carcinoma according to human papillomavirus status. Head Neck. 2019;41(5):1328-1334. [CrossRef] [Medline]
    39. Petrelli F, Sarti E, Barni S. Predictive value of human papillomavirus in oropharyngeal carcinoma treated with radiotherapy: an updated systematic review and meta-analysis of 30 trials. Head Neck. 2014;36(5):750-759. [CrossRef] [Medline]
    40. Liu H, Li J, Zhou Y, Hu Q, Zeng Y, Mohammadreza MM. Human papillomavirus as a favorable prognostic factor in a subset of head and neck squamous cell carcinomas: a meta-analysis. J Med Virol. 2017;89(4):710-725. [CrossRef] [Medline]
    41. Hess CB, Rash DL, Daly ME, Farwell DG, Bishop J, Vaughan AT, et al. Competing causes of death and medical comorbidities among patients with human papillomavirus-positive vs human papillomavirus-negative oropharyngeal carcinoma and impact on adherence to radiotherapy. JAMA Otolaryngol Head Neck Surg. 2014;140(4):312-316. [FREE Full text] [CrossRef] [Medline]
    42. Young L, Cho L. Unique cardiovascular risk factors in women. Heart. 2019;105(21):1656-1660. [CrossRef] [Medline]
    43. He J, Zhu Z, Bundy JD, Dorans KS, Chen J, Hamm LL. Trends in cardiovascular risk factors in US adults by race and ethnicity and socioeconomic status, 1999-2018. JAMA. 2021;326(13):1286-1298. [FREE Full text] [CrossRef] [Medline]
    44. Sturgeon KM, Deng L, Bluethmann SM, Zhou S, Trifiletti DM, Jiang C, et al. A population-based study of cardiovascular disease mortality risk in US cancer patients. Eur Heart J. 2019;40(48):3889-3897. [FREE Full text] [CrossRef] [Medline]
    45. Abramowitz L, Guily JLS, Moyal-Barracco M, Bergeron C, Borne H, Dahlab A, et al. Epidemiological and economic burden of potentially HPV-related cancers in France. PLoS One. 2018;13(9):e0202564. [FREE Full text] [CrossRef] [Medline]
    46. Daniel CL, Atkins A. HPV vaccination disparities and strategies to prevent oropharyngeal cancer in males. Prev Med Rep. 2021;24:101521. [FREE Full text] [CrossRef] [Medline]
    47. Gray JL, Maghlouth AA, Hussain HA, Sheef MA. Impact of oral and oropharyngeal cancer diagnosis on smoking cessation patients and cohabiting smokers. Tob Induc Dis. 2019;17:75. [FREE Full text] [CrossRef] [Medline]
    48. Toman J, Von Larson S, Umeno H, Kurita T, Furusaka T, Hasegawa H, et al. HPV-positive oropharyngeal cancer via p16 immunohistochemistry in Japan. Ann Otol Rhinol Laryngol. 2017;126(2):152-158. [CrossRef] [Medline]
    49. Husain S, Lohia S, Petkov V, Blackwell T, Swisher-McClure S, Mizrachi A, et al. Disparities and guideline adherence for HPV testing among patients with oropharyngeal squamous cell carcinoma, NCDB, and SEER. Head Neck. 2021;43(7):2110-2123. [FREE Full text] [CrossRef] [Medline]

    ‎
    AJCC: American Joint Committee on Cancer
    COD: cause of death
    HNC: head and neck cancer
    HNCSM: head and neck cancer–specific mortality
    HNSCC: head and neck squamous cell carcinoma
    HPV: human papillomavirus
    NCCM: noncancer-caused mortality
    OM: overall mortality
    OPSCC: oropharyngeal squamous cell carcinoma
    SEER: Surveillance, Epidemiology, and End Result
    sHR: subdistribution hazard ratio
    SPCM: second primary cancer mortality

    Edited by A Mavragani, T Sanchez; submitted 25.03.23; peer-reviewed by J Wang, M Khokhar; comments to author 19.04.23; revised version received 14.05.23; accepted 25.07.23; published 29.08.23.

    Copyright

    ©Dong-Dong Zhang, Min Lei, Yue Wang, Pei-Ji Zeng, Yong-Jun Hong, Cheng-Fu Cai. Originally published in JMIR Public Health and Surveillance (https://publichealth.jmir.org), 29.08.2023.

    This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Public Health and Surveillance, is properly cited. The complete bibliographic information, a link to the original publication on https://publichealth.jmir.org, as well as this copyright and license information must be included.