Routine immunization in Pakistan: comparison of multiple data sources and identification of factors associated with vaccination.

Background: Within Pakistan, estimates of vaccination coverage with the pentavalent vaccine, oral polio vaccine (OPV) and measles vaccine (MV) in 2011 were reported to be 74%, 75% and 53%, respectively. These national estimates may mask regional variation. The reasons for this variation have not been explored. Methods: Data from the Multiple Indicator Cluster Surveys (MICS) for Balochistan and Punjab (2010-2011) are analysed to examine factors associated with receiving three or more doses of the pentavalent vaccine and one or more MVs using regression modelling. Pentavalent and OPV estimates from the MICS were compared to vaccine dose histories from surveillance for acute flaccid paralysis (AFP; poliomyelitis) to ascertain agreement. Results: Adjusted coverage of children 12-23 months of age were estimated to be 16.0%, 75.5% and 34.2% in Balochistan and 58.0%, 87.7% and 72.6% in Punjab for the pentavalent vaccine, OPV and MV, respectively. Maternal education, healthcare utilization and wealth were associated with receiving the pentavalent vaccine and the MV. There was a strong correlation of district estimates of vaccination coverage between AFP and MICS data, but AFP estimates of pentavalent coverage in Punjab were biased toward higher values. Conclusions: National estimates mask variation and estimates from individual surveys should be considered alongside other estimates. The development of strategies targeted towards poorly educated parents within low-wealth quintiles that may not typically access healthcare could improve vaccination rates

Tooling-up for infectious disease transmission modelling.

In this introduction to the Special Issue on methods for modelling of infectious disease epidemiology we provide a commentary and overview of the field. We suggest that the field has been through three revolutions that have focussed on specific methodological developments; disease dynamics and heterogeneity, advanced computing and inference, and complexity and application to the real-world. Infectious disease dynamics and heterogeneity dominated until the 1980s where the use of analytical models illustrated fundamental concepts such as herd immunity. The second revolution embraced the integration of data with models and the increased use of computing. From the turn of the century an emergence of novel datasets enabled improved modelling of real-world complexity. The emergence of more complex data that reflect the real-world heterogeneities in transmission resulted in the development of improved inference methods such as particle filtering. Each of these three revolutions have always kept the understanding of infectious disease spread as its motivation but have been developed through the use of new techniques, tools and the availability of data. We conclude by providing a commentary on what the next revoluition in infectious disease modelling may be

The impact of surveillance and other factors on detection of emergent and circulating vaccine derived polioviruses.

Background: Circulating vaccine derived poliovirus (cVDPV) outbreaks remain a threat to polio eradication. To reduce cases of polio from cVDPV of serotype 2, the serotype 2 component of the vaccine has been removed from the global vaccine supply, but outbreaks of cVDPV2 have continued. The objective of this work is to understand the factors associated with later detection in order to improve detection of these unwanted events. Methods: The number of nucleotide differences between each cVDPV outbreak and the oral polio vaccine (OPV) strain was used to approximate the time from emergence to detection. Only independent emergences were included in the analysis. Variables such as serotype, surveillance quality, and World Health Organization (WHO) region were tested in a negative binomial regression model to ascertain whether these variables were associated with higher nucleotide differences upon detection. Results: In total, 74 outbreaks were analysed from 24 countries between 2004-2019. For serotype 1 (n=10), the median time from seeding until outbreak detection was 284 (95% uncertainty interval (UI) 284-2008) days, for serotype 2 (n=59), 276 (95% UI 172-765) days, and for serotype 3 (n=5), 472 (95% UI 392-603) days. Significant improvement in the time to detection was found with increasing surveillance of non-polio acute flaccid paralysis (AFP) and adequate stool collection. Conclusions: cVDPVs remain a risk; all WHO regions have reported at least one VDPV outbreak since the first outbreak in 2000 and outbreak response campaigns using monovalent OPV type 2 risk seeding future outbreaks. Maintaining surveillance for poliomyelitis after local elimination is essential to quickly respond to both emergence of VDPVs and potential importations as low-quality AFP surveillance causes outbreaks to continue undetected. Considerable variation in the time between emergence and detection of VDPVs were apparent, and other than surveillance quality and inclusion of environmental surveillance, the reasons for this remain unclear

The impact of surveillance and other factors on detection of emergent and circulating vaccine derived polioviruses [version 3; peer review: 2 approved, 1 approved with reservations]

Background: Circulating vaccine derived poliovirus (cVDPV) outbreaks remain a threat to polio eradication. To reduce cases of polio from cVDPV of serotype 2, the serotype 2 component of the vaccine has been removed from the global vaccine supply, but outbreaks of cVDPV2 have continued. The objective of this work is to understand the factors associated with later detection in order to improve detection of these unwanted events. Methods: The number of nucleotide differences between each cVDPV outbreak and the oral polio vaccine (OPV) strain was used to approximate the time from emergence to detection. Only independent emergences were included in the analysis. Variables such as serotype, surveillance quality, and World Health Organization (WHO) region were tested in a negative binomial regression model to ascertain whether these variables were associated with higher nucleotide differences upon detection. Results: In total, 74 outbreaks were analysed from 24 countries between 2004-2019. For serotype 1 (n=10), the median time from seeding until outbreak detection was 572 (95% uncertainty interval (UI) 279-2016), for serotype 2 (n=59), 276 (95% UI 172-765) days, and for serotype 3 (n=5), 472 (95% UI 392-603) days. Significant improvement in the time to detection was found with increasing surveillance of non-polio acute flaccid paralysis (AFP) and adequate stool collection. Conclusions: cVDPVs remain a risk; all WHO regions have reported at least one VDPV outbreak since the first outbreak in 2000 and outbreak response campaigns using monovalent OPV type 2 risk seeding future outbreaks. Maintaining surveillance for poliomyelitis after local elimination is essential to quickly respond to both emergence of VDPVs and potential importations as low-quality AFP surveillance causes outbreaks to continue undetected. Considerable variation in the time between emergence and detection of VDPVs were apparent, and other than surveillance quality and inclusion of environmental surveillance, the reasons for this remain unclear

Desirable BUGS in models of infectious diseases.

Bayesian inference using Gibbs sampling (BUGS) is a set of statistical software that uses Markov chain Monte Carlo (MCMC) methods to estimate almost any specified model. Originally developed in the late 1980s, the software is an excellent introduction to applied Bayesian statistics without the need to write a MCMC sampler. The software is typically used for regression-based analyses, but any model that can be specified using graphical nodes are possible. Advanced topics such as missing data, spatial analysis, model comparison and dynamic infectious disease models can be tackled. Three examples are provided; a linear regression model to illustrate parameter estimation, the steps to ensure that the estimates have converged and a comparison of run-times across different computing platforms. The second example describes a model that estimates the probability of being vaccinated from cross-sectional and surveillance data, and illustrates the specification of different models, model comparison and data augmentation. The third example illustrates estimation of parameters within a dynamic Susceptible-Infected-Recovered model. These examples show that BUGS can be used to estimate parameters from models relevant for infectious diseases, and provide an overview of the relative merits of the approach taken

Population sensitivity of acute flaccid paralysis and environmental surveillance for serotype 1 poliovirus in Pakistan: an observational study.

BACKGROUND: To support poliomyelitis eradication in Pakistan, environmental surveillance (ES) of wastewater has been expanded alongside surveillance for acute flaccid paralysis (AFP). ES is a relatively new method of surveillance, and the population sensitivity of detecting poliovirus within endemic settings requires estimation. METHODS: Data for wild serotype 1 poliovirus from AFP and ES from January 2011 to September 2015 from 14 districts in Pakistan were analysed using a multi-state model framework. This framework was used to estimate the sensitivity of poliovirus detection from each surveillance source and parameters such as the duration of infection within a community. RESULTS: The location and timing of poliomyelitis cases showed spatial and temporal variability. The sensitivity of AFP surveillance to detect serotype 1 poliovirus infection in a district and its neighbours per month was on average 30.0% (95% CI 24.8-35.8) and increased with the incidence of poliomyelitis cases. The average population sensitivity of a single environmental sample was 59.4% (95% CI 55.4-63.0), with significant variation in site-specific estimates (median varied from 33.3-79.2%). The combined population sensitivity of environmental and AFP surveillance in a given month was on average 98.1% (95% CI 97.2-98.7), assuming four samples per month for each site. CONCLUSIONS: ES can be a highly sensitive supplement to AFP surveillance in areas with converging sewage systems. As ES for poliovirus is expanded, it will be important to identify factors associated with variation in site sensitivity, leading to improved site selection and surveillance system performance

The cost-effectiveness of controlling dengue in Indonesia using wMel Wolbachia released at scale: a modelling study.

BACKGROUND: Release of virus-blocking Wolbachia-infected mosquitoes is an emerging disease control strategy that aims to control dengue and other arboviral infections. Early entomological data and modelling analyses have suggested promising outcomes, and wMel Wolbachia releases are now ongoing or planned in 12 countries. To help inform government, donor, or philanthropist decisions on scale-up beyond single city releases, we assessed this technology's cost-effectiveness under alternative programmatic options. METHODS: Using costing data from existing Wolbachia releases, previous dynamic model-based estimates of Wolbachia effectiveness, and a spatially explicit model of release and surveillance requirements, we predicted the costs and effectiveness of the ongoing programme in Yogyakarta City and three new hypothetical programmes in Yogyakarta Special Autonomous Region, Jakarta, and Bali. RESULTS: We predicted Wolbachia to be a highly cost-effective intervention when deployed in high-density urban areas with gross cost-effectiveness below $1500 per DALY averted. When offsets from the health system and societal perspective were included, such programmes even became cost saving over 10-year time horizons with favourable benefit-cost ratios of 1.35 to 3.40. Sequencing Wolbachia releases over 10 years could reduce programme costs by approximately 38% compared to simultaneous releases everywhere, but also delays the benefits. Even if unexpected challenges occurred during deployment, such as emergence of resistance in the medium-term or low effective coverage, Wolbachia would remain a cost-saving intervention. CONCLUSIONS: Wolbachia releases in high-density urban areas are expected to be highly cost-effective and could potentially be the first cost-saving intervention for dengue. Sites with strong public health infrastructure, fiscal capacity, and community support should be prioritised