key: cord-261547-8tfbhmzo
authors: Góes, Luiz Gustavo Bentim; Campos, Angélica Cristine de Almeida; Carvalho, Cristiano de; Ambar, Guilherme; Queiroz, Luzia Helena; Cruz-Neto, Ariovaldo Pereira; Munir, Muhammad; Durigon, Edison Luiz
title: Genetic diversity of bats coronaviruses in the Atlantic Forest hotspot biome, Brazil
date: 2016-07-26
journal: Infect Genet Evol
DOI: 10.1016/j.meegid.2016.07.034
sha: 
doc_id: 261547
cord_uid: 8tfbhmzo

Bats are notorious reservoirs of genetically-diverse and high-profile pathogens, and are playing crucial roles in the emergence and re-emergence of viruses, both in human and in animals. In this report, we identified and characterized previously unknown and diverse genetic clusters of bat coronaviruses in the Atlantic Forest Biome, Brazil. These results highlight the virus richness of bats and their possible roles in the public health.

Brazil harbours 15% of the world's bat diversity and carries 178 distinct bats species (Nogueira et al., 2014) . Out of these, a total of 113 species exist only in the Atlantic Forest Biome (AFB), which is the second largest rain forest of the South America, and is one of the unique regions with highest biodiversity in the world (Paglia et al., 2012) .

Bats are historically unique and widespread mammals playing essential roles in the emergence and re-emergence of viruses of both veterinary and public health importance. Although viruses of diverse genetic backgrounds can co-exist asymptomatically in bats, majority of these viruses are single stranded RNA viruses (Calisher et al., 2006) . Coronaviruses (CoVs) are enveloped, positive-sense, single-stranded RNA viruses in the family Coronaviridae, and are usually associated with respiratory, enteric, hepatic and neurological pathologies of varying severity (Woo et al., 2012) . Coronaviruses are classified into four genera; Alphacoronavirus (α-CoV) and Betacoronavirus (β-CoV) have been exclusively identified in mammals, whereas Gammacoronavirus and Deltacoronavirus are mainly detected in avian species (Woo et al., 2012; ICTV, 2015) . Based on the genetic relatedness, the β-CoVs can further be subdivided into four clades: A to D (Drexler et al., 2014) .

All CoVs that can potentially infect human were originated from animal reservoirs, and four of such CoVs are believe to be transmitted to human through bats including α-CoVs (229E and NL63) and highly pathogenic β-CoV (Severe Acute Respiratory Syndrome and Middle East Respiratory Syndrome) (Bolles et al., 2011; Chan et al., 2015; Corman et al., 2013; Huynh et al., 2012) . Recently, a number of novel bats CoVs have been identified, primarily from African, Asian and European bats (Calisher et al., 2006; Chu et al., 2006; Drexler et al., 2014) , as well as from South American countries including Costa Rica, Panama, Ecuador, Mexico and Brazil (Corman et al., 2013; Goes et al., 2013) . Collectively, these studies indicate the co-existence of bats and viruses at the interface of viral evolution and bats ecology.

A limited number of studies have been conducted in Brazil to map the nature and breath of bats in harbouring viral populations. In previous studies, a total of five distinct CoV lineages have been detected in just 10% of local bats (15 species), and most of these are belonging to α-CoV (Brandao et al., 2008; Corman et al., 2013; Goes et al., 2013) . These attributes, and the existence of a large number of human beings (120 million) in the Atlantic Forest Biome, Brazil, clearly highlight the potential of bats in not only carrying zoonotic viruses but also in possible transmission of viruses to human beings.

To ascertain the diversity of bats CoVs circulating in the Brazilian bats, a total of 401 intestine tissues from 17 bat species were collected from 2010 to 2014 encompassing bats with distinct diet habits (Table  1 , Fig. 1A , B). Bats from urban (N = 192) and rural (N = 10) areas were received from 14 municipalities of the Northwestern State of São Paulo by the Rabies Laboratoryof Universidade Estadual Paulista (UNESP), Araçatuba-SP. Additionally, 99 and 100 samples were collected from Iguaçu National Park, and from two distinct disturbed landscape sites, respectively, and provided by our collaborators of the Zoology Department at UNESP, Rio Claro-SP. All activities were Infection, Genetics and Evolution 44 (2016) 

Infection, Genetics and Evolution j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / m e e g i d authorized and approved by the Ethic Committee of Institute of Biomedical Research from University of São Paulo (56-18-03/2014) and bats species were identified based on the morphological characteristics including head-boby, forearm size, and dental arch as described by Vizotto and Tadei (1973) ; Gregorin and Taddei (2002) and Miranda et al. (2011) .

Nucleic acids extracted from 30 mg of intestine tissue of individual bat using NucliSENS® easyMAG® were subjected to first strand cDNA synthesis with random hexamers and RT-PCR High-Capacity cDNA Archive Kit (Applied BioSystems). The cDNA preparations were screened by pancoronavirus Nested PCR assay targeting the RNA-dependent RNA polymerase (RdRp) gene (Chu et al., 2011) . Positive samples were sequenced and final sequences were submitted to the GenBank (KT717380-KT717394) (Supplementary table) . A dataset, consisting of sequence generated in this study and all publically available bats CoVs, was phylogenetically analysed by neighbour-joining method in MEGA6 software using Kimura's two-parameter correction and 10,000 bootstrap values (Tamura et al., 2013) .

Coronaviruses were detected in 15 bat intestines samples from eight bat species with distinct diet habit, demonstrating a marked potential of CoVs distribution among bat species in AFB that harbours 9% of world's bat diversity. All CoVs-positive bat species were geographical distributed in the Neotropics and in anthropogenic area apparently affected by fragmented forests, and bats-abundance urbanized areas. The detection of CoVs varied by bat species; undisturbed forest remnant (4%), fragmented forests (4.5%) and urban areas (3.1%) ( Table 1) .

Phylodynamics analyses indicated the circulation of 13 α-CoV and 2 β-CoV in AFB bats (Fig. 1C) . This is of primary interest as majority of bat CoVs surveillance studies conducted in the New World bats detected only α-CoVs (Corman et al., 2013; Drexler et al., 2014; Osborne et al., 2011) . A higher resolution analysis indicated a distinct distribution and diversity of α-CoV and β-CoV lineages (Fig. 1D) .

α-CoV sequences obtained from bats of same genus presented high nucleotide sequence similarity (e.g. Artibeus, Glossophaga, Carollia, Molossus, Myotis and Sturnira) (Fig. 1D and Supplementary table), even with sequences detected in other studies from bats of geographically distant regions. This relation can be exemplified by the high similarity of the CoV RNA partial sequences (91,9%) detected in Carollia perspicillata species from Fenix-PR, Brazil (KT717385), and Fyzabad, Trinidad and Tobago (EU769557) (Carrington et al., 2008) , located at a distance of at least 3900 km. Similar results were previously reported for a variety of bat CoVs and are taken as evidence of co-evolution of CoV genotypes and specific host genera (Drexler et al., 2014; Corman et al., 2013; Anthony et al., 2013) .

Moreover, two previously uncharacterized α-CoV from Myotis nigricans and Myotis riparius and one β-CoV cluster were identified in studied bats ( Fig. 1C and D) . CoV sequences from Myotis bats genera presented high nt sequence similarity (98.9%) and were most closely related to bat CoVs (KC110771) collected from Brazilian Molossus molossus and Tadarida brasiliensis. A distinct cluster of α-CoVs were also detected in Sturnira bat that grouped with α-CoVs lineage 1, a group with an evolutionary history of recombination and cross-species transmission between domestic and livestock animals, such as feline, canine and swine α-CoVs (Lorusso et al., 2008) . This clade presents nt similarity between 73.9% and 74.4% with transmissible gastroenteritis virus (DQ811789) and feline coronavirus (AF124987) depending on the Sturnira lineage (Supplementary table) . This branching pattern may possibly explain the common-ancestral origin of α-CoVs lineage 1 species from bats and other animal species. However, extensive evolutionary studies on complete genome sequences of these isolates are required to provide information on the virus origin and divergence.

Out of two β-CoVs, the virus from Eumops glaucinus clustered within MERS-CoV containing lineage C. The second β-CoV detected in Artibeus lituratus shown highest similarity (96.4%) with CoV reported from Costa Rica (Corman et al., 2013) and lower similarity (69.3%) with MERS-CoV (Supplementary table) . These sequences showed sequence similarities with β-CoV detected in Pipistrellus bats that poses high homology with human MERS-CoV (JX869059). Notably to observe that the Eumops bat positive for CoV was found on an urban area and was predated by a domesticated cat. It was not possible to rule out any event of virus transmission to human, however, this highlight the potential of domesticated animals in virus transmission and disease dynamics at the virus-animalhuman interface. This is of special importance due to the established role of animals in the transmission of viruses to human being (Johnson et al., 2015) . Taken together, these results represent the first detection of lineage C β-CoV in South American bats. Despite of close relationship between lineage C β-CoVs in Asia, Africa and Europe, cumulative data indicate that this lineage showed a more diversified host family distribution in Americas in Mormoopidae, Phyllostomidae and Molossidae bat species (Anthony et al., 2013; Corman et al., 2013; Goes et al., 2013) . 

We present a great diversity of CoV genotypes and clusters in Brazilian bats, highlighting a biogeographic distribution of bats CoVs in the region. It is indispensable in future to investigate the evolutionary events in genetically diverse bats CoVs using complete genome sequences, and their possible transmission potentials to human being. Although it is not possible to calculate the risk of "spill over" events of Brazilian bats CoVs to humans, our results reinforce the need for expanded and continuing surveillance of CoVs in bat fauna, including those in the AFB regions of Brazil.

Supplementary data to this article can be found online at http://dx. doi.org/10.1016/j.meegid.2016.07.034.

Coronaviruses in bats from Mexico

SARS-CoV and emergent coronaviruses: viral determinants of interspecies transmission

A coronavirus detected in the vampire bat Desmodus rotundus

Bats: important reservoir hosts of emerging viruses

Detection and phylogenetic analysis of group 1 coronaviruses in South American bats

Middle East respiratory syndrome coronavirus: another zoonotic betacoronavirus causing SARS-like disease

Coronaviruses in bent-winged bats (Miniopterus spp

Avian coronavirus in wild aquatic birds

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

Novel bat coronaviruses

Chave Artificial para a Identificação de Molossídeos Brasileiros (Mammalia, Chiroptera)

Evidence supporting a zoonotic origin of human coronavirus strain NL63

Virus Taxonomy

Spillover and pandemic properties of zoonotic viruses with high host plasticity

Chave ilustrada para determinação dos morcegos da Região Sul do Brasil

Checklist of Brazilian bats

Alphacoronaviruses in New World bats: prevalence, persistence, phylogeny, and potential for interaction with humans

Lista Anotada dos Mamíferos do Brasil/Annotated Checklist of Brazilian Mammals

Evolution. An Eocene big bang for bats

MEGA6: molecular evolutionary genetics analysis version 6.0

Chave para determinação de quirópteros brasileiros

Discovery of seven novel mammalian and avian coronaviruses in the genus Deltacoronavirus supports bat coronaviruses as the gene source of Alphacoronavirus and Betacoronavirus and avian coronaviruses as the gene source of Gammacoronavirus and Deltacoronavirus

A: Brazil's map and sites for bats capture/collection. Bats were collected from 17 sites (a-q) in the Atlantic Forest Biome from two adjacent Paraná and São Paulo states, Brazil. B: Based on morphological characteristics, bats were classified into different families (positive families were bold and underlined)

This study was supported by FAPESP (São Paulo Research Foundation) process number 2013/11006-0 and British Council grant number 172710323. Cruz-Neto A.P was supported by FAPESP grant process number 2008/57687-0. We thank Luiz Aurélio de Campos Crispin and Mariana Cristine Pereira de Souza for their contributions to this study.
