Background:

As human papillomavirus positive (HPV+) oral cavity and pharynx cancer (OCPC) incidence increases significantly, our objective was to determine whether selected sociodemographic and clinical factors were associated with HPV+ OCPCs overall and by oropharyngeal and non-oropharyngeal sites.

Methods:

Surveillance, Epidemiology and End Results (SEER) Program data were used in this study. Specifically, univariate and logistic regression models were used to examine the relationships between HPV+ and HPV− OCPC cases and age, sex, race, ethnicity, marital status, factors of neighborhood socioeconomic status (i.e., nSES/Yost index) and rurality/urbanity, first malignancy status, histology, reporting source, stage at diagnosis, and OCPC anatomic site. The same approach was used to identify risk factors for HPV positivity for oropharyngeal and non-oropharyngeal OCPCs separately.

Results:

In all OCPCs, cases that were male, <80 years old, lived in the four highest nSES categories, diagnosed with a non-“gum and other mouth” OCPC (ref = hypopharynx), not locally staged at diagnosis, and a first malignancy had higher odds of being HPV+. Cases that were American Indian/Alaska Native and Asian or Pacific Islander (ref = White), Spanish-Hispanic-Latino ethnicity, non-married/partnered, and not reported by a hospital/clinic had lower odds of being HPV+. Associations were maintained in oropharyngeal OCPCs and only age and race remained significant for non-oropharyngeal OCPCs.

Conclusions:

Sociodemographic and clinical differences in HPV+ and HPV− OCPC, overall and for (non)oropharyngeal, cases exist.

Impact:

Identification of OCPC and (non)oropharyngeal risk factors for HPV positivity may assist in discovering high-risk groups that should receive enhanced public health efforts to reduce the U.S. OCPC burden.

Oral cavity and pharynx cancer (OCPC) incidence and mortality rates have been on upward trends for the past decade (1). It is estimated that there will be 54,000 new OCPC cases identified in the United States and an estimated 11,230 OCPC deaths in 2022 (1). The 5-year relative survival rate for those diagnosed with OCPC is 68.0%, with earlier-stage diagnoses having more favorable survival rates (SEER*Explorer, 2022, SEER RRID: SCR_003293). OCPC disparities are apparent with males having an incidence more than double that of females (17.4 per 100,000 males vs. 6.4 per 100,000 females) and Blacks having a lower 5-year survival rate compared with Whites (52.4% vs. 69.2%, respectively; SEER*Explorer, 2022, SEER RRID: SCR_003293).

The human papillomavirus (HPV) has consistently been identified as a common OCPC risk factor (2–4). OCPC incidence, mortality, and survival rates vary by HPV status. HPV-positive (HPV+) OCPC diagnoses typically occur at later stages and need more drastic treatments as currently there is no standard screening protocol nor pharmacotherapeutic treatment for the earlier presence of an oral HPV infection (5). The number of HPV+ OCPC cases has increased by almost 50% since 1973 (6–9), occurring most often in males who are typically not heavy tobacco users (10). HPV+ OCPC mortality rates have been on a similar trajectory as HPV+ OCPC incidence rates, experiencing equivalent rate increases over time (8). The survival rate from an HPV+ OCPC diagnosis can be as high as 71% to 93% for those in lower risk groups (i.e., less tumor extension, no lymphnode involvement, non-smoking; ref. 5). Surprisingly, higher socioeconomic status (SES) has been associated with higher risk of HPV+ OCPCs (10). This unique disparity may offer some insight into the improved survival rate for HPV+ OCPCs, compared with those diagnosed with HPV negative (HPV−) OCPCs (10).

Only cancers from very well-defined specific OCPC anatomies have been routinely HPV tested and investigated (11–24). Specifically, HPV has been associated with more than 70% of all U.S. cancers found in the oropharyngeal region of the head and neck (25), predominately the tonsils and base of the tongue (3, 5, 6, 24, 26). However, the oral cavity and pharynx anatomic region is large with OCPCs occurring from the nasopharynx through the hypopharynx, with incidence, mortality, and survival rates also differing across anatomic sites (27). OCPC causal pathways are complex and multifactorial, with causes (e.g., HPV, tobacco/alcohol) potentially interacting with one another (i.e., independently, complementarily, or synergistically). Therefore, OCPC sites thought to be related to other risk factors (e.g., alcohol/tobacco) could also be HPV+ (28, 29). Additional OCPC anatomies also remain un(der)investigated and could be associated with HPV positivity (5, 30). For instance, a national study of cancer registry–representative tumors was HPV tested, and HPV was identified in 32.0% of oral cavity tumors—an OCPC site thought to not be HPV-related (24). Non-oropharyngeal OCPC cases may have different risk factors for HPV positivity (11, 18, 23) than those for oropharyngeal OCPCs (13–16, 19, 23), but research is limited as prior studies have concentrated on better understanding the latter.

Much of the existing national-level epidemiologic findings on OCPC cases with known HPV status are outdated and/or restricted in generalizability, relying on specific OCPC anatomic sites or populations (11–23). Moreover, sociodemographic factors reflecting neighborhood SES (nSES) and rurality/urbanity, potentially related to OCPC diagnoses, are untapped risk factors that may hold some of the keys to understanding differences in HPV+ and HPV− OCPCs. Given the significant increases in incidence rates of HPV+ OCPCs (3, 6–9, 31), this lack of relevant epidemiologic information is problematic as these risk factors could elucidate important causal pathways and identify high-risk populations which may aid disease control. Therefore, the objective of this investigation was to evaluate associations between known HPV status (i.e., from HPV testing) and demographic, census tract-level, and clinical factors among OCPC cases, overall, and for oropharyngeal and non-oropharyngeal sites independently, using data from a large, well-described collection of central cancer registries.

The Surveillance, Epidemiology and End Results (SEER) Program is a collection of population-based central cancer registries capturing cases from 22 registries across the United States. With permission from the SEER Program, the SEER Head and Neck with HPV Status and Census Tract-level SES/Rurality Combined Database was used to explore associations between known HPV status and demographic, census-tract level (nSES and rurality/urbanity), and clinical characteristics (SEER*Stat, December 2021, SEER RRID: SCR_003293). Institutional Review Board approval was not required, and ethical consent was not warranted because SEER data are publicly available and deidentified.

Invasive OCPC cases diagnosed and reported to SEER Program registries between 2014 and 2017 were included in this investigation (N = 25,357). Participants were excluded if they had missing socioeconomic data (i.e., nSES level, n = 1,535) or had a reporting source of death certificate only or autopsy only (n = 199). In addition, participants under 30 years old were excluded from analyses (n = 223) as their OCPCs are typically of a unique oral HPV-related disease presentation (i.e., juvenile onset recurrent respiratory papillomatosis) with a nonmodifiable transmission pathway (i.e., vertical; ref. 32).

SEER-defined sociodemographic variables were investigated as potential risk factors for HPV positivity in OCPC cases, overall, and more specifically in oropharyngeal and non-oropharyngeal OCPC sites. With HPV test status as the outcome of interest, OCPC cases were classified into three status groups based on their recorded HPV status—(i) HPV+; (ii) HPV−; and (iii) unknown. (As the SEER Program includes a collection of central cancer registries, HPV positivity was determined independently by reporting sources with varying HPV testing methods.) Sex was “Male” or “Female” and age at diagnosis was described using 10-year incremental categories (i.e., 30–39, 40–49, 50–59, 60–69, 70–79, 80+ years old). Race was described using the variable “Race Recode (W, B, AI, API)”, which classified race as “White”, “Black”, “Asian or Pacific Islander”, “American Indian or Alaskan Native”, or “unknown”. Ethnicity was dichotomized as “Spanish-Hispanic-Latino” or “non-Spanish-Hispanic-Latino”, using the variable “Origin Recode NHIA”. Marital status was defined as “divorced/separated/widowed”, “married/domestic partnered”, “single” or “unknown”. SEER's nSES, or Yost index, is a time-dependent composite score, including census tract-level information about education index, household income, percent below 150% of poverty line, median house value, percent unemployed, median rent, and percent working class variables (33–35). nSES scores were classified into quintiles, with the first and fifth quintiles describing the lowest and highest nSES groups, respectively. SEER census tract-level rurality data come from the U.S. Department of Agriculture's Rural Urban Commuting Area (RUCA) codes (36), using four levels (“All Urban”, “Mostly Urban”, “Mostly Rural”, “All Rural”).

SEER's clinical variables provided additional details on tumor site, histology, stage at diagnosis, first malignancy, diagnostic confirmation, and reporting source. The SEER variable “Site Recode”, based on International Classification of Disease for Oncology version-3 (ICD-O-3) topography codes (37), was used to determine specific anatomies within the category of OCPC (i.e., OCPC site), including “tongue” (ICD-O-3: C019), “tonsil” (ICD-O-3: C024, C090-099), “oropharynx” (ICD-O-3: C100-109), “nasopharynx” (ICD-O-3: C110-119), “pyriform sinus” (ICD-O-3: C129), “hypopharynx” (ICD-O-3: C130-139), “gum and other mouth” (ICD-O-3: C051, 052), and “other oral cavity and pharynx” (ICD-O-3: C140-148). To be comprehensive in the investigation of all potential OCPC risk factors for HPV positivity, all OCPC sites included in the dataset, as defined by the SEER “Site Recode” variable, were included in the current investigation. [It should be noted that SEER has restricted the dataset to only OCPC sites believed to be potentially associated with HPV instead of including all OCPC (sub)sites (e.g., “tongue” only included base of tongue (ICD-O-3: C019), not non-base subsites (ICD-O-3: C020-023, C025-029); therefore, the category is referred to as “base of tongue” hereafter]. For additional analyses, OCPC sites were later recategorized as oropharyngeal (i.e., base of tongue (ICD-O-3: C019), tonsil (ICD-O-3: C024, C090-099), oropharynx (ICD-O-3: C100-109), gum and other mouth (ICD-O-3: C051, C052), other oral cavity and pharynx (ICD-O-3: C140, C142, C148)) and non-oropharyngeal [i.e., nasopharynx (ICD-O-3: C110-119), pyriform sinus (ICD-O-3: C129), hypopharynx (ICD-O-3: C130-139)] (38). Histology was dichotomized as squamous cell neoplasms (ICD-O-3: 8050-8076, 8078, 8083, 8084, 8094) and non-squamous cell neoplasms (ICD-O-3: 8000-8049, 8090-8389, 8430-8579, 8800-8839, 8850-8999, 9120-9169, 9180-9249, 9350-9379, 9490-9529). Stage at diagnosis was characterized using “SEER Combined Summary Stage” and was classified as “Localized,” “Regional,” “Distant,” and “Unknown/unstaged.” As indicated by the name, SEER's “first malignant primary indicator” variable collected data on whether the current OCPC diagnosis was the individual's first malignant tumor or not. Diagnostic confirmation described how diagnoses of OCPC were verified, ranging from clinical diagnosis to positive cytology/histology/radiography [with(out) microscopic confirmation] to unknown methodology. The origin of OCPC diagnosis was described by the SEER variable “type of reporting source” which was dichotomized as either “hospital or clinic” or “other.”

Descriptive statistics [i.e., frequencies (percentages)] for sociodemographic and clinical characteristics were estimated for OCPC cases by HPV status (+, −, or unknown). Differences in missing data (i.e., HPV status) were determined by a 10% difference in the distribution of each characteristic between those with known and unknown HPV status. Unknown HPV status and histologically non-squamous cell neoplasm OCPCs were then removed such that the analysis set only included those with known (i.e., positive or negative) HPV status and squamous cell neoplasms to best assess factors associated with OCPC HPV positivity. Univariate logistic regression models were used to determine whether characteristics were associated with HPV+ status. All variables that were statistically significant (P < 0.05) univariately were included in a multivariable logistic regression model to identify characteristics and populations that continued to be associated with HPV positivity. This same approach was used to identify risk factors for HPV positivity for oropharyngeal and non-oropharyngeal OCPCs separately. ORs (univariate and adjusted) and 95% confidence intervals (CI) were calculated for all models. Variables in which almost all cases were concentrated within one category even after regrouping (e.g., diagnostic confirmation), there was no meaningful contribution to determining HPV status (e.g., year of diagnosis), or there was nonsignificance or collinearity with other variables of interest (e.g., RUCA) were omitted. Analyses were conducted using SAS 9.4 (SAS Institute Inc., SAS RRID: SCR_008567) and Stata 16 (StataCorp, Stata RRID: SCR_012763).

Data availability

The data analyzed in this study were obtained from the SEER Program at www.seer.cancer.gov.

Of the 23,400 OCPCs’ meeting inclusion criteria, 42% (n = 9,832) were HPV+, 18% (n = 4,211) were HPV−, and 40% (n = 9,357) were of unknown HPV status. (Because SEER has restricted the dataset to only OCPC (sub)sites believed to be potentially associated with HPV instead of including all OCPC sites, prevalence of HPV positivity among all U.S. OCPC cases cannot be concluded.) Most tumors were diagnostically confirmed with positive histology (93%). Table 1 shows the breakdown of demographic, census tract-level, and clinical factors according to known HPV status. Overall, most OCPC cases were male, 50 to 69 years old, White, non-Spanish-Hispanic-Latino, and married/domestic partnered. They also tended to reside in urban census tracts and census tracts in the highest (i.e., fifth) nSES quintile. Clinically, most OCPCs were diagnosed on the base of tongue or tonsil, as squamous cell neoplasms, at the regional stage, by a hospital or clinic, and as first malignancies (Table 1). Cases of unknown HPV status had similar sociodemographic and clinical attributes as cases with known HPV status, except for stage at diagnosis which differed by more than 10% (Supplementary Table S1).

Table 1.

Frequencies (percentages) of demographic, socioeconomic, and clinical OCPC data by HPV statusa.

VariableLevelHPV- (n = 4,211) (Col %)HPV+ (n = 9,832) (Col %)Total (n = 14,043) (Col %)
Demographics 
Sex 
 Female 1,051 (25%) 1,358 (14%)  2,409 (17%) 
 Male 3,160 (75%) 8,474 (86%)  11,634 (83%) 
Age 
 30–39 years 62 (1%) 90 (1%)  152 (1%) 
 40–49 years 313 (7%) 944 (10%)  1,257 (9%) 
 50–59 years 1,188 (28%) 3,292 (33%)  4,480 (32%) 
 60–69 years 1,466 (35%) 3,572 (36%)  5,038 (36%) 
 70–79 years 829 (20%) 1,570 (16%)  2,399 (17%) 
 80+ years 353 (8%) 364 (4%)  717 (5%) 
Race 
 White 3,217 (76%) 8,800 (90%)  12,017 (86%) 
 American Indian/Alaska Native 25 (1%) 56 (1%)  81 (1%) 
 Asian or Pacific Islander 375 (9%) 305 (3%)  680 (5%) 
 Black 572 (14%) 618 (6%)  1,190 (9%) 
 Unknown 22 (1%) 53 (1%)  75 (1%) 
Ethnicity 
 Non-Spanish-Hispanic-Latino 3,894 (92%) 9,202 (94%)  13,096 (93%) 
 Spanish-Hispanic-Latino 317 (8%) 630 (6%)  947 (7%) 
Marital status 
 Married/Domestic Partner 2,055 (49%) 6,043 (61%)  8,098 (58%) 
 Divorced/Separate/Widowed 964 (23%) 1,690 (17%)  2,654 (19%) 
 Single 948 (23%) 1,659 (17%)  2,607 (19%) 
 Unknown 244 (6%) 440 (4%)  684 (5%) 
Census tract-based factors 
nSES (Yost) 
 Group 1 (lowest) 915 (22%) 1,183 (12%)  2,098 (15%) 
 Group 2 705 (17%) 1,442 (15%)  2,147 (15%) 
 Group 3 720 (17%) 1,710 (17%)  2,430 (17%) 
 Group 4 836 (20%) 2,352 (24%)  3,188 (23%) 
 Group 5 (highest) 1,035 (25%) 3,145 (32%)  4,180 (30%) 
RUCA 
 All rural 270 (6%) 636 (6%)  906 (6%) 
 All urban 2,819 (67%) 6,259 (64%)  9,078 (65%) 
 Mostly rural 282 (7%) 747 (8%)  1,029 (7%) 
 Mostly urban 840 (20%) 2,190 (22%)  3,030 (22%) 
Clinical data 
Diagnosis year 
 2014 1,009 (24%) 2,027 (21%)  3,036 (22%) 
 2015 1,015 (24%) 2,266 (23%)  3,281 (23%) 
 2016 1,142 (27%) 2,708 (28%)  3,850 (27%) 
 2017 1,045 (25%) 2,831 (29%)  3,876 (28%) 
OCPC site 
 Gum and other mouth 263 (6%) 80 (1%)  343 (2%) 
 Hypopharynx 643 (15%) 233 (2%)  876 (6%) 
 Nasopharynx 474 (11%) 270 (3%)  744 (5%) 
 Oropharynx 439 (10%) 604 (6%)  1,043 (7%) 
 Other oral cavity and pharynx 171 (4%) 193 (2%)  364 (3%) 
 Base of tongue 1,120 (27%) 3,708 (38%)  4,828 (34%) 
 Tonsil 1,101 (26%) 4,744 (48%)  5,845 (42%) 
Histology 
 Squamous cell neoplasms 4,009 (95%) 9,726 (99%)  13,735 (98%) 
 Non-squamous cell neoplasms 202 (5%) 106 (1%)  308 (2%) 
Stage at diagnosis 
 Localized 723 (17%) 787 (8%)  1,510 (11%) 
 Regional 2,425 (58%) 7,405 (75%)  9,830 (70%) 
 Distant 978 (23%) 1,497 (15%)  2,475 (18%) 
 Unknown/unstaged 85 (2%) 143 (1%)  228 (2%) 
Report source 
 Hospital or Clinic 3,997 (95%) 9,458 (96%)  13,455 (96%) 
 Other 214 (5%) 374 (4%)  588 (4%) 
First malignancy 
 No 1,035 (25%) 1,339 (14%)  2,374 (17%) 
 Yes 3,176 (75%) 8,493 (86%)  11,669 (83%) 
VariableLevelHPV- (n = 4,211) (Col %)HPV+ (n = 9,832) (Col %)Total (n = 14,043) (Col %)
Demographics 
Sex 
 Female 1,051 (25%) 1,358 (14%)  2,409 (17%) 
 Male 3,160 (75%) 8,474 (86%)  11,634 (83%) 
Age 
 30–39 years 62 (1%) 90 (1%)  152 (1%) 
 40–49 years 313 (7%) 944 (10%)  1,257 (9%) 
 50–59 years 1,188 (28%) 3,292 (33%)  4,480 (32%) 
 60–69 years 1,466 (35%) 3,572 (36%)  5,038 (36%) 
 70–79 years 829 (20%) 1,570 (16%)  2,399 (17%) 
 80+ years 353 (8%) 364 (4%)  717 (5%) 
Race 
 White 3,217 (76%) 8,800 (90%)  12,017 (86%) 
 American Indian/Alaska Native 25 (1%) 56 (1%)  81 (1%) 
 Asian or Pacific Islander 375 (9%) 305 (3%)  680 (5%) 
 Black 572 (14%) 618 (6%)  1,190 (9%) 
 Unknown 22 (1%) 53 (1%)  75 (1%) 
Ethnicity 
 Non-Spanish-Hispanic-Latino 3,894 (92%) 9,202 (94%)  13,096 (93%) 
 Spanish-Hispanic-Latino 317 (8%) 630 (6%)  947 (7%) 
Marital status 
 Married/Domestic Partner 2,055 (49%) 6,043 (61%)  8,098 (58%) 
 Divorced/Separate/Widowed 964 (23%) 1,690 (17%)  2,654 (19%) 
 Single 948 (23%) 1,659 (17%)  2,607 (19%) 
 Unknown 244 (6%) 440 (4%)  684 (5%) 
Census tract-based factors 
nSES (Yost) 
 Group 1 (lowest) 915 (22%) 1,183 (12%)  2,098 (15%) 
 Group 2 705 (17%) 1,442 (15%)  2,147 (15%) 
 Group 3 720 (17%) 1,710 (17%)  2,430 (17%) 
 Group 4 836 (20%) 2,352 (24%)  3,188 (23%) 
 Group 5 (highest) 1,035 (25%) 3,145 (32%)  4,180 (30%) 
RUCA 
 All rural 270 (6%) 636 (6%)  906 (6%) 
 All urban 2,819 (67%) 6,259 (64%)  9,078 (65%) 
 Mostly rural 282 (7%) 747 (8%)  1,029 (7%) 
 Mostly urban 840 (20%) 2,190 (22%)  3,030 (22%) 
Clinical data 
Diagnosis year 
 2014 1,009 (24%) 2,027 (21%)  3,036 (22%) 
 2015 1,015 (24%) 2,266 (23%)  3,281 (23%) 
 2016 1,142 (27%) 2,708 (28%)  3,850 (27%) 
 2017 1,045 (25%) 2,831 (29%)  3,876 (28%) 
OCPC site 
 Gum and other mouth 263 (6%) 80 (1%)  343 (2%) 
 Hypopharynx 643 (15%) 233 (2%)  876 (6%) 
 Nasopharynx 474 (11%) 270 (3%)  744 (5%) 
 Oropharynx 439 (10%) 604 (6%)  1,043 (7%) 
 Other oral cavity and pharynx 171 (4%) 193 (2%)  364 (3%) 
 Base of tongue 1,120 (27%) 3,708 (38%)  4,828 (34%) 
 Tonsil 1,101 (26%) 4,744 (48%)  5,845 (42%) 
Histology 
 Squamous cell neoplasms 4,009 (95%) 9,726 (99%)  13,735 (98%) 
 Non-squamous cell neoplasms 202 (5%) 106 (1%)  308 (2%) 
Stage at diagnosis 
 Localized 723 (17%) 787 (8%)  1,510 (11%) 
 Regional 2,425 (58%) 7,405 (75%)  9,830 (70%) 
 Distant 978 (23%) 1,497 (15%)  2,475 (18%) 
 Unknown/unstaged 85 (2%) 143 (1%)  228 (2%) 
Report source 
 Hospital or Clinic 3,997 (95%) 9,458 (96%)  13,455 (96%) 
 Other 214 (5%) 374 (4%)  588 (4%) 
First malignancy 
 No 1,035 (25%) 1,339 (14%)  2,374 (17%) 
 Yes 3,176 (75%) 8,493 (86%)  11,669 (83%) 

aSEER*Stat Database: Incidence - SEER Research Plus Data (Specialized Head and Neck Fields with Census Tract SES/Rurality), 18 Registries (excl AK), Nov 2020 Sub (2006–2018) <Vintage 2019 Pops by Tract 2010 Geographies> - Linked To Census Tract Attributes SEER 18 excl AK, Census 2010 Geographies Time Dependent (2006–2018; 2008 copied to 06, 07; 2017 copied to 18). National Cancer Institute, DCCPS, Surveillance Research Program; December 2021. http://www.seer.cancer.gov/. Accessed January 24, 2022. (SEER*Stat, December 2021, SEER RRID: SCR_003293).

When OCPC cases with unknown HPV status (n = 9,357) and histology of non-squamous cell neoplasms (n = 308) were removed to more accurately assess factors associated with OCPC HPV positivity, 13,735 cases with a known HPV status remained [n = 4,009 (HPV−), n = 9,726 (HPV+)]. Because of the high prevalence of HPV+ OCPCs, ORs should not be interpreted as prevalence ratios as they can overestimate the strength of associations. Univariately, the odds of OCPCs being HPV+ were significantly higher for males, those of the five youngest age groups (ref = 80+ years old), those residing in census tracts in each of the four highest nSES quintiles [ref = lowest quintile (i.e., group 1)], those diagnosed with an OCPC in the sites of nasopharynx, oropharynx, other oral cavity and pharynx, base of tongue, and tonsil (ref = hypopharynx), those diagnosed at the regional, distant, and unknown/unstaged stage (ref = local stage), and those for whom this was the first malignancy (Table 2). Most notably, the odds of HPV positivity in OCPCs were two to three times higher for males, those ages 40 to 69 years (vs. those ages 80+ years), those diagnosed at the regional stage (vs. local), and those with first malignancies. The odds of being HPV+ was more than nine times higher for OCPCs on the base of tongue or tonsils than those of the hypopharynx (Table 2). Univariately, odds of OCPCs being HPV+ were significantly lower for American Indians/Alaska Natives and Asians or Pacific Islanders (ref = White), those of Spanish-Hispanic-Latino ethnicity, those with any of the three non-married/domestic partnered marital statuses, and those with an “other” reporting source (ref = hospital or clinic). Interestingly, odds of OCPCs being HPV+ were at least 60% lower among American Indians/Alaska Natives and Asians or Pacific Islanders as compared with Whites (Table 2). On the other hand, neither residence in a rural or urban census tract was significantly associated with OCPCs being HPV+ (Table 2).

Table 2.

Odds of HPV+ oral cavity and pharynx squamous cell neoplasm cancer by demographic, socioeconomic, and clinical dataa.

VariableLevelCrude HPV prevalenceUnadjusted ORb (95% CI)cPdAdjusted ORb (95% CI)cPd
Demographics 
Sex 
 Female (ref) 58% 1.0 — 1.0 — 
 Male 73% 2.04 (1.86–2.24) <0.001 1.65 (1.49–1.84) <0.001 
Age 
 30–39 years 66% 1.79 (1.21–2.64) 0.003 2.33 (1.50–3.62) <0.001 
 40–49 years 77% 3.01 (2.47–3.68) <0.001 2.59 (2.06–3.26) <0.001 
 50–59 years 74% 2.66 (2.26–3.13) <0.001 2.10 (1.74–2.53) <0.001 
 60–69 years 71% 2.32 (1.97–2.72) <0.001 1.92 (1.60–2.30) <0.001 
 70–79 years 66% 1.80 (1.52–2.14) <0.001 1.70 (1.40–2.06) <0.001 
 80+ years (ref) 52% 1.0 — 1.0 — 
Race 
 White (ref) 74% 1.0 — 1.0 — 
 American Indian/Alaska Native 68% 0.39 (0.34–0.44) <0.001 0.55 (0.48–0.64) <0.001 
 Asian or Pacific Islander 48% 0.33 (0.28–0.39) <0.001 0.45 (0.37–0.54) <0.001 
 Black 52% 0.77 (0.48–1.23) 0.275 0.80 (0.47–1.37) 0.423 
 Unknown 72% 0.92 (0.55–1.55) 0.763 1.08 (0.61–1.92) 0.786 
Ethnicity 
 Non-Spanish-Hispanic-Latino (ref) 71% 1.0 — 1.0 — 
 Spanish-Hispanic-Latino 67% 0.84 (0.73–0.97) 0.020 0.75 (0.64–0.88) 0.001 
Marital status 
 Married/Domestic Partnered (ref) 76% 1.0 — 1.0 — 
 Divorced/Separated/Widowed 64% 0.59 (0.53–0.65) <0.001 0.78 (0.70–0.87) <0.001 
 Single 64% 0.58 (0.53–0.64) <0.001 0.68 (0.61–0.76) <0.001 
 Unknown 65% 0.61 (0.51–0.72) <0.001 0.72 (0.59–0.87) 0.001 
Census tract-based factors 
nSES (Yost) 
 Group 1 (lowest) (ref) 57% 1.0 — 1.0 — 
 Group 2 68% 1.60 (1.41–1.82) <0.001 1.38 (1.19–1.59) <0.001 
 Group 3 71% 1.86 (1.64–2.11) <0.001 1.55 (1.34–1.79) <0.001 
 Group 4 75% 2.24 (1.98–2.52) <0.001 1.80 (1.56–2.06) <0.001 
 Group 5 (highest) 76% 2.40 (2.14–2.68) <0.001 1.82 (1.59–2.08) <0.001 
RUCA 
 All rural (ref) 71% 1.0 — — — 
 All urban 70% 0.95 (0.82–1.11) 0.517 — — 
 Mostly rural 73% 1.11 (0.91–1.36) 0.300 — — 
 Mostly urban 73% 1.11 (0.94–1.32) 0.199 — — 
Clinical data 
OCPC site 
 Hypopharynx (ref) 27% 1.0 — 1.0 — 
 Gum and other mouth 23% 0.82 (0.61–1.10) 0.187 0.98 (0.72–1.34) 0.903 
 Nasopharynx 40% 1.77 (1.42–2.22) <0.001 1.77 (1.40–2.25) <0.001 
 Oropharynx 58% 3.74 (3.08–4.55) <0.001 3.47 (2.83–4.25) <0.001 
 Other oral cavity and pharynx 54% 3.18 (2.46–4.12) <0.001 2.61 (1.99–3.42) <0.001 
 Base of tongue 77% 9.15 (7.76–10.79) <0.001 7.16 (6.03–8.50) <0.001 
 Tonsil 81% 11.80 (10.01–13.91) <0.001 9.60 (8.08–11.39) <0.001 
Stage at diagnosis 
 Localized (ref) 53% 1.0 — 1.0 — 
 Regional 76% 2.77 (2.48–3.10) <0.001 2.05 (1.80–2.33) <0.001 
 Distant 61% 1.40 (1.23–1.60) <0.001 1.43 (1.23–1.67) <0.001 
 Unknown/unstaged 64% 1.58 (1.17–2.12) 0.003 1.96 (1.39–2.75) <0.001 
Report source 
 Hospital or Clinic (ref) 71% 1.0 — 1.0 — 
 Other 64% 0.73 (0.61–0.86) <0.001 0.73 (0.60–0.88) 0.002 
First malignancy 
 No (ref) 57% 1.0 — 1.0 — 
 Yes 74% 2.11 (1.92–2.31) <0.001 1.64 (1.47–1.82) <0.001 
VariableLevelCrude HPV prevalenceUnadjusted ORb (95% CI)cPdAdjusted ORb (95% CI)cPd
Demographics 
Sex 
 Female (ref) 58% 1.0 — 1.0 — 
 Male 73% 2.04 (1.86–2.24) <0.001 1.65 (1.49–1.84) <0.001 
Age 
 30–39 years 66% 1.79 (1.21–2.64) 0.003 2.33 (1.50–3.62) <0.001 
 40–49 years 77% 3.01 (2.47–3.68) <0.001 2.59 (2.06–3.26) <0.001 
 50–59 years 74% 2.66 (2.26–3.13) <0.001 2.10 (1.74–2.53) <0.001 
 60–69 years 71% 2.32 (1.97–2.72) <0.001 1.92 (1.60–2.30) <0.001 
 70–79 years 66% 1.80 (1.52–2.14) <0.001 1.70 (1.40–2.06) <0.001 
 80+ years (ref) 52% 1.0 — 1.0 — 
Race 
 White (ref) 74% 1.0 — 1.0 — 
 American Indian/Alaska Native 68% 0.39 (0.34–0.44) <0.001 0.55 (0.48–0.64) <0.001 
 Asian or Pacific Islander 48% 0.33 (0.28–0.39) <0.001 0.45 (0.37–0.54) <0.001 
 Black 52% 0.77 (0.48–1.23) 0.275 0.80 (0.47–1.37) 0.423 
 Unknown 72% 0.92 (0.55–1.55) 0.763 1.08 (0.61–1.92) 0.786 
Ethnicity 
 Non-Spanish-Hispanic-Latino (ref) 71% 1.0 — 1.0 — 
 Spanish-Hispanic-Latino 67% 0.84 (0.73–0.97) 0.020 0.75 (0.64–0.88) 0.001 
Marital status 
 Married/Domestic Partnered (ref) 76% 1.0 — 1.0 — 
 Divorced/Separated/Widowed 64% 0.59 (0.53–0.65) <0.001 0.78 (0.70–0.87) <0.001 
 Single 64% 0.58 (0.53–0.64) <0.001 0.68 (0.61–0.76) <0.001 
 Unknown 65% 0.61 (0.51–0.72) <0.001 0.72 (0.59–0.87) 0.001 
Census tract-based factors 
nSES (Yost) 
 Group 1 (lowest) (ref) 57% 1.0 — 1.0 — 
 Group 2 68% 1.60 (1.41–1.82) <0.001 1.38 (1.19–1.59) <0.001 
 Group 3 71% 1.86 (1.64–2.11) <0.001 1.55 (1.34–1.79) <0.001 
 Group 4 75% 2.24 (1.98–2.52) <0.001 1.80 (1.56–2.06) <0.001 
 Group 5 (highest) 76% 2.40 (2.14–2.68) <0.001 1.82 (1.59–2.08) <0.001 
RUCA 
 All rural (ref) 71% 1.0 — — — 
 All urban 70% 0.95 (0.82–1.11) 0.517 — — 
 Mostly rural 73% 1.11 (0.91–1.36) 0.300 — — 
 Mostly urban 73% 1.11 (0.94–1.32) 0.199 — — 
Clinical data 
OCPC site 
 Hypopharynx (ref) 27% 1.0 — 1.0 — 
 Gum and other mouth 23% 0.82 (0.61–1.10) 0.187 0.98 (0.72–1.34) 0.903 
 Nasopharynx 40% 1.77 (1.42–2.22) <0.001 1.77 (1.40–2.25) <0.001 
 Oropharynx 58% 3.74 (3.08–4.55) <0.001 3.47 (2.83–4.25) <0.001 
 Other oral cavity and pharynx 54% 3.18 (2.46–4.12) <0.001 2.61 (1.99–3.42) <0.001 
 Base of tongue 77% 9.15 (7.76–10.79) <0.001 7.16 (6.03–8.50) <0.001 
 Tonsil 81% 11.80 (10.01–13.91) <0.001 9.60 (8.08–11.39) <0.001 
Stage at diagnosis 
 Localized (ref) 53% 1.0 — 1.0 — 
 Regional 76% 2.77 (2.48–3.10) <0.001 2.05 (1.80–2.33) <0.001 
 Distant 61% 1.40 (1.23–1.60) <0.001 1.43 (1.23–1.67) <0.001 
 Unknown/unstaged 64% 1.58 (1.17–2.12) 0.003 1.96 (1.39–2.75) <0.001 
Report source 
 Hospital or Clinic (ref) 71% 1.0 — 1.0 — 
 Other 64% 0.73 (0.61–0.86) <0.001 0.73 (0.60–0.88) 0.002 
First malignancy 
 No (ref) 57% 1.0 — 1.0 — 
 Yes 74% 2.11 (1.92–2.31) <0.001 1.64 (1.47–1.82) <0.001 

Note: bold text = statistical significance.

aSEER*Stat Database: Incidence - SEER Research Plus Data (Specialized Head and Neck Fields with Census Tract SES/Rurality), 18 Registries (excl AK), Nov 2020 Sub (2006–2018) <Vintage 2019 Pops by Tract 2010 Geographies> - Linked To Census Tract Attributes SEER 18 excl AK, Census 2010 Geographies Time Dependent (2006–2018; 2008 copied to 06, 07; 2017 copied to 18). National Cancer Institute, DCCPS, Surveillance Research Program; December 2021. http://www.seer.cancer.gov/. Accessed January 24, 2022. (SEER*Stat, Dec 2021, SEER RRID: SCR_003293).

bOdds ratio.

c95% confidence interval.

dα ≤ 0.05.

Generally speaking, the univariate associations were maintained in the multivariable model for each sociodemographic and clinical characteristic, when holding all other variables constant. In other words, males, those of the five youngest age groups (i.e., 30–79 year olds), those residing in census tracts in each of the four highest nSES quintiles (i.e., groups 2–5), those diagnosed with an OCPC in the sites of nasopharynx, oropharynx, other oral cavity and pharynx, base of tongue, and tonsil (ref = hypopharynx), those diagnosed at regional, distant, and unknown/unstaged stage at diagnosis (ref = local stage), and those for whom this was the first malignancy retained higher odds of their OCPCs being HPV+. Younger cases (i.e., 30–59 year olds) had at least two times the odds of their OCPCs being HPV+ while older cases (i.e., 60–79 year olds) only had an odds for HPV positivity that was 70% to 92% higher (ref = 80+ year olds). OCPCs of the base of tongue and tonsils were still more than seven and nine times, respectively, likely to be HPV+ than those of the hypopharynx. Compared with those residing in census tracts in the lowest nSES quintile (i.e., group 1), the odds of OCPCs being HPV+ continued to increase through nSES quintiles, suggesting a dose–response relationship (Table 2). On the other hand, American Indians/Alaska Natives and Asians or Pacific Islanders (ref = White), those of Spanish-Hispanic-Latino ethnicity, those with any of the three non-married/domestic partnered marital statuses, and those with an “other” reporting source (ref = hospital or clinic) retained lower odds for their OCPCs being HPV+ (Table 2). More specifically, OCPCs among single and divorced/separated/widowed individuals had a significantly lower odds (32% and 22%, respectively) to be HPV+ than those of their married/partnered counterparts (Table 2).

There were 12,281 oropharyngeal cases [n = 3,021 (HPV−), n = 9,260 (HPV+)] and 1,454 non-oropharyngeal cases [n = 988 (HPV−), n = 466 (HPV+)]. All univariate and multivariable associations for HPV+ OCPCs in general, including directionality and strength of associations, were maintained for oropharyngeal cases (Table 3). For non-oropharyngeal OCPCs, only age (i.e., 40–69 year olds), stage at diagnosis (i.e., regional stage only), and first malignancy were univariately associated with higher odds of HPV positivity; race [i.e., American Indian/Alaska Native (borderline significant), Asian or Pacific Islander] was the only factor with a lower odds of HPV positivity (Table 4). In the multivariable model, age and race were the only significant risk factors. More specifically, non-oropharyngeal OCPC cases ages 40 to 69 years had approximately two times higher odds of HPV positivity than 80+ year olds with non-oropharyngeal OCPCs, after adjusting for all other variables in the model. American Indian/Alaska Native Asian or Pacific Islander cases had 34% and 47% lower odds of HPV positivity compared with Whites with non-oropharyngeal OCPCs, respectively (Table 4).

Table 3.

Odds of HPV-positive oropharyngeal squamous cell neoplasm cancer by demographic, socioeconomic, and clinical dataa.

VariableLevelCrude HPV prevalenceUnadjusted ORb (95% CI)cPdAdjusted ORb (95% CI)cPd
Demographics 
Sex 
 Female (ref) 62% 1.0 — 1.0 — 
 Male 78% 2.15 (1.94–2.38) <0.001 1.81 (1.62–2.02) <0.001 
Age 
 30–39 years 79% 2.82 (1.71–4.64) <0.001 3.38 (2.01–5.70) <0.001 
 40–49 years 82% 3.44 (2.75–4.30) <0.001 2.96 (2.33–3.76) <0.001 
 50–59 years 78% 2.74 (2.30–3.28) <0.001 2.29 (1.89–2.77) <0.001 
 60–69 years 75% 2.35 (1.97–2.80) <0.001 1.98 (1.64–2.38) <0.001 
 70–79 years 72% 1.98 (1.64–2.38) <0.001 1.84 (1.51–2.24) <0.001 
 80+ years (ref) 57% 1.0 — 1.0 — 
Race 
 White (ref) 77% 1.0 — 1.0 — 
 American Indian/Alaska Native 75% 0.39 (0.34–0.44) <0.001 0.52 (0.44–0.60) <0.001 
 Asian or Pacific Islander 66% 0.58 (0.46–0.72) <0.001 0.53 (0.42-0.67) <0.001 
 Black 57% 0.86 (0.50–1.50) 0.596 0.79 (0.44–1.40) 0.415 
 Unknown 78% 1.05 (0.58–1.90) 0.882 1.27 (0.67–2.39) 0.467 
Ethnicity 
 Non-Spanish-Hispanic-Latino (ref) 76% 1.0 — 1.0 — 
 Spanish-Hispanic-Latino 72% 0.83 (0.71–0.97) 0.021 0.79 (0.67–0.93) 0.006 
Marital status 
 Married/Domestic Partnered (ref) 80% 1.0 — 1.0 — 
 Divorced/Separated/Widowed 69% 0.56 (0.50–0.62) <0.001 0.74 (0.66–0.83) <0.001 
 Single 69% 0.57 (0.51–0.63) <0.001 0.65 (0.58–0.73) <0.001 
 Unknown 69% 0.57 (0.48–0.69) <0.001 0.66 (0.54–0.80) <0.001 
Census tract-based factors 
nSES (Yost) 
 Group 1 (lowest) (ref) 61% 1.0 — 1.0 — 
 Group 2 73% 1.68 (1.46–1.93) <0.001 1.46 (1.25–1.69) <0.001 
 Group 3 76% 1.97 (1.72–2.26) <0.001 1.63 (1.40–1.89) <0.001 
 Group 4 79% 2.41 (2.11–2.75) <0.001 1.96 (1.70–2.26) <0.001 
 Group 5 (highest) 81% 2.67 (2.36–3.03) <0.001 2.11 (1.83–2.42) <0.001 
RUCA 
 All rural (ref) 74% 1.0 — — — 
 All urban 75% 1.03 (0.88–1.22) 0.707 — — 
 Mostly rural 77% 1.17 (0.94–1.45) 0.167 — — 
 Mostly urban 77% 1.19 (0.99–1.42) 0.063 — — 
Clinical data 
Stage at diagnosis 
 Localized (ref) 57% 1.0 — 1.0 — 
 Regional 79% 2.89 (2.56–3.27) <0.001 2.34 (2.06–2.67) <0.001 
 Distant 69% 1.67 (1.44–1.93) <0.001 1.52 (1.30–1.77) <0.001 
 Unknown/unstaged 69% 1.65 (1.18–2.30) <0.001 1.61 (1.13–2.30) 0.008 
Report source 
 Hospital or Clinic (ref) 76% 1.0 — 1.0 — 
 Other 68% 0.68 (0.56–0.83) <0.001 0.71 (0.58–0.87) 0.001 
First malignancy 
 No (ref) 62% 1.0 — 1.0 — 
 Yes 78% 2.20 (1.98–2.43) <0.001 1.79 (1.60–2.00) <0.001 
VariableLevelCrude HPV prevalenceUnadjusted ORb (95% CI)cPdAdjusted ORb (95% CI)cPd
Demographics 
Sex 
 Female (ref) 62% 1.0 — 1.0 — 
 Male 78% 2.15 (1.94–2.38) <0.001 1.81 (1.62–2.02) <0.001 
Age 
 30–39 years 79% 2.82 (1.71–4.64) <0.001 3.38 (2.01–5.70) <0.001 
 40–49 years 82% 3.44 (2.75–4.30) <0.001 2.96 (2.33–3.76) <0.001 
 50–59 years 78% 2.74 (2.30–3.28) <0.001 2.29 (1.89–2.77) <0.001 
 60–69 years 75% 2.35 (1.97–2.80) <0.001 1.98 (1.64–2.38) <0.001 
 70–79 years 72% 1.98 (1.64–2.38) <0.001 1.84 (1.51–2.24) <0.001 
 80+ years (ref) 57% 1.0 — 1.0 — 
Race 
 White (ref) 77% 1.0 — 1.0 — 
 American Indian/Alaska Native 75% 0.39 (0.34–0.44) <0.001 0.52 (0.44–0.60) <0.001 
 Asian or Pacific Islander 66% 0.58 (0.46–0.72) <0.001 0.53 (0.42-0.67) <0.001 
 Black 57% 0.86 (0.50–1.50) 0.596 0.79 (0.44–1.40) 0.415 
 Unknown 78% 1.05 (0.58–1.90) 0.882 1.27 (0.67–2.39) 0.467 
Ethnicity 
 Non-Spanish-Hispanic-Latino (ref) 76% 1.0 — 1.0 — 
 Spanish-Hispanic-Latino 72% 0.83 (0.71–0.97) 0.021 0.79 (0.67–0.93) 0.006 
Marital status 
 Married/Domestic Partnered (ref) 80% 1.0 — 1.0 — 
 Divorced/Separated/Widowed 69% 0.56 (0.50–0.62) <0.001 0.74 (0.66–0.83) <0.001 
 Single 69% 0.57 (0.51–0.63) <0.001 0.65 (0.58–0.73) <0.001 
 Unknown 69% 0.57 (0.48–0.69) <0.001 0.66 (0.54–0.80) <0.001 
Census tract-based factors 
nSES (Yost) 
 Group 1 (lowest) (ref) 61% 1.0 — 1.0 — 
 Group 2 73% 1.68 (1.46–1.93) <0.001 1.46 (1.25–1.69) <0.001 
 Group 3 76% 1.97 (1.72–2.26) <0.001 1.63 (1.40–1.89) <0.001 
 Group 4 79% 2.41 (2.11–2.75) <0.001 1.96 (1.70–2.26) <0.001 
 Group 5 (highest) 81% 2.67 (2.36–3.03) <0.001 2.11 (1.83–2.42) <0.001 
RUCA 
 All rural (ref) 74% 1.0 — — — 
 All urban 75% 1.03 (0.88–1.22) 0.707 — — 
 Mostly rural 77% 1.17 (0.94–1.45) 0.167 — — 
 Mostly urban 77% 1.19 (0.99–1.42) 0.063 — — 
Clinical data 
Stage at diagnosis 
 Localized (ref) 57% 1.0 — 1.0 — 
 Regional 79% 2.89 (2.56–3.27) <0.001 2.34 (2.06–2.67) <0.001 
 Distant 69% 1.67 (1.44–1.93) <0.001 1.52 (1.30–1.77) <0.001 
 Unknown/unstaged 69% 1.65 (1.18–2.30) <0.001 1.61 (1.13–2.30) 0.008 
Report source 
 Hospital or Clinic (ref) 76% 1.0 — 1.0 — 
 Other 68% 0.68 (0.56–0.83) <0.001 0.71 (0.58–0.87) 0.001 
First malignancy 
 No (ref) 62% 1.0 — 1.0 — 
 Yes 78% 2.20 (1.98–2.43) <0.001 1.79 (1.60–2.00) <0.001 

Note: bold text = statistical significance.

aSEER*Stat Database: Incidence - SEER Research Plus Data (Specialized Head and Neck Fields with Census Tract SES/Rurality), 18 Registries (excl AK), Nov 2020 Sub (2006–2018) <Vintage 2019 Pops by Tract 2010 Geographies> - Linked To Census Tract Attributes SEER 18 excl AK, Census 2010 Geographies Time Dependent (2006–2018; 2008 copied to 06, 07; 2017 copied to 18). National Cancer Institute, DCCPS, Surveillance Research Program; December 2021. http://www.seer.cancer.gov/. Accessed January 24, 2022. (SEER*Stat, Dec 2021, SEER RRID: SCR_003293).

bOdds ratio.

c95% confidence interval.

dα ≤ 0.05.

Table 4.

Odds of HPV-positive non-oropharyngeal squamous cell neoplasm cancer by demographic, socioeconomic, and clinical dataa.

VariableLevelCrude HPV prevalenceUnadjusted ORb (95% CI)cPdAdjusted ORb (95% CI)cPd
Demographics 
Sex 
 Female (ref) 31% 1.0 — — — 
 Male 32% 1.09 (0.84–1.42) 0.528 — — 
Age 
 30–39 years 21% 0.99 (0.35–2.77) 0.983 1.04 (0.36–2.96) 0.947 
 40–49 years 38% 2.36 (1.29–4.32) 0.006 2.31 (1.24–4.31) 0.009 
 50–59 years 35% 2.02 (1.17–3.49) 0.011 1.85 (1.06–3.24) 0.031 
 60–69 years 35% 2.06 (1.20–3.53) 0.009 1.84 (1.06–3.18) 0.030 
 70–79 years 24% 1.23 (0.69–2.18) 0.479 1.12 (0.63–2.00) 0.703 
 80+ years (ref) 21% 1.0 — 1.0 — 
Race 
 White (ref) 35% 1.0 — 1.0 — 
 American Indian/Alaska Native 33% 0.71 (0.50–1.01) 0.054 0.66 (0.47–0.94) 0.023 
 Asian or Pacific Islander 24% 0.57 (0.41–0.78) 0.001 0.53 (0.38–0.73) <0.001 
 Black 28% 0.93 (0.28–3.10) 0.900 0.93 (0.27–3.17) 0.909 
 Unknown 25% 0.62 (0.12–3.07) 0.556 0.41 (0.08–2.12) 0.290 
Ethnicity 
 Non-Spanish-Hispanic-Latino (ref) 32% 1.0 — — — 
 Spanish-Hispanic-Latino 29% 0.87 (0.56–1.36) 0.543 — — 
Marital status 
 Married/Domestic Partnered (ref) 34% 1.0 — — — 
 Divorced/Separated/Widowed 31% 0.86 (0.65–1.15) 0.313 — — 
 Single 29% 0.80 (0.60–1.06) 0.118 — — 
 Unknown 31% 0.88 (0.52–1.47) 0.618 — — 
Census tract-based factors 
nSES (Yost) 
 Group 1 (lowest) (ref) 30% 1.0 — — — 
 Group 2 35% 1.28 (0.89–1.84) 0.188 — — 
 Group 3 34% 1.23 (0.85–1.77) 0.272 — — 
 Group 4 34% 1.19 (0.84–1.70) 0.324 — — 
 Group 5 (highest) 28% 0.90 (0.64–1.28) 0.564 — — 
RUCA 
 All rural (ref) 38% 1.0 — — — 
 All urban 32% 0.76 (0.47–1.20) 0.238 — — 
 Mostly rural 29% 0.67 (0.35–1.28) 0.229 — — 
 Mostly urban 32% 0.75 (0.45–1.25) 0.266 — — 
Clinical data 
Stage at diagnosis 
 Localized (ref) 26% 1.0 — 1.0 — 
 Regional 34% 1.47 (1.02–2.10) 0.036 1.34 (0.92–1.94) 0.123 
 Distant 31% 1.25 (0.86–1.83) 0.249 1.14 (0.77–1.69) 0.513 
 Unknown/unstaged 0% 1.89 (0.89–4.01) 0.096 1.82 (0.84–3.92) 0.128 
Report source 
 Hospital or Clinic (ref) 32% 1.0 — — — 
 Other 31% 0.96 (0.55–1.67) 0.877 — — 
First malignancy 
 No (ref) 26% 1.0 — 1.0 — 
 Yes 34% 1.42 (1.08–1.88) 0.013 1.32 (0.98–1.77) 0.071 
VariableLevelCrude HPV prevalenceUnadjusted ORb (95% CI)cPdAdjusted ORb (95% CI)cPd
Demographics 
Sex 
 Female (ref) 31% 1.0 — — — 
 Male 32% 1.09 (0.84–1.42) 0.528 — — 
Age 
 30–39 years 21% 0.99 (0.35–2.77) 0.983 1.04 (0.36–2.96) 0.947 
 40–49 years 38% 2.36 (1.29–4.32) 0.006 2.31 (1.24–4.31) 0.009 
 50–59 years 35% 2.02 (1.17–3.49) 0.011 1.85 (1.06–3.24) 0.031 
 60–69 years 35% 2.06 (1.20–3.53) 0.009 1.84 (1.06–3.18) 0.030 
 70–79 years 24% 1.23 (0.69–2.18) 0.479 1.12 (0.63–2.00) 0.703 
 80+ years (ref) 21% 1.0 — 1.0 — 
Race 
 White (ref) 35% 1.0 — 1.0 — 
 American Indian/Alaska Native 33% 0.71 (0.50–1.01) 0.054 0.66 (0.47–0.94) 0.023 
 Asian or Pacific Islander 24% 0.57 (0.41–0.78) 0.001 0.53 (0.38–0.73) <0.001 
 Black 28% 0.93 (0.28–3.10) 0.900 0.93 (0.27–3.17) 0.909 
 Unknown 25% 0.62 (0.12–3.07) 0.556 0.41 (0.08–2.12) 0.290 
Ethnicity 
 Non-Spanish-Hispanic-Latino (ref) 32% 1.0 — — — 
 Spanish-Hispanic-Latino 29% 0.87 (0.56–1.36) 0.543 — — 
Marital status 
 Married/Domestic Partnered (ref) 34% 1.0 — — — 
 Divorced/Separated/Widowed 31% 0.86 (0.65–1.15) 0.313 — — 
 Single 29% 0.80 (0.60–1.06) 0.118 — — 
 Unknown 31% 0.88 (0.52–1.47) 0.618 — — 
Census tract-based factors 
nSES (Yost) 
 Group 1 (lowest) (ref) 30% 1.0 — — — 
 Group 2 35% 1.28 (0.89–1.84) 0.188 — — 
 Group 3 34% 1.23 (0.85–1.77) 0.272 — — 
 Group 4 34% 1.19 (0.84–1.70) 0.324 — — 
 Group 5 (highest) 28% 0.90 (0.64–1.28) 0.564 — — 
RUCA 
 All rural (ref) 38% 1.0 — — — 
 All urban 32% 0.76 (0.47–1.20) 0.238 — — 
 Mostly rural 29% 0.67 (0.35–1.28) 0.229 — — 
 Mostly urban 32% 0.75 (0.45–1.25) 0.266 — — 
Clinical data 
Stage at diagnosis 
 Localized (ref) 26% 1.0 — 1.0 — 
 Regional 34% 1.47 (1.02–2.10) 0.036 1.34 (0.92–1.94) 0.123 
 Distant 31% 1.25 (0.86–1.83) 0.249 1.14 (0.77–1.69) 0.513 
 Unknown/unstaged 0% 1.89 (0.89–4.01) 0.096 1.82 (0.84–3.92) 0.128 
Report source 
 Hospital or Clinic (ref) 32% 1.0 — — — 
 Other 31% 0.96 (0.55–1.67) 0.877 — — 
First malignancy 
 No (ref) 26% 1.0 — 1.0 — 
 Yes 34% 1.42 (1.08–1.88) 0.013 1.32 (0.98–1.77) 0.071 

Note: bold text = statistical significance.

aSEER*Stat Database: Incidence - SEER Research Plus Data (Specialized Head and Neck Fields with Census Tract SES/Rurality), 18 Registries (excl AK), Nov 2020 Sub (2006–2018) <Vintage 2019 Pops by Tract 2010 Geographies> - Linked To Census Tract Attributes SEER 18 excl AK, Census 2010 Geographies Time Dependent (2006–2018; 2008 copied to 06, 07; 2017 copied to 18). National Cancer Institute, DCCPS, Surveillance Research Program; December 2021. http://www.seer.cancer.gov/. Accessed January 24, 2022. (SEER*Stat, Dec 2021, SEER RRID: SCR_003293).

bOdds ratio.

c95% confidence interval.

dα ≤ 0.05.

Both in the univariate and multivariable OCPC models, males, those of the five youngest age groups (i.e., 30–79 year olds), those residing in census tracts of every nSES quintile except the lowest quintile (group 1), those diagnosed with any OCPC site except hypopharynx or gum and other mouth, those diagnosed with any stage at diagnosis except local, and those with first malignancies had significantly higher odds of OCPCs being HPV+. American Indians/Alaska Natives or Asians or Pacific Islanders (vs. Whites), Spanish-Hispanics-Latinos, non-married/partnered cases, and those cases from “other” (vs. hospital or clinic) reporting sources had significantly lower odds of OCPCs being HPV+ compared with their respective referent groups. When investigating as oropharyngeal and non-oropharyngeal subgroups, OCPC risk factors for HPV positivity held for oropharyngeal OCPCs, while age and race were the only significant factors for non-oropharyngeal OCPCs. (Because of the high prevalence of HPV+ OCPCs, ORs should not be interpreted as prevalence ratios as they can overestimate the strength of associations.)

Sociodemographic and clinical factors found to be associated with OCPCs, and more specifically, oropharyngeal OCPCs, being HPV+ were in alignment with previous findings. Earlier investigations concluded that males (9, 10, 12–17, 19, 21–23, 26, 39), non-Hispanics (15–17, 40, 41), Whites (9–13, 15, 17, 21–23, 26, 39, 40, 42, 43), younger ages (13–16, 18, 19, 23, 26, 40), and individuals from higher socioeconomic environments (14, 15, 17, 40) were more likely to have HPV+ OCPCs than HPV− OCPCs. Prior research also determined that HPV+ head and neck cancers were less common among non-White races (14, 19, 20), non-married/single individuals (14, 16, 21, 23), and older ages (14, 15, 19, 22, 23). Clinically, other studies have found that HPV+ OCPCs were first identified at more advanced stages (23, 43) and more likely to occur in the base of tongue/oropharyngeal region (14, 15, 22, 23). Non-oropharyngeal OCPC sites being associated with HPV positivity herein is likely due in part to the (complementary, synergistic) interaction between HPV and tobacco/alcohol use [e.g., predominantly among low SES (10), rural residents (44)] risk factors, often whereby tobacco/alcohol use is an underlying risk factor for HPV-related OCPCs (3, 24, 45). Rurality has also been associated with HPV− OCPCs (46). HPV− OCPCs typically occur more frequently among those with lower SES (10), are diagnosed at earlier stages than HPV+ cases (5, 47), and identified in the oral cavity, mostly the front of the mouth (lips, buccal mucosa, gingiva) where tobacco/alcohol chemical exposure contact is greatest (26, 27).

Findings from the current non-oropharyngeal OCPC analysis were also mostly supported by prior research. Prior studies have also found that younger non-oropharyngeal OCPC cases (<70 years old) were more likely to have HPV+ OCPCs (11, 18, 23). As herein, minority racial status has been shown in the literature to be protective against an HPV+ status with Whites being more likely to have an HPV+ non-oropharyngeal OCPC (11, 18). Sex has generally not been identified as a risk factor for non-oropharyngeal OCPC cases, aligning with our findings (11, 18). Only when studies restricted their investigations to hypopharynx cases did their evidence counter our conclusions—age (no association; ref. 11), race (no association; ref. 18), sex (females: less common; ref. 23). Marital status has been investigated minimally with findings being inconclusive—one showing no association (ref. 11; like the current study) and one finding a positive association (23). Clinically, stage has not been related to HPV status of non-oropharyngeal OCPCs elsewhere either (11, 18).

Findings among sociodemographic and clinical factors being (potentially) associated with OCPC HPV positivity can be different depending on whether HPV status is known or assumed due to anatomic and histologic characteristics of the OCPC when HPV status is unknown. In a prior investigation using SEER datasets which did not contain HPV status, it also observed that HPV-related (vs. non–HPV-related) OCPCs tended to occur among those of younger ages, males, non-Hispanics, Whites, and at advanced stages at diagnosis (48). Without HPV status known and only HPV status anatomically/histologically assumed, the prior study (48) found no association with nSES and that HPV-related OCPCs were more common in non-urban census tracts. Such findings are in direct conflict with our current findings of a positive association between increasing nSES and HPV+ and no association between urban residence and HPV+, among OCPCs in general and more specifically, oropharyngeal OCPCs. Confirmatory HPV testing of OCPCs in the current SEER dataset supply strong confidence in the findings presented herein and reiterate the need for future studies to rely on datasets with known HPV status.

This study was strengthened by the utilization of SEER Program data, which covers a substantial racially and ethnically diverse U.S. population that aids in the generalizability of this investigation's findings to other U.S. OCPCs (49). OCPC HPV status was reported in the dataset, avoiding the need to make HPV status conjectures based on OCPC anatomic and/or histologic characteristics. Inclusion of all SEER-identified OCPCs and further refinement of analyses to oropharyngeal and non-oropharyngeal OCPC sites independently allowed for a more complete investigation on the topic of OCPC risk factors for HPV positivity. Residential census tracts of SEER cases were linked to various nSES measures, including primary SES elements of income, education, and employment. The Yost index, used to assess nSES, has face validity, is clear and easily replicated (50), and, compared with county-level data, census tract-level data more closely align with neighborhoods, allowing for more precise nSES estimations [SEER*Stat, SEER Head and Neck with HPV Status and Census Tract-level SES/Rurality Combined Database (December 2021), SEER RRID: SCR_003293]. Reliance on RUCA Codes, also census tract-level data (vs. county-level data), for rurality/urbanity determination also contributed to a more complete assessment of investigated OCPC cases.

The primary limitation was the amount of missing and/or unknown HPV status data, especially in the earlier years of the dataset. To ensure the accuracy of our analyses, we restricted cases to the diagnosis years of 2014 to 2017 whereby more than 60% of HPV statuses were known. Assessing missing data, stage at diagnosis was the only factor substantively different between cases with known and unknown HPV status. The multivariable model was adjusted for all known possible confounders, including stage at diagnosis; however, bias might be introduced because of unknown confounder(s). Additional cases were excluded because of missing nSES data, but post hoc analyses indicated there were no significant changes in OCPC HPV positivity associations as a relatively large sample size was maintained.

This study also had other limitations. OCPC HPV positivity was evaluated by sources reporting information to the SEER Program; therefore, there was likely heterogeneity in HPV testing approaches across reporting sources. However, SEER has strict reporting guidelines in place and is known for complete and quality data so there is minimal concern with the validity of the HPV status data herein (49). Another restriction of the dataset is SEER's inclusion of only OCPC (sub)sites believed to be potentially associated with HPV instead of including all OCPC (sub)sites. To maximize use of the dataset and produce the most comprehensive results possible, we included all OCPC (sub)sites included in the dataset in our analyses whereas prior studies further limited their investigations to OCPC anatomic sites of interest (11–13, 16–20, 22, 23). Together, categories of combined anatomic sites of the SEER “Site Recode” variable (e.g., other oral cavity and pharynx, gum, and other mouth) coupled with unknown HPV testing methodologies of these sites potentially created an additional weakness of the dataset. In addition, SEER Program data do not include details concerning comorbidity and patient behavior, including OCPC screening behavior, and this is another limitation of the dataset, impeding the adjustment for confounding and/or effect modification by these factors. nSES and RUCA categories were assessed at one point in time; neighborhood status may vary over time. Furthermore, as area-based measures, it is possible that nSES and RUCA categories may not accurately describe individual SES and rurality/urbanity statuses. Finally, reliance on prevalence ORs with a high prevalence of HPV+ OCPCs suggests that ORs should not be interpreted as prevalence ratios. [Prevalence values are contained within (prevalence) ORs and prevalence ratios and therefore would be inaccurate to interpret in place of ORs or prevalence ratios.] In such situations of high disease prevalence, (prevalence) odds ratios can overestimate the prevalence ratio by at least 10%, especially when the intent is to approximate a risk ratio by an OR (51).

In conclusion, there are distinct sociodemographic and clinical differences in OCPC cases by HPV status, especially among oropharyngeal OCPCs, but less so among non-oropharyngeal OCPCs. An accurate understanding of epidemiologic differences in OCPC cases, namely oropharyngeal OCPCs, by HPV status is imperative so that high-risk populations are not ignored and OCPC incidence and mortality rates do not go unchecked. An HPV vaccine exists and can prevent much of the oral HPV infections (i.e., 88%) that cause OCPCs (52), but the prevention measure is underutilized (i.e., Healthy People 2020 goal for 13–17 years old adolescents: 80% vaccinated; United States: 59% up-to-date; refs. 53–55). Given that HPV-related OCPCs tend to occur later in life, and we have not reached herd immunity with current vaccination rates, HPV-related OCPCs may continue to increase until at least 2050, assuming HPV vaccination rates can be rapidly scaled up between now and then (56). In the meantime, health care providers and public health professionals can use this (oropharyngeal) OCPC risk factor for HPV positivity information to more accurately and effectively direct their messages and impact behaviors to reduce the OCPC risk within specific populations. For instance, public health efforts to minimize the oropharyngeal OCPC burden may target males, Whites, younger ages (especially 30–59 year olds), higher nSES groups, and those at higher risk for tonsil and base of tongue OCPCs while also encouraging more frequent OCPC screenings to prevent primary malignancies and later stage diagnoses. They should also consider intervening with awareness messages with Whites and/or those who are 40 to 69 years old to keep non-oropharyngeal OCPCs from increasing significantly. Health care providers can formulate patient risk assessments by evaluating significant sociodemographic and clinical risk factors identified in this study within their patient populations. And although a preliminary finding, clinicians may want to consider testing non-oropharyngeal (e.g., nasopharynx) OCPCs for HPV instead of assuming them to be HPV− as traditionally assumed, especially if patients exhibit additional OCPC risk factors for HPV positivity identified herein (i.e., 40–69 years old). Collectively, public health targeting and provider tailoring of educational efforts on these [(non)oropharyngeal] OCPC risk factors for HPV positivity may work together with repetitive consistent messaging at both the community and individual levels to maximize our potential efforts at preventing these cancers from reaching epidemic proportions in the United States.

K.H. Jordan reports grants from NCI/NIH during the conduct of the study. E.D. Paskett reports grants from Merck Foundation, Pfizer, Genentech, and Guardant Health outside the submitted work. No disclosures were reported by the other authors.

K.H. Jordan: Conceptualization, funding acquisition, investigation, visualization, methodology, writing–original draft, project administration, writing–review and editing, and data interpretation. J.L. Fisher: Conceptualization, data curation, formal analysis, investigation, methodology, writing–original draft, writing–review and editing, and data interpretation. J.A. Stephens: Data curation, formal analysis, investigation, methodology, writing–review and editing, and data interpretation. E.D. Paskett: Conceptualization, resources, supervision, funding acquisition, writing–review and editing.

The NIH's NCI generously supported this research (P01 CA229143-S1 for K.H. Jordan), but the funders did not have any other role in this study.

The publication costs of this article were defrayed in part by the payment of publication fees. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

Note: Supplementary data for this article are available at Cancer Epidemiology, Biomarkers & Prevention Online (http://cebp.aacrjournals.org/).

1.
American Cancer Society
.
Cancer facts & figures 2022
;
2022
.
Available from
: https://www.cancer.org/research/cancer-facts-statistics/all-cancer-facts-figures/cancer-facts-figures-2022.html.
2.
Franceschi
S
,
Muñoz
N
,
Bosch
XF
,
Snijders
PJ
,
Walboomers
JM
.
Human papillomavirus and cancers of the upper aerodigestive tract: a review of epidemiological and experimental evidence
.
Cancer Epidemiol Biomarkers Prev
1996
;
5
:
567
75
.
3.
Gillison
ML
.
Human papillomavirus-associated head and neck cancer is a distinct epidemiologic, clinical, and molecular entity
.
Semin Oncol
2004
;
31
:
744
54
.
4.
Gillison
ML
,
Chaturvedi
AK
,
Anderson
WF
,
Fakhry
C
.
Epidemiology of human papillomavirus-positive head and neck squamous cell carcinoma
.
J Clin Oncol
2015
;
33
:
3235
42
.
5.
Candotto
V
,
Lauritano
D
,
Nardone
M
,
Baggi
L
,
Arcuri
C
,
Gatto
R
, et al
.
HPV infection in the oral cavity: epidemiology, clinical manifestations and relationship with oral cancer
.
Oral Implantol
2017
;
10
:
209
20
.
6.
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
:
299
309
.
7.
American Cancer Society
.
Cancer facts & figures 2021
;
2021
.
Available from
: https://www.cancer.org/research/cancer-facts-statistics/all-cancer-facts-figures/cancer-facts-figures-2021.html.
8.
Howlader
NNA
,
Krapcho
M
,
Miller
D
,
Brest
A
,
Yu
M
,
Ruhl
J
, et al
.
Data from: SEER cancer statistics review, 1975–2018
.
April 2021
;
Based on November 2020 SEER data submission
.
Bethesda, MD
.
Deposited April 2021
.
9.
Javadi
P
,
Sharma
A
,
Zahnd
WE
,
Jenkins
WD
.
Evolving disparities in the epidemiology of oral cavity and oropharyngeal cancers
.
Cancer Causes Control
2017
;
28
:
635
45
.
10.
Dahlstrom
KR
,
Bell
D
,
Hanby
D
,
Li
G
,
Wang
LE
,
Wei
Q
, et al
.
Socioeconomic characteristics of patients with oropharyngeal carcinoma according to tumor HPV status, patient smoking status, and sexual behavior
.
Oral Oncol
2015
;
51
:
832
8
.
11.
Marshall
DC
,
Kao
DD
,
Bakst
R
,
Sharma
S
,
Ferrandino
R
,
Rosenzweig
K
, et al
.
Prognostic role of human papilloma virus status in hypopharyngeal squamous cell carcinoma
.
Laryngoscope Investig Otolaryngol
2020
;
5
:
860
7
.
12.
Li
Y
,
Jiang
B
,
Huang
M
,
Wang
G
,
Bai
Y
,
Huang
W
.
Retrospective study of the clinicopathological characteristics and prognosis of elderly patients with oropharyngeal squamous cell carcinoma
.
J Int Med Res
2021
;
49
:
3000605211016662
.
13.
Pike
LRG
,
Hwang
WL
,
Royce
TJ
,
Sanford
NN
,
Mahal
BA
.
HPV status predicts for improved survival following chemotherapy in metastatic squamous cell carcinoma of the oropharynx
.
Oral Oncol
2018
;
86
:
69
74
.
14.
Pike
LRG
,
Royce
TJ
,
Mahal
AR
,
Kim
DW
,
Hwang
WL
,
Mahal
BA
, et al
.
Outcomes of HPV-associated squamous cell carcinoma of the head and neck: impact of race and socioeconomic status
.
J Natl Compr Canc Netw
2020
;
18
:
177
84
.
15.
Mahal
BA
,
Catalano
PJ
,
Haddad
RI
,
Hanna
GJ
,
Kass
JI
,
Schoenfeld
JD
, et al
.
Incidence and demographic burden of HPV-associated oropharyngeal head and neck cancers in the United States
.
Cancer Epidemiol Biomarkers Prev
2019
;
28
:
1660
7
.
16.
Tewari
P
,
Kashdan
E
,
Walsh
C
,
Martin
CM
,
Parnell
AC
,
O'Leary
JJ
.
Estimating the conditional probability of developing human papilloma virus related oropharyngeal cancer by combining machine learning and inverse Bayesian modelling
.
PLoS Comput Biol
2021
;
17
:
e1009289
.
17.
Sanford
NN
,
Hwang
WL
,
Pike
LRG
,
Lam
AC
,
Royce
TJ
,
Mahal
BA
.
Trimodality therapy for HPV-positive oropharyngeal cancer: a population-based study: trimodality therapy for HPV+ OPC
.
Oral Oncol
2019
;
98
:
28
34
.
18.
Abdel-Rahman
O
.
Prognostic value of HPV status among patients with hypopharyngeal carcinoma: a population-based study
.
Clin Transl Oncol
2020
;
22
:
1645
50
.
19.
Abdel-Rahman
O
.
Sex-based differences in outcomes of patients with HPV-positive versus HPV-negative oropharyngeal carcinoma: a population-based study
.
Expert Rev Anticancer Ther
2020
;
20
:
813
8
.
20.
Wotman
M
,
Oh
EJ
,
Ahn
S
,
Kraus
D
,
Costantino
P
,
Tham
T
.
HPV status in patients with nasopharyngeal carcinoma in the United States: a SEER database study
.
Am J Otolaryngol
2019
;
40
:
705
10
.
21.
Wang
J
,
Tian
Y
,
Huang
H
,
Huang
D
,
Liu
Y
,
Tian
Y
, et al
.
The prognosis of HPV-associated metastatic pharyngeal patients by primary and distant site
.
Oral Oncol
2022
;
125
:
105675
.
22.
Wang
P
,
Zang
S
,
Zhang
M
,
Li
G
,
Wang
Y
,
Qiao
Q
.
Impact of HPV status on metastatic patterns and survival in non-oropharyngeal head and neck cancer with distant metastasis
.
Eur Arch Otorhinolaryngol
2022
;
279
:
3029
39
.
23.
Wu
Q
,
Wang
M
,
Liu
Y
,
Wang
X
,
Li
Y
,
Hu
X
, et al
.
HPV positive status is a favorable prognostic factor in non-nasopharyngeal head and neck squamous cell carcinoma patients: a retrospective study from the surveillance, epidemiology, and end results database
.
Front Oncol
2021
;
11
:
688615
.
24.
Saraiya
M
,
Unger
ER
,
Thompson
TD
,
Lynch
CF
,
Hernandez
BY
,
Lyu
CW
, et al
.
US assessment of HPV types in cancers: implications for current and 9-valent HPV vaccines
.
J Natl Cancer Inst
2015
;
107
:
djv086
.
25.
Berman
TA
,
Schiller
JT
.
Human papillomavirus in cervical cancer and oropharyngeal cancer: one cause, two diseases
.
Cancer
2017
;
123
:
2219
29
.
26.
Chaturvedi
AK
,
Engels
EA
,
Anderson
WF
,
Gillison
ML
.
Incidence trends for human papillomavirus-related and -unrelated oral squamous cell carcinomas in the United States
.
J Clin Oncol
2008
;
26
:
612
9
.
27.
Trotta
BM
,
Pease
CS
,
Rasamny
JJ
,
Raghavan
P
,
Mukherjee
S
.
Oral cavity and oropharyngeal squamous cell cancer: key imaging findings for staging and treatment planning
.
Radiographics
2011
;
31
:
339
54
.
28.
Kumar
R
,
Rai
AK
,
Das
D
,
Das
R
,
Kumar
RS
,
Sarma
A
, et al
.
Alcohol and tobacco increases risk of high risk HPV infection in head and neck cancer patients: study from North-East Region of India
.
PLoS One
2015
;
10
:
e0140700
.
29.
Smith
EM
,
Rubenstein
LM
,
Haugen
TH
,
Hamsikova
E
,
Turek
LP
.
Tobacco and alcohol use increases the risk of both HPV-associated and HPV-independent head and neck cancers
.
Cancer Causes Control
2010
;
21
:
1369
78
.
30.
Johnson
DE
,
Burtness
B
,
Leemans
CR
,
Lui
VWY
,
Bauman
JE
,
Grandis
JR
.
Head and neck squamous cell carcinoma
.
Nat Rev Dis Primers
2020
;
6
:
92
.
31.
Denson
L
,
Janitz
AE
,
Brame
LS
,
Campbell
JE
.
Oral cavity and oropharyngeal cancer: changing trends in incidence in the United States and Oklahoma
.
J Okla State Med Assoc
2016
;
109
:
339
45
.
32.
Amiling
R
,
Meites
E
,
Querec
TD
,
Stone
L
,
Singh
V
,
Unger
ER
, et al
.
Juvenile-onset recurrent respiratory papillomatosis in the United States, epidemiology and HPV types-2015–2020
.
J Pediatric Infect Dis Soc
2021
;
10
:
774
81
.
33.
Yu
M
,
Tatalovich
Z
,
Gibson
JT
,
Cronin
KA
.
Using a composite index of socioeconomic status to investigate health disparities while protecting the confidentiality of cancer registry data
.
Cancer Causes Control
2014
;
25
:
81
92
.
34.
Yost
K
,
Perkins
C
,
Cohen
R
,
Morris
C
,
Wright
W
.
Socioeconomic status and breast cancer incidence in California for different race/ethnic groups
.
Cancer Causes Control
2001
;
12
:
703
11
.
35.
Liu
L
,
Deapen
D
,
Bernstein
L
.
Socioeconomic status and cancers of the female breast and reproductive organs: a comparison across racial/ethnic populations in Los Angeles County, California (United States)
.
Cancer Causes Control
1998
;
9
:
369
80
.
36.
Cromartie
J
,
ERS USDA
.
Rural-urban commuting area codes
.
Available from
: https://www.ers.usda.gov/data-products/rural-urban-commuting-area-codes.aspx.
37.
World Health Organization
.
International classification of diseases for oncology (ICD-O), 3rd ed., 1st revision ed
.
Geneva, Switzerland
:
World Health Organization
;
2013
.
39.
Bean
MB
,
Switchenko
JM
,
Steuer
CE
,
Patel
M
,
Higgins
K
,
McDonald
M
, et al
.
Demographic and socioeconomic factors associated with metastases at presentation in HPV-related squamous cell carcinoma of the head and neck: an NCDB analysis
.
JCO Oncol Pract
2020
;
16
:
e476
87
.
40.
Peterson
CE
,
Khosla
S
,
Jefferson
GD
,
Davis
FG
,
Fitzgibbon
ML
,
Freels
S
, et al
.
Measures of economic advantage associated with HPV-positive head and neck cancers among non-Hispanic black and white males identified through the national cancer database
.
Cancer Epidemiol
2017
;
48
:
1
7
.
41.
Benard
VB
,
Johnson
CJ
,
Thompson
TD
,
Roland
KB
,
Lai
SM
,
Cokkinides
V
, et al
.
Examining the association between socioeconomic status and potential human papillomavirus-associated cancers
.
Cancer
2008
;
113
:
2910
8
.
42.
Freedman
RL
,
Sibley
H
,
Williams
AM
,
Chang
SS
.
Race, not socioeconomic disparities, correlates with survival in human papillomavirus-negative oropharyngeal cancer: a retrospective study
.
Am J Otolaryngol
2021
;
42
:
102816
.
43.
Rotsides
JM
,
Oliver
JR
,
Moses
LE
,
Tam
M
,
Li
Z
,
Schreiber
D
, et al
.
Socioeconomic and racial disparities and survival of human papillomavirus-associated oropharyngeal squamous cell carcinoma
.
Otolaryngol Head Neck Surg
2021
;
164
:
131
8
.
44.
CDC
.
Behavioral risk factor surveillance survey, 2019. Analysis by the american lung association epidemiology and statistics unit using SPSS software
.
45.
Auguste
A
,
Deloumeaux
J
,
Joachim
C
,
Gaete
S
,
Michineau
L
,
Herrmann-Storck
C
, et al
.
Joint effect of tobacco, alcohol, and oral HPV infection on head and neck cancer risk in the French West Indies
.
Cancer Med
2020
;
9
:
6854
63
.
46.
Marks
JA
,
Switchenko
JM
,
Steuer
CE
,
Ryan
M
,
Patel
MR
,
McDonald
MW
, et al
.
Socioeconomic factors influence the impact of tumor HPV status on outcome of patients with oropharyngeal squamous cell carcinoma
.
JCO Oncol Pract
2021
;
17
:
e313
22
.
47.
Carpén
T
,
Sjöblom
A
,
Lundberg
M
,
Haglund
C
,
Markkola
A
,
Syrjänen
S
, et al
.
Presenting symptoms and clinical findings in HPV-positive and HPV-negative oropharyngeal cancer patients
.
Acta Otolaryngol
2018
;
138
:
513
8
.
48.
Jordan
KH
,
Fisher
JL
,
Paskett
ED
.
Distinct sociodemographic differences in incidence and survival rates for anatomically human papillomavirus (HPV)-related, non-HPV-related, and “other” oral cavity and pharynx cancers: an analysis of Surveillance, Epidemiology and End Results (SEER) Program data
.
Front Oncol
2022
;
12
:
980900
.
49.
Duggan
MA
,
Anderson
WF
,
Altekruse
S
,
Penberthy
L
,
Sherman
ME
.
The surveillance, epidemiology, and end results (SEER) program and pathology: toward strengthening the critical relationship
.
Am J Surg Pathol
2016
;
40
:
e94
e102
.
50.
Boscoe
FP
,
Liu
B
,
Lee
F
.
A comparison of two neighborhood-level socioeconomic indexes in the United States
.
Spat Spatiotemporal Epidemiol
2021
;
37
:
100412
.
51.
Tamhane
AR
,
Westfall
AO
,
Burkholder
GA
,
Cutter
GR
.
Prevalence odds ratio versus prevalence ratio: choice comes with consequences
.
Stat Med
2016
;
35
:
5730
5
.
52.
Chaturvedi
AK
,
Graubard
BI
,
Broutian
T
,
Pickard
RKL
,
Tong
ZY
,
Xiao
W
, et al
.
Effect of prophylactic human papillomavirus (HPV) vaccination on oral HPV infections among young adults in the United States
.
J Clin Oncol
2018
;
36
:
262
7
.
53.
NIH CDC
.
State cancer profiles
.
Available from
: https://www.statecancerprofiles.cancer.gov/index.html.
54.
CDC
.
Human papillomavirus (HPV) vaccination coverage among adolescents 13–17 years by State, HHS Region, and the United States, National Immunization Survey-Teen (NIS-Teen)
;
2018
.
Available from
: https://www.cdc.gov/vaccines/imz-managers/coverage/teenvaxview/index.html.
55.
Walker
TY
,
Elam-Evans
LD
,
Yankey
D
,
Markowitz
LE
,
Williams
CL
,
Mbaeyi
SA
, et al
.
National, regional, state, and selected local area vaccination coverage among adolescents aged 13–17 Years - United States, 2017
.
MMWR Morb Mortal Wkly Rep
2018
;
67
:
909
17
.
56.
Simms
KT
,
Steinberg
J
,
Caruana
M
,
Smith
MA
,
Lew
JB
,
Soerjomataram
I
, et al
.
Impact of scaled up human papillomavirus vaccination and cervical screening and the potential for global elimination of cervical cancer in 181 countries, 2020–99: a modelling study
.
Lancet Oncol
2019
;
20
:
394
407
.
This open access article is distributed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license.