Association between pertussis vaccination coverage and other sociodemographic factors and pertussis incidence using surveillance data

Routine vaccination with pertussis vaccines has been successful in driving down pertussis mortality and morbidity globally. Despite high vaccination coverage, countries such as Australia, USA, and UK have experienced increase in pertussis activity over the last few decades. This may be due to local pockets of low vaccination coverage that result in persistence of pertussis in the population and occasionally lead to large outbreaks. The objective of this study was to characterize the association between pertussis vaccination coverage and sociodemographic factors and pertussis incidence at the school district level in King County, Washington, USA. We used monthly pertussis incidence data for all ages reported to the Public Health Seattle and King County between January 1, 2010 and December 31, 2017 to obtain school district level pertussis incidence. We obtained immunization data from the Washington State Immunization Information System to estimate school-district level vaccination coverage as proportion of 19-35 month old children fully vaccinated with ≥4 doses of the Diphtheria-Tetanus-acellular-Pertussis (DTaP) vaccine in a school district. We used two methods to quantify the effects of vaccination coverage on pertussis incidence: an ecological vaccine model and an endemic-epidemic model. Even though the effect of vaccination is modeled differently in the two approaches, both models can be used to estimate the association between vaccination coverage and pertussis incidence. Using the ecological vaccine model, we estimated the vaccine effectiveness of 4 doses of Diphtheria-Tetanus-acellular-Pertussis vaccine to be 83% (95% credible interval: 63%, 95%). In the endemic-epidemic model, under-vaccination was statistically significantly associated with epidemic risk of pertussis (adjusted Relative Risk, aRR: 2.76; 95% confidence interval: 1.44, 16.6). Household size and median income were statistically significantly associated with endemic pertussis risk. The endemic-epidemic model suffers from ecological bias, whereas the ecological vaccine model provides less biased and more interpretable estimates of epidemiological parameters, such as DTaP vaccine effectiveness, for each school district.