key: cord-0797895-osbxayk6
authors: Bittar, Cíntia; Machado, Rafael Rahal Guaragna; Comelis, Manuela Tosi; Bueno, Larissa Mayumi; Beguelini, Mateus Rodrigues; Morielle-Versute, Eliana; Nogueira, Maurício Lacerda; Rahal, Paula
title: Alphacoronavirus Detection in Lungs, Liver, and Intestines of Bats from Brazil
date: 2019-05-29
journal: Microb Ecol
DOI: 10.1007/s00248-019-01391-x
sha: f0867137aaf3f96f4c7b22c41dd0667c007befd5
doc_id: 797895
cord_uid: osbxayk6

Bats are flying mammals distributed worldwide known to host several types of Coronavirus (CoV). Since they were reported as the probable source of spillover of highly pathogenic CoV into the human population, investigating the circulation of this virus in bats around the world became of great importance. We analyzed samples from 103 bats from two distinct regions in Brazil. Coronavirus from the Alphacoronavirus genus was detected in 12 animals, 11 from São José do Rio Preto—SP region and 1 from Barreiras—BA region, resulting in a prevalence of 17.18% and 2.56% respectively. The virus was detected not only in intestines but also in lungs and liver. Phylogenetic analysis based on nsP12 genomic region suggests that the sequences group according to host family and sampling location. Studies on the circulation of these viruses in bats remain important to understand the ecology and evolutionary relationship of these pathogens. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00248-019-01391-x) contains supplementary material, which is available to authorized users.

Coronaviruses (CoVs) are single-stranded positive sense RNA viruses that belong to the Coronaviridae family. The family comprises four genera: Alphacoronavirus (α-CoV), Betacoronavirus (β-CoV), Gammacoronavirus (γ-CoV), and Deltacoronavirus (δ-CoV) [1] . These viruses infect several species of animals and can cause hepatic, enteric, neurological, and acute and chronic respiratory infection [2] . In humans, they generally cause mild upper respiratory tract disease with low mortality [3] . In 2003 and 2012 respectively, the emergence of highly pathogenic β-CoV SARS-CoV (severe acute respiratory syndrome) and MERS-CoV (Middle East respiratory syndrome) revealed that this group of viruses can also cause severe respiratory illness in humans [4, 5] .

Bats are widespread flying mammals that are distributed worldwide. About 1150 bat species have been described and 179 were reported in Brazil [6] [7] [8] . Numerous viruses have been reported to infect bats, and these animals are being studied as potential reservoir for several zoonotic diseases, including Coronaviruses [3, 9, 10] . The majority of Bat Coronavirus (BtCoV) have been identified in bats from Asia, Africa, and Europe [11] [12] [13] [14] [15] [16] [17] [18] [19] ; however, CoV has also been reported in bats from South American countries including Brazil, Trinidad, Costa Rica, Ecuador, Panama, and Mexico [20] [21] [22] [23] [24] [25] .

Although there is no report of CoV pathogenesis in bats, studies suggest that SARS-CoV and MERS-CoV have a zoonotic origin in bats [26] . Therefore, studying the circulation of BtCoV in bats around the world is of great importance to Cíntia Bittar and Rafael Rahal Guaragna Machado contributed equally to this work.

Electronic supplementary material The online version of this article (https://doi.org/10.1007/s00248-019-01391-x) contains supplementary material, which is available to authorized users. detect potential sources of spillover of pathogenic viruses into the human population. Here, we investigate the presence of coronavirus in bats from two distinct regions in Brazil.

Bats were collected in São José do Rio Preto region located in southeastern Brazil at the state of São Paulo (SP) and in Barreiras region, located in the northeastern part of the country at the state of Bahia (BA). The first is a transition zone between the Atlantic Forrest biome, composed by an ecosystem of semi-deciduous forests, and the Cerrado biome characterized by savannahs. While the second region is situated within the Cerrado biome. Both regions are under tropical climates with dry winters (Tropical Aw according to Köppen-Geiger). The animals were collected in urban and peri-urban locations, using mist nets and euthanized according to the Guide for the Care and Use of Laboratory [28] [29] [30] . Intestines, liver, and lungs were removed from bats from São José do Rio Preto-SP and intestines and liver from bats collected in Barreiras-BA. All samples were stored at − 150°C until further processing.

Samples were thawed, immersed in TRIzol (Thermo Fisher Scientific), and tissues were homogenized mechanically by a polytron homogenizer (Marconi, Piracicaba, SP BR). Extraction was carried out according to manufacturer's instructions. RNA pellet was resuspended in water treated with DEPC (Sigma Aldrich), quantified by NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific), and stored at − 150°C. The cDNA synthesis was performed using High Capacity cDNA Reverse Transcription Kit (Applied Biosystems).

The presence of Coronavirus in the organs was tested by a pan-Coronavirus NESTED-PCR as described by Chu et al. targeting the RNA-dependent RNA polymerase (RdRp), nsP12 region [31] . Positive samples were purified using DNA Clean & Concentrator (Zymo Research) and sequenced by Sanger method [32] , using the set of primers from the NESTED reaction, in a 3130 XL Genetic Analyzer (Applied Biosystems). Sequences were submitted to Electropherogram Quality Analysis tool from Embrapa Genetic Resources and Biotechnology website (http://asparagin.cenargen.embrapa.br/ phph/). The platform performed the quality analysis using Phred [33] and contig assembly by CAP3 [34] .

Sequences are identified by BtCoV, followed by the region where the host was collected (BA-Barreiras; SJRP-São José do Rio Preto), the reference of the animal (e.g., 03, 04 …) and by the identification of the organ where the viral RNA was detected (I-intestines; Li-liver; Lu-lungs). Nucleotide sequences were submitted to GenBank under the accession numbers MH974764 to MH974780, and specific identification is available in Online Resource 1.

Two datasets were assembled to perform the phylogenetic analysis. The first one (dataset 1) was composed by all the sequences generated in this study plus sequences representing all four Coronavirus genera (Online Resource 1): α-CoV, β-CoV, γ-CoV and δ-CoV. The second dataset (Dataset 2) is composed by our sequences and sequences available in GenBank of the nsP12 region of BtCoV from the Alphacoronavirus genus (Online Resource 2). Sequences were aligned using Muscle [35, 36] , nested in the Seaview 4.6.1 package [37] . The datasets were edited using BioEdit 7.2.5 package [38] .

A Bayesian Markov Chain Monte Carlo (MCMC) estimation of the phylogenetic trees based on the sequences in both datasets was done using BEAST v1.8 [39] . The analysis was carried out using General Time Reversible model plus gamma distribution (GTR + G) for dataset 1 and GTR + G + invariable sites (GTR + I + G) for dataset 2, as selected by jModelTest [40] , and both runs under a lognormal relaxed clock and constant population size. For dataset 1, the MCMC chain ran for 20,000,000 steps sampling every 2000 steps and for dataset 2 for 100,000,000 steps sampling every 5000 steps. The phylogenetic trees were summarized using Tree Annotator 1.8 and edited using iTOL 4.2.3 [41] [42] [43] .

Association between host traits and tree topology was statistically tested based on the posterior set of trees (PST) generated by BEAST, using the software BaTS 0.9 [44] . The Parsimony Score (PS), Association Index (AI), and Monophyletic Clade (MC) metrics were analyzed, and values of p ≤ 0.01 were considered significant.

The map presented in Fig. 1 was designed by the authors using QGIS 3.2.3 based on maps of cities and states provided by Brazilian Institute of Geography and Statistics (IBGE) [45, 46] . Vectorial images of animals used in Figs. 1 and 2 were provided by vecteezy.com, www.freepik.es, and pixabay.com.

A total of 103 bats were collected, 64 from São José do Rio Preto region and 39 from Barreiras region. In both locations, males were more sampled than females, 60.9% and 66.7% respectively. In São José do Rio Preto, the species that was most sampled was Molossus rufus (21.9%) followed by Molossus molossus and Artibeus literatus (20.3% each). In Barreiras, Molossus molossus (28.2%) was most sampled followed by Artibeus literatus and Artibeus planirostris (23.1% each). Online Resource 3 present the information on sex, family, species, and organs of all animals collected.

The presence of CoV RNA was tested by reverse transcription NESTED-PCR specific for the nsP12 genomic region. From bats collected in São José do Rio Preto region intestines, liver, and lungs were tested. As for the animals from Barreiras region, only intestines and liver were tested, due to availability. Positive samples were sequenced and the resulting sequences were submitted to BLASTn (https://www.ncbi.nlm. nih.gov/BLAST/) in order to confirm CoV identification. The results revealed that 12 animals were infected with coronavirus, 11 from São José do Rio Preto and 1 from Barreiras. Viral RNA was detected in more than one organ in some animals totalizing 17 positive samples. Only one animal, SJRP03, tested positive in all organs. Coronavirus was mostly present in intestine samples, except for SJRP06 and SJRP10. Table 1 displays information on family, species, gender, and the organs where the virus was detected.

Considering the regions separately, the prevalence of CoV in bats from São José do Rio Preto region was 17.18% (11/64) and 2.56% (1/39) in bats from Barreiras region. Also, males 

A total of 17 sequences, from 12 different animals, were generated in this study. The first analysis we performed was to determine the phylogenetic relationship of the sequences identified in both regions with representatives of the four Coronavirus genera. Based on dataset 1, composed by 52 nucleotide sequences of 426 bp (see BMethods^), we ran a Bayesian phylogenetic analysis. The phylogenetic tree indicates that the 17 sequences from this study belong to the α-CoV genus since they group in a monophyletic branch together with all sequences from this genus, with strong branch support (1.00) (Fig. 2) .

We decided to analyze our sequences in the context of Bat Coronaviruses from the α-CoV genus. To this end, we assembled dataset 2 composed by 139 nucleotide sequences of 338 bp of the nsP12 region, which included all sequences from this group available in GenBank, plus our sequences. A Bayesian MCMC analysis using BEAST 1.8 was performed. The summarized tree is presented in Fig. 3 . The analysis of the phylogenetic tree revealed a tendency of viruses isolated from the same host family to group in monophyletic branches with strong branch support, although there were some exceptions. In order to test for this association, we performed a statistical analysis, based on the posterior set of trees using BaTS software. According to the results, there is a positive association between the clustering in the phylogenetic tree and bat families with a p value ≤ 0.01 for PS and AI. As for MC values, all families displayed significant 206 Bittar C. et al. We also tested for an association between phylogeny and the location where the virus was sampled. The results revealed a positive association between this trait and phylogenetic clustering based on AI and PS metrics. The MC metrics was significant for all locations with more than two representatives (Table 3 ) (Online Resource 4B). (Sampling site of each sequence used in the analysis is provided in Online Resource 2).

The outbreaks of the β-CoV SARS-CoV and MERS-CoV evidenced the importance of bats as natural reservoir and as a source of spillover of Coronaviruses into the human population [26] . Understanding the diversity and the ecological aspects of these viruses in bats is important to prevent new outbreaks. Here, we investigated the presence of Coronavirus RNA in bats, living in urban and peri-urban areas, from two geographically distant regions in Brazil, São José do Rio Preto-SP and Barreiras-BA. Our results revealed the presence of Bat Coronavirus from the α-CoV genus in specimens from both regions. None of them are closely related to pathogenic viruses to humans. A prevalence of 17.2% was observed among bats sampled from the São José do Rio Preto region which is similar to the rate of 16.7% reported in bats from Hong Kong [47] . Although one study in the UK reported a comparably high BtCoV prevalence of 23% [48] , reported prevalence rates usually range from 3% to 10% [13, 16, 17, 22, 24, [49] [50] [51] [52] [53] , similar to our current findings in the Barreiras region (2.56%). A global study with over 12,000 samples reported a BtCoV prevalence of 8.6% [54] . The CoV prevalence among bats from the São José do Rio Preto region is therefore higher than estimated for Barreiras as well as compared to other studies. Several factors might influence the estimates, like the number of specimens sampled, host species diversity, number of sites of collection, geographic regions, among others [54] . In this study, 54.5% (6/11) of the infected specimens from São José do Rio Preto region were collected in the same site. The location is a private property known to harbor about six different species of bats in the roof of the house. The fact that these animals were collected in the same location might have influenced the high CoV prevalence in this region, since they live in close contact aiding the transmission. Another factor that might have influenced is that in addition to intestines, we also analyzed lungs and liver for the presence of CoV, while most prevalence studies only test fecal or intestine samples. Although most of the infected animals tested positive in the intestines (81.8%-9/11), in two animals, CoV was only detected in other organs (lungs and liver), raising the chances of viral detection. Other authors have also reported the presence of CoV RNA in tissues and secretions other than intestines or fecal samples, such as kidney, heart, lungs, and oropharyngeal swabs [25, 55, 56] . However, this is the first report of BtCoV detection in the liver.

Since the animals presented no clinical symptoms, the significance of the detection of virus in lungs and liver remains unclear. Experimental infection of bats with CoV from the Betacoronavirus genus has been tried by Watanabe and colleagues [57] . Results revealed that the virus was only detected in the intestines. In our study, the infection of lungs and liver was not restricted to one host species or family since it was found in Molossus rufus (Molossidae), Eptesicus sp. (Vespertilionidae), and Glossophaga soricina (Phyllostomidae). We analyzed the phylogenetic relationship of the virus identified in this study with other BtCoV from the α-CoV genus (Fig. 3) . The sequences generated in our study are distributed in three distinct monophyletic branches. Most sequences, 13 out of 16, of CoV identified in bats from São José do Rio Preto grouped into a monophyletic branch with strong posterior probability (0.77). In the same branch, there are also three sequences of CoV identified in bats from the cities of Araçatuba (KC886321) and Penapolis (KT717388, KT717389), about 100 km distant from São José do Rio Preto. This cluster shares a most recent common ancestor, also with high branch support, with the virus we identified in a Molossus molossus bat in Barreiras . Sequences BtCoV_SJRP08I and BtCoV_SJRP22I, both from Phyllostomus discolor bats, group together in another monophyletic branch sharing a most recent common ancestor with a CoV from a P. discolor (JQ731783) bat collected in Panama [23] . Both clusters include all strains, included in the analysis, identified in the respective species. This indicates that the BtCoV identified in these hosts could be species-specific, being more closely related to the viruses identified in the same host, even if they are geographically distant. In fact, other sequences from our study, from hosts collected in the same sites, were infected with BtCoV that grouped apart from the P. discolor and A. literatus clusters ( Fig. 3 ; Online Resource 2). An alternative explanation is that the occurrence of related strains in distant locations is a consequence of transmission during bat migration. Given the genetic distance between the strains sampled, this would obviously not be due to direct transmission but would presumably involve unsampled intermediate strains. Several bat species are known to migrate seasonally following resources, such as food.

Migration studies report that small-and medium-sized bats can migrate long distances, with a maximum migration distance of 1905 km reported for a Vespertilionidae bat [58] . However, there are no studies of migration patterns for P. discolor and few for A. literatus, none of which report such long distance [59, 60] . It is therefore more likely that the close relatedness of BtCoV from P. discolor and from A. literatus is a reflection of the virus Bhost specificity.Ô verall, the BtCoV phylogenetic tree shows highly supported branches clustering primarily according to host family. Some families have two different clusters like Molossidae, Phyllostomidae, and Rhinolophidae while Miniopteridae and Hipposideridae present only one each. Since there are some exceptions and CoV from Vespertilionidae bats does not seem to group together, we statistically tested the hypothesis of a positive relationship between host family and phylogeny and it was confirmed. However, due to the exceptions, we decided to test for a possible association with location where the hosts were collected. For example, as previously described, most CoV from São José do Rio Preto grouped into a monophyletic and most of them belong to the Molossidae family. Still, among them there are sequences from virus identified in two animals from other families, Vespertilionidae (BtCoV_SJRP06) and Phyllostomidae (BtCoV_SJRP10). So, we statistically tested for an association between the location where the host was collected and phylogeny and there is also a positive association. Our results indicate that despite the association between BtCoV and host family seen here, the close contact between animals from different families may result in virus transmission. We conclude that both traits, host family and location, influence BtCoV phylogenetic relationship.

Here, we report the detection of Alphacoronavirus in bats from Brazil, not only in intestines but also in lungs and liver. We found a higher prevalence of BtCoV in animals collected in São José do Rio Preto-SP region, when compared to Barreiras-BA. The investigation of CoV in bats throughout the world remains important to characterize and understand the circulation of these viruses, identifying potential sources of spillover of pathogenic strains into the human population. 

The authors declare that they have no conflict of interest.

Ethical Approval All applicable international, national, and/or institutional guidelines for the care and use of animals were followed. All procedures performed in studies involving animals were in accordance with the ethical standards of the institution or practice at which the studies were conducted. This article does not contain any studies with human participants performed by any of the authors.

Ratification vote on taxonomic proposals to the International Committee on Taxonomy of Viruses

Coronavirus pathogenesis and the emerging pathogen severe acute respiratory syndrome coronavirus

Coronaviruses: important emerging human pathogens

Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia

Identification of a novel coronavirus in patients with severe acute respiratory syndrome

Mammalogy: adaptation, diversity, ecology

Checklist of Brazilian bats, with comments on original records

Mammal species of the world: a taxonomic and geographic reference

A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special?

Bats: important reservoir hosts of emerging viruses

Genomic characterization of severe acute respiratory syndrome-related coronavirus in European bats and classification of coronaviruses based on partial RNA-dependent RNA polymerase gene sequences

Amplification of emerging viruses in a bat colony

Detection and prevalence patterns of group I coronaviruses in bats

Distant relatives of severe acute respiratory syndrome coronavirus and close relatives of human coronavirus 229E in bats

Identification of a novel coronavirus in bats

Coronaviruses in bent-winged bats (Miniopterus spp

Prevalence and genetic diversity of coronaviruses in bats from China

Detection of novel SARS-like and other coronaviruses in bats from Kenya

Identification of a severe acute respiratory syndrome coronavirus-like virus in a leaf-nosed bat in Nigeria

Coronaviruses in bats from Mexico

Genetic diversity of bats coronaviruses in the Atlantic Forest hotspot biome, Brazil

Highly diversified coronaviruses in neotropical bats

Alphacoronavirus in urban Molossidae and Phyllostomidae bats, Brazil

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

Bat origin of human coronaviruses

Guide for the care and use of laboratory animals

Chave para determinação de quirópteros brasileiros

A new species of Eptesicus (Mammalia: Chiroptera: Vespertilionidae) from the Atlantic Forest

A revision of the Neotropical bats of the genus Myotis

Avian coronavirus in wild aquatic birds

DNA sequencing with chain-terminating inhibitors

Base-calling of automated sequencer traces using phred. II. Error probabilities

CAP3: a DNA sequence assembly program

MUSCLE: a multiple sequence alignment method with reduced time and space complexity

MUSCLE: multiple sequence alignment with high accuracy and high throughput

SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building

BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic acids symposium

Bayesian phylogenetics with BEAUti and the BEAST 1.7

2: more models, new heuristics and parallel computing

Interactive Tree Of Life (iTOL): an online tool for phylogenetic tree display and annotation

Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees

Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy

Correlating viral phenotypes with phylogeny: accounting for phylogenetic uncertainty

QGIS Development Team -QGIS Geographic Information System. Open Source Geospatial Foundation Project

IBGE -Instituto Brasileiro de Geografia e Estatísticabases cartográficas

Complete genome sequence of bat coronavirus HKU2 from Chinese horseshoe bats revealed a much smaller spike gene with a different evolutionary lineage from the rest of the genome

Alphacoronavirus detected in bats in the United Kingdom

Alphacoronaviruses detected in French bats are phylogeographically linked to coronaviruses of European bats

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

Detection of coronaviruses in bats of various species in Italy

Detection of group 1 coronaviruses in bats in North America

Diversity of coronavirus in bats from Eastern Thailand

Global patterns in coronavirus diversity

Longitudinal surveillance of Betacoronaviruses in fruit bats in Yunnan Province

Virome analysis for identification of novel mammalian viruses in bats from Southeast China

Bat coronaviruses and experimental infection of bats, the Philippines

Bat migrations in Europe: a review of banding data and literature. Federal Agency for Nature Conservation

Longdistance movement by a great fruit-eating bat, Artibeus lituratus (Olfers, 1818), in southeastern Brazil (Chiroptera, Phyllostomidae): evidence for migration in Neotropical bats?