J Gen Virol Tips for Better Browsing
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 87 (2006), 375-385; DOI 10.1099/vir.0.81298-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zell, R.
Right arrow Articles by Wutzler, P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zell, R.
Right arrow Articles by Wutzler, P.
Agricola
Right arrow Articles by Zell, R.
Right arrow Articles by Wutzler, P.
© 2006 Society for General Microbiology

Molecular-based reclassification of the bovine enteroviruses

Roland Zell1, Andi Krumbholz1, Malte Dauber2, Elizabeth Hoey3 and Peter Wutzler1

1 Institute for Virology and Antiviral Therapy, Hans-Knöll-Str. 2, 07745 Jena, Germany
2 Institute for Virus Diagnostics, Friedrich Loeffler Institute, Federal Research Institute for Animal Health, Boddenblick 5a, 17493 Insel Riems, Germany
3 School of Biology & Biochemistry, Medical Biology Centre, The Queen's University of Belfast, UK

Correspondence
Roland Zell
Roland.Zell{at}med.uni-jena.de


   ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Bovine enteroviruses are currently classified into two serotypes within the species Bovine enterovirus (BEV). Comparison of the sequences of six American and eleven German BEV isolates with published BEV sequences revealed the necessity to revise the taxonomy of these viruses. Molecular data indicate that the bovine enteroviruses are composed of two clusters (designated BEV-A and -B) each with two and three geno-/serotypes, respectively. Whereas low amino acid identity of the capsid proteins 1C (VP3) and 1D (VP1) is the main criterion for the discrimination of geno-/serotypes, the BEV clusters, presumably representing species, differ in sequence identity of all viral proteins. In addition, characteristic lengths of (i) the capsid proteins 1B, 1C and 1D, (ii) the 2C protein, and (iii) the 3'-non-translated region are observed. The BEVs can be distinguished from the other enteroviruses by sequence identity and unique features of the 5'-non-translated region, i.e. a conserved second cloverleaf and characteristic RNA structures of the internal ribosome entry site. Phylogenetically, the closest relatives of the bovine enteroviruses are the porcine enteroviruses. Incongruent phylogenies of the 5'-non-translated region, the capsid proteins and the 3D polymerase indicate frequent intraserotypic and interserotypic recombination within the non-capsid and the capsid region of the BEV genome.

The GenBank/EMBL/DDBJ accession numbers for the nucleotide sequences reported in the paper are DQ092769–DQ092795.


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Enteroviruses of the family Picornaviridae are small, non-enveloped viruses with an icosahedral virion and a positive-stranded RNA genome. Currently, five human and several animal enterovirus species comprise the genus Enterovirus (Stanway et al., 2005Down). The classified animal enteroviruses include simian, porcine and bovine isolates. The first bovine enteroviruses were collected in the late 1950s (Kunin & Minuse, 1958Down; McFerran, 1958Down; Moll & Davis, 1959Down; see also: Barya et al., 1967Down; Dunne et al., 1974Down). They were isolated from (i) faeces of animals with symptoms of pneumonia, respiratory disease, enteritis, dysentery and fertility disorders, (ii) fetal fluids of an aborted calf, and (iii) the faeces of apparently healthy animals (Table 1Down). Difficulties in reproducing clinical symptoms, following experimental infection of animals, led to the conclusion that bovine enteroviruses were of only minor veterinary medical importance. Early attempts to classify bovine enteroviruses by serological means were reported by Huck & Cartwright (1964)Down, Barya et al. (1967)Down and Dunne et al. (1974)Down. In those studies seven, four and eight serogroups, respectively, were suggested. The proposal of Dunne et al. (1974)Down was approved by the WHO/FAO and prototype strains of seven serotypes were deposited in the ATCC, while the eighth serotype (designated PS 35) was subsequently found to be identical to PS 83, the proposed serotype 5. Later, this classification was revised following a proposal of Knowles & Barnett (1985)Down. According to their work, bovine enteroviruses are classified into two serotypes, which are currently within the species Bovine enterovirus (BEV) of the genus Enterovirus (Stanway et al., 2005Down). However, after accumulation of enterovirus sequence data (Earle et al., 1988Down; McNally et al., 1994Down; McCarthy et al., 1999Down; Goens et al., 2004Down), concerns arose as to whether the present classification correctly represents the genetic diversity of bovine enteroviruses. In the present study, we sequenced six BEV strains that were collected in the USA (between 1957 and 1962) and deposited in the ATCC, as well as 11 recent field isolates from Germany (collected between 1982 and 2003). Among the American isolates are five of the previously suggested prototype strains. The sequence data were compared with published BEV sequences. Our genome data and subsequent phylogenetic analyses support the establishment of two BEV species, one of which contains two and the other three geno-/serotypes.


View this table:
[in this window]
[in a new window]
 
Table 1. BEV isolates

 

   METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Viruses and cells.
The viruses used in this study are described in Table 1Up. BEV strains PS 42 (VR-758), PS 83 (VR-757), PS 87 (VR-774) and PS 89 (VR-755) were obtained from the ATCC. LC-R4 was obtained from the National Veterinary Research Institute, Pulawy, Poland, and BEV-261 from the Institute for Animal Health (Pirbright, UK). Prior to RNA preparation, the viruses were propagated in Madin–Darby bovine kidney (MDBK) cells. This cell line was grown in Dulbecco's modified Eagle's medium supplemented with 10 % fetal bovine serum. Medium and supplements were obtained from Gibco-BRL Life Technologies.

Generation of BEV amplicons by long RT-PCR.
RNA of virus-infected MDBK cells was prepared by using the Perfect RNA kit of Eppendorf. Five micrograms of total RNA was reverse-transcribed with 20 pmol oligo(dT)20 primer and 40 U RevertAid H Minus M-MuLV reverse transcriptase (supplied by MBI Fermentas) in a total volume of 20 µl. Two microlitres of cDNA was used for PCR employing several primer sets. Different protocols were used depending on the expected size of the PCR product. For PCR products of approximately 1 kbp the following protocol was used: 0·2 mM each dNTP, 1 µM each primer, 10 mM Tris/HCl pH 8·3, 1·5 mM MgCl2, 50 mM KCl, 2·5 U Taq polymerase (total volume 50 µl). The PCR cycle was the same as used for routine diagnostics, i.e. 35 cycles of 50 s 94 °C, 50 s 55 °C, 1 min 72 °C. For the amplification of PCR products larger than 1·5 kbp, the Combizyme DNA polymerase mix supplied by InViTek was used. The PCR cycle was modified in the following way: 35 cycles of 30 s 92 °C, 50 s 55 °C, 8 min 68 °C. The 5'- and 3'-ends of the genomes were amplified using the 5'/3' RACE kit of Roche.

Sequencing of the BEV amplicons, sequence alignments and phylogenetic analysis.
PCR products were analysed by electrophoresis on agarose gels stained with ethidium bromide. PCR products of the expected size were purified by using either the QIAquick PCR Purification kit or the QIAquick Gel Extraction kit (Qiagen). Sequencing was performed according to the cycle sequencing protocol of ABI. The products were analysed on a Prism 310 Genetic Analyser of ABI. The sequences were aligned manually or with the help of CLUSTAL W (Thompson et al., 1994Down). Neighbour-joining trees were calculated with the quartet puzzling method (Strimmer & von Haeseler, 1996Down; Strimmer et al., 1997Down) using the JTT substitution model for amino acid sequences (Jones et al., 1992Down) and the Tamura–Nei model for nucleotide sequences (Tamura & Nei, 1993Down). The reliability of the clustering was tested by 10 000–25 000 iterations in the quartet puzzling method depending on the complexity of the dataset. For tree construction, maximum-likelihood branch lengths were computed. Consistency of branching was also tested with the maximum-parsimony algorithm using the PHYLIP program package (Felsenstein, 1995Down) and PROTML program (Adachi & Hasegawa, 1992Down).

RNA secondary predictions and free-energy calculations were performed with the MFOLD program (version 3.0; Zuker et al., 1999Down).

Nucleotide sequence accession numbers.
The nucleotide sequence data reported in this paper were submitted to the GenBank nucleotide sequence database (accession nos DQ092769–DQ092795, see also Table 2Down). For sequence alignments, available sequences of bovine, human, porcine and simian enteroviruses as well as human rhinoviruses were used (for the complete compilation of virus strains and GenBank accession numbers used see Supplementary Table; available in JGV Online).


View this table:
[in this window]
[in a new window]
 
Table 2. Classification of BEV strains

 

   RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Comparison of the capsid protein- and 3D polymerase-encoding gene regions
The RNA sequences and the genome organization of several BEV isolates have been previously reported (Earle et al., 1988Down; McNally et al., 1994Down; McCarthy et al., 1999Down; Goens et al., 2004Down). Briefly, a long open reading frame encodes a polyprotein ranging from 2166 to 2176 aa. The sequences of the polymerase-encoding genome region, which are clearly distinct from the human enteroviruses, suggested the establishment of a discrete species named BEV. However, the low degree of sequence identity between the capsid protein sequences of strains VG(5)27, K 2577 and SL 305 (all serotype 1 according to the classification of Knowles & Barnett, 1985Down) and those of strains PS 87 and RM2 (both serotype 2) prompted us to study bovine enteroviruses in more detail. The aim of this study was to compare the complete sequences of representative strains of the seven serotypes previously proposed by Dunne et al. (1974)Down. With the exception of BEV-165, which has been previously sequenced (M. McIlhatton, personal communication), all strains deposited in the ATCC by C. M. Kunin and H. W. Dunne were sequenced. This collection included five of the suggested prototype strains (BEV-261, PS 42, PS 83, PS 87 and PS 89) and LC-R4, which belongs to the same serogroup as BEV-1 and VG(5)27. Moreover, 11 German isolates collected between 1982 and 2003 were also included in this study. For the latter viruses, a minimum of the 5'-NTR, the P1 coding region, the 3D polymerase coding region and the 3'-NTR were sequenced. The BEV sequences were compared with representative human, porcine and simian enterovirus as well as with human rhinovirus sequences taken from the GenBank database (see Table 2Up).

In order to solve the question of whether there are one or more BEV species, the pairwise sequence relationships for each of the genome-encoded capsid proteins 1A–1D and 3D were calculated. In previous studies, this approach provided useful evidence for species distinction (Van Regenmortel et al., 1997Down; Zell et al., 2001Down). Three frequency peaks of amino acid identity scores are observed for the capsid proteins 1C and 1D, and two peaks for the 3D polymerase (Fig. 1Down). These distribution patterns indicate the existence of heterogeneous geno-/serotypes and a higher order taxon, presumably a species.


Figure 1
View larger version (32K):
[in this window]
[in a new window]
 
Fig. 1. Frequency distribution of pairwise amino acid identity scores of the four capsid proteins (1A–1D) and the 3D polymerase. Amino acid sequences of 22 BEV strains and field isolates were compared with each other in order to calculate amino acid identities. The discontinuous frequency distribution of the plotted amino acid identities indicates the existence of heterologous geno-/serotypes and species.

 
Whereas the protein sequences of the P2/P3 region and 1A (VP4) are highly conserved, the capsid proteins 1B (VP2), 1C (VP3) and 1D (VP1) show characteristic differences suited for the discrimination of geno-/serotypes. The C terminus of the 1B protein of the A cluster is 2 aa longer than that of 1B of the B cluster (overall protein lengths 248 vs 246 aa) and the GH-loop of the 1C protein is 1 aa shorter (overall protein lengths 242 vs 243 aa). The greatest heterogeneity is observed for the 1D protein, which shows differences at the N terminus (deletion of 5 aa in BEV-B), the BC-loop (deletion of 2 aa in BEV-B), the DE-loop (insertion of 1 aa in BEV-B serotype 1, deletion of 1 aa in BEV-A serotype 2), the HI-loop (insertion of 1 aa in BEV-B) and the C terminus (deletion of 1 aa in BEV-B, insertion of 2 aa in strain D 14/3/96). With the exception of 2C, all P2/P3 proteins have a highly conserved primary structure. The 2C protein is thought to be a membrane-associated ATPase and has a conserved nucleotide-binding motif (Möller & Amons, 1985Down). The 2C protein of the BEV cluster B is 1 aa shorter than the 2C of the A cluster.

The phylogenetic analysis of the P1 polyprotein included four HRV strains as outgroups (HRV2, 3, 9 and 14) and at least one member of each recognized enterovirus species. The phylogenetic tree reveals two rhinovirus clusters and nine enterovirus clusters (Fig. 2Down). These clusters represent two BEV clusters (designated BEV-A and BEV-B), as well as the acknowledged enterovirus species SEV-A, HEV-A to -D, PV and PEV, and the rhinovirus species HRV-A and -B. Each BEV cluster is further subdivided, presumably representing geno-/serotypes. In cluster A, group 1 composes strains D 14/3/96, Vir 404/03, LC-R4, VG(5)27, VD 2860/1-9, D 8/01, D 58/96-V2130 and E 6-82, and group 2 includes BEV-165 (M4), PS 42, PS 83, D 3/98, 56/59/1, SD 1182 II, D 14/1/96, SL 305 and K 2577. In cluster B, group 1 consists of BEV-261 (also known as M2) and RM2, group 2 includes PS 89, PS 87/Maryland, Wye-3A and Jena 38/02, whereas PS 87/Belfast is the only member of group 3. Genetic clustering of the BEV geno-/serotypes correlates in part with the former classification of Dunne et al. (1974)Down as LC-R4, BEV-261 (M2), BEV-165 (M4), PS 87/Belfast and PS 89 were previously described as discrete serotypes (see Table 2Up). However, PS 42 and PS 83, which were also described as discrete serotypes, are very similar to each other. One nucleotide insertion within the 5'-NTR and 33 nt substitutions (corresponding to 12 aa substitutions) were observed. Presumably, the substitutions reflect different passage histories of otherwise identical viruses. Moreover, both viruses are closely related to BEV-165, suggesting that they belong to the same geno-/serotype of the BEV-A cluster. BEV-261 (M2) and RM2 have very similar P1 sequences (Fig. 2Down) and are serologically almost indistinguishable (McNally et al., 1994Down). Obviously, these sequences represent different passage histories of the same virus.


Figure 2
View larger version (27K):
[in this window]
[in a new window]
 
Fig. 2. Phylogenetic relationships of the BEVs to the entero- and rhinoviruses. Unrooted neighbour-joining tree of the P1 capsid protein precursor of 58 entero- and rhinoviruses using human rhinovirus 2 as an outgroup. Each acknowledged species is represented by at least two viruses (exception: SEV-A). Amino acid sequences were aligned with the CLUSTAL W program and adjusted manually. Maximum-likelihoodbranch lengths were calculated with the quartet puzzling method. Branch lengths are proportional to the genetic divergence. The scale bar indicates the number of amino acid substitutions per site. Numbers at nodes represent percentages of bipartitions in intermediate trees that have been generated in 25 000 puzzlingsteps. The BEV prototype strains of the previous proposal by Dunne et al. (1974)Down are printed in bold and underlined. The proposed BEV clusters (species) and geno-/serotypes are indicated. Note that the nomenclature of the geno-/serotypes of this study differsfrom the previous proposal.

 
For the German BEV isolate D 14/3/96, phylogenetic analyses yielded incongruent results. Whereas phylogenetic analysis using the complete P1 region indicates a relation to geno-/serotype 1 of the BEV-A cluster, similar analyses of the individual capsid proteins revealed a more complex picture. Comparison of the 1A protein, which has a highly conserved sequence, allowed clustering into BEV-A but not further subdivision into geno-/serotypes (Fig. 3Down). Comparison of 1B and 1C indicates a similarity with geno-/serotype 2 of the BEV-A cluster. However, in the 1D comparison D 14/3/96 clusters with geno-/serotype 1 (Fig. 3Down). These findings are an indication of interserotypic recombination in the evolution of the BEVs. Another example of interserotypic recombination is PS 87/Belfast, which is closely related to BEV-B serotype 2 in the 1A region but is more related to BEV-261 (M2 and RM2) in the capsid proteins 1B–1D. An example of intraserotypic recombination is shown by LC-R4, which groups with Vir 404/03 in the 1B phylogeny and with VG(5)27 in the 1C and 1D phylogenies (Fig. 3Down). Besides these examples, incongruities between the 5'-NTR and 1A phylogenies on the one hand, and between 1D and 3D phylogenies on the other hand, indicate more frequent recombination events outside the capsid-encoding regions.


Figure 3
View larger version (32K):
[in this window]
[in a new window]
 
Fig. 3. Recombination within the BEVs. Unrooted neighbour-joining trees of the 5'-NTR, each of the four capsid proteins and the 3D polymerase. Using porcine enteroviruses 9 and 10 as outgroups, the 5'-NTR, 1A, 1B, 1C and 1D protein of 25 enteroviruses and the 3D protein of 24 enteroviruses were compared. The isolate D 14/3/96 is boxed. Nucleotide and amino acid sequences were aligned manually. For tree construction, programs DNAML and PROTML/Neighbour were employed (Felsenstein, 1995Down; Adachi & Hasegawa, 1992Down). The BEV prototype strains of the previous proposal by Dunne et al. (1974)Down are printed in bold and underlined. The proposed BEV clusters (species) and geno-/serotypes are indicated where applicable.

 
Eight of eleven HEV-C prototype strains (CVA11, 13, 15, 17, 18, 20, 21 and 24) clustered together with the poliovirus sequences (Figs 2 and 3UpUp), indicating the need to revise the present enterovirus taxonomy (see also Brown et al., 2003Down).

The phylogenetic analysis of the 3D polymerase included 72 rhinovirus and enterovirus sequences. The phylogenetic tree based on these sequences reveals genetic clusters representing species, whereas individual serotypes cannot be differentiated (Fig. 4Down). The BEV strains fall into two clusters correlating with clusters A and B of the P1 analysis. Again, CVA1, CVA19 and CVA22 separate from the remaining CVA strains (CVA11, 13, 15, 17, 18, 20, 21 and 24) and the polioviruses.


Figure 4
View larger version (21K):
[in this window]
[in a new window]
 
Fig. 4. Phylogenetic relationships of the BEVs to the entero- and rhinoviruses. Unrooted neighbour-joining tree of the 3D polymerase. Amino acid sequences were aligned with the CLUSTAL W program. For the details of tree construction, refer to the legend of Fig. 2Up. The BEV prototype strains of the previous proposal by Dunne et al. (1974)Down are printed in bold and underlined. The proposed BEV clusters (species) are indicated.

 
Non-translated regions (NTRs)
The open reading frame of the BEV genome is flanked on either side by NTRs. The length of the 5'-NTR ranges from 812 to 822 nt. The overall folding of the BEV 5'-NTR is similar to that of the human enteroviruses but is characterized by sequences that have the potential to form three unique RNA structures, i.e. an additional cloverleaf structure (designated domain I*, Fig. 5Down) separated from the 5'-cloverleaf by a stem–loop, a small domain III and a less conserved domain VI (Zell & Stelzner, 1997Down). These structures are conserved among all sequenced BEV genomes. Phylogenetic analysis of 56 enteroviral 5'-NTR sequences revealed five clades correlating with five folding patterns of entero- and rhinoviral 5'-NTRs (Fig. 5Down and Supplementary Figure in JGV Online). However, three of five groups (i.e. human enteroviruses, human rhinoviruses and bovine enteroviruses) contain more than one species.


Figure 5
View larger version (22K):
[in this window]
[in a new window]
 
Fig. 5. Analysis of enteroviral and rhinoviral 5'-NTR. Unrooted neighbour-joining tree of 5'-NTR sequences. Nucleotide sequences were manually aligned. Sequences representing the second BEV cloverleaf were deleted. For the details of tree construction, refer to the legend of Fig. 2Up. Circles indicate sequence clusters with the same RNA-folding pattern, which are schematically illustrated. Significant differences of the folding patterns are boxed.

 
The 3'-NTR of the BEV-A cluster is 2–4 nt shorter than that of the BEV-B cluster (71 nt vs 73–75 nt). The predicted secondary structures (Fig. 6Down, upper panel) reveal a poliovirus-like folding pattern consisting of two stem–loop structures (domains X and Y), which have the potential to form a pseudoknot-like element (PKLE). Comparison of the tertiary structure models, based on the predicted secondary structures, show that the helices X and K (PKLE) are almost identical, while the Y helices of both clusters differ in length, in sequence and in the central base–base mismatch (Fig. 6Down, lower panel).


Figure 6
View larger version (14K):
[in this window]
[in a new window]
 
Fig. 6. Comparison of the putative 3'-NTR secondary structures of both BEV clusters and corresponding tertiary structure predictions. The putative domains are designated X and Y starting from the 3'-end. Nucleotides involved in the formation of a PKLE (helix K) are underlined.

 

   DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
In an earlier attempt Dunne et al. (1974)Down suggested the existence of eight BEV serotypes. However, only seven representative viruses were deposited in the ATCC, while the eighth strain (designated PS 35) was subsequently found to be identical to PS 83, the proposed serotype 5 isolate. Later, this classification was revised following a proposal that BEVs fell into only two serotypes (Knowles & Barnett, 1985Down). Typing of the BEVs by serological means may be difficult due to significant cross-reaction of type-specific sera with other BEV types. After accumulation of BEV sequence data, concerns arose as to whether the concept of two serotypes within a single BEV species correctly reflects the genetic variability of these viruses. Comparison of the sequences from the capsid protein region revealed relatively low intertypic amino acid identities ranging from 72 to 76 % between the type 1 and type 2 strains [VG(5)27, BEV-165 vs PS 87, RM2] (Zell & Stelzner, 1997Down). In order to establish a molecular-based BEV classification, six of the former prototype strains collected in the USA between 1957 and 1962 as well as 11 German field isolates (collected between 1982 and 2003) were sequenced. These sequences were compared with available BEV sequences. The data revealed that the bovine enteroviruses compose two major clusters that are characterized by structural features and by levels of sequence identity. These clusters (designated BEV-A and -B) correlate with two ‘serotypes' as defined by Knowles & Barnett (1985)Down. However, there are strong arguments supporting the hypothesis that both clusters represent species rather than serotypes: (i) structural features (lengths of proteins 1B, 1C, 1D, 2C and length of 3'-NTR) are correlated with the clusters BEV-A and -B. (ii) The distribution pattern of amino acid identity scores obtained after pairwise comparisons of capsid protein sequences consists of three peaks (Fig. 1Up). (iii) Phylogenetic analyses of the 3D polymerase reveals two clades – an observation that is confirmed by the discontinuous distribution pattern of the 3D amino acid identity scores (Figs 1 and 4UpUp).

Within each cluster, further subgrouping was observed in pairwise sequence comparisons and in the phylogenetic analyses. This subgrouping partly correlates to the serotyping previously proposed by Dunne et al. (1974)Down. While LC-R4, BEV-261 (M2 or RM2), PS 87 and PS 89 represent four discrete geno-/serotypes, BEV-165 (M4), PS 42 and PS 83 are closely related and should be considered as strains of another fifth geno-/serotype. Moreover, PS 42 and PS 83 are almost identical, even though they were deposited in the ATCC as discrete serotypes.

Phylogenetic analyses of the 5'-NTR, single capsid proteins and the 3D polymerase provided evidence for interserotypic and intraserotypic recombination (Fig. 3Up). Intertypic recombination in poliovirus (Cammack et al., 1988Down; Lipskaya et al., 1991Down) and echovirus (Andersson et al., 2002Down) has been described previously. In general, incongruence between phylogenies of different genome regions are considered indicative of recombination events in enteroviruses (Santti et al., 1999Down; Lindberg et al., 2003Down).

It is a generally accepted concept that picornavirus serotypes are molecularly defined by the diversity of the capsid proteins, whereas the less diverse non-structural protein regions define an enterovirus species. Accordingly, amino acid identities of the BEV 1D protein ranged from 50 to 55 % for heterologous species, 70 to 85 % for heterologous serotypes/homologous species and were greater than 90 % for homologous serotypes (compare Fig. 1Up). For the 3D polymerase, the observed amino acid identities were greater than 95 % of heterologous serotypes/homologous species (Fig. 1Up). In this context, difficulties in typing the isolate D 14/3/96 may be explained by an interserotypic recombination event in the evolution of this virus. However, proof of this hypothesis requires the analysis of numerous virus isolates collected in close temporal and spatial proximity. With the exception of D 14/1/96, such virus isolates are not available. However, it is an interesting observation that D 14/1/96 and D 14/3/96 are closely related in their 5'-NTR and 3D gene region (Fig. 3Up). The capsid proteins of both isolates are considerably divergent. This suggests multiple recombination events in the evolution of these epidemiologically related viruses. A prerequisite of enterovirus recombination is the double/multiple infection of an individual host, which is common during vaccination with the live-attenuated polio vaccine but less frequent with non-polio enteroviruses.

Recently, Goens et al. (2004)Down published the genome sequences of PS 87/Maryland and Wye-3A. Both sequences cluster with PS 89 (Figs 2 and 3UpUp). Likewise, PS 87/Pirbright and PS 89 are identical (N. J. Knowles, personal communication). Instead, the sequence of PS 87/Belfast is identical to the previously published partial PS 87 sequence (GenBank accession no. X79368; McNally et al., 1994Down). It is noteworthy that the PS 87/Belfast and PS 89 strains used in the present study were directly received from the ATCC. Thus, one cannot exclude that, during the passage history of PS 87/Maryland and PS 87/Pirbright, this strain was mixed up with PS 89 or that virus stocks were mixtures of both viruses. (Indeed it is possible that some of the confusion presented by previous BEV serological classification studies have arisen by some virus isolates being mixtures of bovine enteroviruses).

Further analysis of the BEV genomes confirmed previous findings concerning the secondary structures of the 5'-NTR. The BEVs possess a second cloverleaf and characteristic RNA structures of the internal ribosome entry site as previously described (Zell & Stelzner, 1997Down; Zell et al., 1999Down). Krumbholz et al. (2002)Down have recently described five groups of entero- and rhinoviruses based on their predicted 5'-NTR folding patterns. Phylogenetic analysis of enteroviral and rhinoviral 5'-NTRs based on a representative number of sequences confirmed this observation (Fig. 5Up). Although characteristic for the BEVs, the 5'-NTR sequences are not a distinctive feature for species and serotype demarcation.

The 3'-NTR of the BEVs is poliovirus-like as suggested from the secondary structure predictions. It consists of two stem–loop structures (domains X and Y) allowing the formation of a PKLE. Tertiary structure models of representative BEV strains suggest that base stacking of helix X protrudes into the loop of the X domain, which forms helix K by base-pairing with the loop of the Y domain. This part of the 3'-NTR is almost identical in all sequenced BEVs (Fig. 6Up) and the porcine enteroviruses (data not shown), while the Y helix differs in these viruses.

Although enterovirus taxonomy was recently revised (King et al., 2000Down; Stanway et al., 2005Down), the present classification of the species Bovine enterovirus, Human enterovirus C and Poliovirus still appears to be inadequate. The sequencing data presented in this study support a revision of the present enterovirus taxonomy. According to our results, the bovine enteroviruses should be classified into two species (BEV-A and -B) containing at least two and three geno-/serotypes, respectively. Sequencing of further bovine, caprine and ovine enterovirus isolates may lead to the identification of additional geno-/serotypes. Like the human enteroviruses, most of the BEV isolates may unambiguously be typed by 1B, 1C or 1D gene regions, with the exception of isolate D 14/3/96. A more comprehensive set of sequence data may help to define subgenomic regions suited for the genotyping of BEV isolates.


   ACKNOWLEDGEMENTS
 
The excellent technical assistance of Veronika Güntzschel, Melanie Harder, Katja Schneider and Sabine Wachsmuth is appreciated. We are indebted to Nick Knowles (Institute for Animal Health, Pirbright, UK) for stimulating discussions, exchange of unpublished data and the communication of useful PCR primers. We thank Rüdiger Wurm for the gift of the isolate Vir 404/03.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Adachi, J. & Hasegawa, M. (1992). Amino acid substitution of proteins coded for in mitochondrial DNA during mammalian evolution. Jpn J Genet 67, 187–197.[CrossRef][Medline]

Andersson, P., Edman, K. & Lindberg, A. M. (2002). Molecular analysis of the echovirus 18 prototype: evidence of interserotypic recombination with echovirus 9. Virus Res 85, 71–83.[CrossRef][Medline]

Barya, M. A., Moll, T. & Mattson, D. E. (1967). Antigenic analysis of bovine enteroviruses through studies of kinetics of neutralization. Am J Vet Res 28, 1283–1294.[Medline]

Brown, B., Oberste, M. S., Maher, K. & Pallansch, M. A. (2003). Complete genomic sequencing shows that polioviruses and members of human enterovirus species C are closely related in the noncapsid coding region. J Virol 77, 8973–8984.[Abstract/Free Full Text]

Cammack, N., Phillips, A., Dunn, G., Patel, V. & Minor, P. D. (1988). Intertypic genomic rearrangements of poliovirus strains in vaccines. Virology 167, 507–514.[Medline]

Dunne, H. W., Huang, C. M. & Lin, W. J. (1974). Bovine enteroviruses in the calf: an attempt at serologic, biologic, and pathologic classification. J Am Vet Med Assoc 164, 290–294.[Medline]

Earle, J. A. P., Skuce, R. A., Fleming, C. S., Hoey, E. M. & Martin, S. J. (1988). The complete nucleotide sequence of a bovine enterovirus. J Gen Virol 69, 253–263.[Abstract/Free Full Text]

Felsenstein, J. (1995). PHYLIP: phylogeny inference package, version 3.57c. University of Washington, Seattle, WA.

Goens, S. D., Botero, S., Zemla, A., Zhou, C. E. & Perdue, M. L. (2004). Bovine enterovirus 2: complete genomic sequence and molecular modelling of a reference strain and a wild-type isolate from endemically infected US cattle. J Gen Virol 85, 3195–3203.[Abstract/Free Full Text]

Huck, R. A. & Cartwright, S. F. (1964). Isolation and classification of viruses from cattle during outbreaks of mucosal or respiratory disease and from herds with reproductive disorders. J Comp Pathol 74, 346–365.[Medline]

Jones, D. T., Taylor, W. R. & Thornton, J. M. (1992). The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci 8, 275–282.[Abstract/Free Full Text]

King, A. M. Q., Brown, F., Christian, P. & 8 other authors (2000). Picornaviridae. In Virus Taxonomy. Seventh Report of the International Committee for the Taxonomy of Viruses, 7th edn, pp. 657–678. Edited by M. H. V. Van Regenmortel, C. M. Fauquet, D. H. L. Bishop, C. H. Calisher, E. B. Carsten, M. K. Estes, S. M. Lemon, J. Maniloff, M. A. Mayo, D. J. McGeoch, C. R. Pringle & R. B. Wickner. New York, San Diego: Academic Press.

Knowles, N. J. & Barnett, I. T. (1985). A serological classification of bovine enteroviruses. Arch Virol 83, 141–155.[CrossRef][Medline]

Krumbholz, A., Dauber, M., Henke, A., Birch-Hirschfeld, E., Knowles, N. J., Stelzner, A. & Zell, R. (2002). Sequencing of porcine enterovirus groups II and III reveals unique features of both virus groups. J Virol 76, 5813–5821.[Abstract/Free Full Text]

Kunin, C. M. & Minuse, E. (1958). The isolation in tissue culture, chick embryo and suckling mice of filtrable agents from healthy dairy cattle. J Immunol 80, 1–11.

Lindberg, A. M., Andersson, P., Savolainen, C., Mulders, M. N. & Hovi, T. (2003). Evolution of the genome of Human enterovirus B: incongruence between phylogenies of the VP1 and 3CD regions indicates frequent recombination within the species. J Gen Virol 84, 1223–1235.[Abstract/Free Full Text]

Lipskaya, G. Y., Muzychenko, A. R., Kutitova, O. K., Maslova, S. V., Equestre, M., Drozdov, S. G., Perez Bercoff, R. P. & Agol, V. I. (1991). Frequent isolation of intertypic poliovirus recombinants with serotype 2 specificity from vaccine-associated polio cases. J Med Virol 35, 290–296.[Medline]

McCarthy, F. M., Smith, G. A. & Mattick, J. S. (1999). Molecular characterisation of Australian bovine enteroviruses. Vet Microbiol 68, 71–81.[CrossRef][Medline]

McFerran, J. B. (1958). ECBO viruses of cattle. Vet Rec 70, 999.

McNally, R. M., Earle, J. A., McIlhatton, M., Hoey, E. M. & Martin, S. J. (1994). The nucleotide sequence of the 5' non-coding and capsid coding genome regions of two bovine enterovirus strains. Arch Virol 139, 287–299.[CrossRef][Medline]

Moll, T. & Davis, A. D. (1959). Isolation and characterization of cytopathogenic enteroviruses from cattle with respiratory disease. Am J Vet Res 20, 27–32.[Medline]

Moll, T. & Ulrich, M. I. (1963). Biologic characteristics of certain bovine enteric viruses. Am J Vet Res 24, 545–550.[Medline]

Möller, W. & Amons, R. (1985). Phosphate-binding sequences in nucleotide-binding proteins. FEBS Lett 186, 1–7.[CrossRef][Medline]

Santti, J., Hyypiä, T., Kinnunen, L. & Salminen, M. (1999). Evidence of recombination among enteroviruses. J Virol 73, 8741–8749.[Abstract/Free Full Text]

Stanway, G., Brown, F., Christian, P. & 9 other authors (2005). Picornaviridae. In Virus Taxonomy, Eighth Report of the International Committee on Taxonomy of Viruses, 8th edn, pp. 757–778. Edited by C. M. Fauquet, M. A. Mayo, J. Maniloff, U. Desselberger & L. A. Ball. London: Elsevier/Academic Press.

Strimmer, K. & von Haeseler, A. (1996). Quartet puzzling: a quartet maximum likelihood method for reconstructing tree topologies. Mol Biol Evol 13, 964–969.

Strimmer, K., Goldman, N. & von Haeseler, A. (1997). Bayesian probabilities and quartet puzzling. Mol Biol Evol 14, 210–211.

Tamura, K. & Nei, M. (1993). Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10, 512–526.[Abstract]

Thompson, J. D., Higgins, D. G. & Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22, 4673–4680.[Abstract/Free Full Text]

Van Regenmortel, M. H. V., Bishop, D. H. L., Fauquet, C. M., Mayo, M. A., Maniloff, J. & Calisher, C. H. (1997). Guidelines to the demarcation of virus species. Arch Virol 142, 1505–1518.[Medline]

Zell, R. & Stelzner, A. (1997). Application of genome sequence information to the classification of bovine enteroviruses: the importance of 5'- and 3'-nontranslated regions. Virus Res 51, 213–229.[CrossRef][Medline]

Zell, R., Sidigi, K., Henke, A., Schmidt-Brauns, J., Hoey, E., Martin, S. & Stelzner, A. (1999). Functional features of the bovine enteroviral 5'-nontranslated region. J Gen Virol 80, 2299–2310.[Abstract/Free Full Text]

Zell, R., Dauber, M., Krumbholz, A., Henke, A., Birch-Hirschfeld, E., Stelzner, A., Prager, D. & Wurm, R. (2001). Porcine teschoviruses comprise at least eleven distinct serotypes: molecular and evolutionary aspects. J Virol 75, 1620–1631.[Abstract/Free Full Text]

Zhang, A. Q. & Burgess, G. W. (1986). Serotyping bovine enteroviruses in Queensland. In Proceedings of the 5th International conference on Livestock Production and Diseases in the Tropics, Kuala Lumpur, Malaysia.

Zuker, M., Mathews, D. H. & Turner, D. H. (1999). Algorithms and thermodynamics for RNA secondary structure prediction: a practical guide. In RNA Biochemistry and Biotechnology. NATO ASI Series, pp. 11–43. Edited by J. Barciszewski & B. F. C. Clark. Kluwer Academic Publishers.

Received 29 June 2005; accepted 12 October 2005.


This article has been cited by other articles:


Home page
J. Virol.Home page
P. Simmonds
Recombination and Selection in the Evolution of Picornaviruses and Other Mammalian Positive-Stranded RNA Viruses
J. Virol., November 15, 2006; 80(22): 11124 - 11140.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zell, R.
Right arrow Articles by Wutzler, P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zell, R.
Right arrow Articles by Wutzler, P.
Agricola
Right arrow Articles by Zell, R.
Right arrow Articles by Wutzler, P.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS