Alphacoronavirus Detection in Lungs, Liver, and Intestines of Bats from Brazil.

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.

Introduction

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-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-19]; however, CoV has also been reported in bats from South American countries including Brazil, Trinidad, Costa Rica, Ecuador, Panama, and Mexico [20-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 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.

Methods

Samples

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 Animals of the Institute for Laboratory Animal Resources Commission on Life Sciences National Research Council [27]. Also, collections were made under environmental licenses from the Brazilian Institute of Environment and Natural Renewable Resources (IBAMA), SISBIO no.: 21707-1 (Process: 02027.001957/2006-02); 46190-1 and the project was approved by the Ethics Committee on the Use of Animals of the Institute of Biosciences, Letters and Exact Sciences (CEUA/IBILCE—Protocol 135/2016). Sex was determined based on Vizotto and Tadei, Miranda et al., and Laval [28-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.

RNA Extraction and cDNA Synthesis

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).

Coronavirus Detection and Sequencing

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].

Sequence ID

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.

Datasets and Phylogenetic Analysis

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-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.

Images Design and Attribution

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 , , and .

Fig. 1

Map of Brazil showing the locations where infected bats where collected. The host families of infected animals are represented in the respective sampling location (see legend). São Paulo state—dark blue; São José do Rio Preto region—light blue; Bahia state—orange; Barreiras region—light orange

Fig. 2

Bayesian phylogenetic tree based on a dataset of 52 nucleotide sequences of 426 bp from the nsP12 genomic region of CoV. Nodes display posterior probabilities higher than 0.7. The analysis includes 17 sequences from this study along with representatives of the four CoV genera, represented by different colors according to the legend. Images beside the tips of the tree are a graphic representation of the host. Accession numbers of each sequence is presented in Online Resource 1. Sequences generated in this study

Results

Samples and Coronavirus Detection

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.

Table 1

Specimen, family, species, gender, and the organs that tested positive for Coronavirus

Bat ID	Family	Species	Gender	Intestines	Liver	Lungs
BtCoV_SJRP03	Molossidae	Molossus rufus	Male	+	+	+
BtCoV_SJRP04	Molossidae	Molossus rufus	Male	+	–	–
BtCoV_SJRP05	Molossidae	Molossus rufus	Male	+	+	–
BtCoV_SJRP06	Vespertilionidae	Eptesicus sp	Female	–	+	+
BtCoV_SJRP08	Phyllostomidae	Phyllostomus discolor	Male	+	–	–
BtCoV_SJRP10	Phyllostomidae	Glossophaga soricina	Male	–	–	+
BtCoV_SJRP16	Molossidae	Molossus molossus	Female	+	–	–
BtCoV_SJRP22	Phyllostomidae	Phyllostomus discolor	Male	+	–	–
BtCoV_SJRP44	Phyllostomidae	Artibeus literatus	Male	+	–	–
BtCoV_SJRP49	Molossidae	Molossus rufus	Male	+	–	–
BtCoV_SJRP52	Molossidae	Molossus rufus	Male	+	+	–
BtCoV_BA19I	Molossidae	Molossus molossus	Male	+	–	–

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 were more infected than females since in São José do Rio Preto, CoV was detected in 23.1% (9/39) of males and in only 8% (2/25) of females. In Barreiras—BA, the only bat infected was male. The locations where infected bats were collected are presented in Fig. 1. It also show the families of bats that hosted CoV found in each location.

Phylogenetic Analysis

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 “Methods”), 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 association except for Mormoopidae, since there is only one virus identified in this family in our analysis (Table 2) (Online Resource 4A).

Fig. 3

Bayesian phylogenetic tree based on a dataset of 139 nucleotide sequences of 338 bp from the nsP12 genomic region of CoV . The analysis includes sequences of BtCoV from the α-CoV genus. Nodes display posterior probabilities higher than 0.7. Tips are colored according to the host family (see legend). Sequences generated in this study

Table 2

Statistical analysis of the association of host family and phylogenetic clustering performed on BaTS software. Parsimony Score (PS), Association Index (AI), Monophyletic Clade (MC). Values of p ≤ 0.01 were considered significant

Statistic	Bat family	p value
AI		0.00
PS		0.00
MC	Miniopteridae	0.01
MC	Molossidae	0.01
MC	Vespertilionidae	0.01
MC	Phyllostomidae	0.01
MC	Rhinolophidae	0.01
MC	Mormoopidae	1.0
MC	Hipposideridae	0.01

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).

Table 3

Statistical analysis of the association of location and phylogenetic clustering performed on BaTS software. Parsimony Score (PS), Association Index (AI), Monophyletic Clade (MC). Values of p ≤ 0.01 were considered significant

Statistic	Sampling location	p value
AI		0.00
PS		0.00
MC	Hong_Kong	1.00
MC	Barreiras_BA	1.0
MC	S.J.Rio Preto_SP	0.01
MC	Bulgaria	0.01
MC	Panama	1.0
MC	Porto Alegre_RS	0.01
MC	CostaRica_Tirimbina	1.0
MC	CostaRica_SanJose	1.0
MC	Araçatuba	1.0
MC	Mexico	1.0
MC	China_Henan	1.0
MC	CostaRica_Parrita	1.00
MC	Mojiang	0.01
MC	Ghana	0.01
MC	Fênix_PR	1.0
MC	Iguaçu_PR	1.00
MC	Ivai_PR	0.01
MC	Andradina_SP	1.0
MC	Penápolis_SP	0.01
MC	SudMennucci_SP	1.0
MC	Campinas_SP	0.01
MC	Ribeirão_Preto_SP	1.0
MC	Valinhos_SP	1.0
MC	Sumaré_SP	1.0
MC	China_Guizhou	0.01
MC	China_Guangdong	0.01
MC	Hainan_China	1.0

Discussion

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–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 (BtCoV_BA19I). The sequence BtCoV_SJRP44I groups apart from all other in a monophyletic branch (0.99) together with sequence KT717382. Both viruses were detected in Artibeus literatus bats. It is interesting that KT717382 was identified in Iguaçu National Park—PR, about 850 km from São José do Rio Preto—SP, in 2012 [22]. 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 “host specificity.”

Overall, 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.

GenBank accession number of the sequences from Dataset 1, including the sequences from this study. (XLSX 9 kb)

GenBank accession number, host family and species, date and place of host sampling of the sequences from Dataset 2. (XLSX 16 kb)

Information on sex, family, species and organs of all animals collected in this study. (XLSX 13 kb)

Results of statistical test of association between traits and phylogeny done on BaTS software. A. Statistical analysis of the association of host family and phylogenetic clustering; B. Statistical analysis of the association of location and phylogenetic clustering. Parsimony Score (PS); Association Index (AI); Monophyletic Clade (MC). Values of p ≤ 0.01 were considered significant. (PDF 976 kb)