JWD Your personal alerts
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Journal of Wildlife Diseases, 40(4), 2004, pp. 617-631
© Wildlife Disease Association  2004
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
Services
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 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.

MOLECULAR EPIDEMIOLOGY OF TERRESTRIAL RABIES IN THE FORMER SOVIET UNION

Ivan V. Kuzmin1,2, Alexandr D. Botvinkin3, Lorraine M. McElhinney4, Jean S. Smith2, Lillian A. Orciari2, Gareth J. Hughes2, Anthony R. Fooks4 and Charles E. Rupprecht2,5

1 Rabies Group, Institute for Natural Foci Infections, Prospekt Mira, 7, Omsk, 644080 Russia
2 Rabies Section, Centers for Disease Control and Prevention, 1600 Clifton Rd., Atlanta, Georgia 30333, USA
3 Epidemiology Chair, State Medical University, Krasnogo Vosstania, 1, Irkutsk 664003, Russia
4 Rabies Research and Diagnostic Group, Veterinary Laboratories Agency, Weybridge, Addlestone, Surrey KT15 3NB, UK

5 Corresponding author (email: cyr5{at}cdc.gov)


   ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Fifty-five rabies virus isolates originating from different regions of the former Soviet Union (FSU) were compared with isolates originating from Eurasia, Africa, and North America according to complete or partial nucleoprotein (N) gene sequences. The FSU isolates formed five distinct groups. Group A represented viruses originating from the Arctic, which were similar to viruses from Alaska and Canada. Group B consisted of "Arctic-like" viruses, originating from the south of East Siberia and the Far East. Group C consisted of viruses circulating in the steppe and forest-steppe territories from the European part of Russia to Tuva and in Kazakhstan. These three phylogenetic groups were clearly different from the European cluster. Viruses of group D circulate near the western border of Russia. Their phylogenetic position is intermediate between group C and the European cluster. Group E consisted of viruses originating from the northwestern part of Russia and comprised a "northeastern Europe" group described earlier from the Baltic region. According to surveillance data, a specific host can be defined clearly only for group A (arctic fox; Alopex lagopus) and for the Far Eastern part of the group B distribution area (raccoon dog; Nyctereutes procyonoides). For other territories and rabies virus variants, the red fox (Vulpes vulpes) is the main virus reservoir. However, the steppe fox (Vulpes corsac), wolf (Canis lupus), and raccoon dog are also involved in virus circulation, depending on host population density. These molecular data, joined with surveillance information, demonstrate that the current fox rabies epizootic in the territory of the FSU developed independently of central and western Europe. No evidence of positive selection was found in the N genes of the isolates. In the glycoprotein gene, evidence of positive selection was strongly suggested in codons 156, 160, and 183. At these sites, no link between amino acid substitutions and phylogenetic placement or specific host species was detected.

  Key words:  Epidemiology, former Soviet Union, rabies, phylogenetics, Russia.


   INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Rabies is enzootic in Eurasia. The Asian tropics are affected mainly by urban dog rabies, and thousands of human cases are reported there annually (Ahuja et al., 1985; Wilde, 1997). In Europe, urban dog rabies has been largely eliminated. The main wildlife reservoir is the red fox (Vulpes vulpes), with only sporadic human cases registered in European countries over the last few decades, mainly because of successfully implemented oral vaccination programs and advanced human postexposure prophylaxis (Toma and Andral, 1977; Bourhy et al., 1999). In the former Soviet Union (FSU), rabies is enzootic from the western borders to the Far East, but the affected areas are discontinuous (Fig. 1Go). General epidemiologic patterns of the disease, surveillance data, and veterinary records suggest that rabies is maintained in Russia by wild canids. During the last several decades, only a few local dog epizootics have occurred sporadically in some territories. Wildlife rabies predominates also in the Ukraine, Belarus, Moldova, and the Baltic states. Dog rabies is predominant in the central Asiatic states of the FSU, especially in Uzbekistan and Tajikistan and in the Caucasus (Cherkasskiy, 1985; Selimov, 1998; Vedernikov et al., 2002).



View larger version (28K):
[in this window]
[in a new window]
 
FIGURE 1. Map of the former Soviet Union showing the terrestrial rabies enzootic region (shaded) and distribution of different phylogenetic groups of rabies virus (A–E). State designation: 1 = Estonia; 2 = Lithuania; 3 = Latvia; 4 = Moldova; 5 = Georgia; 6 = Armenia; 7 = Azerbaijan; 8 = Tajikistan; 9 = Kyrgyzstan. RV305, RV307, and RV308 = un-grouped isolates from Georgia.

 
Certain wild Carnivora species appear to maintain rabies virus circulation in different regions. Arctic rabies covers the tundra and tundra-forest zones in the north, and this enzootic region matches the area of distribution and regular migrations of the arctic fox (Alopex lagopus). A wide zone of conifer taiga forests lying to the south of the tundra-forest zone is considered largely free of rabies. In general, this region is not supportive of wild canids, and their density appears too low to maintain active virus circulation. In the forest-steppe and steppe zones, lying to the south of the taiga, rabies is maintained primarily by the red fox. The raccoon dog (Nyctereutes procyonoides) is involved in virus circulation in the Far East and parts of Europe. The steppe fox (Vulpes corsac) and jackal (Canis aureus) participate in rabies virus circulation in the steppe and desert territories. The wolf (Canis lupus) historically caused rabies outbreaks in different areas. These six canid species were included in the group of likely principal hosts and vectors of rabies virus in Russia and the FSU. Mustelids and other carnivores were considered as additional or occasional hosts (Kantorovich and Reshetnikov, 1968; Selimov, 1978; Malkov and Gribanova, 1980). Rodents and other terrestrial mammals have been diagnosed rabid, but only very rarely (Selimov, 1998; Vedernikov et al., 2002).

At least nine different antigenic variants of rabies virus were detected in the FSU using antinucleocapsid monoclonal antibodies (N-MAbs). A limited relationship of the antigenic patterns to geographic origin and terrestrial host species of the isolate was demonstrated (Botvinkin et al., 1990; Selimov et al., 1994).

Few isolates originating from Russia and states of the FSU were available for genetic comparison. Comparative analyses of some of these isolates from the FSU have been published, but not with specific attention to the FSU (Kissi et al., 1995; McElhinney et al., 2001; Nadin-Davis et al., 2002; Johnson et al., 2003). In this study, we performed phylogenetic analysis of 55 terrestrial rabies virus isolates, originating from different geographic locations of Russia and some surrounding states of the FSU, and compared them with isolates originating from Eurasia, Africa, and North America. Our overall objective was to draw inferences into the origins, reservoirs, and regional epidemiology of the disease, as well as to search for comparative evidence of selection among viral genes.


   MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
All viruses were isolated from rabid animals by intracerebral mouse inoculation (Koprowski, 1996) and were used after a minimum of one to seven mouse brain passages (Table 1Go). The samples were initially identified by the direct fluorescent antibody test (Dean et al., 1996) or the rapid rabies enzyme immunosorbent assay (Bourhy and Perrin, 1996).


View this table:
[in this window]
[in a new window]
 
TABLE 1. Description of viruses used for study of the molecular epidemiology of rabies in the former Soviet Union.
 
Total RNA was extracted from infected mouse brains with TRIzolTM (Gibco-BRL, Inc., Gaithersburg, Maryland, USA) according to the manufacturer’s recommendations. Reverse transcriptase polymerase chain reaction (RT-PCR) was performed with primer sets to the nucleoprotein (N) and glycoprotein (G) genes and sequenced as described elsewhere (Tordo et al., 1996; Johnson et al., 2003).

Two data sets were used for the phylogenetic analysis. The first set (n=74, including 39 isolates from the FSU) consisted of 1,350 nucleotide (nt) sequences of the entire N gene (positions 71–1,420 according to the Pasteur virus [PV] genome; GenBank accession number M13215). The second set (n=92, including 55 isolates from the FSU) consisted of 405 nt sequences of the N gene (positions 71–475). Phylogenetic analysis was performed by the neighbor-joining (NJ) method from the MEGA computer program, version 2.1 (Kumar et al., 2001). Maximum likelihood (ML) analysis, with empirical base frequencies, gamma distribution of rate variation among sites, and the hidden Markov model of inferring different rates of evolution at different sites, was performed by the method implemented in PHYLIP, version 3.6-alpha (Felsenstein, 1993). Initially, the data set was subjected to multiple replications in the SEQBOOT module. The obtained file was processed subsequently in the DNAPARS, PROT-PARS, or DNAML modules, respectively, and finally, the consensus tree was generated with the CONSENSE module. Bootstrap values were determined for 1,000 replicates by the NJ method and for 100 replicates by the ML method.

Three data sets were used for positive selection analysis: 1) entire N gene sequences (n=67); 2) 852 nt G gene sequences (n=36; 284 codons, positions 3,408–4,259); and 3) 251 nt G gene sequences (n=52; 117 codons, positions 4,009–4,259). Positive selection analysis was performed with various models of codon substitution (Yang et al., 2000) implemented in the CODEML program of the PAML package (Yang, 1997). Various ML models are applied with varying constraints on the values of synonymous (dS) and nonsynonymous (dN) substitution rates, and their ratio ({omega}). Models allowing for positive selection (i.e., {omega}>1) are nested within models that do not allow for positive selection. This allows the significance of the fit for positive selection models to be tested by the likelihood ratio test. Positive selection is inferred if the positive selection model has a significantly higher likelihood than the null model, and a value of {omega}>1 is estimated. If evidence of positive selection is suggested, Bayesian methods are used to identify which individual codons fall into the {omega}>1 class.

Sequence alignments were trimmed to include only complete nonstop codons, and identical sequences were removed. For each set of sequences, a ML tree was generated with PAUP* 4.0, beta 10 (Swofford, 2000) by a heuristic search incorporating tree bisection-recombination (TBR) branch swapping. This tree was then used for positive selection analysis. In each case, the model of nucleotide substitution was selected with Modeltest (Posada and Crandall, 1998). Model testing with CODEML was performed as described by Woelk et al. (2002). Phylogenetic trees used in this analysis are available from the authors on request.


   RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Phylogenetic analysis of the entire N gene sequences, performed by either ML or NJ methods, revealed trees of an identical topology (Fig. 2Go). The sequences were subdivided into two main groups. The first group consisted of two tropical dog viruses from Thailand (8738THA) and India (AF374721) and might be used as an out-group to the second group, which consisted of all other sequences. The second group was further subdivided into two distinct clusters. One of them (cluster I) joined Arctic viruses originating from Eurasia and North America and two branches of "Arctic-like" viruses, originating from India, Pakistan, East Siberia, and the Far East (bootstrap support of 100%). The second (cluster II) contained viruses originating from different territories of Eurasia and Africa (bootstrap support of 99%). Included in this cluster were representatives of all European groups previously identified (Bourhy et al., 1999) and other viruses of Old World origin (Kissi et al., 1995; McElhinney et al., 2001). The topology of the tree for all these sequences matched that of earlier reports.



View larger version (21K):
[in this window]
[in a new window]
 
FIGURE 2. Neighbor joining phylogenetic tree of 74 rabies virus isolates according to the entire nucleoprotein gene. Genbank accession numbers are given when strain names are unavailable. The tree is rooted with the isolates AF374721 and 8738THA. Bootstrap values are presented for key nodes, and branch lengths are drawn to scale. Isolates from the former Soviet Union are italicized. Abbreviations for the main phylogenetic groups from the former Soviet Union (A–E) are described in the text, as well as big clusters (I and II); abbreviations for earlier described European groups (NEE, CE, EE, WE) are given according to Bourhy et al. (1999).

 
Isolates originating from the FSU could be clearly subdivided into five main groups on the basis of phylogenetic placement. Diversity varied primarily with geographic location rather than host species, even when a number of isolates from different species were available for comparison. Group A consisted of isolates originating from the Arctic zone of Eurasia and North America. This cluster was further split into two clades, and the Canadian isolate "1991 fox 6199" was placed ancestral to both. No difference was registered among these Arctic viruses with regard to continent of origin or host species, with one exception: the isolate RV250, obtained from a rodent, Citellus undulatus, in the south of East Siberia, distant from the Arctic zone.

Group B consisted of Arctic-like viruses originating from East Siberia and the Far East, including the Japanese isolate, Komatsugawa. Another Arctic-like group included viruses originating from India and Pakistan (isolates RV61, 196p, and 277p). Group C consisted of viruses originating from the steppe and forest-steppe territories from the European part of Russia to the Altai and Sayan mountains. Differences were well supported from the other group of viruses, which joined isolates from different parts of Europe (Kissi et al., 1995; Bourhy et al., 1999), with bootstrap support of 100%.

Group D joined viruses originating from the center of the European part of Russia (RV234, RV299, and RV262) and an isolate from Hungary, 9215HON, which was not related to any clade from an earlier study (Bourhy et al., 1999). This group was linked to a cluster of European viruses with limited bootstrap support (52%). Group E consisted of viruses isolated in the northwestern part of Russia. Together with isolates from Estonia, these viruses belonged to the earlier described "northeastern Europe" group (Bourhy et al., 1999). The isolate RV308 from Georgia was not classified into any cluster.

When the extended set of 405 nt N gene sequences was analyzed, tree topology was identical. All virus groups determined for the entire N gene remained and were extended with the new representatives, whereas bootstrap support was reduced (Fig. 3Go). In this tree, the newly incorporated isolates RV305 from Georgia and RV1133 from Turkey were joined to isolate RV308 and formed a solid group with substantial bootstrap support. Another newly incorporated isolate from Georgia, RV307, was joined to the cluster of Middle East viruses.



View larger version (21K):
[in this window]
[in a new window]
 
FIGURE 3. Neighbor joining phylogenetic tree of 92 rabies virus isolates according to 405 nucleotide fragment of the nucleoprotein gene. For explanations see legend to Figure 2Go.

 
No evidence of positive selection was found in any set of N gene sequences. In a number of cases, the null model was significantly less likely than one that allowed positive selection, but in no case was {omega}>1. The results of positive selection analysis for G gene sequences are shown in Table 2Go. For both G gene data sets, a number of null models were less favored than models that allowed for positive selection. In both cases, this includes the most precise model (M8) for detecting positive selection (Anisimova et al., 2002). For the sequence set covering the larger area of the G gene (data set 2), two positive selection models (M3 and M8) selected three codons (156, 160, 183) as having a {omega}=2.2709 and 1.7427, respectively. For the smaller sequence set (data set 3), positive selection was inferred at just one codon (183) by M3 and M8, with {omega}=2.3301 and 2.3266, respectively.


View this table:
[in this window]
[in a new window]
 
TABLE 2. Results of positive selection analysis
 

   DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Geographic distribution of the genetic variants of rabies virus is presented in Figure 1Go. The isolates belonging to groups A–E originated from distinct geographic regions. Arctic rabies is distributed in a circumpolar fashion. Absence of significant distinctions between North American and Siberian Arctic rabies virus isolates suggests that the virus population is intermixed. This is in concordance with the results of N-MAbs application (Botvinkin et al., 1990) and previous genetic data (Kissi et al., 1995). Our set of sequences was limited, and we could not make a species comparison. However, surveillance data suggests that the arctic fox is the principal host of rabies virus in the Arctic region of Russia. Rare human rabies cases caused by arctic foxes, dogs, and wolves were described in northern Russian (Kantorovich and Reshetnikov, 1968; Botvinkin et al., 1995; Kuzmin, 1999). Additional attention should be paid to isolate RV250. The isolation point of this virus is separated from the native area of Arctic rabies virus circulation by at least 1,000 km of the "rabies-free" zone of the taiga forests. Further active surveillance of C. undulatus and another rodent species did not bring positive results. According to antigenic typing, this isolate reacted with N-MAb P-41, like the Arctic rabies viruses, but did not react with some other N-MAbs, being different from all other rabies viruses (n>300) available for comparison (Botvinkin et al., 1990). Isolates with similar antigenic patterns were described from Thailand and Madagascar (Sureau et al., 1983). We have no satisfactory explanation for the origin of this virus.

Circulation of group B viruses in the Russian Far East is currently maintained by raccoon dogs, red foxes, and wolves. According to available surveillance data, rabid raccoon dogs have been identified in the Far East since 1931 (Mirolubov, 1934), and these animals have been known as sources of human rabies since at least 1951. During the outbreak of 1979–80, active surveillance demonstrated rabies virus in 13.5±4.6% of raccoon dogs and in 3.5±2.0% of red foxes trapped randomly (Botvinkin et al., 1981). In an experimental study, susceptibility of raccoon dogs to indigenous rabies virus was significantly greater than susceptibility of foxes. During the clinical period, the animals were agitated and aggressive, and the virus was detected in their salivary glands more frequently than in the salivary glands of foxes (Botvinkin et al., 1983). These data provide an indirect suggestion that the raccoon dog might be the principal host of this rabies virus variant. However, about 80% of human rabies cases were caused by dogs, with raccoon dogs the next most important source. Only one human case was reported after a red fox bite (Savitsky et al., 1981; Yanovich, 2003).

In the Russian Far East, the rabies enzootic region is limited and fragmented. It lies in the wet lowlands along the Amur and Ussuri rivers. Most probably, a major part of the enzootic region is situated in the northern provinces of China. Unfortunately, we had no samples from this area to make a comparison. The Komatsugawa virus, isolated from a dog in Tokyo shortly after the end of World War II (Ito et al., 1999), belonged to the same group.

Another lineage of group B (isolates 248c and 304c) originated from a local rabies focus situated in the highland steppes of East Siberia near the border with Mongolia (Transbaikal region). During the first half of the 20th century, dogs and wolves were known as the main reservoirs and vectors of rabies virus in this area. Later, from 1950 to 1977, wolves were a source of 8.7% of human rabies cases and domestic animals (mainly dogs) were a source of the remaining 91.3%. Rabies in red fox was first diagnosed in 1961. From 1971 to 1977, rabid red foxes and steppe foxes were identified by veterinary laboratories almost annually, but foxes never caused a documented human case (Botvinkin et al., 1980). Both isolates, reported in this study, were found in 1977 as a result of active surveillance. Their Arctic-like properties were demonstrated in reactions with N-Mabs; both viruses reacted with N-MAb P-41 (Selimov et al., 1994). Rabies was not reported from this territory after 1984. It was not possible to discern whether this focus has been maintained. However, rabies is still occasionally reported from Mongolia (Angar, 2001).

Both group B lineages, as well as other Arctic-like viruses originating from India and Pakistan, are separated from the northern area of Arctic rabies virus circulation by thousands of kilometers of allegedly rabies-free territory. It appears that these viruses originate from one progenitor, and current diversity could be the result of further independent evolution. However, it was not possible to make an appropriate age estimation because of lack of samples belonging to one lineage separated by a number of years.

Group C viruses formed a solid cluster with considerable bootstrap support. They were isolated in the huge steppe and forest-steppe territory of southeastern Europe, West Siberia, Kazakhstan, and Tuva. The last focus in Tuva is separated from the continuous western area of virus circulation by the Altay mountains. It appears to be a marginal zone of the rabies enzootic territory of the Mongolian highland steppes.

According to veterinary records and field surveillance data, the disease in the steppe zone of the FSU is maintained mainly by the red fox. The steppe fox participates in virus circulation concurrently with the red fox, particularly in territories of high population density (Sansizbaev, 1975; Malkov and Gribanova, 1980; Kuzmin et al., 2001). We did not register any genetic distinctions between the red fox and steppe fox isolates.

In the area of group C circulation, the first wildlife rabies outbreak was reported among raccoon dogs, foxes, and wolves in 1945 in the region of the Volga River (Isakov, 1949). Since 1951, the red fox has been considered the main rabies reservoir (Amitrov, 1956; Nazarov, 1961). In 1949, a fox rabies outbreak was described in Kazakhstan (Sludsky, 1954). In West Siberia, fox rabies has been registered in Novosibirsk province since 1958 and in Altay province since 1961 (Vyazhevich, 1959; Kiryanov, 1962). Dogs (with sporadic cases caused by wolves) were the major source of human rabies in the first half of the 20th century, but foxes predominated as a source of human rabies in the second half of the century (Selimov, 1978). Group D viruses originated from the western European part of Russia, including the junction of the Russian, Ukrainian, and Belorussian borders (so-called Polesie). This group retained an intermediate phylogenetic position between the lineages originating from other parts of Europe and group C, which was concordant with geographic distribution of the viruses. Particular antigenic patterns of these viruses were also suggested with N-MAbs (Selimov et al., 1988). Red fox rabies, with involvement of other wildlife (the raccoon dog, badger (Meles meles), and sometimes small mustelids), predominated in this territory after World War II, and a high incidence of human rabies has been found there during the last few decades (Selimov, 1978; Cherkasskiy et al., 1995). We had only two human isolates originating from this region, one of them belonged to group D (RV239), but the second belonged to group C (RV241). Because of limited samples, we could not establish which virus lineage caused human disease more frequently.

Group E in our study was the same as the "northeastern Europe" group described earlier in Bourhy et al. (1999). The isolates from the Baltic region (northwestern part of Russia, Estonia, and Finland) reacted with N-MAb P-41 (Kulonen and Boldina, 1993; Selimov et al., 1994). Wildlife rabies predominated there after World War II (Selimov, 1978, 1998; Cherkasskiy et al., 1995). Because raccoon dogs were frequently involved in virus circulation in this region and because the viruses formed a separate lineage, the raccoon dog was proposed to be an intermediary host in the switch of rabies virus from dogs and wolves to foxes (Bourhy et al., 1999). However, there is no support for this statement in the epidemiologic surveillance data and historical veterinary records (Selimov, 1978, 1998; Cherkasskiy et al., 1995). The raccoon dog population density in the region is estimated at 3–20 times less than that of the red fox population (Nasymovich and Isakov, 1985). It appears that the red fox is more affected in the Russian part of this area. However, raccoon dog cases predominated during the rabies outbreak in Finland in 1988–89 (Nyberg et al., 1992). Phylogenetically, this group was related to other European lineages, and there is no reason to consider it as specific to the raccoon dog.

The isolates originating from Georgia and Turkey tended to form another lineage that probably represented viruses circulating south of the Black Sea region (Johnson et al., 2003). One of the Georgian isolates (RV307) was joined to viruses originating from the Middle East. According to epidemiologic surveillance, dog rabies predominated in the Transcaucasian states of the FSU with only occasional incidence among wild canids (Selimov, 1978, 1998; Cherkasskiy et al., 1995). We have very few isolates from this territory to make appropriate comparisons. Different lineages of rabies virus appear to circulate in southeastern Europe and the Middle East (Bourhy et al., 1999; Johnson et al., 2003), although further investigations are needed. Because the viruses of group C were not identified south of the Caucasus Mountains, this mountain range could be a physical barrier separating these rabies virus populations.

All defined virus groups corresponded to certain geographic regions. Conversely, an association with host species was less obvious. The principal host might be identified only in Arctic regions (arctic fox) and in the Far East (raccoon dog). In other territories, rabies virus is apparently associated with foxes (red fox and steppe fox).

As has been suggested by others, the current European epizootic of fox rabies started in eastern Prussia or at the former Russia-Polish border during 1935–45 (Toma and Andral, 1977; Bourhy et al., 1999). Our data clearly demonstrate that viruses of groups B and C could not derive from that area. The distance from the western borders of Russia to the Volga delta is about 1,500 km, and to Siberia and Kazakhstan about 3,000–3,500 km. The carrying capacity of biotopes and properties of fox populations in the western part of European Russia are similar to those described for central and western Europe. The speed of the epizootic front movement was estimated in Europe as 20–60 km per year (Toma and Andral, 1977; Pastoret and Brochier, 1999). With this velocity, the epizootic could not have reached the Volga River region by the time of the first outbreak there in 1945. In contrast, the estimated carrying capacity of Asian steppe biotopes appears limited. Fox populations here are mobile, and the animals seem to maintain rabies epizootics at a population density of 0.3–0.8 foxes per km2 (Malkov and Gribanova, 1980). The speed of the epizootic front movement in Western Siberia was estimated as 160–513 km/yr (Rybak et al., 1992). There are no natural barriers to rabies progression, even during winter. Hence, the viruses of group C could spread along the steppe zone in Asia rapidly. It is not clear where the epizootic started. Additionally, small canids maintain viruses of group B, as well as other lineages of rabies virus in the Middle East (David et al., 2000; Johnson et al., 2003). Perhaps different populations of rabies virus independently switched to small canids in distinct areas of Eurasia during 1930–70.

Previous analyses of rabies virus sequences for evidence of positive selection have suggested that rabies virus evolves largely under purifying or neutral selection (Badrane and Tordo, 2001; Badrane et al., 2001; Holmes et al., 2002). Application of methods identical to those applied here for rabies virus isolates from different geographic regions suggested that positive selection was only occurring at a single codon in the ectodomain of the G gene (position 183). Both G gene data sets used in this study also selected position 183 as being under positive selection. Additionally, in this study, we found a further two codons to be under positive selection, corresponding to ~1% of the codons analyzed. Although this value is low, it remains in the range estimated for genes of other negative-strand RNA viruses (Woelk et al., 2002).

In nature, rabies virus exists in distinct lineages with little interaction of variants following division (Bourhy et al., 1999; Smith, 2002). We suggest that the detection of a greater number of codons under positive selection here might be a result of the restriction of this analysis to a defined subset of virus lineages, rather than an amalgamation of every available complete gene sequence (Holmes et al., 2002). In reality, processes of positive selection are likely to be lineage-dependent, and an assessment of positive selection across the entire rabies virus genotype could distort and dilute individual processes, such that they are no longer detectable. Ideally, the analysis here of Russian rabies virus isolates would have been further divided into individual variants. However, the number of sequences available limits adequate subdivision at this time.

The rabies virus glycoprotein is responsible for binding to host cell receptors and is a primary target of the immune response (Dietzschold et al., 1983). In contrast, the nucleoprotein is largely hidden from immune surveillance. The lack of detectable positive selection in the N gene concurs with the results of others (Holmes et al., 2002). All three of the G gene positive selection sites occur within the neurotoxin-like region (Lentz et al., 1984), a region thought to participate in host cell binding at neuromuscular junctions. These sites might be subject to adaptive evolution through host cell receptor binding, although, using the sequences from this study, we saw no obvious link between amino acid substitutions at these sites and phylogenetic placement. Although the majority of rabies virus replication occurs in neurons that are under relatively weak immune surveillance (Miller, 1999), it might be over simplistic to state that rabies virus evolves according to a neutral model (Badrane and Tordo, 2001). Subtle positive selection at a handful of codons might be sufficient to ensure that limited immune pressure continues to force the rabies virus to adapt to changing host environments, especially in a diversity of small-bodied carnivores.


   ACKNOWLEDGMENTS
 
We are deeply indebted to the late Professor M. Selimov (Institute of Polio and Viral Encephalitis RAMS, Moscow, Russia) and the late A. A. King (Central Veterinary Laboratory, Weybridge, UK), who collected and stored the main part of the samples and made valuable input into their study with monoclonal antibodies. We thank Y. T. Arai (National Institute of Infectious Diseases, Tokyo, Japan) for providing the cDNA of isolate Komatsugawa. We appreciate the help of G. N. Sidorov (State Pedagogical University, Omsk, Russia), E. M. Poleschuk (Institute for Natural Foci Infections, Omsk, Russia), and many staff in the veterinary and sanitary surveillance system who helped us to collect the samples. This study was supported by the Association of Public Health Laboratories (international EID fellowship program of 2002–03) and by the Royal Society. The work of the Veterinary Laboratories Agency is funded by the UK Department for Environment, Food, and Rural Affairs (Defra ROAME SE0417). G.J.H. was funded by an American Society for Microbiology/National Centers for Infectious Disease postdoctoral fellowship.


   LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
AHUJA, S., K. K. TRIPATHI, S. M. SAHA, AND S. N. SAXENA. 1985. Epidemiology of rabies in India. In Rabies in the Tropics, E. Kuwert, C. Merieux, H. Koprowski, and K. Bogel (eds.). Springer-Ver-lag, Berlin, pp. 571–582.

AMITROV, V. E. 1956. Foxes as rabies carriers. Veterinary 6: 33–34. [In Russian.]

ANGAR, D. 2001. Morbidity trend of plague, anthrax and rabies [in Mongolia]. Scientific Journal of the Center of Infectious Diseases with Natural Foci 9: 20–24. [In Mongolian.]

ANISIMOVA, M., J. P. BIELAWSKI, AND Z. YANG. 2002. Accuracy and power of Bayes prediction of amino acid sites under positive selection. Molecular Biology and Evolution 19: 950–958.[Abstract/Free Full Text]

BADRANE, H., AND N. TORDO. 2001. Host switching in lyssavirus history from the Chiroptera to the Carnivora orders. Journal of Virology 75: 8096–8104.[Abstract/Free Full Text]

———, C. BAHLOUL, P. PERRIN, AND N. TORDO. 2001. Evidence of two lyssavirus phylogroups with distinct pathogenicity and immunogenicity. Journal of Virology 75: 3268–3276.[Abstract/Free Full Text]

BOTVINKIN, A. D., V. P. SAVITSKY, G. N. SIDOROV, G. D. MALKOV, L. YA. GRIBANOVA, AND V. A. BUTINA. 1980. About wildlife rabies in Eastern Transbaikal region. In Modern methods in the study of natural foci infections. Institute for Natural Foci Infections, Omsk, Russia, pp. 52–61. [In Russian.]

———, ———, ———, AND V. G. YUDIN. 1981. The role of raccoon dog in the epidemiology of rabies in the Far East. Journal of Microbiology, Epidemiology and Immunobiology 1: 50–54. [In Russian.]

———, L. YA. GRIBANOVA, AND T. A. NIKIPHOROVA. 1983. Experimental rabies in raccoon dog. Journal of Microbiology, Epidemiology and Immunobiology 12: 37–40. [In Russian.]

———, M. A. SELIMOV, E. V. KLUEVA, L. YA. GRIBANOVA, AND N. A. KHISMATULLINA. 1990. Antigenic characteristic of field rabies virus isolates with the help of antinucleocapsid monoclonal antibodies from different regions of the USSR. Journal of Microbiology, Epidemiology and Immunobiology 1: 50–54. [In Russian.]

———, V. F. CHERNIAVSKI, AND I. YA. EGOROV. 1995. Rabies in Siberian subarctic. East Siberian Journal of Infectious Pathology 1: 5–10. [In Russian.]

BOURHY, H., AND P. PERRIN. 1996. Rapid rabies enzyme immunodiagnosis (RREID) for rabies antigen detection. In Laboratory techniques in rabies, 4th Edition, F.-X. Meslin, M. M. Kaplan, and H. Koprowski (eds.). WHO, Geneva, Switzerland, pp. 105–113.

———, B. KISSI, L. AUDRY, M. SMRECZAK, M. SAD-KOWSKA-TODYS, K. KULONEN, N. TORDO, J. F. ZMUDZINSKI, AND E. C. HOLMES. 1999. Ecology and evolution of rabies virus in Europe. Journal of General Virology 80: 2545–2557.[Abstract/Free Full Text]

CHERKASSKIY, B. L. 1985. Epidemiology and prevention of rabies. Medicine, Moscow, USSR. [In Russian.]

———, A. G. KNOP, V. A. VEDERNIKOV, V. A. SEDOV, A. E. KHAIRUSHEV, AND S. A. CHERNICHENKO. 1995. Epidemiology and epizootology of rabies in the territory of the former USSR. Journal of Microbiology, Epidemiology and Immunobiology 1: 21–26. [In Russian.]

CONZELMANN, K. K., J. H. COX, L. G. SCHNEIDER, AND H. J. THIEL. 1990. Molecular cloning and complete nucleotide sequence of the attenuated rabies virus SAD B19. Virology 175: 485–499.[Medline]

DAVID, D., B. YAKOBSON, J. S. SMITH, AND Y. STRAM. 2000. Molecular epidemiology of rabies virus isolates from Israel and other Middle- and Near-Eastern countries. Journal of Clinical Microbiology 38: 755–762.[Abstract/Free Full Text]

DEAN, D. J., M. K. ABELSETH, AND P. ATANASIU. 1996. The fluorescent antibody test. In Laboratory techniques in rabies, 4th Edition, F.-X. Meslin, M. M. Kaplan, and H. Koprowski (eds.). World Health Organizatin, Geneva, Switzerland, pp. 88–93.

DIETZSCHOLD, B., W. H. WUNNER, T. J. WIKTOR, A. D. LOPES, M. LAFON, C. L. SMITH, AND H. KOPROWSKI. 1983. Characterization of an antigenic determinant of the glycoprotein that correlates with pathogenicity of rabies virus. Proceedings of the National Academy of Sciences USA 80: 70–74.[Abstract/Free Full Text]

FELSENSTEIN, J. 1993. PHYLIP: Phylogenetic inference package (version 3.5c). Department of Genetics, University of Washington, Seattle, Washington.

HOLMES, E. C., C. H. WOELK, R. KASSIS, AND H. BOURHY. 2002. Genetic constraints and the adaptive evolution of rabies virus in Nature. Virology 292: 247–257.[Medline]

ISAKOV, YU. A. 1949. Rabies of wild animals inhabiting the region of Volga delta. Issues on the Regional Experimental and Descriptive Parasitology 6: 82–86. [In Russian.]

ITO, N., M. SUGIYAMA, K. ORAVEERAKUL, P. PIYAVIRIYAKUL, B. LUMLERTDACHA, Y. T. ARAI, Y. TAMURA, Y. MORI, AND N. MINAMOTO. 1999. Molecular epidemiology of rabies in Thailand. Microbiology and Immunology 43: 551–559.[Medline]

JOHNSON, N., C. BLACK, J. SMITH, H. UN, L. M. MCELHINNEY, O. AYLAN, AND A. R. FOOKS. 2003. Rabies emergence among foxes in Turkey. Journal of Wildlife Diseases 39: 262–270.[Abstract]

KANTOROVICH, R. A., AND P. M. RESHETNIKOV. 1968. Questions of hydrophobia geography in the USSR. In Proceedings of Scientific Meeting on the Problems of Medical Geography 18–23 September 1968, Moscow, Russia, pp. 188–189. [In Russian.]

KIRYANOV, G. I. 1962. Animal rabies in the Altay region. Proceedings of Irkutsk Plague Control Institute, Irkutsk, Russia, 3: 96–97. [In Russian.]

KISSI, B., N. TORDO, AND H. BOURHY. 1995. Genetic polymorphism in the rabies virus nucleoprotein gene. Virology 209: 526–537.[Medline]

KOPROWSKI, H. 1996. The mouse inoculation test. In Laboratory techniques in rabies, 4th Edition, F.-X. Meslin, M. M. Kaplan, and H. Koprowski (eds.). World Health Organization, Geneva, Switzerland, pp. 80–86.

KULONEN, K., AND I. BOLDINA. 1993. Differentiation of two rabies strains in Estonia with reference to recent Finnish isolates. Journal of Wildlife Diseases 29: 209–213.[Abstract]

KUMAR, S., K. TAMURA, I. B. JAKOBSEN, AND M. NEI. 2001. MEGA2: Molecular evolutionary genetics analysis software. Bioinformatics 17: 1244–1245.[Abstract/Free Full Text]

KUZMIN, I. V. 1999. An arctic fox rabies virus strain as the cause of human rabies in Russian Siberia. Archives of Virology 144: 627–629.[Medline]

———, A. D. BOTVINKIN, G. N. SIDOROV, AND E. I. RECHOV. 2001. Epizootic situation and prospectives of rabies elimination in the south of West Siberia. Journal of Microbiology, Epidemiology and Immunobiology 3: 28–35. [In Russian.]

LENTZ, T. L., P. T. WILSON, E. HAWROT, AND D. W. SPEICHER. 1984. Amino acid sequence similarity between rabies virus glycoprotein and snake venom curaremimetic neurotoxins. Science 226: 847–848.[Abstract/Free Full Text]

MALKOV, G. B., AND L. YA. GRIBANOVA. 1980. Main results of the wildlife rabies study in Siberia and Far East. In Modern methods in study of natural foci infections. Institute for Natural Foci Infections, Omsk, Russia, pp. 52–61. [In Russian.]

MCELHINNEY, L. M., N. JOHNSON, K. L. MANSFIELD, C. J. FINNEGAN, J. SMITH, J.-P. LOWINGS, AND A. R. FOOKS. 2001. Molecular characterization of rabies and rabies-related viruses. In Rabies control in Asia, B. Dodet and F.-X. Meslin (eds.). John Libbey Eurotext, Paris, France, pp. 178–184.

MILLER, D. W. 1999. Immunobiology of the blood-brain barrier. Journal of Neurovirology 5: 570–578.[Medline]

MIROLUBOV, I. N. 1934. Japanese deer diseases. In Questions of deer-breeding. Nauka, Moscow, USSR, pp. 1–22. [In Russian.]

NADIN-DAVIS, S. A., G. A. CASEY, AND A. I. WAN-DELER. 1994. A molecular epidemiological study of rabies virus in central Ontario and western Quebec. Journal of General Virology 75: 2575–2583.[Abstract/Free Full Text]

———, M. ABDEL-MALIK, J. ARMSTRONG, AND A. I. WANDELER. 2002. Lyssavirus P gene characterization provides insight into the phylogeny of the genus and identifies structural similarities and diversity within the encoded phosphoprotein. Virology 298: 286–305.[Medline]

NASYMOVICH, F. F., AND U. A. ISAKOV. 1985. Arctic fox, red fox, raccoon dog. Nauka, Moscow, USSR, pp. 1–160. [In Russian.]

NAZAROV, V. P. 1961. Animal rabies. Selhozgiz, Moscow, USSR, 160 pp. [In Russian.]

NYBERG, M., K. KULONEN, E. NEUVONEN, AND C. EK-KOMMONEN. 1992. An epidemic of sylvatic rabies in Finland. Descriptive epidemiology and results of oral vaccination. Acta Veterinaria Scandinavian 33: 43–57.

PASTORET, P. P., AND B. BROCHIER. 1999. Epidemiology and control of fox rabies in Europe. Vaccine 17: 1750–1754.[Medline]

POSADA, D., AND K. A. CRANDALL. 1998. MODEL-TEST: Testing the model of DNA substitution. Bioinformatics 14: 817–818.[Abstract/Free Full Text]

PREHAUD, C., S. LAY, B. DIETZSCHOLD, AND M. LA-FON. 2003. Glycoprotein of nonpathogenic rabies viruses is a key determinant of human cell apoptosis. Journal of Virology 77: 10537–10547.[Abstract/Free Full Text]

RYBAK, I. V., G. N. SIDOROV, AND I. V. KUZMIN. 1992. Properties of an existence and distribution of the rabies epizootics in Omsk province during 1962–1991. In Proceedings of the International Symposium "The problems of wild animals pathology and health protection, ecological interactions of wild and agricultural animal diseases," Astrachan, Moscow, Russia, 5–9 October 1992, pp. 37–39. [In Russian.]

SANSIZBAEV, B. K. 1975. Some epizootological and epidemiological aspects of rabies in Kazakchstan. Journal of Microbiology, Epidemiology and Immunobiology 5: 139–140. [In Russian.]

SAVITSKY, V. P., A. D. BOTVINKIN, V. I. BELKO, S. P. MAJOROV, N. F. SIDELNIKOVA, AND L. E. GORKOVENKO. 1981. Epidemiological peculiarities of rabies in Far East. In Modern methods in the study of natural foci infections. Institute for Natural Foci Infections, Omsk, USSR, pp. 31–40. [In Russian.]

SELIMOV, M. A. 1978. Rabies. Medicine, Moscow, USSR, 336 pp. [In Russian.]

———. 1998. Modern epizootic situation and perspectives of rabies elimination. Voprosi Virusologii 5: 196–198. [In Russian.]

———, A. A. KING, E. V. KLUEVA, AND A. D. BOTVINKIN. 1988. Antigenic variation in rabies viruses of the USSR. Rabies Information Exchange 17: 40–42.

———, A. D. BOTVINKIN, V. V. KHOZINSKI, AND L. YA. GRIBANOVA. 1994. New data about distribution of P-41 positive strains of rabies virus in arctic and nonarctic regions. Journal of Microbiology, Epidemiology and Immunobiology 2: 53–56. [In Russian.]

SLUDSKY, A. A. 1954. The role of wild mammals in distribution of infectious diseases of domestic animals. In Natural foci of infectious diseases in Kazakhstan. Medicine, Alma-Ata, USSR, pp. 67–98. [In Russian.]

SMITH, J. S. 2002. Molecular epidemiology. In Rabies, A. C. Jackson and W. H. Wunner (eds.). Academic Press, New York, pp. 79–113.

SUREAU, P., P. ROLLIN, AND T. J. WIKTOR. 1983. Epidemiologic analysis of antigenic variations of street rabies virus: Detection by monoclonal antibodies. American Journal of Epidemiology 117: 605–609.[Abstract/Free Full Text]

SWOFFORD, D. L. 2000. Phylogenetic analysis using parsimony (*and other methods), Version 4.0b6. Sinauer Associates, Sunderland, Massachusetts.

TOMA, B., AND L. ANDRAL. 1977. Epidemiology of fox rabies. Advances in Virus Research 21: 1–36.[Medline]

TORDO, N., O. POCH, A. ERMINE, G. KEITH, AND F. ROUGEON. 1986. Walking along the rabies genome: Is the large G-L intergenic region a remnant gene? Proceedings of National Academy of Sciences USA 83: 3914–3918.[Abstract/Free Full Text]

———, D. SACRAMENTO, AND H. BOURHY. 1996. The polymerase chain reaction (PCR) technique for diagnosis, typing and epidemiological studies of rabies. In Laboratory techniques in rabies, 4th Edition, F.-X. Meslin, M. M. Kaplan, and H. Koprowski (eds.). World Health Organization, Geneva, Switzerland, pp. 157–174.

VEDERNIKOV, V. A., A. A. SHABEIKIN, A. A. KHAR-KEVICH, A. M. GULUKIN, N. A. KLEMENTEVA, V. A. SEDOV, AND S. A. KOLOMITSEV. 2002. A review on the rabies epizootic situation in Russian Federation in 2000 and its prediction for 2001. Veterinary Pathology 1: 52–58. [In Russian.]

VYAZHEVICH, V. K. 1959. The case of human rabies caused by wild fox. Journal of Microbiology, Epidemiology and Immunobiology 10: 146–147. [In Russian.]

WILDE, H. 1997. Rabies, 1996. International Journal of Infectious Diseases 1: 135–142.

WOELK, C. H., O. G. PYBUS, L. JIN, D. W. G. BROWN, AND E. C. HOLMES. 2002. Increased positive selection pressure in persistent (SSPE) versus acute measles virus infections. Journal of General Virology 83: 1419–1430.[Abstract/Free Full Text]

YANG, Z. 1997. PAML: A program package for phylogenetic analysis by maximum likelihood. Computer Applications in the Biosciencies 13: 555–556.

———, NIELSEN, N. GOLDMAN, AND A. M. PEDERSEN. 2000. Codonsubstitution models for heterogeneous selection pressure at amino acid sites. Genetics 155: 431–449.[Abstract/Free Full Text]

YANOVICH, V. A. 2003. Epidemiological situation and prophylaxis of rabies in Jewish autonomous region. All-Russian Medical Scientific-Practical Journal 3: 32. [In Russian.]

Received for publication 3 November 2003.



This article has been cited by other articles:


Home page
J. Gen. Virol.Home page
H. Bourhy, J.-M. Reynes, E. J. Dunham, L. Dacheux, F. Larrous, V. T. Q. Huong, G. Xu, J. Yan, M. E. G. Miranda, and E. C. Holmes
The origin and phylogeography of dog rabies virus
J. Gen. Virol., November 1, 2008; 89(11): 2673 - 2681.
[Abstract] [Full Text] [PDF]


Home page
J. Gen. Virol.Home page
D. David, G. J. Hughes, B. A. Yakobson, I. Davidson, H. Un, O. Aylan, I. V. Kuzmin, and C. E. Rupprecht
Identification of novel canine rabies virus clades in the Middle East and North Africa
J. Gen. Virol., March 1, 2007; 88(3): 967 - 980.
[Abstract] [Full Text] [PDF]


Home page
J. Gen. Virol.Home page
I. V. Kuzmin, G. J. Hughes, and C. E. Rupprecht
Phylogenetic relationships of seven previously unclassified viruses within the family Rhabdoviridae using partial nucleoprotein gene sequences.
J. Gen. Virol., August 1, 2006; 87(Pt 8): 2323 - 2331.
[Abstract] [Full Text] [PDF]


Home page
J. Gen. Virol.Home page
G. J. Hughes, L. A. Orciari, and C. E. Rupprecht
Evolutionary timescale of rabies virus adaptation to North American bats inferred from the substitution rate of the nucleoprotein gene
J. Gen. Virol., May 1, 2005; 86(5): 1467 - 1474.
[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
Services
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 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.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS