A Population-Based Study of Sociodemographic and Geographic Variation in HPV Vaccination

Background: Human papillomavirus (HPV) vaccination rates in the United States remain low and lag behind other recommended adolescent vaccines. Studies evaluating the association of geographic and area-level characteristics with HPV vaccination rates provide a valuable resource for public health planning. Method: We used the Rochester Epidemiology Project data linkage system to ascertain HPV vaccination rates between 2010 and 2015 in a 7-county region of southern Minnesota. Geocoded individual patient data were spatially linked to socioeconomic data from the American Community Survey at the census block group level. Bayesian hierarchical logistic regression was used to model incident vaccination rates, adjusting for individual- and area-level sociodemographic characteristics, and geolocation. Geolocation was modeled as an approximated Gaussian field using a Stochastic Partial Differential Equations approach. All models were estimated using Integrated Nested Laplace Approximations. Results: In adjusted models, increasing age and female sex were associated with increased HPV vaccination. Lower socioeconomic status was associated with decreased rates of initiation [adjusted odds ratio (AOR); 95% confidence interval = 0.90 (0.86–0.95)], completion of the second dose [AOR = 0.88 (0.83–0.93)], and completion of the third dose [AOR = 0.85 (0.80–0.92)]. Geographic spatial analysis demonstrated increased odds of vaccination for the eastern region and in the greater Rochester metropolitan area, showing significant spatial variation not explained by individual level characteristics and ACS block group-level data. Conclusions: HPV vaccination rates varied geographically and by individual and geographically indexed sociodemographic characteristics. Impact: Identifying geographic regions with low HPV vaccination rates can help target clinical and community efforts to improve vaccination rates. Cancer Epidemiol Biomarkers Prev; 26(4); 533–40. ©2017 AACR. See all the articles in this CEBP Focus section, “Geospatial Approaches to Cancer Control and Population Sciences.”

[1]  S. Jacobsen,et al.  Correlates for human papillomavirus vaccination of adolescent girls and young women in a managed care organization. , 2010, American journal of epidemiology.

[2]  Arthur F. Pillsbury,et al.  The Salinity of Rivers , 1981 .

[3]  H. Rue,et al.  An explicit link between Gaussian fields and Gaussian Markov random fields; The SPDE approach , 2010 .

[4]  Jenny Jeyarajah,et al.  Human Papillomavirus Vaccination Coverage Among Adolescents, 2007–2013, and Postlicensure Vaccine Safety Monitoring, 2006–2014 — United States , 2014, MMWR. Morbidity and mortality weekly report.

[5]  S. Stokley,et al.  National, Regional, State, and Selected Local Area Vaccination Coverage Among Adolescents Aged 13–17 Years — United States, 2018 , 2019, MMWR. Morbidity and mortality weekly report.

[6]  A. Giuliano,et al.  Impact of baseline covariates on the immunogenicity of a quadrivalent (types 6, 11, 16, and 18) human papillomavirus virus-like-particle vaccine. , 2007, The Journal of infectious diseases.

[7]  A. D. Diez Roux,et al.  Neighborhoods and health , 2010, Annals of the New York Academy of Sciences.

[8]  Finn Lindgren,et al.  Bayesian Spatial Modelling with R-INLA , 2015 .

[9]  M. Schootman,et al.  Geographic disparity, area poverty, and human papillomavirus vaccination. , 2010, American journal of preventive medicine.

[10]  S. Stokley,et al.  National, Regional, State, and Selected Local Area Vaccination Coverage Among Adolescents Aged 13–17 Years — United States, 2014 , 2015, MMWR. Morbidity and mortality weekly report.

[11]  H. Lawson,et al.  Quadrivalent Human Papillomavirus Vaccine: Recommendations of the Advisory Committee on Immunization Practices (ACIP). , 2007, MMWR. Recommendations and reports : Morbidity and mortality weekly report. Recommendations and reports.

[12]  P. Diggle Applied Spatial Statistics for Public Health Data , 2005 .

[13]  J. Snow On the Mode of Communication of Cholera , 1856, Edinburgh medical journal.

[14]  X. Castellsagué,et al.  Comparison of the Immunogenicity and Reactogenicity of a Prophylactic Quadrivalent Human Papillomavirus (Types 6, 11, 16, and 18) L1 Virus-Like Particle Vaccine in Male and Female Adolescents and Young Adult Women , 2006, Pediatrics.

[15]  J. Trogdon,et al.  Geospatial Patterns in Human Papillomavirus Vaccination Uptake: Evidence from Uninsured and Publicly Insured Children in North Carolina , 2015, Cancer Epidemiology, Biomarkers & Prevention.

[16]  Behavioral Domains,et al.  Capturing Social and Behavioral Domains and Measures in Electronic Health Records: Phase 2 , 2015 .

[17]  Beth A. Glenn,et al.  Exploring the Role of Neighborhood Socio-Demographic Factors on HPV Vaccine Initiation Among Low-Income, Ethnic Minority Girls , 2013, Journal of Immigrant and Minority Health.

[18]  K. Koch Introduction to Bayesian Statistics , 2007 .

[19]  Eduardo Lazcano-Ponce,et al.  Safety and Persistent Immunogenicity of a Quadrivalent Human Papillomavirus Types 6, 11, 16, 18 L1 Virus-Like Particle Vaccine in Preadolescents and Adolescents: A Randomized Controlled Trial , 2007, The Pediatric infectious disease journal.

[20]  N. Krieger Historical roots of social epidemiology: socioeconomic gradients in health and contextual analysis. , 2001, International journal of epidemiology.

[21]  Kevin A. Henry,et al.  Geographic Factors and Human Papillomavirus (HPV) Vaccination Initiation among Adolescent Girls in the United States , 2016, Cancer Epidemiology, Biomarkers & Prevention.

[22]  L. Melton,et al.  History of the Rochester Epidemiology Project. , 1996, Mayo Clinic proceedings.

[23]  J. Ebbert,et al.  Younger age at initiation of the human papillomavirus (HPV) vaccination series is associated with higher rates of on-time completion. , 2016, Preventive medicine.

[24]  J. Macneil,et al.  National, Regional, State, and Selected Local Area Vaccination Coverage Among Adolescents Aged 13-17 Years - United States, 2015. , 2016, MMWR. Morbidity and mortality weekly report.

[25]  FDA licensure of bivalent human papillomavirus vaccine (HPV2, Cervarix) for use in females and updated HPV vaccination recommendations from the Advisory Committee on Immunization Practices (ACIP). , 2010, MMWR. Morbidity and mortality weekly report.

[26]  L. Kurland,et al.  The patient record in epidemiology. , 1981, Scientific American.

[27]  S. Stokley,et al.  Barriers to Human Papillomavirus Vaccination Among US Adolescents: A Systematic Review of the Literature , 2015 .

[28]  Jarvis T. Chen,et al.  Geocoding and monitoring of US socioeconomic inequalities in mortality and cancer incidence: does the choice of area-based measure and geographic level matter?: the Public Health Disparities Geocoding Project. , 2002, American journal of epidemiology.

[29]  M. Schootman,et al.  A Comparison of Parent- and Provider-Reported Human Papillomavirus Vaccination of Adolescents. , 2017, American journal of preventive medicine.

[30]  A. D. Diez Roux,et al.  A review of spatial methods in epidemiology, 2000-2010. , 2012, Annual review of public health.

[31]  S. Martino Approximate Bayesian Inference for Latent Gaussian Models , 2007 .

[32]  B. Yawn,et al.  American Journal of Epidemiology Practice of Epidemiology Use of a Medical Records Linkage System to Enumerate a Dynamic Population over Time: the Rochester Epidemiology Project , 2022 .

[33]  J. Pankow,et al.  Geospatial patterns of human papillomavirus vaccine uptake in Minnesota , 2015, BMJ Open.

[34]  S. Stokley,et al.  Human Papillomavirus Vaccination Coverage Among Adolescent Girls, 2007–2012, and Postlicensure Vaccine Safety Monitoring, 2006–2013 — United States , 2013, MMWR. Morbidity and mortality weekly report.