Previous Article | Next Article ![]()
Journal of Virology, May 2004, p. 4876-4883, Vol. 78, No. 9
0022-538X/04/$08.00+0 DOI: 10.1128/JVI.78.9.4876-4883.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Microbiology, Enterovirus Laboratory, National Public Health Institute (KTL), 00300 Helsinki, Finland,1 Central Laboratory of Virology, Health Protection Inspectorate, 11315 Tallinn, Estonia2
Received 24 September 2003/ Accepted 2 January 2004
|
|
|---|
|
|
|---|
The evolution rate of polioviruses is very high, partly due to the high error frequency in RNA synthesis: roughly 104 per base pair per replication cycle. For wild polioviruses circulating in human populations, nucleotide substitutions accumulate at a rate of approximately 1% per year and consist primarily of changes at synonymous codon positions thus not leading to amino acid substitutions (22). The genetic diversity of poliovirus strains is exploited in molecular epidemiology, a key component of poliovirus surveillance (33), currently based on sequence analysis of the VP1 coding region of the genome. The replication of the attenuated poliovirus strains (Sabin vaccine strains) in the human gut also leads to genetic variation and may result in reversions of the vaccine strains to pathogenic phenotypes reminiscent of wild polioviruses.
In addition to the high mutation frequency, the divergence of poliovirus strains is increased by recombination. Intertypic recombination is a frequent phenomenon in poliovirus vaccinees, and strains with a recombinant genome have been isolated from both healthy vaccinees and from patients with vaccine-associated poliomyelitis (5, 15-17, 24). In natural intertypic recombinant poliovirus strains, the recombination junctions are usually located in the genomic region encoding the nonstructural proteins. However, in some strains, the recombination site has been shown to reside in the capsid protein VP1 coding region (2, 27).
OPV recipients are known to excrete OPV-derived polioviruses for various periods but usually not more than for a few months. During extended replication in immunodeficient individuals, the virus may accumulate point mutations and the modified virus, referred to as immunodeficiency-associated vaccine-derived poliovirus (iVDPV), has occasionally caused paralytic disease in the vaccinee (4, 18, 23, 26). VDPVs may also be generated through circulation of OPV-derived virus strains in human populations with deficient herd immunity (circulating VDPV [cVDPV]). Some VDPVs characterized in the literature cannot be classified in these two main categories, and their origin remains open to question (10, 11, 37). Recently, cVDPVs have been associated with four outbreaks of poliomyelitis. Type 2 VDPV circulated in Egypt for over 10 years (1983 to 1993) and was isolated from 30 patients (7, 41). An outbreak of poliomyelitis on the island of Hispaniola was associated with type 1 cVDPV (21). In the Philippines, type 1 cVDPV was involved in three poliomyelitis cases in 2001 (6, 8), and in Madagascar in 2002, type 2 cVDPV was the causative agent of four paralytic cases (9, 34). In each of these outbreaks, the cVDPV strains were recombinants, originating from either a type 1 or 2 OPV strain and in which most of the genomic region encoding the nonstructural proteins was suggested to be derived from an unknown non-polio enterovirus.
We report here the complete genomic sequence of a highly divergent variant of the Sabin strain of type 3 poliovirus (PV3) isolated from the main sewage channel of Tallinn, Estonia, in October 2002. Estonia has been free of poliomyelitis cases since early 1960s (20). The current virus was isolated during routine screening for wild-type poliovirus circulation. It was not linked to any epidemics, and the origin of the virus remains unresolved.
|
|
|---|
Virus isolation. Wastewater samples for environmental poliovirus surveillance were regularly collected in Tallinn by the "trap principle" essentially as described previously (3). Tightly bound multilayer pads of cotton gauze (about 10 by 10 by 2 cm) were fixed with a string to the sewer wall and immersed into the sewage stream in the inlet sewer entering the main sewage treatment plant or similar places elsewhere in the sewage network. After 6 days of exposure, the pads were removed and transported to the laboratory in a container with 10 ml of cold Hanks' balanced salt solution (pH 7.5). The container was shaken vigorously for 5 min, followed by squeezing the pad several times. All liquid was finally squeezed from the pad and poured into a sterile centrifuge tube. pH was adjusted to 8.0 to 8.4 with 1 M HCl or NaOH, and the sample was clarified by low-speed centrifugation for 20 min. Supernatant was poured into a new sterile tube and extracted with 0.5 ml of chloroform per 10 ml of supernatant. The chloroform-extracted elute was examined for the presence of poliovirus as recommended for stool extracts in the WHO Polio Laboratory Manual (40).
During the intensified sampling period, two 1-liter samples of raw sewage were drawn from the inlet of the sewage treatment plant and transported at 4°C within a few hours to the National Public Health Institute (KTL) in Helsinki, Finland. At KTL, the two-phase separation method was used essentially as described previously (31), except that the entire sample was concentrated for analysis in cell culture.
Serotyping and ITD. The serotype of the isolate was confirmed at KTL by neutralization with polyclonal antisera obtained from the National Institute of Public Health and the Environment (RIVM), Bilthoven, The Netherlands, according to the protocol recommended by the WHO (40).
Antigenic intratypic differentiation (ITD) was performed with two validated assays. First, an enzyme immuno assay (EIA) was carried out with cross-absorbed antisera for PV3/Sabin and PV3/non-Sabin-like (RIVM) (39). Second, a neutralization assay was carried out with a mixture of Sabin-specific monoclonal antibodies (MAbs) for the three poliovirus serotypes as described by Martin et al. (27).
Reverse transcription-PCR restriction fragment length polymorphism (RT-PCR-RFLP) was used for molecular ITD (1). The sequence of the VP1 capsid protein coding region was determined as previously described (12).
Plaque purification. The isolate was purified by using a plaque assay in monolayer cultures of GMK cells and passaged twice in GMK cells at 36°C. The same virus stock (from now on referred to as PV3/EST02/E252, or E252) was used in all subsequent studies.
Sequencing and sequence analysis. Nearly complete genome sequencing of the isolate E252 was carried out as previously described (2) with primers designed by a primer-walking strategy. Electropherograms obtained from cycle sequencing (ABI Prism BigDye Terminator Cycle Sequencing Ready Reaction kit; Applied Biosystems, Espoo, Finland) in an automated sequencer (model 310; Applied Biosystems) were analyzed by Sequencing Analysis, v3.1, and Sequence Navigator, v.1 (Applied Biosystems). At least two corresponding electropherograms were performed and analyzed for each genomic region. Multiple sequence alignments were made using Clustal X (1.8). Distance matrices were estimated using the DNADIST and PROTDIST programs (PHYLIP, version 3.572c) (13) with the maximum-likelihood model of nucleotide substitution. Dendrograms were drawn with NEIGHBOR (PHYLIP) using the neighbor-joining option and visualized with NJPlot. Bootstrap analysis was performed using SEQBOOT (PHYLIP) with 1,000 replicates. Similarity analysis of the complete genome alignments was performed with the SimPlot software package (version 2.5) (25). The complete genome sequences of PV3/Sabin and PV1/Sabin used in comparisons were from GenBank (accession no. X00925 and V01150, respectively). The illustrations of the molecular structure of E252 were generated by the RasWin Molecular Visualization program (version 2.6; Glaxo Wellcome Research and Development, Stevenage, Herefordshire, United Kingdom) (36). The PV3/Sabin coordinates were from the Protein Data Bank, under identification no. 1PVC (14).
Assay for temperature sensitivity. Temperature sensitivity of the virus E252 was assayed by comparison of plaque formation in GMK and RD cells at 36 and 39.5°C. Serial 10-fold dilutions were prepared from stock viruses PV3/Sabin, E252, and PV3/Leon. Monolayer cultures of both cell lines on six-well plates were inoculated with 100 µl of virus dilutions. After adsorption for 30 min at either 36 or 39.5°C, the unabsorbed virus inoculums were removed. Cell monolayers were washed twice with Hanks' balanced salt solution and overlayed with 2 ml of 0.5% carboxymethyl cellulose in minimal essential medium supplemented with 1% fetal calf serum. The plates were incubated at 36 or 39.5°C and stained with crystal violet at 1 (RD) or 2 (GMK) days postinfection.
Microneutralization assay. A standard microneutralization assay was used to analyze neutralizing antibodies against the isolate E252 in human sera. Briefly, 50 µl of fourfold serum dilutions beginning at 1:4 was mixed with a standard amount of pretitrated virus (30 to 300 50% tissue culture infective doses [TCID50]/50 µl) in 96-well microtiter plates and incubated for 1 h at 36°C. After a further overnight incubation at 4°C, 2 x 104 GMK cells were added to each well, and the incubation was continued for 6 days at 36°C. The wells were stained with crystal violet, and the virus-induced cytopathic effect (CPE) was determined by photometry.
Human sera. The first study group included 15 OPV-vaccinated Estonian children ages 3 to 7 years. One serum specimen was drawn from each child between November 2002 and February 2003. The second study group consisted of 23 Finnish children vaccinated with three doses of inactivated polio vaccine (IPV) at ages 4, 6, and 12 months. Serum samples for antibody analysis were drawn at the age of 18 months (29).
Neutralizing MAbs. PV3-specific MAbs 204, 472, 495, 520, 838, 875, 879, 881, and 882, corresponding to antigenic sites 1 to 3, were kindly provided by Philip D. Minor (NIBSC).
Nucleotide sequence accession number. The nucleotide sequence of the complete poliovirus genome determined in this study has been submitted to GenBank under accession no. AY421739.
|
|
|---|
Environmental surveillance comprised monthly sampling of two locations in Tallinn. Five of the 13 samples examined during the second half of 2002 revealed echovirus 13, and 1 of them (22 October 2002), in addition, a PV3 labeled PV3/EST/02/E252 (E252), is characterized in detail below. As a response to this observation, four other locations in Tallinn were sampled in December 2002. One PV3 was isolated from the sample collected on 16 December and was shown to be typically Sabin-like. Furthermore, two 1-liter samples of raw sewage collected on 18 December 2002 and 22 January 2003 were concentrated at KTL in Helsinki. The latter sample revealed Sabin-like PV3. Samples collected and analyzed in Tallinn in early 2003 have revealed two isolates of PV2 with Sabin-like characteristics.
Preliminary characterization of the isolate. Serotype 3 of poliovirus E252 was confirmed by complete neutralization with polyclonal antisera specific for type 3. RT-PCR-RFLP revealed an atypical restriction pattern with one of the three enzymes (Hpa), while the other two (Dde and Hae) gave a result typical of PV3/Sabin. In EIA with cross-absorbed rabbit antisera to Sabin-like and non-Sabin-like PV3, the isolate E252 was not recognized by either antiserum. However, retesting under more sensitive assay conditions suggested the double-reactive nature of the isolate. In the second antigenic ITD assay based on MAbs (27), E252 was not neutralized with Sabin-specific MAbs. The nucleotide sequence of the VP1 capsid protein was 86.7% identical to the PV3/Sabin sequence (GenBank accession no. X00925).
Complete genomic sequence. The (almost) complete genomic sequence of isolate E252 was determined from nucleotide (nt) 19 to the end of the coding region and compared to that of PV3/Sabin and other sequence data available in GenBank. According to a similarity analysis with a 200-nt sliding window, the closest relative of E252 from the 5' noncoding region (NCR) close to the 3' end of the 2C coding region was PV3/Sabin, the rest of the genome being more closely related to PV1/Sabin (Fig. 1). The most divergent part of the genome of E252 was the approximately 100-nt-long hypervariable region in the 5' NCR. This region included two deletions (3 and 8 nt, respectively) compared to PV3/Sabin, and the identity to PV3/Sabin was only approximately 50%. Nearly complete genomic sequencing revealed multiple nucleotide substitutions throughout the genome as compared to PV3/Sabin or PV1/Sabin. The phylogenetic relationships of the E252 with other polioviruses and human enterovirus C (HEV-C) in VP1 and three-dimensional (3D) polymerase coding regions are shown in Fig. 2a and b, respectively.
![]() View larger version (36K): [in a new window] |
FIG. 1. Similarity analysis of complete PV genomes. Analysis was done with a sliding window of 200 nt moving in steps of 20 nt, with all positions with a gap deleted, and the nucleotide similarity was plotted by using the JC model of nucleotide substitution. (a) E252 is compared against PV1/Sabin and PV3/Sabin (b) E252 is compared against 12 different poliovirus or HEV-C genomes. The GenBank accession numbers of the sequences used are as follows: PV3/Sabin, X00925; PV3/23127, X04468; PV3/Usol-d-bac, AJ293918; PV1/Sabin, V01150; PV1/DOR00013, AF405690; PV1/RUS1161-96, AF462419; PV2/Sabin, X00595; CV-A21, D00538; CV-A24, D90457; and PV1/66months, AJ132961.
|
![]() View larger version (22K): [in a new window] |
FIG. 2. Neighbor-joining tree of the nucleotide sequence relationships of the isolate E252 and other polioviruses and HEV-C in VP1 (a) and 3D polymerase (b) protein coding regions. The location of E252 is indicated by an arrow. Numbers at nodes represent reliability values for the internal branch (= bootstrap values). Values over 70% are shown. The scale bar indicates the evolutionary distance (0.01 = 1%). The GenBank accession numbers for the sequences are as follows: CV-A1, AF081293 and AF499635; CV-A11, AF081301 and AF499636; CV-A13, AF081303 and AF499637; CV-A15, AF081305 and AF499638; CV-A17, AF081306 and AF499639; CV-A18, AF081307 and AF499640; CV-A19, AF081308 and AF499641; CV-A20, AF081309 and AF499642; and CV-A22, AF081310 and AF499643. For information on other sequences, see the legend to Fig. 1.
|
|
View this table: [in a new window] |
TABLE 1. Amino acid sequences of different antigenic sites of PV3
|
![]() View larger version (87K): [in a new window] |
FIG. 3. Predicted molecular structure of the protomer E252. A space-fill model of a protomer (14) is shown at three different angles. The picture was generated with the RasWin Molecular Visualization Program (version 2.6) (36). (a) Front view. (b) Left-side view. (c) Right-side view. VP1, blue; VP2, yellow; VP3, red; VP4, green. A changed amino acid is shown by cyan, an antigenic site is in white, and a changed amino acid in an antigenic site is in violet. Identification of amino acids by the one-letter code: T, threonine; K, lysine. The first number after the one-letter code indicates the capsid protein: 1, VP1; 2, VP2; and 3, VP3. The three-number code indicates the position of the amino acid in a given capsid protein.
|
Antigenic properties of the isolate. Isolate E252 was neutralized completely with polyclonal PV3-specific antisera, but not with a mixture of Sabin-specific MAbs used in the ITD assay. To further characterize the antigenic properties of the isolate, neutralizing activities of nine antigenic site 1- to 3-specific MAbs were determined with a standard microneutralization assay. The isolate E252 was not neutralized with any of the tested PV3-specific MAbs. In contrast, PV3/Sabin was neutralized with eight and wild-type PV3/Saukett was neutralized with six of the nine tested MAbs (Table 2).
|
View this table: [in a new window] |
TABLE 2. Reactivity of isolate E252 with MAbs specific for different antigenic sites (1, 2, and 3) of PV3
|
![]() View larger version (10K): [in a new window] |
FIG. 4. Serum neutralizing antibodies to PV3 strains E252 and Sabin in children vaccinated with IPV (A) or OPV (B). One serum sample from each child was analyzed for neutralizing antibodies by using a standard microneutralization assay with fourfold serum dilutions. The negative sera (titers of <4) were blotted, as they were 1:2.
|
|
View this table: [in a new window] |
TABLE 3. Temperature sensitivity of E252
|
|
|
|---|
The initial sequencing of the capsid protein VP1 coding region suggested that the closest relative to the isolate E252 was the oral vaccine strain PV3/Sabin, but the two viruses shared only 86.7% of the VP1 nucleotides. The attenuated Sabin strain and its parental strain PV3/Leon belong to an extinct genotype. Putative reemergence of related wild strains is highly unlikely. Complete genomic sequencing revealed multiple nucleotide substitutions, compared to Sabin 3, throughout the genome. These included all the genetic markers of reverted, neurovirulent PV3/Sabin-derived virus strains. Therefore, we did not test the neurovirulence of the E252 virus in experimental animals. There is plenty of evidence for increased neurovirulence in strains that are genetically much closer to the vaccine strain than E252.
A probable Sabin 3/Sabin 1-recombination junction was found in the 2C coding region. It is noteworthy, however, that the sequences of the P3 regions are not serotype specific and those of HEV-C are not necessarily separable from the corresponding sequences found in polioviruses. Therefore, theoretically it is also possible that the P3 region of the E252 virus was derived from an unidentified genetic cluster C virus.
Amino acid substitutions were seen in most of the viral proteins. They were relatively less frequent in the nonstructural proteins than in the major capsid proteins. In the capsid proteins VP1, VP2, and VP3, most of the substitutions were located near the virion surface, where the degree of freedom of polypeptide conformation is greater than in the beta barrels. Several substitutions were seen at the known neutralization antigenic sites. However, we do not know if this was due to immune selection in the host or had occurred by successive random selections associated with putative bottleneck transmission events over the years. The numerous substitutions at the antigenic sites were reflected in the strikingly altered antigenic properties of the virus. None of the nine tested MAbs recognizing one or both of the tested PV3 reference strains was able to neutralize the E252 virus. Furthermore, the capacity of human sera from groups of IPV- or OPV-immunized children to neutralize the virus was significantly reduced compared to that of the corresponding vaccine strains. From this point of view, the virus should have had a potential for transmission in both Estonian and Finnish populations. No evidence for this type of circulation was obtained.
The route of the virus to the sewerage remains unresolved. There have not been any "polio-like" epidemics or even suspected cases of poliomyelitis in Estonia during the time the virus was isolated. All the fecal specimens collected from patients with meningitis or a similar diagnosis during autumn 2002 were subjected to reisolation attempts. In spite of these duplicate analyses of the specimens, no sign of poliovirus circulation was seen before or after sewage isolate E252. The environmental surveillance was intensified by collecting more wastewater samples from distinct areas of Tallinn, but this did not lead to a second isolation of a similar virus. It seems highly unlikely that the numerous substitutions in the genome of E252 would have been generated by widespread circulation of the virus in the local population.
Importation of the evolved virus from abroad cannot be excluded. Wild-type PV3, most probably after importation, was isolated from Tallinn sewage in 1987 and 1988 in the absence of reported cases of poliomyelitis (20). Tallinn is a busy harbor city with numerous visiting ships every week. Occasionally the ships empty their waste containers into the sewers of the City of Tallinn, and one could speculate about this route of introduction of the virus into the city sewers. However, according to the harbor records, no ship emptied its waste containers into the sewerage of Tallinn in October 2002.
The origin of the virus E252 is unknown, but one likely reason for emergence of this virus is a long-term replication in an immunodeficient individual. A Sabin 3/Sabin 1 recombination is typical of most recipients of OPV, whereas the P3 region in discovered cVDPVs has been derived from an unidentified genetic group C virus (including the HEV-C and poliovirus species) (6, 8, 21, 41). There are two problems with this hypothesis, however. First, no poliovirus-excreting immunodeficient patients are known in Tallinn. This does not, of course, exclude their existence, because humoral immunodeficiency sufficient to allow chronic poliovirus infection sometimes can be asymptomatic. Second, the virus was isolated by a trap sampling method. Theoretically, the likelihood of detecting virus excreted by a single person in a sewage system where excreta from a given person are mixed with sewage derived from half a million nonexcretors, is very small (19, 32). The nonselective adsorbent immersed in the sewage is likely to be rapidly saturated with unrelated suspended materials abundant in a city's sewage. We have to assume that part of the E252 virus excreted by the putative immunodeficient individual happened to be in the very location and at the very time where and when the gauze pad was immersed in the sewage.
The AFP surveillance is the "gold standard" in the WHO program aiming at the global eradication of wild-type poliovirus. Environmental surveillance, recognized as a supplementary tool, has shown to be extremely powerful in some situations. For example, a recent survey in Egypt revealed prolonged and extensive circulation of wild-type PV1 at a time when the poliovirus-positive AFP cases were already rare (12). Environmental surveillance has also previously yielded isolation of diverged vaccine-derived strains with a wild-type phenotype without any concurrent reported cases of paralytic poliomyelitis. For example, a PV2 strain recovered from sewage in Israel, in the absence of any polio cases, had 91.7% nucleotide identity with PV2/Sabin in the P1 capsid region and had regained some wild-type characteristics (37).
In conclusion, we have described genetic and phenotypic features of a PV3/PV1 recombinant VDPV that appears to have characteristics of an iVDPV, although the putative chronically infected individual is not known. No evidence for circulation of this virus was obtained, although the observed antigenic changes would help the virus to escape vaccine-induced immunity in human populations. The existence of this kind of viruses in human populations is a serious risk, considering the option of abruptly stopping the use of OPV after the desired and continuing global eradication of wild-type poliovirus-dependent poliomyelitis.
This work was partially supported by WHO TSA I8/181/526 to KTL.
|
|
|---|
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»