Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention24274823384089112-150610.3201/eid1912.121506DispatchDispatchDistinct Lineage of Vesiculovirus from Big Brown Bats, United StatesVesiculovirus from Big Brown BatsNgTerry Fei FanDriscollCindyCarlosMaria PazPrioleauAlgernonSchmiederRobertDwivediBhaktiWongJakkChaYunheeHeadStevenBreitbartMyaDelwartEricBlood Systems Research Institute, San Francisco, California, USA (T.F.F. Ng, J. Wong, Y. Cha, E. Delwart); University of California, San Francisco (T.F.F. Ng, E. Delwart); University of South Florida, St. Petersburg, Florida, USA (T.F.F. Ng, B. Dwivedi, M. Breitbart); Maryland Department of Natural Resources, Oxford, Maryland, USA (C. Driscoll); Maryland Department of Health and Mental Hygiene, Baltimore, Maryland, USA (M.P. Carlos, A. Prioleau); San Diego State University, San Diego, California, USA (R. Schmieder); The Scripps Research Institute, La Jolla, California, USA (S. Head)Address for correspondence: Terry Fei Fan Ng, Blood Systems Research Institute, Virology, 270 Masonic Ave, San Francisco, CA 94116, USA; email: terryfeifan@gmail.com122013191219781980

We identified a novel rhabdovirus, American bat vesiculovirus, from postmortem tissue samples from 120 rabies-negative big brown bats with a history of human contact. Five percent of the tested bats were infected with this virus. The extent of zoonotic exposure and possible health effects in humans from this virus are unknown.

Keywords: Rhabdoviridaevesiculovirusviral metagenomicbatsvirus discoverybig brown batsvirusesNorth AmericaMarylandUnited StatesrabiesEptesicus fuscus

Bats are reservoirs for many emerging viral pathogens, including Ebola viruses, Marburg viruses, henipaviruses, and severe acute respiratory syndrome coronaviruses; >80 bat virus species have been characterized (1,2). The diversity of these viruses and their high infection rates in bats may be attributed to multiple factors might that facilitate virus transmission and maintenance, including bats’ large social group size, high species diversity, long life, long-distance migration, roost sharing by multiple species, and social habits such as mutual grooming and biting (1,2).

Rabies virus (family Rhabdoviridae, genus Lyssavirus) is commonly detected in bats from the United States. Analyses of several cases of human rabies infections have reported insectivorous bats as the source (3). The Rhabdoviridae family contains 6 formally approved genera, but most bat rhabdoviruses belong to the Lyssavirus genus (Figure). Nonrabies lyssaviruses have been characterized from bats in other parts of the world, including Australia, Europe, Africa, and Asia (49). In contrast to the known diversity in bats of the extensively analyzed Lyssavirus genus, the diversity of other Rhabdoviridae genera in bats remains largely undetermined. Vesiculoviruses (genus Vesiculovirus), such as vesicular stomatitis virus, cause fever and vesicular diseases in animals such as cattle, horses, and pigs. Some vesiculoviruses, including Chandipura virus and vesicular stomatitis virus, are also zoonotic and cause acute diseases in humans.

Analyses of American bat vesiculovirus (ABVV) compared with other members of the family Rhabdoviridae. A) Genome organization of ABVV; B) Bayesian inference tree of the ABVV N gene; C) Bayesian inference tree of the 5 concatenated ABVV genes (N, P, M, G, L). For the Bayesian analyses, sequences from the entire gene were used, except for a few partially sequenced genomes for which only ≈100 aa were publicly available. Posterior probabilities (>75%) of the Bayesian analysis are shown next to each node. Formally classified vesiculoviruses are labeled with “V,” whereas potential vesiculoviruses not formally recognized by the International Committee on Taxonomy of Viruses are labeled with “V?.” Distinct clades of bat rhabdoviruses are labeled with triangles of different colors: black, vesiculovirus; gray, lyssavirus; white, unclassified. N, nucleoprotein; P, phosphoprotein; M, matrix protein; G, glycoprotein; L, polymerase protein. Scale bar indicates amino acid substitutions per site.

The bat virome has not been fully characterized. Most bat virome studies have been conducted by analyzing fecal, anal swab, or pharyngeal swab specimens from healthy bats (1013). These studies have revealed a variety of viruses but no new rhabdoviruses. However, viruses in fecal and pharyngeal samples could include ingested and inhaled viruses that originated from insects and plants (1013). To focus specifically on viruses infecting the bats themselves, we performed unbiased metagenomic sequencing of RNA viruses purified from the lungs and livers of 120 rabies-negative big brown bats (Eptesicus fuscus) collected in Maryland, USA.

The Study

During 2008, more than 500 bats associated with possible human exposure were submitted to the Maryland Department of Health and Mental Hygiene State Laboratory for postmortem diagnosis of rabies by direct fluorescent antibody assay (Technical Appendix). For this study, virus particles were purified from the lungs and livers of 120 rabies-negative bats with good carcass condition, and viral nucleic acids were extracted, randomly amplified, and sequenced by using 454 pyrosequencing (Roche, Mannheim, Germany) and Solexa Illumina sequencing (Illumina, San Diego, CA, USA). A total of ≈100,000 pyrosequences and 13.5 million Solexa Illumina sequences were generated and then assembled to form contigs. More than 30 contigs showed low protein identities to known vesiculoviruses (BLASTx; http://blast.ncbi.nlm.nih.gov/Blast.cgi), indicating a possible novel virus. PCR and rapid amplification of cDNA ends were performed to obtain the complete genome of this virus (primers shown in Technical Appendix Figures 1, 2). We proposed the name American bat vesiculovirus (ABVV) for this virus.

The full-length ABVV genome (GenBank accession no. JX569193) consists of 10,692 nt of negative-sense, single-stranded RNA beginning with a 103-nt 5′ untranslated region, followed by open reading frames encoding for the nucleoprotein, phosphoprotein, matrix protein, glycoprotein, and polymerase protein genes (Figure, panel A). Phylogenetic analyses showed that ABVV is related to vesiculoviruses including Chandipura virus and Isfahan virus, both of which are associated with encephalitic illness in humans. ABVV is located close to the root of vesiculoviruses in the Bayesian analysis of the nucleoprotein gene (Figure, panel B) and shares 41%–49% aa identity with known vesiculoviruses, similar to the vesiculovirus interspecies identities reported (47.9%–72.5%) and higher than the intergenera identities between vesiculoviruses, lyssaviruses, and ephemeroviruses (17.0%–33.1%) (14). Analyses of the polymerase gene alone (Technical Appendix Figure 3) and of a concatenation of all 5 genes (Figure, panel C) suggested that ABVV lies within the vesiculovirus clade, more closely related to mammalian than fish vesiculoviruses. Combined, these analyses indicate that ABVV is likely to belong to the Vesiculovirus genus, rather than representing a novel genus. The basal phylogenetic position of ABVV suggests early divergence from other mammalian vesiculovirus species.

Lung and liver tissues from 60 of the bats used for the pooled metagenomic analyses were screened individually for ABVV by using reverse transcription PCR targeting the polymerase gene (Technical Appendix). Three (5%) bats tested positive for ABVV: 1 adult female, 1 adult male, and 1 juvenile male. Viral RNA was found in liver tissue from the 2 male bats and in lung and liver tissues from the female adult bat.

Considering the extensive lyssavirus diversity in bats, we hypothesize that bat vesiculoviruses are similarly diverse. To facilitate characterization of diverse vesiculoviruses in bats, we designed 2 pairs of degenerate PCR primers (VesiConAF-KCDGAYAARAGYCAYTCVATGA; VesiConAR-TGNGCNACDGTNARDGCATT; VesiConBF-GGNMGRTTYTTYTCHYTDATGTC; VesiConBR-TCHGCNGAYTGCATNGTYTCA) on the basis of a sequence alignment of the polymerase gene of ABVV and the formally classified mammalian vesiculoviruses. When the ABVV-positive bat liver cDNA was used as a control, the nested PCR yielded an amplicon of 704 bases, and its sequence was confirmed by cloning and Sanger sequencing. Future studies may use these pan-vesiculovirus PCR primers to investigate vesiculovirus diversity in other bat species and in other regions.

Conclusions

Big brown bats are prevalent in North America, where their geographic range overlaps extensively with that of humans, and considerable interactions occur between big brown bats and humans and their pets. Big brown bats from this region are a known reservoir of rabies virus; our analysis shows that these bats also constitute a sylvatic mammalian reservoir of vesiculoviruses.

The characterization of ABVV sheds light on vesiculovirus diversity in bats. The other bat rhabdoviruses—Mount Elgon bat virus, Oita rhabdovirus, and Kern Canyon virus—do not belong to the vesiculovirus clade but cluster together in a separate clade (14) (Figure, panel B). A recent report described several rhabdoviruses in oropharyngeal swab specimens from Spanish bats, but the short reads (100 bases) precluded a detailed phylogenetic analysis (15). Of the bats tested in our study, 5% were infected with ABVV, a finding that suggests vesiculoviruses are prevalent in bats. The characterization of a novel rhabdovirus in bats with a history of human contact raises questions for further research, including health effects on the virus’ hosts, seroprevalence, possible transmission by insect vectors, and the extent of zoonotic exposure in humans. ABVV-specific and vesiculovirus-consensus PCRs, as well as future endeavors to culture this virus, will help address these questions.

<supplementary-material content-type="local-data" id="SD1"><caption><title>Technical Appendix

Detailed methods for postmortem diagnosis of rabies by direct fluorescent antibody assay in >500 bats associated with possible human exposure submitted to the Maryland Department of Health and Mental Hygiene State Laboratory; location of the primers used in this study relative to the American bat vesiculovirus; sequences of the primers used in this study; and Bayesian inference tree based on the polymerase gene depicting relationships among the members of the family Rhabdoviridae.

Suggested citation for this article: Ng TFF, Driscoll C, Carlos MP, Prioleau A, Schmieder R, Dwivedi B, et al. Distinct lineage of vesiculovirus from big brown bats, United States. Emerg Infect Dis [Internet]. 2013 Dec [date cited]. http://dx.doi.org/10.3201/eid1912.121506

Acknowledgments

We thank Robert A. Myers, Katherine Feldman, Dana Limpert, and the laboratory scientists at the Rabies Laboratory of the Maryland Department of Health and Mental Hygiene for sample collection and discussion. We also thank John Shimashita, Lana Schaffer, and the staff at the Next Generation Sequencing Core at the Scripps Research Institute for providing DNA sequencing support.

T.F.F.N., J.W., Y.C., and E.L.D. were supported by the Blood Systems Research Institute and the National Heart, Lung, and Blood Institute (R01HL083254). M.B. was funded by National Science Foundation Biodiversity Inventories grant DEB-1025915.

Dr Ng is a postdoctoral fellow at the Blood System Research Institute and University of California, San Francisco. His current research focuses on metagenomic discovery of novel human and animal viruses using deep sequencing, including rhabdoviruses, caliciviruses, picornaviruses, picaliviruses, hepeviruses, and astroviruses.

ReferencesWong S, Lau S, Woo P, Yuen KY. Bats as a continuing source of emerging infections in humans. Rev Med Virol. 2007;17:6791 and. 10.1002/rmv.52017042030Calisher CH, Childs JE, Field HE, Holmes KV, Schountz T. Bats: important reservoir hosts of emerging viruses. Clin Microbiol Rev. 2006;19:53145 and. 10.1128/CMR.00017-0616847084Lyles D, Rupprecht C. Rhabdoviridae. In: Knipe D, Howley P, editors. Fields virology, 5th ed. Philadelphia: Williams & Wilkins; 2006Kuzmin IV, Hughes GJ, Botvinkin AD, Orciari LA, Rupprecht CE. Phylogenetic relationships of Irkut and West Caucasian bat viruses within the Lyssavirus genus and suggested quantitative criteria based on the N gene sequence for lyssavirus genotype definition. Virus Res. 2005;111:2843 and. 10.1016/j.virusres.2005.03.00815896400Bourhy H, Kissi B, Tordo N. Molecular diversity of the Lyssavirus genus. Virology. 1993;194:7081 and. 10.1006/viro.1993.12368386891Arai YT, Kuzmin IV, Kameoka Y, Botvinkin AD. New lyssavirus genotype from the lesser mouse-eared bat (Myotis blythi), Kyrghyzstan. Emerg Infect Dis. 2003;9:3337 and. 10.3201/eid0903.02025212643828Botvinkin AD, Poleschuk EM, Kuzmin IV, Borisova TI, Gazaryan SV, Yager P, Novel lyssaviruses isolated from bats in Russia. Emerg Infect Dis. 2003;9:16235 and. 10.3201/eid0912.03037414720408Kuzmin IV, Orciari LA, Arai YT, Smith JS, Hanlon CA, Kameoka Y, Bat lyssaviruses (Aravan and Khujand) from Central Asia: phylogenetic relationships according to N, P and G gene sequences. Virus Res. 2003;97:6579 and. 10.1016/S0168-1702(03)00217-X14602198Kuzmin IV, Mayer AE, Niezgoda M, Markotter W, Agwanda B, Breiman RF, Shimoni bat virus, a new representative of the Lyssavirus genus. Virus Res. 2010;149:197210 and. 10.1016/j.virusres.2010.01.01820138934Ge X, Li Y, Yang X, Zhang H, Zhou P, Zhang Y, Metagenomic analysis of viruses from bat fecal samples reveals many novel viruses in insectivorous bats in China. J Virol. 2012;86:462030 and. 10.1128/JVI.06671-1122345464Wu Z, Ren X, Yang L, Hu Y, Yang J, He G, Virome analysis for identification of novel mammalian viruses in bat species from Chinese provinces. J Virol. 2012;86:109991012 and. 10.1128/JVI.01394-1222855479Donaldson EF, Haskew AN, Gates JE, Huynh J, Moore CJ, Frieman MB. Metagenomic analysis of the viromes of three North American bat species: viral diversity among different bat species that share a common habitat. J Virol. 2010;84:1300418 and. 10.1128/JVI.01255-1020926577Li L, Victoria JG, Wang C, Jones M, Fellers GM, Kunz TH, Bat guano virome: predominance of dietary viruses from insects and plants plus novel mammalian viruses. J Virol. 2010;84:695565 and. 10.1128/JVI.00501-1020463061Kuzmin IV, Hughes GJ, Rupprecht CE. Phylogenetic relationships of seven previously unclassified viruses within the family Rhabdoviridae using partial nucleoprotein gene sequences. J Gen Virol. 2006;87:232331 and. 10.1099/vir.0.81879-016847128Aznar-Lopez C, Vazquez-Moron S, Marston DA, Juste J, Ibáñez C, Berciano JM, Detection of rhabdovirus viral RNA in oropharyngeal swabs and ectoparasites of Spanish bats. J Gen Virol. 2013;94:6975 and. 10.1099/vir.0.046490-023100368