Email updates

Keep up to date with the latest news and content from BMC Research Notes and BioMed Central.

Open Access Short Report

Epidemiology of vampire bat-transmitted rabies virus in Goiás, central Brazil: re-evaluation based on G-L intergenic region

Shinji Hirano1, Takuya Itou1*, Adolorata AB Carvalho2, Fumio H Ito3 and Takeo Sakai1

Author Affiliations

1 Nihon University Veterinary Research Center, 1866 Kameino, Fujisawa 252-0880, Kanagawa, Japan

2 Faculty of Agriculture and Veterinary Science, Department of Preventive Veterinary Medicine, UNESP, Via de Acesso Prof. Paulo Donato Castellane, Jaboticabal, São Paulo 14884-900, Brazil

3 Department of Preventive Veterinary Medicine and Animal Health, Faculty of Veterinary Medicine and Zootechny, University of São Paulo, Av. Prof. Dr. Orlando Marques de Paiva, 87, Cidade Universtiátria, São Paulo 05508-000, Brazil

For all author emails, please log on.

BMC Research Notes 2010, 3:288  doi:10.1186/1756-0500-3-288

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1756-0500/3/288


Received:26 May 2010
Accepted:8 November 2010
Published:8 November 2010

© 2010 Itou et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Vampire bat related rabies harms both livestock industry and public health sector in central Brazil. The geographical distributions of vampire bat-transmitted rabies virus variants are delimited by mountain chains. These findings were elucidated by analyzing a high conserved nucleoprotein gene. This study aims to elucidate the detailed epidemiological characters of vampire bat-transmitted rabies virus by phylogenetic methods based on 619-nt sequence including unconserved G-L intergenic region.

Findings

The vampire bat-transmitted rabies virus isolates divided into 8 phylogenetic lineages in the previous nucleoprotein gene analysis were divided into 10 phylogenetic lineages with significant bootstrap values. The distributions of most variants were reconfirmed to be delimited by mountain chains. Furthermore, variants in undulating areas have narrow distributions and are apparently separated by mountain ridges.

Conclusions

This study demonstrates that the 619-nt sequence including G-L intergenic region is more useful for a state-level phylogenetic analysis of rabies virus than the partial nucleoprotein gene, and simultaneously that the distribution of vampire bat-transmitted RABV variants tends to be separated not only by mountain chains but also by mountain ridges, thus suggesting that the diversity of vampire bat-transmitted RABV variants was delimited by geographical undulations.

Background

Rabies is a zoonosis that kills infected mammals, including humans, and is mainly transmitted by carnivores. In the Americas, chiropterans (insectivorous, frugivorous and hematophagous bat) are another reservoir of this disease. Although dog-transmitted rabies in central Brazil has been reduced by aggressive vaccination programs [1], chiroptera (particularly the common vampire bat, Desmodus rotundas)-transmitted rabies remains endemic in this region, and harms both the livestock industry and the public health sector [2,3].

To date, vampire bat-transmitted rabies in livestock has been controlled by reducing the population of vampire bats and by vaccinating livestock [3,4]. However, the depopulation of vampire bats has limitations and the effects are temporary, while vaccination of livestock is only carried out for some animals and is ineffective in decreasing rabies levels in vampire bats.

For the sustainable and effective control of vampire bat rabies, further knowledge of epidemiological features, such as vampire bat ecology and the dynamics of vampire bat-transmitted rabies, is necessary. Molecular epidemiological analysis of vampire bat-transmitted cattle rabies cases using the partial nucleoprotein gene, which is the most conserved gene in the rabies virus (RABV) genome, has suggested that the distribution of variants in Brazil is delimited by mountain chains and clustered in tens of thousands of square kilometers [5]. However, the vampire bats migrate several kilometers from their nests [6]. To elucidate a more detailed genetic clustering and geo-distribution of genetic clades of vampire bat-transmitted RABV, 204 isolates from Goiás, which includes the 185 isolates analyzed previously, were employed and analyzed by a phylogenetic method based on a nucleotide sequence encompassing the G-L intergenic region locating between glycoprotein (G) and polymerase (L) gene loci, which is the most divergent region in the RABV genome and is used for monitoring epidemiological changes in the evolution of RABV [7,8].

Results

Phylogenetic analysis

The 204 RABV isolates divided into 8 phylogenetic lineages in the previous nucleoprotein gene analysis were divided into 10 phylogenetic lineages with significant bootstrap values (Figure 1; details shown in Table 1). Isolates of the C-5 and C-6 lineages designated by Kobayashi et al. belonged to the A-lineage, while the B-lineage consisted of some isolates of C-5. The C-, D-, F- and J-lineages included isolates belonging to C-12, C-1, C-22 and C-3, respectively. Isolates of C-20 belonged to the I-lineage, and C-21 was divided into three lineages; G-, H- and I-lineages. Two isolates not belonging to any lineages in previous studies were assigned to the E-lineage.

thumbnailFigure 1. NJ tree based on 619-nt sequence of the partial G and the G-L intergenic region. Full circles on internal branches indicate significant bootstrap values (>70%). The number of isolates retaining 100% nucleotide similarity with the representative isolate is shown in parentheses. Lineages of the isolates defined in a previous study [5] are shown in italics after a diagonal (Ud: undefined).

Table 1. Isolates from Goiás

Geographical plotting

In this study, geographical areas in Goiás were divided by mountain chains into the Northwest, North central, Northeast and South regions (Figure 2). Most isolates of the A-lineage were distributed in the South region. The B-lineage was likely to exist in the Northeast and North central regions. The isolates belonging to the C-, E-, F-, G-, H- and I-lineages were distributed in the Northwest region. The G-, H- and I-lineages have narrow distributions and are apparently separated by mountain ridges (Figure 2; Area I), while the A-lineage was distributed in a wide range throughout a Southeastern basin (Figure 2; Area II). The isolates of the D-lineage were plotted on an eastern edge of the North central region. The J-lineage had a wide geographical distribution in Goiás.

thumbnailFigure 2. Distribution of vampire bat-transmitted RABV isolates in Goiás. Area I: distribution of G-, H- and I-lineages in mountain ridges. Area II: distribution of A-lineage in the basin. Symbols are same as those used in Fig. 1. Numbers inside the symbols indicate the number of isolates obtained in the same region. Small and open symbols indicated that only one isolate was obtained.

Discussion

Previous studies had elucidated that the distributions of vampire bat and transmitted RABV are delimited by mountain chains [5]. In the present study, it was reconfirmed that mountain chains divide the distribution patterns of each viral lineage. Furthermore, the isolates belonging to C-5 having a wide range in Goiás were divided into the A- and B-lineages, and were found to be distributed in the South and North regions on either side of a mountain chain. This finding supports Kobayashi's hypothesis that distribution of vampire bat-transmitted RABV is affected by mountain chains.

The same variants of vampire bat-transmitted RABV were spread widely in flat low lands (< 800 m), but at higher elevations (800-1600 m), they had a narrower distribution [9]. However, the G-, H- and I-lineages were found to be separated by mountain ridges in low areas (400-800 m) located in a southern undulating area of the Northwestern region (Figure 2; Area I). Furthermore, the A-lineage was located in an eastern basin of the South region (> 800 m; Figure 2; Area II). Considering that the higher lands showed an undulating landscape, the results suggest that the distribution patterns of vampire bat-transmitted RABV variants depend on such undulations. On the other hand, the distribution of common vampire bats in a valley is limited by the ridges that form the valley [10], thus supporting the notion that the distribution of RABV variants is affected by smaller topography than mountain chains.

Conclusions

The present study analyzed the epidemiology of vampire bat-transmitted RABV using a 619-nt region containing the partial glycoprotein gene and the G-L intergenic region, and indicated that the isolates can be further divided into several phylogenetic lineages with significant bootstrap values when compared to characterization based on the 203-nt partial N gene. Furthermore, the phylogenetic lineages were divided by both mountain chains and mountain ridges. In future studies, it will be important to analyze samples from different time points and to elucidate the dynamics of vampire bat-transmitted rabies in order to establish effective and sustainable control measures for preventing rabies circulation among vampire bats.

Methods

Samples

A total of 204 samples obtained from 192 cattle, 10 horses and 2 sheep in Goiás from October 2001 to August 2002 were employed, which had been confirmed as rabies positive through fluorescent antibody test and mouse inoculation test (Table 1). Viral RNA was extracted from the brain as described previously [11]. Lineages of the 164 cattle isolates, C-1, C-3, C-5, C-6, C-12, C-20, C-21 and C-22, were previously characterized based on a 203-nt sequence of the nucleoprotein gene [5], and are shown in Table 1.

Determination of nucleotide sequences

RT-PCR and direct sequencing with the HmG5-1302 (4619TGTTGAAGTTCACCTYCCMGATGT4642, positions relative to PV strain genome (Accession No. M13215) and RVLa-1 (5435ATRGGGTCATCATAAACCTC5416) primer pair were performed as described previously [11]. The target sequence includes the partial glycoprotein gene and G-L intergenic region. Nucleotide sequences were determined using the ATGC program version 4.02 (GENETYX Co., Tokyo, Japan).

Phylogenetioc analysis

Multiple nucleotide sequence alignments of the partial glycoprotein gene and G-L intergenic region were generated by the ClustalW package in MEGA ver. 4.0 [12]. A phylogenetic tree was constructed by the neighbor joining (NJ) method with bootstrap analysis (1000 psuedoreplicates) under the p-distance model. Phylogenetic clustering supported by a bootstrap value exceeding 70% was regarded as a reliable lineage [13]. Results were validated by the maximum likelihood method using PhyML [14]. In order to reconfirm the shape of the NJ tree, the ML tree was constructed under HKY substitution model justified by MODELTEST packaged in Hyphy program [15].

Geographical plotting

The 204 RABV isolates were plotted onto a geographical map described using the DIVA-GIS program [16] with GIS data from Instituto Brasileiro de Geografia e Estatística [17] and DIVA-GIS gData [18].

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

SH participated in the design of the study, performed the experimental procedures and the data analysis, and wrote the manuscript. TI, AABC, FHI and TS elaborated the study design, management, coordination, and drafting the manuscript. The authors have read and approved the final manuscript.

Acknowledgements

This work was partly supported by the "Strategic Research Base Development" Program for Private Universities subsidized by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan.

References

  1. Romijn PC, van der Heide R, Cattaneo CA, Silva RC, van der Poel WH: Study of lyssaviruses of bat origin as a source of rabies for other animal species in the State of Rio De Janeiro, Brazil.

    Am J Trop Med Hyg 2003, 69:81-86. PubMed Abstract | Publisher Full Text OpenURL

  2. Goncalves MA, Sa Neto RJ, Brazil TK: Outbreak of aggressions and transmission of rabies in human beings by vampire bats in northeastern Brazil.

    Rev Soc Bras Med Trop 2002, 35:461-464. PubMed Abstract | Publisher Full Text OpenURL

  3. Mayen F: Haematophagous bats in Brazil, their role in rabies transmission, impact on public health, livestock industry and alternatives to an indiscriminate reduction of bat population.

    J Vet Med B Infect Dis Vet Public Health 2003, 50:469-472. PubMed Abstract | Publisher Full Text OpenURL

  4. Massad E, Coutinho FA, Burattini MN, Sallum PC, Lopez LF: A mixed ectoparasite--microparasite model for bat-transmitted rabies.

    Theor Popul Biol 2001, 60:265-279. PubMed Abstract | Publisher Full Text OpenURL

  5. Kobayashi Y, Sato G, Mochizuki N, Hirano S, Itou T, Carvalho AA, Albas A, Santos HP, Ito FH, Sakai T: Molecular and geographic analyses of vampire bat-transmitted cattle rabies in central Brazil.

    BMC Vet Res 2008, 4:44. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  6. Arellano-Sota C: Biology, ecology, and control of the vampire bat.

    Rev Infect Dis 1988, 10:S615-619. PubMed Abstract OpenURL

  7. Nel LH, Sabeta CT, von Teichman B, Jaftha JB, Rupprecht CE, Bingham J: Mongoose rabies in southern Africa: a re-evaluation based on molecular epidemiology.

    Virus Res 2005, 109:165-173. PubMed Abstract | Publisher Full Text OpenURL

  8. Ngoepe CE, Sabeta C, Nel L: The spread of canine rabies into Free State province of South Africa: A molecular epidemiological characterization.

    Virus Res 2009, 142:175-180. PubMed Abstract | Publisher Full Text OpenURL

  9. Kobayashi Y, Ogawa A, Sato G, Sato T, Itou T, Samara SI, Carvalho AA, Nociti DP, Ito FH, Sakai T: Geographical distribution of vampire bat-related cattle rabies in Brazil.

    J Vet Med Sci 2006, 68:1097-1100. PubMed Abstract | Publisher Full Text OpenURL

  10. Trajano E: Movements of Cave Bats in Southeastern Brazil, with Emphasis on the Population Ecology of the Common Vampire Bat, Desmodus rotundus (Chiroptera).

    Biotropica 1996, 28:121-129. Publisher Full Text OpenURL

  11. Sato G, Itou T, Shoji Y, Miura Y, Mikami T, Ito M, Kurane I, Samara SI, Carvalho AA, Nociti DP, Ito FH, Sakai T: Genetic and phylogenetic analysis of glycoprotein of rabies virus isolated from several species in Brazil.

    J Vet Med Sci 2004, 66:747-753. PubMed Abstract | Publisher Full Text OpenURL

  12. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0.

    Mol Biol Evol 2007, 24:1596-1599. PubMed Abstract | Publisher Full Text OpenURL

  13. Hills DM, Bull JJ: An empirical test of boostrapping as a method for assessing confidence in phylogenetic analysis.

    Syst Biol 1993, 42:182-192. OpenURL

  14. Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML Online--a web server for fast maximum likelihood-based phylogenetic inference.

    Nucleic Acids Res 2005, 33:W557-559. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Pond SL, Frost SD, Muse SV: HyPhy: hypothesis testing using phylogenies.

    Bioinformatics 2005, 21:676-679. PubMed Abstract | Publisher Full Text OpenURL

  16. Hijmans RJ, Guarino L, Cruz M, Rojas E: Computer tools for spatial analysis of plant genetic resources data: 1. DIVA-GIS.

    Plant Genet Res Newsl 2002, 127:15-19. OpenURL

  17. Instituto Brasileiro de Geografia e Estatística (IBGE) [http://www.ibge.gov.br/servidor_arquivos_geo] webcite

  18. DIVA-GIS gData [http://www.diva-gis.org/gData] webcite