Sabin Vaccine Reversion in the Field: a Comprehensive Analysis of Sabin-Like Poliovirus Isolates in Nigeria
Authors: Famulare Michael , Chang Stewart , Iber Jane , Zhao Kun , Adeniji Johnson A. , Bukbuk David , Baba Marycelin , Behrend Matthew , Burns Cara C. , Oberste M. Steven 
Published: October 2015

To assess the dynamics of genetic reversion of live poliovirus vaccine in humans, we studied molecular evolution in Sabin-like poliovirus isolates from Nigerian acute flaccid paralysis cases obtained from routine surveillance. We employed a novel modeling approach to infer substitution and recombination rates from whole-genome sequences and information about poliovirus infection dynamics and the individual vaccination history. We confirmed observations from a recent vaccine trial that VP1 substitution rates are increased for Sabin-like isolates relative to the rate for the wild type due to increased nonsynonymous substitution rates. We also inferred substitution rates for attenuating nucleotides and confirmed that reversion can occur in days to weeks after vaccination. We combine our observations for Sabin-like virus evolution with the molecular clock for VP1 of circulating wild-type strains to infer that the mean time from the initiating vaccine dose to the earliest detection of circulating vaccine-derived poliovirus (cVDPV) is 300 days for Sabin-like virus type 1, 210 days for Sabin-like virus type 2, and 390 days for Sabin-like virus type 3. Phylogenetic relationships indicated transient local transmission of Sabin-like virus type 3 and, possibly, Sabin-like virus type 1 during periods of low wild polio incidence. Comparison of Sabin-like virus recombinants with known Nigerian vaccine-derived poliovirus recombinants shows that while recombination with non-Sabin enteroviruses is associated with cVDPV, the recombination rates are similar for Sabin isolate-Sabin isolate and Sabin isolate–non-Sabin enterovirus recombination after accounting for the time from dosing to the time of detection. Our study provides a comprehensive picture of the evolutionary dynamics of the oral polio vaccine in the field.

IMPORTANCE
The global polio eradication effort has completed its 26th year. Despite success in eliminating wild poliovirus from most of the world, polio persists in populations where logistical, social, and political factors have not allowed vaccination programs of sustained high quality. One issue of critical importance is eliminating circulating vaccine-derived polioviruses (cVDPVs) that have properties indistinguishable from those of wild poliovirus and can cause paralytic disease. cVDPV emerges due to the genetic instability of the Sabin viruses used in the oral polio vaccine (OPV) in populations that have low levels of immunity to poliovirus. However, the dynamics responsible are incompletely understood because it has historically been difficult to gather and interpret data about evolution of the Sabin viruses used in OPV in regions where cVDPV has occurred. This study is the first to combine whole-genome sequencing of poliovirus isolates collected during routine surveillance with knowledge about the intrahost dynamics of poliovirus to provide quantitative insight into polio vaccine evolution in the field.