Skip to main content

Sabin polio virus protein 1 (VP1) evolution in patients with acute flaccid paralysis from 2010 to 2016 in Uganda

This article has been updated


Acute flaccid paralysis (AFP) is a rare side effect of the oral polio vaccine but can be associated with outbreaks and permanent disability in patients harboring circulating vaccine-derived polioviruses (cVDPVs). With the advancement of polio abolition in a glimpse, cVDPVs are causing outbreaks and slowing the polio eradication process. The polio virus protein 1 (VP1) contains the binding site that is key for virus transmission. Understanding the evolution of VP1 among AFP patients could yield more insight into the early events of cVDPVs. Polioviruses were identified from stool specimens of AFP patients using cell culture; and confirmed by the real time RT PCR intra-typic differentiation and vaccine-derived poliovirus assays. Seventy-nine (79) Sabin-like poliovirus 1 (SL1) and 86 Sabin-like poliovirus 3 (SL3) were sequenced. The VP1 amino acid substitutions T106A in Sabin poliovirus 1 and A54V in Sabin poliovirus 3 were common among the AFP patients as has been found in previous studies. Other substitutions that were associated with AFP were: T290A and A54T in SL1 and SL3 respectively. Nucleotide mutations that were common among the AFP patients included T402C, C670A, and T816C in SL1, and G22A, C375Y, A472R, and A694T in SL3 polioviruses. Characterizing mutations that are associated with AFP could contribute to efforts pursued to mitigate the risk of vaccine-derived polioviruses and promote development of safer vaccines.


The poliovirus belongs to Enterovirus C species and has 3 serotypes; serotype 1, 2 and 3. Transmission occurs through the fecal–oral route. The virus multiplies in the gut and in 2% of infected persons, it enters and replicates in the nervous system resulting in muscle paralysis that manifests as acute flaccid paralysis [1]. Acute flaccid paralysis (AFP) is clinically characterized by sudden onset of weakness and fever. The patient presents with asymmetrical paralysis, reduced muscle tone however, sensation remains intact [2].

The genome of the poliovirus is 7.5 kb long with one open reading frame which is translated into four structural proteins (VP1–VP4) and non-structural viral proteins (2A–2C, 3A–3D). The VP1 protein plays a role in virus binding and has been used to type and track poliovirus transmission [3]. Furthermore, it has been shown to predict phylogenetic inference for the complete genome thus the small protein has been used to estimate the evolutionary relationship for polioviruses [4]. Virus RNA packaging and release has been mapped to the amino end of the VP1 region [5]. Lastly, the VP1 forms part of the major neutralization site that seldom contribute to virus escape [6, 7].

The evolution of VP1 of the poliovirus occurs because of the lack of 3′–5′ exonuclease proof reading mechanism for the virus coded RNA dependent RNA polymerase [8, 9]. During virus replication, the RNA polymerase introduces mutations in the virus genome, some of which occur in the VP1 region and are occasionally associated with reversion to a neuro pathogenic trait. The rate of misincorporation of bases ranges 10−5–10−3 per replication cycle [10, 11]. VP1 substitutions accumulate at a higher rate in Sabin-like polioviruses than the wild polioviruses [12] and occur within the first or second months following the administration of oral polio vaccine (OPV) [13, 14]

Numerous mutations have been identified among patients of AFP. These are usually silent synonymous mutations, and a selective constraint is applied on the non-synonymous ones [15, 16]. Yan et al. [17] identified one nucleotide mutation C161U that resulted in substitution A54V in 3 out of 4 Sabin-like type 3 polioviruses collected from a cluster of high-risk AFP patients. In another study, amino acid changes in VP1 of Sabin polioviruses 1 namely H149Y, T106A or I90L have been reported among AFP patient with residual paralysis [18].

De-attenuating mutations like T6I in Sabin poliovirus 3 [19] and single mutations can also alter pathogenesis; a single mutation (I143T) in Sabin poliovirus 2 VP1 region has been associated with neuro-pathogenesis [20]. Numerous mutations among AFP patients result in vaccine-derived polioviruses (VDPV) and are classified as VDPV type 1, 2, and 3 [21].

The aim of this study was to characterize VP1 mutations of Sabin poliovirus serotype 1 and 3 isolates that are associated with acute flaccid paralysis in Uganda. This study could provide additional insight into the early events of vaccine-derived polioviruses and could open avenues for intervention.


It is a retrospective cross-sectional study designed to investigate mutations that are common in patients of acute flaccid paralysis. The permission to undertake this study was obtained from the Ministry of Health, Uganda, the owner of the specimens and the Uganda Virus Research Institute, Research and Ethics Committee approved the study.

Archived poliovirus isolates from AFP patients aged 1 day–180 months were identified through the routine national AFP surveillance system. The Sabin poliovirus 1 and 3 that were investigated were collected between 2010 and 2016 before the switch from trivalent OPV (OPV 1, 2, and 3) to bivalent OPV (OPV1 and 3). The incidence of AFP patients during this period ranged from 2.3 to 8.5% (data not published). The percentage described represents the number of AFP patients who were shedding the poliovirus compared to the total number of AFP patients who were reported annually. The highest rate was observed in 2016 and corresponded to the period of enhanced polio campaigns prior to the switch.

Virus isolation

Poliovirus isolation was primarily characterized through cell culture using the rhabdomyosarcoma (RD) and mouse-derived L cell lines (L20B) [22]. The poliovirus type was then confirmed using real time Reverse Transcriptase Intra-Typic Differentiation and Vaccine-Derived Poliovirus assays (rRTPCR ITD and VDPV assays) [23].

RNA extraction

RNA extraction was performed from archived virus isolates using QIAamp Viral RNA mini kit (Qiagen) according to the manufacturer’s protocol.

Sabin RT-PCR procedure

The procedure was adopted and optimized using the standard protocol for VP1 gene characterization (CDC Laboratory). The Qiagen One step RT-PCR was used according to the manufacturer’s protocol. Briefly the RT-PCR amplifying a 1.1 kb encompassing the VP1 region was performed in a 50 µl reaction that consisted of: 30.5 µl RNase free water; 10.0 µl 5 × Qiagen One Step RT-PCR buffer; 2.0 µl dNTPs (10 mM each); 2.0 µl 5 × Qiagen One Step RT-PCR enzyme mix; 1.0 µl Y7R primer (40 pmol/µl), 5′GGTTTTGTGTCAGCITGYAAYGA-3′; 1.0 µl Q8 primer (10 pmol/µl), 5′AAGAGGTCTCTRTTCCACAT-3′; 0.5 µl RNase Inhibitor (40 U/µl), and 3.0 µl of the template (RNA extract). The thermo-profile condition used for the RT-PCR reactions was as follows: reverse transcription 50 °C for 30 min, RT enzyme inactivation 95 °C for 15 min; amplification: 94 °C for 30 s, 45 °C for 30 s, 72 °C for 1 min × 35 cycles; final extension 72 °C for 10 min; and finally left to hold at 4 °C.

Gel electrophoresis

One gram (1.00gm) of agarose powder was weighed and added to 100mls of 0.5% TBE buffer which was heated and cooled. Four micro-litres (4.0 µl) of ethidium bromide were added to the gel and a comb positioned in the casting tray to create wells for the samples. After the gel had solidified, the PCR products were then loaded onto the wells. A 100 bp molecular DNA ladder marker and a positive control specimen were added for amplicon sizing to ensure that the correct virus fragment was amplified. The PCR products were run at 130 V for 45 min in the electrophoresis tank. After gel electrophoresis, the 1.1 kb DNA products were viewed under ultraviolet light and the images were taken for documentation. The DNA products were stored at 20 °C till further testing.

Cleaning DNA

The Invitrogen Charge Switch PCR Clean-up Kit was used according to the manufacturer’s protocol. Briefly, binding of cDNA was achieved when 50 µl of the purification buffer and 50 µl PCR product together with 10 µl of the Charge Switch magnetic beads were mixed in a microcentrifuge tube and incubated at room temperature for one minute. The mixture was then placed on a Magna rack and the supernatant was removed and discarded. The complementary DNA (cDNA) bound to the beads was washed twice using 150 µl wash buffer for each sample. The sample was again removed from the rack and 50 µl of the elution buffer was added. The beads and buffer were mixed by gentle pipetting. The mixture was placed back on the Magna Rack and incubated for 1 min. The supernatant which contained the purified cDNA was collected, quantified using nanodrop, and stored at − 20 °C.

VP1 region sequencing

The Big Dye Terminator v3.1 cycle sequencing kit was used according to the manufacturer’s protocol. Briefly, a 10.0 µl sequencing reaction was set up consisting of: 5.0 µl RNase free H20, 1.0 µl sequencing buffer × 5, 1.0 µl primer (3.2 pmol) (primarily 246S/249S/249A/Q8 for SL1 or 248S/251S/ 261A/Q8 for SL3), 2.0 µl Big dye, and 1.0 µl [cDNA] (20-40 ng). The sequencing positive control reaction consisted of: 4.5 µl dH20, 1.0 µl Seq. buffer, 1.25 µl primer M13, 2.0 µl Big dye, and 1.25 µl pGEM DNA. The following thermo-profile condition was used for the sequencing reactions: amplification 95 °C for 15 s, 42 °C for 15 s, 60 °C for 4 min for 25 cycles and left to hold at 4 °C. The DNA was then stored at − 20 °C.

Cleaning sequencing products

Sequencing product cleanup was performed using the Agencourt ‘Clean SEQ according to the manufacturer’s protocol with a modification. Briefly, 10 µl of the magnetic beads were added to each specimen. Forty-two (42) µl of 85% ethanol was added to the specimen and mixed to form a homogeneous mixture. To enhance DNA binding the plate was left for 5 min before placing on the magnet (modification). The plate was then placed on the Magnetic Agencourt block for 5 min. The supernatant was removed, and each well was washed twice using 85% ethanol while the plate was on the magnet. The plate was left to dry for 10 min at room temperature. Thereafter, 42 µl of 0.1 mM EDTA (elution buffer) was added, and the plate was incubated at room temperature for 5 min to elute the DNA. Thirty (30) µl of the clear solution that contained the eluted DNA was carefully removed with a pipette ensuring that there was no bead carry over. The DNA product was transferred into the new wells ready for loading onto the detector. The plate was sealed with a septae and stored at − 20 °C until it was run on the 3500 ABI Genetic analyzer.

Data analysis

The full VP1 contigs were assembled using Sequencher software 4.10.1. Alignment of the consensus sequences and phylogenetic inference were performed using Mega 7.0 sequence analysis software [24]. Conservative and non-conservative mutations in VP1 were then identified in each independent alignment. The identification of nucleotide and amino acid mutations was performed by tally when changes in nucleotide and amino acid residues were compared to the reference Sabin poliovirus 1 and 3. Mutations were captured as point nucleotide and amino acid differences in the study VP1 sequence alignments.

Mutations from previous circulating VDPV outbreaks

Poliovirus VP1 sequences from outbreaks linked to circulating VDPVs [25] were downloaded from the National Center for Biotechnology Information (NCBI) and analyzed. Fifty-nine (59) and 24 sequences of SL1 and SL3 were available. The study VP1 mutation sites described were compared with those identified in viruses from the reported outbreaks.

AFP surveillance system

The study specimens investigated were kindly obtained from the routine national AFP surveillance system and the data reported was obtained from the “poliomyelitis/ AFP investigation data collection tool for acute illness”. The clinical presentation, immunization history, together with the preliminary specimen tracking were obtained from the collection tool. Sabin follow-up was performed for the AFP cases in whom the Sabin-like viruses were isolated.


The VP1 sequencing protocol was primarily optimized to detect polio sequences. A total of 79 SL1 and 86 SL3 were sequenced from the archived poliovirus isolates. In total, there were 65 SL1s, 81 SL3s and 16 mixtures of SL1 and SL3 polioviruses. Two (2) SL1 and 7 SL3 sequences failed on sequencing and 4 SL3 were missing. Twenty-nine and 50 mutations were observed in Sabin poliovirus 1 and 3 respectively. The nucleotide and amino acid mutations identified are respectively shown below in Tables 1 and 2.

Table 1 Nucleotide and corresponding amino acid changes among Sabin polioviruses type 1 isolates
Table 2 Nucleotide and corresponding amino acid change among Sabin polioviruses type 3 isolates

In this study, conservative and non-conservative changes have been described. A nucleotide substitution that changes the corresponding amino acid property in the protein is denoted non- conservative substitution (‘NC’), whereas a nucleotide substitution that does not change the amino acid property in the protein is denoted conservative substitution (‘C’) (see Tables 1 and 2 below).

There were more VP1 mutations for Sabin poliovirus 3 than Sabin poliovirus 1 (SL1). The SL1 viruses contained 29 mutations with 4 non-conservative substitution compared to SL3 that contained 50 mutations with 3 non-conservative substitutions.

A complementary presentation of VP1 nucleotide mutations for SL1 and SL3 sequences is shown in Table 3 placed after the text.

Table 3 Complementary presentation of VP1 nucleotide mutations for Sabin-like poliovirus 1 and 3

A total of 24 SL1 and 33 SL3 VP1 mutants have been identified in Table 3. The frequency of the nucleotide mutations for SL1 and SL3 can also be deduced from Table 3. The common nucleotide mutations sites (appearing more than once) in SLI were: A316R/G, T402C, C670A, T816C, and G870R, and those in SL3 were : G22A/R, G160A, C161T, C375Y, A472R, and A694T/W. The estimated VP1 mutation rate for SL1 and SL3 viruses in this study has been 4.1 × 10−4 (29 × 100/79 × 906), and 5.2 × 10−4 (50 × 100/86 × 900) respectively.

A map of Uganda showing the districts of origin for the VP1 non-conservative mutations among AFP patients is shown in Fig. 1 below.

Fig. 1
figure 1

Map of Uganda showing the districts of origin of the VP1 non-conservative mutations among AFP patients

The non-conservative VP1 mutants were identified in districts across the country at specified times. To understand VP1 sequence variation better, phylogeny inference to reference Sabin polioviruses 1 and 3 was performed (data not shown). There was minimal nucleotide divergence among the Sabin 1 and 3 sequences however, isolate UGA-16-3603/4 was divergent from the reference Sabin poliovirus 3. Study VP1 mutation sites described were compared with those identified in viruses from the reported cVDPV outbreaks. Mutation T106A in SL1 viruses; and A54V and M105T in SL3 viruses existed in the VDPVs as well [25, 26].

Table 4 displaying the relationship of the OPV dose to the mutations identified in SL1 and SL3 is shown below.

Table 4 Oral polio vaccine doses corresponding to VP1 mutations with complete vaccination, clinical and specimen documentation

Sixty-seven percent (67%) of AFP patients selected in Table 4 above had received adequate OPV doses. Acute flaccid paralysis was linked to multiple VP1 mutations in SL3 compared to SL1 isolates. The children who have been investigated above recovered with no residual paralysis.

The relationship between the ‘time interval from vaccination to specimen collection’, and  ‘VP1 mutations’

The Mann–Whitney test was used to establish the relationship between the ‘interval from vaccination to sample collection,’ and ‘VP1 mutations’. The assessment was constrained by the missing data; the vaccination date was missing for many children. Two samples are expected per child; to avoid the bias of repeated measure, one sample was considered per child and 42% of the samples were analyzed. There was no difference in the ‘vaccination interval’ for children who had a mutation (s) and those who did not have (p-value 0.867).


There was no vaccine-derived poliovirus nor wild type poliovirus detected among specimens that were collected from 2010 to 2016. No poliovirus type 2 was identified because virus sequencing was performed after the switch from trivalent to bivalent OPV when work on type 2 polioviruses was ceased under the WHO Global Action Plan (GAP III) containment. All polio virus type 2 isolates were destroyed and only SL1 and SL3 isolates were investigated. There were 29 mutations versus 4 non-conservative mutations for Sabin poliovirus 1 and 50 mutations versus 3 non-conservative mutations for Sabin poliovirus 3 which implicates a robust negative selection for Sabin poliovirus 3. Several amino acid substitutions were common among the SL1 and SL3 isolates and these included T106A, T290A, A54V and A54T.

For SL1, T106A was identified in 3 isolates from children residing in different districts namely Soroti, Kyankwanzi and Bududa in 2011, 2014, and 2015 respectively. Kyankwazi and Bududa are 530 km apart and had an 8 months’ time interval of virus detection. The distance and time interval of processing the samples did not predict a possible epidemiological linkage. Mutation T290A in SL1 was also reported in three AFP patients from Tororo district in 2012; and Kween and Bushenyi districts in 2015. Bushenyi and Kween are 631 km apart and sample processing time interval for the two viruses was 1 month. There was a study limitation in that samples from other countries that were processed around the same time of identifying the mutation could not be included and sequenced due to ethical issues.

Another mutation A54V was identified in SL3 and was detected in 3 isolates from Soroti, Iganga, and Rubirizi districts. Isolates from Soroti and Iganga were identified in 2012 whereas the isolate from Rubirizi was identified in 2014. The specimens identified in the same year had a 5 months’ difference between the processing time. Again, the time interval was not agreeable to an epidemiological linkage. Another mutation was detected at the same position A54T. This mutation was detected in three AFP patients from Maracha, Oyam, and Mityana districts in 2011, 2014, and 2015 as well. Oyam district is 320 km distant from Mityana district and the 2 viruses from both districts were detected within a time interval of 8 months. VP1 amino acid substitution at position 106 for SL1, and 54 for the SL3 viruses existed in both AFP and reported VDPV sequences. Mutations at these sites might be linked to vaccine-derived polioviruses.

The virus sequence UGA-16-3603/4 was obtained from a 2 months’ old child who developed AFP one day after receiving OPV1. The detected isolate contained 4 mutations and was divergent from the SL3 VP1 reference sequence. The virus is unlikely to have evolved in the recipient within one day after receiving the vaccine. There are two possible scenarios; the mutant could have evolved from the first OPV dose or, it could have been a case of environmental spread. The child resided in Bududa district which had sub-optimal tOPV + 3 coverage during the study period [27]. This district is one of the hard-to-reach districts with natural calamities of landslides.

It would be interesting to know whether there are immune deficient children among the AFP children who were shedding the virus. It is not possible to identify such children in this study because specimens for AFP surveillance are obtained at two time points within an interval of 24–48 h.

Overall, the AFP patients from the different regions of Uganda did not progress to vaccine-derived polioviruses implying that the immunity of the population was adequate and able to interrupt transmission. The national OPV3+ coverage ranged from 80 to 90% during the study period [28]. Polio campaigns were conducted in the whole country in 2012, 2015, and 2016 and in addition, targeted campaigns were also performed in the high-risk regions for VDPV and wild poliovirus transmission.


Several VP1 mutations common among the AFP patients have been described and supported by this study. The data generated could support innovations to design safer attenuated vaccines.


In this study, virus isolates were used, and culture adapted viruses might pose a bias. It was not possible to adequately relate ‘OPV doses’ and ‘vaccination interval’ to ‘VP1 mutations’ because of a gap in documentation.

Availability of data and materials

All data generated during this study appears in the published article.

Change history

  • 27 August 2023

    Some values from Tables were missed to be in bold



Virus protein 1


Acute flaccid paralysis


Vaccine-derived poliovirus


Sabin-like poliovirus 1


Sabin-like poliovirus 3


Ribonucleic acid


Oral polio vaccine


Rhabdomyosarcoma cells


Mouse L cells specific for poliovirus


Polymerase chain reaction


Real-time reverse transcriptase PCR


Intra-typic differentiation


Centers for disease control and prevention




Deoxyribonucleic acid


Complementary DNA


Ethylenediamine tetra acetic acid


Conservative amino acid substitution


Non-conservative amino acid substitution


National center for biotechnology information


World Health Organization


Global action plan III


  1. Racaniello VR. One hundred years of poliovirus pathogenesis. Virology. 2006;344:9–16.

    Article  CAS  PubMed  Google Scholar 

  2. Nomoto A. No title eradication of poliomyelitis. Nat Immunol. 2002;3(3):205–8.

    Article  CAS  PubMed  Google Scholar 

  3. Kew OM, Mulders MN, Lipskaya GY, da Silva EE, Patlansch MA. Molecular epidemiology of polioviruses. Semin Virol. 1995;6(6):401–14.

    Article  Google Scholar 

  4. Liu Y, Ma T, Liu J, Zhao X, Cheng Z, Guo H, et al. Bioinformatics analysis and genetic diversity of the poliovirus. J Med Microbiol. 2014;63:1724–31.

    Article  PubMed  Google Scholar 

  5. Kirkegaard K. Mutations in VP1 of poliovirus specifically affect both encapsidation and release of viral RNA. J Virol. 1990;64(1):195–206.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Shaw J, Jorba J, Zhao K, Iber J, Chen Q, Adu F, et al. Dynamics of evolution of poliovirus neutralizing antigenic sites and other capsid functional domains during a large and prolonged outbreak. J Virol. 2018;92(9):10–1128.

    Article  Google Scholar 

  7. Chen Z, Fischer ER, Kouiavskai D, Hansen BT, Ludtke SJ, Bidzhieva B, et al. Cross-neutralizing human anti-poliovirus antibodies bind the recognition site for cellular receptor. Proc Natl Acad Sci U S A. 2013;110(50):20242–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Jorba J, Campagnoli R, De L, Kew O. Calibration of multiple poliovirus molecular clocks covering an extended evolutionary range. J Virol. 2008;82(9):4429–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Kew OM, Sutter RW, De Gourville EM, Dowdle WR, Pallansch MA. Vaccine-derived polioviruses and the endgame strategy for global polio eradication. Ann Rev Microbiol. 2005;59:587–635.

    Article  CAS  Google Scholar 

  10. Ward CD, Flanegan JB. Determination of the poliovirus RNA polymerase error frequency at eight sites in the viral genome. J Virol. 1992;66(6):3784–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. De la Torre JC, Giachetti C, Semler BL, Holland JJ. High frequency of single-base transitions and extreme frequency of precise multiple-base reversion mutations in poliovirus. Proc Natl Acad Sci U S A. 1992;89(7):2531–5.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Famulare M, Chang S, Iber J, Zhao K, Adeniji JA, Bukbuk D, et al. Sabin vaccine reversion in the field: a comprehensive analysis of sabin-like poliovirus isolates in Nigeria. J Virol. 2016;90(1):317–31.

    Article  CAS  PubMed  Google Scholar 

  13. Minor PD, John A, Ferguson M, Icenogle JP. Antigenic and molecular evolution of the vaccine strain of type 3 poliovirus during the period of excretion by a primary vaccinee. J Gen Virol. 1986;67(4):693–706.

    Article  CAS  PubMed  Google Scholar 

  14. Valesano AL, Taniuchi M, Fitzsimmons WJ, Islam MO, Ahmed T, Zaman K, et al. The early evolution of oral poliovirus vaccine is shaped by strong positive selection and tight transmission bottlenecks. Cell Host Microbe. 2021;29(1):32–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Jorgensen D, Pons-Salort M, Shaw AG, Grassly NC. The role of genetic sequencing and analysis in the polio eradication programme. Virus Evol. 2020;6(2):veaa040.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Wong W, Gauld J, Famulare M. From vaccine to pathogen: modelling sabin 2 vaccine virus reversion and evolutionary epidemiology. 2020

  17. Yan DM, Zhang Y, Wang DY. Analysis of 4 clustered high risk acute flaccid paralysis cases in Shanxi Province in 2006. Zhongguo ji = hua mian yi Chin J Vacc Immun. 2010;16(2):127–31.

    Google Scholar 

  18. Rahimi P, Tabatabaie H, Gouya MM, Zahraie M, Mahmudi M, Ziaie A, et al. Characterization of mutations in the VP(1) region of Sabin strain type 1 polioviruses isolated from vaccine-associated paralytic poliomyelitis cases in Iran. J Clin Virol. 2007;39(4):304–7.

    Article  CAS  PubMed  Google Scholar 

  19. Minor PD. The molecular biology of poliovaccines. J Gen Virol. 1992;73:3065–77.

    Article  CAS  PubMed  Google Scholar 

  20. De Melo Cassemiro KMS, Burlandy FM, Barbosa MRF, Chen Q, Jorba J, Hachich EM, et al. Molecular and phenotypic characterization of a highly evolved type 2 vaccine-derived poliovirus isolated from seawater in Brazil, 2014. PLoS ONE. 2016;11(3):e0152251.

    Article  Google Scholar 

  21. World Health Organization(WHO). Classification and reporting of vaccine-derived polioviruses (VDPV) GPEI guidelines. 2016;

  22. World Health Organization. Polio Laboratory Manual [Internet]. 4th ed. 2004.

  23. Kilpatrick DR, Ching K, Iber J, Chen Q, Yang SJ, De L, et al. Identification of vaccine-derived polioviruses using dual-stage real-time RT-PCR. J Virol Meth. 2014;197:25–8.

    Article  CAS  Google Scholar 

  24. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Kew O, Morris-Glasgow V, Landaverde M, Burns C, Shaw J, Garib Z, et al. Outbreak of poliomyelitis in hispaniola associated with circulating type 1 vaccine-derived poliovirus. Science (80-). 2002;296(5566):356–9.

    Article  CAS  Google Scholar 

  26. Sadeuh-Mba SA, Kavunga-Membo H, Joffret ML, Yogolelo R, Endegue-Zanga MC, Bessaud M, et al. Genetic landscape and macro-evolution of cocirculating coxsackieviruses A and vaccinederived polioviruses in the democratic Republic of Congo, 2008–2013. PLoS Negl Trop Dis. 2019;13(4):e0007335.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Unpublished data. Aggregated Country Data for AFP Surveillance, Uganda National Expanded Programme on Immunisation/Ministry of Health, Uganda.

  28. WHO and UNICEF. Uganda: WHO and UNICEF estimates of Immunization coverage: 2019 revision [Internet]. 2019 Accessed 15 Sept 2021

Download references


The primers used in this study were kindly donated by Cara Burns. She critically reviewed the manuscript and together with Qi Chen assisted in the optimization of the VP1 sequencing protocol. The WHO Country office, and Government of Uganda offered logistical support. The appreciation is extended further to the Director General Health Services, Ministry of Health, Dr. Henry G. Mwebesa, and Management, Uganda Virus Research Institute for the continued administrative support; and to Mr. Samuel Ofori Gyasi for his assistance with data presentation. We are also indebted to the parents and children who donated specimens and made it possible to perform this study.


The World Health Organization and Government of Uganda supported this work.

Author information

Authors and Affiliations



Conceptualization MBN; methodology MBN, HB, PN, PT, MB, RD, JPE, MT, PK, BB, and JB, formal analysis HB, MBN; writing original draft preparation MBN; writing—review and editing PT, CRB, MBN, BB, JB; supervision CRB, BB, JB; and administration BB, JB, CRB. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Mary Bridget Nanteza.

Ethics declarations

Ethics approval and consent to participate

Ethical approval and waiver of consent was obtained from the Uganda Virus Research Institute, Research and Ethics Committee, Reference: GC/127/878 and GC/127/002 respectively.

Consent for publication

The publication of this data has been authorized by the Ministry of Health, Uganda, Reference: ADM.105/261/80.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nanteza, M.B., Bakamutumaho, B., Tushabe, P. et al. Sabin polio virus protein 1 (VP1) evolution in patients with acute flaccid paralysis from 2010 to 2016 in Uganda. Virol J 20, 172 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: