Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
npj Vaccines volume 8, Article number: 137 (2023)
2
Metrics details
There is an increasing burden of circulating vaccine-derived polioviruses (cVDPVs) due to the continued use of oral poliovirus vaccine (OPV). However, the informativeness of routine OPV VP1 sequencing for the early identification of viruses carrying virulence-associated reversion mutations has not been directly evaluated in a controlled setting. We prospectively collected 15,331 stool samples to track OPV shedding from children receiving OPV and their contacts for ten weeks following an immunization campaign in Veracruz State, Mexico and sequenced VP1 genes from 358 samples. We found that OPV was genetically unstable and evolves at an approximately clocklike rate that varies across serotypes and by vaccination status. Overall, 61% (11/18) of OPV-1, 71% (34/48) OPV-2, and 96% (54/56) OPV-3 samples with available data had evidence of a reversion at the key 5’ UTR attenuating position and 28% (13/47) of OPV-1, 12% (14/117) OPV-2, and 91% (157/173) OPV-3 of Sabin-like viruses had ≥1 known reversion mutations in the VP1 gene. Our results are consistent with previous work documenting rapid reversion to virulence of OPV and underscores the need for intensive surveillance following OPV use.
The global use of oral poliovirus vaccine (OPV), a live, attenuated virus, has been essential in reducing the burden of poliomyelitis caused by wild poliovirus by 99.9% since 19881,2. The vaccine covers all three wild poliovirus serotypes, is inexpensive, easy to deliver, and—because it is a live virus—is transmissible from vaccinated children onwards, extending the reach of vaccination campaigns. However, OPV vaccination is risky—the vaccine virus is genetically unstable and continues to evolve in vaccinated individuals and their contacts. If OPV mutates at positions responsible for attenuation of the vaccine virus, it can regain virulence. OPV can regain virulence within a vaccinated individual, causing vaccine-associated paralytic poliomyelitis, though rarely3. OPV may also revert to virulence over transmission chains, generating outbreaks of circulating vaccine derived polioviruses (cVDPVs).
To reduce the risks of cVDPVs, in May 2016, trivalent OPV vaccines (containing serotypes 1, 2, and 3) were replaced with a bivalent OPV (containing serotypes 1 and 3), following World Health Organization recommendations. The WHO recommended routine immunization with at least one dose of inactivated poliovirus vaccine (IPV) prior to the switch to bivalent OPV, to increase protection against paralytic disease4. However, access to IPV was limited, leaving large populations of children immunologically naïve to poliovirus type 2. Since the global switch to bivalent OPV, the majority of paralytic poliovirus cases are now caused by cVDPVs, not wild poliovirus, the majority of which are serotype 2 viruses (cVDPV2). From 2016 to 2022, type 2 cVDPVs emerged 68 times in 34 countries, leading to more than 1500 cases of paralytic polio5 (Fig. 1).
a Confirmed cases of paralytic polio were reported to the World Health Organization from 2000-2022. Color indicates source: circulating vaccine derived poliovirus (cVDPV) and wild poliovirus cases. b Map of countries reporting cVDPV cases from 2020-2022; fill indicates number of cVDPV cases identified. Data from the World Health Organization.
The GPEI recommends vaccination with OPV as part of an outbreak response to induce mucosal immunity and interrupt transmission. While in the past, monovalent OPV-2 (mOPV-2) was used to respond to cVDPV2 outbreaks, a more genetically stable replacement, type 2 novel OPV (nOPV-2)6,7, was recommended by the WHO under an Emergency Use Listing in 20208. Since March 2021, nOPV-2 has been the primary vaccine used to respond to cVDPV2 outbreaks. Yet even with the development of a new vaccine, mOPV-2 has not been completely replaced and the risk of cVDPVs remains a serious global health concern.
Outbreak responses are guided by routine testing of cases of acute flaccid paralysis in children less 15, environmental surveillance in more than 30 high-risk countries, and sequencing of the poliovirus VP1 gene9. The Global Polio Eradication Initiative (GPEI) defines VDPVs by their divergence from OPV in the VP1 gene9 and cVDPVs as VDPVs with evidence of transmission10,11. The identification of new cVPDVs results in an epidemiological investigation and an immunization response9.
Previous genomic studies of type 2 cVDPVs and OPV-2 shedding in vaccinated children and contacts have identified strong within-host selection and tight transmission bottlenecks that govern the early evolution of type 2 viruses away from the Sabin vaccine strain along a predictable pathway to virulence12,13. Yet previous studies have not yet integrated detailed molecular and epidemiological data to assess the use of VP1 sequencing for the early identification of viruses posing a public health risk. Because mOPV2 as well as tOPV continue to be used in outbreak responses, we need a better understanding of the risks associated with the reversion to virulence and transmission of viruses with virulence-associated mutations, in particular, in the months directly following vaccination campaigns. Given that global surveillance and outbreak responses rely on VP1 sequencing, we need a better understanding of how informative VP1 diversity is for characterizing the risk of reversion to virulence and transmission of virulent viruses.
Here, we tracked the evolution and transmission of the three OPV serotypes during a prospective study of OPV shedding following a vaccine campaign in semi-rural communities in Veracruz State, Mexico14. We sequenced OPV VP1 genes from stool sampled from vaccinated children, their household contacts, and community members in unvaccinated households for 10 weeks following the vaccination campaign. This longitudinal collection provides a rare opportunity to observe the early evolution and transmission of the three OPV serotypes following vaccination campaigns and to evaluate the current use of VP1 sequencing as an epidemiological tool.
During a prospective study of OPV shedding following a National Health Week (NHW) vaccination campaign in Veracruz State, Mexico, we collected 15,331 stool samples over 10 weeks from children receiving the OPV vaccine, household contacts, and community members (Fig. 2). In our previous study of OPV shedding and transmission, we found that among 148 vaccinated children, 380 household contacts, and 1124 unvaccinated community contacts, 78%, 18%, and 7%, respectively, shed OPV15. We found that transmission of OPV began as soon as 1 day following OPV vaccination and, while most community members shed transmitted OPV for less than a week, we observed shedding among community members as long as 71 days after the vaccination campaign. We also observed more frequent shedding of OPV by community members in communities with the highest vaccination rates15.
a Map of Mexico with Veracruz State highlighted and Orizaba city, where the study took place, indicated with a red point. b We conducted a prospective observational study of OPV viral shedding following a vaccination campaign. 148 children were enrolled and vaccinated with trivalent OPV. The stool was collected from NHW vaccinees, unvaccinated household members, and community members in unvaccinated households. We extracted RNA from stool samples, tested samples with OPV serotype-specific RT-PCR, and Sanger sequenced positive samples.
From these samples, 551 samples were positive for at least one OPV serotype, including 267 OPV-1, 402 OPV-2, and 317 OPV-3 positive samples. We sequenced 358 high quality VP1 genes from 174 individuals, including 18.4% (49/267) positive OPV-1 samples, 33.8% (136/402) positive OPV-2 samples, and 54.5% (173/317) OPV-3 samples. Sequence data was available for 99 children who were vaccinated with OPV during the NHW vaccination campaign, whom we refer to as NHW vaccinees; and 72 community members, including 41 household contacts of NHW vaccinees. Three individuals did not have a known vaccination status. For 84 infections (unique individuals and serotypes), longitudinal samples were available at 2–5 sampling points. The prospective study was conducted in three semi-rural villages; samples include 106 individuals from Capoluca, where 70% of children less than five in participating households were vaccinated during the NHW campaign; 49 from Campo Grande, where 30% were vaccinated; and 19 individuals from Tuxpanguillo, where 10% were vaccinated.
We observed a rapid rise of VP1 diversity in all three serotypes over weeks following OPV vaccination. The most common haplotype for both OPV-1 and OPV-2 in each sampling location represents the Sabin vaccine strain (Fig. 3) and we observed a star-like pattern in the haplotype networks, where most haplotypes were directly connected to the Sabin haplotype. Outside of the dominant haplotype, most haplotypes are represented by a single sample.
For each OPV serotype and sampling site, haplotype networks of OPV VP1 represent sampled genetic diversity. Nodes indicate identical haplotypes (VP1 consensus sequences) and node size indicates the number of samples sharing a haplotype. Points on edges between nodes indicate the SNP distance between haplotypes and node colors indicate an individual’s vaccination status. Labels indicate OPV serotype and study site. Site C Capoluca, had 70% vaccination coverage of eligible children, site CG Campo Grande: 30%, site T Tuxpanguillo: 10%. One sample with a distant haplotype (>14 SNPs distant to the Sabin vaccine) was removed from OPV-2 site C, to aid visualization.
In contrast, for OPV-3, the dominant haplotype differs by a single nucleotide polymorphism (SNP) from the vaccine strain, a previously described attenuating site, C2493U (Fig. 3). Again, there is a star-like pattern of OPV-3 diversity, with many branches stemming directly from the dominant haplotype. As with OPV-2 and OPV-3, most haplotypes outside the dominant haplotype are represented by a single sample.
The Global Polio Eradication Initiative (GPEI) defines VDPVs by their divergence from the Sabin vaccine in the VP1 gene (>1% divergent or ≥10 substitutions for serotypes 1 and 3, 0.6% divergent of ≥6 substitutions for serotype 2)9. The majority, 99% (338/340), of all samples sequenced following vaccination fell under the thresholds for a VDPV. However, one sample from an NHW vaccinee, collected 4 days after vaccination met the definition of a VDPV2 (VP1 genetic distance 14 from Sabin 2), and one sample from a household contact sampled 30 days following vaccination met the definition of a VDPV1 (10 SNPs distant from Sabin 1) (Fig. 4).
Genetic distance of the OPV VP1 gene to the Sabin vaccine strain versus time from vaccination for NHW vaccinees (blue) and not vaccinated study participants, including household contacts and unvaccinated community members (red), for each OPV serotype (a-c). Black lines and grey shading of a linear model for distance to the Sabin strain. Days following vaccination is measured as days from the first vaccination (for NHW vaccinees), day from the first household vaccination (for household members), and days from the first community vaccination (for other community members). Red dashed lines indicate the genetic distance threshold for a VDPV: 10 SNPs distant from the Sabin vaccine for serotypes 1 and 3 and 6 SNPs distant for serotype 2.
We then tested whether OPV accumulated substitutions in a clock-like and, therefore, if OPV divergence from Sabin vaccine strains could be used to estimate the duration of viral circulation in the months following a vaccination campaign. We were not able to calculate a true molecular clock rate because of the limited observation period and because of our inclusion of sequences that do not represent a single genetic lineage. Genetic divergence of OPV-1 was significantly associated with time following vaccination in both NHW vaccinees (p = 0.0099, Pearson’s r = 0.42) and community members (p < 0.001, Pearson’s r = 0.90) (Fig. 4). When considering all samples together, the observed mutations in OPV-1 corresponded to a VP1 gene evolutionary rate of 1.2 × 10-4 substitutions per site per day. Genetic divergence of OPV-2 was not significantly associated with time following vaccination in NHW vaccinees (p = 0.49, Pearson’s r = −0.075) and was positively correlated with time among community members, though not significantly so (p = 0.075, Pearson’s r = 0.31) (Fig. 4). OPV-3 divergence was positively associated with time following vaccination for community members (p = 0.0058, Pearson’s r = 0.48), but not NHW vaccinees (p = 0.39, Pearson’s r = −0.073)
We reasoned that mutations at known attenuating positions may be under strong selection and could affect an estimate of the molecular clock. However, excluding the known attenuating sites in the VP1 gene did not substantially change observed associations between genetic divergence and time (Supplementary Fig. 1).
VP1 genetic diversity is also used by the GPEI to determine if a VDPV is genetically linked to previously collected environmental or clinical samples, which makes it a cVDPV, and whether viruses are linked to a previously identified cVDPV emergence or constitute a new emergence.
To evaluate the genetic structure of VP1, we compared pairwise genetic distances across samples. As expected with the relatively short duration of sampling following vaccination campaigns, 27.9% (644/2304) of OPV-1; 41.9% (6693/16,616) of OPV-2; and 41.4% (12,272/29,628) of OPV-3 VP1 pairs of sequences were identical. The mean pairwise distance was 2.5 SNPs for OPV-1, 1.4 SNPs for OPV-2, and 1.3 SNPs for OPV-3 (Fig. 5).
The proportion of pairs of OPV VP1 sequences within a given, binned pairwise genetic distance for pairs of samples collected longitudinally from the same individual, from the same household, different households within the same community, or outside the community, i.e. pairs of samples from different communities. Facets indicate OPV serotype.
In a general linear model, the pairwise genetic distance between VP1 samples was elevated for individuals from the same household (aOR: 1.31; 95% CI: 1.12–1.52), from the same community (aOR: 1.38; 95% CI: 1.26–1.53), and from different communities (aOR: 1.53; 95% CI: 1.39-1.769), compared to samples collected from an individual. As expected, longitudinal samples collected from the same individual were largely identical. Pairwise genetic distances differed across serotypes: compared to OPV-1, we observed decreased pairwise genetic distances for both OPV-2 (aOR: 0.59; 95% CI: 0.57–0.61) and OPV-3 (aOR: 0.54; 95% CI: 0.52–0.55).
We then examined known mutations at key attenuating sites in the vaccine virus, associated with the reversion of vaccine viruses to virulence. First, as stool may contain viruses with both revertant and non-revertant alleles at the same locus, we assessed the proportion of viral RNA in each sample with a revertant mutation at the major attenuating position in the 5’ UTR, or the revertant proportion, with an allele-specific qPCR test16. Of the 122 samples with available results, 99 (81%) had a signal of reversion to virulence (revertant proportion > 0), and in 82 samples (67%), the reversion proportion was ≥ 50%. Among samples with complete assay information, we found evidence of a 5’ UTR reversion in 82% (87/106) of vaccines, 85% of household contacts (11/13), and 33% (1/3) of other community members. We identified transmission of the 5’ UTR reversion mutation in six independent households in which the revertant mutation was first identified in NHW vaccinees and subsequently identified in household members (4 households) or concurrently identified in both NHW vaccinees and household members (2 households). Transmission of the 5’ UTR reversion mutation occurred for all three serotypes.
We then examined additional mutations within the VP1 gene associated with reversion to virulence. Of three known virulence attenuating substitutions in the OPV-1 VP1 gene, 16% (8/49) of samples had a A2749G (Sabin to virulent wildtype) mutation, including 4 NHW vaccinees and 2 household contacts. The mutation first appeared 14 days following vaccination in a vaccinated child and contact in one house and in another vaccinated child independently and rose in frequency, although there were limited OPV-1 samples available 4+ weeks after vaccination. It was observed in stool up to 52 days after the vaccination campaign. Three samples with the A2749G mutation, including two NHW vaccinees and one household member, also had evidence of fixation of the reversion mutation in the 5’ UTR (reversion proportion of 100%). A single sample shed by a vaccinated child harbored a A2795G substitution, and no U2879C attenuating mutations occurred (Fig. 6).
For each serotype, points indicate the revertant proportion, the proportion of samples with a mutation at a key attenuating position in OPV in the weeks following vaccination. Bars indicate the total number of sequences available for that sampling week. Each serotype has unique attenuating sites.
Among OPV-2 samples, the known VP1 reversion mutation, U2909C/A, occurred in 4.4% (6/136) of samples, in four NHW vaccinees. The substitution was first observed 7 to 14 days after vaccination; two children shed OPV-2 carrying this reversion mutation on multiple days. The mutation was observed in stool up to 21 days following the vaccination campaign, though again, there were limited OPV-2 samples available 4+ weeks after vaccination. No individuals had mutations present at the adjacent site 2908, also associated with reversion7. Four of these samples also had evidence of a reversion mutation in the 5’ UTR.
Among OPV-3 samples, the majority, 96.5% (167/173), had a C2493U eversion substitution and an additional 11 samples had a mixed allele or no allele call at that site, possibly evidence of mixed infection. Samples with a C2493/U mutation included 95% (89 of 94) of individuals with OPV-3 sequences available, including 62 NHW vaccinees, 19 household members, and 8 individuals outside of vaccinated households. The mutation first appeared one day following the vaccination in NHW vaccinees, household members, and individuals outside vaccinated households and it was observed in stool up to 72 days after the vaccination campaign from a vaccinated child and 70 and 54 days from virus shed in a household member and community member, respectively. 35% (58/167) of these samples also had evidence of a reversion mutation in the 5’ UTR.
While rapid reversion to virulence has been documented in OPV, the risk of transmission of OPV with virulence-associated mutations following vaccination campaigns is less clear. Further, global poliovirus surveillance uses molecular epidemiology to directly guide outbreak responses. However, the informativeness of the routinely sequenced marker gene, VP1, in the early identification of VDPVs posing a public health risk has not been evaluated in a controlled setting.
Here, we leveraged OPV sampled for 10 weeks following an immunization campaign in a population previously vaccinated with IPV to assess the epidemiological information within VP1 samples. We found that all three OPV serotypes were genetically unstable in the weeks following vaccination campaigns, consistent with previous reports. We identified a linear relationship between genetic divergence from vaccine strains and sampling time (though this signal differed by vaccine status) and significant genetic structure. Consistent with prior reports, we found evidence of frequent reversion substitutions among children who had received the OPV vaccine. Importantly, reversion mutations were frequently also found among community members who had not received OPV as part of the National Immunization Week, evidence that revertant viruses may be frequently transmitted onwards.
Polioviruses are among the fastest evolving RNA viruses, and we observed that the rapid generation of genetic diversity in NHW vaccinees and their household and community contacts contains epidemiological information. Robust molecular clock estimates require samples from a well-defined virus lineage and over a significant sampling period, which were not available here, and therefore, we report trends observed in the genetic divergence from the Sabin vaccine over time rather than precise estimates of the OPV molecular clock. We found that the early evolution of OPV is largely clocklike, with variation in mutation rate across serotypes. The lack of a signal of clocklike evolution of OPV-2 and OPV-3 in NHW vaccinees likely reflects the shorter sampling period for NHW vaccinees; unvaccinated household and community contacts were infected later in the study by viruses that had already passed through at least one host.
Though we were not able to calculate a true molecular clock rate, our observed OPV-1 VP1 substitution rate of 1.2 × 10−4 substitutions/site/day (including all samples) was similar to previous estimates of OPV-1 evolution 0.83 × 10−4 substitutions/site/day17 (estimated from three years of surveillance for acute flaccid paralysis). Previously, the molecular clock for wild poliovirus type 1 was estimated to be 0.28 × 10−4 substitutions/site/day, from a set of 31 viruses sampled over ten years18 3.12 × 10−5 substitutions/site/day in samples shed from a chronically infected immune-compromised individual19. The differences in substitution rate may reflect the different sampling durations, as shorter sampling durations often yield higher pathogen substitution rate estimates, likely because purifying selection has not yet removed transient deleterious mutations that are still present prior to the effect of purifying selection20,21, in addition to different selective pressures in these different sampling contexts.
We observed two samples with unexpectedly high levels of genetic divergence that met the definition of VDPVs. Likely, these samples represent OPV infection from sources outside of the focal NHW vaccination campaign, indicating prolonged replication or circulation elsewhere.
We observe moderate genetic structure in VP1 sequences, with increasing genetic divergence between samples from different communities compared to that sampled within the same community, household, or longitudinally within an individual. Whole genome sequences and additional epidemiological information about the timing of infection and contacts could improve the inference of transmission networks.
The rapid reversion to virulence of Sabin vaccines is well documented22,23. While each serotype has key attenuating mutations in the 5’ UTR, additional mutations along the genome, including in the VP1 gene, stabilize the attenuated phenotype. In OPV-2, a predictable evolutionary route towards virulence includes gatekeeper mutations in the 5’UTR and the VP1 gene, followed by recombination events, followed by additional stabilizing mutations12,13. Most samples had evidence of the major 5’ UTR reversion mutation, including samples 82% of vaccinee samples, 85% from household contacts, and 33% from other community members. Importantly, we also found evidence of household transmission of the 5’ UTR mutation for all three serotypes, indicating that samples that do not meet the genetic distance threshold for a VDPV are both virulent and transmissible.
We additionally found a high frequency of attenuating substitutions in VP1, with substitutions appearing early after vaccination. The higher frequency of OPV-2 and OPV-3 reversion mutations we observed could reflect the greater observed genetic stability of OPV-1, which differs from the virulent wild-type progenitor virus by 59 substitutions, compared to OPV-2 and OPV-3 which are more closely related to the virulent ancestor24. Our group previously found that 97% of OPV-3 recipients shed viruses that had a site 472 reversion mutation within two weeks following vaccination25,26. The instability of OPV-3 has epidemiological consequences: OPV-3 causes the majority of vaccine-associated paralytic poliomyelitis (VAPP), paralysis attributable to the vaccine itself, followed by OPV-2 and then OPV-13. However, previous studies have reported heterogeneity at the OPV-3 C2493U attenuating mutation in vaccine stocks27. We did not have access to vaccine stocks and therefore cannot determine if the near-fixation of an attenuating mutation we observed was due to heterogeneity in the vaccine or de novo evolution in NHW vaccinees and subsequent onward transmission.
OPV remains the primary tool used to control poliovirus outbreaks. Novel OPV2 (nOPV-2), a genetically stable replacement for mOPV-2, was granted Emergency Use Listing in 2020 and first used in outbreak vaccination responses in 20218. Since then, over 600 million doses of nOPV2 have been administered and the WHO Strategic Advisory Group of Experts on Immunization recommended the transition from its initial use period to wider use because of its documented safety and effectiveness in inducing strong immune responses in a cVDPV2 outbreak response8,28,29. The increasing replacement of mOPV-2 is expected to reduce the risk of cVDPVs in the future; however, several questions remain about how nOPV-2 can be optimally used for outbreak control30.
Although the majority of OPV samples in our study collected were defined as Sabin-like viruses, rather than VDPVs, many of the Sabin-like viruses had rapidly lost attenuating mutations at previously described positions. While VP1 sequence divergence may be useful in estimating the duration of circulation or predicting transmission clusters, circulating OPV that does not yet meet GPEI definitions for a cVDPV, can be virulent and can pose a public health risk in populations without prior immunity. For example, Sabin-like OPV-2 viruses, which were 3 SNPs distant to the Sabin virus VP1 sequence, were previously found to have caused a cluster of acute flaccid paralysis cases in an orphanage in the Altai Region of Russia31. All outbreak genomes shared 15 substitutions compared to the Sabin virus, including the loss of attenuation mutations A481G and U2909C.
The global increase of cVDPV2 outbreaks has increased dramatically since 2016, when Sabin 2 was removed from the trivalent poliovirus vaccine5. Our results are consistent with previous work on the genetic instability of OPV and rapid reversion to virulence and underscore the urgency of heightened surveillance during and following vaccination campaigns.
We conducted a prospective study of OPV shedding and transmission following a trivalent OPV vaccination campaign during Mexico’s February 2015 National Health Week (NHW) in three semi-rural populations in Orizaba, Veracruz State14 (Fig. 2). Wild poliovirus has not been reported since 1994 in Mexico. Mexico’s national vaccine program has included routine IPV as well as two annual OPV vaccination campaigns since 2007 and therefore provides a model for transmission of OPV in a setting with routine IPV immunization, an epidemiological context that will increasingly exist as IPV replaces OPV and OPV is only used as a tool for cVDPV or wild poliovirus outbreak response.
We randomized 466 households in three communities to the vaccine and no-vaccine arm of our study. To measure the impact of vaccination rate on the frequency of household and community shedding, we randomized households so that community-level OPV vaccination rates were 10%, 30%, and 70% in the three communities14. In total, 148 children (8.5% of the population) were vaccinated with trivalent OPV during the NHW. 148 children receiving the vaccine during the NHW (NHW vaccinees), 380 household contacts, and 1124 community members from households that did not receive OPV during the NHW were followed prospectively and stool was collected from all members of sampled households 1 day prior to the National Health Week OPV immunization campaign, and 1, 4, 7, 10, 14, 21, 28, 51, and 71 days following, for a total of 10 samples per person14.
Study nurses collected immunization histories and additional demographic and clinical metadata from all study participants. We report rates of OPV shedding and transmission among child vaccines, household contacts, and community contacts in prior studies15; the focus of this analysis was on the molecular evolution of tOPV following vaccination campaigns.
The study was approved by the Stanford University Institutional Review Board (Protocol #31546), the Comité de Etica, Bioseguridad e Investigación of the Instituto Nacional de Salud Pública (CI: 1260, No. 1581), and the Instituto Veracruzano para la Formación e Investigación en Salud (SESVER/IVEFIS//SIS/DIB/0109/02014, classification 15 S).
We extracted RNA from frozen stool samples using an RNeasy minikit (Qiagen, Valencia, CA) according to the manufacturer’s protocol. Positive and negative controls were included in each set of extractions. Extracted RNA preparations were stored at −80°C and tested samples for OPV with serotype-specific RT-qPCR14,16. We Sanger sequenced the VP1 gene using the Y7 and Q8R primers32 and generated consensus sequences by aligning forward and reverse sequences to serotype-specific reference genomes in Geneious v7.
We additionally used a previously described allele-specific RT-qPCR assay to assess the proportion of viral RNA in each sample with the major 5’ UTR reversion mutation, or the revertant proportion16. We combined five microliters of the reverse transcription reaction mixture with 10 μl TaqMan universal PCR master mix without uracil-N-glycosylase (Applied Biosystems, Foster City, CA), 0.05 μl TaqMan probe (Applied Biosystems, Foster City, CA), 0.9 μl each of forward and reverse primers, and 3.15 μl distilled water. Cycling conditions were as described previously25. Samples that returned average cycle threshold (CT) values of >35 from duplicate assays were considered negative, due to an increased number of false-positive results above this threshold.
We calculated the revertant proportion from the allele-specific RT-qPCR cycle threshold values as previously described16: 2−revCT/(2−revCT + 2−nonrevCT), where revCT and nonrevCT are revertant and nonrevertant CT values, respectively.
We aligned VP1 gene sequences of each OPV serotype to the corresponding Sabin vaccine sequence (GenBank Accession Numbers: AY184219, AY184220, AY184221) with MAFFT v.733, with the algorithm recommended for closely related viral sequences (https://mafft.cbrc.jp/alignment/software/closelyrelatedviralgenomes.html), preserving the alignment length of the reference sequence. We used the default alignment strategy and a nucleotide sequence scoring matrix of 1PAM / κ = 2, recommended for closely related DNA sequences.
We quantified pairwise distances between samples with the R package ape (pairwise deletion = TRUE)34. For household contacts and community members who were not vaccinated with OPV, we measured the time since vaccination as the time since vaccination of the household member or the earliest vaccination given in the community, respectively. We visualized haplotype networks built from minimum spanning trees with the R package pegas35.
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
The full-length VP1 sequences used in our study are publicly available on GenBank (accession numbers OQ714877 – OQ715012, OQ715013 – OQ715185, OQ715186 – OQ715234). The references sequences used for sequence alignment are publicly available on GenBank (AY184219, AY184220, AY184221). Multiple sequence alignments and metadata to reproduce figures and analyses are publicly available on GitHub (https://github.com/ksw9/opv-consensus). We accessed World Health Organization poliovirus case data from the WHO extranet (https://extranet.who.int/polis/public/CaseCount.aspx), including both global and country-level reported cases from 2000–2022.
Code for analysis and figures is publicly available on GitHub (https://github.com/ksw9/opv-consensus).
Khan, F. et al. Progress toward polio eradication — Worldwide, January 2016–March 2018. Morb. Mortal. Wkly. Rep. 67, 524 (2018).
Article Google Scholar
Jorba, J. et al. Update on vaccine-derived poliovirus outbreaks — Worldwide, January 2018–June 2019. Morb. Mortal. Wkly. Rep. 68, 1024 (2019).
Article Google Scholar
Platt, L. R., Estivariz, C. F. & Sutter, R. W. Vaccine-associated paralytic poliomyelitis: a review of the epidemiology and estimation of the global burden. J. Infect. Dis. 210, S380–S389 (2014).
Article PubMed Google Scholar
Chumakov, K. M., Brechot, C., Gallo, R. C. & Plotkin, S. Choosing the right path toward polio eradication. NEJM 388, 577–579 (2023).
Article PubMed Google Scholar
Macklin, G. R. et al. Epidemiology of type 2 vaccine-derived poliovirus outbreaks between 2016 and 2020. Vaccine (2022) https://doi.org/10.1016/j.vaccine.2022.08.008.
Wahid, R. et al. Assessment of genetic changes and neurovirulence of shed Sabin and novel type 2 oral polio vaccine viruses. NPJ Vaccines 6, 1–11 (2021).
Article Google Scholar
Wahid, R. et al. Genetic and phenotypic stability of poliovirus shed from infants who received novel type 2 or Sabin type 2 oral poliovirus vaccines in Panama: an analysis of two clinical trials. Lancet Microbe 3, e912–e921 (2022).
Article CAS PubMed PubMed Central Google Scholar
The World Health Organization. Standard operating procedures: responding to a poliovirus event or outbreak. (2022).
Global Polio Eradication Initiative. Classification and reporting of vaccine-derived polioviruses (VDPV) – GPEI guidelines. (2016).
Jorgensen, D., Pons-Salort, M., Shaw, A. G. & Grassly, N. C. The role of genetic sequencing and analysis in the polio eradication programme. Virus Evol. 6, veaa040 (2020).
Article PubMed PubMed Central Google Scholar
Jorba, J. Phylogenetic analysis of poliovirus sequences. Methods Mol. Biol. 1387, 227–237 (2016).
Article CAS PubMed Google Scholar
Valesano, A. L. et al. The early evolution of oral poliovirus vaccine is shaped by strong positive selection and tight transmission bottlenecks. Cell Host Microbe 29, 32–43.e4 (2021).
Article CAS PubMed PubMed Central Google Scholar
Stern, A. et al. The evolutionary pathway to virulence of an RNA Virus. Cell 169, 35–46.e19 (2017).
Article CAS PubMed PubMed Central Google Scholar
Sarnquist, C. et al. Protocol paper: oral poliovirus vaccine transmissibility in communities after cessation of routine oral poliovirus vaccine immunization. Clin. Infect. Dis. 67, S115–S120 (2018).
Article PubMed PubMed Central Google Scholar
Altamirano, J. et al. Characterization of household and community shedding and transmission of oral polio vaccine in mexican communities with varying vaccination coverage. Clin. Infect. Dis. 67, S4–S17 (2018).
Article CAS PubMed PubMed Central Google Scholar
Troy, S. B. et al. Use of a novel real-time PCR assay to detect oral polio vaccine shedding and reversion in stool and sewage samples after a mexican national immunization day. J. Clin. Microbiol 49, 1777–1783 (2011).
Article PubMed PubMed Central Google Scholar
Famulare, M. et al. Sabin vaccine reversion in the field: a comprehensive analysis of sabin-like poliovirus isolates in Nigeria. J. Virol. 90, 317–331 (2016).
Article CAS PubMed Google Scholar
Jorba, J., Campagnoli, R., De, L. & Kew, O. Calibration of multiple poliovirus molecular clocks covering an extended evolutionary range. J. Virol. 82, 4429–4440 (2008).
Article CAS PubMed PubMed Central Google Scholar
Yang, C.-F. et al. Intratypic recombination among lineages of type 1 vaccine-derived poliovirus emerging during chronic infection of an immunodeficient patient. J. Virol. 79, 12623–12634 (2005).
Article CAS PubMed PubMed Central Google Scholar
Duchêne, S., Holmes, E. C. & Ho, S. Y. W. Analyses of evolutionary dynamics in viruses are hindered by a time-dependent bias in rate estimates. Proc. R. Soc. B: Biol. Sci. 281, 20140732 (2014).
Article Google Scholar
Holmes, E. C. Molecular clocks and the puzzle of RNA virus origins. J. Virol. 77, 3893 (2003).
Article CAS PubMed PubMed Central Google Scholar
Duintjer Tebbens, R. J. et al. Oral poliovirus vaccine evolution and insights relevant to modeling the risks of circulating vaccine-derived polioviruses (cVDPVs). Risk Anal. 33, 680–702 (2013).
Article PubMed PubMed Central Google Scholar
Burns, C. C., Diop, O. M., Sutter, R. W. & Kew, O. M. Vaccine-derived polioviruses. J. Infect. Dis. 210, S283–S293 (2014).
Article PubMed Google Scholar
Kew, O. M., Sutter, R. W., De Gourville, E. M., Dowdle, W. R. & Pallansch, M. A. Vaccine-derived polioviruses and the endgame strategy for global polio eradication. Annu. Rev. Microbiol 59, 587–635 (2005).
Article CAS PubMed Google Scholar
Gnanashanmugam, D. et al. Shedding and reversion of oral polio vaccine Type 3 in Mexican vaccinees: comparison of mutant analysis by PCR and enzyme cleavage to a real-time PCR assay. J. Clin. Microbiol 45, 2419 (2007).
Article CAS PubMed PubMed Central Google Scholar
Martinez, C. V. et al. Shedding of Sabin poliovirus type 3 containing the nucleotide 472 uracil-to-cytosine point mutation after administration of oral poliovirus vaccine. J. Infect. Dis. 190, 409–416 (2004).
Article CAS PubMed Google Scholar
Neverov, A. & Chumakov, K. Massively parallel sequencing for monitoring genetic consistency and quality control of live viral vaccines. Proc. Natl Acad. Sci. USA 107, 20063–20068 (2010).
Article CAS PubMed PubMed Central Google Scholar
Mirzoev, A. et al. Assessment of serological responses following vaccination campaigns with Type 2 novel oral polio vaccine: a population-based study in Tajikistan in 2021. Lancet Glob. Health 10, e1807–e1814 (2022).
Article CAS PubMed PubMed Central Google Scholar
GPEI. Two years since rollout of novel oral polio vaccine type 2 (nOPV2): How’s it all working out? https://polioeradication.org/news-post/two-years-since-rollout-of-novel-oral-polio-vaccine-type-2-nopv2-hows-it-all-working-out/ (2023).
Thompson, K. M. Effectiveness of a new vaccine for outbreak response and the increasingly complicated polio endgame. Lancet Glob. Health 10, e1697–e1698 (2022).
Article CAS PubMed Google Scholar
Korotkova, E. A. et al. A cluster of paralytic poliomyelitis cases due to transmission of slightly diverged Sabin 2 vaccine poliovirus. J. Virol. 90, 5978 (2016).
Article CAS PubMed PubMed Central Google Scholar
Kilpatrick, D. R. et al. Poliovirus serotype-specific VP1 sequencing primers. J. Virol. Methods 174, 128–130 (2011).
Article CAS PubMed Google Scholar
Nakamura, T., Yamada, K. D., Tomii, K. & Katoh, K. Parallelization of MAFFT for large-scale multiple sequence alignments. Bioinformatics 34, 2490–2492 (2018).
Article CAS PubMed PubMed Central Google Scholar
Paradis, E. & Schliep, K. Ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35, 526–528 (2019).
Article CAS PubMed Google Scholar
Paradis, E. pegas: an {R} package for population genetics with an integrated-modular approach. Bioinformatics 26, 419–420 (2010).
Article CAS PubMed Google Scholar
Download references
This project was funded by the National Institutes of Health 5R21AI148810 to Y.M; K.S.W. was funded by a Thrasher Early Career Award and a Stanford Child Health Research Institute postdoctoral fellowship.
Division of Epidemiology, University of Utah, Salt Lake City, UT, 84105, USA
Katharine S. Walter
Department of Epidemiology and Population Health, Stanford University School of Medicine, Stanford, CA, 94305, USA
Jonathan Altamirano & Yvonne Maldonado
Department of Pediatrics, Stanford University School of Medicine, Stanford, CA, USA
ChunHong Huang, Yuan J. Carrington, Frank Zhou & Yvonne Maldonado
Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine, Stanford, CA, 94305, USA
Jason R. Andrews
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Y.M., J.R.A., and K.S.W. designed the study. J.A., C.H., Y.C., and F.Z. conducted DNA extractions, sample preparation, Sanger sequencing, and data management. K.W. performed data analysis and drafted the paper. Y.M. and J.R.A. reviewed the paper and supervised the project.
Correspondence to Katharine S. Walter.
The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Walter, K.S., Altamirano, J., Huang, C. et al. Rapid emergence and transmission of virulence-associated mutations in the oral poliovirus vaccine following vaccination campaigns. npj Vaccines 8, 137 (2023). https://doi.org/10.1038/s41541-023-00740-9
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41541-023-00740-9
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.