Epidemiology of the silent polio outbreak in Rahat, Israel, based on modeling of environmental surveillance data

Significance The 2013–2014 silent polio epidemic in Israel was a setback to global eradication efforts because Israel had previously been certified as polio-free by the World Health Organization. Fortunately, Israel has a robust environmental surveillance program that detected the epidemic and allowed rapid mobilization of a vaccine campaign before any cases of acute flaccid paralysis. This kind of silent (caseless) epidemic will be increasingly common as we approach global eradication, demonstrating the need for both enhanced environmental surveillance and an accompanying inference framework to translate environmental data into public health metrics. We incorporate environmental data into a population-level disease transmission model, generating insights into the epidemiology of the outbreak. This framework can be used to guide future interventions. Israel experienced an outbreak of wild poliovirus type 1 (WPV1) in 2013–2014, detected through environmental surveillance of the sewage system. No cases of acute flaccid paralysis were reported, and the epidemic subsided after a bivalent oral polio vaccination (bOPV) campaign. As we approach global eradication, polio will increasingly be detected only through environmental surveillance. We developed a framework to convert quantitative polymerase chain reaction (qPCR) cycle threshold data into scaled WPV1 and OPV1 concentrations for inference within a deterministic, compartmental infectious disease transmission model. We used this approach to estimate the epidemic curve and transmission dynamics, as well as assess alternate vaccination scenarios. Our analysis estimates the outbreak peaked in late June, much earlier than previous estimates derived from analysis of stool samples, although the exact epidemic trajectory remains uncertain. We estimate the basic reproduction number was 1.62 (95% CI 1.04–2.02). Model estimates indicate that 59% (95% CI 9–77%) of susceptible individuals (primarily children under 10 years old) were infected with WPV1 over a little more than six months, mostly before the vaccination campaign onset, and that the vaccination campaign averted 10% (95% CI 1–24%) of WPV1 infections. As we approach global polio eradication, environmental monitoring with qPCR can be used as a highly sensitive method to enhance disease surveillance. Our analytic approach brings public health relevance to environmental data that, if systematically collected, can guide eradication efforts.

[1]  Dani Gamerman,et al.  Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference , 1997 .

[2]  M. Pallansch,et al.  Modeling options to manage type 1 wild poliovirus imported into Israel in 2013. , 2015, The Journal of infectious diseases.

[3]  M. Koopmans,et al.  Metagenomic Sequencing for Surveillance of Food- and Waterborne Viral Diseases , 2017, Front. Microbiol..

[4]  R. Sutter,et al.  Vaccine-associated paralytic poliomyelitis: a review of the epidemiology and estimation of the global burden. , 2014, The Journal of infectious diseases.

[5]  Michael Famulare,et al.  Assessing the stability of polio eradication after the withdrawal of oral polio vaccine , 2017, bioRxiv.

[6]  Joel C. Miller,et al.  A Note on the Derivation of Epidemic Final Sizes , 2012, Bulletin of mathematical biology.

[7]  K. Thompson,et al.  Managing population immunity to reduce or eliminate the risks of circulation following the importation of polioviruses. , 2015, Vaccine.

[8]  Shihab U. Sobuz,et al.  Development and assessment of molecular diagnostic tests for 15 enteropathogens causing childhood diarrhoea: a multicentre study. , 2014, The Lancet. Infectious diseases.

[9]  M. Pallansch,et al.  Characterizing Poliovirus Transmission and Evolution: Insights from Modeling Experiences with Wild and Vaccine‐Related Polioviruses , 2013, Risk analysis : an official publication of the Society for Risk Analysis.

[10]  Y. Manor,et al.  Intensified environmental surveillance supporting the response to wild poliovirus type 1 silent circulation in Israel, 2013. , 2014, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[11]  Dani Gamerman,et al.  Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference, Second Edition , 2006 .

[12]  Sheng Li,et al.  Dynamics and control of infections transmitted from person to person through the environment. , 2009, American journal of epidemiology.

[13]  Young Ran Kim,et al.  Characterization of Prohibitin 1 as a Host Partner of Vibrio vulnificus RtxA1 Toxin. , 2016, The Journal of infectious diseases.

[14]  Y. Manor,et al.  Advanced Environmental Surveillance and Molecular Analyses Indicate Separate Importations Rather than Endemic Circulation of Wild Type 1 Poliovirus in Gaza District in 2002 , 2007, Applied and Environmental Microbiology.

[15]  E. Ionides,et al.  Inference for dynamic and latent variable models via iterated, perturbed Bayes maps , 2015, Proceedings of the National Academy of Sciences.

[16]  S. Cairncross,et al.  Water, sanitation, hygiene and enteric infections in children , 2013, Archives of Disease in Childhood.

[17]  A. Vincent,et al.  Rapid Group-, Serotype-, and Vaccine Strain-Specific Identification of Poliovirus Isolates by Real-Time Reverse Transcription-PCR Using Degenerate Primers and Probes Containing Deoxyinosine Residues , 2009, Journal of Clinical Microbiology.

[18]  C. Cobelli,et al.  Parameter and structural identifiability concepts and ambiguities: a critical review and analysis. , 1980, The American journal of physiology.

[19]  Y. Manor,et al.  Molecular epidemiology of silent introduction and sustained transmission of wild poliovirus type 1, Israel, 2013. , 2014, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[20]  Y. Manor,et al.  Development and validation of a real time quantitative reverse transcription-polymerase chain reaction (qRT-PCR) assay for investigation of wild poliovirus type 1-South Asian (SOAS) strain reintroduced into Israel, 2013 to 2014. , 2014, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[21]  Lerato Seakamela,et al.  Diagnostic Assay Development for Poliovirus Eradication , 2017, Journal of Clinical Microbiology.

[22]  N. Grassly The final stages of the global eradication of poliomyelitis , 2013, Philosophical Transactions of the Royal Society B: Biological Sciences.

[23]  K.S.P. Kumar,et al.  On the observability of nonlinear systems: I , 1971 .

[24]  R. Hudson,et al.  A model for the genetic architecture of quantitative traits under stabilizing selection , 2017, 1704.06707.

[25]  Yves F. Atchad'e,et al.  Iterated filtering , 2009, 0902.0347.

[26]  L. Roberts Infectious disease. Israel's silent polio epidemic breaks all the rules. , 2013, Science.

[27]  A. Sabin,et al.  Immunization with live attenuated poliovirus vaccine. , 1961, American journal of diseases of children.

[28]  Haley R Pipkins,et al.  Polyamine transporter potABCD is required for virulence of encapsulated but not nonencapsulated Streptococcus pneumoniae , 2017, PloS one.

[29]  R. Haque,et al.  Household Environmental Conditions Are Associated with Enteropathy and Impaired Growth in Rural Bangladesh , 2013, The American journal of tropical medicine and hygiene.

[30]  Y. Manor,et al.  Neurovirulent Vaccine-Derived Polioviruses in Sewage from Highly Immune Populations , 2006, PloS one.

[31]  J. Gog,et al.  A status-based approach to multiple strain dynamics , 2002, Journal of mathematical biology.

[32]  E. Hargreaves Epidemiology of poliomyelitis. , 1949, Lancet.

[33]  S. Wassilak,et al.  Using Acute Flaccid Paralysis Surveillance as a Platform for Vaccine-Preventable Disease Surveillance , 2017, The Journal of infectious diseases.

[34]  Y. Manor,et al.  Estimation of polio infection prevalence from environmental surveillance data , 2017, Science Translational Medicine.

[35]  S. Tegegne,et al.  Contribution of Environmental Surveillance Toward Interruption of Poliovirus Transmission in Nigeria, 2012–2015 , 2016, Journal of Infectious Diseases.

[36]  D. Shmueli,et al.  Bedouin Communities in the Negev , 2011 .

[37]  Ursula Klingmüller,et al.  Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood , 2009, Bioinform..

[38]  Daisuke Sano,et al.  Microfluidic Quantitative PCR for Simultaneous Quantification of Multiple Viruses in Environmental Water Samples , 2014, Applied and Environmental Microbiology.

[39]  Suzanne L Robertson,et al.  Identifiability and estimation of multiple transmission pathways in cholera and waterborne disease. , 2013, Journal of theoretical biology.

[40]  Andrew W. Brooks,et al.  Genetic signatures for Helicobacter pylori strains of West African origin , 2017, PloS one.

[41]  M. Pallansch,et al.  Review and Assessment of Poliovirus Immunity and Transmission: Synthesis of Knowledge Gaps and Identification of Research Needs , 2013, Risk analysis : an official publication of the Society for Risk Analysis.

[42]  Y. Manor,et al.  Quantitative multiplex one-step RT-PCR assay for identification and quantitation of Sabin strains of poliovirus in clinical and environmental specimens. , 2018, Journal of virological methods.

[43]  Yann S. Dufour,et al.  Hook length of the bacterial flagellum is optimized for maximal stability of the flagellar bundle , 2018, PLoS biology.

[44]  Y. Manor,et al.  Laboratory challenges in response to silent introduction and sustained transmission of wild poliovirus type 1 in Israel during 2013. , 2014, The Journal of infectious diseases.

[45]  Marisa C Eisenberg,et al.  Determining identifiable parameter combinations using subset profiling. , 2014, Mathematical biosciences.

[46]  Stephen R. Lindemann,et al.  Predicting Species-Resolved Macronutrient Acquisition during Succession in a Model Phototrophic Biofilm Using an Integrated ‘Omics Approach , 2017, Front. Microbiol..

[47]  R. Bassal,et al.  Field study of fecal excretion as a decision support tool in response to silent reintroduction of wild-type poliovirus 1 into Israel. , 2015, Journal of clinical virology : the official publication of the Pan American Society for Clinical Virology.

[48]  Marisa C. Eisenberg,et al.  Dose-response relationships for environmentally mediated infectious disease transmission models , 2017, PLoS Comput. Biol..

[49]  D. Earn,et al.  Multiple Transmission Pathways and Disease Dynamics in a Waterborne Pathogen Model , 2010, Bulletin of mathematical biology.

[50]  Michael H. Cortez,et al.  Distinguishing between Indirect and Direct Modes of Transmission Using Epidemiological Time Series , 2013, The American Naturalist.

[51]  R. Blattner Epidemiology of poliomyelitis , 1955 .

[52]  S. Lowther,et al.  Environmental surveillance for polioviruses in the Global Polio Eradication Initiative. , 2014, The Journal of infectious diseases.

[53]  Y. Manor,et al.  Insidious reintroduction of wild poliovirus into Israel, 2013. , 2013, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[54]  E. Kopel,et al.  Lessons from a public health emergency--importation of wild poliovirus to Israel. , 2014, The New England journal of medicine.

[55]  E. Ionides,et al.  Successes and shortcomings of polio eradication: a transmission modeling analysis. , 2013, American journal of epidemiology.

[56]  Y. Manor,et al.  Modeling the spread of polio in an IPV-vaccinated population: lessons learned from the 2013 silent outbreak in southern Israel , 2016, BMC Medicine.

[57]  R. Sutter,et al.  Vaccine-derived polioviruses. , 2014, The Journal of infectious diseases.

[58]  K. Hanson,et al.  Candida auris: an Emerging Fungal Pathogen , 2017, Journal of Clinical Microbiology.

[59]  K. Thompson,et al.  Health and economic consequences of different options for timing the coordinated global cessation of the three oral poliovirus vaccine serotypes , 2015, BMC Infectious Diseases.

[60]  Walter A Orenstein,et al.  The role of supplementary environmental surveillance to complement acute flaccid paralysis surveillance for wild poliovirus in Pakistan – 2011–2013 , 2017, PloS one.

[61]  M. Tashiro,et al.  A community cluster of influenza A(H1N1)pdm09 virus exhibiting cross-resistance to oseltamivir and peramivir in Japan, November to December 2013. , 2014, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[62]  John H. Seinfeld,et al.  Observability of nonlinear systems , 1972 .

[63]  L. Potena,et al.  Monitoring of Cytomegalovirus (CMV)-Specific Cell-Mediated Immunity in Heart Transplant Recipients: Clinical Utility of the QuantiFERON-CMV Assay for Management of Posttransplant CMV Infection , 2018, Journal of Clinical Microbiology.

[64]  N. Nathanson,et al.  From Emergence to Eradication: The Epidemiology of Poliomyelitis Deconstructed , 2010, American journal of epidemiology.

[65]  D. Earn,et al.  Generality of the Final Size Formula for an Epidemic of a Newly Invading Infectious Disease , 2006, Bulletin of mathematical biology.