Resolving pertussis immunity and vaccine effectiveness using incidence time series

Pertussis vaccines, epidemiology and evolution

Pertussis, which is caused by Bordetella pertussis, has plagued humans for at least 800 years, is highly infectious and can be fatal in the unvaccinated, especially very young infants. Although the rollout of whole-cell pertussis (wP) vaccines in the 1940s and 1950s was associated with a drastic drop in incidence, concerns regarding the reactogenicity of wP vaccines led to the development of a new generation of safer, acellular (aP) vaccines that have been adopted mainly in high-income countries. Over the past 20 years, some countries that boast high aP coverage have experienced a resurgence in pertussis, which has led to substantial debate over the basic immunology, epidemiology and evolutionary biology of the bacterium. Controversy surrounds the duration of natural immunity and vaccine-derived immunity, the ability of vaccines to prevent transmission and severe disease, and the impact of evolution on evading vaccine immunity. Resolving these issues is made challenging by incomplete detection of pertussis cases, the absence of a serological marker of immunity, modest sequencing of the bacterial genome and heterogeneity in diagnostic methods of surveillance. In this Review, we lay out the complexities of contemporary pertussis and, where possible, propose a parsimonious explanation for apparently incongruous observations.

Immunological heterogeneity informs estimation of the durability of vaccine protection

Deciphering the properties of vaccines against an emerging pathogen is essential for optimizing immunization strategies. Early after vaccine roll-out, however, uncertainties about vaccine immunity raise the question of how much time is needed to estimate these properties, particularly the durability of vaccine protection. Here we designed a simulation study, based on a generic transmission model of vaccination, to simulate the impact of a breadth of vaccines with different mean (range: 10 months-2 years) and variability (coefficient of variation range: 50-100%) of the duration of protection. Focusing on the dynamics of SARS-CoV-2 in the year after start of mass immunization in Germany as a case study, we then assessed how confidently the duration of protection could be estimated under a range of epidemiological scenarios. We found that lower mean and higher heterogeneity facilitated estimation of the duration of vaccine protection. Across the vaccines tested, rapid waning and high heterogeneity permitted complete identification of the duration of protection; by contrast, slow waning and low heterogeneity allowed only estimation of the fraction of vaccinees with rapid loss of immunity. These findings suggest that limited epidemiological data can inform the duration of vaccine immunity. More generally, they highlight the need to carefully consider immunological heterogeneity when designing transmission models to evaluate vaccines.

Dissecting recurrent waves of pertussis across the boroughs of London

Pertussis has resurfaced in the UK, with incidence levels not seen since the 1980s. While the fundamental causes of this resurgence remain the subject of much conjecture, the study of historical patterns of pathogen diffusion can be illuminating. Here, we examined time series of pertussis incidence in the boroughs of Greater London from 1982 to 2013 to document the spatial epidemiology of this bacterial infection and to identify the potential drivers of its percolation. The incidence of pertussis over this period is characterized by 3 distinct stages: a period exhibiting declining trends with 4-year inter-epidemic cycles from 1982 to 1994, followed by a deep trough until 2006 and the subsequent resurgence. We observed systematic temporal trends in the age distribution of cases and the fade-out profile of pertussis coincident with increasing national vaccine coverage from 1982 to 1990. To quantify the hierarchy of epidemic phases across the boroughs of London, we used the Hilbert transform. We report a consistent pattern of spatial organization from 1982 to the early 1990s, with some boroughs consistently leading epidemic waves and others routinely lagging. To determine the potential drivers of these geographic patterns, a comprehensive parallel database of borough-specific features was compiled, comprising of demographic, movement and socio-economic factors that were used in statistical analyses to predict epidemic phase relationships among boroughs. Specifically, we used a combination of a feed-forward neural network (FFNN), and SHapley Additive exPlanations (SHAP) values to quantify the contribution of each covariate to model predictions. Our analyses identified a number of predictors of a borough's historical epidemic phase, specifically the age composition of households, the number of agricultural and skilled manual workers, latitude, the population of public transport commuters and high-occupancy households. Univariate regression analysis of the 2012 epidemic identified the ratio of cumulative unvaccinated children to the total population and population of Pakistan-born population to have moderate positive and negative association, respectively, with the timing of epidemic. In addition to providing a comprehensive overview of contemporary pertussis transmission in a large metropolitan population, this study has identified the characteristics that determine the spatial spread of this bacterium across the boroughs of London.

Durability of protection after 5 doses of acellular pertussis vaccine among 5-9 year old children in King County, Washington

PURPOSE

Waning of immunity after vaccination with the acellular Pertussis (aP) vaccine has been proposed as one of the main reasons for pertussis resurgence in the US. In this study, we estimated time-varying vaccine effectiveness after 5 doses of aP vaccine.

METHODS

We conducted a retrospective cohort study among children 5-9 years old (born between 2008 and 2012) living in King County, Washington, USA, who participated in the Washington State Immunization Information System. We estimated time-varying vaccine effectiveness after 5 doses of aP using smoothed scaled Schoenfeld residuals obtained from fitting Cox proportional hazards models to the data as well as piecewise constant Poisson regression.

RESULTS

There were 55 pertussis cases in this cohort, of whom 22 (40%) were fully-vaccinated and 33 (60%) were under-vaccinated. Vaccine effectiveness (VE) remained high for up to 42 months after the fifth dose (VE(t) = 89%; 95% CI: 64%, 97%) as estimated using survival analysis methods and up to 4 years (VE(t) = 93%; 95% CI: 67%, 98%) as estimated using Poisson regression.

CONCLUSION

We did not find evidence for waning of vaccine effectiveness for up to four years after 5 doses of aP among 5 -9 years old children in King County, WA.

Pertussis seroprevalence among adults of reproductive age (20-39 years) in fourteen European countries

The reported incidence of pertussis in European countries varies considerably. We aimed to study specific Bordetella pertussis seroprevalence in Europe by measuring serum IgG antibody levels to pertussis toxin (anti-PT IgG). Fourteen national laboratories participated in this study including Belgium, Denmark, Finland, Greece, Hungary, Italy, Lithuania, Malta, Norway, Poland, Portugal, Romania, Spain, and Sweden. Each country collected approximately 250 samples (N = 7903) from the age groups 20-29 years (N = 3976) and 30-39 years (N = 3927) during 2010-2013. Samples were anonymous residual sera from diagnostic laboratories and were analyzed at the national laboratories by a Swedish reference method, a commercial ELISA kit, or were sent to Sweden for analysis. The median anti-PT IgG concentrations ranged from 4 to 13.6 IU/mL. The proportion of samples with anti-PT IgG ≥100 IU/mL, indicating a recent infection ranged from 0.2% (Hungary) to 5.7% (Portugal). The highest proportion of sera with anti-PT IgG levels between 50 and <100 IU/mL, indicating an infection within the last few years, was found in Portugal (12.3%) and Italy (13.9%). This study shows that the circulation of B. pertussis is quite extensive in adults, aged 20-39 years, despite well-established vaccination programs in Europe.

Long-term changes in endemic threshold populations for pertussis in England and Wales: A spatiotemporal analysis of Lancashire and South Wales, 1940-69

Metapopulation dynamics play a critical role in driving endemic persistence and transmission of childhood infections. The endemic threshold concept, also referred to as critical community size (CCS), is a key example and is defined as the minimum population size required to sustain a continuous chain of infection transmission. The concept is fundamental to the implementation of effective vaccine-based disease control programmes. Vaccination serves to increase endemic threshold population size, promoting disease fadeout and eventual elimination of infection. To date, empirical investigations of the relationship between vaccination and endemic threshold population size have tended to focus on isolated populations in island communities. Very few studies have examined endemic threshold dynamics in 'mainland' regional populations with complex hierarchical spatial structures and varying levels of connectivity between subpopulations. The present paper provides the first spatially explicit analysis of the temporal changes in endemic threshold populations for one vaccine-preventable childhood infection (pertussis) in two dynamic regions of England and Wales: Lancashire and South Wales. Drawing upon weekly disease records of the Registrar-General of England and Wales over a 30-year period (January 1940-December 1969) regression techniques were used to estimate the endemic threshold size for pertussis in the two study regions. Survival analyses were performed to compare disease fadeout duration and probability for both regions in the pre-vaccine and vaccine eras, respectively. Our findings reveal the introduction of mass vaccination led to a considerable increase in threshold size for both Lancashire (~387,333) and South Wales (~1,460,667). Significant growth in fadeout duration was observed in the vaccine era for pertussis non-hotspots in both regions, consistent with geographical synchronisation of epidemic activity. Regional differences in endemic threshold populations reflect significant regional variations in spatial connectivity, population dispersion and level of geographical isolation.

Commentary: resolving pertussis resurgence and vaccine immunity using mathematical transmission models

The epidemiology of pertussis-a vaccine-preventable respiratory infection typically caused by the bacterium Bordetella pertussis-remains puzzling. Indeed, the disease seems nowhere close to eradication and has even re-emerged in certain countries-such as the US-that have maintained high vaccination coverage. Because the dynamics of pertussis are shaped by past vaccination and natural infection rates, with the relevant timescale spanning decades, the interpretation of such unexpected trends is not straightforward. In this commentary, we propose that mathematical transmission models play an essential role in helping to interpret the data and in closing knowledge gaps in pertussis epidemiology. We submit that recent advances in statistical inference methods now allow us to estimate key parameters, such as the nature and duration of vaccinal immunity, which have to date been difficult to quantify. We illustrate these points with the results of a recent study based on data from Massachusetts (Domenech de Cellès, Magpantay, King, and Rohani, Sci. Transl. Med. 2018;10: eaaj1748. doi:10.1126/scitranslmed.aaj1748), in which we used such methods to elucidate the mechanisms underlying the ongoing resurgence of pertussis. In addition, we list a number of safety checks that can be used to critically assess mathematical models. Finally, we discuss the remaining uncertainties surrounding pertussis vaccines, in particular the acellular vaccines used for teenage booster immunizations.

A second-order iterated smoothing algorithm

Simulation-based inference for partially observed stochastic dynamic models is currently receiving much attention due to the fact that direct computation of the likelihood is not possible in many practical situations. Iterated filtering methodologies enable maximization of the likelihood function using simulation-based sequential Monte Carlo filters. Doucet et al. (2013) developed an approximation for the first and second derivatives of the log likelihood via simulation-based sequential Monte Carlo smoothing and proved that the approximation has some attractive theoretical properties. We investigated an iterated smoothing algorithm carrying out likelihood maximization using these derivative approximations. Further, we developed a new iterated smoothing algorithm, using a modification of these derivative estimates, for which we establish both theoretical results and effective practical performance. On benchmark computational challenges, this method beat the first-order iterated filtering algorithm. The method's performance was comparable to a recently developed iterated filtering algorithm based on an iterated Bayes map. Our iterated smoothing algorithm and its theoretical justification provide new directions for future developments in simulation-based inference for latent variable models such as partially observed Markov process models.

Effects of reactive social distancing on the 1918 influenza pandemic

The 1918 influenza pandemic was characterized by multiple epidemic waves. We investigated reactive social distancing, a form of behavioral response where individuals avoid potentially infectious contacts in response to available information on an ongoing epidemic or pandemic. We modelled its effects on the three influenza waves in the United Kingdom. In previous studies, human behavioral response was modelled by a Power function of the proportion of recent influenza mortality in a population, and by a Hill function, which is a function of the number of recent influenza mortality. Using a simple epidemic model with a Power function and one common set of parameters, we provided a good model fit for the observed multiple epidemic waves in London boroughs, Birmingham and Liverpool. We further applied the model parameters from these three cities to all 334 administrative units in England and Wales and including the population sizes of individual administrative units. We computed the Pearson's correlation between the observed and simulated for each administrative unit. We found a median correlation of 0.636, indicating that our model predictions are performing reasonably well. Our modelling approach is an improvement from previous studies where separate models are fitted to each city. With the reduced number of model parameters used, we achieved computational efficiency gain without over-fitting the model. We also showed the importance of reactive behavioral distancing as a potential non-pharmaceutical intervention during an influenza pandemic. Our work has both scientific and public health significance.

The pertussis enigma: reconciling epidemiology, immunology and evolution

Pertussis, a highly contagious respiratory infection, remains a public health priority despite the availability of vaccines for 70 years. Still a leading cause of mortality in developing countries, pertussis has re-emerged in several developed countries with high vaccination coverage. Resurgence of pertussis in these countries has routinely been attributed to increased awareness of the disease, imperfect vaccinal protection or high infection rates in adults. In this review, we first present 1980-2012 incidence data from 63 countries and show that pertussis resurgence is not universal. We further argue that the large geographical variation in trends probably precludes a simple explanation, such as the transition from whole-cell to acellular pertussis vaccines. Reviewing available evidence, we then propose that prevailing views on pertussis epidemiology are inconsistent with both historical and contemporary data. Indeed, we summarize epidemiological evidence showing that natural infection and vaccination both appear to provide long-term protection against transmission and disease, so that previously infected or vaccinated adults contribute little to overall transmission at a population level. Finally, we identify several promising avenues that may lead to a consistent explanation of global pertussis epidemiology and to more effective control strategies.

Pertussis immunity and epidemiology: mode and duration of vaccine-induced immunity

The resurgence of pertussis in some countries that maintain high vaccination coverage has drawn attention to gaps in our understanding of the epidemiological effects of pertussis vaccines. In particular, major questions surround the nature, degree and durability of vaccine protection. To address these questions, we used mechanistic transmission models to examine regional time series incidence data from Italy in the period immediately following the introduction of acellular pertussis (aP) vaccine. Our results concur with recent animal-challenge experiments wherein infections in aP-vaccinated individuals proved as transmissible as those in naive individuals but much less symptomatic. On the other hand, the data provide evidence for vaccine-driven reduction in susceptibility, which we quantify via a synthetic measure of vaccine impact. As to the precise nature of vaccine failure, the data do not allow us to distinguish between leakiness and waning of vaccine immunity, or some combination of these. Across the range of well-supported models, the nature and duration of vaccine protection, the age profile of incidence and the range of projected epidemiological futures differ substantially, underscoring the importance of the remaining unknowns. We identify key data gaps: sources of data that can supply the information needed to eliminate these remaining uncertainties.

An evidence synthesis approach to estimating the incidence of symptomatic pertussis infection in the Netherlands, 2005-2011

BACKGROUND

Despite high vaccination coverage, infection with Bordetella pertussis is a current public health concern in the Netherlands and other European Union member states. Because surveillance data are subject to extensive under-ascertainment and under-reporting, incidence is difficult to determine. Our objective was to estimate the age-group specific incidence of symptomatic pertussis infection in the Netherlands over the period 2005-2011, using multi-parameter evidence synthesis.

METHODS

Age-specific seroconversion probabilities were estimated for 2007 using Netherlands population data stratified by age-group and cross-sectional population-wide serosurvey (PIENTER-2) data, with a sero-diagnostic cut-off of 125 EU/ml as a proxy for recent infection. Symptomatic probabilities were derived from a study of household contacts and from PIENTER-2. The annual number of symptomatic infected (SI) persons was estimated using evidence synthesis methods in a Bayesian framework, by combining the estimated incidence of infection with notification data and symptomatic probabilities.

RESULTS

An incidence rate of 128 SI cases per 10,000 population (95 % credible interval [CrI]: 110-150) was estimated for 2005, which decreased to 107 per 10,000 (95 % CrI: 91-126) for 2011. The degree of underestimation in statutory notified cases was age-dependent, ranging from 10-fold (10-19 years) to 69-fold (60+ years). The largest annual decreases in SI incidence rate over the study period were in the 1-4 and 5-9 years age-groups (24.3 %, 15.9 % per year, respectively).

CONCLUSIONS

By synthesising all available data, the incidence of symptomatic pertussis and the extent to which SI is underrepresented by notification data can be estimated. Such estimates are essential for disease burden computation and for informing public health priority-setting.

Demographic buffering: titrating the effects of birth rate and imperfect immunity on epidemic dynamics

Host demography can alter the dynamics of infectious disease. In the case of perfectly immunizing infections, observations of strong sensitivity to demographic variation have been mechanistically explained through analysis of the susceptible–infected–recovered (SIR) model that assumes lifelong immunity following recovery from infection. When imperfect immunity is incorporated into this framework via the susceptible–infected–recovered–susceptible (SIRS) model, with individuals regaining full susceptibility following recovery, we show that rapid loss of immunity is predicted to buffer populations against the effects of demographic change. However, this buffering is contrary to the dependence on demography recently observed for partially immunizing infections such as rotavirus and respiratory syncytial virus. We show that this discrepancy arises from a key simplification embedded in the SIR(S) framework, namely that the potential for differential immune responses to repeat exposures is ignored. We explore the minimum additional immunological information that must be included to reflect the range of observed dependencies on demography. We show that including partial protection and lower transmission following primary infection is sufficient to capture more realistic reduced levels of buffering, in addition to changes in epidemic timing, across a range of partially and fully immunizing infections. Furthermore, our results identify key variables in this relationship, including R₀.

Combating pertussis resurgence: One booster vaccination schedule does not fit all

Pertussis has reemerged as a major public health concern in many countries where it was once considered well controlled. Although the mechanisms responsible for continued pertussis circulation and resurgence remain elusive and contentious, many countries have nevertheless recommended booster vaccinations, the timing and number of which vary widely. Here, using a stochastic, age-stratified transmission model, we searched for cost-effective booster vaccination strategies using a genetic algorithm. We did so assuming four hypothesized mechanisms underpinning contemporary pertussis epidemiology: (I) insufficient coverage, (II) frequent primary vaccine failure, (III) waning of vaccine-derived protection, and (IV) vaccine "leakiness." For scenarios I-IV, successful booster strategies were identified and varied considerably by mechanism. Especially notable is the inability of booster schedules to alleviate resurgence when vaccines are leaky. Critically, our findings argue that the ultimate effectiveness of vaccine booster schedules will likely depend on correctly pinpointing the causes of resurgence, with misdiagnosis of the problem epidemiologically ineffective and economically costly.