<figure><img src="image_1.png" /></figure>

HAL
open science

# Genomic Diversity and Evolution of the Lyssaviruses

Olivier Delmas, Edward Holmes, Chiraz Talbi, Florence Larrous, Laurent Dacheux, Christiane Bouchier, Hervé Bourhy

► To cite this version:

Olivier Delmas, Edward Holmes, Chiraz Talbi, Florence Larrous, Laurent Dacheux, et al.. Genomic Diversity and Evolution of the Lyssaviruses. PLoS ONE, 2008, 3 (4), pp.e2057. 10.1371/journal.pone.0002057. pasteur-03695822

HAL Id: pasteur-03695822

https://pasteur.hal.science/pasteur-03695822

Submitted on 15 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Genomic Diversity and Evolution of the Lyssaviruses

Olivier Delmas¹, Edward C. Holmes²,⁴, Chiraz Talbi¹, Florence Larrous¹, Laurent Dacheux¹, Christiane Bouchier³, Hervé Bourhy¹*

1 Institut Pasteur, UPRE Lyssavirus Dynamics and Host Adaptation, World Health Organization Collaborating Centre for Reference and Research on Rabies, Paris, France,
2 Mueller Laboratory, Center for Infectious Disease Dynamics, Department of Biology, The Pennsylvania State University, University Park, Pennsylvania, United States of America, 3 Institut Pasteur, Plate-forme Génomique - Pasteur Genopole® Ile de France, Paris, France, 4 Fogarty International Center, National Institutes of Health, Bethesda, Maryland, United States of America

## Abstract

Lyssaviruses are RNA viruses with single-strand, negative-sense genomes responsible for rabies-like diseases in mammals. To date, genomic and evolutionary studies have most often utilized partial genome sequences, particularly of the nucleoprotein and glycoprotein genes, with little consideration of genome-scale evolution. Herein, we report the first genomic and evolutionary analysis using complete genome sequences of all recognised lyssavirus genotypes, including 14 new complete genomes of field isolates from 6 genotypes and one genotype that is completely sequenced for the first time. In doing so we significantly increase the extent of genome sequence data available for these important viruses. Our analysis of these genome sequence data reveals that all lyssaviruses have the same genomic organization. A phylogenetic analysis reveals strong geographical structuring, with the greatest genetic diversity in Africa, and an independent origin for the two known genotypes that infect European bats. We also suggest that multiple genotypes may exist within the diversity of viruses currently classified as 'Lagos Bat'. In sum, we show that rigorous phylogenetic techniques based on full length genome sequence provide the best discriminatory power for genotype classification within the lyssaviruses.

**Citation:** Delmas O, Holmes EC, Talbi C, Larrous F, Dacheux L, et al. (2008) Genomic Diversity and Evolution of the Lyssaviruses. PLoS ONE 3(4): e2057.
doi:10.1371/journal.pone.0002057

**Editor:** Oliver G. Pybus, University of Oxford, United Kingdom

**Received** January 28, 2008; **Accepted** March 17, 2008; **Published** April 30, 2008

**Copyright:** © 2008 Delmas et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** O. Delmas was supported by the European commission through the "VIZIER" Integrated project (LSHG-CT-2004-511966). This work was also funded by the VIZIER program, and by the grant of Genomic Platform Ile de France "Analyse de la Diversité Génétique des Rhabdoviridae" 2005-2006, and in part by NIH grant number GM080533-01.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: hbourhy@pasteur.fr

## Introduction

Lyssaviruses (LYSSAV) are RNA viruses with single-stranded, negative-sense genomes of the family *Rhabdoviridae* [1] that infect a variety of mammals causing rabies-like diseases. Rabies is an ancient disease that may have been reported in the Old World before 2300 B.C. [2]. However, the absence of effective control measures in animal reservoir populations combined with a widespread lack of human access to vaccination means that more than 50,000 people annually die of rabies, particularly in Asia and Africa [3,4]. Currently, there are seven recognised genotypes (GT) of LYSSAV defined on the basis of their genetic similarity [5,6]: rabies virus (RABV, GT1) responsible for classical rabies in terrestrial mammals globally and in bats on the American continent, as well as the cause of most rabies-related human deaths worldwide [3]; Lagos bat virus (LBV, GT2); Mokola virus (MOKV, GT3); Duvenhage virus (DUVV, GT4); European bat lyssavirus type 1 (EBLV-1, GT5); European bat lyssavirus type 2 (EBLV-2, GT 6); and Australian bat lyssavirus (ABLV, GT7). All genotypes except MOKV (where the host species is unknown) have bat reservoirs, hinting that lyssaviruses originated in these mammals [7]. Additionally, four new lyssavirus genotypes that infect bats in central and southeast Asia have been proposed: Aravan virus, Khujand virus, Irkut virus and West Caucasian Bat virus [8,9]. The negative-sense LYSSAV genome encodes five proteins: the nucleoprotein (N), phosphoprotein (P), matrix protein

(M), glycoprotein (G) and RNA polymerase (L) in the order 3'-N-P-M-G-L-5' [10].

Despite the importance of LYSSAV for human and wildlife populations, the number of complete genome sequences of field isolates of LYSSAV is sparse, with only eight currently available for limited type species [11–15]. Herein, we present the first genomic and evolutionary analysis of the seven known genotypes of LYSSAV, therein significantly increasing the extent of available genome sequence data available for these important mammalian pathogens.

## Materials and Methods

### Viruses and RNA isolation

Total RNA (Table 1) was isolated from original specimens or from suckling mice brain after early passage using Tri-Reagent (Euromedex). The only exception was the 8743THA isolate that was adapted on BSR cells (passage 22). For this isolate, total RNA was isolated from infected BSR cells infected at a low multiplicity of infection (0,1). Reverse transcription was performed with random hexamer primer (Roche Boehringer) using Superscript II (Invitrogen) following the manufacturer instructions.

### PCR and sequence determination

Long-range PCR products were obtained using ExTaq (Takara) and specific primers (Table S1) using manufacturer recommen-

<figure><img src="image_3.png" /></figure>
<table><caption>**Table 1.** Isolates of lyssavirus analysed in this study.</caption><thead><tr><th>Genus and name</th><th>Reference no.</th><th>Host species/vector</th><th>Origin</th><th>Year of first isolation</th><th>GenBank accession no.</th></tr></thead><tbody><tr><td colspan="6"><strong>Lyssavirus Genotype 1</strong></td></tr><tr><td>Rabies virus</td><td>8743THA</td><td>Human</td><td>Thailand</td><td>1983</td><td>EU293121</td></tr><tr><td>Rabies virus</td><td>8764THA</td><td>Human</td><td>Thailand</td><td>1983</td><td>EU293111</td></tr><tr><td>Rabies virus</td><td>9147FRA</td><td>Fox</td><td>France</td><td>1991</td><td>EU293115</td></tr><tr><td>Rabies virus</td><td>9001FRA</td><td>Dog bitten by a bat</td><td>French Guyana</td><td>1990</td><td>EU293113</td></tr><tr><td>Rabies virus</td><td>9704ARG</td><td>Bat (Tadarida brasiliensis)</td><td>Argentina</td><td>1997</td><td>EU293116</td></tr><tr><td>Rabies virus</td><td>SHBRV-18</td><td>Bat (L. noctivagans)</td><td>USA</td><td>1983</td><td>AY705373</td></tr><tr><td>Rabies virus</td><td>NNV-RAB-H</td><td>Human</td><td>India</td><td>2006</td><td>EF437215</td></tr><tr><td>Rabies virus</td><td>RABV</td><td>Human</td><td>India</td><td>2004</td><td>AY956319</td></tr><tr><td>Rabies virus</td><td>SADB19</td><td>Vaccine</td><td></td><td></td><td>M31046</td></tr><tr><td>Rabies virus</td><td>PV</td><td>Vaccine</td><td></td><td></td><td>NC_001542</td></tr><tr><td colspan="6"><strong>Genotype 2</strong></td></tr><tr><td>Lagos bat virus</td><td>8619NGA</td><td>Bat (Eidolon helvum)</td><td>Nigeria</td><td>1956</td><td>EU293110</td></tr><tr><td>Lagos bat virus</td><td>0406SEN</td><td>Bat (Eidolon helvum)</td><td>Senegal</td><td>1985</td><td>EU293108</td></tr><tr><td colspan="6"><strong>Genotype 3</strong></td></tr><tr><td>Mokola virus</td><td>MOKV</td><td>Cat</td><td>Zimbabwe</td><td>1981</td><td>NC_006429</td></tr><tr><td>Mokola virus</td><td>86100CAM</td><td>Shrew</td><td>Cameroun</td><td>1974</td><td>EU293117</td></tr><tr><td>Mokola virus</td><td>86101RCA</td><td>Rodent</td><td>RCA</td><td>1981</td><td>EU293118</td></tr><tr><td colspan="6"><strong>Genotype 4</strong></td></tr><tr><td>Duvenhage virus</td><td>94286SA</td><td>Bat (Minopterus sp.)</td><td>South Africa</td><td>1981</td><td>EU293120</td></tr><tr><td>Duvenhage virus</td><td>86132SA</td><td>Human</td><td>South Africa</td><td>1971</td><td>EU293119</td></tr><tr><td colspan="6"><strong>Genotype 5</strong></td></tr><tr><td>European bat lyssavirus 1</td><td>8918FRA</td><td>Bat (Eptesicus serotinus)</td><td>France</td><td>1989</td><td>EU293112</td></tr><tr><td>European bat lyssavirus 1</td><td>03002FRA</td><td>Bat (Eptesicus serotinus)</td><td>France</td><td>2003</td><td>EU293109</td></tr><tr><td>European bat lyssavirus 1</td><td>RV9</td><td>Bat (Eptesicus serotinus)</td><td>Germany</td><td>1968</td><td>EF157976</td></tr><tr><td colspan="6"><strong>Genotype 6</strong></td></tr><tr><td>European lyssavirus 2</td><td>9018HOL</td><td>Bat (Myotis dasycneme)</td><td>Holland</td><td>1986</td><td>EU293114</td></tr><tr><td>European lyssavirus 2</td><td>RV1333</td><td>Human</td><td>Scotland</td><td>2002</td><td>EF157977</td></tr><tr><td colspan="6"><strong>Genotype 7</strong></td></tr><tr><td>Australian bat lyssavirus</td><td>ABLh</td><td>Human</td><td>Australia</td><td>1986</td><td>AF418014</td></tr><tr><td>Australian bat lyssavirus</td><td>ABLb</td><td>Bat (Pteropus species)</td><td>Australia</td><td>1996</td><td>NC_003243</td></tr></tbody></table>

Genbank accession numbers for the newly acquired sequences are designated EU293108-EU293121.

doi:10.1371/journal.pone.0002057.t001

For sequence determination we used a shotgun base approach called LoPPS (Long PCR Product Sequencing) [16,17]. 3' genomic ends were generated by RACE protocol [14] using a 5' phosphorylated reverse complementary T7 primer. T7 cDNAs were further used for heminested-PCR with ExTaq using T7 and two strain specific primers designed in the N coding region (supplementary Table 1). To determine the 5' sequence of the genomic RNA we used a 5'RACE version 2.0 kit from Invitrogen following manufacturer instructions. The PCR products (5' or 3' RACE) were then purified on gel using Qiaquick gel extraction kit (Qiagen) and cloned in PCR 2.1 TOPO T/A (Invitrogen) for sequencing. Each position of the consensus nucleotide sequence was determined from at least three independent sequences. All consensus sequences obtained using Sequencher 4.7 (Gene Codes) software were aligned using ClustalX 1.83.1 [18]. The untranslated regions were further aligned manually using the SE-AL program (http://tree.bio.ed.ac.uk/). GenBank accession numbers for the sequences newly acquired here are designated EU293108-EU293121.

## Phylogenetic analysis

Phylogenetic analysis of LYSSAV genomes was based on a multiple alignment of concatenated coding region sequences (12105 nt). A maximum likelihood (ML) phylogenetic analysis of these data was undertaken using PAUP* [19] employing the best-fit GTR+I+γ₄ model of nucleotide substitution inferred by ModelTest [20]. To determine the extent of support for different groupings on the tree a bootstrap resampling analysis was undertaken employing 1000 replicate neighbor-joining trees estimated under the ML substitution model.

## Results and Discussion

In total we determined 14 new complete genome sequences of field isolates representing six (GT1, GT2, GT3, GT4, GT5 and GT6) of the seven genotypes of LYSSAV, with complete genome from GT4 obtained for the first time. These genomes were combined with eight genomes described previously (with the exception of one Australian bat lyssavirus for which leader and

<figure><img src="image_4.png" /></figure>
<table><caption>**Table 2.** Coding potential, genome size (in nucleotides) and G+C content of 24 genomes representing the 7 genotypes of the lyssavirus genus.</caption><thead><tr><th>Genotype</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>All genotypes</th></tr></thead><tbody><tr><td>3' UTR</td><td>70</td><td>70</td><td>70</td><td>70</td><td>70</td><td>70</td><td>70</td><td>70</td></tr><tr><td>N protein</td><td>1353</td><td>1353</td><td>1353</td><td>1356</td><td>1356</td><td>1356</td><td>1353</td><td>1353-1356</td></tr><tr><td>N-P</td><td>90-94</td><td>101</td><td>100-102</td><td>90</td><td>90-96</td><td>101</td><td>93-94</td><td>90-102</td></tr><tr><td>P protein</td><td>894</td><td>918</td><td>912</td><td>897</td><td>894</td><td>894</td><td>894</td><td>894-918</td></tr><tr><td>P-M</td><td>87-90</td><td>75</td><td>80-83</td><td>83</td><td>83</td><td>88</td><td>87</td><td>75-90</td></tr><tr><td>M protein</td><td>609</td><td>609</td><td>609</td><td>609</td><td>609</td><td>609</td><td>609</td><td>609</td></tr><tr><td>M-G</td><td>211-215</td><td>204</td><td>203-204</td><td>191</td><td>211</td><td>205-210</td><td>207-209</td><td>191-215</td></tr><tr><td>G protein</td><td>1575</td><td>1569</td><td>1569</td><td>1602</td><td>1575</td><td>1575</td><td>1578-1581</td><td>1569-1602</td></tr><tr><td>G-L</td><td>515-525</td><td>578-588</td><td>546-562</td><td>562-563</td><td>560</td><td>511-512</td><td>508-509</td><td>508-588</td></tr><tr><td>L protein</td><td>6384-6429</td><td>6384</td><td>6381-6384</td><td>6384</td><td>6384</td><td>6384</td><td>6384-6387</td><td>6381-6429</td></tr><tr><td>5' UTR</td><td>86-131</td><td>145</td><td>112-114</td><td>131</td><td>130-131</td><td>131</td><td>131</td><td>86-145</td></tr><tr><td>Genome</td><td>11923-11928</td><td>12006-12016</td><td>11940-11957</td><td>11975-11976</td><td>11966-11971</td><td>11924-11930</td><td>11918</td><td>11918-12016</td></tr><tr><td>G+C%</td><td>44.9-45.4</td><td>40.9-43.5</td><td>44.1-44.9</td><td>44.1-44.2</td><td>44.6-45.0</td><td>44.8</td><td>43.4-44.2</td><td>40.9-45.4</td></tr></tbody></table>

Upper and lower size ranges are indicated.
doi:10.1371/journal.pone.0002057.t002

trailer sequences are unavailable). Eight field isolates of viruses isolated from humans, canids and bats were chosen as representative of the diversity of GT1. Two vaccine strains (SAD-B19 and PV) were included in all sequence comparisons but not in the phylogenetic analysis. Our study also represent the first analysis of the intrinsic genetic diversity of GT2, GT3, GT4, GT5 and GT6 based on full length genomes.

All genomes have the same structural organization although their lengths varied between 11918 nt. (GT7) and 12016 nt. (GT2) (Table 2). The predicted size of the coding regions is similar among genotypes, with the M protein identical in length across all genotypes and the P protein the most variable [14,21,22]. As observed in other RNA viruses, all genotypes show a bias toward G+C richness [23], with the lowest G+C content observed in GT2 and the highest in GT1 (Table 2). All genomes have a polycistronic genome organization surrounded by untranslated regions (Table S2) similar to that already described [10,14]. The extent of genetic diversity, reflected in percentage identity, varies within and among proteins (Figure 1), in the order N>L>M>G>P (95.2, 94.2, 92.3, 85.8, 81.5% amino acid identity, respectively). A similar pattern was previously observed using more limited data sets [14]. This same order was also observed in terms of overall selection pressure, measured as the mean ratio of nonsynonymous ($d_N$) to synonymous substitutions ($d_S$) per site ($d_N/d_S$), estimated using the maximum likelihood SLAC (Single Likelihood Ancestor Counting; http://www.datamonkey.org/) method [24]: $N=0.048$; $L=0.055$; $M=0.078$; $G=0.119$; $P=0.187$. This approximately four-fold difference in mean $d_N/d_S$ reflects major differences in selective constraint among proteins. This trend was also reflected in previous analyses of full length genomes of vaccine strains [22] and through partial gene comparisons [25,26].

Our study represents the largest analysis of the 3' and 5' UTR of the lyssavirus genomes undertaken to date. The 3' UTR comprises 70 nt and includes the leader regions potentially transcribed into the leader RNA. The 5' UTR region comprises 86-145 nt and contains the trailer regions of size 68-69 nt. Both the 3' and 5' UTR have conserved signals that play a role to modulate replication and transcription (Figure S1) [27]. Our data also reveals a strict complementarity limited to the 9 terminal

nucleotides as well as nucleotide positions 14 and 16 from both ends of the genome [14,28,29].

There have been several attempts to estimate the evolutionary relationships among lyssaviruses, with most utilizing only one or two genes [1,7,21,26,30-34]. We therefore undertook a phylogenetic analysis of 22 genomes representative of the seven genotypes of LYSSAV based on a multiple alignment of concatenated coding sequences. Our phylogenetic analysis reveals the separation of LYSSAV into two major branches previously defined as different 'phylogroups' [7] and 7 component lineages defined as genotypes [5,35]. Phylogroup 1 comprised GT1, 4, 5, 6 and 7, while phylogroup 2 contains only GT2 and GT3 (Figure 2). Notably, phylogroup 2 contains viruses of sampled exclusively from Africa – LBV and MOKV – while a third African genotype (DUVV) is found within phylogroup 1 [32,36]. Also of note was the observation that although GT5 and GT6 both circulate in European insectivorous bats [32], the former is more closely related to the African GT4 viruses [32]. Hence, there has clearly been an independent origin of genotypes 5 and 6 in European bats, as previously documented in analyses of the N and G genes in isolation [32]. Finally, that bats appear as the principle host species across such a large phylogeographic range indicates that the association between lyssaviruses and bats is likely to be the ancestral condition (with a secondary loss of bat transmission in GT3), such that the movement of bats is likely to be responsible for the global dissemination of these viruses [7].

Notably, our study represents the first analysis of the genetic diversity of four complete genomes of GT2 and GT4, both of which are African in origin. While little variation is seen within GT4, the degree of divergence among the two GT2 isolates (Lagos bat virus – LBV) is striking (23.7% and 12.1% at the nucleotide and amino acid levels, respectively) and greater than that seen within any other genotype. Hence, although 0406SEN and 8619NGA are related according to the arbitrary classification system based on nucleotide identity between N coding regions (80.3% between 0406SEN and 8619NGA compared to a cut-off of 80%) [6,21], this classification system will likely need to be revised as expanded surveys of LYSSAV in Africa (this study and [37]) and in Eurasia [8,9,21] reveal greater genetic diversity. More

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 1.** Schematic representation of lyssavirus genome organization and sequence similarity among 24 aligned genomes. A. The 3' leader, N-, P-, M-, G- and L-coding regions and the 5' trailer region are shown. B. Sequence similarity is calculated by moving a window of 60 nucleotides along the aligned sequences. C. Sequence similarity is calculated by moving a window of 20 amino acids along the aligned sequences. Within each window, the similarity of any one position is taken to be the average of all the possible pairwise scores at that position and is calculated using PLOTCON (available at http://bioweb.pasteur.fr/seqanal/interfaces/plotcon.html). doi:10.1371/journal.pone.0002057.g001</figcaption></figure>

fundamentally, if, as we suggest, complete genomes represent the best tools for genotyping, we propose that 0406SEN should constitute a new GT8 different from GT2 (8619NGA) and that the genotype division should be set at 76.4 to 81.6% nucleotide identity at coding sequences for all five viral proteins. Such a cut-

off would provide more discriminatory power than systems that utilize the N gene in isolation (Table 3).

Finally, we suggest that the phylogenetic methods used here – based on a realistic model of nucleotide substitution, a robust phylogenetic method, and rigorous bootstrap resampling –

<figure><img src="image_7.png" /><figcaption>**Figure 2.** Phylogenetic relationships of 22 complete coding regions of LYSSAV genomes representatives of the 7 genotypes. The phylogeny was inferred using an ML procedure, and all horizontal branches are scaled according to the number of substitutions per site. Boot strap values (>95%) are shown for key nodes. The tree is mid-point rooted for purposes of clarity only. doi:10.1371/journal.pone.0002057.g002</figcaption></figure>

<figure><img src="image_8.png" /></figure>
<table><caption>**Table 3.** Minimum intra-genotype and maximum inter-genotype sequence similarities among 24 lyssaviruses.</caption><thead><tr><th>Coding regions</th><th>Number of genotypes</th><th>Minimum intragenotype similarity</th><th>Maximum intergenotype similarity</th></tr></thead><tbody><tr><td rowspan="2">N</td><td>8*</td><td>83,3</td><td>80,3</td></tr><tr><td>7</td><td>80,3</td><td>79,8</td></tr><tr><td rowspan="2">N,P,M,G,L</td><td>8*</td><td>81,6</td><td>76,4</td></tr><tr><td>7</td><td>76,3</td><td>76,4</td></tr></tbody></table>

*when 04065EN is considered as the representative isolate of a new GT8, with 8619 the representative isolate of GT2.
doi:10.1371/journal.pone.0002057.t003

represent a more powerful method of lyssavirus classification than those based on pairwise genetic diversity alone, particularly as they account for any lineage-specific rate variation that will compromise all distance-based approaches used to date. This method has also been proposed for HIV to try to standardize viral classification [38] confirming the interest of this method for viral classification.

## Supporting Information

**Table S1** List of primers

Found at: doi:10.1371/journal.pone.0002057.s001 (0.06 MB PDF)

**Table S2** Transcription and termination signals for all lyssavirus genotypes.

Found at: doi:10.1371/journal.pone.0002057.s002 (0.05 MB PDF)

**Figure S1** Comparison of the 5' and the reverse complementary 3' genomic termini of the antigenomic (+) sense RNA of

lyssaviruses. Identical nucleotides are indicated by a vertical line. A, 23 lyssaviruses representing the 7 genotypes. B, consensus sequences. Only regions corresponding to 3' and 5' UTR sequences are shown. TTS: transcription termination signal.

Found at: doi:10.1371/journal.pone.0002057.s003 (0.27 MB DOC)

## Acknowledgments

We are grateful to Dr Alma Palazollo for providing isolate 9704ARG from Argentina and to Magali Tichit for its technical help.

## Author Contributions

Conceived and designed the experiments: EH HB LD OD. Performed the experiments: EH CB HB FL OD CT. Analyzed the data: EH HB FL OD CT. Contributed reagents/materials/analysis tools: EH CB HB LD OD. Wrote the paper: EH HB OD.

## References

1. Bourhy H, Cowley JA, Larrous F, Holmes EC, Walker PJ (2005) Phylogenetic relationships among rhabdoviruses inferred using the L polymerase gene. *J Gen Virol* 86: 2849–2858.

2. Steele JH, Fernandez PJ (1991) History of rabies and global aspects. In: Baer GM, ed. The Natural History of Rabies. 2nd edition ed. Boca Raton, USA: CRC Press. pp 1–26.

3. Knobel DL, Cleaveland S, Coleman PG, Fevre EM, Meltzer MI, et al. (2005) Re-evaluating the burden of rabies in Africa and Asia. Bull World Health Organ 83: 360–368.

4. Warrell MJ, Warrell DA (2004) Rabies and other lyssavirus diseases. Lancet 363: 959–969.

5. Gould AR, Hyatt AD, Lunt R, Kattenbelt JA, Hengstberger S, et al. (1998) Characterisation of a novel lyssavirus isolated from Pteropid bats in Australia. Virus Res 54: 165–187.

6. Kissi B, Tordo N, Bourhy H (1995) Genetic polymorphism in the rabies virus nucleoprotein gene. Virology 209: 526–537.

7. Badrane H, Tordo N (2001) Host switching in Lyssavirus history from the Chiroptera to the Carnivora orders. *J Virol* 75: 8096–8104.

8. Arai YT, Kuzmin IV, Kameoka Y, Botvinkin AD (2003) New lyssavirus genotype from the Lesser Mouse-eared Bat (Myotis blythi), Kyrghyzstan. Emerg Infect Dis 9: 333–337.

9. Botvinkin AD, Poleschuk EM, Kuzmin IV, Borisova TI, Gazaryan SV, et al. (2003) Novel lyssaviruses isolated from bats in Russia. Emerg Infect Dis 9: 1623–1625.

10. Tordo N, Poch O, Ermine A, Keith G, Rougeon F (1986) Walking along the rabies genome: is the large G-L intergenic region a remnant gene? Proc Natl Acad Sci U S A 83: 3914–3918.

11. Faber M, Pulmanausahakul R, Nagao K, Prosniak M, Rice AB, et al. (2004) Identification of viral genomic elements responsible for rabies virus neuroinvasiveness. Proc Natl Acad Sci U S A 101: 16328–16332.

12. Gould AR, Kattenbelt JA, Gumley SG, Lunt RA (2002) Characterisation of an Australian bat lyssavirus variant isolated from an insectivorous bat. Virus Res 89: 1–28.

13. Le Mercier P, Jacob Y, Tordo N (1997) The complete Mokola virus genome sequence: structure of the RNA-dependent RNA polymerase. *J Gen Virol* 78 (Pt 7): 1571–1576.

14. Marston DA, McElhinney LM, Johnson N, Muller T, Conzelmann KK, et al. (2007) Comparative analysis of the full genome sequence of European bat lyssavirus type 1 and type 2 with other lyssaviruses and evidence for a conserved

transcription termination and polyadenylation motif in the G-L 3' non-translated region. *J Gen Virol* 88: 1302–1314.

15. Warrilow D, Smith IL, Harrower B, Smith GA (2002) Sequence analysis of an isolate from a fatal human infection of Australian bat lyssavirus. Virology 297: 109–119.

16. Emonet S, Grard G, Brisbarre N, Moureau G, Temmam S, et al. (2006) LoPPS: a long PCR product sequencing method for rapid characterisation of long amplicons. Biochem Biophys Res Commun 344: 1080–1085.

17. Emonet SF, Grard G, Brisbarre NM, Moureau GN, Temmam S, et al. (2007) Long PCR Product Sequencing (LoPPS): a shotgun-based approach to sequence long PCR products. Nat Protoc 2: 340–346.

18. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25: 4876–4882.

19. Swofford DL (2003) PAUP*. Phylogenetic Analysis Using Parsimony (*and other methods) version 4. Sunderland MA: Sinauer Associates.

20. Posada D, Crandall KA (1998) MODELTEST: testing the model of DNA substitution. Bioinformatics 14: 817–818.

21. Kuzmin IV, Hughes GJ, Botvinkin AD, Orciari LA, Rupprecht CE (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.

22. Wu X, Franka R, Velasco-Villa A, Rupprecht CE (2007) Are all lyssavirus genes equal for phylogenetic analyses? Virus Res 129: 91–103.

23. Auewarakul P (2005) Composition bias and genome polarity of RNA viruses. Virus Res 109: 33–37.

24. Pond SL, Frost SD (2005) Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics 21: 2531–2533.

25. Kuzmin IV, Orciari LA, Arai YT, Smith JS, Hanlon CA, et al. (2003) Bat lyssaviruses (Aravan and Khujand) from Central Asia: phylogenetic relationships according to N, P and G gene sequences. Virus Res 97: 65–79.

26. Nadin-Davis SA, Abdel-Malik M, Armstrong J, Wandeler AI (2002) Lyssavirus P gene characterisation provides insights into the phylogeny of the genus and identifies structural similarities and diversity within the encoded phosphoprotein. Virology 298: 286–305.

27. Finke S, Conzelmann KK (1999) Virus promoters determine interference by defective RNAs: selective amplification of mini-RNA vectors and rescue from cDNA by a 3' copy-back ambisense rabies virus. *J Virol* 73: 3818–3825.

<figure><img src="image_9.png" /></figure>
28. Bourhy H, Tordo N, Lafon M, Sureau P (1989) Complete cloning and molecular organization of a rabies-related virus, Mokola virus. J Gen Virol 70 (Pt 8): 2063-2074.

29. Tordo N, Poch O, Ermine A, Keith G, Rougeon F (1988) Completion of the rabies virus genome sequence determination: highly conserved domains among the L (polymerase) proteins of unsegmented negative-strand RNA viruses. Virology 165: 565-576.

30. Bourhy H, Kissi B, Audry L, Smreczak M, Sadkowska-Todys M, et al. (1999) Ecology and evolution of rabies virus in Europe. J Gen Virol 80 (Pt 10): 2545-2557.

31. Davis PL, Bourhy H, Holmes EC (2006) The evolutionary history and dynamics of bat rabies virus. Infect Genet Evol 6: 464-473.

32. Davis PL, Holmes EC, Larrous F, Van der Poel WH, Tjornehoj K, et al. (2005) Phylogeography, population dynamics, and molecular evolution of European bat lyssaviruses. J Virol 79: 10487-10497.

33. Davis PL, Rambaut A, Bourhy H, Holmes EC (2007) The evolutionary dynamics of canid and mongoose rabies virus in Southern Africa. Arch Virol 152: 1251-1258.

34. Holmes EC, Woelk CH, Kassis R, Bourhy H (2002) Genetic constraints and the adaptive evolution of rabies virus in nature. Virology 292: 247-257.

35. Bourhy H, Kissi B, Tordo N (1993) Molecular diversity of the Lyssavirus genus. Virology 194: 70-81.

36. Shope RE (1982) Rabies-related viruses. Yale J Biol Med 55: 271-275.

37. Kuzmin IV, Niezgoda M, Franka R, Agwanda B, Markotter W, et al. (2008) Lagos bat virus in Kenya. J Clin Microbiol.

38. Gifford R, de Oliveira T, Rambaut A, Myers RE, Gale CV, et al. (2006) Assessment of automated genotyping protocols as tools for surveillance of HIV-1 genetic diversity. Aids 20: 1521-1529.

<figure><img src="image_10.png" /></figure>