key: cord-0912611-52u2sct7
authors: Moreira‐Soto, A.; Taylor‐Castillo, L.; Vargas‐Vargas, N.; Rodríguez‐Herrera, B.; Jiménez, C.; Corrales‐Aguilar, E.
title: Neotropical Bats from Costa Rica harbour Diverse Coronaviruses
date: 2015-02-04
journal: Zoonoses Public Health
DOI: 10.1111/zph.12181
sha: eac77cf93387fc5c58daa8f6b99ba59a772de504
doc_id: 912611
cord_uid: 52u2sct7

Bats are hosts of diverse coronaviruses (CoVs) known to potentially cross the host–species barrier. For analysing coronavirus diversity in a bat species‐rich country, a total of 421 anal swabs/faecal samples from Costa Rican bats were screened for CoV RNA‐dependent RNA polymerase (RdRp) gene sequences by a pancoronavirus PCR. Six families, 24 genera and 41 species of bats were analysed. The detection rate for CoV was 1%. Individuals (n = 4) from four different species of frugivorous (Artibeus jamaicensis, Carollia perspicillata and Carollia castanea) and nectivorous (Glossophaga soricina) bats were positive for coronavirus‐derived nucleic acids. Analysis of 440 nt. RdRp sequences allocated all Costa Rican bat CoVs to the α‐CoV group. Several CoVs sequences clustered near previously described CoVs from the same species of bat, but were phylogenetically distant from the human CoV sequences identified to date, suggesting no recent spillover events. The Glossophaga soricina CoV sequence is sufficiently dissimilar (26% homology to the closest known bat CoVs) to represent a unique coronavirus not clustering near other CoVs found in the same bat species so far, implying an even higher CoV diversity than previously suspected.

Bats are hosts of diverse coronaviruses (CoVs) known to potentially cross the host-species barrier. For analysing coronavirus diversity in a bat species-rich country, a total of 421 anal swabs/faecal samples from Costa Rican bats were screened for CoV RNA-dependent RNA polymerase (RdRp) gene sequences by a pancoronavirus PCR. Six families, 24 genera and 41 species of bats were analysed. The detection rate for CoV was 1%. Individuals (n = 4) from four different species of frugivorous (Artibeus jamaicensis, Carollia perspicillata and Carollia castanea) and nectivorous (Glossophaga soricina) bats were positive for coronavirus-derived nucleic acids. Analysis of 440 nt. RdRp sequences allocated all Costa Rican bat CoVs to the a-CoV group. Several CoVs sequences clustered near previously described CoVs from the same species of bat, but were phylogenetically distant from the human CoV sequences identified to date, suggesting no recent spillover events. The Glossophaga soricina CoV sequence is sufficiently dissimilar (26% homology to the closest known bat CoVs) to represent a unique coronavirus not clustering near other CoVs found in the same bat species so far, implying an even higher CoV diversity than previously suspected.

Approximately 71.8% of emerging infectious diseases have originated from wildlife (Jones et al., 2008) . Bats have been involved in outbreaks of emerging diseases around the world and are known reservoirs of coronaviruses (CoVs) capable of crossing the host-species barrier (Calisher et al., 2006) . One example is the 2002-2003 severe acute respiratory syndrome (SARS) CoV outbreak. After the subsequent finding that bats constitute this virus' natural reservoir (Lau et al., 2005) , several research groups have detected aand b-CoVs in a variety of bat species throughout the world (Drexler et al., 2014) . Also the newly emerged Middle East respiratory syndrome (MERS) CoV depicts limited sequence similarities with CoVs found in bats (Annan et al., 2013) . Around 33% of overall bat species and 75% of the known 113 bat genera are found in the neotropics (Osborne et al., 2011) . Costa Rica exhibits an extraordinary bat diversity with more than 109 recorded species so far, representing 10% of the total bat diversity worldwide (Timm and LaVal, 2000) . In this study, the diversity of bat coronaviruses in neotropical bats was assessed by analysing putative phylogenetic relationships with previously detected RdRp sequences around the world.

For an intensive screening of bats potentially infected with CoV, permission for bat capturing and sampling was obtained from the Institutional Committee of Care and Use of Animals of the University of Costa Rica (CICUA-124-12) according to international animal welfare standards and by University of Costa Rica project B2-285. Bats were sampled between 2012 and early 2014 at 11 distinct locations in Costa Rica (Fig. 1a) (Platt et al., 2000) set at each location in open areas near drinking or feeding sites, to maximize capture success (Newman et al., 2011) . Once captured, bats were collected and transported to the field laboratory. Bats were examined for signs of disease and taxonomically classified (Timm and LaVal, 2000) . Anal swabs or faecal pellets for a total of 421 bats were collected and dissolved in 400 ll of RNA stabilizing reagent (RNAlater â , Life Technologies, Carlsbad, CA, USA) and (Tamura et al., 2013) . To attempt coronavirus isolation, CoV RNA-positive faecal samples were inoculated after debris and bacterial clarification in Desmodus rotundus embryo primary kidney-, fibroblasts-and intestinal-cell-derived confluent monolayers. Briefly, low passage cells were seeded in 24-well plates with D-MEM 10% NCS and incubated at 37°C in a 5% CO 2 atmosphere until 70% confluence. Cell culture medium was then removed, and cells were inoculated with Log 10 diluted samples. After a 2-hour virus adsorption period, D-MEM 2% NCS was added. Cells were daily observed during 22 days to detect any changes in cell morphology (CPE). Cells and supernatants were analysed for CoVs presence by RT-PCR at days 7, 14 and 21 post-inoculation. We sampled 41 neotropical bat species, comprising 24 genera and 6 of the existing 9 neotropical bat families (Supplementary Table S1 ) and representing a wide variety of trophic guilds. The overall finding for the targeted RdRp sequence was approximately 1% (4 of 421 bats). Four bat species Artibeus jamaicensis, Carollia perspicillata, Carollia castanea and Glossophaga soricina were found to be CoV positive, but none showed any obvious signs of disease, consistent with other reports (Tang et al., 2006; Dominguez et al., 2007) . The finding of CoV positivity in these species may be related to the fact that these bats are highly adaptable, resistant to habitat fragmentation and widespread throughout Costa Rica (Medell ın et al., 2000; Timm and LaVal, 2000) . Genera found previously positive such as Mollossus sp (Lima et al., 2013) were negative. This suggests that bat CoV prevalence might be low, although relative sampling for this species in this study was high.

Sequences found (one each from C. perspicillata, A. jamaicensis, G. soricina and C. castanea) were submitted to Genbank (Accession Numbers KC779225, KC779226, KM215146 and KM215147, respectively) .

Phylogenetic analysis of the nucleotide sequences clustered all of the CoVs detected to the a-CoV group (Fig. 1b) . Sequences found did not cluster near so far known pathogenic human CoVs suggesting that the risk of spillover of known human pathogenic CoVs carried by bats is low. Our C. perspicillata CoV and C. castanea sequences (KC779225 and KM215147, respectively) clustered with previously published Carollia CoVs from Brazil, Costa Rica and Trinidad with 99% posterior probability. The sequence from C. perspicillata has 96.8-97.8% nucleotide sequence identity and shows 1 or 2 amino acid differences (98.6-99.3% amino acid sequence identity for the RdRp) with the Carollia sp. sequences described by Corman et al. (2013) . The sequence of C. castanea CoV has a nucleotide sequence identity of 92.5-92.9% with the other Carollia sp. CoVs detected in the country by us and by Corman et al. (2013) , which represents an amino acid identity of 98.6% with only two amino acid substitutions. This is the first detection of CoV in C. castanea so far reported. All Carollia sp. CoV sequences found to date cluster together in a genus-specific monophyletic clade. The A. jamaicensis CoV sequence (KC779226) was most closely related to A. jamaicensis CoV II found in Panama by Corman et al. (2013) (JQ731784), having 97.3% nucleotide sequence identity and clustering together on a branch with 100% posterior probability. The observed nucleotide sequence difference translates into 4 amino acid substitutions resulting in an amino acid identity of 96.8%. Finding the same species of bats carrying similar CoVs even when geographically distant suggests for species specificity of bat CoVs . Therefore, these findings support the theory that CoVs coevolve with individual bat species or genera.

Interestingly, the G. soricina CoV RdRp sequence (KM2 15146) represents a divergent sequence. In fact, the highest homology (76% sequence identity) corresponds to an Australian Rhinolophus megaphyllus bat CoV (bat coronavirus R.meg/Australia/CoV100/2007 EU834953.1) and one found by us in a Costa Rican A. jamaicensis bat (KC779226). Its sequence identity to CoV sequences from the same species of bat [G. soricina in Trinidad EU769558 (Carrington et al., 2008) ] is only 70% at the nucleotide level and 79.4% at the amino acid level. Posterior probability of 73% makes the placing of this bat CoV sequence in the phylogeny uncertain. Also, our G. soricina bat CoV sequence did not cluster near previously reported CoVs from the same species of bat in Trinidad (Carrington et al., 2008) . Nonetheless, the detection of such a divergent sequence of CoV in G. soricina suggests that the diversity of bat CoVs in a given species may be greater than previously known. It is tempting to speculate that nectivorous bats could be more exposed to different CoVs because of their feeding habits . Additionally, G. soricina is known to roost with bats of eight different families (Alvarez et al., 1991) , which presents a plausible scenario of virus exchange among different species.

Attempts for virus isolation were unsuccessful. No CPE was observed in any of the primary embryonic cell lines, and none of the cell cultures were positive by the CoV-specific RT-PCR. We hypothesize that at least these embryonic bat cell lines originating from D. rotundus are not susceptible to bat CoV infection or that the sample contained noninfectious virus.

Our study emphasizes that even if sampling has been already performed in a given country, continuous surveillance aids in finding novel viruses which may expand the information regarding CoV diversity. The phylogenetic analysis performed in the current study show that the derived bat CoVs are unrelated to coronaviruses from other animals or humans, suggesting that recent zoonotic spillover events have not occurred. For an accurate understanding of the events involving an emergence of a nondescribed pathogen, for example a bat coronavirus, more profound ecological and taxonomical studies in their natural hosts have to be performed. These studies will support a more accurate surveillance and proper conservation management for these mammals. Moreover, because Latin America has been suggested as an emerging disease hotspot (Jones et al., 2008) , more studies on these rich biodiversity countries are encouraged, as it might aid in the forecasting readiness of emerging infectious diseases.

Glossophaga soricina. Mamm. species

Human betacoronavirus 2c EMC/2012-related viruses in bats

Bats: important reservoir hosts of emerging viruses

Detection and Phylogenetic Analysis of Group 1 Coronaviruses in South American Bats

Highly diversified coronaviruses in neotropical bats

Detection of group 1 coronaviruses in bats in North America

Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS

Global trends in emerging infectious diseases

Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats

Detection of Alphacoronavirus in velvety free-tailed bats (Molossus molossus) and Brazilian free-tailed bats (Tadarida brasiliensis) from urban area of Southern Brazil

Bat Diversity and Abundance as Indicators of Disturbance in Neotropical Rainforests

Investigating the Role of Bats in Emerging Zoonoses: Balancing Ecology, Conservation and Public Health Interest. FAO Manual No. 12

Alphacoronaviruses in New World Bats: Prevalence, Persistence, Phylogeny, and Potential for Interaction with Humans

Detection of dengue virus neutralizing antibodies in bats from Costa Rica and Ecuador

Detection of Group 1 Coronaviruses in Bats Using Universal Coronavirus Reverse Transcription Polymerase Chain Reactions. SARS-Other Coronaviruses

MEGA6: Molecular Evolutionary Genetics Analysis version 6.0

Prevalence and genetic diversity of coronaviruses in bats from

Clave de campo para los murcielagos de Costa Rica

We thank all the students who helped with capturing the bats, Carlos Vargas for valuable technical assistance and Emmanuel Rojas of Reserva Biologica Tirimbina. We thank Dr. Victor Hugo Sancho from the Rabies Control Program from the Department of Agriculture of Costa Rica (MAG). We also thank Dr. Mirko Trilling from the Institute for Virology, Essen, Germany, Dr. Cesar Rodriguez Sanchez from the University of Costa Rica and Dr. Christine Carrington from the University of West Indies, Trinidad and Tobago for critical reading of the manuscript. This work was supported by the University of Costa Rica Project VI-803-B2-285.

Conceived and designed experiments: AMS, LTC, BRH, CJ and ECA. Performed the experiments: AMS, NVV and BRH. Analysed the data: AMS, LTC and ECA. Wrote the manuscript: AMS and ECA.

The authors disclose no conflict of interest.

Additional Supporting Information may be found in the online version of this article: Table S1 . Bat species collected, tested and listed according to location and CoV positivity.