- Short Report
- Open access
- Published:
Epidemiology of vampire bat-transmitted rabies virus in Goiás, central Brazil: re-evaluation based on G-L intergenic region
BMC Research Notes volume 3, Article number: 288 (2010)
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.
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.
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].
References
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.
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. 10.1590/S0037-86822002000500006.
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.
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. 10.1006/tpbi.2000.1494.
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-10.1186/1746-6148-4-44.
Arellano-Sota C: Biology, ecology, and control of the vampire bat. Rev Infect Dis. 1988, 10: S615-619.
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. 10.1016/j.virusres.2004.12.003.
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. 10.1016/j.virusres.2009.02.012.
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. 10.1292/jvms.68.1097.
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. 10.2307/2388777.
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. 10.1292/jvms.66.747.
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. 10.1093/molbev/msm092.
Hills DM, Bull JJ: An empirical test of boostrapping as a method for assessing confidence in phylogenetic analysis. Syst Biol. 1993, 42: 182-192.
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. 10.1093/nar/gki352.
Pond SL, Frost SD, Muse SV: HyPhy: hypothesis testing using phylogenies. Bioinformatics. 2005, 21: 676-679. 10.1093/bioinformatics/bti079.
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.
Instituto Brasileiro de Geografia e Estatística (IBGE). [http://www.ibge.gov.br/servidor_arquivos_geo]
DIVA-GIS gData. [http://www.diva-gis.org/gData]
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.
Author information
Authors and Affiliations
Corresponding author
Additional information
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.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to 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.
About this article
Cite this article
Hirano, S., Itou, T., Carvalho, A.A. et al. Epidemiology of vampire bat-transmitted rabies virus in Goiás, central Brazil: re-evaluation based on G-L intergenic region. BMC Res Notes 3, 288 (2010). https://doi.org/10.1186/1756-0500-3-288
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1756-0500-3-288