<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2023.1207482</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>The difference in the composition of gut microbiota is greater among bats of different phylogenies than among those with different dietary habits</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Guo</surname>
<given-names>Min</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/611056/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xie</surname>
<given-names>Siwei</given-names>
</name>
<xref rid="aff2" ref-type="aff"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Junhua</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Yuzhi</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>He</surname>
<given-names>Xiangyang</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Luo</surname>
<given-names>Pengfei</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Deng</surname>
<given-names>Jin</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
<xref rid="aff3" ref-type="aff"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhou</surname>
<given-names>Chunhui</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Qin</surname>
<given-names>Jiao</given-names>
</name>
<xref rid="aff4" ref-type="aff"><sup>4</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/640171/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Huang</surname>
<given-names>Chen</given-names>
</name>
<xref rid="aff5" ref-type="aff"><sup>5</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1410926/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhang</surname>
<given-names>Libiao</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
<xref rid="c001" ref-type="corresp"><sup>&#x002A;</sup></xref>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Guangdong Key Laboratory of Animal Conservation and Resource Utilization, Guangdong Public Laboratory of Wild Animal Conservation and Utilization, Institute of Zoology, Guangdong Academy of Sciences</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>College of Mathematics and Informatics, South China Agricultural University</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>College of Life Sciences, South China Normal University</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country></aff>
<aff id="aff4"><sup>4</sup><institution>Institute of Plant Protection, Guangdong Academy of Agricultural Sciences</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country></aff>
<aff id="aff5"><sup>5</sup><institution>Dr. Neher&#x2019;s Biophysics Laboratory for Innovative Drug Discovery, State Key Laboratory of Quality Research in Chinese Medicine, Macau University of Science and Technology</institution>, <addr-line>Taipa</addr-line>, <country>Macao SAR, China</country></aff>
<author-notes>
<fn fn-type="edited-by" id="fn0006">
<p>Edited by: Decai Jin, Chinese Academy of Sciences (CAS), China</p>
</fn>
<fn fn-type="edited-by" id="fn0007">
<p>Reviewed by: Anbu Poosakkannu, University of Turku, Finland; Sunil Banskar, University of Arizona, United States</p>
</fn>
<corresp id="c001">&#x002A;Correspondence: Libiao Zhang, <email>zhanglb@giz.gd.cn</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>28</day>
<month>07</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1207482</elocation-id>
<history>
<date date-type="received">
<day>17</day>
<month>04</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>10</day>
<month>07</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2023 Guo, Xie, Wang, Zhang, He, Luo, Deng, Zhou, Qin, Huang and Zhang.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Guo, Xie, Wang, Zhang, He, Luo, Deng, Zhou, Qin, Huang and Zhang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<p>Bats have a very long evolutionary history and are highly differentiated in their physiological functions. Results of recent studies suggest effects of some host factors (e.g., phylogeny and dietary habit) on their gut microbiota. In this study, we examined the gut microbial compositions of 18 different species of bats. Results showed that Firmicutes, Gammaproteobacteria, and Actinobacteria were dominant in all fecal samples of bats. However, the difference in the diversity of gut microbiota among bats of different phylogenies was notable (<italic>p</italic>&#x2009;=&#x2009;0.06). Various species of Firmicutes, Actinobacteria, and Gammaproteobacteria were found to contribute to the majority of variations in gut microbiota of all bats examined, and <italic>Aeromonas</italic> species were much more abundant in bats that feed on both insects and fish than in those of insectivores. The abundance of various species of <italic>Clostridium</italic>, <italic>Euryarchaeota</italic>, and ancient bacterial phyla was found to vary among bats of different phylogenies, and various species of <italic>Vibrio</italic> varied significantly among bats with different dietary habits. No significant difference in the number of genes involved in various metabolic pathways was detected among bats of different phylogenies, but the abundance of genes involved in 5 metabolic pathways, including transcription; replication, recombination, and repair; amino acid transport and metabolism; and signal transduction mechanisms, was different among bats with different dietary habits. The abundance of genes in 3 metabolic pathways, including those involved in stilbenoid, diarylheptanoid, and gingerol biosynthesis, was found to be different between insectivorous bats and bats that feed on both insects and fish. Results of this study suggest a weak association between dietary habit and gut microbiota in most bats but a notable difference in gut microbiota among bats of different phylogenies.</p>
</abstract>
<kwd-group>
<kwd>bat</kwd>
<kwd>phylogeny</kwd>
<kwd>dietary habit</kwd>
<kwd>gut microbiota</kwd>
<kwd>dietary shift</kwd>
</kwd-group>
<contract-num rid="cn1">2021GDASYL-20210103052</contract-num>
<contract-num rid="cn2">2021B1212110003</contract-num>
<contract-num rid="cn2">2021B1212050021</contract-num>
<contract-num rid="cn3">32172435</contract-num>
<contract-sponsor id="cn1">GDAS Special Project of Science and Technology Development</contract-sponsor>
<contract-sponsor id="cn2">Guangdong Provincial Science and Technology Program</contract-sponsor>
<contract-sponsor id="cn3">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content></contract-sponsor>
<counts>
<fig-count count="4"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="44"/>
<page-count count="9"/>
<word-count count="5706"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Microbial Symbioses</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="sec1">
<title>Introduction</title>
<p>Gut microbiota is composed of many bacteria, fungi, viruses, and other organisms in the gastrointestinal system. The composition of gut microbiota in animals varies greatly in their early life but becomes stable in adulthood (<xref ref-type="bibr" rid="ref27">Lim et al., 2015</xref>; <xref ref-type="bibr" rid="ref41">Yin et al., 2020</xref>). Many large-scale studies have revealed effects of gut microbiota on the well-being of the host (<xref ref-type="bibr" rid="ref38">Wan et al., 2020</xref>; <xref ref-type="bibr" rid="ref21">Kuang et al., 2022</xref>). Host factors such as phylogeny, physiology, behavior, and dietary habit have also been shown to affect the composition of gut microbiota (<xref ref-type="bibr" rid="ref4">Benson et al., 2010</xref>; <xref ref-type="bibr" rid="ref23">Kurilshikov et al., 2021</xref>).</p>
<p>A recent large-scale study of the gut bacteria of bats, amphibians, and birds showed that bats and birds share similar profiles of gut bacteria (<xref ref-type="bibr" rid="ref35">Song et al., 2020</xref>), suggesting that an innate ability, such as flying, may have similar effects on gut bacteria in animals of different phylogenies. In contrast, other studies disclosed a link between specific microorganisms and their hosts (<xref ref-type="bibr" rid="ref8">Brucker and Bordenstein, 2012</xref>; <xref ref-type="bibr" rid="ref14">Easson and Thacker, 2014</xref>; <xref ref-type="bibr" rid="ref32">Sanders et al., 2014</xref>; <xref ref-type="bibr" rid="ref6">Brooks et al., 2016</xref>). Such relationship is referred to as phylosymbiosis that predicts host clades harbor distinguishable microbial communities (<xref ref-type="bibr" rid="ref8">Brucker and Bordenstein, 2012</xref>) and is common in some animals, such as sponges (<xref ref-type="bibr" rid="ref14">Easson and Thacker, 2014</xref>), Nasonia (<xref ref-type="bibr" rid="ref8">Brucker and Bordenstein, 2012</xref>), ants, and apes (<xref ref-type="bibr" rid="ref32">Sanders et al., 2014</xref>). These findings suggest a relationship between the composition of gut microbiota and host evolutionary history, behaviors, and ecological factors (<xref ref-type="bibr" rid="ref42">Zepeda Mendoza et al., 2018</xref>; <xref ref-type="bibr" rid="ref28">Lutz et al., 2019</xref>).</p>
<p>Bat is the second largest order of mammals and have a long evolutionary history, distinct physiology, diverse dietary habits, echolocations, and flying ability (<xref ref-type="bibr" rid="ref1">Adams and Pedersen, 2013</xref>). Similar to carnivores, the gut microbiota of bats are dominated by Firmicutes and Proteobacteria (<xref ref-type="bibr" rid="ref36">Sun et al., 2020</xref>). It has been shown that the compositions of gut bacteria of bats vary with their diet, sex, age, mass, geographic location, physiology, reproduction, and symbiotic parasites (<xref ref-type="bibr" rid="ref29">Phillips et al., 2012</xref>; <xref ref-type="bibr" rid="ref10">Carrillo-Araujo et al., 2015</xref>; <xref ref-type="bibr" rid="ref28">Lutz et al., 2019</xref>). Although the effects of host factors on gut bacteria have been described (<xref ref-type="bibr" rid="ref29">Phillips et al., 2012</xref>; <xref ref-type="bibr" rid="ref10">Carrillo-Araujo et al., 2015</xref>; <xref ref-type="bibr" rid="ref28">Lutz et al., 2019</xref>), the mechanisms of such effects are unknown. We hypothesized that host phylogeny and dietary habits could affect the diversity of gut microbiota and therefore conducted this study to compare the compositions of gut microbiota in bats of different phylogenies and those of insect-eating (insectivores) bats, fruit-eating (frugivores) bat, and bats that feed on both insects and fish, aiming to uncover the impact of various host factors on the gut microbiota of bats.</p>
</sec>
<sec sec-type="materials|methods" id="sec2">
<title>Materials and methods</title>
<sec id="sec3">
<title>Samples</title>
<p>A total of 24 bats of 18 different species belonging to 8 genera and 3 families (Pteropodidae, Rhinolophidae, and Vespertilionidae) were captured from Yanyan Cave, Huizhou City, China. These bats included 6 <italic>Monascus pilosus</italic> bats, 2 <italic>Monascus chinensis</italic> bats, and one each of the following: <italic>Monascus altarium</italic>, <italic>Monascus davidii</italic>, <italic>Monascus longipes</italic>, <italic>Rhinolophus pearsonii</italic>, <italic>Rhinolophus pusillus</italic>, <italic>Rhinolophus siamensis</italic>, <italic>Rhinolophus affinis</italic>, <italic>Rhinolophus macrotis</italic>, <italic>Hypsugo pulveratus</italic>, <italic>Eptesicus serotinus</italic>, <italic>Pipistrellus abramus</italic>, <italic>Murina huttoni</italic>, <italic>Murina aurata</italic>, <italic>Miniopterus pusillus</italic>, <italic>Macropus fuliginosus</italic>, <italic>Rousettus leschenaultii</italic> (<xref ref-type="bibr" rid="ref40">Wilson and Mittermeier, 2019</xref>). These bat species were selected because they have different dietary habits in addition to different phylogenies. Among the 18 bat species, <italic>M. pilosus</italic> bats feed on both insects and fish, <italic>R. leschenaultii</italic> bat is frugivores, and the remaining bats are insectivores (<xref ref-type="bibr" rid="ref40">Wilson and Mittermeier, 2019</xref>).</p>
<p>To collect fecal samples, each bat was kept in a clean cloth bag for about 5&#x2009;h. Fecal particles were then collected and placed in 2&#x2009;mL tubes. For those with no fecal particles, anal swabs were used. A total of 24 fecal samples were collected (one sample from each bat). The samples were immediately snap frozen in liquid nitrogen and then stored in a &#x2212;80&#x00B0;C freezer until used.</p>
</sec>
<sec id="sec4">
<title>DNA extraction, metagenomic assembly, and taxonomical annotation</title>
<p>DNA was extracted from fecal samples using the QIAamp<sup>&#x00AE;</sup> Fast DNA Stool Mini kit following manufacturer&#x2019;s instructions without any modifications. DNA samples were then subjected metagenomic analysis by next generation sequencing (NGS) with a NovaSeq 6000 sequencer (read length: 150 bases; sequencing depth: 10 Gb; paired-end sequencing by Shanghai Origingene Corporation, China).</p>
<p>To filter out low-quality and adapter contaminated reads, the raw reads were subjected to quality control by SOAPnuke (1.5.6; <xref ref-type="bibr" rid="ref12">Chen et al., 2018</xref>). The high-quality reads were assembled by Megahit<xref rid="fn0001" ref-type="fn"><sup>1</sup></xref> (<xref ref-type="bibr" rid="ref26">Li et al., 2015</xref>). MetaGeneMark<xref rid="fn0002" ref-type="fn"><sup>2</sup></xref> (<xref ref-type="bibr" rid="ref44">Zhu et al., 2010</xref>) was employed to search for coding sequences (CDSs) in assembled contigs longer than 100&#x2009;bp. To construct a non-redundant gene set, genes of all samples were merged and clustered with 95% identity by Cd-hit (v4.6.4)<xref rid="fn0003" ref-type="fn"><sup>3</sup></xref> (<xref ref-type="bibr" rid="ref25">Li and Godzik, 2006</xref>). Taxonomic identification of genes in the non-redundant gene set was performed with Blastp (v2.2.28)<xref rid="fn0004" ref-type="fn"><sup>4</sup></xref> (<xref ref-type="bibr" rid="ref9">Camacho et al., 2009</xref>) by aligning the genes against those in the Nucleotide Sequence Database (Nr database; v20170519) at an <italic>e</italic>-value &#x2264;1e&#x2212;5. When the gene length coverage was less than 50%, the gene was assigned an &#x201C;Unclassified&#x201D; status. Genes were annotated by aligning high quality reads to the non-redundant gene set using Bowtie2 (v2.2.5; <xref ref-type="bibr" rid="ref24">Langmead and Salzberg, 2012</xref>), and gene abundance was calculated by PathoScope2 (v2.0.6; <xref ref-type="bibr" rid="ref16">Hong et al., 2014</xref>). Species abundance was determined by the sum of assigned gene counts. The taxonomic tree of microbes was built according to the NCBI taxonomy. Divergence time of bat species in the phylogenetic tree was determined by TimeTree 5<xref rid="fn0005" ref-type="fn"><sup>5</sup></xref> (<xref ref-type="bibr" rid="ref22">Kumar et al., 2022</xref>).</p>
</sec>
<sec id="sec5">
<title>Gene function and metabolic pathway annotations</title>
<p>The non-redundant gene set was annotated by blasting translated proteins against the database of Clusters of Orthologous Groups of proteins (COG) (v20090331) and that of Kyoto Encyclopedia of Genes and Genomes (KEGG) (v81). COG-annotated gene functions were visualized using the boxplot function of the Lattice package in the R software. Significantly changed gene functions and pathways between groups were detected by the Wilcoxon test of the R software (adjusted <italic>p</italic>&#x2009;&#x003C;&#x2009;0.05, by the method of Holm).</p>
</sec>
<sec id="sec6">
<title>Evaluation of microbial diversity</title>
<p>To assess the compositional homogeneity of the microbial communities in samples, Shannon&#x2013;Wiener index (<xref ref-type="bibr" rid="ref34">Shannon, 2001</xref>) was calculated using vegan, ggplot2, RColorBrewer, and reshape2 packages of the R software. Weighted principal component analysis (PCA) was executed using ggbiplot, ggrepel, and ggplot2 packages. <italic>Z</italic>-score normalization (<xref ref-type="bibr" rid="ref33">Shalabi et al., 2006</xref>) was performed with the prcomp package prior to calculation of the principal components. The ggbiplot package was used to visualize the distribution of each sample in PC1 and PC2. Significant difference between groups was determined using the adonis package of the R software (permutations&#x2009;=&#x2009;999, by the method of Euclidean; <xref ref-type="bibr" rid="ref2">Anderson, 2001</xref>).</p>
</sec>
<sec id="sec7">
<title>Identification of significant gut microorganisms and host metabolic pathways</title>
<p>To identify the significant gut microorganisms, 500 tree models were constructed by random forest (RF) modeling. Each tree model was trained by bootstrapping samples from all samples; those samples unemployed in the training (~36.8%) were used to test the accuracy of the tree model (<xref ref-type="bibr" rid="ref5">Breiman, 2001</xref>). RF models were built by the package randomForest of the R software (with parameters: ntree&#x2009;=&#x2009;500, mtry&#x2009;=&#x2009;sqrt(<italic>p</italic>) where <italic>p</italic> is number of variables, replace&#x2009;=&#x2009;TRUE, maxnodes&#x2009;=&#x2009;NULL, oob.prox&#x2009;=&#x2009;proximity, nPerm&#x2009;=&#x2009;1, maxnodes&#x2009;=&#x2009;NULL) to analyze data of bats with different dietary habits and those of different phylogenies by random permutation of species-level abundance profiles. The out-of-bag (OOB) error rate was used to assess the accuracy of modeling. Microbial influence value was determined by the mean decrease in the Gini index of each node of the constructed model, and a higher value indicated a more significant effect on the composition of gut microbiota. The number of microbial species was determined by repeating RF modeling 5 times with 10-fold cross validation.</p>
<p>The significantly changed (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.01) abundance of bacterial species between insectivorous bats and bats that eat both insects and fish were determined by Metastats non-parametric multiple tests with <italic>p</italic>-values adjusted for false discovery rate (FDR) to aviod errors (<xref ref-type="bibr" rid="ref39">White et al., 2009</xref>). The significant difference in the number of genes involved in certain functions or pathways between groups were detected by the Wilcoxon test (adjusted <italic>p</italic>&#x2009;&#x003C;&#x2009;0.05 by the method of Holm). The abundance profiles of various significant microbes were visualized in heat map by packages gplots, fields, and akima of the R software.</p>
</sec>
</sec>
<sec sec-type="results" id="sec8">
<title>Results</title>
<sec id="sec9">
<title>Sequencing data</title>
<p>The NGS performed with the NovaSeq 6000 platform on the 24 bat fecal samples generated approximately 7.9&#x2009;GB raw data per sample. After cleaning, an average of 7.8 Gb clean data were obtained from each sample with 22,213&#x2013;1,050,926 contigs and 25,969&#x2013;688,291 genes. The average N50 of all genes was 570&#x2009;bp (<xref rid="fig1" ref-type="fig">Figure 1</xref>), and 7,525,231 unique genes were obtained from all samples.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p>Sequencing data. <bold>(A)</bold> Length distribution of raw data. <bold>(B)</bold> Number distribution of contigs and genes. <bold>(C)</bold> Length distribution of N50 contigs and genes.</p>
</caption>
<graphic xlink:href="fmicb-14-1207482-g001.tif"/>
</fig>
</sec>
<sec id="sec10">
<title>Variations in gut microbiota among bats of different phylogenies or with different diets</title>
<p>After taxonomic annotation, 394 archaea species, 18,271 bacterial species, 757 fungal species, and 1,912 viral species were detected. Except for <italic>P. abramus</italic>, all samples showed similar microbial profiles at the phylum level (<xref rid="fig2" ref-type="fig">Figure 2A</xref>) with Proteobacteria being most predominant (average relative abundance: 0.658), followed by Firmicutes (0.204), norank_d__Viruses (0.054), Chlamydiae (0.018), and Deinococcus-Thermus (0.010). The predominant bacterial genera were <italic>Enterobacter</italic> (0.054), <italic>Escherichia</italic> (0.047), norank_f__Enterobacteriaceae (0.041), <italic>Serratia</italic> (0.040), <italic>Klebsiella</italic> (0.032), <italic>Bacillus</italic> (0.031), <italic>Morganella</italic> (0.030), <italic>Lactococcus</italic> (0.022), <italic>Citrobacter</italic> (0.021), <italic>Aeromonas</italic> (0.019), <italic>Kluyvera</italic> (0.017), and <italic>Enterococcus</italic> (0.015) (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S1</xref>).</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>Diversity of gut microbiota in various species of bats. <bold>(A)</bold> Top 18 microbial phyla in gut microbiota of bats. The bar below the phylogenetic tree indicates divergence time of the bat species. <bold>(B)</bold> Diversity of gut microbiota in bats with different diet habits (the two upper panels) and those of different phylogenies (at genus level, the two lower panels).</p>
</caption>
<graphic xlink:href="fmicb-14-1207482-g002.tif"/>
</fig>
<p>There was no significant difference in microbial diversity was detected among bats with insect (IB) or both insect and fish (IFB) diets by statistical (Wilcoxon test adjusted by the method of Holm, all <italic>p</italic>&#x2009;&#x003E;&#x2009;0.05) and PCA analyses (<italic>p</italic>&#x2009;=&#x2009;0.69; <xref rid="fig2" ref-type="fig">Figure 2B</xref>). Meanwhile, in PCA analysis, there was no significant difference at the family level of microbes among bats of different phylogenies (<italic>R</italic><sup>2</sup>&#x2009;=&#x2009;0.049, <italic>p</italic>&#x2009;=&#x2009;0.743) (<xref ref-type="supplementary-material" rid="SM2">Supplementary Figure S2</xref>), but a high variation in microbial diversity at the genus level was observed (<italic>R</italic><sup>2</sup>&#x2009;=&#x2009;0.54, <italic>p</italic>&#x2009;=&#x2009;0.06; <xref rid="fig2" ref-type="fig">Figure 2B</xref>).</p>
</sec>
<sec id="sec11">
<title>Identification of microbes responsible for major variations in gut microbiota and major metabolic pathways in bats</title>
<p>The random forest (RF) method was used to identify key microbes responsible for major variations in the composition of gut microbiota, and 53 different species of microbes were found among bats of different phylogenies (PG bats; OOB error rate: 41.67%). For bats with different dietary habits (DG bats), 38 different species of microbes were found (OOB error rate: 16.67%) (<xref rid="fig3" ref-type="fig">Figure 3A</xref>). By constructing a phylogenetic tree of all 91 microbes, species of Firmicutes, Actinobacteria, and Gammaproteobacteria were found to contribute the most variations in gut microbiota in all bats. For bats of different phylogenies, species of <italic>Clostridium</italic>, <italic>Euryarchaeota</italic>, and ancient bacterial phyla contributed the major variation. In contrast, species of <italic>Vibrio</italic> contributed almost solely the major variation in gut microbiota in bats with different dietary habits (<xref rid="fig3" ref-type="fig">Figure 3B</xref>).</p>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption>
<p>Major microbial species in gut microbiota of bats. <bold>(A)</bold> Gray lines show results of 5 repeats of RF modeling with 10-fold cross validations. Dark lines indicate average values of the 5 repeats. Dashed lines denote the lowest point of the dark line, where the model had the lowest out-of-bag (OOB) error. The bar length in MeanDecreaseGini shows the extent of microbes affecting node impurity in RF modeling, and a microbe with a longer bar has a higher effect on the composition of gut microbiota. <bold>(B)</bold> Phylogenetic tree of the 91 most significant microbes among bats. Red branches represent the 53 most significant microbes among bats of different phylogenies, and green branches denote the 38 most significant microbes among bats with different dietary habits. Scale bar length indicates the number of amino acid substitutions per site.</p>
</caption>
<graphic xlink:href="fmicb-14-1207482-g003.tif"/>
</fig>
<p>For all bats, the predominant orthologous gene functions were Transcription, Replication, recombination, and repair, Amino acid transport and metabolism, and Signal transduction mechanisms (<xref ref-type="supplementary-material" rid="SM3">Supplementary Figure S3</xref>). The predominant metabolic pathways were Transporters for insectivorous bats (IB), Transporters and Carbon metabolism for bats that feed on both insects and fish (IFB), and Transporters and Biosynthesis of amino acids for frugivorous bat (FB) (<xref ref-type="supplementary-material" rid="SM4">Supplementary Figure S4</xref>). No significant difference in the number of these orthologous genes or metabolic pathways was detected among all bats. However, the number of genes of metabolic pathways was found to vary among bats with different dietary habits as follows: Photosynthesis - antenna proteins between IB and IFB (<italic>p</italic>&#x2009;=&#x2009;0.036); Flavonoid biosynthesis between IB and FB (<italic>p</italic>&#x2009;=&#x2009;0.0002) and between IFB and FB (<italic>p</italic>&#x2009;=&#x2009;0.041); Isoflavonoid biosynthesis between IB and IFB (<italic>p</italic>&#x2009;=&#x2009;0.035); Flavone and flavonol biosynthesis between IB and IFB (<italic>p</italic>&#x2009;=&#x2009;0.032); and Stilbenoid, diarylheptanoid, and gingerol biosynthesis between IB and FB (<italic>p</italic>&#x2009;=&#x2009;0.0002) and between IFB and FB (<italic>p</italic>&#x2009;=&#x2009;0.041).</p>
</sec>
<sec id="sec12">
<title>Significantly changed microbes during dietary transition</title>
<p>By performing Metastats non-parametric multiple tests, the abundance of 62 bacterial species was found to vary (<xref rid="fig4" ref-type="fig">Figure 4A</xref>). Most of them were also observed in the RF modeling of bats with different dietary habits, e.g., members of <italic>Vibrio</italic>, <italic>Aeromonas</italic>, <italic>Streptomyces</italic>, <italic>Cytophagaceae</italic>, <italic>Pseudoalteromonas</italic>, <italic>Alishewanella</italic>, <italic>Pseudoalteromonas</italic>, <italic>Marinobacter</italic>, <italic>Neisseria</italic>, and <italic>Gammaproteobacteria bacterium</italic> TMED95. A large number of <italic>Aeromonas</italic> species were much more abundant in bats that feed on both insects and fish than in those with insect diet, and <italic>Micromonospora</italic> sp. RP3T was found to be enriched in bats that feed on both insects and fish (<xref rid="fig4" ref-type="fig">Figure 4B</xref>).</p>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption>
<p>Significantly changed bacterial species and metabolic pathways between IB and IFB. <bold>(A)</bold> The 62 significantly changed bacterial species between IB and IFB. The phylogenetic relationships of these species are shown in the right panel. <bold>(B)</bold> Significantly altered KEGG pathways between IB and IFB.</p>
</caption>
<graphic xlink:href="fmicb-14-1207482-g004.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussions" id="sec13">
<title>Discussion</title>
<p>The compositions of gut microbiota vary among different species of animals and may also vary under different physiological and behavioral conditions. As bats are highly differentiated animals with various dietary preferences and some bat species such as <italic>M. pilosus</italic> are in transition from insect-eating to fish-eating, they are ideal for investigation of the effects of diets on gut microbiota. In this study, we found that the gut microbiota profile of the one frugivorous bat (<italic>R. leschenaultii</italic>) examined was greatly different from that of insectivorous bats and bats that feed on both insects and fish, but those of the latter two dietary types of bats were similar. However, the abundance of 62 different species of bacteria was found to vary among different species of bats. These observations suggest that the difference in the profiles of gut microbiota in bats of different phylogenies was higher than those with different dietary habits as previously described (<xref ref-type="bibr" rid="ref28">Lutz et al., 2019</xref>). We also found a much less difference in gut microbiota among bats of different phylogenies at the genus (<italic>R</italic><sup>2</sup>&#x2009;=&#x2009;0.54, <italic>p</italic>&#x2009;=&#x2009;0.06) level than at the family level (<italic>R</italic><sup>2</sup>&#x2009;=&#x2009;0.049, <italic>p</italic>&#x2009;=&#x2009;0.743) by PCA analysis. A possible reason for such is that closely related bats have similar behaviors, physiologies, and environmental factors and therefore have similar gut microbial compositions (<xref ref-type="bibr" rid="ref3">Archie, 2019</xref>).</p>
<p>Many of the physiological or behavioral characteristics of bats have a unique evolutionary history (<xref ref-type="bibr" rid="ref7">Brualla et al., 2023</xref>). For example, the basic immune systems of bats (e.g., innate and adaptive immune systems) are more tolerant than those of other mammals (<xref ref-type="bibr" rid="ref30">Randolph and Barreiro, 2018</xref>). A study of bat rabies virus found that bat antibodies failed to penetrate the blood-brain barrier and therefore were unable to defend the infection in the central nervous system (<xref ref-type="bibr" rid="ref31">Roy et al., 2007</xref>). However, bats are known to carry many potentially lethal viruses, such as filoviruses, henipaviruses, lyssaviruses, paramyxoviruses, and coronaviruses (<xref ref-type="bibr" rid="ref20">Karunarathna et al., 2020</xref>; <xref ref-type="bibr" rid="ref36">Sun et al., 2020</xref>). Such symbiosis of various viruses may affect their physiologies, behaviors, and adaption to the virulent gut ecology. Bats also harbor many bacterial pathogens such as various species of <italic>Bacillus</italic>, <italic>Enterobacter</italic>, <italic>Enterococcus</italic>, <italic>Escherichia</italic>, <italic>Klebsiella</italic>, <italic>Pantoea</italic>, <italic>Pseudomonas</italic>, and <italic>Serratia</italic> (<xref ref-type="bibr" rid="ref13">Daniel et al., 2013</xref>; <xref ref-type="bibr" rid="ref36">Sun et al., 2020</xref>). We postulate that the symbiosis of various viruses and pathogenic bacteria is a survival strategy for bats to evade predators and that the host-microbe interactions may be a driving force for the evolution of bats to survive in various ecologic conditions. Dietary habits and other behaviors may also affect animal evolution. For instance, the bamboo diet promotes the adaptive development of strong teeth and paws in giant panda (<xref ref-type="bibr" rid="ref37">Vallittu et al., 2021</xref>). Another example is that the flying ability of birds or bats is endowed with light bone, high degree of skull healing, shoulder strap, sternum with keel protuberance, and degradation of posterior tibia and fibula (<xref ref-type="bibr" rid="ref11">Chang et al., 2019</xref>).</p>
<p>Different from most terrestrial mammals, the gut microbiota of bats with different dietary habits have few Bacteroidetes, but similar to carnivores, bat gut microbiota is enriched in Firmicutes and Proteobacteria (<xref ref-type="bibr" rid="ref36">Sun et al., 2020</xref>). A previous study of bat gut bacteria showed that species of Gammaproteobacteria and Firmicutes were dominant in non-breeding adult bats; this observation is also consistent with ours in this study. We also did not observe a significant correlation between dietary habits and gut microbiota in bats that feed on insets or both insects and fish as that reported by other investigators (<xref ref-type="bibr" rid="ref29">Phillips et al., 2012</xref>; <xref ref-type="bibr" rid="ref10">Carrillo-Araujo et al., 2015</xref>; <xref ref-type="bibr" rid="ref28">Lutz et al., 2019</xref>). We noticed that <italic>Aeromonas</italic> species were much more abundant in bats that feed on both insects and fish than those with insect diet. As <italic>Aeromonas</italic> is often detected in fish (<xref ref-type="bibr" rid="ref15">Hanninen et al., 1997</xref>), it may alter the composition of gut microbiota in bats that feed on fish.</p>
<p>As observed in previous studies (<xref ref-type="bibr" rid="ref42">Zepeda Mendoza et al., 2018</xref>; <xref ref-type="bibr" rid="ref17">Ingala et al., 2021</xref>), we found that the predominant metabolic pathways were Transporters for insectivorous bats; Transporters and Carbon metabolism for bats that feed on both insects and fish; and Transporters and Biosynthesis of amino acids for frugivorous bat. A recent study showed that the numbers of genes involved in 37 metabolic pathways were different between herbivorous and animalivorous bats (<xref ref-type="bibr" rid="ref17">Ingala et al., 2021</xref>). However, we found only 3 metabolic pathways involved in stilbenoid, diarylheptanoid and gingerol biosynthesis that were different in the number of genes between insectivorous bats and bats (<italic>M. pilosus</italic>) that feed on both insects and fish. A probable reason for this discrepancy is that we examined a very small dietary shift, i.e., from insect to both insect and fish diet. As <italic>M. pilosus</italic> bats have a high viral load, the enrichment of genes in these 3 metabolic pathways may be a symbiotic strategy with the viruses. The finding that the number of genes involved in photosynthesis of antenna proteins varied between IB and IFB (<italic>p</italic>&#x2009;=&#x2009;0.036) suggests that photosynthesis also plays a role in gut microbial variation among bats with different dietary habits as previously postulated (<xref ref-type="bibr" rid="ref19">Jin et al., 2021</xref>).</p>
<p>By investigating the effects of host speciation and dietary habits on gut microbiota, this study disclosed the effect of host&#x2013;microbe interactions on bat evolution. A major weakness in this study is in bat sampling. To avoid the influence of geographic variation on gut microbiota, bats used in this study were captured from the same cave. As a result, we only captured bats from 3 families (Pteropodidae, Rhinolophidae, Vespertilionidae), and the number of bats is very limited with only one bat in most species. Such problem may introduce deviations in the comparison of microbial diversity. Furthermore, a previous study surveyed the variations in gut microbiota among 9 bat families, and nearly half of the families showed significant differences in gut bacterial composition (<xref ref-type="bibr" rid="ref28">Lutz et al., 2019</xref>). We investigated 3 bat families and detected no significant difference in gut microbiota among them, indicating that investigations with more families are warranted in future studies. Although we have investigated the effects of host speciation and diets on the composition of gut microbiota, there may be other host factors that also play a role in the formation of gut microbiota in bats. We also noticed that more than 50% of gut microbial genes in each sample were annotated as &#x201C;function unknown.&#x201D; It has been postulated that these genes may be involved in echolocation, longevity, flying, and immune system (<xref ref-type="bibr" rid="ref43">Zhang et al., 2013</xref>; <xref ref-type="bibr" rid="ref18">Jebb et al., 2020</xref>).</p>
</sec>
<sec sec-type="conclusions" id="sec14">
<title>Conclusion</title>
<p>In this study, we found a very weak correlation between dietary habits and gut microbiota in bats that feed on insects or both insects and fish. However, the compositions of gut microbiota were found to vary among different species of bats, suggesting that host speciation plays a major role in the formation of gut microbiota. We also found differences in the number of genes related to certain metabolic pathways among bats with different dietary habits, suggesting the effect of diet on physiological functions of bats.</p>
</sec>
<sec sec-type="data-availability" id="sec15">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: <ext-link xlink:href="https://www.ncbi.nlm.nih.gov/" ext-link-type="uri">https://www.ncbi.nlm.nih.gov/</ext-link>, PRJNA954561.</p>
</sec>
<sec id="sec16">
<title>Ethics statement</title>
<p>The animal study was reviewed and approved by the Animal Ethics Committee of Institute of Zoology, Guangdong Academy of Sciences (No. GIABR20200810).</p>
</sec>
<sec id="sec17">
<title>Author contributions</title>
<p>MG and LZ designed the research. MG performed the research. JW, YZ, and XH collected samples and extracted fecal DNA. MG, SX, JQ, PL, CZ, and CH provided analytic tools. MG, SX, JD, and JQ analyzed data. MG and LZ wrote the paper. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec sec-type="funding-information" id="sec18">
<title>Funding</title>
<p>This work was supported by grants from the GDAS Special Project of Science and Technology Development (2021GDASYL-20210103052), the Guangdong Provincial Science and Technology Program (2021B1212110003 and 2021B1212050021), and the National Natural Science Foundation of China (32172435).</p>
</sec>
<sec sec-type="COI-statement" id="sec19">
<title>Conflict of interest</title>
<p>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.</p>
<p>The handling editor DJ declared a shared parent affiliation with some of the authors MG, JW, YZ, XH, PL, JD, CZ, and LZ at the time of review.</p>
</sec>
<sec id="sec100" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<ack>
<p>The authors thank Feifan Chen (Zhongkai University of Agriculture and Engineering) and Liwei Liu (Guangzhou University) for their advice in the study.</p>
</ack>
<sec sec-type="supplementary-material" id="sec20">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2023.1207482/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fmicb.2023.1207482/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Image_1.TIF" id="SM1" mimetype="image/tif" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image_2.TIF" id="SM2" mimetype="image/tif" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image_3.TIF" id="SM3" mimetype="image/tif" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image_4.TIF" id="SM4" mimetype="image/tif" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Adams</surname> <given-names>R. A.</given-names></name> <name><surname>Pedersen</surname> <given-names>S. C.</given-names></name></person-group> (<year>2013</year>). <source>Bat evolution, ecology and conservation</source>. <publisher-loc>New York</publisher-loc>: <publisher-name>Springer</publisher-name>.</citation></ref>
<ref id="ref2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname> <given-names>M. J.</given-names></name></person-group> (<year>2001</year>). <article-title>A new method for non-parametric multivariate analysis of variance</article-title>. <source>Austral Ecol.</source> <volume>26</volume>, <fpage>32</fpage>&#x2013;<lpage>46</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1442-9993.2001.01070.pp.x</pub-id></citation></ref>
<ref id="ref3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Archie</surname> <given-names>E. A.</given-names></name></person-group> (<year>2019</year>). <article-title>Bat microbiomes are socially synchronized</article-title>. <source>Nat. Ecol. Evol.</source> <volume>3</volume>, <fpage>18</fpage>&#x2013;<lpage>19</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41559-018-0749-2</pub-id>, PMID: <pub-id pub-id-type="pmid">30532041</pub-id></citation></ref>
<ref id="ref4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benson</surname> <given-names>A. K.</given-names></name> <name><surname>Kelly</surname> <given-names>S. A.</given-names></name> <name><surname>Legge</surname> <given-names>R.</given-names></name> <name><surname>Ma</surname> <given-names>F.</given-names></name> <name><surname>Low</surname> <given-names>S. J.</given-names></name> <name><surname>Kim</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2010</year>). <article-title>Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>107</volume>, <fpage>18933</fpage>&#x2013;<lpage>18938</lpage>. doi: <pub-id pub-id-type="doi">10.1073/pnas.1007028107</pub-id>, PMID: <pub-id pub-id-type="pmid">20937875</pub-id></citation></ref>
<ref id="ref5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Breiman</surname> <given-names>L.</given-names></name></person-group> (<year>2001</year>). <article-title>Random forests</article-title>. <source>Mach. Learn.</source> <volume>2001</volume>, <volume>45</volume>, <fpage>5</fpage>&#x2013;<lpage>32</lpage>. doi: <pub-id pub-id-type="doi">10.1023/A:1010933404324</pub-id></citation></ref>
<ref id="ref6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brooks</surname> <given-names>A. W.</given-names></name> <name><surname>Kohl</surname> <given-names>K. D.</given-names></name> <name><surname>Brucker</surname> <given-names>R. M.</given-names></name> <name><surname>van Opstal</surname> <given-names>E. J.</given-names></name> <name><surname>Bordenstein</surname> <given-names>S. R.</given-names></name></person-group> (<year>2016</year>). <article-title>Phylosymbiosis: relationships and functional effects of microbial communities across host evolutionary history</article-title>. <source>PLoS Biol.</source> <volume>14</volume>:<fpage>e2000225</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pbio.2000225</pub-id>, PMID: <pub-id pub-id-type="pmid">27861590</pub-id></citation></ref>
<ref id="ref7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brualla</surname> <given-names>N. L. M.</given-names></name> <name><surname>Wilson</surname> <given-names>L. A. B.</given-names></name> <name><surname>Doube</surname> <given-names>M.</given-names></name> <name><surname>Carter</surname> <given-names>R. T.</given-names></name> <name><surname>McElligott</surname> <given-names>A. G.</given-names></name> <name><surname>Koyabu</surname> <given-names>D.</given-names></name></person-group> (<year>2023</year>). <article-title>The vocal apparatus: an understudied tool to reconstruct the evolutionary history of echolocation in bats?</article-title> <source>J. Mamm. Evol.</source> <volume>30</volume>, <fpage>79</fpage>&#x2013;<lpage>94</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s10914-022-09647-z</pub-id></citation></ref>
<ref id="ref8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brucker</surname> <given-names>R. M.</given-names></name> <name><surname>Bordenstein</surname> <given-names>S. R.</given-names></name></person-group> (<year>2012</year>). <article-title>The roles of host evolutionary relationships (genus: <italic>Nasonia</italic>) and development in structuring microbial communities</article-title>. <source>Evolution</source> <volume>66</volume>, <fpage>349</fpage>&#x2013;<lpage>362</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1558-5646.2011.01454.x</pub-id>, PMID: <pub-id pub-id-type="pmid">22276533</pub-id></citation></ref>
<ref id="ref9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Camacho</surname> <given-names>C.</given-names></name> <name><surname>Coulouris</surname> <given-names>G.</given-names></name> <name><surname>Avagyan</surname> <given-names>V.</given-names></name> <name><surname>Ma</surname> <given-names>N.</given-names></name> <name><surname>Papadopoulos</surname> <given-names>J.</given-names></name> <name><surname>Bealer</surname> <given-names>K.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>BLAST+: architecture and applications</article-title>. <source>BMC Bioinformatics</source> <volume>10</volume>:<fpage>421</fpage>. doi: <pub-id pub-id-type="doi">10.1186/1471-2105-10-421</pub-id>, PMID: <pub-id pub-id-type="pmid">20003500</pub-id></citation></ref>
<ref id="ref10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carrillo-Araujo</surname> <given-names>M.</given-names></name> <name><surname>Ta&#x015F;</surname> <given-names>N.</given-names></name> <name><surname>Alc&#x00E1;ntara-Hern&#x00E1;ndez</surname> <given-names>R. J.</given-names></name> <name><surname>Gaona</surname> <given-names>O.</given-names></name> <name><surname>Schondube</surname> <given-names>J. E.</given-names></name> <name><surname>Medell&#x00ED;n</surname> <given-names>R. A.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Phyllostomid bat microbiome composition is associated to host phylogeny and feeding strategies</article-title>. <source>Front. Microbiol.</source> <volume>6</volume>:<fpage>447</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2015.00447</pub-id>, PMID: <pub-id pub-id-type="pmid">26042099</pub-id></citation></ref>
<ref id="ref11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname> <given-names>Y.</given-names></name> <name><surname>Song</surname> <given-names>S.</given-names></name> <name><surname>Li</surname> <given-names>A.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>Z.</given-names></name> <name><surname>Xiao</surname> <given-names>Y.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>The roles of morphological traits, resource variation and resource partitioning associated with the dietary niche expansion in the fish-eating bat <italic>Myotis pilosus</italic></article-title>. <source>Mol. Ecol.</source> <volume>28</volume>, <fpage>2944</fpage>&#x2013;<lpage>2954</lpage>. doi: <pub-id pub-id-type="doi">10.1111/mec.15127</pub-id>, PMID: <pub-id pub-id-type="pmid">31063664</pub-id></citation></ref>
<ref id="ref12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>Y.</given-names></name> <name><surname>Shi</surname> <given-names>C.</given-names></name> <name><surname>Huang</surname> <given-names>Z.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data</article-title>. <source>Gigascience</source> <volume>7</volume>, <fpage>1</fpage>&#x2013;<lpage>6</lpage>. doi: <pub-id pub-id-type="doi">10.1093/gigascience/gix120</pub-id>, PMID: <pub-id pub-id-type="pmid">29220494</pub-id></citation></ref>
<ref id="ref13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Daniel</surname> <given-names>D. S.</given-names></name> <name><surname>Ng</surname> <given-names>Y. K.</given-names></name> <name><surname>Chua</surname> <given-names>E. L.</given-names></name> <name><surname>Arumugam</surname> <given-names>Y.</given-names></name> <name><surname>Wong</surname> <given-names>W. L.</given-names></name> <name><surname>Kumaran</surname> <given-names>J. V.</given-names></name></person-group> (<year>2013</year>). <article-title>Isolation and identification of gastrointestinal microbiota from the short-nosed fruit bat <italic>Cynopterus brachyotis</italic></article-title>. <source>Microbiol. Res.</source> <volume>168</volume>, <fpage>485</fpage>&#x2013;<lpage>496</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.micres.2013.04.001</pub-id></citation></ref>
<ref id="ref14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Easson</surname> <given-names>C. G.</given-names></name> <name><surname>Thacker</surname> <given-names>R. W.</given-names></name></person-group> (<year>2014</year>). <article-title>Phylogenetic signal in the community structure of host-specific microbiomes of tropical marine sponges</article-title>. <source>Front. Microbiol.</source> <volume>5</volume>:<fpage>532</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2014.00532</pub-id>, PMID: <pub-id pub-id-type="pmid">25368606</pub-id></citation></ref>
<ref id="ref15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hanninen</surname> <given-names>M. L.</given-names></name> <name><surname>Oivanen</surname> <given-names>P.</given-names></name> <name><surname>Hirvela-Koski</surname> <given-names>V.</given-names></name></person-group> (<year>1997</year>). <article-title>Aeromonas species in fish, fish-eggs, shrimp and freshwater</article-title>. <source>Int. J. Food Microbiol.</source> <volume>34</volume>, <fpage>17</fpage>&#x2013;<lpage>26</lpage>. doi: <pub-id pub-id-type="doi">10.1016/s0168-1605(96)01163-4</pub-id>, PMID: <pub-id pub-id-type="pmid">9029253</pub-id></citation></ref>
<ref id="ref16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hong</surname> <given-names>C.</given-names></name> <name><surname>Manimaran</surname> <given-names>S.</given-names></name> <name><surname>Shen</surname> <given-names>Y.</given-names></name> <name><surname>Perez-Rogers</surname> <given-names>J. F.</given-names></name> <name><surname>Byrd</surname> <given-names>A. L.</given-names></name> <name><surname>Castro-Nallar</surname> <given-names>E.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>PathoScope 2.0: a complete computational framework for strain identification in environmental or clinical sequencing samples</article-title>. <source>Microbiome</source> <volume>2</volume>:<fpage>33</fpage>. doi: <pub-id pub-id-type="doi">10.1186/2049-2618-2-33</pub-id>, PMID: <pub-id pub-id-type="pmid">25225611</pub-id></citation></ref>
<ref id="ref17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ingala</surname> <given-names>M. R.</given-names></name> <name><surname>Simmons</surname> <given-names>N. B.</given-names></name> <name><surname>Dunbar</surname> <given-names>M.</given-names></name> <name><surname>Wultsch</surname> <given-names>C.</given-names></name> <name><surname>Krampis</surname> <given-names>K.</given-names></name> <name><surname>Perkins</surname> <given-names>S. L.</given-names></name></person-group> (<year>2021</year>). <article-title>You are more than what you eat: potentially adaptive enrichment of microbiome functions across bat dietary niches</article-title>. <source>Anim. Microbiome</source> <volume>3</volume>:<fpage>82</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s42523-021-00139-8</pub-id>, PMID: <pub-id pub-id-type="pmid">34906258</pub-id></citation></ref>
<ref id="ref18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jebb</surname> <given-names>D.</given-names></name> <name><surname>Huang</surname> <given-names>Z.</given-names></name> <name><surname>Pippel</surname> <given-names>M.</given-names></name> <name><surname>Hughes</surname> <given-names>G. M.</given-names></name> <name><surname>Lavrichenko</surname> <given-names>K.</given-names></name> <name><surname>Devanna</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Six reference-quality genomes reveal evolution of bat adaptations</article-title>. <source>Nature</source> <volume>583</volume>, <fpage>578</fpage>&#x2013;<lpage>584</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41586-020-2486-3</pub-id>, PMID: <pub-id pub-id-type="pmid">32699395</pub-id></citation></ref>
<ref id="ref19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jin</surname> <given-names>H.</given-names></name> <name><surname>Wang</surname> <given-names>Y.</given-names></name> <name><surname>Zhao</surname> <given-names>P.</given-names></name> <name><surname>Wang</surname> <given-names>L.</given-names></name> <name><surname>Zhang</surname> <given-names>S.</given-names></name> <name><surname>Meng</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Potential of producing flavonoids using cyanobacteria as a sustainable chassis</article-title>. <source>J. Agric. Food Chem.</source> <volume>69</volume>, <fpage>12385</fpage>&#x2013;<lpage>12401</lpage>. doi: <pub-id pub-id-type="doi">10.1021/acs.jafc.1c04632</pub-id>, PMID: <pub-id pub-id-type="pmid">34649432</pub-id></citation></ref>
<ref id="ref20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Karunarathna</surname> <given-names>S. C.</given-names></name> <name><surname>Dong</surname> <given-names>Y.</given-names></name> <name><surname>Karasaki</surname> <given-names>S.</given-names></name> <name><surname>Tibpromma</surname> <given-names>S.</given-names></name> <name><surname>Hyde</surname> <given-names>K. D.</given-names></name> <name><surname>Lumyong</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Discovery of novel fungal species and pathogens on bat carcasses in a cave in Yunnan Province, China</article-title>. <source>Emerg. Microbes. Infect.</source> <volume>9</volume>, <fpage>1554</fpage>&#x2013;<lpage>1566</lpage>. doi: <pub-id pub-id-type="doi">10.1080/22221751.2020.1785333</pub-id>, PMID: <pub-id pub-id-type="pmid">32573334</pub-id></citation></ref>
<ref id="ref21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kuang</surname> <given-names>Z.</given-names></name> <name><surname>Li</surname> <given-names>F.</given-names></name> <name><surname>Duan</surname> <given-names>Q.</given-names></name> <name><surname>Tian</surname> <given-names>C.</given-names></name> <name><surname>Nevo</surname> <given-names>E.</given-names></name> <name><surname>Li</surname> <given-names>K.</given-names></name></person-group> (<year>2022</year>). <article-title>Host diet shapes functionally differentiated gut microbiomes in sympatric speciation of blind mole rats in upper Galilee, Israel</article-title>. <source>Front. Microbiol.</source> <volume>13</volume>:<fpage>1062763</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2022.1062763</pub-id>, PMID: <pub-id pub-id-type="pmid">36458196</pub-id></citation></ref>
<ref id="ref22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumar</surname> <given-names>S.</given-names></name> <name><surname>Suleski</surname> <given-names>M.</given-names></name> <name><surname>Craig</surname> <given-names>J. M.</given-names></name> <name><surname>Kasprowicz</surname> <given-names>A. E.</given-names></name> <name><surname>Sanderford</surname> <given-names>M.</given-names></name> <name><surname>Li</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>TimeTree 5: an expanded resource for species divergence times</article-title>. <source>Mol. Biol. Evol.</source> <volume>39</volume>:<fpage>msac174</fpage>. doi: <pub-id pub-id-type="doi">10.1093/molbev/msac174</pub-id>, PMID: <pub-id pub-id-type="pmid">35932227</pub-id></citation></ref>
<ref id="ref23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kurilshikov</surname> <given-names>A.</given-names></name> <name><surname>Medina-Gomez</surname> <given-names>C.</given-names></name> <name><surname>Bacigalupe</surname> <given-names>R.</given-names></name> <name><surname>Radjabzadeh</surname> <given-names>D.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Demirkan</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Large-scale association analyses identify host factors influencing human gut microbiome composition</article-title>. <source>Nat. Genet.</source> <volume>53</volume>, <fpage>156</fpage>&#x2013;<lpage>165</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41588-020-00763-1</pub-id>, PMID: <pub-id pub-id-type="pmid">33462485</pub-id></citation></ref>
<ref id="ref24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Langmead</surname> <given-names>B.</given-names></name> <name><surname>Salzberg</surname> <given-names>S. L.</given-names></name></person-group> (<year>2012</year>). <article-title>Fast gapped-read alignment with Bowtie 2</article-title>. <source>Nat. Methods</source> <volume>9</volume>, <fpage>357</fpage>&#x2013;<lpage>359</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.1923</pub-id>, PMID: <pub-id pub-id-type="pmid">22388286</pub-id></citation></ref>
<ref id="ref25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>W.</given-names></name> <name><surname>Godzik</surname> <given-names>A.</given-names></name></person-group> (<year>2006</year>). <article-title>Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences</article-title>. <source>Bioinformatics</source> <volume>22</volume>, <fpage>1658</fpage>&#x2013;<lpage>1659</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btl158</pub-id></citation></ref>
<ref id="ref26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>D.</given-names></name> <name><surname>Liu</surname> <given-names>C. M.</given-names></name> <name><surname>Luo</surname> <given-names>R.</given-names></name> <name><surname>Sadakane</surname> <given-names>K.</given-names></name> <name><surname>Lam</surname> <given-names>T. W.</given-names></name></person-group> (<year>2015</year>). <article-title>MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph</article-title>. <source>Bioinformatics</source> <volume>31</volume>, <fpage>1674</fpage>&#x2013;<lpage>1676</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btv033</pub-id>, PMID: <pub-id pub-id-type="pmid">25609793</pub-id></citation></ref>
<ref id="ref27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lim</surname> <given-names>E. S.</given-names></name> <name><surname>Zhou</surname> <given-names>Y.</given-names></name> <name><surname>Zhao</surname> <given-names>G.</given-names></name> <name><surname>Bauer</surname> <given-names>I. K.</given-names></name> <name><surname>Droit</surname> <given-names>L.</given-names></name> <name><surname>Ndao</surname> <given-names>I. M.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Early life dynamics of the human gut virome and bacterial microbiome in infants</article-title>. <source>Nat. Med.</source> <volume>21</volume>, <fpage>1228</fpage>&#x2013;<lpage>1234</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nm.3950</pub-id>, PMID: <pub-id pub-id-type="pmid">26366711</pub-id></citation></ref>
<ref id="ref28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lutz</surname> <given-names>H. L.</given-names></name> <name><surname>Jackson</surname> <given-names>E. W.</given-names></name> <name><surname>Webala</surname> <given-names>P. W.</given-names></name> <name><surname>Babyesiza</surname> <given-names>W. S.</given-names></name> <name><surname>Kerbis Peterhans</surname> <given-names>J. C.</given-names></name> <name><surname>Demos</surname> <given-names>T. C.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Ecology and host identity outweigh evolutionary history in shaping the bat microbiome</article-title>. <source>mSystems</source> <volume>4</volume>:<fpage>e00511</fpage>. doi: <pub-id pub-id-type="doi">10.1128/mSystems.00511-19</pub-id>, PMID: <pub-id pub-id-type="pmid">31719140</pub-id></citation></ref>
<ref id="ref29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Phillips</surname> <given-names>C. D.</given-names></name> <name><surname>Phelan</surname> <given-names>G.</given-names></name> <name><surname>Dowd</surname> <given-names>S. E.</given-names></name> <name><surname>McDonough</surname> <given-names>M. M.</given-names></name> <name><surname>Ferguson</surname> <given-names>A. W.</given-names></name> <name><surname>Delton hanson</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Microbiome analysis among bats describes influences of host phylogeny, life history, physiology and geography</article-title>. <source>Mol. Ecol.</source> <volume>21</volume>, <fpage>2617</fpage>&#x2013;<lpage>2627</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1365-294X.2012.05568.x</pub-id>, PMID: <pub-id pub-id-type="pmid">22519571</pub-id></citation></ref>
<ref id="ref30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Randolph</surname> <given-names>H. E.</given-names></name> <name><surname>Barreiro</surname> <given-names>L. B.</given-names></name></person-group> (<year>2018</year>). <article-title>Holy immune tolerance, batman!</article-title> <source>Immunity</source> <volume>48</volume>, <fpage>1074</fpage>&#x2013;<lpage>1076</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.immuni.2018.05.016</pub-id>, PMID: <pub-id pub-id-type="pmid">29924972</pub-id></citation></ref>
<ref id="ref31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Roy</surname> <given-names>A.</given-names></name> <name><surname>Phares</surname> <given-names>T. W.</given-names></name> <name><surname>Koprowski</surname> <given-names>H.</given-names></name> <name><surname>Hooper</surname> <given-names>D. C.</given-names></name></person-group> (<year>2007</year>). <article-title>Failure to open the blood-brain barrier and deliver immune effectors to central nervous system tissues leads to the lethal outcome of silver-haired bat rabies virus infection</article-title>. <source>J. Virol.</source> <volume>81</volume>, <fpage>1110</fpage>&#x2013;<lpage>1118</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JVI.01964-06</pub-id>, PMID: <pub-id pub-id-type="pmid">17108029</pub-id></citation></ref>
<ref id="ref32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sanders</surname> <given-names>J. G.</given-names></name> <name><surname>Powell</surname> <given-names>S.</given-names></name> <name><surname>Kronauer</surname> <given-names>D. J. C.</given-names></name> <name><surname>Vasconcelos</surname> <given-names>H. L.</given-names></name> <name><surname>Frederickson</surname> <given-names>M. E.</given-names></name> <name><surname>Pierce</surname> <given-names>N. E.</given-names></name></person-group> (<year>2014</year>). <article-title>Stability and phylogenetic correlation in gut microbiota: lessons from ants and apes</article-title>. <source>Mol. Ecol.</source> <volume>23</volume>, <fpage>1268</fpage>&#x2013;<lpage>1283</lpage>. doi: <pub-id pub-id-type="doi">10.1111/mec.12611</pub-id>, PMID: <pub-id pub-id-type="pmid">24304129</pub-id></citation></ref>
<ref id="ref33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shalabi</surname> <given-names>L. A.</given-names></name> <name><surname>Shaaban</surname> <given-names>Z.</given-names></name> <name><surname>Kasasbeh</surname> <given-names>B.</given-names></name></person-group> (<year>2006</year>). <article-title>Data mining: a preprocessing engine</article-title>. <source>J. Comput. Sci.</source> <volume>2</volume>, <fpage>735</fpage>&#x2013;<lpage>739</lpage>. doi: <pub-id pub-id-type="doi">10.3844/jcssp.2006.735.739</pub-id></citation></ref>
<ref id="ref34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shannon</surname> <given-names>C. E.</given-names></name></person-group> (<year>2001</year>). <article-title>A mathematical theory of communication</article-title>. <source>ACM SIGMOBILE Mob. Comput. Commun. Rev.</source> <volume>5</volume>, <fpage>3</fpage>&#x2013;<lpage>55</lpage>. doi: <pub-id pub-id-type="doi">10.1145/584091.584093</pub-id></citation></ref>
<ref id="ref35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>S. J.</given-names></name> <name><surname>Sanders</surname> <given-names>J. G.</given-names></name> <name><surname>Delsuc</surname> <given-names>F.</given-names></name> <name><surname>Metcalf</surname> <given-names>J.</given-names></name> <name><surname>Amato</surname> <given-names>K.</given-names></name> <name><surname>Taylor</surname> <given-names>M. W.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Comparative analyses of vertebrate gut microbiomes reveal convergence between birds and bats</article-title>. <source>mBio</source> <volume>11</volume>:<fpage>e02901</fpage>. doi: <pub-id pub-id-type="doi">10.1128/mBio.02901-19</pub-id>, PMID: <pub-id pub-id-type="pmid">31911491</pub-id></citation></ref>
<ref id="ref36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sun</surname> <given-names>D. L.</given-names></name> <name><surname>Gao</surname> <given-names>Y. Z.</given-names></name> <name><surname>Ge</surname> <given-names>X. Y.</given-names></name> <name><surname>Shi</surname> <given-names>Z. L.</given-names></name> <name><surname>Zhou</surname> <given-names>N. Y.</given-names></name></person-group> (<year>2020</year>). <article-title>Special features of bat microbiota differ from those of terrestrial mammals</article-title>. <source>Front. Microbiol.</source> <volume>11</volume>:<fpage>1040</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2020.01040</pub-id>, PMID: <pub-id pub-id-type="pmid">32582057</pub-id></citation></ref>
<ref id="ref37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vallittu</surname> <given-names>P. K.</given-names></name> <name><surname>Varrela</surname> <given-names>J.</given-names></name> <name><surname>Salo</surname> <given-names>J.</given-names></name> <name><surname>Rengui</surname> <given-names>L.</given-names></name> <name><surname>Shanshan</surname> <given-names>L.</given-names></name> <name><surname>Shan</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Temporomandibular joint and Giant Panda&#x2019;s (<italic>Ailuropoda melanoleuca</italic>) adaptation to bamboo diet</article-title>. <source>Sci. Rep.</source> <volume>11</volume>:<fpage>14252</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-021-93808-2</pub-id>, PMID: <pub-id pub-id-type="pmid">34244613</pub-id></citation></ref>
<ref id="ref38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wan</surname> <given-names>Y.</given-names></name> <name><surname>Tang</surname> <given-names>J.</given-names></name> <name><surname>Li</surname> <given-names>J.</given-names></name> <name><surname>Yuan</surname> <given-names>J.</given-names></name> <name><surname>Wang</surname> <given-names>F.</given-names></name> <name><surname>Li</surname> <given-names>D.</given-names></name></person-group> (<year>2020</year>). <article-title>Contribution of diet to gut microbiota and related host cardiometabolic health: diet-gut interaction in human health</article-title>. <source>Gut Microbes</source> <volume>11</volume>, <fpage>603</fpage>&#x2013;<lpage>609</lpage>. doi: <pub-id pub-id-type="doi">10.1080/19490976.2019.1697149</pub-id>, PMID: <pub-id pub-id-type="pmid">31964203</pub-id></citation></ref>
<ref id="ref39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>White</surname> <given-names>J. R.</given-names></name> <name><surname>Nagarajan</surname> <given-names>N.</given-names></name> <name><surname>Pop</surname> <given-names>M.</given-names></name></person-group> (<year>2009</year>). <article-title>Statistical methods for detecting differentially abundant features in clinical metagenomic samples</article-title>. <source>PLoS Comput. Biol.</source> <volume>5</volume>:<fpage>e1000352</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pcbi.1000352</pub-id>, PMID: <pub-id pub-id-type="pmid">19360128</pub-id></citation></ref>
<ref id="ref40"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Wilson</surname> <given-names>D. E.</given-names></name> <name><surname>Mittermeier</surname> <given-names>R. A.</given-names></name></person-group> (<year>2019</year>). <source>Handbook of the mammals of the world volume 9: bats</source>. <publisher-loc>Barcelona</publisher-loc>: <publisher-name>Lynx Edicions</publisher-name>.</citation></ref>
<ref id="ref41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yin</surname> <given-names>Z.</given-names></name> <name><surname>Sun</surname> <given-names>K.</given-names></name> <name><surname>Li</surname> <given-names>A.</given-names></name> <name><surname>Sun</surname> <given-names>D.</given-names></name> <name><surname>Li</surname> <given-names>Z.</given-names></name> <name><surname>Xiao</surname> <given-names>G.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Changes in the gut microbiota during Asian particolored bat (<italic>Vespertilio sinensis</italic>) development</article-title>. <source>PeerJ</source> <volume>8</volume>:<fpage>e9003</fpage>. doi: <pub-id pub-id-type="doi">10.7717/peerj.9003</pub-id>, PMID: <pub-id pub-id-type="pmid">32435532</pub-id></citation></ref>
<ref id="ref42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zepeda Mendoza</surname> <given-names>M. L.</given-names></name> <name><surname>Xiong</surname> <given-names>Z.</given-names></name> <name><surname>Escalera-Zamudio</surname> <given-names>M.</given-names></name> <name><surname>Runge</surname> <given-names>A. K.</given-names></name> <name><surname>Th&#x00E9;z&#x00E9;</surname> <given-names>J.</given-names></name> <name><surname>Streicker</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Hologenomic adaptations underlying the evolution of sanguivory in the common vampire bat</article-title>. <source>Nat. Ecol. Evol.</source> <volume>2</volume>, <fpage>659</fpage>&#x2013;<lpage>668</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41559-018-0476-8</pub-id>, PMID: <pub-id pub-id-type="pmid">29459707</pub-id></citation></ref>
<ref id="ref43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>G.</given-names></name> <name><surname>Cowled</surname> <given-names>C.</given-names></name> <name><surname>Shi</surname> <given-names>Z.</given-names></name> <name><surname>Huang</surname> <given-names>Z.</given-names></name> <name><surname>Bishop-Lilly</surname> <given-names>K. A.</given-names></name> <name><surname>Fang</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>Comparative analysis of bat genomes provides insight into the evolution of flight and immunity</article-title>. <source>Science</source> <volume>339</volume>, <fpage>456</fpage>&#x2013;<lpage>460</lpage>. doi: <pub-id pub-id-type="doi">10.1126/science.1230835</pub-id>, PMID: <pub-id pub-id-type="pmid">23258410</pub-id></citation></ref>
<ref id="ref44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname> <given-names>W.</given-names></name> <name><surname>Lomsadze</surname> <given-names>A.</given-names></name> <name><surname>Borodovsky</surname> <given-names>M.</given-names></name></person-group> (<year>2010</year>). <article-title><italic>Ab initio</italic> gene identification in metagenomic sequences</article-title>. <source>Nucleic Acids Res.</source> <volume>38</volume>:<fpage>e132</fpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkq275</pub-id>, PMID: <pub-id pub-id-type="pmid">20403810</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn id="fn0001">
<p><sup>1</sup><ext-link xlink:href="http://i.cs.hku.hk/~alse/hkubrg/projects/idba_ud/" ext-link-type="uri">http://i.cs.hku.hk/~alse/hkubrg/projects/idba_ud/</ext-link>
</p>
</fn>
<fn id="fn0002">
<p><sup>2</sup><ext-link xlink:href="http://exon.gatech.edu/meta_gmhmmp.cgi" ext-link-type="uri">http://exon.gatech.edu/meta_gmhmmp.cgi</ext-link>
</p>
</fn>
<fn id="fn0003">
<p><sup>3</sup><ext-link xlink:href="http://www.bioinformatics.org/cd-hit/" ext-link-type="uri">http://www.bioinformatics.org/cd-hit/</ext-link>
</p>
</fn>
<fn id="fn0004">
<p><sup>4</sup><ext-link xlink:href="http://blast.ncbi.nlm.nih.gov/Blast.cgi" ext-link-type="uri">http://blast.ncbi.nlm.nih.gov/Blast.cgi</ext-link>
</p>
</fn>
<fn id="fn0005">
<p><sup>5</sup><ext-link xlink:href="http://www.timetree.org/" ext-link-type="uri">http://www.timetree.org/</ext-link>
</p>
</fn>
</fn-group>
</back>
</article>