Desmodus rotundus and Artibeus spp. bats might present distinct rabies virus lineages

ORIGINAL ARTICLE

Desmodus rotundus and Artibeus spp. bats might present distinct rabies virus lineages

Willian Oliveira Fahl^I,*; Pedro Carnieli Jr^I; Juliana Galera Castilho^I; Maria Luiza Carrieri^I; Ivanete Kotait^I; Keila Iamamoto^I; Rafael Novaes Oliveira^I; Paulo Eduardo Brandão^II

^IInstituto Pasteur, São Paulo, SP, Brazil

^IIDepartment of Preventive Veterinary Medicine and Animal Health, School of Veterinary Medicine, Universidade de São Paulo, São Paulo, SP, Brazil

ABSTRACT

In Brazil, bats have been assigned an increasing importance in public health as they are important rabies reservoirs. Phylogenetic studies have shown that rabies virus (RABV) strains from frugivorous bats Artibeus spp. are closely associated to those from the vampire bat Desmodus rotundus, but little is known about the molecular diversity of RABV in Artibeus spp. The N and G genes of RABV isolated from Artibeus spp. and cattle infected by D. rotundus were sequenced, and phylogenetic trees were constructed. The N gene nucleotides tree showed three clusters: one for D. rotundus and two for Artibeus spp. Regarding putative N amino acid-trees, two clusters were formed, one for D. rotundus and another for Artibeus spp. RABV G gene phylogeny supported the distinction between D. rotundus and Artibeus spp. strains. These results show the intricate host relationship of RABV's evolutionary history, and are invaluable for the determination of RABV infection sources.

Keywords: Rabies Bat Artibeus spp.; Desmodus rotundus; Phylogeny; Molecular epidemiology

Introduction

Rabies is a worldwide neglected fatal encephalitis,^1,2 which is listed amongst the ten major infectious causes of human deaths worldwide, estimated at 55,000 per year.³ This infection is efficiently preventable by vaccination,⁴ but treatment costs of rabid or exposed patients, diagnostic procedures, and vaccines make it a significant challenge for public health systems in endemic countries.^5,6

The disease is caused by rabies virus (RABV) (family Rhabdoviridae, genus Lyssavirus),⁷ an enveloped virus with a length of 110-250 nm and a diameter of 75 nm, usually transmitted by the saliva of an infected mammal. The RABV genome has a negative-stranded non-segmented ssRNA with 11,932-nucleotides that encodes the five structural proteins N (nucleoprotein), P(phosphoprotein), M (matrix), G (envelope glycoprotein), and L (large, RNA-dependent RNA-polymerase).^8,9 Between 2004 and 2005, 62 people died in the Brazilian Amazon of rabies transmitted by Desmodus rotundus vampire bat, a primary reservoir of rabies in Latin America.^10,11 Regarding frugivorous bats of the genus Artibeus, rabies has been reported in the species A. fimbriatus, A. jamaicensis, A. lituratus, and A. planirostris.¹²Artibeus spp. bats have been assigned an increasing importance in public health, as they are considered a rabies reservoir for humans in urban areas in Brazil, which is aggravated by the increasing population of these bats and the fact that they share roosts with D. rotundus.^12-14

From 2003 to 2008, the Instituto Pasteur in Brazil tested 18,007 non-hematophagous bats for rabies, 252 of which were found to be positive. Also, from 2005 and 2007, 56 out of 160 non-hematophagous bats that tested positive for rabies were classified as A. lituratus.

Phylogenetic studies based on the N gene have suggested that RABV lineages from Artibeus sp. are not divergent from those from D. rotundus,¹⁵ all belonging to the same genic¹⁶ and antigenic variant 3 (AgV3).¹⁷ Nonetheless, these studies have provided only inconclusive results, as they were based on a very restricted sampling regarding geographic area and sample number, mainly in the case of Artibeus spp. In addition, these studies were based on a single gene sequences for phylogenetic reconstructions.

The ability to determine the source of infection and the epidemiology of rabies cycles are paramount for accurate decision-making in public health, mainly regarding vaccination strategies and animal population control. This study aimed to evaluate the possibility of distinction between RABV genetic lineages related to D. rotundus and Artibeus spp. bats based on N and G genes sequences.

Materials and methods

Rabies virus strains

Twenty Artibeus spp. RABV strains were obtained from firstpassage isolates in mice inoculated with 20% suspensions of A. lituratus and Artibeus spp. central nervous systems (CNS); 15 D. rotundus-related strains were obtained directly from naturally infected cattle brain tissues (Table 1). These 35 strains were collected in nine municipalities from São Paulo State, Southeastern Brazil (Fig. 1), between 2004 and 2005. All samples were diagnosed positive for rabies by direct immunofluorescence test (DIFT) targeted to the viral nucleoprotein.¹⁸

The nucleoprotein and the glycoprotein sequences generated in this study have been assigned the GenBank accession numbers JF682392 -JF682426 and JF682427-JF682461, respectively.

Amplification and sequencing of the nucleoprotein and glycoprotein genes

Total RNA from the 35 RABV isolates CNS samples tested and the from positive and negative controls were extracted with TRIzol^TM (Invitrogen - Carlsbad, CA, USA) method, following the manufacturer's instructions. The challenge virus standard (CVS) fixed strain of RABV isolated in mice brain and nuclease free-water were used as positive and negative controls, respectively.

Reverse transcription polymerase chain reaction (RT-PCR) to partially amplify the N and G genes was performed according to the protocol described by Carnieli et al.,¹⁹ using the primers described by Orciari et al.²⁰ f or the N gene and those described by Sato et al.²¹ for the G gene (Table 2). The PCR products were purified from the PCR reactions using the QIAquick^TM Gel Extraction Kit^TM (Qiagen - Valencia, CA, USA), according to the manufacturer's instructions. The products with nonspecific bands were purified using 1% agarose gel and the QIAquick^® kit. After the purification step, the DNA samples were visually quantified in 2% agarose gel with Low Mass DNA Ladder (Invitrogen - Carlsbad, CA, USA), following the manufacturer's instructions.

The DNA sequencing reaction mixture consisted of 4 fL of BigDye 3.1^TM (Applied Biosystems - Foster City, CA, USA), 3.2 pmol of both sense and antisense primer for each gene in separate reactions, 30-60 ng of target DNA and DNase-free water to a final reaction of 10 fL. The reaction was performed in a Mastercycler Gradient thermal cycler^TM (Eppendorf, NY, USA) with 35 cycles at 96ºC for 10s, 50 ºC for 5s, and 60ºC for 4 min, with a ramp of 1ºC/s between each temperature. Sequencing reaction products were purified with Sephadex^TMG-50 fine beads (GE Healthcare Biosciences) in 96well multiscreen HV plates. After purification, the sequences were resolved in an ABI-3130^TM Automatic Sequencer (Applied Biosystems - Foster City, CA, USA).

Phylogenetic analyses

To construct the genealogic trees, nucleotide and putative amino acids sequences for the 35 RABV strains of N and G genes were aligned by the CLUSTAL/W multiple alignment algorithm method using the BioEdit program,²² and then by manually checking the alignments for each set of aligned sequences. A score was assigned to each of the nucleotides shown on the electropherograms for each of the sequencing reactions using the online Phred application. Only positions with a Phred score > 20 were used.²³ The final sequence for each strain was obtained using the Contig Assembly Program (CAP) in BioEdit v.5.0.9²² and submitted to BLASTn for homology confirmation.

Phylogenetic trees of the RABV isolates were built using the Neighbor-joining algorithm and the maximum composite likelihood (MCL) evolutionary model implemented in Mega 4.1 (^© 1993-2008)²⁴ with 1,000 bootstrap replicates. Additionally, 38 homologous sequences recovered from Gen-Bank were included in the phylogenetic trees of the N and G genes (29 and none, respectively) and European bat Lyssavirus 1 (another specie in the Lyssavirus) was used as an outgroup.

The minimum, maximum, and mean nucleotide (with the MCL model) and amino acids (with the Poisson correction) identities for the clusters for the N and G gene sequences were calculated using Excel (^©1985-2003 Microsoft Corporation) based on the identity matrices calculated with the BioEdit program. The changes in the amino acids observed in the samples analyzed were studied using Mega 4.1 (^© 1993-2008)²⁴ and BioEdit v.7.0.0²² programs.

Results

After editing, the N gene was 1,281-nucleotides long, located between nucleotides 68 and 1,350 in relation to the N gene CVS reference strain (GenBank number AF406696) and had a putative amino acid sequence with 427 amino acids. Regarding the G gene, sequences were 1,571-nucleotides long, located between nucleotides 8 and 1,579 in relation to the G gene CVS reference strain (GenBank number FJ979833) and had a putative amino acid sequence with 520 amino acids.

The N gene nucleotide phylogenetic tree showed three clusters for the 35 RABV strains included in the study; one cluster for D. rotundus strains and two for the Artibeus spp. strains. Regarding N protein amino acids tree, the 35 strains segregated in only two clusters, one for D. rotundus and one formed mainly by Artibeus spp. strains (Fig. 2). Nucleotides and amino acids identities for the N region under analysis between D. rotundus and Artibeus spp. sequences ranged from 97.4% to 98.7% and 98.1% to 99.7%, respectively.

The G gene nucleotide phylogenetic tree showed two clusters for the 35 RABV strains included in the study; one cluster for D. rotundus strains and one for the Artibeus spp. strains. Regarding G amino acids tree, the 35 strains segregated in two clusters, one for D. rotundus and one formed by Artibeus spp. strains (Fig. 3). Nucleotides and amino acids identities for the G region under analysis between D. rotundus and Artibeus spp. sequences ranged from 97.0% to 99.1% and 96.1% to 99.0%, respectively.

Discussion

Frugivorous bats of Artibeus spp. have an emerging importance for rabies epidemiology in Brazil, mainly in urban centers, and have been reported as carriers of RABV lineages close to those found in the vampire bat D. rotundus.^12,16 In the present investigation, RABV N and G phylogenies of strains recovered from these bats showed the existence of viral lineages that can be accurately attributed to D. rotundus or Artibeus spp. bats, a previously unknown fact.

RABV lineages heterogeneity expressed phylogenetically as host-specific lineages is a widely documented epidemiological phenomenon,^16,19,25,26 and is influenced by geographic barriers rather than by species barriers only.²⁷ Accordingly, the two Artibeus spp. clusters found for N gene sequences might also represent regional sub lineages of RABV, a fact already described for D. rotundus RABV lineages in Brazil.^28,29 These observations agree with the proposition that, regarding rabies in bats, host species are as important as geographic variations.³⁰ Regarding viral lineages, genetic variations that occur within a host species are different from those that occur in another, and this variation, coupled with the host's geographical isolation, may explain the RABV genetic differences reported in this work.³¹

Seven strains (05/Art7270, 05/Art8456, 05/Art7436, 05/Art8921, 05/Art8688, 05/Art3738, and 08IacriSP3577B) were classified in the Artibeus spp. group for the N nucleotides tree; however they segregated into the D. rotundus cluster of the N amino acids tree. This fact suggests the possibility that Artibeus spp. strains are still under an adaptation process after the spill-over event from D. rotundus, as already suggested by Kissi,³¹ who have already experimentally reported this type of stepped adaptation of RABV. This phenomenon can occur as a consequence of the predominance of synonymous mutations over non-synonymous nucleotides mutations, which leads to greater differences among the nucleotide sequences than among the amino acids sequences.^32,33

Regarding the Bov7525 and Bov6314 strains, in cattle, both segregated in the Artibeus spp. RABV cluster; this was unexpected, as cattle rabies is related to that of D. rotundus and not to that of frugivorous bats.^34,35 The most plausible explanation is that these two strains have been transmitted from an Artibeus spp. to a D. rotundus, and then from this bat to the cattle in a rare class of spill-over transmission.

Kobayashi et al.,¹⁶ analyzing the N gene of the RABV isolated from frugivorous bats, insectivorous bats, and D. rotundus, reported lineages associated with Artibeus spp. bats (frugivorous), D. rotundus, and insectivorous bats, suggesting that there are species-specific lineages.²⁶ However, they have not provided significant results to distinguish between RABV isolated from D. rotundus and Artibeus spp., possibly due to a restricted number of samples and sampling area. In the present study, this problem was compensated by the inclusion of 35 RABV isolates of bats (20 from Artibeus spp. and 15 from D. rotundus) from a broad geographical area of São Paulo. In addition, this study also analyzed the complete G gene, which provides a more specific distinction between the genetic sequences.

Different hosts pose different challenges for rabies control. This is more complex in bats due to the large number of species, the different ecologic niches that they occupy, and the impossibility to vaccinate these animals. For instance, the population of D. rotundus in Latin America can be legally controlled with vampiricide anticoagulants applied on cattle or on the bats themselves.³⁴ Conversely, nonhematophagous bats are under legal protection, and only now the knowledge of rabies epidemiology in these bats species is increasing.^35,36

In this context, the results obtained in this study are valuable, because based on the partial amino acid sequences for the N gene it is possible to differentiate RABV strains from Artibeus spp. and D. rotundus for the purpose of defining the infection sources in molecular epidemiology. These results show the close host relationship of RABV transitions, and have an invaluable application for determining the sources of rabies infections transmitted mainly to dogs and cats in urban centers.

In conclusion, for rabies virus isolates related to frugivorous bats of the Artibeus spp. and to the vampire bat D. rotundus, the phylogeny based on sequences of the N and G genes shows segregation patterns in genus-specific agreement in each of these bats. Data from this study suggest that a lineage of RABV is possibly being established in the Artibeus spp. genus.

Conflict of interest

All authors declare to have no conflict of interest.

Acknowledgements

The authors are grateful to CAPES (W.O. Fahl's PhD fellowship) and CNPq (P.E. Brandão's PQ-2 fellowship).

REFERENCES

1. Fooks AR, Brookes SM, Johnson N, Mcelhinney LM, Hutson AM. European bat lyssaviruses: an emerging zoonosis. Epidem and Infect. 2003;131:1029-39.

2. Dodet B. Preventing the incurable: Asian rabies, experts advocate rabies control. Vaccine. 2006;24:3045-9.

3. World Health Organization. World survey of rabies: n. 34 for the year 1998. Geneva: WHO; 2000.

4. Briggs D, Hanlon CA. World rabies day 7: focusing attention on a neglected disease. Vet Rec. 2007;161:288-9.

5. Kaplan G, Turner GS, Warrel D. Rabies: the facts. 2^nd ed. Oxford; New York: Oxford University Press; 1986.

6. Who Expert Consultation On Rabies, 2004. Geneva: Switzerland. WHO Expert consultation on rabies: First report. Geneva: WHO; 2005. 87 p.(Technical report series, 931).

7. Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA. Virus taxonomy classification and nomenclature of viruses Eighth Report of the International Committee on the Taxonomy of Viruses. Amsterdam: Elsevier Acad Press; 2005.

8. Tordo N, Poch O, Ermine A, Keith G. Primary structure of leader RNA and nucleoprotein genes of the rabies genome: segmented homology with VSV. Nucl Assim Res. 1986;14:2671-83.

9. Tordo N, Poch O, Ermine A, Keith G, Rougeon F. Walking along the rabies genome: is the large G-L intergenic region a remnant gene? Proc Natl Acad of Sci. 1986;83:3914-8.

10. Da Rosa ES, Kotait I, Barbosa TF, et al. Bat-transmitted human rabies outbreaks, Brazilian Amazon. Emerg Infect Dis. 2006;12:1197-202.

11. Barbosa TF, Medeiros DB, Travassos da Rosa ES, et al. Molecular epidemiology of rabies virus isolated from different sources during a bat-transmitted human outbreak occurring in Augusto Correa municipality, Brazilian Amazon. Virol. 2008;370:228-36.

12. Sodré MM, Gama AR, Almeida MF. Updated list of bat species positive for rabies in Brazil. Rev Inst Med Trop. 2010;52:75-81.

13. Uieda W, Hayashi MM, Gomes LH, Silva MMS. Espécies de quirópteros diagnosticados com raiva no Brasil. Bull Inst Pasteur. 1996;1:17-35.

14. Cunha EMS, Lara MCCSH, Nassar AFC, Sodré M, Amaral LVF. Isolamento do vírus da raiva em Artibeus fimbriatus no estado de São Paulo, Brasil. Rev Saude Publica. 2005;39:683-4.

15. Shoji Y, Kobayashi Y, Sato G, et al. Genetic characterization of rabies viruses isolated from frugivorous bat (Artibeus spp.) in Brazil. J Vet Med Sci. 2004;66:1271-3.

16. Kobayashi Y, Sato G, Shoji Y, et al. Molecular epidemiological analysis of bat rabies viruses in Brazil. J Vet Med Sci. 2005;67:647-52.

17. Albas A, Souza EA, Lourenço RA, Favoretto SR, Sodré MM. Antigen profile of rabies virus isolated from different species of non-hematophagous bats in the region of Presidente Prudente, State of São Paulo. Rev Soc Bras Med Trop. 2009;42:15-7.

18. Dean DJ, Abelseth MK, Atanasiu P. The fluorescent antibody test. In: Meslin FX, Kaplan MM, Koprowski H, editors. Laboratory techniques in rabies. 4^th ed. Geneva: World Health Organization; 1996. p. 88-95.

19. Carnieli Jr P, Fahl WO, Castilho JG, et al. Characterization of rabies virus isolated from canids and identification of the main wild canid host in northeastern Brazil. Virus Res. 2008;131:33-46.

20. Orciari LA, Niezgoda M, Hanlon CA, et al. Rapid clearance of SAG-2 rabies virus from dogs after oral vaccination. Vaccine. 2001;19:4511-8.

21. Sato G, Itou T, Shoji Y, et al. Genetic and phylogenetic analysis of glycoprotein of rabies virus isolated from several species in Brazil. J Vet Med Sci. 2004;66:747-53.

22. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Series. 1999;41:95-8.

23. Ewing B, Green P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Gen Res. 1998;8:186-94.

24. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596-9.

25. Favoretto SR, Carrieri ML, Cunha EMS, et al. Antigenic typing of Brazilian rabies virus samples isolated from animals and humans, 1989-2000. Rev Inst Med Trop. 2002;44:91-5; Kobayashi Y, Sato G, Kato M, et al. Genetic diversity of bat rabies viruses in Brazil. Arch Virol. 2007;152:1995-2004.

26. Holmes EC, Woelk CH, Kassis R, Bourhy H. Genetic constraints and adaptive evolution of rabies virus in nature. Virology. 2002;292:247-57.

27. Carnieli Jr P, Castilho JG, Fahl WO, Véras NM, Timenetsky MCST. Genetic characterization of rabies virus isolated from cattle between 1997 and 2002 in an epizootic area in the state of São Paulo, Brazil. Virus Res. 2009;144:215-24.

28. Castilho JG, Carnieli Jr P, Oliveira RN, et al. A comparative study of rabies virus isolates from hematophagous bats in Brazil. J Wildl Dis. 2010;46:1335-9.

29. Oliveira RN, De Souza SP, Lobo RS, et al. Rabies virus in insectivorous bats: implications of the diversity of the nucleoprotein and glycoprotein genes for molecular epidemiology. Virology. 2010;405:352-60.

30. Domingo E, Holland JJ. RNA virus mutations and fitness for survival. Annu Ver Microbiol. 1997;51:151-78.

31. Kissi B, Badrane H, Lavenu A, Tordo N, Brahimi M, Bourhy H. Dynamics of virus quasispecies during serial passages in heterologous hosts. J Genl Virol. 1999;80:2041-50.

32. Wunner HW. Rabies virus. In: Jackson AC, Wunner HW, editors. Rabies. San Diego: Academic Press; 2007. p. 23-68.

33. chneider MC, Santos-Burgoa C, Aron J, Munoz B, Ruiz-Velazco S, Uieda W. Potential force of infection of human rabies transmitted by vampire bats in the Amazonian region of Brazil. Am J Trop Med Hyg. 1996;55:680-4.

34. Ito M, Arai YT, Ito UT, et al. Genetic characterization and geographic distribution of rabies virus isolates in Brazil: identification of two reservoirs, dogs and vampire bats. Virology. 2001;284:214-22.

35. Kotait I, Nogueira Filho VS, Souza MCAM, Carrieri ML, Gomes MN, Peres NF. Manual de controle da raiva dos herbívoros. 9ed. São Paulo/SP: Instituto Pasteur; 2010.

36. Reis NR, Shibatta AO, Peracchi AL, Pedro WA, Lima IP. Sobre os morcegos brasileiros. In: Reis NR, Peracchi AL, Pedro WA, Lima IP, editors. Morcegos do Brasil. Londrina: Editora da Universidade Estadual de Londrina; 2007. p. 17-24.

Received 18 May 2012

Accepted 26 July 2012

Available online 10 November 2012

Publication Dates

• Publication in this collection
  11 Dec 2012
• Date of issue
  Dec 2012

History

• Received
  18 May 2012
• Accepted
  26 July 2012