|
|
||||||||
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 |
|---|
|
|
|---|
The GenBank/EMBL/DDBJ accession numbers for the nucleotide sequences reported in the paper are DQ092769DQ092795.
| INTRODUCTION |
|---|
|
|
|---|
|
| METHODS |
|---|
|
|
|---|
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., 1994
). Neighbour-joining trees were calculated with the quartet puzzling method (Strimmer & von Haeseler, 1996
; Strimmer et al., 1997
) using the JTT substitution model for amino acid sequences (Jones et al., 1992
) and the TamuraNei model for nucleotide sequences (Tamura & Nei, 1993
). The reliability of the clustering was tested by 10 00025 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, 1995
) and PROTML program (Adachi & Hasegawa, 1992
).
RNA secondary predictions and free-energy calculations were performed with the MFOLD program (version 3.0; Zuker et al., 1999
).
Nucleotide sequence accession numbers.
The nucleotide sequence data reported in this paper were submitted to the GenBank nucleotide sequence database (accession nos DQ092769DQ092795, see also Table 2
). 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).
|
| RESULTS |
|---|
|
|
|---|
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 1A1D and 3D were calculated. In previous studies, this approach provided useful evidence for species distinction (Van Regenmortel et al., 1997
; Zell et al., 2001
). 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. 1
). These distribution patterns indicate the existence of heterogeneous geno-/serotypes and a higher order taxon, presumably a species.
|
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. 2
). 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)
as LC-R4, BEV-261 (M2), BEV-165 (M4), PS 87/Belfast and PS 89 were previously described as discrete serotypes (see Table 2
). 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. 2
) and are serologically almost indistinguishable (McNally et al., 1994
). Obviously, these sequences represent different passage histories of the same virus.
|
|
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. 4
). 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.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
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)
. 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. 3
). Intertypic recombination in poliovirus (Cammack et al., 1988
; Lipskaya et al., 1991
) and echovirus (Andersson et al., 2002
) has been described previously. In general, incongruence between phylogenies of different genome regions are considered indicative of recombination events in enteroviruses (Santti et al., 1999
; Lindberg et al., 2003
).
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. 1
). For the 3D polymerase, the observed amino acid identities were greater than 95 % of heterologous serotypes/homologous species (Fig. 1
). 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. 3
). 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)
published the genome sequences of PS 87/Maryland and Wye-3A. Both sequences cluster with PS 89 (Figs 2 and 3![]()
). 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., 1994
). 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, 1997
; Zell et al., 1999
). Krumbholz et al. (2002)
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. 5
). 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 stemloop 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. 6
) and the porcine enteroviruses (data not shown), while the Y helix differs in these viruses.
Although enterovirus taxonomy was recently revised (King et al., 2000
; Stanway et al., 2005
), 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 |
|---|
| REFERENCES |
|---|
|
|
|---|
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, 7183.[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, 12831294.[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, 89738984.
Cammack, N., Phillips, A., Dunn, G., Patel, V. & Minor, P. D. (1988). Intertypic genomic rearrangements of poliovirus strains in vaccines. Virology 167, 507514.[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, 290294.[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, 253263.
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, 31953203.
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, 346365.[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, 275282.
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. 657678. 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, 141155.[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, 58135821.
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, 111.
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, 12231235.
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, 290296.[Medline]
McCarthy, F. M., Smith, G. A. & Mattick, J. S. (1999). Molecular characterisation of Australian bovine enteroviruses. Vet Microbiol 68, 7181.[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, 287299.[CrossRef][Medline]
Moll, T. & Davis, A. D. (1959). Isolation and characterization of cytopathogenic enteroviruses from cattle with respiratory disease. Am J Vet Res 20, 2732.[Medline]
Moll, T. & Ulrich, M. I. (1963). Biologic characteristics of certain bovine enteric viruses. Am J Vet Res 24, 545550.[Medline]
Möller, W. & Amons, R. (1985). Phosphate-binding sequences in nucleotide-binding proteins. FEBS Lett 186, 17.[CrossRef][Medline]
Santti, J., Hyypiä, T., Kinnunen, L. & Salminen, M. (1999). Evidence of recombination among enteroviruses. J Virol 73, 87418749.
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. 757778. 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, 964969.
Strimmer, K., Goldman, N. & von Haeseler, A. (1997). Bayesian probabilities and quartet puzzling. Mol Biol Evol 14, 210211.
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, 512526.[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, 46734680.
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, 15051518.[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, 213229.[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, 22992310.
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, 16201631.
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. 1143. 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:
![]() |
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] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |