Geographical disparities in human papillomavirus herd protection

Abstract Background Human papillomavirus (HPV) vaccination has occurred unequally across the United States, potentially contributing to uneven vaccine‐type HPV prevalence between regions. We examined whether emerging vaccine‐related herd protection exhibits regional differences among unvaccinated girls and women. Methods We evaluated the prevalence of vaginal HPV among women 14‐59 years of age from 2003 to 2014 using repeated cross‐sectional data from the National Health and Nutrition Examination Survey (NHANES). Women who provided an adequate vaginal swab sample were included. Vaginal prevalence of vaccine‐type HPV (types 6, 11, 16, 18) were examined in four regions of the United States between 2003 and 2014. We examined vaccine‐type HPV prevalence in 2007‐2014 in each US census region among younger participants (14‐34 years old) stratified by vaccination status to determine whether one or both groups contributed to uneven HPV prevalence. Results A total of 12 175 participants 14‐59 years of age met inclusion criteria. Vaccine‐type HPV prevalence decreased in all regions. Vaccine‐type HPV varied by region only among unvaccinated 14‐34 year olds, with a higher prevalence in the Midwest (13.8%, 95% confidence interval (CI): 10.7‐17.0) and South (12.5%, 95% CI: 10.2‐14.8) compared to the Northeast (8.9%, 95% CI: 6.5‐11.2). No regional variation in vaccine‐type HPV prevalence was observed among vaccinated participants. Conclusions Higher prevalence of vaccine‐type HPV among unvaccinated women in the South and Midwest may contribute to regional disparities in HPV‐related cancer incidence, as emerging herd immunity may not be as strong in those regions.


| 5273
BERENSON Et al. 11-12 years of age, 4 but as of 2018, only 68% of [13][14][15][16][17] year olds had initiated the HPV vaccine series, and only 51% had completed it. 5 National surveillance studies indicate that HPV vaccine initiation and completion rates among females are lower in the South and Midwest compared to other regions of the United States. 3,6,7 This regional inequality among girls and women is a major concern because cervical cancer rates are higher in the South and in some Midwestern states. 8 Continued low vaccination rates in these regions could further contribute to geographic disparities in HPV infection rates and HPVrelated cancers.
Recently, regional variation in HPV prevalence following the launch of the HPV vaccine in the United States in 2006 was described, with a higher prevalence of vaccine-type HPV occurring in the South and Midwest than in the West and Northeast in post-licensure years according to data from the National Health and Nutrition Examination Survey (NHANES). 9 Additionally, herd immunity is emerging in the US population. These studies demonstrated reduced HPV, genital wart, and precancerous cervical lesion prevalence among unvaccinated groups eligible to receive the vaccine in the United States. [10][11][12][13][14] Age groups not eligible, or that may have already been exposed before catch-up vaccination could occur, did not exhibit similar decreases in HPV-related disease which indicates that HPV vaccination is likely responsible for the emerging herd immunity. 10 Given these findings and prior reports of regional disparities in HPV vaccination, it is likely that emergent herd protection is also geographically heterogeneous. Therefore, we sought to examine variation in HPV prevalence among vaccinated and unvaccinated women in four regions of the United States during post-licensure years (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014) and explore the distribution of herd protection among unvaccinated girls and women in each region using NHANES data.

| MATERIALS AND METHODS
For this study, we utilized the NHANES dataset to conduct a secondary data analysis. NHANES consists of a series of nationally representative, cross-sectional surveys of the noninstitutionalized US population, conducted in two-year cycles by the US Centers for Disease Control and Prevention. 15 Detailed methods have been previously described. 16 Briefly, NHANES participants complete a household survey and undergo physical exams in a mobile examination center. Females 14-59 years of age were also asked to self-collect a vaginal swab for HPV testing. All NHANES swab samples were extracted and tested for types of HPV DNAusing Linear Array HPV genotypying assay (Roche Diagnostics). This data analysis, conducted using NHANES data, was exempted from review by the University of Texas Medical Branch Institutional Review Board. Restricted geographic variables were accessed through the Research Data Center at the National Center for Health Statistics.
The prevalence of HPV by region for six biennial cycles was determined between 2003 and 2014. Inclusion in this study required the following: (a) female, (b) 14-59 years of age, (c) adequate vaginal swab sample, (d) state of residence, and (e) self-reported HPV vaccination status for those participating 2007-2014. HPV results were grouped by type: any of 37 HPV types, vaccine-types (6,11,16,18), and non-vaccine types (any HPV types except vaccine-types). Each HPV type group was transformed into a binary variable, with 0 = no type in the group detected and 1 = at least one type in the group detected.
Participants before the 2007-2008 cycle were considered unvaccinated, as the vaccine was approved for sale in mid-2006. Among those who responded to a question asking if they had received at least 1 dose of the HPV vaccine, the proportion of participants who responded "yes" was calculated. Demographics included race/ethnicity, marital status, and health insurance coverage. Behavioral characteristics included cigarette use, sexual behaviors, and alcohol and marijuana use in the past 12 months among 20-59 year olds. History of gonorrhea and chlamydia were based on self-report while history of herpes simplex type 2 (HSV2) was based on lab-confirmed testing and self-report.
The region of residence for participants in the NHANES dataset was defined by four US census regions, including the Northeast, Midwest, South, and West. Bivariate comparisons were conducted using Rao-Scott chi-square tests. For 2007-2008 through 2013-2014 cycles, analyses were restricted to participants who responded to whether they had received at least 1 HPV vaccine dose. After examining comparisons among 14-59 year olds, we grouped the sample into younger (14-34-year-old) and older (35-59-year-old) participants. We included older participants who would not have been eligible for HPV vaccination during the periods examined to determine whether observations were similar or different in both groups. We anticipated that women in the 35 -59 year old group would not exhibit signs of herd immunity, as most would not have received HPV vaccine doses, and would therefore have different results than the younger age group. Unadjusted logistic regression models were used to compare HPV population adjusted odds ratios (PaORs), which were weighted to produce national estimates as well as adjusted for the complex survey design. Each model was stratified by region. A multivariable logistic regression model was used to calculate adjusted PaORs (aPaORs) for vaccine-type HPV during post-licensure cycles, after adjusting for demographics and HPV risk  factors. MEC weights were used to account for nonresponse and complex survey methods. 17 Analyses on the NHANES dataset were conducted using SAS statistical software version 9.4 (SAS Institute, Inc).

| RESULTS
A total of 14 033 females 14-59 years of age participated across the 6 cycles of NHANES sampling from 2003 to 2014.
Of these, 12 182 (86.8%) had valid HPV genotyping results. 7 After excluding 7 participants without known region of residence, 12 175 participants were eligible for final analyses. For post-vaccine licensure cycles (2007-2014), 8025 girls and women 14-59 years old reported their HPV vaccination status. There were regional differences in several demographic characteristics, including age group, race/ethnicity, and health insurance coverage (Table 1). Regional differences were also observed in behavioral characteristics including cigarette use,    Figure 1B). Among 14-59-year-old participants in the post-vaccine licensure cycles, we found significant regional differences in the prevalence of non-vaccine-type HPV (Table 2). However, after stratifying by vaccination status, there were significant regional differences in HPV prevalence among unvaccinated participants only. Regional differences among the unvaccinated were observed for the prevalence of any type, vaccine type, and non-vaccine type HPV. No regional differences were observed among vaccinated participants in the post-licensure cycles.
To determine whether emerging herd immunity is apparent, we examined the prevalence of HPV between vaccinated and unvaccinated women by region. Among 14-34 year olds in post-licensure cycles, we observed regional differences in the prevalence of vaccine-type HPV, but not the other groupings of HPV types (Table 3). HPV vaccination among the 14-34 year olds increased from 10.9% in 2007-2008 to 34.9% in 2013-2014, with the greatest increase occurring in the Northeast. However, the vaccination rate within each cycle and in all cycles combined did not vary by US region (Table S2). Conversely, among women 35-59 years of age, of whom only 1.55% have received the HPV vaccine during 2007-2014, we observed regional differences for any type, for high-risk type, for high-risk non-vaccine type, and non-vaccine type HPV in the post-licensure cycles (Table S1).
We next sought to determine if emerging herd protection, indicated by declines in vaccine-type HPV prevalence among the unvaccinated, varies across the United States. We found no regional differences in HPV prevalence for any grouping of HPV types among vaccinated 14-34-year-old participants ( Table 3). Among unvaccinated 14-34 year olds, however, we observed regional variations in vaccine-type HPV, which was significantly higher among unvaccinated young women compared to vaccinated young women in every region. Among Northeastern participants, those who were unvaccinated had lower prevalence rates of any type HPV compared to those who were vaccinated.
To determine whether these findings were the result of confounding factors, we used multivariable models to calculate aPaORs for demographics and behavioral characteristics among unvaccinated 14-to 34-year-old women. Vaccine-type HPV prevalence in unvaccinated participants was lower in the  Table 4). In a model that included only region and number of lifetime same sex partners, we found that the time effects were magnified compared to the model adjusted with other confounders. (Table  S2).

| DISCUSSION
Human papillomavirus vaccination has been shown to vary across the United States, with the lowest HPV vaccination rates occurring in the South. 7,18 This study confirms that geographic disparities in HPV vaccination may be problematic, as there were significant geographic differences in vaccinetype HPV prevalence among unvaccinated women. Regions with lower vaccination rates had the highest prevalence of vaccine-type HPV among 14-to 34-year-old unvaccinated women compared to vaccinated women. Similar findings have been reported in a Switzerland-based modeling study, in which higher HPV type 16 prevalence was associated with regions of lower vaccine uptake. 19 While this effect may be somewhat mitigated by sexual activity occurring between occupants in different regions of a geographically small nation such as Switzerland, in the larger United States, the majority of residents tend to remain within a single census region. Moves between regions are relatively rare, particularly among populations from lower socioeconomic status groups. 20 Thus, uneven vaccination rates could contribute to continuing geographic disparities in HPV-related cancers. [21][22][23] Our study adds further evidence that herd immunity is emerging among unvaccinated women. A reduction in vaccine-type HPV prevalence was seen among unvaccinated young women across time in the post-licensure years, even after controlling for several factors associated with HPV acquisition and persistence. Our findings fit with several  prior studies identifying a decrease in HPV prevalence among unvaccinated people in the United States and other countries including Australia and Scotland. [24][25][26] Although herd immunity in the United States has previously been reported, regional variations were not. 14 Understanding emerging herd immunity effects at the regional or even local level is important as seemingly strong protection at the national level may belie vulnerability in populations experiencing weaker protection. Recently, clustered outbreaks of measles and pertussis in localized areas with low vaccination rates have caused concerns about how they contribute to perpetuation of these diseases. [27][28][29][30] It may take more than a decade for the consequences of disparities in HPV immunity to become apparent due to the slow progression from HPV infection to cancer. Nonetheless, the risk of higher regional HPV-related cancer morbidity and mortality warrants greater attention to regional disparities in HPV prevalence, particularly for regions already presenting with greater incidence of HPV-related cancers.
The primary strength of this study is that it used repeated cross-sectional survey data representative of the United States. However, these data also create limitation as they may not be representative of each region. Moreover data on HPV vaccination in NHANES are collected by self-report, which is subject to recall bias. Furthermore, accuracy of adolescent HPV vaccination reports have been found to vary by race/ethnicity, and could potentially vary by region, which could introduce bias within and between regions. 31 It is possible that lower vaccine-type HPV in Northeast and West regions may not be due to herd protection related to higher HPV vaccination rates, but rather due to overall lower circulation of HPV in those regions as a result of differing risk behaviors or other factors.
In summary, this study sheds light on the origin of recently reported regional variation in HPV prevalence in the T A B L E 3 Human papillomavirus (HPV) prevalence detected through vaginal swab by region and by vaccination status among women aged United States. Only among unvaccinated girls and women did vaccine-type HPV prevalence vary by region. Notably, odds of vaccine-type HPV prevalence were higher among unvaccinated young women in the Midwest and South where vaccination rates are not increasing as quickly as in the Northeast. Given the higher prevalence of HPV in the South and Midwest, interventions are needed to increase vaccination rates and to maintain or improve cervical cancer screening in these areas.