Previous Article | Next Article ![]()
Journal of Virology, August 2003, p. 8366-8377, Vol. 77, No. 15
0022-538X/03/$08.00+0 DOI: 10.1128/JVI.77.15.8366-8377.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Division of Viral and Rickettsial Diseases, National Center for Infectious Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia 30333,1 Virology Department, The Egyptian Organization for Biological Products and Vaccine Production (VACSERA), Cairo, Egypt,2 National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands,3 Department of Virology II, National Institute of Infectious Diseases, Tokyo, Japan4
Received 13 March 2003/ Accepted 16 May 2003
|
|
|---|
|
|
|---|
Some OPV-derived isolates, however, are more divergent (
99% VP1 sequence identity to the corresponding Sabin strain), indicative of prolonged replication of the vaccine virus. These isolates are described as vaccine-derived polioviruses (VDPVs). Two categories of VDPV isolates have been identified. First to be recognized and extensively characterized were the immunodeficient VDPV (iVDPV) isolates, obtained from patients with defects in antibody production (generally common variable immunodeficiency or X-linked agammaglobulinemia) (4, 25, 60, 68, 79). Some iVDPV isolates are highly divergent (
90% VP1 sequence identity to the parental OPV strain), suggesting that the chronic poliovirus infections had persisted for 10 years or more (25, 68). So far, all reports of persistent iVDPV infections have been from countries with high or intermediate levels of development, where the rates of OPV coverage are high and where the survival times of immunodeficient patients may be extended by their access to appropriate clinical management (N. A. Halsey, personal communication). Chronic iVDPV infections are independent events (57, 60), and the isolates obtained from such infections trace separate pathways of divergence from the original OPV strains (25, 35; R. Park, unpublished results). Currently, there is no clear evidence of spread of iVDPVs from immunodeficient patients to the wider community (57, 63, 68).
The second, more recently recognized category of divergent OPV-derived isolates includes circulating VDPVs (cVDPVs). The immediate public health importance of cVDPVs was underscored by the occurrence of an outbreak of 21 confirmed polio cases (including 2 fatal cases) associated with type 1 cVDPV on the Caribbean island of Hispaniola in 2000 to 2001 (22). Person-to-person transmission of VDPVs was initially suspected when it was found that the first two outbreak isolates were 2 to 3% divergent in VP1 sequence from the parental Sabin 1 OPV strain but shared numerous nucleotide differences from Sabin 1. Phylogenetic analysis of VP1 sequences of the 31 outbreak isolates suggested that the outbreak originated from a single OPV dose given in 1998 to 1999 (22). A separate, apparently much more limited outbreak of polio in the Philippines in 2001 was also found to be associated with type 1 cVDPV (8, 64). A third outbreak, in Madagascar, involving type 2 cVDPV, is currently under investigation (9). The critical risk factors for all three outbreaks were declining rates of OPV coverage and the prior eradication of the corresponding serotype of indigenous wild poliovirus (22).
In this report, we describe the widespread circulation of type 2 VDPV in Egypt in the 1980s and early 1990s. Like the outbreaks in Hispaniola, the Philippines, and Madagascar, VDPV circulation occurred after the apparent elimination of the corresponding serotype of indigenous wild poliovirus. However, unlike the recent outbreaks, where VDPV circulation was estimated to have occurred for less than 3 years (8, 22), type 2 VDPV circulation probably continued in Egypt for approximately a decade. Circulation occurred at a time when the rates of OPV coverage in Egypt were insufficient to eliminate the indigenous wild poliovirus types 1 and 3 (70, 71). VDPV circulation ceased when OPV coverage rates increased. Here we describe key genetic and biological properties of the Egyptian type 2 cVDPV isolates and discuss the implications of our findings to the overall strategy for global polio eradication.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Type 2 poliovirus isolates
|
Nucleic acid sequencing. Conditions for reverse transcription-PCR amplification and cycle sequencing were as described previously (32), using the primers listed in Table 2. Sequencing was performed in both directions, and every nucleotide position was sequenced at least once from each strand. Terminal sequences were determined by using the 5' and 3' rapid amplification of cDNA ends system kits (Life Technologies, Gaithersburg, Md.) according to the manufacturer's instructions.
|
View this table: [in a new window] |
TABLE 2. PCR and sequencing primers
|
Trees of VP1 sequence relationships among Egyptian type 2 cVDPV isolates were constructed by using the maximum likelihood method implemented in the fastDNAml program (40). A heuristic search was performed with 1,000 bootstrap replicates, with the input order of the sequences randomized for each replicate. The Ts/Tv ratio was initially set to 1, and the base frequencies were determined from the input data to simulate the F81 proportional model (13). The candidate tree was further analyzed using the program PAML (77) to compare other models of evolution in likelihood ratio tests (19). The TN93 model of nucleotide substitution (62) gave a significantly higher likelihood score than the F81 model; more complex models did not significantly improve the likelihood score. Parameters in TN93 were set to a Ts/Tv ratio of 10, and the shape parameter of the gamma distribution was calculated by PAML as 0.5 with eight categories of discrete gamma rates. The tree with the highest log likelihood score was rooted to the VP1 sequence of Sabin 2 and displayed using the program TreeView (41).
Estimation of the time of the initiating OPV dose. The time of the initiating OPV dose was estimated from the rate of fixation into VP1 of synonymous nucleotide substitutions among the 28 cVDPV isolates for which the dates of sample collection were known (Minya/92 and Sharqiya/92 were accordingly excluded). The number of VP1 substitutions at synonymous sites (Ks) that accumulated from the Sabin 2 sequence were computed according to the method of Yang and Nielsen (78) implemented in the program YN00 included in the PAML package (77). This method corrects for the Ts/Tv rate bias, the codon frequency bias, and multiple substitutions at a site. Corrected Ks values relative to the Sabin 2 root sequence (for which Ks is set to zero) for each cVDPV isolate were plotted as a function of the date of sample collection. The rate of accumulation of synonymous substitutions was estimated by weighted linear regression (where the weight for each data point was proportional to the reciprocal of the error variance for the corrected Ks value) using statistical applications within the SAS system, version 8 (SAS Institute, Inc., Cary, N.C.). The date of the initiating OPV dose was estimated from the intercept on the abscissa at Ks = 0, and the 95% confidence interval (CI) for the estimated date was calculated following procedures described by Sokal and Rohlf (56).
Estimation of nucleotide diversity. Nucleotide diversity was calculated from VP1 sequence differences and corrected for superimposed substitutions by use of the K2P algorithm (Ts/Tv ratio set at 10) by using the MEGA program, version 2 (30).
Comparison of complete genomic sequences of with those of Sabin 2. Complete genomic sequences of the recombinant cVDPV isolates, Qalyubiya/88-1 and Beni Suef/93, were compared with those of the Sabin 2 master seed stock. Sequence relationships for each viral gene (i.e., the 5' untranslated region [5' UTR] and each interval encoding a major cleavage product from the poliovirus polyprotein) were summarized in neighbor-joining trees constructed using the parameters described above for diverse wild type 2 poliovirus genotypes. Insertion-deletion 5' UTR differences were treated as single nucleotide substitutions.
Numbering of nucleotide and amino acid positions. VP1 sequences of all type 2 isolates were colinear. The coding sequences of Qalyubiya/88-1 and Beni Suef/93 were colinear with those of Sabin 2, but their 5' UTRs were not. To facilitate comparisons, numbering of nucleotide positions of all three isolates followed that described for Sabin 2 (49), with serial letters assigned to insertions.
Neurovirulence testing in PVR-Tg21 mice. Neurovirulence tests on cVDPV isolates were carried out on PVR-Tg21 mice as previously described (29, 31). The mice were purchased from the Central Laboratories for Experimental Animals (Kanagawa, Japan). The type 2 reference strains were MEF-1/EGY42 (neurovirulent) and Sabin 2 (attenuated). Six mice (equal numbers of males and females) were inoculated (30 µl/mouse) intracerebrally for each virus dilution (in 10-fold increments; range: 1.5 to 7.5 log 50% tissue culture infective doses/mouse). Mice were examined daily for 14 days postinoculation, and the times of paralysis or death were recorded. The virus titer that induced paralysis or death in 50% of inoculated mice was calculated by the method of Kärber (20) and expressed as 50% tissue culture infective doses/mouse.
Single-step growth curves. The growth rates and plaque yields at 39.5°C of Qalyubiya/88-1 and Beni Suef/93 were compared with those of Sabin 2, Beni Suef/98 (OPV-like), and MEF-1/EGY42 in single-step growth experiments (input multiplicity, 10 PFU/cell) in S3 HeLa cell suspension cultures (human cervical carcinoma, ATCC CCL-2.2). Virus titers were determined in plaque assays on HeLa cell monolayer cultures (ATCC CCL-2) incubated at 37°C essentially as described previously (51).
Nucleotide sequence accession numbers. VP1 sequences of the type 2 polioviruses described in this study were submitted to the GenBank library under accession numbers AF551795 to AF551841 (Table 1). Complete genomic sequences of isolates Qalyubiya/88-1 and Beni Suef/93 were submitted under accession numbers AF448782 and AF448783.
|
|
|---|
From 1988 to 1993, type 2 poliovirus was frequently isolated from patients with clinically confirmed poliomyelitis in Egypt (Fig. 1; Table 1). Because many of the patients were unimmunized and had no history of recent exposure to OPV, it was assumed at the time that these cases signaled the continued circulation of the indigenous wild poliovirus type 2. This view was reinforced by the finding that 30 of the type 2 isolates had DR antigenic properties when tested in the standard ELISA intratypic differentiation test using cross-absorbed antisera (65). The 30 DR case isolates were selected for further characterization as part of a broader study monitoring the eradication of wild poliovirus in Egypt.
![]() View larger version (34K): [in a new window] |
FIG. 1. Geographic distribution of polio cases in Egypt associated with type 2 cVDPV from 1988 to 1993. Individual cVDPV cases, represented by filled circles, are mapped by governorate (2 of the 30 cVDPV isolates were not mapped because their source governorates were not recorded). cVDPV cases occurring in different years are differentiated by increasing thickness of circle rims from 1988 to 1993. OPV-like isolates (>99% VP1 sequence identity to Sabin 2) from AFP cases with onset from 1995 to 2000 are identified by open circles. Governorate abbreviations: AL, Alexandria; AW, Aswan; AY, Asyut; BH, Behaira; BS, Beni Suef; CR, Cairo; DQ, Daqahliya; DM, Damietta; FY, Fayyum; GB, Gharbiya; GZ, Giza; IS, Ismailiya; KS, Kafr El-Sheikh; MF, Menufiya; MT, Matrouh; MY, Minya; NS, North Sinai; NV, New Valley; PS, Port Said; QL, Qalyubiya; QN, Qina; RS, Red Sea; SH, Sohag; SQ, Sharqiya; SS, South Sinai; SZ, Suez.
|
33% of total VP1 nucleotides, primarily at positions of codon degeneracy) saturate within
30 years. However, synonymous transversions are fixed into variable poliovirus sites at an approximately 10-fold lower rate than synonymous transitions (J. Jorba, unpublished results). Therefore, closely related isolates can be better resolved from more divergent viruses by comparing only transversional differences at fourfold degenerate sites. By this method of comparison, Sabin 2, Qalyubiya/88-1, and Beni Suef/93 are indeed very closely related (0.2 to 0.4% transversional differences in VP1) and are only distantly related to any of the wild type 2 polioviruses (3.7 to 10.6% transversional differences in VP1).
![]() View larger version (17K): [in a new window] |
FIG. 2. Unrooted neighbor-joining tree summarizing VP1 (nucleotides 2482 to 3384) sequence relationships among the Sabin 2 OPV strain, the Egyptian cVDPV isolates Qalyubiya/88-1 and Beni Suef/93, three Egyptian wild type 2 isolates (MEF-1/EGY42, 0317/EGY55, and 0297/EGY79), and nine wild type 2 poliovirus isolates from cases occurring in Africa, Asia, Europe, and South America from 1977 to 1991.
|
![]() View larger version (31K): [in a new window] |
FIG. 3. Maximum likelihood tree of VP1 sequence relationships among the Sabin 2 OPV strain (root of tree) and 30 type 2 cVDPV isolates from Egypt (1988 to 1993). Brackets enclose separate genetic groups (I to III) and their component clusters (a to k).
|
A single node connected the sequences of the 30 cVDPV isolates to the root Sabin 2 sequence. This node represents the VP1 sequence of the inferred most recent common ancestor to all 30 cVDPV isolates. Although the bootstrap values (14) supporting the node shown in Fig. 3 were only moderate (
44%), replicate trees constructed from the same data set by using the neighbor-joining (53) and quartet puzzling (58) algorithms had similarly placed ancestral nodes. The location of the ancestral node on the tree and the short length of the branch joining it to the Sabin 2 root sequence suggest that the most recent ancestral infection common to all subsequent lineages occurred about a year after the initiating OPV dose (see next section) and well before 1988, the year of the cases yielding the earliest available cVDPV isolates (Giza/88, Qalyubiya/88-1, and Qalyubiya/88-2) (Fig. 3). Soon after the ancestral infection, the virus radiated along separate chains of transmission, observed as separate branches on the tree. We have no further information about the early events in the spread of cVDPVs in Egypt because OPV-derived isolates dating prior to 1988 are not available.
The isolates within group I had the most diverse VP1 sequences (0.071 nucleotide differences per site), with the eight isolates distributing into five clusters (a to e). Most group I sequences were connected to each other by long branches. Long branches indicate missing sequences, as many polio cases were likely missed and specimens were not available for virologic testing. At least half (four of eight) of the group I isolates (the source governorates for two of the isolates were not recorded) were from Giza. Isolates from Giza were found in all three groups and in seven different clusters (a, b, c, d, g, i, and k). The high sequence diversity of the Giza isolates (0.070 nucleotide differences per site) does not necessarily indicate that multiple cVDPV lineages originated from that governorate. It is more likely that different cVDPV lineages were introduced into Giza from other communities in Egypt, as Giza is host to migrants from other parts of the country. Only one group I isolate (Minya/92, cluster b) was from the southern governorates, where OPV coverage and AFP surveillance has been lower than in the rest of Egypt. It is possible that several of the Giza cVDPV isolates (such as those comprising clusters b, c, and d) were derived from viruses circulating in areas remote from the Poliomyelitis Institute in Cairo.
Group II isolates were more closely related to each other (0.046 nucleotide differences per site). Most of the isolates were from the adjoining governorates of Qalyubiya, Menufiya, Giza, and Cairo (Fig. 1 and 3). The sequence of Qalyubiya/88-1 is close to the ancestral node connecting all group II isolate sequences. Some of the genetic clusters correlated with the geographic source of the virus (Fig. 1 and 3). For example, all three isolates of cluster h were from Qalyubiya, and most (six of nine) isolates of cluster i were from Beni Suef. Cluster g represents a lineage that circulated over four governorates and then apparently faded out after 1990. Cluster h, apparently localized to Qalyubiya, was also not detected after 1990. Isolates of cluster i were the last type 2 cVDPVs to be found in Egypt (last isolate: Beni Suef/93, 9 February 1993). There is clear evidence of sustained localized cVDPV transmission during a single polio season (which in Egypt runs from about the beginning of April to the end of March of the following calendar year). For example, isolates of cluster h were obtained from polio patients in Qalyubiya from 26 June 1990 to 24 November 1990, and isolates of cluster i were obtained from polio patients in Beni Suef from 8 August 1992 to 9 February 1993. There is also evidence of continued transmission of cVDPV lineages through one or more winter low-transmission seasons for poliovirus, as is seen in clusters g (Qalyubiya/88-2 [29 October 1988] to Giza/90-2 [8 August 1990]) and i (Giza/91-1 [18 October 1991] to Beni Suef/93 [9 February 1993]).
The two isolates forming group III (cluster k), Qalyubiya/89 and Giza/90-1, represent a lineage that diverged from the others well before 1988 and may have established a pocket of localized circulation in 1989 to 1990 (Fig. 3).
Genetic distance from the Sabin 2 VP1 sequence correlated with time of specimen collection (Fig. 3). Sequences of the earliest (1988 to 1989) isolates were located at the tips of branches closest to the root (left side) of the tree, sequences of isolates from 1990 to 1991 were located centrally, and sequences of isolates from 1992 to 1993 were located at the tips of branches at the right side of the tree.
Estimated time of occurrence of initiating OPV dose. The approximate time of the initiating OPV infection was estimated from the VP1 sequence differences from Sabin 2 among the cVDPV isolates, the dates of specimen collection, and the rate of VP1 evolution (22). The rate of substitution into synonymous sites of VP1 was approximately linear from 1988 to 1993 (R2 = 0.670) (Fig. 4). When the regression line was extrapolated to zero VP1 substitutions, the intercept was August 1983 (95% CI: July 1979 to March 1986). The CI was wide because of the stochastic nature of base substitution, because the earliest available cVDPV isolates were obtained from specimens taken about 5 years after the initiating OPV dose and because of the bias in the sources of isolates available for analysis. The VP1 evolution rate estimated from the slope of the regression line was (2.47 ± 0.69) x 10-2 substitutions per synonymous site per year, approximately 80% of the rate observed for circulating wild type 1 polioviruses (16, 25, 32).
![]() View larger version (18K): [in a new window] |
FIG. 4. Estimate of the date of initiating OPV dose from the rate of accumulation of synonymous substitutions into VP1 among the 28 type 2 cVDPV isolates for which the dates of sample collection are known. Abscissa: date of sample collection for each isolate. Ordinate: Ks, the number of substitutions (Sabin 2 sequence set to zero substitutions) at synonymous sites in VP1. The evolution rate was estimated by weighted linear regression (R2 = 0.670) as described in the Materials and Methods. The 95% CI for the estimated date of the initiating OPV dose, July 1979 to March 1986, is bounded by parentheses along the abscissa.
|
Tripartite recombinant structure of Qalyubiya/88-1 and Beni Suef/93 genomes.
The complete genomic sequences of an early type 2 cVDPV isolate, Qalyubiya/88-1, and the last isolate, Beni Suef/93 (Table 1), were compared with that of the Sabin type 2 OPV strain (Fig. 5). Qalyubiya/88-1 was selected as the early cVDPV isolate for comparison because it appeared to be closely related to the ancestor of clusters f to j (which includes Beni Suef/93, cluster i) (Fig. 3). While the capsid sequences were homologous to those of Sabin 2 (97.4% [Qalyubiya/88-1] and 94.4% [Beni Suef/93] nucleotide sequence identity) and to each other (95.1% nucleotide sequence identity), the flanking sequences were dissimilar. The 5' UTR sequences of both Qalyubiya/88-1 and Beni Suef/93 were apparently derived from viruses other than Sabin 2 (having <85% nucleotide sequence identity to the 5' UTR of Sabin 2). The conserved 5' UTR sequences at positions 1 to 646 (maximum observed divergence among poliovirus isolates:
89%) of both cVDPV isolates were not closely related for that interval (93.1% nucleotide sequence identity). By contrast, the highly variable and rapidly evolving (H.-M. Liu, unpublished results) 5' UTR sequences at positions 647 to 747 (maximum observed divergence among poliovirus isolates:
46%) appeared to be relatively more closely related (90.1% nucleotide sequence identity). Similarly, the noncapsid region and 3' UTR sequences of the two cVDPV isolates differed from Sabin 2 (82.7% [Qalyubiya/88-1] and 83.2% [Beni Suef/93] nucleotide sequence identity) and from each other (86.4% nucleotide sequence identity) (Fig. 5). The sequence relationships among the genomes of Sabin 2, Qalyubiya/88-1, and Beni Suef/93 suggest that at least two rounds of recombination occurred in both the 5' UTR and noncapsid region during cVDPV evolution. The heterologous sequences were not derived from the Sabin 1 or Sabin 3 OPV strains or any other characterized poliovirus and were most likely derived from unidentified species C enteroviruses. A type 2 OPV-like isolate having a similar tripartite recombinant structure has been described previously (17).
![]() View larger version (36K): [in a new window] |
FIG. 5. (A) Nucleotide (upper bars) and amino acid (lower bars) substitutions into the genomes of cVDPV isolates Qalyubiya/88-1 (above) and Beni Suef/93 (below). The reference sequence is that of Sabin 2. The substitution maps are aligned with a schematic of the poliovirus genome; the single open reading frame is indicated by a rectangle, flanked by the 5' and 3' UTRs. Asterisks indicate attenuating mutations replaced by recombination (5' UTR) or nucleotide substitution (VP1). (B) Neighbor-joining trees summarizing sequence relationships in the 5' UTR and in each interval encoding a major cleavage product from the poliovirus polyprotein; all trees are drawn to the same scale. Isolate sequence abbreviations: S2, Sabin 2; QL, Qalyubia/88-1; BS, Beni Suef/93.
|
G revertants, restoring the consensus nucleotide at that position (47). Both Qalyubiya/88-1 and Beni Suef/93 had the consensus G at position 481, but unlike in most OPV-like isolates, the surrounding 5' UTR sequences were probably derived from other enteroviruses. Isolates Qalyubiya/88-1 and Beni Suef/93 shared the same VP1-143 Ile
Thr amino acid replacement (restoring the consensus residue for wild type 2 polioviruses) encoded by a U
C substitution at nucleotide position 2909. Neurovirulence of Qalyubiya/88-1 and Beni Suef/93 to PVR-Tg21 transgenic mice. All 30 of the cVDPVs were isolated from patients with paralytic polio (Table 1), demonstrating their neurovirulence for humans under natural conditions. Moreover, the sequence properties of isolates Qalyubiya/88-1 and Beni Suef/93 suggested that they would be neurovirulent when introduced into the central nervous systems of experimental animals (34, 48). When the neurovirulence of these isolates was measured quantitatively in PVR-Tg21 transgenic mice expressing the human receptor for poliovirus, they were found to be as neurovirulent as the prototype wild type 2 poliovirus strain, MEF-1/EGY42 (Table 3), isolated in 1942 from patients with paralytic polio.
|
View this table: [in a new window] |
TABLE 3. Neurovirulence of Qalyubiya/88-1 and Beni Suef/93 in PVR-Tg21 transgenic mice
|
The growth rates and virus yields at 39.5°C of Qalyubiya/88-1 and Beni Suef/93 were compared with those of Sabin 2, Beni Suef/98, and MEF-1/EGY42 in single-step growth experiments in S3 HeLa cells (Fig. 6). Beni Suef/98 is an OPV-like isolate, having vaccinelike antigenic properties, one nucleotide substitution in VP1 (a U
A substitution at nucleotide position 2909 encoding a VP1-143 Ile
Asn replacement), and the consensus G nucleotide at position 481. The wild isolate, MEF-1/EGY42, grows to high yields in cell culture, a property that favored its selection as the type 2 component of the inactivated poliovirus vaccine (IPV) (46). MEF-1/EGY42 grew more rapidly and produced
90-fold higher titers than Sabin 2 during infection of S3 HeLa cells incubated at 39.5°C (Fig. 6). Interestingly, the late cVDPV isolate, Beni Suef/93, had growth properties at 39.5°C similar to those of MEF-1/EGY42. By contrast, the OPV-like isolate, Beni Suef/98, grew to only slightly higher titers (
twofold) than Sabin 2. The early cVDPV isolate, Qalyubiya/88-1, grew more rapidly and to much higher titers than Sabin 2 and Beni Suef/98, but the virus yield at 8 h was about one-third of the yields for MEF-1/EGY42 and Beni Suef/93 (Fig. 6).
![]() View larger version (22K): [in a new window] |
FIG. 6. Single-step growth curves at 39.5°C in S3 HeLa cells of MEF-1/EGY42 (wild), Qalyubiya/88-1 (cVDPV), Beni Suef/93 (cVDPV), Beni Suef/98 (OPV-like), and Sabin 2.
|
|
|
|---|
10 years), the high nucleotide diversity of the cVDPV isolates (consistent with multiple chains of transmission), the generally very low paralytic attack rate for type 2 poliovirus infections (39), and the likelihood that the large majority of polio cases that had occurred during the period of cVDPV endemicity in Egypt had not been investigated. VDPV circulation in Egypt apparently ceased after 1993 as the rates of OPV coverage increased (70, 71). We do not know what levels of OPV coverage actually achieved cessation of VDPV circulation in Egypt because of the absence of complete national coverage data for the early 1990s. It is likely that the OPV coverage rates necessary to interrupt virus transmission would be comparatively high in Egypt because of the prevalence of conditions favorable to poliovirus circulation (15, 39, 42). Nonetheless, the threshold rates of coverage necessary to stop poliovirus transmission are lower for type 2 than for types 1 and 3, probably because seroconversion to type 2 antigen is more efficient (43). The lower OPV coverage threshold required to control type 2 poliovirus may explain why circulation of the indigenous type 2 wild poliovirus and the subsequent type 2 cVDPV ceased in Egypt several years before elimination of the indigenous wild poliovirus types 3 (in 2000) and 1 (circulation continued throughout 2002) (71). We also do not know when the last case associated with type 2 cVDPV occurred in the country because surveillance gaps persisted into the 1990s, especially in the southern governorates where the rates of OPV coverage were lowest. However, as the sensitivity of AFP surveillance improved, no evidence has been found for VDPV circulation after 1993.
The endemic circulation of VDPV in Egypt differs in several important respects from the recent cVDPV outbreaks in Hispaniola (Dominican Republic and Haiti) (22), the Philippines (8, 64), and Madagascar (9). The most significant difference was the much longer period of VDPV circulation in Egypt than in the three recent outbreaks (estimated at
3 years). Moreover, the outbreaks in Egypt and Madagascar were associated with type 2 cVDPV, whereas the outbreaks in Hispaniola and the Philippines were associated with type 1 cVDPV. Genetic evidence suggests that cVDPV spread along multiple independent chains of transmission in both Egypt and Haiti (22). Two separate lineages of type 2 cVDPV have been found in Madagascar, with the first lineage arising about a year before the second (9). By contrast, the more localized outbreak in the Dominican Republic appears to have resulted from the point source introduction of cVDPV from neighboring Haiti (22). cVDPV spread appears to have been even more limited in the Philippines, as the genetic properties of the outbreak isolates suggest that transmission occurred along a single, minimally branched chain (8, 64).
The early virologic events leading to cVDPV outbreaks have not been observed directly. Such events would be difficult to observe even under ideal conditions because of the inherently low paralytic attack rates of poliovirus infections (39). However, the prospects for early detection of cVDPV outbreaks are further reduced in the areas of highest risk, because gaps in surveillance usually accompany gaps in OPV coverage. Nonetheless, a reasonable model for the genesis of cVDPV outbreaks may be proposed from the known properties of the Sabin OPV strains. Because the base substitutions contributing to the attenuated phenotype of the Sabin strains (48) readily revert during replication in the human intestine (34, 36), the virus excreted by healthy vaccine recipients is frequently less attenuated than the original OPV strains. Reversion is very rapid for the Sabin 2 strain because attenuation is primarily determined by only two highly unstable substitutions (34, 36, 48, 80). The intense selection against the attenuating mutations suggests that the revertants replicate more efficiently in human intestine (36). This may be expressed as an overall increase in the titers of excreted virus, an extended period of virus excretion, and possibly an increased potential for person-to-person transmission. Introduction of OPV at low rates of coverage could trigger the transmission of excreted vaccine virus. The extent of spread would depend on several factors, including the number and density of nonimmune susceptible persons in the population. It therefore appears likely that Sabin 2 revertants with increased potentials for neurovirulence and transmissibility are constantly selected in all communities where OPV is used, but their spread is normally restricted by high population immunity.
The observation that all cVDPV isolates were derived from a single initiating OPV dose does not necessarily indicate that only one cVDPV lineage emerged in Egypt during the 1980s. The conditions leading to the observed lineages may have also favored the independent generation of other cVDPV lineages, as has recently been found in Madagascar (9). However, the other lineages might have escaped detection in several different ways, namely, (i) by circulating in communities outside of those sampled, (ii) by spontaneously fading out before 1988, or (iii) by being displaced before 1988 by the observed lineages.
All outbreak-associated cVDPV isolates described thus far have been recombinants with other species C enteroviruses (8, 22, 64). However, this observation does not necessarily imply that recombination plays an obligatory mechanistic role in the phenotypic reversion of OPV, because the main determinants of attenuation of all three Sabin strains map to 5' UTR and capsid region sites (21, 34, 48, 66), and most of the observed recombination sites map to the noncapsid region (9, 22, 64; J. Iber, unpublished results). Although we observed that the interval containing the critical attenuating sites in the 5' UTR had been exchanged out by recombination in the Egyptian cVDPV isolates, it is likely that back mutation preceded recombination, because selection of 5' UTR revertants is very rapid during replication of OPV in the human intestine (33, 36, 80). Moreover, the successive recombination events involving noncapsid region (22) and 5' UTR sequences may not confer any progressive increase in replicative fitness. Rather, poliovirus recombination with other species C enteroviruses may at times occur during mixed infection, with the frequency of recombination being a function of the enterovirus carriage rate and the total number of mixed infections. Thus, if a vaccine-related isolate has significant divergence in its capsid nucleotide sequences (>1% from the parental OPV strains) and has noncapsid or 5' UTR sequences of non-Sabin strain origin, it should be considered a possible cVDPV isolate and the associated case should be investigated further.
The episodes of VDPV circulation in Egypt and other developing countries challenge the assumption that poliovirus endemicity can be restored only by the reintroduction of wild poliovirus and underscore the urgency of reaching the goal of global polio eradication as quickly as possible. The first priority is to eliminate the remaining reservoirs of wild poliovirus endemicity in South Asia and Sub-Saharan Africa (24, 63, 72). At the same time, it is essential to maintain high levels of poliovaccine coverage in all countries, both to prevent the spread of imported wild polioviruses and to suppress the emergence of cVDPVs. Areas at highest risk for emergence of cVDPVs are those where poliovaccine coverage rates have declined, the competing wild polioviruses have been eliminated, and epidemiologic conditions had previously favored wild poliovirus transmission. In response to the recent cVDPV outbreaks, the World Health Organization (WHO) has recommended reinstatement of mass immunization campaigns to close the immunity gap in areas where the rates of routine OPV coverage have been insufficient to prevent a growing susceptibility to polio (63). The optimal frequency of the mass campaigns follows from the rate of accumulation of susceptible persons in the population and the basic reproduction number for poliovirus in the highest risk populations in each area (15, 42). No additional measures have been recommended for the remaining countries of polio endemicity, because activities currently planned for elimination of the last pockets of types 1 and 3 poliovirus circulation will also effectively prevent dissemination of cVDPVs.
The occurrence of cVDPV cases also highlights the need to maintain sensitive poliovirus surveillance into the foreseeable future. Intensive screening of recent poliovirus isolates for cVDPVs has been initiated by laboratories within the WHO Global Polio Laboratory Network (5). Beginning in 2001, all vaccine-related poliovirus isolates from AFP cases have been screened for evidence of prolonged replication or circulation (6, 69). Poliovirus isolates are identified according to their genetic properties by probe hybridization (11), diagnostic PCR assays (26, 75), or PCR-restriction fragment polymorphism analysis (3). All isolates are also tested for antigenic change by using specific cross-absorbed sera in an ELISA format (65) or panels of monoclonal antibodies in neutralization tests (65). Alternatively, isolates have been screened for recombinant noncapsid sequences using PCR primers targeting P2 and P3 region sequences characteristic for each Sabin strain (D. R. Kilpatrick, unpublished results). Any isolate having non-vaccinelike, DR, or NR antigenic properties or having a recombinant genome is characterized further by VP1 sequencing (6). The combined strategy of molecular and antigenic characterization of poliovirus isolates led to the recognition of the cVDPV outbreak in the Philippines (8, 64) and to the identification of several candidate iVDPV isolates from developing countries (6).
Several observations suggest that the risk for emergence of cVDPVs may be highest for poliovirus type 2 (15). The type 2 OPV strain appears to spread more readily to unimmunized people than the other Sabin strains, as shown by the more frequent association of contact cases of vaccine-associated paralytic poliomyelitis with poliovirus type 2 (57) and by the much higher seroprevalence to poliovirus type 2 (relative to types 1 and 3) found among unvaccinated individuals in the United States (10, 15, 59) and Europe (2). Type 2 VDPVs resembling the cVDPVs from Egypt and Madagascar had been isolated from polio cases in other developing countries during times when OPV coverage rates were low (R. Park, unpublished results). However, because of the low paralytic attack rates for type 2 infections, circulation of type 2 VDPVs is the most difficult to detect by AFP surveillance. Consequently, early detection of any future type 2 cVDPV outbreaks will require maintenance, and possible augmentation in some countries (54), of the present very high global standard for AFP and poliovirus surveillance (63, 72).
The recognition that cVDPVs can cause polio outbreaks and potentially reestablish polio endemicity to polio-free areas has prompted a reassessment of global strategies for maintaining polio-free status after wild poliovirus circulation has ceased (12, 15, 38, 63, 68). The number of viable options for the endgame strategy using the existing poliovirus vaccines now appears to be quite limited (12, 38, 63). OPV use should continue in countries of polio endemicity and in developing countries where polio is not endemic but which are at high risk for outbreaks caused by imported wild polioviruses or emerging cVDPVs. Transition from OPV to IPV is likely to continue in polio-free countries, especially in developed countries in temperate zones where IPV efficacy is highest (67) and where high rates of IPV coverage can be maintained through routine immunization (63). Soon after certification of global eradication, countries still using OPV may discontinue its use immediately after completion of the last round of synchronous mass OPV campaigns (12, 63). In addition, IPV use might be implemented in all countries for an interim period until evidence is compelling that the risk of recurrence of poliovirus circulation is very low (12, 63). Whatever the strategy, implementation will be most challenging in low-income, tropical developing countries where the risks of cVDPV emergence are highest. In the meantime, universally high rates of polio vaccine coverage must be maintained at least until implementation of the WHO endgame strategy (63).
Hiroyuki Shimizu, Tetsuo Yoneyama, and Tatsuo Miyamura were supported by a grant in aid for Promotion of Polio Eradication from the Ministry of Health, Welfare and Labor, Japan.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»