Limited or no epidemiological information has been reported for rabies viruses (RABVs) isolated from livestock in the northeastern Brazilian states of Paraíba (PB) and Pernambuco (PE). The aim of this study was to clarify the molecular epidemiology of RABVs circulating in livestock, especially cattle, in these areas between 2003 and 2009.
Findings
Phylogenetic analysis based on 890 nt of the nucleoprotein (N) gene revealed that the 52 livestock-derived RABV isolates characterized here belonged to a single lineage. These isolates clustered with a vampire bat-related RABV lineage previously identified in other states in Brazil; within PB and PE, this lineage was divided between the previously characterized main lineage and a novel sub-lineage.
Conclusions
The occurrences of livestock rabies in PB and PE originated from vampire bat RABVs, and the causative RABV lineage has been circulating in this area of northeastern Brazil for at least 7 years. This distribution pattern may correlate to that of a vampire bat population isolated by geographic barriers.
Background
Rabies is a fatal infectious disease that causes encephalomyelitis. In Brazil, various rabies viruses (RABVs) have been isolated from numerous animal species, including dogs, foxes, cats, and cattle, as well as from hematophagous, insectivorous, and frugivorous bats. Vampire bats, particularly Desmodus rotundus, are an important rabies vector in Latin America. Transmission from vampire bats to humans has been reported, primarily in the Amazon regions of Brazil and Peru, and a large number of cases of cattle rabies transmitted by vampire bats also have been reported in Brazil [1]. Since the introduction of a regional elimination program, the incidence of human and canine rabies in Latin America has fallen by 90% over the past 20 years. However, northeastern Brazil remains a "hotspot" for human rabies because of circulation of the virus among the dog population [2]. Carnieli Jr. et al. reported the molecular characterization and epidemiology of RABVs isolated from canids in northeastern Brazil [3–5], and Shoji et al. reported the genetic and phylogenetic characterization of RABVs isolated from wild fox, insectivorous bats, and livestock in Paraíba (PB) [6]. In Olinda, a city in Pernambuco (PE), 7,062 patients underwent prophylactic antirabies treatment between 2002 and 2006 [7]. Molecular and geographic analyses of livestock rabies in central and southeast Brazil revealed that RABVs isolated from livestock were related to the virus found in vampire bat populations, and this epidemiological pattern was maintained over time and space in these areas [8–11]. However, little or no information on the molecular epidemiology of RABVs has been reported for the virus isolated from livestock in PB and PE. The aim of the present study was to analyse the molecular epidemiology of RABVs circulating among livestock, especially cattle, in these areas between 2003 and 2009.
Results
The sequences of 890 nt PCR products, corresponding to nucleotides 89-978 of the Pasteur vaccine (PV) strain, were determined for all 52 RABV isolates. Among the 52 isolates, the nucleotide and amino acid sequence identities were 97.7-100% and 97.9-100%, respectively.
Phylogenetic analysis based on the sequences of 890 nt of the N gene revealed that the 52 RABV isolates included in this study clustered with a vampire bat-related RABV lineage; these isolates did not cluster with the dog-, fox-, or insectivorous bat-related RABVs (Figure 1). Comparison with RABVs isolated from other states in Brazil indicated that these 52 RABV isolates belonged to a single lineage; furthermore, this lineage was divided between a previously characterized main lineage (B) and a novel sub-lineage (A) consisting of several isolates located in PE. Geographical plotting showed that RABV isolates of the novel sub-lineage were derived from neighbouring areas (Figure 2). The topographical distributions of isolates from lineages A and B could not be distinguished in the areas covered by this study.
Figure 1
A phylogenetic tree based on the nucleotide sequences of the N gene. A phylogenetic tree based on the nucleotide sequences of 890 nt (bases 89-978) of the N gene was constructed using the method devised by Saitou and Nei [12]; the bootstrap probabilities of each node were calculated using 1,000 replicates. The designations BRbv, BRsp, BRgt, BRhr, BR-Pfx, and BR-DR indicate samples from Brazilian cattle, sheep, goats, horses, foxes, and vampire bats, respectively. MOKV and ABLV denote the Mokola virus and Australian bat lyssavirus, respectively. State abbreviations are as follows: PB, Paraíba; PE, Pernambuco; GO, Goiás; SP, São Paulo; RJ, Rio de Janeiro; MT, Mato Grosso; TO, Tocantins; MA, Maranhão; PA, Pará; MS, Mato Grosso do Sul. The triangle and star symbols represent the new sub-lineage (group A) and the previously reported lineage (group B), respectively.
Figure 2
Geographic distribution of rabies virus isolates in the Brazilian states of Paraíba and Pernambuco. (a) Map of Brazil indicating location of states pertinent to this study. State abbreviations are as follows: PB, Paraíba; PE Pernambuco; GO, Goiás; SP, São Paulo; RJ, Rio de Janeiro; MT, Mato Grosso; TO, Tocantins; MA, Maranhão; PA, Pará; MS, Mato Grosso do Sul. (b, c) Detailed geographic distribution of livestock isolates classified as genetic variants in the states of PB (b) and PE (c). The two symbols (triangle and star) correspond to the new sub-lineage and previously reported lineage, respectively, as used in Figure 1. Samples for which the geographic origin and the genetic variant are identical are illustrated using the same symbol. Brazilian maps were obtained from Brasil em Relevo - Embrapa Monitoramento por Satélite http://www.relevobr.cnpm.embrapa.br/.
Discussion
Segments of the N gene from the 52 RABV isolates, collected from PB and PE between 2003 and 2009, were sequenced and phylogenetically analysed. This N gene segment displayed greater than 97.7% nucleotide and amino acid sequence identity among these 52 RABVs. This correlation implies that the lineage has been maintained during transmission in PB and PE.
The phylogenetic analysis described here indicates that all of the livestock RABVs in PB and PE were derived from vampire bat rabies. Furthermore, the present study reveals that there are two RABV lineages that are separate from other regional Brazilian vampire bat-related RABV lineages. The first, sub-lineage A, consists of BRbv1173 (collected in PB in 2004) and another seven RABV isolates collected in PE in 2008-2009; the second, sub-lineage B, consists of isolates collected in PB in 2003-2009 and in PE in 2007-2009. Comparison among RABV isolates from multiple northern states of Brazil revealed that lineages from PB and PE were distinct from those obtained from the Maranhão (MA) state. Because of geographical barriers (mountains and rivers) between PB/PE and MA, it would be difficult for vampire bats (the prime vector for RABV) to move freely between PB/PE and MA. Geographical mapping (Figure 2) demonstrates that sub-lineage A was located mainly in high-altitude areas of PE, while sub-lineage B was widely distributed and present in both PE and PB. Thus, these two lineages seem to correlate with geographic factors and/or vampire bat populations, but do not appear to correlate with the year of isolation. Using a 203 nt segment (bases 109 to 311) of the N gene, Kobayashi et al. [8, 9] reported the existence of at least 24 RABV genetic variants among vampire bat-transmitted cases of rabies in cattle in Brazil; the distribution of several of these RABV genetic variants was found to be delimited by geographic boundaries, including mountain ranges and rivers. Our analysis (data not shown) using the same 203 nt segment reveals that the 52 RABV isolates (from PB and PE) of the present study correlate with the same PB lineage described by Kobayashi et al. These results indicate that livestock rabies has been transmitted by vampire bats in PB and PE during this study period. Furthermore, this RABV lineage seems to have been circulating in this area for at least 7 years, with transmission affected by geographic factors and resulting in dispersion of vampire bats among regional populations.
Previous reports have shown human exposure to vampire bat-transmitted rabies in northeastern Brazil, including PE [1, 7]. The present phylogenetic analysis suggests that rabies epidemics that occurred in cattle in PB and PE were transmitted by vampire bats. These RABV isolates comprised a lineage independent from that of other Brazilian isolates, with the distinction reflecting isolation from neighbouring regions by geographic factors. Thus, the vampire bat-derived rabies in this area represents an endemic disease, suggesting that the regional control of vampire bat rabies in this area may be a workable model for local elimination of human and livestock rabies.
Conclusions
The present study indicates that occurrences of livestock rabies in PB and PE were caused by vampire bat RABVs, and that this RABV lineage has been circulating in this area of northeastern Brazil for at least 7 years. This pattern of distribution may correlate to that of a vampire bat population isolated by geographic barriers.
Methods
The 52 RABV isolates used in this study were obtained from cattle (46), sheep (3), goats (2), and horse (1), and were collected in PB and PE between 2003 and 2009 (Table 1). Brain specimens from these livestock were diagnosed as RABV-positive by an immunofluorescent antibody test and a mouse inoculation test. These study procedures were implemented in accordance with the Institutional Guidelines for Animal Experiments at the Campina Grande University under the permission (number 129/2009) of the Committee for Experimental Animals of this College. Viral RNA was extracted from the brains of livestock using the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany). Nucleoprotein (N) gene sequences from the Brazilian RABV isolates were amplified using RT-PCR with primers JW12 (5'-ATGTAACACCYCTACAATG-3') (position: 55-73 of PV) and N8 (5'-AGTTTCTTCAGCCATCTC-3') (position: 1585-1568 of PV), followed by hemi-nested PCR with primer pairs as follows. Primer pair A consisted of JW12 and P2 (5'-CCCATATAACATCCAACAAAGTG-3') (position: 1029-1007 of PV), and generated a 975-nt amplicon. Primer pair B consisted of P1 (5'-CTACAATGGATGCCGACAAGA-3') (position: 66-86 of PV) and N8, and generated a 1,520 nt amplicon. The nucleotide sequences of RABVs from Brazilian foxes, livestock, and hematophagous and insectivorous bats were obtained from GenBank (Table 1). Cycle sequencing, nucleotide and amino acid sequence alignments, and phylogenetic analyses were performed as previously described [13, 14]. The geographic origins of the RABV isolates sequenced from the Brazilian livestock were plotted at the municipal level of the respective federal states using MapInfo Professional GIS software (ver. 8.0, MapInfo Japan K.K., Tokyo, Japan). Brazilian maps were obtained from Brasil em Relevo - Embrapa Monitoramento por Satélite [15].
Table 1
Brazilian rabies virus isolates used in this study
*1 Names in parentheses are designations from Shoji et al. (2006)
*2 State abbreviations are as follows: PB Paraíba, PE Pernambuco, GO Goiás, SP São Paulo, RJ Rio de Janeiro, MT Mato Grosso, TO Tocantins, MA Maranhão, PA Pará, MS Mato Grosso do Sul
*3 Lineages are based on the phylogenetic tree of Figure 1
*4 We determined extra sequences in this study
Declarations
Acknowledgements
This work was supported in part by a Grant-in-Aid for the Academic Frontier Project for Private Universities from the Ministry of Education, Culture, Sports, Science, and Technology of Japan.
Authors’ Affiliations
(1)
Nihon University Veterinary Research Center, Nihon University
(2)
Academic Unit of Veterinary Medicine, Center of Health and Rural Technology, Federal University of Campina Grande
(3)
Bovine Clinic, Federal Rural University of Pernambuco, Campus of Garanhuns
(4)
Department of Preventive Veterinary Medicine and Animal health, Faculty of Veterinary Medicine and Zootechny, University of São Paulo
References
Schneider MC, Romijn PC, Uieda W, Tamayo H, da Silva DF, Belotto A, da Silva JB, Leanes LF: Rabies transmitted by vampire bats to humans: an emerging zoonotic disease in Latin America?. Rev Panam Salud Publica. 2009, 25 (3): 260-269.PubMedView ArticleGoogle Scholar
Schneider MC, Aguilera XP, Barbosa da Silva J, Ault SK, Najera P, Martinez J, Requejo R, Nicholls RS, Yadon Z, Silva JC, et al: Elimination of neglected diseases in latin america and the Caribbean: a mapping of selected diseases. PLoS Negl Trop Dis. 2011, 5 (2): e964-10.1371/journal.pntd.0000964.PubMedPubMed CentralView ArticleGoogle Scholar
Carnieli P, Brandão PE, Carrieri ML, Castilho JG, Macedo CI, Machado LM, Rangel N, de Carvalho RC, de Carvalho VA, Montebello L, et al: Molecular epidemiology of rabies virus strains isolated from wild canids in Northeastern Brazil. Virus Res. 2006, 120 (1-2): 113-120. 10.1016/j.virusres.2006.02.007.PubMedView ArticleGoogle Scholar
Carnieli P, Castilho JG, Fahl Wde O, Veras NM, Carrieri ML, Kotait I: Molecular characterization of rabies virus isolates from dogs and crab-eating foxes in Northeastern Brazil. Virus Res. 2009, 141 (1): 81-89. 10.1016/j.virusres.2008.12.015.PubMedView ArticleGoogle Scholar
Carnieli P, Fahl Wde O, Castilho JG, Oliveira Rde N, Macedo CI, Durymanova E, Jorge RS, Morato RG, Spindola RO, Machado LM, et al: Characterization of rabies virus isolated from canids and identification of the main wild canid host in Northeastern Brazil. Virus Res. 2008, 131 (1): 33-46. 10.1016/j.virusres.2007.08.007.PubMedView ArticleGoogle Scholar
Shoji Y, Kobayashi Y, Sato G, Gomes AA, Itou T, Ito FH, Sakai T: Genetic and phylogenetic characterization of rabies virus isolates from wildlife and livestock in Paraiba, Brazil. Acta Virol. 2006, 50 (1): 33-37.PubMedGoogle Scholar
Dantas-Torres F, Oliveira-Filho EF: Human exposure to potential rabies virus transmitters in Olinda, State of Pernambuco, between 2002 and 2006. Rev Soc Bras Med Trop. 2007, 40 (6): 617-621.PubMedView ArticleGoogle Scholar
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.PubMedView ArticleGoogle Scholar
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.PubMedPubMed CentralView ArticleGoogle Scholar
Carnieli P, Castilho JG, Fahl Wde O, Veras NM, Timenetsky Mdo C: Genetic characterization of Rabies virus isolated from cattle between 1997 and 2002 in an epizootic area in the state of Sao Paulo, Brazil. Virus Res. 2009, 144 (1-2): 215-224. 10.1016/j.virusres.2009.05.001.PubMedView ArticleGoogle Scholar
Macedo CI, Carnieli P, Fahl Wde O, Lima JY, Oliveira Rde N, Achkar SM, Castilho JG, Carrieri ML, Kotait I: Genetic characterization of rabies virus isolated from bovines and equines between 2007 and 2008, in the States of Sao Paulo and Minas Gerais. Rev Soc Bras Med Trop. 2010, 43 (2): 116-120. 10.1590/S0037-86822010000200002.PubMedView ArticleGoogle Scholar
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMedGoogle Scholar
Mochizuki N, Kobayashi Y, Sato G, Hirano S, Itou T, Ito FH, Sakai T: Determination and Molecular Analysis of the Complete Genome Sequence of Two Wild-Type Rabies Viruses Isolated from a Haematophagous Bat and a Frugivorous Bat in Brazil. J Vet Med Sci. 2011, 73 (6): 759-766. 10.1292/jvms.10-0238.PubMedView ArticleGoogle Scholar
Mochizuki N, Kobayashi Y, Sato G, Itou T, Gomes AA, Ito FH, Sakai T: Complete genome analysis of a rabies virus isolate from Brazilian wild fox. Arch Virol. 2009, 154 (9): 1475-1488. 10.1007/s00705-009-0475-9.PubMedView ArticleGoogle Scholar
Kobayashi Y, Sato G, Kato M, Itou T, Cunha EM, Silva MV, Mota CS, Ito FH, Sakai T: Genetic diversity of bat rabies viruses in Brazil. Arch Virol. 2007, 152 (11): 1995-2004. 10.1007/s00705-007-1033-y.PubMedView ArticleGoogle Scholar
Ito M, Itou T, Shoji Y, Sakai T, Ito FH, Arai YT, Takasaki T, Kurane I: Discrimination between dog-related and vampire bat-related rabies viruses in Brazil by strain-specific reverse transcriptase-polymerase chain reaction and restriction fragment length polymorphism analysis. J Clin Virol. 2003, 26 (3): 317-330. 10.1016/S1386-6532(02)00048-3.PubMedView ArticleGoogle Scholar
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.
