|
|
||||||||
1 Rabies Union, Centers for Disease Control and Prevention, 1600 Clifton Road, MS G-33, Atlanta, GA 30333, USA
2 Laboratory for Clinical and Molecular Virology, The University of Edinburgh, Summerhall, Edinburgh EH9 1QH, UK
Correspondence
I. V. Kuzmin
ibk3{at}cdc.gov
| ABSTRACT |
|---|
|
|
|---|
The GenBank/EMBL/DDBJ accession numbers of the rhabdovirus partial N gene sequences obtained in this study are DQ457098DQ457104, as outlined in Table 1
.
| INTRODUCTION |
|---|
|
|
|---|
Early taxonomy of the Rhabdoviridae, as well as for other viruses, was based on virion morphology and serological cross-reactivity. Indeed, many as-yet unclassified rhabdoviruses were believed to exhibit cross-reactivity with members of the genus Lyssavirus (Calisher et al., 1989
). When the terms rabies-related viruses' and rabies serogroup were introduced (Shope et al., 1970
), these included Rabies virus (RABV), Mokola virus (MOKV) and Lagos bat virus (LBV). Shortly thereafter, new members of the serogroup were described, such as Duvenhage virus (DUVV) [including at that time European bat lyssavirus 1 (EBLV-1)], Kotonkon virus (KOTV) and Obodhiang virus (OBOV) (Meredith et al., 1971
; Kemp et al., 1973
; Tignor et al., 1977
). The latter two viruses, isolated from mosquitoes, demonstrated limited cross-reactivity with MOKV, but not with other members of the serogroup. Additionally, MOKV can be propagated in insect cell lines (Buckley, 1975
), although no insect-derived isolates have been obtained. As was speculated, the proposed evolutionary pathway for the rabies-related viruses included OBOV and KOTV as progenitors, through MOKV, to other highly neurotropic mammalian viruses, such as LBV, DUVV and RABV (Shope, 1982
). Later, other rhabdoviruses, isolated from both arthropods and vertebrates, predominantly in Africa, were described as related serologically to those listed above (Calisher et al., 1989
). Currently, the evolutionary history of the Rhabdoviridae remains unclear.
Gene sequencing and phylogenetic relationships replaced the initial basis of serological and antigenic relationships for virus taxonomy. European bat lyssaviruses were differentiated from DUVV and separated into two species (EBLV-1 and EBLV-2) (King et al., 1990
). A number of new lyssaviruses were isolated in different parts of the world (Gould et al., 1998
; Kuzmin et al., 2003
, 2005
). Phylogenetically, lyssaviruses segregate into a single monophyletic clade, distant from other rhabdovirus genera and unclassified rhabdoviruses (Kuzmin & Rupprecht, 2005
; Springfeld et al., 2005
). A phylogenetic study based on limited polymerase (L) gene sequences was recently described for 38 rhabdoviruses (Bourhy et al., 2005
). Although no new relatives of lyssaviruses were detected in that study, only KOTV from the number of serologically rabies-related viruses was analysed. No other information is available to date on the molecular sequences of other rhabdoviruses that demonstrated serological cross-reactivity with lyssaviruses.
In the present study, we generated partial nucleoprotein (N) gene sequences of OBOV, KOTV, Sandjimba virus (SJAV), Kolongo virus (KOLV), Mount Elgon bat virus (MEBV), Kern Canyon virus (KCV) and Rochambeau virus (RBUV) and attempted to determine their phylogenetic relationships within the family Rhabdoviridae. Estimates of non-synonymous/synonymous (dN/dS) substitution ratios and detection of positively and negatively selected sites were used to assess evolutionary selection pressures along a fragment of the N gene for the three genera. The objective of this study was to provide further insight into the evolutionary history of the Rhabdoviridae.
| METHODS |
|---|
|
|
|---|
|
Primary assembly, alignment and consensus sequence generation, as well as DNA translation, were performed in BioEdit software (Hall, 1999
). Multiple datasets of both nucleotide and amino acid sequences were subjected to phylogenetic analysis. Initially, the fragment of the N gene obtained was aligned and compared with whole-length N gene sequences of all rhabdovirus representatives available in GenBank (Table 2
). To avoid the influence of multiple alignment gaps, only a set of the most closely related sequences was aligned with the newly generated sequences. This alignment was truncated from both the 3' and 5' ends to avoid gaps at the ends of the alignment and to ensure maximum compatibility of the sequences compared.
|
For positive selection analysis, 241 codons of the N gene were selected for a number of representative rhabdoviruses, producing alignments without gaps (positions 3861108 relative to the PV genome): lyssaviruses, 11 taxa (seven recognized species, ARAV, KHUV, IRKV and WCBV); vesiculoviruses, 10 taxa [ISFV, SVCV, PIRYV, CHPV, VSINV (three lineages) and VSNJV (three lineages)]; and ephemeroviruses, six taxa (BEFV, ARV, FLAV, OBOV, KOTV and RBUV). All three datasets were analysed using the Datamonkey online positive selection interface (Kosakovsky Pond & Frost, 2005a
). For each alignment, an NJ tree and ML model were estimated and used as the input for the single likelihood ancestor counting (SLAC) model of positive selection (Kosakovsky Pond & Frost, 2005b
). Estimation of the ratio of non-synonymous substitutions per non-synonymous site (dN) to synonymous substitutions per synonymous site (dS) provides an indication of selection pressures acting on a gene. The majority of codons within a viral gene are likely to be under negative or purifying selection (dN/dS
1). Although average dN/dS values for genes are uninformative for potential positively selected codons, they do provide a means of assessing adaptive evolutionary pressures at the gene level. A value of P of 0.05 was used for inference of positive and negative selection for individual codons. A global estimate of dN/dS was generated with 95 % profile likelihood intervals by SLAC for each dataset.
| RESULTS |
|---|
|
|
|---|
Phylogenetic analysis recognized all seven newly sequenced viruses as members of the Dimarhabdovirus supergroup (Bourhy et al., 2005
), without particular relatedness to the genus Lyssavirus (Fig. 1
). OBOV, KOTV and RBUV were placed in the genus Ephemerovirus, independently of the phylogenetic method (for both nucleotide and amino acid sequences). Support for the joining of RBUV was high (7697 % bootstrap values) only for the nucleotide NJ tree, whereas for the amino acid NJ trees and for all ML and MP trees, bootstrap support was limited (4167 %). In all amino acid trees, FLAV joins to the ephemeroviruses as well. However, in these trees, the position of FLAV was not consistently supported by high bootstrap values (3267 % depending on the method used). Moreover, in the nucleotide trees, FLAV was always moved to the cluster of KOLV, SJAV and TRV (yet ancestrally and with limited bootstrap support; not shown). KOLV and SJAV formed a monophyletic clade, separated from the other genera of rhabdoviruses. TRV was joined to the same clade in all phylogenetic constructions, quite distantly but with consistent bootstrap support. MEBV and KCV reliably formed a monophyletic clade with OITAV, distinct from other rhabdovirus genera and unclassified representatives.
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
None of the seven viruses sequenced in this study were placed phylogenetically within the genus Lyssavirus, as has been suggested previously by serological cross-reactivity (Shope, 1982
; Calisher et al., 1989
). Instead, they were placed within the Dimarhabdovirus supergroup defined in a recent study (Bourhy et al., 2005
). Representatives of this group have been isolated from arthropods (predominantly dipterans) and mammals and are proposed to be vector-borne. Thus, there is currently no documented support for the hypothesis that lyssaviruses evolved from rhabdoviruses of plants and arthropods.
OBOV and KOTV were placed within the genus Ephemerovirus. Ephemeroviruses are distributed widely in the Old World tropics and subtropics, including Australia, Africa, Asia and the Middle East. Mosquitoes are the only insects from which ephemeroviruses have been isolated. Wind-borne spread across the tropics and subtropics has been suggested for viruses of this group, and cattle translocation between continents has been suggested to explain their broad distribution in the Eastern hemisphere (Nandi & Negi, 1999
; Walker, 2005
). The identity between Australian ephemeroviruses ARV and BEFV was less than the identity between ARV and OBOV, which was isolated in Africa. Too few ephemerovirus isolates are available for accurate calculations of the time of divergence and, at present, such global evolutionary questions are unanswerable. The presence in this cluster of the more distantly related KOTV and RBUV suggests that intercontinental translocations of ephemeroviruses or their progenitors may have occurred long ago. In taxonomic terms, RBUV and FLAV are not placed with confidence into the genus Ephemerovirus, but their genetic relatedness is clear. Further nucleotide sequencing and phylogenetic analysis may allow this relationship to be investigated more substantially.
KOLV and SJAV form a monophyletic clade which tentatively should be considered as a new genus of the Rhabdoviridae. Since these viruses were isolated from birds, there is no reason to suggest that their distribution is limited to central Africa. Possibly, these viruses are distributed widely by migratory birds, as occurs with other avian viruses (e.g. West Nile virus). TRV demonstrates distant relatedness to this group. TRV was isolated from a tree shrew (Tupaia belangeri) imported from Thailand and has been shown to replicate only in tree shrew cells (Springfeld et al., 2005
). These mammals are indigenous to tropical Asia, but not to Africa. Taken together, these data suggest that representatives of this group could be more numerous and circulate throughout different continents among different hosts. Further surveillance is needed to evaluate the phylogenetic and epidemiological patterns of these viruses.
MEBV and KCV, together with the recently described OITAV (Iwasaki et al., 2004
), constituted another monophyletic clade at the genus level. Remarkably, all three viruses were isolated from insectivorous bats in distant regions of the world: Japan, central Africa and western North America. If these viruses are associated specifically with bats, it is unlikely that they are transmitted by broad-ranging sanguivorous insects, such as mosquitoes. More likely, they may be transmitted by some specific bat ectoparasites (which still may be dipterans, such as nycteribiids) or they may circulate directly among bats without the requirement for the participation of arthropods.
The functional roles of conserved residues and motifs along the deduced N protein (Fig. 2
) are unclear. They do not overlap with any antigenic epitopes described for RABV (Dietzschold et al., 1987
; Ertl et al., 1991
; Fu et al., 1994
; da Cruz et al., 2001
). They may be involved in RNA binding or conformational patterns, as was shown for proline, glycine and alanine residues in the rhabdovirus N protein (Chou & Fasman, 1978
; Wang et al., 1995
).
Lyssaviruses constitute the most homogeneous clade within the phylogenetic tree of the Rhabdoviridae. Some authors have suggested that this may be a consequence of their youngest evolutionary age (Bourhy et al., 2005
). If so, this would imply that the predominance of purifying selection along the fragment of the N gene tested is due to very rapid (and successful) adaptation to mammalian transmission and constraints placed upon lyssaviruses by their unique pathobiology. Such a situation may not to be true of the ephemeroviruses and vesiculoviruses. Estimation of divergence times for rhabdovirus genera is highly problematic. Although it is possible to estimate substitution rates for individual species, e.g. RABV (Hughes et al., 2005
) and EBLV-1 (Davis et al., 2005
), application of these rates to longer trees (i.e. above the species level) is confounded by substitution saturation and selection pressures that may have occurred during adaptation of progenitor viruses for emerging species.
However, our estimations clearly demonstrate that lyssaviruses are subject to strong constraints against amino acid substitutions. Positive selection analysis using SLAC detected no codons under positive selection for any dataset. Considering the conservative nature of the SLAC method (Kosakovsky Pond & Frost, 2005b
) and previous findings of largely purifying selection in lyssavirus N genes (Holmes et al., 2002
; Kuzmin et al., 2004
; Hughes et al., 2005
; Davis et al., 2005
), this is not surprising. Such constraints appear to be much greater than those observed among ephemeroviruses and vesiculoviruses.
In general, the greatest diversity in the Rhabdoviridae is seen among the insect-borne plant viruses (genera Cytorhabdovirus and Nucleorhabdovirus), followed by Novirhabdovirus, isolated from fish and other aquatic animals (including invertebrates), the dimarhabdoviruses (able to replicate in vertebrate and invertebrate hosts) and finally the lyssaviruses, which replicate exclusively in mammals (Johnson et al., 1999
; Tordo et al., 2004
; Fu, 2005
; Bourhy et al., 2005
).
No particular phylogenetic relatedness was demonstrated between lyssaviruses and the other rhabdoviruses sequenced to date. Limited serological cross-reactivity detected in earlier studies between lyssaviruses and the seven viruses which were sequenced during the current work (Shope, 1982
; Calisher et al., 1989
) could be attributed to common antigenic determinants, since they are all members of one family. We are clearly missing some important chains in the evolutionary pathway of the Rhabdoviridae. With more than 40 unassigned animal rhabdoviruses in the list of the International Committee on Taxonomy of Viruses (Tordo et al., 2004
), additional sequencing is essential to clarify phylogenetic relationships. Such studies with known isolates, together with enhanced field surveillance for new pathogens, should provide greater insight into the biology and evolution of this family.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Buckley, S. M. (1975). Arbovirus infection of vertebrate and insect cell cultures, with special emphasis on Mokola, Obodhiang, and kotonkan viruses of the rabies serogroup. Ann N Y Acad Sci 266, 241250.[Medline]
Calisher, C. H., Karabatsos, N., Zeller, H., Digoutte, J. P., Tesh, R. B., Shope, R. E., Travassos da Rosa, A. P. & St George, T. D. (1989). Antigenic relationships among rhabdoviruses from vertebrates and hematophagous arthropods. Intervirology 30, 241257.[Medline]
Chou, P. Y. & Fasman, D. G. (1978). Prediction of the secondary structure of proteins from their amino acid sequences. Adv Enzymol Relat Areas Mol Biol 47, 45148.[Medline]
da Cruz, F. W., McBride, A. J. A., Conceição, F. R., Dale, J. W., McFadden, J. & Dellagostin, O. A. (2001). Expression of the B-cell and T-cell epitopes of the rabies virus nucleoprotein in Mycobacterium bovis BCG and induction of an humoral response in mice. Vaccine 20, 731736.[CrossRef][Medline]
Davis, P. L., Holmes, C. E., Larrous, F., Van der Poel, W. H. M., Tjørnehøj, K., Alonso, W. J. & Bourhy, H. (2005). Phylogeography, population dynamics, and molecular evolution of European bat lyssaviruses. J Virol 79, 1048710497.
Dietzschold, B., Lafon, M., Wang, H., Otvos, L., Jr, Celis, E., Wunner, W. H. & Koprowski, H. (1987). Localization and immunological characterization of antigenic domains of the rabies virus internal N and NS proteins. Virus Res 8, 103125.[Medline]
Ertl, H., Dietzchold, B. & Otvos, L., Jr (1991). T helper cell epitope of rabies virus nucleoprotein defined by tri- and tetrapeptides. Eur J Immunol 21, 110.[Medline]
Felsenstein, J. (1993). PHYLIP phylogeny inference package (version 3.5c). Distributed by the author. Department of Genome Sciences, University of Washington, Seattle, USA.
Fu, Z. F. (2005). Genetic comparisons of rhabdoviruses from animals and plants. Curr Top Microbiol Immunol 292, 124.[Medline]
Fu, Z. F., Wunner, W. H. & Dietzschold, B. (1994). Immunoprotection by rabies virus nucleoprotein. In Lyssaviruses, pp. 207218. Edited by C. E. Rupprecht, B. Dietzchold & H. Koprowski. Berlin: Springer.
Gould, A. R., Hyatt, A. D., Lunt, R., Kattenbelt, J. A., Hengstberger, S. & Blacksell, S. D. (1998). Characterisation of a novel lyssavirus isolated from pteropid bats in Australia. Virus Res 54, 165187.[CrossRef][Medline]
Hall, T. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 41, 9598.
Holmes, E. C., Woelk, C. H., Kassis, R. & Bourhy, H. (2002). Genetic constraints and the adaptive evolution of rabies virus in nature. Virology 292, 247257.[CrossRef][Medline]
Hughes, G. J., Orciari, L. A. & Rupprecht, C. E. (2005). Evolutionary timescale of rabies virus adaptation to North American bats inferred from the substitution rate of the nucleoprotein gene. J Gen Virol 86, 14671474.
Iwasaki, T., Inoue, S., Tanaka, K. & 8 other authors (2004). Characterization of Oita virus 296/1972 of Rhabdoviridae isolated from a horseshoe bat bearing characteristics of both lyssavirus and vesiculovirus. Arch Virol 149, 11391154.[CrossRef][Medline]
Jeanmougin, F., Thompson, J. D., Gouy, M., Higgins, D. G. & Gibson, T. J. (1998). Multiple sequence alignment with Clustal X. Trends Biochem Sci 23, 403405.[CrossRef][Medline]
Johnson, M. C., Maxwell, J. M., Loh, P. C. & Leong, J. C. (1999). Molecular characterization of the glycoproteins from two warm water rhabdoviruses: snakehead rhabdovirus (SHRV) and rhabdovirus of penaeid shrimp (RPS)/spring viremia of carp virus (SVCV). Virus Res 64, 95106.[CrossRef][Medline]
Kemp, G. E., Lee, V. H., Moore, D. L., Shope, R. E., Causey, O. R. & Murphy, F. A. (1973). Kotonkan, a new rhabdovirus related to Mokola virus of the rabies serogroup. Am J Epidemiol 98, 4349.
King, A., Davis, P. & Lawrie, A. (1990). The rabies viruses of bats. Vet Microbiol 23, 165174.[CrossRef][Medline]
Koprowski, H. (1996). The mouse inoculation test. In Laboratory Techniques in Rabies, 4th edn, pp. 8086. Edited by F.-X. Meslin, M. M. Kaplan & H. Koprowski. Geneva: World Health Organization.
Kosakovsky Pond, S. L. & Frost, S. D. W. (2005a). Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics 21, 25312533.
Kosakovsky Pond, S. L. & Frost, S. D. W. (2005b). Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol 22, 12081222.
Kumar, S., Tamura, K., Jakobsen, I.-B. & Nei, M. (2001). MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17, 12441245.
Kuzmin, I. V. & Rupprecht, C. E. (2005). Lyssaviruses. In Encyclopedia of Medical Genomics and Proteomics, online edition. Edited by J. Fuchs & M. Podda. New York: Marcel Dekker. http://www.dekker.com/sdek/issues
content=t713172965
Kuzmin, I. V., Orciari, L. A., Arai, Y. T., Smith, J. S., Hanlon, C. A., Kameoka, Y. & Rupprecht, C. E. (2003). Bat lyssaviruses (Aravan and Khujand) from Central Asia: phylogenetic relationships according to N, P and G gene sequences. Virus Res 97, 6579.[CrossRef][Medline]
Kuzmin, I. V., Botvinkin, A. D., McElhinney, L. M., Smith, J. S., Orciari, L. A., Hughes, G. J., Fooks, A. R. & Rupprecht, C. E. (2004). Molecular epidemiology of terrestrial rabies in the former Soviet Union. J Wildl Dis 40, 617631.
Kuzmin, I. V., Hughes, G. J., Botvinkin, A. D., Orciari, L. A. & Rupprecht, C. E. (2005). Phylogenetic relationships of Irkut and West Caucasian bat viruses within the Lyssavirus genus and suggested quantitative criteria based on the N gene sequence for lyssavirus genotype definition. Virus Res 111, 2843.[CrossRef][Medline]
Meredith, C. D., Prossouw, A. R. & Van Praag Koch, H. (1971). An unusual case of human rabies thought to be of chiropteran origin. S Afr Med J 45, 767769.[Medline]
Nandi, S. & Negi, B. S. (1999). Bovine ephemeral fever: a review. Comp Immunol Microbiol Infect Dis 22, 8191.[CrossRef][Medline]
Sacramento, D., Bourhy, H. & Tordo, N. (1991). PCR technique as an alternative method for diagnosis and molecular epidemiology of rabies virus. Mol Cell Probes 5, 229240.[CrossRef][Medline]
Shope, R. E. (1982). Rabies-related viruses. Yale J Biol Med 55, 271275.[Medline]
Shope, R. E., Murphy, F. A., Harrison, A. K., Causey, O. R., Kemp, G. E., Simpson, D. I. & Moore, D. L. (1970). Two African viruses serologically and morphologically related to rabies virus. J Virol 6, 690692.
Springfeld, C., Darai, G. & Cattaneo, R. (2005). Characterization of the Tupaia rhabdovirus genome reveals a long open reading frame overlapping with P and a novel gene encoding a small hydrophobic protein. J Virol 79, 67816790.
Tignor, G. H., Murphy, F. A., Clark, H. F., Shope, R. E., Madore, P., Bauer, S. P., Buckley, S. M. & Meredith, C. D. (1977). Duvenhage virus: morphological, biochemical, histopathological and antigenic relationships to the rabies serogroup. J Gen Virol 37, 595611.
Tordo, N., Benmansour, A., Calisher, C. & 7 other authors (2004). Rhabdoviridae. In Virus Taxonomy, Eighth Report of the International Committee on Taxonomy of Viruses, pp. 623644. Edited by C. M. Fauquet, M. A. Mayo, J. Maniloff, U. Desselberger & L. A. Ball. London: Elsevier/Academic Press.
Walker, P. J. (2005). Bovine ephemeral fever in Australia and the world. Curr Top Microbiol Immunol 292, 5780.[Medline]
Wang, Y., Cowley, J. A. & Walker, P. J. (1995). Adelaide River virus nucleoprotein gene: analysis of phylogenetic relationships of ephemeroviruses and other rhabdoviruses. J Gen Virol 76, 995999.
Received 27 January 2006;
accepted 11 April 2006.
This article has been cited by other articles:
![]() |
B. Longdon, D. J. Obbard, and F. M. Jiggins Sigma viruses from three species of Drosophila form a major new clade in the rhabdovirus phylogeny Proc R Soc B, January 7, 2010; 277(1678): 35 - 44. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. V. Kuzmin, M. Niezgoda, R. Franka, B. Agwanda, W. Markotter, J. C. Beagley, O. Y. Urazova, R. F. Breiman, and C. E. Rupprecht Lagos Bat Virus in Kenya J. Clin. Microbiol., April 1, 2008; 46(4): 1451 - 1461. [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 | |