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


     


J Gen Virol 87 (2006), 2323-2331; DOI 10.1099/vir.0.81879-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
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 Kuzmin, I. V.
Right arrow Articles by Rupprecht, C. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kuzmin, I. V.
Right arrow Articles by Rupprecht, C. E.
Agricola
Right arrow Articles by Kuzmin, I. V.
Right arrow Articles by Rupprecht, C. E.
© 2006 Society for General Microbiology

Phylogenetic relationships of seven previously unclassified viruses within the family Rhabdoviridae using partial nucleoprotein gene sequences

I. V. Kuzmin1, G. J. Hughes2 and C. E. Rupprecht1

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
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Partial nucleoprotein (N) gene sequences of the rhabdoviruses Obodhiang (OBOV), Kotonkon (KOTV), Rochambeau (RBUV), Kern canyon (KCV), Mount Elgon bat (MEBV), Kolongo (KOLV) and Sandjimba (SJAV) were generated and their phylogenetic positions within the family Rhabdoviridae were determined. Both OBOV and KOTV were placed within the genus Ephemerovirus. RBUV was joined to the same cluster, but more distantly. MEBV and KCV were grouped into a monophyletic cluster (putative genus) with Oita virus (OITAV). These three viruses, originating from different regions of the world, were all isolated from insectivorous bats and may be specific for these mammals. African avian viruses KOLV and SJAV were joined to each other and formed another clade at the genus level. Further, they were grouped with the recently characterized rhabdovirus Tupaia virus (TRV). Although the genetic distance was great, the grouping was supported by consistent bootstrap values. This observation suggests that viruses of this group may be distributed widely in the Old World. Non-synonymous/synonymous substitution ratio estimations (dN/dS) using a partial N gene fragment (241 codons) for the three rhabdovirus genera revealed contrasting patterns of evolution, where dN/dS values follow the pattern Ephemerovirus > Vesiculovirus > Lyssavirus. The magnitude of this ratio corresponds well with the number of negatively selected codons. The accumulation of dS appears evenly distributed along the gene fragment for all three genera. These estimations demonstrated clearly that lyssaviruses are subjected to the strongest constraints against amino acid substitutions, probably related to their particular niche and unique pathobiology.

The GenBank/EMBL/DDBJ accession numbers of the rhabdovirus partial N gene sequences obtained in this study are DQ457098–DQ457104, as outlined in Table 1Up.


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
The family Rhabdoviridae includes representatives of six defined genera and more than 100 unassigned rhabdoviruses isolated worldwide (Tordo et al., 2004Down). In this family, members of the genus Lyssavirus are the most important in terms of public health and economic losses. Members of the genera Vesiculovirus and Ephemerovirus are of particular veterinary concern because of livestock diseases, whereas representatives of the genus Novirhabdovirus circulate among fish and other aquatic animals and members of the genera Nucleorhabdovirus and Cytorhabdovirus are plant parasites.

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., 1989Down). When the terms ‘rabies-related viruses' and ‘rabies serogroup’ were introduced (Shope et al., 1970Down), 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., 1971Down; Kemp et al., 1973Down; Tignor et al., 1977Down). 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, 1975Down), 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, 1982Down). Later, other rhabdoviruses, isolated from both arthropods and vertebrates, predominantly in Africa, were described as related serologically to those listed above (Calisher et al., 1989Down). 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., 1990Down). A number of new lyssaviruses were isolated in different parts of the world (Gould et al., 1998Down; Kuzmin et al., 2003Down, 2005Down). Phylogenetically, lyssaviruses segregate into a single monophyletic clade, distant from other rhabdovirus genera and unclassified rhabdoviruses (Kuzmin & Rupprecht, 2005Down; Springfeld et al., 2005Down). A phylogenetic study based on limited polymerase (L) gene sequences was recently described for 38 rhabdoviruses (Bourhy et al., 2005Down). 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
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Viruses (Table 1Down) were obtained from the reference collection of the Centers for Disease Control and Prevention (Ft Collins, CO, USA) and were subjected to one intracerebral passage in suckling mice (Koprowski, 1996Down). Total RNA was extracted from infected mouse brains, using TRIzol (Gibco-BRL) according to the manufacturer's instructions. Two methods were used for the amplification of the rhabdovirus N gene. (i) Conventional RT-PCR was carried out with primers designed to conserved regions of the rhabdovirus N gene, RHNF631 [3'-GAYATGTTYTTYTCAAG-5'; messenger sense; positions 701–717 relative to the Pasteur virus RABV (PV) genome (GenBank accession no. M13215 [GenBank] )], RHNB1520 (5'-AGTCTCYTCTGCCATYTC-3'; genomic sense; positions 1568–1585) and RHNB2134 (5'-AGGATCTRGARGGGAACTTRT-3'; genomic sense; positions 2151–2171). The reactions were performed as described previously (Sacramento et al., 1991Down) with modifications (Kuzmin et al., 2003Down). Once a partial sequence had been obtained and the preliminary phylogenetic placement of the virus was evaluated, additional primers were constructed based on the alignment of the closest relatives of the newly sequenced virus. (ii) BD SMART PCR cDNA synthesis, followed by BD Advantage long-distance PCR (BD Bioscience Clontech), was performed according to the manufacturer's recommendations.


View this table:
[in this window]
[in a new window]
 
Table 1. Rhabdoviruses sequenced in the present study

 
PCR products of the expected size were separated on low-melting-temperature agarose gel, excised and purified using the Wizard PCR Preps DNA purification system (Promega) according to the manufacturer's instructions and then inserted into the pGEM T-Easy vector (Promega) and cloned in E. coli JM109 competent cells (Promega). Between 36 and 48 clones of each library were screened by limited sequencing using the Big Dye Terminator cycle sequencing ready reaction kit (Applied Biosystems) and a primer complementary to the pGEM T-Easy vector. The sequencing products were purified with Centrisep spin columns (Princeton Separations) and processed on an ABI Prism 377 DNA Sequencer (Applied Biosystems) according to the manufacturers' instructions. If cDNA to virus genome was detected, the clones were sequenced completely and a consensus sequence was obtained for each set of such clones. Additionally, primers complementary to the 3' and 5' ends of the consensus sequence were constructed. Conventional RT-PCR with these primers was implemented with the original mouse brain samples. Products of these RT-PCRs were sequenced as well to confirm the results obtained by cloning.

Primary assembly, alignment and consensus sequence generation, as well as DNA translation, were performed in BioEdit software (Hall, 1999Down). 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 2Down). 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.


View this table:
[in this window]
[in a new window]
 
Table 2. Full-length rhabdovirus N gene sequences used in the analysis

 
Alignments were produced using the CLUSTAL X package (Jeanmougin et al., 1998Down). Neighbour-joining (NJ) analysis using either p-distances or Kimura's two-parameter and maximum-parsimony (MP) analysis were performed in MEGA software, version 2.1 (Kumar et al., 2001Down). Maximum-likelihood (ML) analysis using a transition/transversion ratio of 2.0 with empirical base frequencies, gamma distribution of rate variation among sites and the Hidden Markov model was performed by the method implemented in the PHYLIP package, version 3.61 (Felsenstein, 1993Down). At the first step, the dataset was subjected to multiple replication in the SEQBOOT module. The file obtained was subsequently processed in the DNAML module and a consensus tree was generated using the CONSENSE module. Bootstrap values were determined for 1000 replicates in the NJ and MP methods and 100 replicates in the ML method.

For positive selection analysis, 241 codons of the N gene were selected for a number of representative rhabdoviruses, producing alignments without gaps (positions 386–1108 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, 2005aDown). 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, 2005bDown). 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
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Partial sequences of OBOV, KOLV and SJAV were obtained by conventional RT-PCR using primers to conserved regions of the N gene. Further elongation of the sequences was possible once general phylogenetic placement of these viruses was established, and new RT-PCR primers were constructed to the most closely related rhabdovirus representatives. BD SMART PCR cDNA synthesis, followed by BD Advantage long-distance PCR, yielded amplification of partial N gene sequences of MEBV, KCV, RBUV and KOTV. Entire N gene sequences were not generated; however, the 5' end of the gene was present in all four sequences obtained by the latter method. Truncation of the 3' end was probably caused by the quality of the RNA extracted by TRIzol reagent. The length of the sequences obtained (947–1287 nt) was adequate for phylogenetic comparisons.

Phylogenetic analysis recognized all seven newly sequenced viruses as members of the Dimarhabdovirus supergroup (Bourhy et al., 2005Down), without particular relatedness to the genus Lyssavirus (Fig. 1Down). 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 (76–97 % 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 (41–67 %). 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 (32–67 % 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.


Figure 1
View larger version (20K):
[in this window]
[in a new window]
 
Fig. 1. Neighbour-joining phylogenetic trees of the Rhabdoviridae based on N gene sequences (left, alignment of 915 nucleotides, positions 439–1329 relative to the PV genome; right, alignment of 305 deduced amino acids). Branch lengths are drawn to scale and bootstrap values for 1000 replicates are given as percentages for nodes. Gaps at the 3' and 5' ends of the alignment have been truncated to ensure maximum compatibility of the compared sequences; gaps within the alignment were treated by pairwise deletion. I, Vesiculovirus; II, Ephemerovirus; III, Lyssavirus; IV, Novirhabdovirus; V, Cytorhabdovirus; VI, Nucleorhabdovirus. Bars, 0.05 changes per nucleotide position (left) and 0.2 changes per amino acid position (right).

 
Limited amino acid sequence alignment of the deduced nucleoprotein (Fig. 2Down) demonstrated the presence of highly conserved residues and motifs (particularly SPYS at positions 171–174 of the alignment).


Figure 2
View larger version (105K):
[in this window]
[in a new window]
 
Fig. 2. Alignment of 305 amino acids of N gene of dimarhabdoviruses (for positions see Fig. 1Up). Grey background indicates residues conserved in 80 % of sequences; black background indicates residues conserved in 100 % of sequences. Residues that are also conserved in lyssaviruses and SIGMAV are marked by dashes above the alignment.

 
Estimations of dN/dS using a partial N gene fragment (241 codons) for the three rhabdovirus genera revealed contrasting patterns of evolution. The two genera from the Dimarhabdovirus supergroup (ephemeroviruses and vesiculoviruses) and lyssaviruses (Table 3Down; Fig. 3Down) demonstrated global dN/dS values which were quite different, in the order Ephemerovirus > Vesiculovirus > Lyssavirus (with an overlap of 95 % likelihood profile limits for vesiculoviruses and ephemeroviruses). The magnitude of this ratio corresponded to the number of negatively selected codons. The accumulation of dS appears evenly distributed along the gene fragment for all three genera. No evidence of positive selection was found at any codon for any of the datasets tested.


View this table:
[in this window]
[in a new window]
 
Table 3. dN/dS ratios and numbers of positively and negatively selected sites in three datasets as inferred from the SLAC substitution model

Datasets are represented by partial N gene sequences (241 codons, positions 386–1108 relative to the PV RABV genome). The three datasets included the following viruses. Ephemerovirus: BEFV, ARV, FLAV, OBOV, KOTV and RBUV. Vesiculovirus: ISFV, SVCV, PIRYV, CHPV, VSINV (three lineages) and VSNJV (three lineages). Lyssavirus: RABV, LBV, MOKV, DUVV, EBLV-1, EBLV-2, ABLV, ARAV, KHUV, IRKV and WCBV.

 

Figure 3
View larger version (66K):
[in this window]
[in a new window]
 
Fig. 3. Results of SLAC analysis for a 241 codon fragment of the N gene of three rhabdovirus genera. Horizontal scale shows codon position. Bars above the horizontal axis show dN counts, bars below the axis show dS counts.

 
The identities of a longer (305 codons) nucleotide and amino acid alignment (Table 4Down) of the same rhabdovirus genera, together with the tentatively joined ‘Oita group’ viruses and Sigma virus (SIGMAV), demonstrate that amino acid identities are greater than nucleotide identities only for lyssaviruses. The lyssaviruses also demonstrated a higher degree of identity to each other than representatives of other genera, except ephemeroviruses, where high identity was registered between ARV and OBOV. The identities between genera were always significantly lower than within the genera, and nucleotide identities between genera were always greater than amino acid identities between them.


View this table:
[in this window]
[in a new window]
 
Table 4. Identities within and between the genera Lyssavirus, Ephemerovirus and Vesiculovirus, ‘Oita-group’ viruses and SIGMAV

Values are percentage identity inferred from a truncated N gene alignment (915 nucleotides, positions 439–1329 relative to PV genome; 305 deduced amino acids), with pairwise gap exclusion applied. The groups contained the following viruses. Lyssavirus: RABV, LBV, MOKV, DUVV, EBLV-1, EBLV-2, ABLV, ARAV, KHUV, IRKV and WCBV. Vesiculovirus: ISFV, SVCV, PIRYV, CHPV, VSINV and VSNJV. Ephemerovirus: BEFV, ARV, OBOV, KOTV and RBUV. ‘Oita-group’ viruses: OITAV, MEBV and KCV.

 

   DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
We focused our study on the N gene because this gene has been sequenced for the largest number of rhabdoviruses available for comparison in GenBank (Tordo et al., 2004Down). Additionally, this gene provides the best resolution for lyssavirus phylogeny (Kuzmin et al., 2005Down). Rhabdovirus phylogeny based on limited L gene sequence (450 nt) gave inconsistent branching within lyssaviruses and poor discrimination for some other rhabdovirus lineages (Bourhy et al., 2005Down).

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, 1982Down; Calisher et al., 1989Down). Instead, they were placed within the Dimarhabdovirus supergroup defined in a recent study (Bourhy et al., 2005Down). 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, 1999Down; Walker, 2005Down). 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., 2005Down). 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., 2004Down), 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. 2Up) are unclear. They do not overlap with any antigenic epitopes described for RABV (Dietzschold et al., 1987Down; Ertl et al., 1991Down; Fu et al., 1994Down; da Cruz et al., 2001Down). 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, 1978Down; Wang et al., 1995Down).

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., 2005Down). 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., 2005Down) and EBLV-1 (Davis et al., 2005Down), 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, 2005bDown) and previous findings of largely purifying selection in lyssavirus N genes (Holmes et al., 2002Down; Kuzmin et al., 2004Down; Hughes et al., 2005Down; Davis et al., 2005Down), 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., 1999Down; Tordo et al., 2004Down; Fu, 2005Down; Bourhy et al., 2005Down).

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, 1982Down; Calisher et al., 1989Down) 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., 2004Down), 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
 
We thank Dr C. Calisher (retired) for suggestions on the list of viruses which demonstrated serologic cross-reactivity with lyssaviruses. We are grateful to Drs A. Powers and S. Bartlet (DVBID, CDC, Ft Collins, CO, USA) for providing viruses from the CDC Arbovirus Reference Collection. We thank Dr X. Wu (DVRD, CDC, Atlanta, GA, USA) for his expertise and help in the molecular cloning. We also thank Sergei Kosakovsky Pond for assistance with the positive selection analysis. The work was supported in part by the Association of Public Health Laboratories (Washington, DC, USA; International EID fellowship program of 2002–2005). The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the funding agency.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Bourhy, H., Cowley, J. A., Larrous, F., Holmes, E. C. & Walker, P. J. (2005). Phylogenetic relationships among rhabdoviruses inferred using the L polymerase gene. J Gen Virol 86, 2849–2858.[Abstract/Free Full Text]

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, 241–250.[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, 241–257.[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, 45–148.[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, 731–736.[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, 10487–10497.[Abstract/Free Full Text]

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, 103–125.[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, 1–10.[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, 1–24.[Medline]

Fu, Z. F., Wunner, W. H. & Dietzschold, B. (1994). Immunoprotection by rabies virus nucleoprotein. In Lyssaviruses, pp. 207–218. 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, 165–187.[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, 95–98.

Holmes, E. C., Woelk, C. H., Kassis, R. & Bourhy, H. (2002). Genetic constraints and the adaptive evolution of rabies virus in nature. Virology 292, 247–257.[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, 1467–1474.[Abstract/Free Full Text]

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, 1139–1154.[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, 403–405.[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, 95–106.[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, 43–49.[Abstract/Free Full Text]

King, A., Davis, P. & Lawrie, A. (1990). The rabies viruses of bats. Vet Microbiol 23, 165–174.[CrossRef][Medline]

Koprowski, H. (1996). The mouse inoculation test. In Laboratory Techniques in Rabies, 4th edn, pp. 80–86. 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, 2531–2533.[Abstract/Free Full Text]

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, 1208–1222.[Abstract/Free Full Text]

Kumar, S., Tamura, K., Jakobsen, I.-B. & Nei, M. (2001). MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17, 1244–1245.[Abstract/Free Full Text]

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, 65–79.[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, 617–631.[Abstract/Free Full Text]

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, 28–43.[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, 767–769.[Medline]

Nandi, S. & Negi, B. S. (1999). Bovine ephemeral fever: a review. Comp Immunol Microbiol Infect Dis 22, 81–91.[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, 229–240.[CrossRef][Medline]

Shope, R. E. (1982). Rabies-related viruses. Yale J Biol Med 55, 271–275.[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, 690–692.[Abstract/Free Full Text]

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, 6781–6790.[Abstract/Free Full Text]

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, 595–611.[Abstract/Free Full Text]

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. 623–644. 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, 57–80.[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, 995–999.[Abstract/Free Full Text]

Received 27 January 2006; accepted 11 April 2006.


This article has been cited by other articles:


Home page
Proc R Soc BHome page
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]


Home page
J. Clin. Microbiol.Home page
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]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
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 Kuzmin, I. V.
Right arrow Articles by Rupprecht, C. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kuzmin, I. V.
Right arrow Articles by Rupprecht, C. E.
Agricola
Right arrow Articles by Kuzmin, I. V.
Right arrow Articles by Rupprecht, C. E.


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