Printer Friendly

Novel sylvatic rabies virus variant in endangered golden palm civet, Sri Lanka.

Rabies has been eliminated from domestic animals in industrialized countries, but sylvatic rabies remains an endemic disease. The ecology of rabies in wildlife populations and natural ecosystems is poorly understood (1), and, as a result, eliminating rabies from the wild is difficult. Little is known about sylvatic rabies in developing countries, where rabies takes its biggest toll on humans. Rabies is endemic to Sri Lanka and has been identified in different wild animals. However, all documented cases of rabies in wildlife in Sri Lanka have been considered a consequence of spillover from dogs. Rabies viruses circulating in this country are distinctly highly homogeneous (2,3).

Two species of palm civet are commonly found in Sri Lanka: the common palm civet, Paradoxurus hermaphroditus, which is widespread in southern Asia and Southeast Asia, and the golden palm civet, P. zeylonensis, which is indigenous to Sri Lanka. This species is closely related to the brown palm civet (P. jerdoni), which lives only in southern India (4). Moreover, 3 additional new species have been identified in Sri Lanka: the golden wetzone palm civet (P. aureus), the golden dry-zone palm civet (P. stenocephalus), and the Sri Lankan brown palm civet (P. montanus) (4). Palm civets in Sri Lanka are, however, endangered because of hunting, parasitic diseases, and dwindling habitat. We report on a sylvatic rabies virus variant detected in a golden palm civet in Sri Lanka.

The Study

On a November morning in 2009, a "wild cat" appeared in the garden of a basic health clinic in Moneragala district, Uva Province, Sri Lanka. The animal, which showed aggressive behavior, was suspected to be rabid and was thus killed to prevent transmission of rabies virus to humans. The animal's head was packed in ice to avoid decomposition and sent to the Medical Research Institute (Colombo, Sri Lanka) for testing. We detected rabies virus in the animal's brain by using the fluorescent antibody test and extracted viral RNA and DNA by using Trizol (Invitrogen, Carlsbad, CA, USA). The rabies virus from this sample was designated as H-1413-09. The whole genome of the virus was sequenced directly from the sample as described (5).

To confirm the species of the rabid animal, we determined the nucleotide sequence of the mitochondrial cytochrome b (cytb) gene and performed a BLAST search ( for similarity with other sequences. By aligning nucleotide sequences of the cytb gene of mitochondrial DNA of domestic cat, jungle cat, fishing cat, Asiatic golden cat, marbled pole cat, European pole cat, lynx, puma, leopard, African lion, tiger, jaguar, civet, and palm civet with ClustalW2 (, we designed primer Felis cytb-F, 5'-ATGACCAACATTCGAAAATCACACC-3' (nt 1-25), and primer Felis cytb-R, 5'-CAATAAT GCCTGAGATGGGTATTAG-3' (nt 1093-1,117). Using these primers, we performed PCR as follows: initial denaturation at 94[degrees]C for 2 min, followed by 94[degrees]C for 30 s, 55[degrees]C for 30 s, and 72[degrees]C for 1 min for 35 cycles, followed by a final extension at 72[degrees]C for 5 min. PCR generated a 1,117-bp fragment from which a 1,004-nt sequence was determined. Analysis showed that the sequence has 100% identity with the partial (224-nt) sequence of the cytb gene of P. zeylonensis (GenBank accession no. FJ881681); this is the only sequence available for P. zeylonensis. The sequence also has 95% identity with P. jerdoni and 90%-92% identity with P. hermaphroditus.

We performed an evolutionary analysis by using the N gene. We inferred a maximum clade credibility phylogenetic tree by using the Bayesian Markov chain Monte Carlo method available in BEAST version 1.6.1 (6). The analysis used a relaxed (uncorrelated lognormal) molecular clock and GTR + [GAMMA] + I model of nucleotide substitution. We selected the model on the basis of Akaike Information Criterion by using jModelTest software (7). All chains were run for 9 x [10.sup.7] generations and sampled every 3,000 steps. This procedure resulted in an effective sample size of >2,000 for all estimated parameters. The posterior densities were calculated with 10% burn-in and checked for convergence by using Tracer version 1.5 ( The mean rate of nucleotide substitution estimated for the N gene was 2.2 x [10.sup.4] substitutions/site/year (95% highest posterior density [HPD] values 1.3-3.2 x [10.sup.4] substitutions/site/ year). This rate is in agreement with previous findings (8). Approximately 155.5 years ago (95% HPD 91.3-249.5 years)--that is, circa 1854 (95% HPD range 1760-1918)--rabies viruses from Sri Lanka and southern India diverged from their most recent common ancestor. Approximately 76.4 years ago (95% HPD 46.4-123.0 years)--that is, in [approximately equal to] 1933 (95% HPD range 1886-1963)--strain H-1413-09 diverged from canine rabies virus in Sri Lanka.

Compared with the genome sequence of rabies virus strain H-08-1320 from Sri Lanka, the genome sequence of strain H-1413-09 had a nucleotide deletion at residue 2,417 and an addition at nt 11,807. As a result, the start and stop signals of the mRNA and the start and stop positions of coding sequences of different genes were altered (Table 1). These altered start and stop codons advanced the coordinates of the intergenic signal of relevant intergenic regions. The P-M region and terminal sequence of strain H-08-1320 were 1 nt shorter (91 nt) and longer (136 nt), respectively. The nucleotide and amino acid identities between the coding regions of strains H-1413-09 and H-081320 are shown in Table 1. The substitutions detected in the deduced amino acid sequences of H-1413-09 were compared with the genomic sequence of H-08-1320 (Table 2). By using the MEGA5 Tamura-Nei model (www., we determined the genetic distance of the N gene to determine whether strain H-1413-09 is more diverse than other rabies viruses in Sri Lanka (Figure). The rate of variation among sites was modeled with a gamma distribution (shape parameter = 0.5). The genetic distances between strain H-1413-09 and other rabies viruses from Sri Lanka (0.027-0.036) were greater than the genetic distances among other rabies viruses (0.001-0.011). These results support our finding that strain H-1413-09 differs from other rabies viruses circulating in Sri Lanka.



Rabies virus probably survives favorably in the wild because it can infect a large spectrum of animals, thereby maximizing replication and dispersal opportunities (9). Most viruses replicate poorly when transferred to new hosts, but greater genetic variation assists in such species adaptation (10). Increased mutation in an RNA virus like rabies virus can give rise to variants with altered levels of fitness to persist and spread. A large number of substitutions were found in strain H-1413-09 compared with strain H-081320; these substitutions might represent changes that resulted from species adaptation. Phylogenetic analysis and comparative sequence data indicated that strain H-1413-09 is a variant rabies virus.

Palm civets are facing extinction in Sri Lanka because the species is losing its habitat, being hunted for its meat, and dying of parasitic diseases (www.sundaytimes. lk/090118/Plus/sundaytimesplus_01.html). Our study indicates that rabies might be another risk factor for extinction of these animals. Identification of a variant rabies virus in wildlife has serious implications for rabies control in Sri Lanka. Identification of such a virus would help provide epidemiologic data about the spread of rabies and its incursion into new geographic regions and would justify allocation of increased resources to help control rabies (11,12).

Several rabies virus variants associated with wildlife are known in the Americas and Africa (1,13-15), and this report identified classical sylvatic rabies in Asia. Whether P. zeylonensis is a reservoir of rabies virus or represents spillover from another animal deserves extensive investigation. The detection of rabies in wildlife indicates that much remains to be discovered in the tropical ecosystem of Sri Lanka. The circulation of a sylvatic variant rabies virus may be another hurdle in the rabies-control effort in Sri Lanka.

This study was supported in part by a Grant-in-Aid for Scientific Research (B) from the Japan Society for the Promotion of Sciences (grant no. 20406026) and by the Research Fund at the Discretion of the President, Oita University (grant no. 610000N5010).

Dr Matsumoto is a project research associate in the Department of Microbiology, Faculty of Medicine, Oita University. His research interest is molecular epidemiology of rabies virus.

DOI: 10.3201/eid1712.110811


(1.) East ML, Hofer H, Cox JH, Wulle U, Wiik H, Pitra C. Regular exposure to rabies virus and lack of symptomatic disease in Seren geti spotted hyenas. Proc Natl Acad Sci U S A. 2001;98:15026-31. doi:10.1073/pnas.261411898

(2.) Arai YT, Takahashi H, Kameoka Y, Shiino T, Wimalaratne O, Lodmell DL. Characterization of Sri Lanka rabies virus isolates using nucleotide sequence analysis of nucleoprotein gene. Acta Virol. 2001;45:327-33.

(3.) Nanayakkara S, Smith JS, Rupprecht CE. Rabies in Sri Lanka: splendid isolation. Emerg Infect Dis. 2003;9:368-71.

(4.) Groves CP, Rajapaksha C, Manemandra-Arachchi K. The taxonomy of the endemic golden palm civet of Sri Lanka. Zool J Linn Soc. 2009;155:238-51. doi:10.1111/j.1096-3642.2008.00451.x

(5.) Matsumoto T, Ahmed K, Wimalaratne O, Yamada K, Nanayakkara S, Perera D, et al. Whole-genome analysis of a human rabies virus from Sri Lanka. Arch Virol. 2011;156:659-69. doi:10.1007/s00705 010-0905-8

(6.) Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214. doi:10.1186/14712148-7-214

(7.) Posada D, Crandall KA. MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998;14:817-8. doi:10.1093/ bioinformatics/14.9.817

(8.) Bourhy H, Reynes JM, Dunham EJ, Dacheux L, Larrous F, Huong VT, et al. The origin and phylogeography of dog rabies virus. J Gen Virol. 2008;89:2673-81. doi:10.1099/vir.0.2008/003913-0

(9.) Keesing F, Holt RD, Ostfeld RS. Effects of species diversity on disease risk. Ecol Lett. 2006;9:485-98. doi:10.1111/j.14610248.2006.00885.x

(10.) Cutler SJ, Fooks AR, van der Poel WH. Public health threat of new, reemerging, and neglected zoonoses in the industrialized world. Emerg Infect Dis. 2010;16:1-7. doi:10.3201/eid1601.081467

(11.) Nadin-Davis SA, Turner G, Paul JP, Madhusudana SN, Wandeler AI. Emergence of Arctic-like rabies lineage in India. Emerg Infect Dis. 2007;13:111-6. doi:10.3201/eid1301.060702

(12.) Slate D, Rupprecht CE, Rooney JA, Donovan D, Lein DH, Chip man RB. Status of oral rabies vaccination in wild carnivores in the United States. Virus Res. 2005;111:68-76. doi:10.1016/j. virusres.2005.03.012

(13.) De Mattos CC, De Mattos CA, Loza-Rubio E, Aguilar-Setien 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:587-97.

(14.) Carnieli P Jr, 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:81-9. doi:10.1016/j.virusres.2008.12.015

(15.) Van Zyl N, Markotter W, Nel LH. Evolutionary history of African mongoose rabies. Virus Res. 2010;150:93-102. doi:10.1016/j. virusres.2010.02.018

Takashi Matsumoto, Kamruddin Ahmed, Omala Wimalaratne, Susilakanthi Nanayakkara, Devika Perera, Dushantha Karunanayake, and Akira Nishizono

Author affiliation: Oita University, Yufu, Japan (T. Matsumoto, K. Ahmed, A. Nishizono); and Medical Research Institute, Colombo, Sri Lanka (O. Wimalaratne, S. Nanayakkara, D. Perera, D. Karunanayake)

Address for correspondence: Kamruddin Ahmed, Research Promotion Project, Oita University, Yufu 879-5593, Oita, Japan; email: ahmed@
Table 1. Percentage identity shared between
genes of 2 rabies virus strains from Sri Lanka,
H-08-1320 and H-1413-09, by gene coding region *

Gene H-1413-09 H-08-1320

 Coding region Start codon, Coding region Start codon,
 stop codon stop codon

N 71-1,423 ATG, TGA 71-1,423 ATG, TGA
P 1,514-2,404 ATG, TGA 1,514-2,404 ATG, TGA
M 2,496-3,104 ATG, TAA 2,497-3,105 ATG, TAA
G 3,317-4,891 ATG, TGA 3,318-4,892 ATG, TGA
L 5,407-11,793 ATG, TGA 5,408-11,794 ATG, TAA

Gene % Identity

 Nucleotide Amino acid

N 99.1 97.1
P 99.1 97.5
M 98.0 96.9
G 98.7 97.4
L 99.0 97.0

* H-08-1320 is a human strain typical of canine
rabies virus circulating in Sri Lanka; H-1413-09
is a novel sylvatic rabies virus variant from a
golden palm civet in Sri Lanka.

Table 2. Substitutions in genome sequence of
rabies virus strain H-1413-09 from Sri Lanka,
compared with genome sequence of strain
H-08-1320 *

Protein, amino acid substitution Site/domain/region
 of protein ([dagger])

 [Leu.sub.80] [right arrow] [Phe.sub.80]
 [Glu.sub.110] [right arrow] [Asp.sub.110]
 [Ile.sub.246] [right arrow] [Val.sub.246]
 [Ala.sub.372] [right arrow] [Val.sub.372] Antigenic site I
 [Gln.sub.167] [right arrow] [Arg.sub.167] N protein binding
 site in
 [Ile.sub.16] [right arrow] [Ala.sub.16]
 [Pro.sub.19] [right arrow] [Ser.sub.19]
 [Ile.sub.55] [right arrow] [Val.sub.55]
 [Lys.sub.77] [right arrow] [Arg.sub.77]
 [Val.sub.193] [right arrow] [Ile.sub.193]
 [Arg.sub.264] [right arrow] [His.sub.264]
 [Ile.sub.449] [right arrow] [Thr.sub.449] Transmembrane region
 [Thr.sub.459] [right arrow] [Ile.sub.459] Transmembrane region
 [Ala.sub.467] [right arrow] [Thr.sub.467] Cytoplasmic domain
 [Glu.sub.475] [right arrow] [Gly.sub.457] Cytoplasmic domain
 [Asn.sub.499] [right arrow] [Ser.sub.499] Cytoplasmic domain
 [Ser.sub.26] [right arrow] [Pro.sub.26]
 [Ile.sub.49] [right arrow] [Leu.sub.49]
 [Cys.sub.137] [right arrow] [Tyr.sub.137]
 [Leu.sub.222] [right arrow] [Ile.sub.222]
 [Ser.sub.312] [right arrow] [Gln.sub.312] Conserved domain I
 [Glu.sub.313] [right arrow] [Lys.sub.313] Conserved domain I
 [Ser.sub.314] [right arrow] [Ala.sub.314] Conserved domain I
 [Arg.sub.315] [right arrow] [Glu.sub.315] Conserved domain I
 [Val.sub.317] [right arrow] [Phe.sub.317] Conserved domain I
 [Lys.sub.1056] [right arrow] [Arg.sub.1056] Conserved domain IV
 [Thr.sub.1137] [right arrow] [Val.sub.1137] Conserved domain V
 [Ala.sub.1520] [right arrow] [Glu.sub.1520]
 [Ile.sub.1555] [right arrow] [Val.sub.1555]
 [Leu.sub.1570] [right arrow] [Met.sub.1570]
 [Met.sub.1577] [right arrow] [Leu.sub.1577]
 [Lys.sub.1625] [right arrow] [Arg.sub.1625]
 [Asn.sub.1763] [right arrow] [Asp.sub.1763]
 [Arg.sub.1876] [right arrow] [His.sub.1876]
 [Asn.sub.2023] [right arrow] [Asp.sub.1763]
 [Gly.sub.2098] [right arrow] [Arg.sub.2098]
 [Leu.sub.2113 [right arrow] [Phe.sub.2113]

* H-08-1320 is a human strain typical of canine
rabies virus circulating in Sri Lanka; H-1413-09
is a novel sylvatic rabies virus variant from a
golden palm civet in Sri Lanka.

([dagger]) Blank spaces indicate no site/domain/
region has been identified in that portion of
the protein.
COPYRIGHT 2011 U.S. National Center for Infectious Diseases
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2011 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:DISPATCHES
Author:Matsumoto, Takashi; Ahmed, Kamruddin; Wimalaratne, Omala; Nanayakkara, Susilakanthi; Perera, Devika;
Publication:Emerging Infectious Diseases
Geographic Code:9SRIL
Date:Dec 1, 2011
Previous Article:Characterization of African swine fever virus Caucasus isolate in European wild boars.
Next Article:Rickettsia parkeri in Amblyomma maculatum ticks, North Carolina, USA, 2009-2010.

Terms of use | Copyright © 2018 Farlex, Inc. | Feedback | For webmasters