Emerg Infect DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention17326958329134906-042910.3201/eid1212.060429DispatchRabies Virus Maintained by Dogs in Humans and Terrestrial Wildlife, Ceará State, BrazilRunning head: Rabies Virus, BrazilFavorettoSilvana R.*de MattosCecília C.de MoraisNélio B.§CarrieriMaria Luíza*RolimBenedito N.§SilvaLucia M.§RupprechtCharles E.DurigonEdison L.de MattosCarlos A.Instituto Pasteur de São Paulo, São Paulo, Brazil;Universidade de São Paulo, São Paulo, Brazil;Centers for Disease Control and Prevention, Atlanta, Georgia, USA;Secretaria Estadual de Saúde do Ceará, Fortaleza, BrazilAddress for correspondence: Silvana R. Favoretto, Instituto Pasteur of São Paulo, Brazil, Av: Paulista, 393, Paraíso, São Paulo, SP – Brasil, CEP: 01311000; email: srfavoretto@ig.com.br122006121219781981

Rabies viruses circulating in Ceará, Brazil, were identified by molecular analysis to be related to variants maintained by dogs, bats, and other wildlife. Most of these viruses are associated with human rabies cases. We document the emergence of a rabies virus variant responsible for an independent epidemic cycle in the crab-eating fox (Cerdocyon thous).

Keywords: RabiesBrazilzoonosisepidemiologyinfectious diseasesvirologyphylogenywildlifedomestic animals

After dog rabies control programs were implemented in Ceará State, Brazil, a the number of human cases decreased (1,2). Thereafter, the epidemiologic importance of rabies in wildlife became evident. From 1990 through 2005, a total of 173 cases of rabies were reported in Cerdocyon thous (crab-eating fox), 25 in Callithrix. j. jacchus (common marmoset) and 6 in Procyon cancrivorous (crab-eating raccoon). During this period, in 13 of 40 human cases reported in Ceará, wildlife was the source of infection (2).

In 1996, because of this new epidemiologic situation, public health authorities launched an educational program, and no human cases due to wildlife were recorded in 1999, despite 84 cases in wildlife registered that year (2). The objective of this study was to elucidate some of the epidemiologic events involved in rabies emergence among wildlife in Ceará.

The Study

We studied 22 samples, from dogs, cattle, wildlife, and humans in Ceará, obtained from 1997 to 2003 (Table). Samples were antigenically characterized by using a monoclonal antibody (MAb) panel against the viral nucleoprotein (35). Isolates were injected into the brains of suckling mice, and brain impressions were made for MAb typing (3,5). Characterization of the samples identified 4 antigenic variants. Antigenic variant-2 (AgV2), maintained by dogs, was found in all C. thous, P. crancrivourous, and human cases and in all dog isolates with the exception of brdg5360, which was positive with all the MAbs. Antigenic variant-3 (AgV3), epidemiologically associated with vampire bats, Desmodus rotundus, was identified in 3 bovine samples. A previously reported profile, representing an AgV that circulates in marmosets in Ceará (5), was detected in sample brsg5696.

Identification, origin, and antigenic and genetic variant of 22 rabies virus samples isolated from Ceará State, Brazil
IdentificationAnimal speciesYear of isolationOriginAntigenic variantGenetic variantGroupGenBank accession no.
Brhm4531Human1997FortalezaAgV2DogBDQ447947
Brcth4122Cerdocyon thous1998-AgV2DogBDQ447948
Brhm5325Human2000CaucaiaAgV2DogCDQ447949
Brdg5360Dog2000CaucaiaAll+DogCDQ447950
Brcth5361C. thous2000ParacuruAgV2DogBDQ447951
Brbv5339Bovine2000Antonina do NorteAgV3Vampire batDDQ447952
Brbv5374Bovine2000QuixereAgV3Vampire batDDQ447953
Brhm5691Human2001CaucaiaAgV2DogCDQ447954
Brdg5693Dog2001MaranguapeAgV2DogADQ447955
Brcth5695C. thous2001BarroquinhaAgV2DogBDQ447956
Brcth5697C. thous2001CanindeAgV2DogBDQ447957
Brcth5692C. thous2001MaranguapeAgV2DogBDQ447958
Brpcr5698P. cancrivorous2001MaranguapeAgV2DogBDQ447959
Brbv5694Bovine2001AquirazAgV3Vampire batDDQ447960
Brsg5696C. j. jacchus2001CaucaiaAgV new*MarmosetEDQ447961
Brhmu138Human2002FortalezaAgV2DogCDQ447962
Brhmu142Human2003FortalezaAgV2DogCDQ447963
Brhmu130Human2003UmirimAgV2DogCDQ447964
Brhmu146Human2003FortalezaAgV2DogCDQ447965
Brhmu145Human2003TururuAgV2DogCDQ447966
Brhmu129Human2003FortalezaAgV2DogCDQ447967
Brhmu131Human2003MaracanauAgV2DogCDQ447968

*Pattern related to isolates from marmosets.

The Ceará viruses were analyzed genetically through a comparative phylogenetic study based on a 320-bp fragment of the nucleoprotein gene, from position 1157 to 1476, as compared with SADB19 (58). These isolates were also compared with rabies virus variants circulating among domestic animals and wildlife from the Americas. The viral RNA was extracted from infected tissues, and the cDNA was obtained by reverse transcription–PCR techniques, using primers 21 g and 304, and was sequenced with primer 304 (7,9). The phylogenetic analyses were made by using the PileUp program of the Wisconsin Package Version 10.1 (10) and the programs DNADIST, NEIGHBOR, SEQBOOT, and CONSENSE of the PHYLIP package (11). The expressed percentages of identity refer to the nucleotide sequences. The trees were obtained with the TREEVIEW program (12). The phylogenetic analyses showed segregation in 5 lineages, A–E (Figure 1), which was statistically supported by high bootstrap values.

Rabies virus isolates by geographic localization and neighbor-joining tree showing a comparison of the groups formed by Ceará State, Brazil, samples isolated from 1997 to 2003. Bootstrap values of >50% obtained from 100 resamplings of the data using distance matrix methods are shown in the nods.

Lineage A was represented by a sample from a dog from Maranguape, which was obtained in 2001. This virus showed its closest genetic relationship with lineage B (identity 92.4%–94.2%). Lineage B was formed by all the C. thous isolates, a sample from a human bitten by a P. cancrivorous raccoon in Fortaleza in 1997, and a virus from a P. crancrivorous raccoon collected in Maranguape during 2001 (intrinsic identity 96.5%–100%). This lineage showed its highest percentage of identity with lineage C (intrinsic identity 90.6%–92.8%). Lineage C consisted of 9 human samples collected in 5 different counties from 2000 to 2003 and an isolate obtained from a dog in 2000. The samples were highly homologous (intrinsic identity 99.1%–100%). Lineage D included 3 bovines collected in 3 geographically distant counties during 2000 and 2001 (intrinsic identity 97.5%–98.4%). Lineage E was represented by the only sample collected from a C. j jacchus marmoset, These last 2 lineages were related distantly to all the others.

When compared with representatives of rabies variants maintained by terrestrial and bat species in the Americas (Figure 2), lineages A, B, and C continue to segregate as independent lineages with high statistical support. Ceará bovine samples representing lineage D clustered with D. rotundus and D. rotundus–related cases from Latin America (intrinsic identity 94%–97.8%).

Neighbor-joining tree showing a comparison of Ceará samples (groups A, B, C, D, E) with isolates obtained from the Americas. Bootstrap values of >50% obtained from 100 resamplings of the data using distance matrix methods are shown in the nods. The sequences from Latin America used in the comparison were identified as as follows: group I, dogs and terrestrial wildlife from Mexico, Venezuela, Colombia, Dominican Republic, and Peru (mxsk, skunk from Mexico; mxdg and mxmx, dog from Mexico; vedg, dog from Venezuela; codg, dog from Colombia; drmg, mongoose from Dominican Republic; mxgm, bobcat from Mexico; pefx, fox from Peru) and terrestrial wildlife from the United States (caussk, skunk from California; txuscy, coyote from Texas; wiussk, striped skunk from Wisconsin; arussk, striped skunk from Arkansas; azusfox, gray fox from Arizona; txusfx, gray fox from Texas); and group II, terrestrial wildlife from the United States (flrac, raccoon from Florida; parac, raccoon from Philadelphia; ksussk, striped skunk from Kansas; arussk, striped skunk from Arkansas). The antigenic variant and endemic cycle to which it belongs are shown in the tree. (GenBank accession no. AB201803 is a vampire bat from Brazil and nos. AY654585, AY654587, AY654586 are humans and a marmoset from Brazil). EC, endemic cycle.

The only sample representative of lineage E segregated with 2 isolates from humans bitten by C. j. jacchus and a sample collected from a marmoset kept as a pet (5). The isolates were highly homologous to each other (intrinsic identity 98.7%–100%).

Conclusions

A thorough description of rabies epidemiology depends on a comprehensive surveillance program and application of accurate molecular methods to discriminate among different variants and the emergence of new foci. Antigenic and limited sequencing analyses were used to better understand the emergent epidemiologic events in wildlife in Ceará, Brazil. These analyses allowed identification of 5 potential cycles in this region, despite antigenic homogeneity.

Lack of antigenic and genetic relationships of sample brdg5693, representing lineage A, with the rest of the isolates from Ceará and the known terrestrial rabies vectors from the Americas shows that this virus is a variant not previously described. This virus was geographically and temporarily associated with samples brpcr5698 and brcth5692, obtained in Maranguape during 2001. These circumstances demonstrated the existence of at least 2 overlapping endemic cycles in this area. Lineage B was formed mainly by isolates from C. thous, which indicates the existence of an emerging rabies cycle in this species.

The epidemiologic situation in Ceará was complicated because of overlapping distributions of dog and C. thous rabies cases (Figure 1). Tree topology and genetic relationships between dog and C. thous variants suggested that the canine virus was introduced in C. thous populations because of spillover events, which gave rise to an emergent cycle. A similar event was described between domestic dogs and Canis adustus (jackal) in Zimbawe. In this case, the variant circulating in dogs was introduced into the C. adustus population by spillover events, with the consequent emergence of an independent cycle (13). Recently, the hoary fox has been identified as a rabies reservoir in Brazil (14).

Inclusion in lineage B of an isolate obtained from a human bitten by a P. cancrivorous raccoon and another sample collected from this species suggested the risk of establishing C. thous variant in P. cancrivorous. The niches of these 2 species overlap, which facilitates their encounters. Additional surveillance is necessary to clarify this situation.

Epidemiologic data which indicates that humans had been exposed to dog bites, results of molecular characterization, and inclusion of a dog isolate in the C lineage strongly incriminate the dog as the reservoir of this variant. Identification of the source of infection by using classic surveillance alone is complicated by the presence of multiple cycles of transmission. Genetic comparison of samples from lineage D with viruses representing bats viruses from the Americas helped to identify D. rotundus as the source of livestock infection.

The close genetic relationship of sample brsg5696 with rabies isolates obtained from C. j. jacchus and human cases bitten by marmosets further supported C. j. jaccuss as the most important vector of this variant. This finding indicates that this species plays an important role for disease maintenance in nature.

Methods for antigenic and genetic identification of rabies samples isolated in the Americas have contributed effectively to the development of health programs, as well as recognition of possible wild reservoirs of urban rabies. The emergence of new cycles in Latin American wildlife indicates the need to strengthen surveillance programs in these species and research development for the evaluation of the feasibility of oral vaccination interventions.

Suggested citation for this article: Favoretto SR, de Mattos CC, de Morais NB, Carrieri ML, Rolim BN, Silva LM, et al. Rabies virus maintained by dogs in humans and terrestrial wildlife, Ceará State, Brazil. Emerg Infect Dis [serial on the Internet]. 2006 Dec. [date cited]. http://dx.doi.org/10.3201/eid1212.060429

Acknowledgments

We thank Débora R. V. Sacramento for her technical collaboration, Fabiana L. A. da Silveira Lebrun for structural support, Maria das Graças Silva for help with the bibliography and illustrations, and all staff at Instituto Pasteur of São Paulo, Brazil, and Laboratory of Rabies Diagnosis of Ceará, Brazil, for collaboration and incalculable incentives.

Dr Favoretto is a scientific researcher and specialist in public health at the Instituto Pasteur of São Paulo. Her research interests include diagnosis and molecular epidemiology of rabies virus.

ReferencesPan-American Health Organization/World Health Organization Prevention and control of diseases. Veterinary Public Health Unit. Elimination of the human rabies transmitted by dogs in Latin America: analysis of the situation, 2004. Washington. Pan American Health Organization. 2005;37:1029Núcleo de Controle das Endemias Transmissíveis por vetores–NUEND, Coordenadoria de Apoio ao Desenvolvimento da Atenção à Saúde–CODAS, 1998–2003. Boletim da Secretaria da Saúde do Estado do Ceará, Brasil, 2004Favoretto SR, Carrieri ML, Cunha EMS, Aguiar EAC, Silva LHQ, Sodré M, Antigenic typing of Brazilian rabies virus samples isolated from animals and humans, 1989–2000. Rev Inst Med Trop Sao Paulo. 2002;44:915 10.1590/S0036-4665200200020000712048546Pan-American Health Organization/World Health Organization Los anticuerpos monoclonales en la caracterización y vigilancia de los virus de la rábia en América Latina y el Caribe. Pan American Journal of Public Health. 2000;8:2147 10.1590/S1020-49892000000800016Favoretto SR, Mattos CC, Morais NB, Araújo FAA, Mattos CA Rabies in marmosets (Callithrix jacchus) from the state of Ceará, Brazil. Emerg Infect Dis. 2001;7:10625 10.3201/eid0706.01063011747745Conzelmann KK, Cox JK, Schneider LG, Thiel HJ Molecular cloning and complete sequence of the attenuated rabies virus SABB19. Virology. 1990;175:48599 10.1016/0042-6822(90)90433-R2139267Smith JS Rabies virus. In.Murray PR, Baron EJ,Pfaller MA, Tenover FC, Yolken R, editors. Manual of clinical microbiology, 6th ed. Washington: American Society for Microbiology Press; 1995 p. 997–1003.de Mattos CC, de Mattos CA, Loza-Rubio E, Aguilar-Setién A, Orciari LA, Smith JS Molecular characterization of rabies virus isolates from Mexico: implications for transmission dynamics and human risk. Am J Trop Med Hyg. 1999;61:5879710548293de Mattos CA, de Mattos CC, Smith JS, Miller E, Papo S, Utrera A, Genetic characterization of rabies field isolates from Venezuela. J Clin Microbiol. 1996;34:155388735118Version WP 10.1 Madison (WI): Genetics Computer Group; 2000Felsenstein J Phylip: Phylogeny inference pakage. Version 3,5c. Seattle (WA): University of Washington; 1993Page R D. TREEVIEW: An application to display phylogenetic trees on personal computers. Computer Applied Biocience. 1996;12:3578Bingham J, Foggim CM, Wandeler AI, Hill FW The epidemiology of rabies in Zimbabwe. 2. Rabies in jackals (Canis adustus and Canis mesomelas). Onderstepoort J Vet Res. 1999;1:112310396757Bernardi F, Nadin Davis SA, Wandeler AI, Armstrong J, Gomes A, Lima FS, Antigenic and genetic characterization of rabies viruses isolated from domestic and wild animals of Brazil identifies the hoary fox as a rabies reservoir. J Gen Virol. 2005;86:315362 10.1099/vir.0.81223-016227239