Previous Article | Next Article ![]()
Journal of Virology, September 2007, p. 8944-8952, Vol. 81, No. 17
0022-538X/07/$08.00+0 doi:10.1128/JVI.00180-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Dorsaf Nasri,1,2,
Lionel Chollet,3
Khaoula Belguith,1
Thomas Bourlet,2
Mahjoub Aouni,1
Bruno Pozzetto,2* and
Sylvie Pillet2
Laboratory of Transmissible Diseases and Biologically Active Substances, Faculty of Pharmacy, Monastir, Tunisia,1 Laboratory of Bacteriology-Virology, GIMAP EA3064, Faculty of Medicine of Saint-Etienne, Saint-Etienne, France,2 Laboratory of Bacteriology-Virology, Hospital of Toulon-La Seyne sur Mer, Toulon-La Seyne sur Mer, France3
Received 26 January 2007/ Accepted 21 May 2007
|
|
|---|
|
|
|---|
The RNA genome, 7.5 kb long, is constituted by a single open reading frame flanked by 5' and 3' untranslated regions. The coding region is translated into a single polyprotein of 2,200 amino acids and is then processed to generate four structural proteins (VP1 to VP4) and seven nonstructural proteins (2A to 2C and 3A to 3D). The four structural proteins constitute the icosahedral capsid that contains the major antigenic determinants (22) and the principal attachment sites to the cellular receptors (for a review, see reference 10).
Recombination between two strains is a known phenomenon for enteroviruses replicating within a single cell. Although recombination has long been recognized as an important property of PV (9, 17), several recent publications have demonstrated that it is also extremely frequent in non-PV enteroviruses (2, 6-8, 15, 17-20, 24, 27, 28, 31, 35). Despite this high level of genetic instability, the occurrence of intra- and interserotypic recombination events in the VP2-VP3-VP1 structural coding region has been shown to be a rare phenomenon and up to now restricted to PV strains (5, 14, 21, 40). Structural requirements of the virion shell or of receptor binding were thought to be involved in the restriction of recombination within this region (17, 30, 36).
In this paper, we describe an interserotypic recombination event occurring in the VP3 coding region of a clinical strain of HEV-B and leading to a chimeric coxsackievirus B3 (CV-B3)/CV-B4 type.
|
|
|---|
Identification of strain SE-03-78616 was performed by neutralization of the cytopathic effect using A to H and J to P Lim-Benyesh Melnick intersecting pools of specific monovalent antisera (Statens Serum Institute, Copenhagen, Denmark) (23). Further neutralization experiments were conducted with specific polyclonal antisera directed toward the six CV-B serotypes (Statens Serum Institute, Copenhagen, Denmark). Clinical isolates of CV-B3, CV-B4, and CV-B5 were used as positive controls.
Molecular typing. The molecular typing of strain SE-03-78616 was performed by partial sequencing of the VP1 gene, as described previously (29), and of the VP2 gene using original primers (AM12 and AM32) (Table 1). Viral RNA was extracted from 140 µl of tissue culture supernatant and then reverse transcribed using SuperScriptIII reverse transcriptase (Invitrogen, Cergy Pontoise, France). Amplification reactions were done with a proofreading Taq DNA polymerase (Platinum Pfx; Invitrogen, Cergy Pontoise, France). The amplified products were purified with the QIAquick gel extraction kit (QIAGEN, Courtaboeuf, France) and sequenced with the Genome Lab Dye Terminator Cycle Sequencing Quick Start kit (Beckman Coulter, Villepinte, France) according to the manufacturer's recommendations. Sequence alignment was performed with BLAST software (1).
|
View this table: [in a new window] |
TABLE 1. Primers used for amplification of the P1 region
|
Sequence, recombination, and phylogenetic analyses. The P1 region (VP4 to VP1) sequence of the strain was compared to that of reference and clinical CV-B strains (http://www.picornaviridae.com/enterovirus/hev-b/hev-b_seqs.htm) with Clustal W (version 1.81) (39). Similar alignments were performed for the P2 (2A to 2C) and P3 (3A to 3D) regions by comparison to reference and clinical HEV-B strains.
The location of recombination events was determined by similarity plotting of the resulting alignment and bootscanning analysis (34) using the SimPlot program, version 2.5 (http://sray.med.som.jhmi.edu/SCRoftware/SimPlot) (16), with a sliding window of 200 nucleotides moving in 20-nucleotide steps.
Phylogenetic trees were designed by imputing the aligned sequences into the MEGA program (version 3.1) (13) and constructed with the neighbor-joining algorithm (33). Genetic distances were calculated with the Kimura-2 parameter model (12) with a transition/transversion ratio of 2.0, and the reliability of the trees was determined by bootstrap analysis with 1,000 pseudoreplicate data sets (11).
Nucleotide sequence accession numbers. The GenBank accession numbers of the sequences used in the molecular analyses are provided in Fig. 1, 2, 3, and 5.
![]() View larger version (33K): [in a new window] |
FIG. 1. Similarity plot and bootscanning analysis of the P1 region comparing strain SE-03-78616 and CV-B strains with complete genomes available in databases. (Top) Similarity plot of strain SE-03-78616 versus CV-B strains constructed by the SimPlot program, version 2.5, with a sliding window of 200 nucleotides moving in 20-nucleotide steps. (Bottom) Bootscan graph of strain SE-03-78616 versus CV-B strains constructed by using the neighbor-joining tree algorithm, the Kimura-2 distance model, and 100 pseudoreplicates. The dotted line shows an arbitrary 70% reliability threshold. The following sequences (with GenBank accession numbers given in parentheses) were used: CV-B1 (AY186745 to AY186748, M16560); CV-B2 (AF081485, EF174468, EF174469); CVB-3 (AF231763 to AF231765, AY752944 to AY752946, M16572, M33854, M88483), CV-B4 (AF311939, DQ480420, X05690), CV-B5 (AF114383, AY875692, X67706), CV-B6 (AF039205).
|
![]() View larger version (38K): [in a new window] |
FIG. 2. Location of the breakpoint site of strain SE-03-78616 by nucleotide sequence alignment with the CV-B3 and CV-B4 strains listed in Fig. 1 (1B3, AY752949; 2B3, AF231765; 3B3, AF231762; 4B3, M88483; 5B3, M33854; 6B3, AY752945; 7B3, AY752944; 8B3, AF231763; 9B3, M16572; 10B3, U57056; 1B4, DQ480420; 2B4, X05690; 3B4, AF311939). Regarding positions exhibiting base homology for all tested strains of a single serotype but divergence between the two serotypes, the nucleotide shared with strain SE-03-78616 is shown against a gray background. The breakpoint site is represented by a black background.
|
![]() View larger version (19K): [in a new window] |
FIG. 3. Phylogenetic relationships in the P1 region between strain SE-03-78616 and the CV-B strains listed in Fig. 1. (A) Phylogenetic tree of the genomic region situated upstream of the breakpoint site (1,206 nucleotides, position 745 to 1950). (B) Phylogenetic tree of the genomic region situated downstream of the breakpoint site (1,357 nucleotides, position 1950 to 3306). Only bootstrap values over 75% are shown at tree nodes. (C) Nucleotide divergence (derived from the MEGA program [13] and expressed as a percentage) between strain SE-03-78616 and CV-B3/CV-B4 strains (abbreviated according to the legend to Fig. 2), upstream and downstream of the breakpoint site.
|
![]() View larger version (25K): [in a new window] |
FIG. 5. Phylogenetic tree showing relationships between strain SE-03-78616 and HEV-B strains in the P2-P3 region. Only bootstrap values over 75% are shown at tree nodes. The HEV-B strains are indicated by serotype and accession number in the GenBank database. Asterisks identify the strains used in Fig. 3.
|
|
|
|---|
Molecular characterization of strain SE-03-78616. To further explore the discrepant results obtained by two different molecular typing strategies, the entire P1 region of the viral genome was sequenced. As depicted in Fig. 1, the obtained sequence (2,562 nucleotides) was shown to be close to that of CV-B3 in the 5' part of the genome until the first third of VP3 and to that of CV-B4 after this position, suggesting the occurrence of an interserotypic recombination event in the capsid region. To locate precisely the breakpoint, the sequence was compared by alignment with other CV-B3 and CV-B4 strains. Figure 2 shows that the recombination site corresponds to nucleotide 1207 of the coding region (nucleotide 1950 with reference to the entire genome). For all positions shared by strains of the same serotype but diverging between the two serotypes, strain SE-03-78616 clustered with CV-B3 and CV-B4 strains upstream and downstream of the breakpoint, respectively (Fig. 2). One of the primer pairs used to sequence the structural coding region (RC3 and RC4 [Table 1]) covered unambiguously the recombinant region.
Phylogenetic analyses (Fig. 3) confirmed that strain SE-03-78616 is the result of a recombination event in the VP3 coding region between two CV-B types. Despite the fact that, in the region located upstream of the breakpoint, this strain is relatively distant from the other strains of CV-B3 used to construct the tree depicted in Fig. 3A, the nucleotide divergence is less than 25% (Fig. 3C), confirming that the strain belongs to the CV-B3 serotype in this region. By contrast, downstream of the breakpoint, the phylogenetic profile and nucleotide divergence data identify the strain as CV-B4 (Fig. 3B and 3C, respectively).
Functional analysis of the breakpoint site. One hundred 50% tissue culture infective doses of three clones of strain SE-03-78616 were completely neutralized by a horse polyclonal anti-CV-B3 monovalent antiserum diluted 1:50 to 1:200 for at least 72 h; in contrast, no neutralization was observed under the same experimental conditions by antisera directed toward CV-B1, CV-B2, CV-B4, CV-B5, or CV-B6.
In order to further explore the putative relationship between the antigenic sites and the location of the breakpoint, an amino acid alignment was performed in the whole capsid region with reference to the available sequences of CV-B3 and CV-B4 strains (Fig. 4). The topography of the antigenic sites was represented according to the mapping described by Auvinen et al. (3). As shown in Fig. 4, the recombination site was found to be located in the central part of the 3b antigenic site in the VP3 protein. In addition, this alignment illustrates the analogy of amino acid sequences with CV-B3 in the 5' part of the recombinant strain and CV-B4 in the 3' part, confirming that this strain is the result of a recombinant event between these two CV-B serotypes. This finding is reinforced by the identification of significant differences in the amino acid sequence (deletions or insertions) between the two serotypes, notably in antigenic site 2b of VP2, in antigenic site 1 of VP1, and in the 3'-terminal part of the latter protein (Fig. 4).
![]() View larger version (57K): [in a new window] |
FIG. 4. Amino acid sequence alignment of the structural proteins of strain SE-03-78616 with the CV-B3 and CV-B4 strains listed in Fig. 1 and abbreviated according to the legend to Fig. 2. Regarding positions exhibiting amino acid homology for all tested strains of a single serotype but divergence between the two serotypes, the amino acid shared with strain SE-03-78616 is shown against a gray background. Deletions or insertions having the same characteristics are shown against a black background. The solid vertical lines delineate the four structural viral proteins (VP). The dotted vertical line locates the CV-B3/CV-B4 switch on the amino acid sequence of the recombinant strain. The antigenic sites described for CV-B3 by Auvinen et al. (3) are depicted in bold type and designated according to the nomenclature of those authors. The amino acids are numbered according to the sequence of the reference strain of CV-B3 (accession number M16572). The positions of the main tertiary structures described for CV-B3 (25) are indicated by double-headed horizontal arrows above the sequences.
|
![]() View larger version (21K): [in a new window] |
FIG. 6. Similarity plot (Top) and bootscanning analysis (Bottom) of the P2-P3 region comparing strain SE-03-78616 with CV-B3 (U57056), CV-B4 (X05690), and E-30 (AM237325) strains. The graphs were constructed as described in the legend to Fig. 1.
|
|
|
|---|
In this study, we analyzed a strain of HEV-B that represents an exception to the above rule, as it was clearly found to exhibit a recombination event within the VP3 gene, with a CV-B3 genome at the 5' end and a CV-B4 genome at the 3' end, as evidenced by bootscanning, alignment, nucleotide divergence, and phylogenetic analyses. A technical artifact of the polymerases was ruled out, since the recombination event was confirmed on three different clones. Except for PV serotypes that have been shown to recombine within the VP1 gene (5, 21, 40) and the VP2 genes (14), this is the first report of an intertypic recombination event within the capsid region of HEV.
The rarity of this naturally occurring phenomenon could be due to the fact that the recombinant strain would not be able to recognize its cellular receptor; however, in the present case, the two serotypes participating in the recombination have the same receptor (10). Another explanation could be linked to capsid constraints within the VP2-VP3-VP1 block that prevent the occurrence of recombination events. For strain SE-03-78616, the precise location of the recombination site is position 1950 of the VP3 gene. As shown in Fig. 4, the breakpoint is located in the ßB barrel at the junction between the knob and the BC loop previously described in the tertiary structure of CV-B3 (25), a region known to harbor the 3b antigenic site of this serotype (3).
With regard to neutralization epitopes, it is striking to note that the recombinant strain was neutralized by a CV-B3 monovalent polyclonal antiserum but not by the Lim-Benyesh-Melnick pools of antisera. This finding could be explained by the fact that the two anti-CV-B3 neutralizing antisera are not directed toward the same epitope; alternatively, the dilution inherent to the constitution of pooled antisera could account for the lower sensitivity of this reagent. More interestingly, although strain SE-03-78616 exhibited a CV-B4-like VP1 region, it was not recognized by the CV-B4 monovalent antiserum. This finding suggests that the major neutralization antigenic site(s) is located in the proximal part of the VP2-VP3-VP1 block. In fact, Beatrice et al. reported that the VP2 peptide was the sole immunogen capable of inducing neutralizing antibodies against CV-B3 (4). These findings support the fact that typing methods based only on the VP1 region may underestimate the weight of the VP2/VP3 region in the antigenic variability of enteroviruses.
In complement to the description of a recombination event within the capsid region of strain SE-03-78616, the occurrence of another breakpoint was documented in its nonstructural coding region. This observation illustrates that recombination in the latter part of the enterovirus genome, which is very common in HEV-B circulating strains (36), is also possible in an isolate exhibiting a recombined capsid region.
Finally, from the taxonomic point of view, this finding suggests the putative existence of chimeric viruses within the HEV-B species which could be considered equivalents to the recombinant forms between different clades of human immunodeficiency virus type 1 strains belonging to subtype M (26), adding to the complexity of type identification within the HEV-B species. The extensive use of typing methods targeting two different parts of the capsid region, VP1 and VP2, will contribute to an evaluation of the impact of interserotypic recombination in terms of enterovirus identification and functional activity.
Published ahead of print on 30 May 2007. ![]()
Lamjed Bouslama and Dorsaf Nasri contributed equally to this work. ![]()
|
|
|---|
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»