Molecular epidemiological tracing of a cattle rabies outbreak lasting less than a month in Rio Grande do Sul in southern Brazil
BMC Research Notes volume 9, Article number: 87 (2016)
Vampire bat-transmitted cattle rabies cases are typically encountered in areas where the disease is endemic. However, over the period of a month in 2009, an outbreak of cattle rabies occurred and then ended spontaneously in a small area of the Rio Grande do Sul State in southern Brazil. To investigate the epidemiological characteristics of this rabies outbreak in Rio Grande do Sul, 26 nucleotide sequences of rabies virus (RABV) genomes that were collected in this area were analyzed phylogenetically.
Nucleotide sequence identities of the nucleoprotein gene and G–L intergenic region of the 26 RABVs were greater than 99.6 %. Phylogenetic analysis showed that all RABVs clustered with the vampire bat-related cattle RABV strains and that the RABVs were mainly distributed in southern Brazil.
The findings of the present study suggested that a small population of rabid vampire bats carrying a single RABV strain produced a spatiotemporally restricted outbreak of cattle rabies in southern Brazil.
Rabies is a fatal infectious disease that causes encephalomyelitis. Genetic diversity among rabies viruses (RABVs) depends on the host species and its geographic distribution. Rabies is endemic in Brazil, and RABVs have been isolated from a variety of mammalian hosts, including dogs, foxes, cats, cattle, and hematophagous, insectivorous and frugivorous bats [1–3]. Vampire bats, particularly Desmodus rotundus, are important vectors of rabies in Latin America. Rabies transmission from vampire bats to humans occurs frequently, primarily in the Amazonian regions of Brazil and Peru, and cattle rabies transmitted by vampire bats is also common in Brazil . Vampire bat-related rabies is thus a serious problem from both animal and public health standpoints.
Molecular and geographic epidemiological analyses of livestock rabies in northeastern, central and southeastern Brazil revealed that RABVs isolated from livestock were closely related to the virus strains found in vampire bat populations, and that the RABVs in these areas were stable within a spatiotemporal context . However, relatively little information currently exists on the molecular epidemiology of RABVs isolated from livestock in southern Brazil.
The incidence of vampire bat-transmitted cattle rabies cases is typically restricted to areas where the disease is endemic [5–7]. However, between February and March 2009, an outbreak of cattle rabies occurred and ended spontaneously in a small region of Rio Grande do Sul State (RS) in southern Brazil. To investigate the epidemiological characteristics of this rabies outbreak in RS, the RABVs obtained in this area were subjected to phylogenetic analysis.
An outbreak of cattle rabies occurred in a hilly area (<5 km radius) in the city of Nova Roma do Sul in Rio Grande do Sul State, southern Brazil (Fig. 1). This outbreak occurred on a total of 19 cattle farms between February and March 2009, and ended spontaneously without any mitigation measures being implemented. Outbreaks of cattle rabies in this region had not been observed for at least 3 years before or after this outbreak. Brain specimens from these livestock were diagnosed as RABV-positive by an immunofluorescent antibody test and 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 samples were amplified by RT-PCR and determined using applied biosystems 3130 genetic analyzer (ABI) using primers described previously [8, 9]. Information of primers used for RT-PCR and sequencing are shown in Additional file 1: Table S1. Multiple alignment for the datasets of N and G–L genes were made by MAFFT program . Neighbor-joining tree of N gene with 1000 bootstrap replicates was constructed with the p-distance using MEGA ver 6.10 . The nucleotide sequences retrieved from GenBank to construct the phylogenetic tree are listed in Additional file 2: Table S2.
The geographic origins of the RABV sequences obtained from the Brazilian livestock were plotted at the municipal level in the federal states in which infection was observed using MapInfo Professional GIS software (ver. 8.0, MapInfo Japan K.K., Tokyo, Japan) and Google Earth (ver.6.0, Google, Mountain View, CA, USA). Brazilian maps were obtained from Brasil em Relevo—Embrapa Monitoramento por Satélite .
A total of 26 RNA samples were extracted from brains of rabid cattle from the 19 farms surveyed in this study (Fig. 1; Table 1). The PCR amplified sequences of the 1424-nt N-mRNA region and 578-nt G–L region were determined for all 26 RABV samples [GenBank: AB685222–AB685247 and LC057591–LC057616] (Table 1). The results from N gene sequences showed that, with the exception of BRbv1289, the nucleotide and amino acid sequence identities of the RS RABVs were 100 %. Between BRbv1285, which was selected among the 25 identical RABV RNAs as the representative sequence, and BRbv1289, nucleotide and amino acid sequence identities were more than 99.9 and 100 %, respectively. To validate the genomic homogeneity among the RABV samples, the variable G–L intergenic 578-bp region (position: 4929–5514 of PV) was additionally analyzed. The nucleotide sequence identity was more than 99.6 % among all 26 RABVs, showing the quite high homogeneity of epidemic RABVs in this outbreak.
Phylogenetic analysis based on a partial N sequence (position: 203–1420 of PV) revealed that 26 RABV samples formed a monophyletic cluster in vampire bat-related RABV lineage (Fig. 2). Comparison with RABVs obtained from neighboring countries, indicated that these 26 RABVs formed a sub-cluster with RABVs collected in the states of Sao Paulo (SP) and Rio de Janeiro (RJ) in Brazil. In addition to this sub-cluster, there are two another sub-clusters consisted of RABV samples collected in the state of SP and neighboring areas. These lineages were named the lineages of RD1, RD2/3, and Old strains by Carnieli et al. . RS strains belonged to the lineage of Old strains and shared a common ancestor with RD1 lineage. RABVs in RD1 and RD2/3 lineages were distributed in mountain range located in the boundary of the states of SP, Minas Gerais, and RJ, while Old strains were isolated from towns in the interior of the state of SP .
Using a 203-bp segment (position: 109–311 of PV) of the N gene, Kobayashi et al. reported the existence of at least 24 RABV genetic variants among the vampire bat-transmitted rabies cases in cattle in Brazil . Using the same 203-bp segment, our analysis revealed that the 26 RS RABVs of the present study correlate with the C-19 lineage described by Kobayashi et al.  (Additional file 3: Figure S1). Samples belonging to the C-19 lineage were distributed in the states of RS, Goiás, Mato Grosso do Sul, SP and RJ, which are the southeastern regions of Brazil .
Nucleoprotein gene and G–L intergenic segments from the 26 RABV samples collected in early spring of 2009 from separate farms in a restricted region of Nova Roma do Sul, RS, were sequenced and subjected to phylogenetic analysis. The genetic analysis revealed that the nucleotide and amino acid sequences were almost identical among the 26 RABVs and shared a common ancestor, suggesting the RABV strain was transmitted from a small population of rabid vampire bats infected with RABV, which may be derived from a single introduction, to a cattle population grazing in Nova Roma do Sul.
Generally, vampire bat-transmitted RABVs were restricted to specific geographic areas because of the vampire bat ecology. The phylogenetic analysis described here revealed that all of the livestock RABVs in RS were derived from vampire bat rabies. Although we were unable to investigate the prevalence of rabid vampire bats in this area, the findings of the present study revealed that the 26 RABVs correspond with the “Old strains” lineage reported by Carnieli Jr et al. and C-19 lineage identified by Kobayashi et al. [5, 13]. The distribution of RABV lineages in Brazil is known to be delimited by geomorphological boundaries, such as mountain ranges and rivers . The reason why the C-19 lineage consists of RABVs from Brazil and none from Argentina may be thus because a very large river, the Uruguay River, separates the Brazil and Argentina. However, the C-19 lineage is distributed over larger areas. Since the “Old strains” lineage, which was made up basically of the isolates collected before 1998, has been endemic in the state of SP before epizootic of RD2/RD3 were observed , suggesting that RABVs in Old strains lineages may have gradually spread in southern Brazil, including the state of RS.
The typical distribution range of D. rotundas activity around a roost is about 10 km2 ; the range of this outbreak was approximately 80 km2. In extensive random surveys conducted in Trinidad and Brazil, the prevalence of the virus among vampire bats in the field was 0.46–0.75 % (see Baer) . Further, no rabies outbreaks have been reported in the past 3 years in this region, and the outbreak ended without any specific prevention measures being implemented for cattle rabies. Although the non-migratory, vampire bats have distribution ranges of between 10 and 20 km2 , it has been suggested that the vampire bats may need to fly over longer distances to feed and that they share roosts with other bat species; indeed, this may explain how the RABV could spread over a short period of time . However, the cattle rabies cases described in the present study ended suddenly over the period of a month. There are several possible reasons why the activities of the vampire bats may have been limited and why this outbreak of cattle rabies ended suddenly. The affected area was geographically isolated by a river valley that may have limited the feeding activities and distribution of the bats. It is also possible that the size of the vampire bat population in the region may have been too small to act as reservoir for RABV.
Taken together, it is considered that this outbreak was caused by small population of vampire bats that were vectors of the RABV lineage that is distributed mainly in southern Brazil. It is considered that the relatively restricted range of bats contributed to the sudden end of the outbreak before it could spread to other regions.
Segments of N gene and G–L intergenic region sequences from 26 RABV samples that were collected in early spring from separate farms in small region of Nova Roma do Sul, Rio Grande do Sul, in 2009, were subjected to phylogenetic analysis. The N gene and G–L region sequence analysis revealed that the 26 RABVs had almost identical nucleotide and amino acid sequences. The outbreak was likely caused by a small population of vampire bats that were infected with a single strain of the RABV lineage found mainly in southern Brazil that invaded the area of the study site.
Ito M, Arai YT, Itou T, Sakai T, Ito FH, Takasaki T, Kurane I. Genetic characterization and geographic distribution of rabies virus isolates in Brazil: identification of two reservoirs, dogs and vampire bats. Virology. 2001;284:214–22.
Sato G, Kobayashi Y, Shoji Y, Sato T, Itou T, Ito FH, Santos HP, Brito CJ, Sakai T. Molecular epidemiology of rabies from Maranhão and surrounding states in the northeastern region of Brazil. Arch Virol. 2006;151:2243–51.
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:1995–2004.
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–100.
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.
Condori-Condori RE, Streicker DG, Cabezas-Sanchez C, Velasco-Villa A. Enzootic and epizootic rabies associated with vampire bats, peru. Emerg Infect Dis. 2013;19(9):1463–9.
Torres C, Lema C, Dohmen FG, Beltran F, Novaro L, Russo S, Freire MC, Velasco-Villa A, Mbayed VA, Cisterna DM. Phylodynamics of vampire bat-transmitted rabies in Argentina. Mol Ecol. 2014;23:2340–52.
Hirano S, Itou T, Carvalho AA, Ito FH, Sakai T. Epidemiology of vampire bat-transmitted rabies virus in Goiás, central Brazil: re-evaluation based on G–L intergenic region. BMC Res Notes. 2012;3:288.
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:759–66.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Brasil em Relevo—Embrapa Monitoramento por Satélite http://www.relevobr.cnpm.embrapa.br/. Accessed 13 May 2015.
Carnieli P Jr, Castilho JG, Fahl Wde O, Véras NM, Timenetsky Mdo C. 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.
Mochizuki N, Kawasaki H, Silva ML, Afonso JA, Itou T, Ito FH, Sakai T. Molecular epidemiology of livestock rabies viruses isolated in the northeastern Brazilian states of Paraíba and Pernambuco from 2003–2009. BMC Res Notes. 2012;5:32.
Arellano-Sota C. Biology, ecology, and control of the vampire bat. Rev Infect Dis. 1988;10(Suppl 4):S615–9.
Baer GM. Vampire bat and bovine paralytic rabies. In: Baer GM, editor. Natural History of Rabies. 2nd ed. Boston: CRC Press; 1991. p. 389–403.
Lord RD. Seasonal reproduction of vampire bats and its relation to seasonality of bovine rabies. J Wildl Dis. 1992;28:292–4.
Romijn PC, van der Heide R, Cattaneo CA, Silva Rde C, 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–6.
ERD and AASBC collected epidemiological samples and information in the field. TI, TF and NM drafted the manuscript. TF conducted the RT-PCR and sequencing analyses. TI, NM, YK, FHI and TS contributed to study design, management, coordination, and drafting of the manuscript. All authors read and approved the final manuscript.
This study was supported in part by a Grant-in-Aid for Scientific Research (24580453) from the Japan Society for the Promotion of Science and the Strategic Research Base Development Program, “International joint research and training of young researchers for zoonosis control in a globalized world”, and a matching fund subsidy from the Ministry of Education, Culture, Sports, Science and Technology of Japan (S0991023 and S1491007).
The authors declare that they have no competing interests.
Additional file 3. Phylogenetic tree based on nucleoprotein gene 203 nts (corresponding 107–311 nts of PV strain) using Brazilian RABV isolates. The 24 lineages classified by Kobayashi et al.  and RS RABV samples in this study were compared. All RS samples belonged to C-19 lineage maintained in Southwestern Brazil. The symbols for each lineage correspond to the results by Kobayashi et al. .
About this article
Cite this article
Itou, T., Fukayama, T., Mochizuki, N. et al. Molecular epidemiological tracing of a cattle rabies outbreak lasting less than a month in Rio Grande do Sul in southern Brazil. BMC Res Notes 9, 87 (2016). https://doi.org/10.1186/s13104-016-1898-5