Dong-Lei Sun
Yi-Zhou Gao
Xing-Yi Ge
Zheng-Li Shi
Ning-Yi Zhou- 1Wuhan Institute of Virology, Chinese Academy of Sciences, Wuhan, China
- 2State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
- 3College of Biology, Hunan University, Changsha, China
Bats (order Chiroptera) are one of the most diverse and widely distributed group of mammals with a close relationship to humans. Over the past few decades, a number of studies have been performed on bat viruses; in contrast, bacterial pathogens carried by bats were largely neglected. As more bacterial pathogens are being identified from bats, the need to study their natural microbiota is becoming urgent. In the current study, fecal samples of four bat species from different locations of China were analyzed for their microbiota composition. Together with the results of others, we concluded that bat microbiota is most commonly dominated by Firmicutes and Proteobacteria; the strict anaerobic phylum Bacteroidetes, which is dominant in other terrestrial mammals, especially humans and mice, is relatively rare in bats. This phenomenon was interpreted as a result of a highly specified gastrointestinal tract in adaptation to the flying lifestyle of bats. Further comparative study implied that bat microbiota resemble those of the order Carnivora. To discover potential bacterial pathogens, a database was generated containing the 16S rRNA gene sequences of known bacterial pathogens. Potential bacterial pathogens belonging to 12 genera were detected such as Salmonella, Shigella, and Yersinia, among which some have been previously reported in bats. This study demonstrated high resolution and repeatability in detecting organisms of rare existence, and the results could be used as guidance for future bacterial pathogen isolation.
Introduction
Bats (common term for order Chiroptera) are one of the most diverse and widely distributed groups of mammals which are capable of sustained flight and have gained important scientific interest in recent decades. Studies on bat-associated viral pathogens regarded bats as natural reservoirs for many emerging viruses (Calisher et al., 2006; Ge et al., 2012), including some of potential threat to human health (Halpin et al., 2000; Leroy et al., 2005; Li et al., 2005). Besides viruses, pathogenic parasites were also found in bats (Murakoshi et al., 2018). In contrast to previous reports on bat viral pathogens (Zhou et al., 2018), studies on bat-associated bacterial pathogens and their natural commensal microbiota are less common (Muhldorfer et al., 2010; Muhldorfer, 2013). Although limited, several potential pathogenic bacterial strains have been isolated or identified in bats over the past few decades, emphasizing the importance of bat bacteria investigations (Muhldorfer, 2013). These potential pathogens include those from genera Salmonella (Adesiyun et al., 2009;Reyes et al., 2011), Shigella (Arata et al., 1968), Yersinia (Muhldorfer et al., 2010), Vibrio (Muhldorfer et al., 2011), Clostridium (Hajkova and Pikula, 2007), Campylobacter (Hazeleger et al., 2010), Listeria (Hohne et al., 1975), Leptospira (Matthias et al., 2005), and Pasteurella (Helmick et al., 2004). Among these species, some can cause diseases to their hosts or other animals while the others did not show obvious pathogenicity and may be considered as opportunistic pathogens (Muhldorfer, 2013).
Besides bacterial pathogens, studies on the natural microbial composition of bats are also scant. A review on bacteria of bats emphasized this point and called for more research efforts on this area (Muhldorfer, 2013). Earlier microbial ecological studies of bat microbiota using culture-based methods have been performed decades ago (Klite, 1965; Heard et al., 1997; Di Bella et al., 2003). The microbial composition of bats may be influenced by their feeding strategies (Carrillo-Araujo et al., 2015; Li et al., 2018) or their ages (Vengust et al., 2018). However, previous studies did not focus on the differences in gut microbiome between bats and land mammals. Moreover, it is now well known that a large number of microbes are unculturable under laboratory conditions, and culture-independent methods have been developed intensively in the last two decades (Muyzer et al., 1993; Ley et al., 2008). Using cloning and sequencing technology, Ley et al. (2008) studied the gut microbiota of 59 mammalian species including two Chiroptera species (Carollia perspicillata and Pteropus giganteus). Due to the low throughput of this technology, they were only able to acquire 274 and 228 sequences from each bat species. Later, a more comprehensive sequencing-based study of bat microbiota was performed by Phillips et al. (2012) on bats from North America. However, they did not compare bat microbiota with other mammals or focus on bacterial pathogens.
In the current study, gut microbiota of four bat species belonging to four different families were investigated. The four species of bats come from different habitats of China, including a cave in mountains in Yunnan (Eonycteris spelaea), a cave in mountains located at Xianning city periphery in Hubei (Hipposideros armiger), a cave in a mountain in Guangxi (Taphozous melanopogon), and a historic building in the city of Tianjin (Myotis petax). Bats of E. spelaea are herbivores feeding on banana, longan, etc., while the other three species are insectivorous. This study describes the general features of bat natural microbiota and compared with microbiota of other terrestrial mammals. Moreover, a database containing 16S rRNA gene sequences of known bacterial pathogens was generated to predict potential bacterial pathogens.
Materials and Methods
Sample Collection and Phylogenetic Determination
Fecal samples of four wild bat species located at different habitats of China were collected in summer without disturbing their natural life (Figure 1). All animals included this study displayed no apparent clinical signs of illness based on their normal behavior as well as smooth fur appearance. The sampling process, in brief, was performed by spreading an aseptic plastic cloth under bat colonies and waiting until sufficient fecal samples to be collected (normally 1–2 h) (Fofanov and Furstenau, 2018). Samples were frozen in liquid nitrogen immediately after collection. Feces from the same bat species were eventually mixed for two reasons. (1) Sampling of the gregarious animal individually without capturing is difficult. (2) Variations within same species can be minimized by mixing, and repeatability of techniques can be measured. Before this experiment, long-term viral monitor research for these bat colonies had been done, and few bats for each colony were captured and morphologically identified. Two hundred pellets of the same bat species were mixed in the lab.
Figure 1. A map showing the four bat sample collection locations in China.
The phylogenetic relationship was determined by sequencing of cytochrome b (Cyto b) gene (please see below). A phylogenetic dendrogram demonstrating the relationship of bat species was constructed using MEGA 5.2 (Tamura et al., 2013). The alignment data sets were used to generate the phylogenetic trees under the maximum likelihood method with a bootstrap value of 1,000 replicates.
DNA Extraction
Genomic DNA from fecal samples was extracted using a QIAamp DNA stool mini kit (Qiagen, Germany) following the manufacturer’s instructions; a bead-beating step was introduced to facilitate bacteria lysis (Yu and Morrison, 2004). Concentration and purity of extracted DNA were tested on NanoDrop (Thermo, United States). For each species, three aliquots of samples (0.2 g each) were processed simultaneously.
DNA Sequencing and Data Analysis
Briefly, a 3-mm wing membrane biopsy specimen from wild bats was obtained for DNA extraction by the DNeasy Tissue Kit (QIAGEN). The primers L14724ag-F (5′-ATGATATGAAAAACCATCGTTG-3′) and H15915ag-R (5′-TTTCCNTTTCTGGTTTACAAGAC-3′) were used for amplification of the mitochondrial Cyto b gene of bat (Guillén-Servent and Francis, 2006). The PCR amplification was performed as follows: 94°C for 4 min followed by 40 cycles consisting of 94°C for 30 s, 60°C for 30 s, 68°C for 1.5 min, and a final extension of 68°C for 10 min. Degenerate primers F515 (5′-GTGNCAGCMGCCGCGGTAA-3′) and R926 (5′-CCGYCAATTYMTTTRAGTTT-3′) (Quince et al., 2011) targeting the V4–V5 hypervariable region of 16S rRNA genes was used as this region has been shown to have the least intragenomic heterogeneity (Sun et al., 2013) and thus is less likely to overestimate microbial diversity. A 10-nucleotide barcode attached to primer F515 was designed to distinguish different samples in multiplexed sequencing (Table 1). Every two barcodes contained at least two different nucleotides. A 50 μl PCR mixture contained 2.5 U rTaq DNA polymerase (Takara, China), 0.2 mM dNTPs, 0.2 μM each primer, and 50 ng of total genomic DNA. PCR was performed under the following conditions: initial denaturation at 95°C for 5 min, followed by 29 cycles at 95°C for 30 s, 55°C for 30 s, and 72°C for 30 s, with a final extension at 72°C for 10 min. PCR products were further purified using Cycle Pure Kit (OMEGA, United States). Sequencing was performed on the ROCHE 454 FLX Titanium platform (Roche, Switzerland) at the National Human Genome Center of China at Shanghai, China. Sequencing began from forward primer F515 so as to capture the V4 and most of the V5 region.
Table 1. Bat species included for pyrosequencing, barcodes and number of reads for each sample.
Sequencing raw data (.sff file) was first denoised using the PyroNoise algorithm (Quince et al., 2009) and assigned to each sample according to the 10-nt barcode. Reads without a valid barcode or shorter than 300 bp were excluded from further study. Chimeras were identified and removed using the UCHIME program embedded in MOTHUR (Schloss et al., 2009). Operational taxonomic units (OTUs) were clustered at 3% and 5% dissimilarity levels using the furthest neighbor method and were further used to generate rarefaction curves and estimate richness and diversity indices. The taxonomic assignment of reads was determined by the RDP classifier (Wang et al., 2007) at a confidence threshold of 70%. To compare microbial communities between samples, both the OTU-based Bray–Curtis index and phylogenetic-based unweighted UniFrac (Lozupone et al., 2006) were used. Dendrograms were visualized either by Figtree or EvolView (Zhang et al., 2012). Principal coordinate analysis (PCoA) was also performed using the resulting UniFrac distance matrix. A comparison of the difference of communities was performed by analysis of molecular variance (AMOVA, p < 0.05 was considered as significant). The raw data were also transformed into.fastq files with MOTHUR (Schloss et al., 2009) and submitted into GenBank (accession number PRJNA576557). Those who are interested in this can extract effective information with the corresponding barcodes listed in Table 1 from any.fastq files.
Bacterial Pathogen Detection
A database containing 16S rRNA genes of well-identified bacterial pathogens was constructed from existing pathogen databases including the National Microbial Pathogen Database Resource (NMPDR1) and Virulence Factors Database2 (Supplementary Table S1). The genes for 16S rRNA in bacterial pathogens were retrieved from the Ribosomal Database Project database3. Local BLAST was performed using our sequencing data as the database and the pathogen database as the query. To ensure maximum reliability, only results with sequence identity higher than 99% were recorded.
Results
Phylogenetic Analysis of Bat Species
Phylogenetic information deduced from Cyto b gene (Figure 2) showed four bat species each located in a separate cluster. The only frugivorous species E. spelaea (YN), although varies significantly in body size and feeding habit, is more closely related to the insectivorous H. armiger (XN). M. petax (TJ) is more distant from the aforementioned two. The last species T. melanopogon (GX) which was located in a distant clade indicated the furthest evolutionary distance from the others.
Figure 2. Phylogenetic dendrogram constructed from Cyto b gene sequences of bats. According to traditional classification, Eonycteris spelaea belong to the Megachiroptera suborder and the others belong to the Microchiroptera suborder.
Diversity and Richness Estimators Based on Sequencing Data
A total of 126,809 barcoded raw sequences were produced in this study, the ROCHE 454 FLX Titanium platform allowed a mean read length of 421 bp. After strict quality control, trimming of primers and barcodes, and removal of chimeric sequences, 92,714 high-quality sequences with an average of 7,726 reads per sample (range 4,687–8,596, SD = 1,014) proceeded into the downstream process. To minimize its effect on the calculation of diversity and richness estimators, a subsampling procedure was introduced using the smallest sample size and repeated 1,000 times to calculate the average (Schloss et al., 2011). Since a majority of microbial ecological studies have used 3% and 5% dissimilarity levels to define OTUs (McKenna et al., 2008; Nam et al., 2011), this criterion was also used in the current study. At the 3% and 5% dissimilarity levels, the observed OTUs for each species averaged 427 (3% dissimilarity level) (range 346–525, SD = 58) and 197 (5% dissimilarity level) (range 148–256, SD = 36), respectively (Table 2), among which T. melanopogon (GX) showed the greatest number of observed OTUs and M. petax (TJ) samples showed the least number of observed OTUs at both levels. Chao1 and abundance-based coverage estimator (ACE) estimates of richness were consistent with observed OTUs (Table 2); however, the Shannon diversity index indicated that the frugivorous species of E. spelaea of the YN samples (averaged 4.33) was the highest. The result is not surprising because Shannon diversity is a reflection of both richness and evenness. Although YN samples did not show the highest richness, it did have the highest evenness (Table 2). Analysis of rank abundance showed the bat gut microbiota consisted of few highly abundant species and a majority of rare species (Supplementary Figure S1).
Table 2. Alpha diversity of each species at 3% and 5% dissimilarity levels.
To study whether the sequencing effort is sufficient to recover most of the colony members, a rarefaction analysis was performed based on OTUs clustered at 3% and 5% cutoff (Figures 3A,B). These curves indicated that the sequencing depth of the current study well represented the microbial communities as curves of all samples were approaching plateaus. This conclusion was further consolidated by another parameter reflecting the sequencing completeness good’s coverage which ranged 94–96% among all samples (Table 2). Coverage of this level means no more than six new OTUs were expected for every 100 additional sequencing efforts.
Figure 3. Diversity of bacterial communities as revealed by a rarefaction curve constructed from OTUs at 3% (A) and 5% (B) dissimilarity levels.
Taxonomic Assignment of Bacterial Pyrosequencing Reads
At a confidence threshold of 70%, 92,525 out of 92,714 (99.7%) high-quality sequences can be assigned to 12 known phyla (Figure 4A). Surprisingly, Firmicutes and Proteobacteria dominated the gut microbiota of all four bat species. Firmicutes, which was also a common composition of human gut flora, took up around one half of the microbial community in GX and YN samples; however, this group was even more dominant in TJ and XN samples with a percentage of more than 90%. Proteobacteria is also commonly discovered in guts of terrestrial mammals (mice, Nava et al., 2011; macaques, McKenna et al., 2008; and humans, Nam et al., 2011), but this phylum is rarely dominant in those mammals (usually less than 5%). In this study, Proteobacteria occupied 44% sequences in GX and 54% in YN and thus can be regarded as dominant. The proportions of this phylum in TJ and XN samples (6% and 9%, respectively) were also higher than those in other terrestrial mammals. To our surprise, another dominant phylum other than Firmicutes in higher mammals, Bacteroidetes, was nearly absent in bats. Only GX samples detected around 2% Bacteroidetes sequences; its proportion in other samples was less than 0.01%. Similarly, TJ samples detected around 1% sequences from Fusobacteria; however, it was hardly detected in the others. Phyla of Actinobacteria, Synergistetes, Tenericutes, Cyanobacteria, Deinococcus–Thermus, Gemmatimonadetes, and TM7 were found to be minority groups.
Figure 4. Taxonomic classification of pyrosequencing reads at phylum level (A) and family level (B) using the RDP classifier at a confidence threshold of 70%.
A more detailed community composition at family level was demonstrated in Figure 4B. GX, XN, and YN samples shared a common feature as they all consisted of a large proportion of sequences from Streptococcaceae, Enterobacteriaceae, and Enterococcaceae. In TJ samples, Clostridiaceae seemed to have taken the place of Streptococcaceae and made up more than 60% of total microbiota. Some species-specific groups can be distinguished at this taxonomic level, for example, Neisseriaceae for GX, Peptostreptococcaceae for TJ, Ruminococcaceae for XN, and Moraxellaceae for YN. Further study of this result revealed, at least in some bat species, that their gut flora tends to be dominated by facultative (e.g., Streptococcaceae, Enterobacteriaceae, and Enterococcaceae) or even aerobic (e.g., Neisseriaceae) microorganisms rather than strict anaerobic groups like Bacteroidetes. A heatmap constructed on the major 30 genera (Supplementary Figure S2) showed that Lactococcus was consistently detected in all samples, which is reminiscent of a previous report about the slight acidity of bat intestinal contents (Klite, 1965). Another genus detected with considerable proportion among all samples was Enterococcus, which is also a common commensal organism of animal guts.
Community Similarity Analysis
Multiple approaches were employed to compare the microbial community similarity of bat species. Dendrogram constructed from the Yue and Clayton measure of dissimilarity and Bray–Curtis index based on OTU information (3% cutoff) both indicated TJ samples were the most divergent from the others; on the other hand, GX and XN samples consistently showed the closest similarity (Figures 5A,B). We also want to question whether the host phylogenetic relationship is in accordance with that of its gut microbiota. However, a comparison between the host evolutionary relationship (Figure 2) and gut flora dendrogram (Figure 5) showed no obvious correlation (congruence test, de Vienne et al., 2007, p > 0.05). Since the above OTU-based similarity analysis did not take into consideration the phylogenetic information, unweighted UniFrac was performed, and PCoA was used to study the communities along axes of maximal variance (Figure 6). It is interesting that four species each located in a separate quadrant, and AMOVA concluded that every between-group difference was significant (p < 0.05). Three replicates strongly clustered together, indicating good repeatability in this study.
Figure 5. Dendrogram showing the similarity of bacterial communities based on the (A) Yue and Clayton measure of dissimilarity and (B) Bray–Curtis index.
Figure 6. Principal coordinate analysis of UniFrac metric for all samples: (A) PC1 and PC2 and (B) PC1 and PC3.
To address the question about the position of bat microbiota in the mammalian world, we took advantage of Ley’s data of mammalian gut microbes (Ley et al., 2008). The gut genome of bats is compared with that of terrestrial animals here. Unweighted UniFrac analysis showed bat microbiota was generally restricted to a cluster which contained species almost entirely from order Carnivora (Figure 7). These species include black bear (BB), polar bear (PB), giant panda (GP), red panda (RP and RPSD), and so on. It should be noted that although some species (e.g., giant panda and red panda) are herbivorous, they evolutionary belong to Carnivora and share a similar gut physiology (simple gut compared to foregut or hindgut) with their carnivorous cousins (Ley et al., 2008). The only exception that was not from Carnivora in the cluster was hedgehog (HH) from order Insectivora; however, it also has a simple gut type similar to that of the above Carnivora species.
Figure 7. Comparison of bat microbiota with that of other mammals (unweighted UniFrac analysis). Reference sequence data from Ley et al. (2008) were trimmed to match the same position and length of the 16S rRNA gene. Different colors represent different feeding habits (green: herbivorous, red: carnivorous, gray: omnivorous) except blue which marked our sample. Abbreviations of species: orangutans, Orang1,2; western lowland gorillas, GOR and GORSD; chimpanzees, CHIMP1,12; black rhinoceros, BNO; flying fox, FF; bighorn sheep, BHSD, BH1, and BH2; horses, HORSEJ and HORSEM; Somali wild ass, WA; Grevy’s zebra, GZ; Asian elephants, AE1-3; rock hyraxes, RH and RHSD; Capybara, CAP; Bwindi gorilla, BG; giraffe, gir; springboks, SBK; Okapi, OK1-3; Transcaspian Urial sheep, TU1,2; Speke’s gazelles, SP3; Argali sheep, AS1-3; takin, TAK; François langur, FL; springboks, SB and SBSD; douc langur, DL; eastern black and white colobus, COL; red kangaroos, KO1,2; red river hog, RRH; baboons, AFBAB; Indian rhinoceros, IR; naked mole rat, MOLERAT; Visayan warty pig, VWP; babirusa, BARB; spotted hyenas, HY1,2; East Angolan colobus, EAC; ring-tailed lemur, RT; white-faced saki, Saki; Prevost’s squirrel, SQ; spider monkey, SPIM; mongoose lemur, ML; black lemur, BKLE; calimicos (Goeldi’s marmoset), CAL; Geoffrey’s marmoset, MAR; bush dogs, bdog1,3; armadillo, arma; lions, Li1-3; cheetahs, CE2,3; African elephants, AFEL, AFYEL, AFEL2,3; Hartmann’s mountain zebra, AFZEB; red panda, RP and RPSD; polar bears, PB1,2; hedgehog, HH; spectacled bear, SB; giant panda, GP; black bears, BB1,2; Seba’s short-tailed bat, bat; echidna, ECH; rabbit, RA.
Potential Bacterial Pathogen Exploration
Since a considerable number of bacterial pathogens carried by bats have been reported (Muhldorfer, 2013), a 16S rRNA gene database constructed from well-characterized bacterial pathogens was used in an effort to search for potential bacterial pathogens in bats. To maximize reliability, a threshold of 99% sequence identity was applied. Potential pathogens from 12 genera were detected with different proportions (Figure 8), among which Salmonella, Staphylococcus, and Bacillus were consistently detected in all bat species. Brucella was detected in three species except TJ. Shigella, Enterococcus, and Escherichia were detected in XN and YN; on the other hand, GX and XN samples showed the existence of Yersinia. Other rarely detected but potentially pathogenic genera include Bartonella, Bordetella, Streptococcus, and Vibrio.
Figure 8. Heatmap of detected potential bacterial pathogens. The strains in brackets represent the pathogenic bacteria used for retrieval and comparison.
Discussion
The order Chiroptera is one of the most diverse and geographically distributed groups of mammals capable of sustained flight (Muhldorfer, 2013). Their wide distribution and close relationship to humans have drawn much research attention over the past few decades (Ciminski et al., 2019; Fan et al., 2019; Laing et al., 2019). A variety of viral pathogens with clinical significance have been identified in bats (Halpin et al., 2000; Leroy et al., 2005; Ge et al., 2012), especially a recently discovered SARS-like coronavirus that might account for the 2003 SARS pandemic in China (Ge et al., 2013). Bats are able to both tolerate viral infections and live well beyond expectations (Hayman, 2019). Compared to the number of studies on bat virus, knowledge about bat bacteria seemed to have lagged behind. Early researches on bat commensal bacteria can be dated back to the 1960s (Klite, 1965); this and later studies of bat microbiota were mostly based on cultivation (Heard et al., 1997). Due to the wide existence of uncultivable microorganisms, those studies were obviously limited. The development of sequencing techniques provided a convenient and powerful approach in studying the composition of extremely complex environments such as the gut (Nam et al., 2011), soil (Roesch et al., 2007), and seawater (Lee et al., 2011).
Traditional opinions classified Chiroptera into Megachiroptera (bigger in size and mostly feed on fruits) and Microchiroptera (echolocating and mostly feeding on insects). From this view, E. spelaea belongs to Megachiroptera while the other three belong to Microchiroptera. However, some evidences suggested Chiroptera be classified into Yinpterochiroptera (including Eonycteris and Hipposideros) and Yangochiroptera suborders (including Taphozous and Myotis) (Springer et al., 2001; Teeling et al., 2002, 2005). Our phylogenetic evidence based on the Cyto b gene sequence supported this theory as E. spelaea first joined H. armiger; these two together clustered with M. petax and T. melanopogon.
The current study revealed that strict anaerobic bacteria from Bacteroidetes were almost absent in bats; on the contrary, facultative or even aerobic species were more abundant. This result is surprising because Bacteroidetes has been shown as a dominating group of commensal community in a majority of terrestrial mammals including humans, mice, and many other species (Suchodolski et al., 2008; Nam et al., 2011; Nava et al., 2011; Shanks et al., 2011). Although it has not been pointed out, the gut microbiota of two bat species in Ley’s study were also significantly dominated by Firmicutes or Proteobacteria and lacking Bacteroidetes (Ley et al., 2008). The result is also in consistency with Phillips’ study as Bacteroidetes only made up 1% of the total microbiota (Phillips et al., 2012). We hereby argue this phenomenon is a coevolution of gut flora and host gut physiology because the gut system of bats has many unique features compared to other mammals. The key features to bat gut include, firstly, shorter intestinal length (about one third to one fifth of that of mouse of comparable weight, Klite, 1965; Caviedes-Vidal et al., 2007); secondly, a unique structure of the gastrointestinal system, usually lacking the cecum and appendix (Klite, 1965; Gadelha-Alves et al., 2008) and in some species even the colon (Makanya and Maina, 1994); and thirdly, significantly reduced transit time of food digest through the intestine (as short as 15 min in some species, Klite, 1965, compared to several hours in mice of comparable size). In addition, as the only mammals with a true flying ability, which is the most energy-consuming activity, bats have well adapted to carrying increased oxygen demands of flight and elevated oxygen diffusion capacity (Canals et al., 2005). The metabolic rate of a flying bat is three to five times higher than that of similar-sized terrestrial mammals in maximum exercise (Shen et al., 2010). Further efforts are required to verify whether the lack of a strict anaerobic Bacteroidetes community is universal to other flying animals.
It has been well established that commensal microbiota is affected by a variety of factors including host diet, phylogeny, and gut type (Ley et al., 2008; Zhang et al., 2010). Using culture-dependent methods, Klite (1965) concluded that food habits alone would not seem to account for the structural changes of microbial communities between bat species. A recently published mouse model reported that host diet explained 57% of total variation compared to less than 12% explained by genetic mutation (Zhang et al., 2010). Our result also indicated that the phylogenic relationship did not show clear consistency with microbiota similarity. However, if more species of the same family were included to minimize within-family variation, the influence of a phylogenic relationship might be seen as shown by a previous study (Phillips et al., 2012). Beyond host diet, other geographical and environmental factors may have more significant effects on microbiota. Further investigation may focus on more bat species of diversified feeding habits as well as their living environments to reveal the relationship of microbiota to diet, age, sex, and environmental factors. It should be noted that a bat virus study detected a considerable number of sequences belonging to bacteriophages (Ge et al., 2012); whether bacteriophages have any effect on bat microbiota is also an untouched area.
The studies on bat microbiota have become epidemiologically significant as more potential bacterial pathogens are being identified in bats, including Salmonella enteritidis (Muhldorfer et al., 2011), Yersinia enterocolitica (Muhldorfer et al., 2010), Shigella boydii (Arata et al., 1968), and Streptococcus dentirousetti (Takada and Hirasawa, 2008). Some of these species were reported to be able to cause local or systemic diseases in bats or other animals while some others did not show obvious pathogenicity and were regarded as opportunistic pathogens (Muhldorfer, 2013).
Earlier attempts to isolate bacterial pathogens usually require enormous sampling efforts and poor repeatability as pathogens were only isolated in few individuals (Hohne et al., 1975; Hazeleger et al., 2010; Muhldorfer et al., 2011). With higher resolution and better repeatability, the culture-independent technology can be applied to perform preliminary pathogen surveys. Although partial 16S rRNA genes were not able to assign sequences down to species level, the result is still meaningful as a guidance to targeted identification of pathogens. Since each pathogen carries unique DNA or RNA signatures that differentiate it from other organisms, if the database of these signatures can be curated, future metagenomic or transcriptomic studies will provide more detailed and conclusive results about bacterial pathogens (Leon et al., 2018). As an example, the identification of Bartonella mayotimonensis from metagenomic research proved the feasibility of this approach (Veikkolainen et al., 2014).
In conclusion, the pyrosequencing analysis of the commensal microbiota of four Chinese bat species revealed key features unique to bat gut flora, including the absence of Bacteroidetes group, which is a dominant group in other terrestrial mammals. Besides that, general bat microbiota is more similar to that of the Carnivora animals. Potential bacterial pathogens from 12 genera were detected, which can be further investigated.
Data Availability Statement
The datasets generated for this study can be accessed from GenBank, PRJNA576557.
Ethics Statement
This study was approved by the Animal Ethics Committee of the Wuhan Institute of Virology (Animal Ethics Approval Number: WIVA05201202).
Author Contributions
D-LS acquired and analysed data, and wrote the manuscript. Y-ZG analysed data and wrote the manuscript. X-YG and Z-LS collected samples and provided technical guidance. N-YZ contributed research ideas and modified the manuscript.
Funding
This work was funded by the National Natural Science Foundation of China (81290341).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01040/full#supplementary-material
Footnotes
References
Adesiyun, A. A., Stewart-Johnson, A., and Thompson, N. N. (2009). Isolation of enteric pathogens from bats in Trinidad. J. Wildl. Dis. 45, 952–961. doi: 10.7589/0090-3558-45.4.952
Arata, A. A., Vaughn, J. B., Newell, K. W., Barth, R. A. J., and Gracian, M. (1968). Salmonella and Shigella infections in bats in selected areas of Colombia. Am. J. Trop. Med. Hyg. 17, 92–95.
Calisher, C. H., Childs, J. E., Field, H. E., Holmes, K. V., and Schountz, T. (2006). Bats: important reservoir hosts of emerging viruses. Clin. Microbiol. Rev. 19, 531–545. doi: 10.1128/Cmr.00017-16
Canals, M., Atala, C., Olivares, R., Guajardo, F., Figueroa, D. P., Sabat, P., et al. (2005). Functional and structural optimization of the respiratory system of the bat Tadarida brasiliensis (Chiroptera, Molossidae): does airway geometry matter? J. Exp. Biol. 208, 3987–3995. doi: 10.1242/jeb.01817
Carrillo-Araujo, M., Taş, N., Alcántara-Hernández, R. J., Gaona, O., Schondube, J. E., Medellín, R. A., et al. (2015). Phyllostomid bat microbiome composition is associated to host phylogeny and feeding strategies. Front. Microbiol. 6:447. doi: 10.3389/fmicb.2015.00447
Caviedes-Vidal, E., McWhorter, T. J., Lavin, S. R., Chediack, J. G., Tracy, C. R., and Karasov, W. H. (2007). The digestive adaptation of flying vertebrates: high intestinal paracellular absorption compensates for smaller guts. Proc. Natl. Acad. Sci. U.S.A. 104, 19132–19137. doi: 10.1073/pnas.0703159104
Ciminski, K., Ran, W., Gorka, M., Lee, J., Malmlov, A., Schinkothe, J., et al. (2019). Bat influenza viruses transmit among bats but are poorly adapted to non-bat species. Nat. Microbiol. 4, 2298–2309. doi: 10.1038/s41564-019-0556-559
de Vienne, D. M., Giraud, T., and Martin, O. C. (2007). A congruence index for testing topological similarity between trees. Bioinformatics 23, 3119–3124. doi: 10.1093/bioinformatics/btm500
Di Bella, C., Piraino, C., Caracappa, S., Fornasari, L., Violani, C., and Zava, B. (2003). Enteric microflora in Italian Chiroptera. J. Mt. Ecol. 7, 221–224.
Fan, H., Walker, A. P., Carrique, L., Keown, J. R., Serna Martin, I., Karia, D., et al. (2019). Structures of influenza A virus RNA polymerase offer insight into viral genome replication. Nature 573, 287–290. doi: 10.1038/s41586-019-1530-1537
Fofanov, V. Y., and Furstenau, T. N. (2018). Guano exposed: impact of aerobic conditions on bat fecal microbiota∗. 8, 5563–5574. doi: 10.1002/ece3.4084
Gadelha-Alves, R., Rozensztranch, A. M. D., and Rocha-Barbosa, O. (2008). Comparative intestinal histomorphology of five species of Phyllostomid Bats (Phyllostomidae, Microchiroptera): ecomorphological relations with alimentary habits. Int. J. Morphol. 26, 591–602.
Ge, X. Y., Li, J. L., Yang, X. L., Chmura, A. A., Zhu, G., Epstein, J. H., et al. (2013). Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature 503:12711. doi: 10.1038/nature12711
Ge, X. Y., Li, Y., Yang, X. L., Zhang, H. J., Zhou, P., Zhang, Y. Z., et al. (2012). Metagenomic analysis of viruses from bat fecal samples reveals many novel viruses in insectivorous bats in China. J. Virol. 86, 4620–4630. doi: 10.1128/Jvi.06671-6611
Guillén-Servent, A., and Francis, C. M. (2006). A new species of bat of the Hipposideros bicolor group (Chiroptera: Hipposideridae) from Central Laos, with evidence of convergent evolution with Sundaic taxa. Acta Chiropterol. 8, 39–61. doi: 10.3161/1733-5329(2006)8[39:ansobo]2.0.co;2
Hajkova, P., and Pikula, J. (2007). Veterinary treatment of evening bats (Vespertilionidae) in the Czech Republic. Vet. Rec. 161, 139–140. doi: 10.1136/vr.161.4.139
Halpin, K., Young, P. L., Field, H. E., and Mackenzie, J. S. (2000). Isolation of Hendra virus from pteropid bats: a natural reservoir of Hendra virus. J. Gen. Virol. 81, 1927–1932. doi: 10.1099/0022-1317-81-8-1927
Hayman, D. T. S. (2019). Bat tolerance to viral infections. Nat. Microbiol. 4, 728–729. doi: 10.1038/s41564-019-0430-439
Hazeleger, W. C., Jacobs-Reitsma, W. F., van Bergen, M., Lina, P. H., van den Bosch, T., and Beumer, R. R. (2010). “Bats may be carriers of Campylobacter,” in Proceedings of the 15th Bat Research Conference Abstract book, Prague, 72.
Heard, D. J., De Young, J. L., Goodyear, B., and Ellis, G. A. (1997). Comparative rectal bacterial flora of four species of flying fox (Pteropus sp.). J. Zoo Wildl. Med. 28, 471–475.
Helmick, K. E., Heard, D. J., Richey, L., Finnegan, M., Ellis, G. A., Nguyen, A., et al. (2004). A pasteurella-like bacterium associated with pneumonia in captive megachiropterans. J. Zoo Wildl. Med. 35, 88–93. doi: 10.1638/01-083
Hohne, K., Loose, B., and Seeliger, H. P. R. (1975). Isolation of Listeria-monocytogenes in slaughter animals and bats of togo (West-Africa). Ann. Microbiol. A 126, 501–507.
Klite, P. D. (1965). Intestinal bacterial flora and transit time of 3 neotropical bat species. J. Bacteriol. 90, 375–379. doi: 10.1128/jb.90.2.375-379.1965
Laing, E. D., Navaratnarajah, C. K., Cheliout Da Silva, S., Petzing, S. R., Xu, Y., Sterling, S. L., et al. (2019). Structural and functional analyses reveal promiscuous and species specific use of ephrin receptors by Cedar virus. Proc. Natl. Acad. Sci. U.S.A. 116, 20707–20715. doi: 10.1073/pnas.1911773116
Lee, O. O., Wang, Y., Yang, J. K., Lafi, F. F., Al-Suwailem, A., and Qian, P. Y. (2011). Pyrosequencing reveals highly diverse and species-specific microbial communities in sponges from the Red Sea. ISME J. 5, 650–664. doi: 10.1038/ismej.2010.165
Leon, M. P. D., Montecillo, A. D., Pinili, D. S., Siringan, M. A. T., and Park, D. S. (2018). Bacterial diversity of bat guano from Cabalyorisa Cave, Mabini, Pangasinan, Philippines: a first report on the metagenome of Philippine bat guano. PLoS One 13:e0200095. doi: 10.1371/journal.pone.0200095
Leroy, E. M., Kumulungui, B., Pourrut, X., Rouquet, P., Hassanin, A., Yaba, P., et al. (2005). Fruit bats as reservoirs of Ebola virus. Nature 438, 575–576. doi: 10.1038/438575a
Ley, R. E., Hamady, M., Lozupone, C., Turnbaugh, P. J., Ramey, R. R., Bircher, J. S., et al. (2008). Evolution of mammals and their gut microbes. Science 320, 1647–1651. doi: 10.1126/science.1155725
Li, J., Li, L., Jiang, H., Yuan, L., Zhang, L., Ma, J. E., et al. (2018). Fecal bacteriome and mycobiome in bats with diverse diets in South China. Curr. Microbiol. 75, 1352–1361. doi: 10.1007/s00284-018-1530-1530
Li, W. D., Shi, Z. L., Yu, M., Ren, W. Z., Smith, C., Epstein, J. H., et al. (2005). Bats are natural reservoirs of SARS-like coronaviruses. Science 310, 676–679. doi: 10.1126/science.1118391
Lozupone, C., Hamady, M., and Knight, R. (2006). UniFrac - An online tool for comparing microbial community diversity in a phylogenetic context. BMC Bioinformatics 7:371. doi: 10.1186/1471-2105-7-371
Makanya, A. N., and Maina, J. N. (1994). The morphology of the intestine of the insectivorous horseshoe bat (Rhinolophus-Hildebrandti, Peters) - a scanning electron and light-microscopic study. Afr. J. Ecol. 32, 158–168. doi: 10.1111/j.1365-2028.1994.tb00566.x
Matthias, M. A., Diaz, M. M., Campos, K. J., Calderon, M., Willig, M. R., Pacheco, V., et al. (2005). Diversity of bat-associated Leptospira in the Peruvian Amazon inferred by Bayesian phylogenetic analysis of 16S ribosomal DNA sequences. Am. J. Trop. Med. Hyg. 73, 964–974. doi: 10.4269/ajtmh.2005.73.964
McKenna, P., Hoffmann, C., Minkah, N., Aye, P. P., Lackner, A., Liu, Z., et al. (2008). The macaque gut microbiome in health, lentiviral infection, and chronic enterocolitis. PLoS Pathog. 4:e20. doi: 10.1371/journal.ppat.0040020
Muhldorfer, K. (2013). Bats and bacterial pathogens: a review. Zoonoses Public Health 60, 93–103. doi: 10.1111/j.1863-2378.2012.01536.x
Muhldorfer, K., Speck, S., Kurth, A., Lesnik, R., Freuling, C., Muller, T., et al. (2011). Diseases and causes of death in European bats: dynamics in disease susceptibility and infection rates. PLoS One 6:e29773. doi: 10.1371/journal.pone.0029773
Muhldorfer, K., Wibbelt, G., Haensel, J., Riehm, J., and Speck, S. (2010). Yersinia species isolated from bats, Germany. Emerging Infect. Dis. 16, 578–580. doi: 10.3201/eid1603.091035
Murakoshi, F., Koyama, K., Akasaka, T., Horiuchi, N., and Kato, K. (2018). Molecular and histopathological characterization of Cryptosporidium and Eimeria species in bats in Japan. J. Vet. Med. Sci. 80, 1395–1399. doi: 10.1292/jvms.18-0130
Muyzer, G., Dewaal, E. C., and Uitterlinden, A. G. (1993). Profiling of complex microbial-populations by Denaturing Gradient Gel-Electrophoresis analysis of polymerase chain reaction-amplified genes-coding for 16S Ribosomal-RNA. Appl. Environ. Microbiol. 59, 695–700. doi: 10.1128/aem.59.3.695-700.1993
Nam, Y. D., Jung, M. J., Roh, S. W., Kim, M. S., and Bae, J. W. (2011). Comparative analysis of Korean human gut microbiota by barcoded pyrosequencing. PLoS One 6:e22109. doi: 10.1371/journal.pone.0022109
Nava, G. M., Friedrichsen, H. J., and Stappenbeck, T. S. (2011). Spatial organization of intestinal microbiota in the mouse ascending colon. ISME J. 5, 627–638. doi: 10.1038/ismej.2010.161
Phillips, C. D., Phelan, G., Dowd, S. E., McDonough, M. M., Ferguson, A. W., Delton Hanson, J., et al. (2012). Microbiome analysis among bats describes influences of host phylogeny, life history, physiology and geography. Mol. Ecol. 21, 2617–2627. doi: 10.1111/j.1365-294X.2012.05568.x
Quince, C., Lanzen, A., Curtis, T. P., Davenport, R. J., Hall, N., Head, I. M., et al. (2009). Accurate determination of microbial diversity from 454 pyrosequencing data. Nat. Methods 6, 639–641. doi: 10.1038/nmeth.1361
Quince, C., Lanzen, A., Davenport, R. J., and Turnbaugh, P. J. (2011). Removing noise from pyrosequenced amplicons. BMC Bioinformatics 12:38. doi: 10.1186/1471-2105-12-38
Reyes, A. W. B., Rovira, H. G., Masangkay, J. S., Ramirez, T. J., Yoshikawa, Y., and Baticados, W. N. (2011). Polymerase chain reaction assay and conventional isolation of Salmonella spp. from Philippine bats∗. Acta Sci. Vet. 39.
Roesch, L. F., Fulthorpe, R. R., Riva, A., Casella, G., Hadwin, A. K. M., Kent, A. D., et al. (2007). Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 1, 283–290. doi: 10.1038/ismej.2007.53
Schloss, P. D., Gevers, D., and Westcott, S. L. (2011). Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS One 6:e27310. doi: 10.1371/journal.pone.0027310
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/Aem.01541-1549
Shanks, O. C., Kelty, C. A., Archibeque, S., Jenkins, M., Newton, R. J., McLellan, S. L., et al. (2011). Community structures of fecal bacteria in cattle from different animal feeding operations. Appl. Environ. Microbiol. 77, 2992–3001. doi: 10.1128/Aem.02988-2910
Shen, Y. Y., Liang, L., Zhu, Z. H., Zhou, W. P., Irwin, D. M., and Zhang, Y. P. (2010). Adaptive evolution of energy metabolism genes and the origin of flight in bats. Proc. Natl. Acad. Sci. U.S.A. 107, 8666–8671. doi: 10.1073/pnas.0912613107
Springer, M. S., Teeling, E. C., Madsen, O., Stanhope, M. J., and de Jong, W. W. (2001). Integrated fossil and molecular data reconstruct bat echolocation. Proc. Natl. Acad. Sci. U.S.A. 98, 6241–6246. doi: 10.1073/pnas.111551998
Suchodolski, J. S., Camacho, J., and Steiner, J. M. (2008). Analysis of bacterial diversity in the canine duodenum, jejunum, ileum, and colon by comparative 16S rRNA gene analysis. FEMS Microbiol. Ecol. 66, 567–578. doi: 10.1111/j.1574-6941.2008.00521.x
Sun, D. L., Jiang, X., Wu, Q. L., and Zhou, N. Y. (2013). Intragenomic heterogeneity in 16S rRNA genes causes overestimation of prokaryotic diversity. Appl. Environ. Microbiol. 79:AEM.01282-13. doi: 10.1128/AEM.01282-1213
Takada, K., and Hirasawa, M. N. (2008). Streptococcus dentirousetti sp nov., isolated from the oral cavities of bats. Int. J. Syst. Evol. Micr. 58, 160–163. doi: 10.1099/ijs.0.65204-65200
Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729. doi: 10.1093/molbev/mst197
Teeling, E. C., Madsen, O., Van den Bussche, R. A., de Jong, W. W., Stanhope, M. J., and Springer, M. S. (2002). Microbat paraphyly and the convergent evolution of a key innovation in Old World rhinolophoid microbats. Proc. Natl. Acad. Sci. U.S.A. 99, 1431–1436. doi: 10.1073/pnas.022477199
Teeling, E. C., Springer, M. S., Madsen, O., Bates, P., O’Brien, S. J., and Murphy, W. J. (2005). A molecular phylogeny for bats illuminates biogeography and the fossil record. Science 307, 580–584. doi: 10.1126/science.1105113
Veikkolainen, V., Vesterinen, E. J., Lilley, T. M., and Pulliainen, A. T. (2014). Bats as reservoir hosts of human bacterial pathogen, Bartonella mayotimonensis. Emerging Infect. Dis. 20, 960–967. doi: 10.3201/eid2006.130956
Vengust, M., Knapic, T., and Weese, J. S. (2018). The fecal bacterial microbiota of bats; Slovenia. PLoS One 13:e0196728. doi: 10.1371/journal.pone.0196728
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-67
Yu, Z. T., and Morrison, M. (2004). Improved extraction of PCR-quality community DNA from digesta and fecal samples. Biotechniques 36, 808–812. doi: 10.2144/04365st04
Zhang, C., Zhang, M., Wang, S., Han, R., Cao, Y., Hua, W., et al. (2010). Interactions between gut microbiota, host genetics and diet relevant to development of metabolic syndromes in mice. ISME J. 4, 232–241. doi: 10.1038/ismej.2009.112
Zhang, H., Gao, S., Lercher, M. J., Hu, S., and Chen, W. H. (2012). EvolView, an online tool for visualizing, annotating and managing phylogenetic trees. Nucleic Acids Res. 40, W569–W572. doi: 10.1093/nar/gks576
Keywords: bat, bacterial pathogen, gut microbiota, aerobic microbes, sequencing
Citation: Sun D-L, Gao Y-Z, Ge X-Y, Shi Z-L and Zhou N-Y (2020) Special Features of Bat Microbiota Differ From Those of Terrestrial Mammals. Front. Microbiol. 11:1040. doi: 10.3389/fmicb.2020.01040
Received: 13 October 2019; Accepted: 27 April 2020;
Published: 03 June 2020.
Edited by:
Spyridon Ntougias, Democritus University of Thrace, GreeceReviewed by:
Michelle Baker, Australian Animal Health Laboratory (CSIRO), AustraliaKonstantinos Papadimitriou, University of Peloponnese, Greece
Copyright © 2020 Sun, Gao, Ge, Shi and Zhou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ning-Yi Zhou, bmluZ3lpLnpob3VAc2p0dS5lZHUuY24=
†These authors have contributed equally to this work