Risk factors and short-term projections for serotype-1 poliomyelitis incidence in Pakistan: A spatiotemporal analysis

Natalie A. Molodecky, Isobel M. Blake, Kathleen M. O’Reilly, Mufti Zubair Wadood, Rana M. Safdar, Amy Wesolowski, Caroline O. Buckee, Ananda S. Bandyopadhyay, Hiromasa Okayasu, Nicholas C. Grassly

Research output: Contribution to journalArticle

Abstract

Background: Pakistan currently provides a substantial challenge to global polio eradication, having contributed to 73% of reported poliomyelitis in 2015 and 54% in 2016. A better understanding of the risk factors and movement patterns that contribute to poliovirus transmission across Pakistan would support evidence-based planning for mass vaccination campaigns. Methods and findings: We fit mixed-effects logistic regression models to routine surveillance data recording the presence of poliomyelitis associated with wild-type 1 poliovirus in districts of Pakistan over 6-month intervals between 2010 to 2016. To accurately capture the force of infection (FOI) between districts, we compared 6 models of population movement (adjacency, gravity, radiation, radiation based on population density, radiation based on travel times, and mobile-phone based). We used the best-fitting model (based on the Akaike Information Criterion [AIC]) to produce 6-month forecasts of poliomyelitis incidence. The odds of observing poliomyelitis decreased with improved routine or supplementary (campaign) immunisation coverage (multivariable odds ratio [OR] = 0.75, 95% confidence interval [CI] 0.67–0.84; and OR = 0.75, 95% CI 0.66–0.85, respectively, for each 10% increase in coverage) and increased with a higher rate of reporting non-polio acute flaccid paralysis (AFP) (OR = 1.13, 95% CI 1.02–1.26 for a 1-unit increase in non-polio AFP per 100,000 persons aged <15 years). Estimated movement of poliovirus-infected individuals was associated with the incidence of poliomyelitis, with the radiation model of movement providing the best fit to the data. Six-month forecasts of poliomyelitis incidence by district for 2013–2016 showed good predictive ability (area under the curve range: 0.76–0.98). However, although the best-fitting movement model (radiation) was a significant determinant of poliomyelitis incidence, it did not improve the predictive ability of the multivariable model. Overall, in Pakistan the risk of polio cases was predicted to reduce between July–December 2016 and January–June 2017. The accuracy of the model may be limited by the small number of AFP cases in some districts. Conclusions: Spatiotemporal variation in immunization performance and population movement patterns are important determinants of historical poliomyelitis incidence in Pakistan; however, movement dynamics were less influential in predicting future cases, at a time when the polio map is shrinking. Results from the regression models we present are being used to help plan vaccination campaigns and transit vaccination strategies in Pakistan.

Original languageEnglish (US)
Article numbere1002323
JournalPLoS Medicine
Volume14
Issue number6
DOIs
StatePublished - Jun 1 2017
Externally publishedYes

Fingerprint

Spatio-Temporal Analysis
Pakistan
Poliomyelitis
Incidence
Poliovirus
Radiation
Paralysis
Immunization Programs
Odds Ratio
Confidence Intervals
Immunization
Serogroup
Logistic Models
Mass Vaccination
Cell Phones
Gravitation
Population Density
Population
Area Under Curve
Vaccination

ASJC Scopus subject areas

  • Medicine(all)

Cite this

Molodecky, N. A., Blake, I. M., O’Reilly, K. M., Wadood, M. Z., Safdar, R. M., Wesolowski, A., ... Grassly, N. C. (2017). Risk factors and short-term projections for serotype-1 poliomyelitis incidence in Pakistan: A spatiotemporal analysis. PLoS Medicine, 14(6), [e1002323]. https://doi.org/10.1371/journal.pmed.1002323

Risk factors and short-term projections for serotype-1 poliomyelitis incidence in Pakistan : A spatiotemporal analysis. / Molodecky, Natalie A.; Blake, Isobel M.; O’Reilly, Kathleen M.; Wadood, Mufti Zubair; Safdar, Rana M.; Wesolowski, Amy; Buckee, Caroline O.; Bandyopadhyay, Ananda S.; Okayasu, Hiromasa; Grassly, Nicholas C.

In: PLoS Medicine, Vol. 14, No. 6, e1002323, 01.06.2017.

Research output: Contribution to journalArticle

Molodecky, NA, Blake, IM, O’Reilly, KM, Wadood, MZ, Safdar, RM, Wesolowski, A, Buckee, CO, Bandyopadhyay, AS, Okayasu, H & Grassly, NC 2017, 'Risk factors and short-term projections for serotype-1 poliomyelitis incidence in Pakistan: A spatiotemporal analysis', PLoS Medicine, vol. 14, no. 6, e1002323. https://doi.org/10.1371/journal.pmed.1002323
Molodecky, Natalie A. ; Blake, Isobel M. ; O’Reilly, Kathleen M. ; Wadood, Mufti Zubair ; Safdar, Rana M. ; Wesolowski, Amy ; Buckee, Caroline O. ; Bandyopadhyay, Ananda S. ; Okayasu, Hiromasa ; Grassly, Nicholas C. / Risk factors and short-term projections for serotype-1 poliomyelitis incidence in Pakistan : A spatiotemporal analysis. In: PLoS Medicine. 2017 ; Vol. 14, No. 6.
@article{6af4deb88f67491ea822cd895ec661c9,
title = "Risk factors and short-term projections for serotype-1 poliomyelitis incidence in Pakistan: A spatiotemporal analysis",
abstract = "Background: Pakistan currently provides a substantial challenge to global polio eradication, having contributed to 73{\%} of reported poliomyelitis in 2015 and 54{\%} in 2016. A better understanding of the risk factors and movement patterns that contribute to poliovirus transmission across Pakistan would support evidence-based planning for mass vaccination campaigns. Methods and findings: We fit mixed-effects logistic regression models to routine surveillance data recording the presence of poliomyelitis associated with wild-type 1 poliovirus in districts of Pakistan over 6-month intervals between 2010 to 2016. To accurately capture the force of infection (FOI) between districts, we compared 6 models of population movement (adjacency, gravity, radiation, radiation based on population density, radiation based on travel times, and mobile-phone based). We used the best-fitting model (based on the Akaike Information Criterion [AIC]) to produce 6-month forecasts of poliomyelitis incidence. The odds of observing poliomyelitis decreased with improved routine or supplementary (campaign) immunisation coverage (multivariable odds ratio [OR] = 0.75, 95{\%} confidence interval [CI] 0.67–0.84; and OR = 0.75, 95{\%} CI 0.66–0.85, respectively, for each 10{\%} increase in coverage) and increased with a higher rate of reporting non-polio acute flaccid paralysis (AFP) (OR = 1.13, 95{\%} CI 1.02–1.26 for a 1-unit increase in non-polio AFP per 100,000 persons aged <15 years). Estimated movement of poliovirus-infected individuals was associated with the incidence of poliomyelitis, with the radiation model of movement providing the best fit to the data. Six-month forecasts of poliomyelitis incidence by district for 2013–2016 showed good predictive ability (area under the curve range: 0.76–0.98). However, although the best-fitting movement model (radiation) was a significant determinant of poliomyelitis incidence, it did not improve the predictive ability of the multivariable model. Overall, in Pakistan the risk of polio cases was predicted to reduce between July–December 2016 and January–June 2017. The accuracy of the model may be limited by the small number of AFP cases in some districts. Conclusions: Spatiotemporal variation in immunization performance and population movement patterns are important determinants of historical poliomyelitis incidence in Pakistan; however, movement dynamics were less influential in predicting future cases, at a time when the polio map is shrinking. Results from the regression models we present are being used to help plan vaccination campaigns and transit vaccination strategies in Pakistan.",
author = "Molodecky, {Natalie A.} and Blake, {Isobel M.} and O’Reilly, {Kathleen M.} and Wadood, {Mufti Zubair} and Safdar, {Rana M.} and Amy Wesolowski and Buckee, {Caroline O.} and Bandyopadhyay, {Ananda S.} and Hiromasa Okayasu and Grassly, {Nicholas C.}",
year = "2017",
month = "6",
day = "1",
doi = "10.1371/journal.pmed.1002323",
language = "English (US)",
volume = "14",
journal = "PLoS Medicine",
issn = "1549-1277",
publisher = "Nature Publishing Group",
number = "6",

}

TY - JOUR

T1 - Risk factors and short-term projections for serotype-1 poliomyelitis incidence in Pakistan

T2 - A spatiotemporal analysis

AU - Molodecky, Natalie A.

AU - Blake, Isobel M.

AU - O’Reilly, Kathleen M.

AU - Wadood, Mufti Zubair

AU - Safdar, Rana M.

AU - Wesolowski, Amy

AU - Buckee, Caroline O.

AU - Bandyopadhyay, Ananda S.

AU - Okayasu, Hiromasa

AU - Grassly, Nicholas C.

PY - 2017/6/1

Y1 - 2017/6/1

N2 - Background: Pakistan currently provides a substantial challenge to global polio eradication, having contributed to 73% of reported poliomyelitis in 2015 and 54% in 2016. A better understanding of the risk factors and movement patterns that contribute to poliovirus transmission across Pakistan would support evidence-based planning for mass vaccination campaigns. Methods and findings: We fit mixed-effects logistic regression models to routine surveillance data recording the presence of poliomyelitis associated with wild-type 1 poliovirus in districts of Pakistan over 6-month intervals between 2010 to 2016. To accurately capture the force of infection (FOI) between districts, we compared 6 models of population movement (adjacency, gravity, radiation, radiation based on population density, radiation based on travel times, and mobile-phone based). We used the best-fitting model (based on the Akaike Information Criterion [AIC]) to produce 6-month forecasts of poliomyelitis incidence. The odds of observing poliomyelitis decreased with improved routine or supplementary (campaign) immunisation coverage (multivariable odds ratio [OR] = 0.75, 95% confidence interval [CI] 0.67–0.84; and OR = 0.75, 95% CI 0.66–0.85, respectively, for each 10% increase in coverage) and increased with a higher rate of reporting non-polio acute flaccid paralysis (AFP) (OR = 1.13, 95% CI 1.02–1.26 for a 1-unit increase in non-polio AFP per 100,000 persons aged <15 years). Estimated movement of poliovirus-infected individuals was associated with the incidence of poliomyelitis, with the radiation model of movement providing the best fit to the data. Six-month forecasts of poliomyelitis incidence by district for 2013–2016 showed good predictive ability (area under the curve range: 0.76–0.98). However, although the best-fitting movement model (radiation) was a significant determinant of poliomyelitis incidence, it did not improve the predictive ability of the multivariable model. Overall, in Pakistan the risk of polio cases was predicted to reduce between July–December 2016 and January–June 2017. The accuracy of the model may be limited by the small number of AFP cases in some districts. Conclusions: Spatiotemporal variation in immunization performance and population movement patterns are important determinants of historical poliomyelitis incidence in Pakistan; however, movement dynamics were less influential in predicting future cases, at a time when the polio map is shrinking. Results from the regression models we present are being used to help plan vaccination campaigns and transit vaccination strategies in Pakistan.

AB - Background: Pakistan currently provides a substantial challenge to global polio eradication, having contributed to 73% of reported poliomyelitis in 2015 and 54% in 2016. A better understanding of the risk factors and movement patterns that contribute to poliovirus transmission across Pakistan would support evidence-based planning for mass vaccination campaigns. Methods and findings: We fit mixed-effects logistic regression models to routine surveillance data recording the presence of poliomyelitis associated with wild-type 1 poliovirus in districts of Pakistan over 6-month intervals between 2010 to 2016. To accurately capture the force of infection (FOI) between districts, we compared 6 models of population movement (adjacency, gravity, radiation, radiation based on population density, radiation based on travel times, and mobile-phone based). We used the best-fitting model (based on the Akaike Information Criterion [AIC]) to produce 6-month forecasts of poliomyelitis incidence. The odds of observing poliomyelitis decreased with improved routine or supplementary (campaign) immunisation coverage (multivariable odds ratio [OR] = 0.75, 95% confidence interval [CI] 0.67–0.84; and OR = 0.75, 95% CI 0.66–0.85, respectively, for each 10% increase in coverage) and increased with a higher rate of reporting non-polio acute flaccid paralysis (AFP) (OR = 1.13, 95% CI 1.02–1.26 for a 1-unit increase in non-polio AFP per 100,000 persons aged <15 years). Estimated movement of poliovirus-infected individuals was associated with the incidence of poliomyelitis, with the radiation model of movement providing the best fit to the data. Six-month forecasts of poliomyelitis incidence by district for 2013–2016 showed good predictive ability (area under the curve range: 0.76–0.98). However, although the best-fitting movement model (radiation) was a significant determinant of poliomyelitis incidence, it did not improve the predictive ability of the multivariable model. Overall, in Pakistan the risk of polio cases was predicted to reduce between July–December 2016 and January–June 2017. The accuracy of the model may be limited by the small number of AFP cases in some districts. Conclusions: Spatiotemporal variation in immunization performance and population movement patterns are important determinants of historical poliomyelitis incidence in Pakistan; however, movement dynamics were less influential in predicting future cases, at a time when the polio map is shrinking. Results from the regression models we present are being used to help plan vaccination campaigns and transit vaccination strategies in Pakistan.

UR - http://www.scopus.com/inward/record.url?scp=85021784358&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85021784358&partnerID=8YFLogxK

U2 - 10.1371/journal.pmed.1002323

DO - 10.1371/journal.pmed.1002323

M3 - Article

C2 - 28604777

AN - SCOPUS:85021784358

VL - 14

JO - PLoS Medicine

JF - PLoS Medicine

SN - 1549-1277

IS - 6

M1 - e1002323

ER -