Differences in Diet and Gut Microbiota Between Lactating and Non-lactating Asian Particolored Bats (Vespertilio sinensis): Implication for a Connection Between Diet and Gut Microbiota

In mammals, lactation is considered the most energetically costly phase for females. To meet nutritional and energy demands, lactating females usually change feeding patterns by eating food that is higher in protein and calories. Their gut microbes respond accordingly to help adapt to the changes in diet. In this study, we examined differences in diet and gut microbial composition between lactating and non-lactating Asian particolored bats (Vespertilio sinensis) using COI and 16S amplicon sequencing. When compared with non-lactating bats, we found that the diversity and composition of lactating bats’ diets differed; the proportion of Diptera increased and Coleoptera and Orthoptera decreased significantly. This could be attributed to the easy availability and high protein content of Diptera. Comparative analysis of the gut microbiota of lactating and non-lactating females showed that although the diversity of gut microbiota did not change, the relative abundance of specific gut microbiota associated with a particular diet did change. For example, when the consumption of Coleoptera decreased in lactating bats, the relative abundance of Lactobacillaceae was also reduced. Lactobacillaceae are thought to be involved in the digestion of Coleopteran exoskeletons. This study suggests that during lactation, Asian particolored bats eat a diet that yields higher levels of protein, and at the same time, the abundance of specific gut microbes change to help their hosts adapt to these changes in diet.

Introduction

In the animal kingdom, females typically play a larger role in the care of young offspring. To maximize the chances of their offspring’s survival, they invest more energy and pay a higher parental cost by providing larger sex cells and they invest more energy into rearing them (Taylor et al., 2000). This phenomenon is most prominent in lactating female mammals. For example, in a study of lemurs, the females’ energy expenditure increased as their babies grew (Tarnaud, 2006). Similarly, in non-human primates, females showed higher dietary protein requirements during lactation (Oftedal, 1991). Therefore, lactation is the most costly investment in offspring by female mammals. To meet their infants’ needs for energy and nutrition, females face greater energetic and nutritional pressure during lactation (Aiello and Wells, 2002). This is especially true for the only flying mammals, bats. During reproduction, males and females bats of many species inhabit separate locations and females form separate breeding colonies to carry out the young offspring duties alone (Rydell and Baagøe, 1994; Safi et al., 2007). And during lactation, females require the high energy for milk production as well as for carrying their young while flying (Kunz, 1987; Sigsgaard et al., 2020).

Optimal foraging theory (OFT) suggests that an animal will attempt to gain the greatest energetic benefit for the lowest energetic cost, while foraging to maximize fitness (Schoener, 1971). Although, an individual can consume a wide diversity of prey items, it may adopt different diets and select specific foods depending on energetic benefits relative to handling time costs (Araujo et al., 2011). On the basis of the energetic demands of lactation in previous studies, lactating females take in more energy than females in other reproductive states (Herrera and Heymann, 2004; Wichert et al., 2009). They may accomplish this by eating for longer (Denryter et al., 2020), eating faster (McCabe and Fedigan, 2007), eating more nutritious food (Herrera and Heymann, 2004) or all three. Dietary changes are one of the most effective ways to cope with these energy and nutrition needs (Denryter et al., 2020). Many studies have emphasized that selection of high-quality foods, rather than bulk feeding, is the most important component of nutrient acquisition during lactation (Vasey, 2000, 2002). State-dependent diets, especially reproductive state-dependent diets, have been demonstrated in many mammal species (Thompson and Wrangham, 2008). To our knowledge, no study has yet examined whether the composition of lactating bats’ diets differs from that of non-lactating bats, and whether lactating bats choose to eat more high-protein foods in response to higher energy and nutritional requirements.

The gut microbiota is the most abundant microbial system co-occurring with mammals and is largely a consequence of both the nutrient-rich environment of mammalian digestive systems and the beneficial functions these communities provide to their hosts (Muegge et al., 2011). In mammals, the gut microbiota is expected to adapt to the varying energetic and nutritional pressures of females in different reproductive states. Recent studies implied that the community composition and diversity of the gut microbiota can covary with dietary changes, and that these changes are beneficial to the host, such that they can meet the energy and nutrition requirements of lactation (Phillips et al., 2017). So far, most research on lactating females has focused on the effects of the altered reproductive status and hormones on gut microbiota (Mallott et al., 2020; Sun et al., 2020); few studies have examined whether changes in diet composition during lactation cause corresponding changes in gut microbiota. To examine the food composition and gut microbial structure of mammals during lactation, and try to establish the relationship between diet and gut microbiota during their most energetically costly phase, we focused on an insectivorous bat species in northeast China: the Asian particolored bat (Vespertilio sinensis). In this study, COI and 16S rRNA amplicon sequencing was used to investigate the diet and gut microbiota of adult females in and out of lactation periods. We hypothesized that lactating bats, compared with non-lactating bats, change their diet composition to meet higher energy and nutritional requirements and their gut microbial diversity and composition will change accordingly. We predict that (1) there are significant differences in the composition of diet between lactating and non-lactating bats; (2) gut microbial diversity and composition of lactating bats are different from that in non-lactating bats; and (3) there is a strong correlation between diet and gut microbes.

Materials and Methods

Sample Collection

To avoid the influence of age and season (Ma et al., 2007), all female bats in this study were randomly collected from July to August 2020 in the Asian particolored bats habitat under the highway bridge of Acheng District, Harbin City, Heilongjiang Province, China (45°55'N, 126°94'E). Asian particolored bats mainly feed on insects and inhabit the roofs or eaves of bridges or old buildings. During reproduction, males and females inhabit separate locations, which provides ideal conditions for sampling. Bats were captured with a mist net as they returned to their habitat after hunting. We gently caught the bats from the net by hand with sterile gloves. We released the bats from the net as soon as possible and none were injured. We collected lactating (n=10) and non-lactating (n=14) Asian particolored bats in two periods. The determination of reproductive status was based on external morphological characteristics: the nipples of female bats become enlarged during lactation and the hairs around the nipples of female bats disappear because of suckling by their infants (Haarsma, 2008; Racey, 2009). Each captured individual was put into sterilized kraft bags to collect their excrement. After defecation, fecal samples were immediately put into an empty cryopreservation tube with sterile forceps and placed into dry ice for preservation until they were transported to the laboratory, where they were frozen at −80°C until DNA extraction. On completion of the sampling, the bats were immediately released back to their habitat.

DNA Extraction and PCR Amplification

FastDNA™ Spin Kit (Santa Ana, California, United States) was used to extract total insect and bacterial genomic DNA samples following the manufacturer’s instructions, and the extracted DNA was tested. The quantity of extracted DNA was measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States) and the quality was assessed using agarose gel electrophoresis.

Forward primer ZBJ-ArtF1cF and reverse primer ZBJ-ArtR2cR were used for Cytochrome oxidase fragments of PCR amplification, and forward primers 338F and reverse 806R were used for PCR amplification of the V3–V4 region of bacterial 16S rRNA gene. The PCR reaction system and conditions are described in detail in the Supplementary Material. PCR amplicons were purified with Agencourt AMPure Beads (Beckman Coulter, Indianapolis, IN, United States) and quantified using the PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, United States). After the individual quantification step, amplicons were pooled in equal amounts and pair-end 2×300-bp sequencing was performed using the Illlumina MiSeq platform with MiSeq Reagent Kit v3.

Bioinformatics Analysis

Raw fastq files were quality-filtered by Trimmomatic (Bolger et al., 2014) and merged by FLASH (Magoč and Salzberg, 2011). After performing quality control to obtain the optimized sequences, the operational taxonomic units (OTUs) were clustered using Usearch (Edgar, 2010) with singletons removed and chimera filtering. All optimized sequences were mapped to representative sequences, and sequences with >97% similarity with the representative sequence were selected to generate an OTU table. All analyses of the sequencing filtering and normalization were performed in Usearch. Sequencing effort coverage was visually assessed using alpha diversity rarefaction curves of number of OTUs.

Diet Data Analysis

Taxonomic identification of insect species was made by comparing a representative sequence from each OTU to reference sequences in the Barcode of Life Database (BOLD)¹ and the Genbank database.² By doing so, the insects in the bats’ diet could be identified to the species level. Sequence data analyses were mainly performed using QIIME v1.8.0 (Caporaso et al., 2010) and R packages (v3.5.1). OTU-level ranked abundance curves were generated to compare the richness and evenness of OTUs among samples. The Shannon diversity index and Chao1 richness index were calculated using the OTU table in QIIME to determine the diet alpha diversity for the groups of bats. We then used the Wilcoxon rank-sum test to calculate significant differences. Beta diversity was analyzed to investigate variation in the dietary composition of the different group samples using unweighted UniFrac distance metrics. Differences in group dispersion among bat groups were assessed using ADONIS, using R package, vegan (Jari Oksanen et al., 2019). Beta diversity was then visualized with principal coordinate analysis (PCoA; Ramette, 2007). The species composition and relative abundance of bat diets in each group were calculated at order taxonomic levels, and the composition of dominant species in different groups was visualized by community Pie plots. The relative abundance of insect species at order, family, and genus consumed by lactating and non-lactating bats were analyzed using Wilcoxon rank-sum test in R package, stats (R Core Team, 2019).

Gut Microbiota Data Analysis

QIIME software, Chao1 and Shannon indices were used to evaluate the alpha diversity of each sample. We then used the Wilcoxon rank-sum test to calculate significant differences between lactating and non-lactating bats. Beta diversity analysis was performed to investigate structural variation in microbial communities of the different group samples on the basis of unweighted UniFrac and weighted UniFrac distance metrics. Differences in group dispersion among bat groups were assessed using ADONIS. Beta diversity was then visualized with PCoA. Taxonomies were grouped at the phylum, family, and genus levels. Taxa abundances at the phylum, family, and genus levels were statistically compared among groups with using Wilcoxon rank-sum test, in R package, stats. To further identify the highly-dimensional gut microbes and characterize the differences between two or more biological conditions, we combined the linear discriminant analysis effect size (LEfSe) analysis.

Diet and Gut Microbiota Correlation Analysis

We performed canonical correspondence analysis (CCA) to test the relative abundance of different insect species in the diet on microbial community at different taxonomic levels in lactating and non-lactating bats, applying forward selection and the Monte Carlo permutation test with 999 random permutations. In addition, we used Procrustes analysis to assess congruence between different insect species relative abundance dietary composition in different groups on microbial community structure and the significance of the Procrustes statistic was tested by 999 permutations with “protest” function.

Results

Diet

After quality processing, a total of 293,784 effective sequences were obtained, with an average effective sequence of 12,241 per sample. We identified 276 OTUs from reads, 99 of which were only found in the diet of lactating bats and 138 of which were only found in the diet of non-lactating bats (Supplementary Table 1).

Dietary diversity differed significantly between the two reproductive periods. Compared with non-lactating bats, the alpha diversity of the diet based on the Shannon index was higher in lactating bats (Wilcoxon rank-sum test, p=0.004, Q=0.009; Figure 1A; Supplementary Figure 1A). However, the alpha diversity based on the Chao1 index was no different for lactating and non-lactating bats (Wilcoxon rank-sum test, p=0.23, Q=0.23; Figure 1B). The unweighted UniFrac metric data showed distinct clustering by lactating bat and non-lactating bat values, such that individuals from the two groups clustered separately. The results of ADONIS showed significant differences between the two groups (ADONIS: R²=0.14, p=0.001; Figure 2A).

FIGURE 1

Figure 1. Alpha diversity of different groups. (A,B) The diet diversity estimated by Chao1 and Shannon indexes. (C,D) Gut microbial diversity estimated by Chao1 and Shannon indexes. L refers to lactating bats; A refers to non-lactating bats. ^*0.01 < p ≤ 0.05, ^**0.001 < p ≤ 0.01, and ^***p ≤ 0.001.

FIGURE 2

Figure 2. Principal coordinate analysis (PCoA) analysis based on unweighted UniFrac metrics. (A) Differential insect composition in diet between lactating and non-lactating bats. (B) Differential gut microbiota communities between lactating and non-lactating bats.

Ten orders were found in the diet of lactating bats, with Diptera accounting for 67.64% and Lepidoptera for 23.74% (Figure 3A). Thirteen orders were found in the diet of non-lactating bats, with Diptera accounting for 40.24% and Lepidoptera for 26.94% (Figure 3B). We compared the relative read abundance of more than 1% of insect species and found that compared with non-lactating bats, lactating bats consumed more Diptera. The consumption of Coleoptera was significantly reduced (p=0.03), and the lactating bats even stopped consuming the larger Orthoptera (p=0.001; Figure 3C; Supplementary Table 3). At the family level, the proportion of Limoniidae (p=0.0002), Lasiocampidae (p=0.00002), and Limacodidae (p=0.04) were significantly increased and Carabidae (p=0.009) was significantly reduced in the diet of lactating bat diets (Figure 3D; Supplementary Table 3). At the genus level, the proportion of Rhipidia (p=0.000009), Symplecta (p=0.0002), Dicranomyia (p=0.02), and Dendrolimus (p=0.00003) were significantly increased and unclassified_f__Tipulidae (p=0.012), Amara (p=0.12), Hylaea_f__Geometridae (p=0.043; Figure 3E; Supplementary Table 3).

FIGURE 3

Figure 3. (A,B) Composition at the order level of insect species in the diet: (A) Lactating bats and (B) Non-lactating bats. (C–E) Analysis of insect species composition differences: (C) Analysis of order-level differences, (D) Analysis of family-level differences, and (E) Analysis of genus-level differences (^*0.01<p≤0.05, ^**0.001<p≤0.01, and ^***p≤0.001).

Gut Microbiota

After quality inspection, 730,728 effective sequences were obtained (average of 30,447 reads/sample). We identified 995 OTUs from reads, 276 of which were only found in the gut microbiota of lactating bats and 380 of which were only found in the gut microbiota of non-lactating bats (Supplementary Table 1).

There was no significant difference in the alpha diversity of gut microbiotas between lactating bats and non-lactating bats (Wilcoxon rank-sum test, Chao1 index, p=0.62, Q=0.98, Shannon index, p=0.98, Q=0.98; Figures 1C,D; Supplementary Figure 1B). For the comparison of individual bat microbiotas, the PCoA results based on unweighted UniFrac metric data showed that individuals from the lactating and non-lactating groups clustered separately. The results of ADONIS showed significant differences between the two groups (ADONIS: R²=0.1387, p=0.001; Figure 2B). But weighted UniFrac measurements showed no difference between the lactating and non-lactating groups (ADONIS: R²=0.0755, p=0.137; Supplementary Figure 2).

Three major bacterial phyla with relative abundance greater than 1% were identified in Asian particolored bat microbiotas, and most sequences were classified as Firmicutes (93.14%), followed by Proteobacteria (3.32%) and Actinobacteria (3.02%; Figure 4A).

FIGURE 4

Figure 4. (A–C) Relative abundance of major taxa in microbiotas of the bat gut: (A) Phylum level, (B) Family level, and (C) Genus level. (D–F) Analysis of gut microbial community differences: (D) Analysis of phylum-level differences, (E) Analysis of family-level differences, and (F) Analysis of genus-level differences (^*0.01<p≤0.05, ^**0.001<p≤0.01, and ^***p≤0.001).

In terms of microbial composition, the relative abundances of some dominant microbial groups at the level of phylum, family, and genus changed significantly. At the phylum level, compared with the gut microbiotas of non-lactating bats, Actinobacteria significantly increased (p=0.004) and Firmicutes decreased (p=0.01) in lactating bats (Figures 4A,D; Supplementary Figure 3; Supplementary Table 3). At the family level, Corynebacteriaceae (p=0.0001), Staphylococcaceae (p=0.0001), Mycoplasmataceae (p=0.001), and Rhizobiaceae (p=0.02) significantly increased and Lactobacillaceae (p=0.001) significantly decreased in lactating bats (Figures 4B,E; Supplementary Figure 3; Supplementary Table 3). At the genus level, Mycoplasma (p=0.001), Corynebacterium (p=0.05), Jeotgalicoccus (p=0.0004), unclassified_f__Corynebacteriaceae (p=0.0001) significantly increased and Weissella (p=0.03) and Lactobacillus (p=0.001) significantly decreased in lactating bats (Figures 4C,F; Supplementary Figure 3; Supplementary Table 3).

Correlation Between Diet and Gut Microbiota

Because effects of dietary composition on gut microbiome community may vary with taxonomic level, we summarized gut microbiome community and dietary datasets at multiple taxonomic levels, and performed correlation score for each rank. CCA suggests that there was a significant correlation between diet and gut microbiota in lactating and non-lactating bats at OTU and family levels (p<0.05; Figure 5; Supplementary Table 4). At the genus level, diet was found to be only correlated with the OTU level of gut microbiome (p=0.004), but no significant correlations were observed between genus and family level of gut microbiome (p>0.05; Figure 5; Supplementary Table 4). The Procrustes analysis showed that there was a strong correlation between the diet and gut microbiome of lactating bats compared with non-lactating bats (p<0.05; Supplementary Table 5).

FIGURE 5

Figure 5. Correlations between the gut microbial datasets and the dietary datasets summarized at different taxonomic levels.

Discussion

In this study, we examined the differences in diets and gut microbiotas between lactating and non-lactating bats V. sinensis and explored the correlations between the bats’ diets and the microbial composition of their guts during lactation and non-lactation periods. Our main findings are as follows: (1) Compared with non-lactating bats, the composition of lactating bats’ diets differed significantly, which suggests that lactating bats are selective in their diet. (2) Diversity of the gut microbiome did not differ between lactating and non-lactating bats, but relative abundance of gut microbes did vary. (3) Changes in the abundance of certain groups of gut microbes were associated with some particular types of foods, which confirms the role of diet in shaping gut microbiota. The results supported our hypothesis and prediction 1 and 3, while results comparing gut microbial diversity did not support our hypothesis and prediction 2.

Diet Changes

Lactation is the most energetically costly activity during reproduction in mammalian females (Vanslambrouck et al., 2021). Bats are the only truly flying mammals and during lactation, in addition to the high energy expenditure required to produce milk, females have to carry their young while flying. This has been shown to require staggering energy expenditure (Kunz et al., 1995; McLean and Speakman, 1999). Our study suggests that to meet the high energy requirements of lactation, female V. sinensis have a significantly higher diversity of insect species in their diet than non-lactating females, and the proportion of specific groups in their diet changed significantly. For example, lactating bats relied heavily on Diptera and decreased their consumption of Orthoptera and Coleoptera. The consumption of Lasiocampidae and Limacodidae in lactating bats was also significantly higher than that in non-lactating bats. There are several reasons for these results. First, Diptera have a higher protein content than Lepidoptera (Churchward-Venne et al., 2017). Lactating bats increased their dietary diversity and chose high-protein foods to ensure their nutritional needs were met and that lactation functioned properly to provide nutrient factors that are essential for infant growth and development (Lee et al., 2016; Lee and Kelleher, 2016; Marangoni et al., 2016). Thus, it is evident that high-quality food selection, rather than bulk feeding, is the most important component of nutrient acquisition when nursing (Vesterinen et al., 2016). Second, OFT suggests that an animal will attempt to gain the greatest energetic benefit for the lowest energetic cost while foraging to maximize fitness. Coleoptera have large exoskeleton that take time to handle, so bats tend to eat less of them during lactation (Dufour and Sauther, 2002). Although Orthoptera has a higher protein content, most species are active during the day. Few nocturnal species are active on the ground, which makes them more difficult for bats to hunt (Whitby et al., 2020). Finally, lactating females consumed significantly more species of certain insect families of the order Lepidoptera than did non-lactating females. We speculated that these particular species may have higher quality proteins that are more effective in supplementing the nutritional and energetic needs of lactating females, but evidence for this is lacking. On the basis of these findings, the OFT is supported, such that during lactation, Asian particolored bats support their energy demands by eating more nutritious food, which increases their nutritional intake.

The Gut Microbiota Changes

The gut microbiota plays an important role in the nutrition and health of the host (Al-Asmakh et al., 2014; Gentile and Weir, 2018). Dietary components have the capability to modulate the composition of this biota (Wu et al., 2011; Barrett et al., 2018). Although, we found differences in dietary composition between lactating and non-lactating bats, we did not detect significant differences in gut microbial diversity. The results suggest that healthy mammals tend to stabilize the diversity of their gut microbes in adulthood (Albenberg and Wu, 2014). Our results are consistent with studies in a variety of primates that confirmed that the alpha diversity of the gut microbial community remained stable throughout lactation. Studies on the human gut microbiota have shown that markers of microbial stability, such as richness and diversity, are often used as indicators of gut health (Gentile and Weir, 2018). In a few cases, they confirmed that the diversity of the intestinal microbiota of mammalian females changed with reproductive status (Gaona et al., 2019; Sun et al., 2020). However, these studies did not emphasize the role of dietary changes during lactation in shaping the structural diversity of gut microbiota (Phillips et al., 2012). On the contrary, they suggested that the diversity of gut microbiota in lactating females may be more susceptible to other factors besides diet, such as hormones and metabolism (Liu et al., 2019; Coleman et al., 2021).

Although, the overall diversity of gut microbes was unrelated to diet, we found that changes in diet had influences on specific, related microbiota. In human studies, it has been confirmed that dietary fat consumption is positively correlated with abundant Actinobacteria and negatively correlated with abundant Firmicutes (Abecia et al., 2007; Wu et al., 2011; Jin et al., 2021). Compared with non-lactating bats, the lactating bats had a higher relative abundance of Actinobacteria and lower relative abundance of Firmicutes, which would be useful for consuming more fat and reducing fat storage, thus improving the survival rate of offspring (Tarnaud, 2006; Hamady et al., 2008). Lactobacillus and Weissella mainly degrade macromolecules with complex structures to promote gut absorption (Famularo et al., 2005; Safika et al., 2019). Because of the reduced consumption of insects with large bodies and hard shells during lactation, it is speculated that the significant decrease in the relative abundance of Lactobacillus and Weissella in the intestines of V. sinensis during lactation is related to the decrease in the proportion of Coleoptera in the diet. Lactobacillus and Weissella participates in the production of short-chain fatty acids and maintains the balance of the gut and maintaining normal host health (Markowiak-Kopec and Slizewska, 2020). The relative abundance of Staphylococcaceae increased in lactating bats, which increased microflora nutrient acquisition by improving the hydrolysis of indigestible proteins and polysaccharides in the diets, and improving the energy efficiency of females to breastfeed their infants (Collado et al., 2008; Hunt et al., 2012; Rowland et al., 2018). In conclusion, by analyzing the correlation between diet and the gut microbiota at different classification levels, we confirmed that diet does have a regulatory effect on gut microbiota of lactating bats. Although there is no correlation between insect species at the genus level and gut microbes, this may be related to some unclassified species, but diet is only one of the factors regulating gut microbiota during lactation.

Mammals have evolved many behavioral and physiological strategies to meet the energy and nutrient requirements of their altered reproductive state. These strategies include eating longer, eating faster, eating more nutritious food, and enhancing metabolic efficiency (Bell and Bauman, 1997), as well as using seasonal reproduction (Bronson, 2009). In this study, we found that during lactation, Asian particolored bats did eat more nutritious food by consuming a diet that yields higher levels of protein. In addition, the abundance of specific microbes in their gut changed to help them adapt to changes to the composition of their diet. Because of limitations of sample size and the variables we considered, we were not able to fully reveal other behavioral and physiological changes caused by changes in mammalian reproductive status. Future studies should incorporate more factors to reveal the characteristics of adaptive evolution in different reproductive states of mammals, especially the adaptive mechanism of lactation, which is the reproductive period requiring the highest energy investments.

Data Availability Statement

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 in the article/Supplementary Material.

Ethics Statement

The animal study was reviewed and approved by the Laboratory Animal Welfare and Ethics Committee of Jilin Agricultural University.

Author Contributions

JL and HW contributed to the design of the study. WY, and YC collected the feces samples in all period. JF and HW contributed to the writing of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 32071492 and 32171489) and the Jilin Provincial Natural Science Foundation (Grant No. 20200201186JC).

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.

Publisher’s Note

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.

Acknowledgments

We are grateful to Tinglei Jiang and Lei Feng for their great contributions to samples collection. Sequencing service was provided by Ji Mei Biotechnology Co., Ltd., Shanghai, China. We thank Elizabeth Kelly from Liwen Bianji (Edanz; www.liwenbianji.cn/) for editing the English text of a draft of this manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.735122/full#supplementary-material

Footnotes

1. ^http://www.boldsystems.org/

2. ^http://www.ncbi.nlm.nih.gov/GenBank

References

Abecia, L., Fondevila, M., Balcells, J., and McEwanz, N. R. (2007). The effect of lactating rabbit does on the development of the caecal microbial community in the pups they nurture. J. Appl. Microbiol. 103, 557–564. doi: 10.1111/j.1365-2672.2007.03277.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Aiello, L. C., and Wells, J. C. K. (2002). Energetics and the evolution of the genusHOMO. Annu. Rev. Anthropol. 31, 323–338. doi: 10.1146/annurev.anthro.31.040402.085403

CrossRef Full Text | Google Scholar

Al-Asmakh, M., Stukenborg, J. B., Reda, A., Anuar, F., Strand, M. L., Hedin, L., et al. (2014). The gut microbiota and developmental programming of the testis in mice. PLoS One 9:e103809. doi: 10.1371/journal.pone.0103809

PubMed Abstract | CrossRef Full Text | Google Scholar

Albenberg, L. G., and Wu, G. D. (2014). Diet and the intestinal microbiome: associations, functions, and implications for health and disease. Gastroenterology 146, 1564–1572. doi: 10.1053/j.gastro.2014.01.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Araujo, M. S., Bolnick, D. I., and Layman, C. A. (2011). The ecological causes of individual specialisation. Ecol. Lett. 14, 948–958. doi: 10.1111/j.1461-0248.2011.01662.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrett, H. L., Gomez-Arango, L. F., Wilkinson, S. A., McIntyre, H. D., Callaway, L. K., Morrison, M., et al. (2018). A vegetarian diet is a major determinant of gut microbiota composition in early pregnancy. Nutrients 10:890. doi: 10.3390/nu10070890

PubMed Abstract | CrossRef Full Text | Google Scholar

Bell, A. W., and Bauman, D. E. (1997). Adaptations of glucose metabolism during pregnancy and lactation. J. Mammary Gland Biol. Neoplasia 2, 265–278. doi: 10.1023/A:1026336505343

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Bronson, F. (2009). Climate change and seasonal reproduction in mammals. Philos. Trans. R. Soc. B Biol. Sci. 364, 3331–3340. doi: 10.1098/rstb.2009.0140

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Churchward-Venne, T. A., Pinckaers, P. J. M., van Loon, J. J. A., and van Loon, L. J. C. (2017). Consideration of insects as a source of dietary protein for human consumption. Nutr. Rev. 75, 1035–1045. doi: 10.1093/nutrit/nux057

PubMed Abstract | CrossRef Full Text | Google Scholar

Coleman, D. N., Alharthi, A. S., Liang, Y., Lopes, M. G., Lopreiato, V., Vailati-Riboni, M., et al. (2021). Multifaceted role of one-carbon metabolism on immunometabolic control and growth during pregnancy, lactation and the neonatal period in dairy cattle. J. Anim. Sci. Biotechnol. 12:27. doi: 10.1186/s40104-021-00547-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Collado, M. C., Isolauri, E., Laitinen, K., and Salminen, S. (2008). Distinct composition of gut microbiota during pregnancy in overweight and normal-weight women. Am. J. Clin. Nutr. 88, 894–899. doi: 10.1093/ajcn/88.4.894

PubMed Abstract | CrossRef Full Text | Google Scholar

Denryter, K., Cook, R. C., Cook, J. G., Parker, K. L., and Gillingham, M. P. (2020). State-dependent foraging by caribou with different nutritional requirements. J. Mammal. 101, 544–557. doi: 10.1093/jmammal/gyaa003

PubMed Abstract | CrossRef Full Text | Google Scholar

Dufour, D. L., and Sauther, M. L. (2002). Comparative and evolutionary dimensions of the energetics of human pregnancy and lactation. Am. J. Hum. Biol. 14, 584–602. doi: 10.1002/ajhb.10071

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461

PubMed Abstract | CrossRef Full Text | Google Scholar

Famularo, G., De Simone, C., Pandey, V., Sahu, A. R., and Minisola, G. (2005). Probiotic lactobacilli: an innovative tool to correct the malabsorption syndrome of vegetarians? Med. Hypotheses 65, 1132–1135. doi: 10.1016/j.mehy.2004.09.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaona, O., Gomez-Acata, E. S., Cerqueda-Garcia, D., Neri-Barrios, C. X., and Falcon, L. I. (2019). Fecal microbiota of different reproductive stages of the central population of the lesser-long nosed bat, Leptonycteris yerbabuenae. PLoS One 14:e0219982. doi: 10.1371/journal.pone.0219982

PubMed Abstract | CrossRef Full Text | Google Scholar

Gentile, C. L., and Weir, T. L. (2018). The gut microbiota at the intersection of diet and human health. Science 362, 776–780. doi: 10.1126/science.aau5812

PubMed Abstract | CrossRef Full Text | Google Scholar

Haarsma, A.-J. (2008). Manual for assessment of reproductive status, age and health in European Vespertilionid bats. Electronic publication, version 2.

Google Scholar

Hamady, M., Walker, J. J., Harris, J. K., Gold, N. J., and Knight, R. (2008). Error-correcting barcoded primers for pyrosequencing hundreds of samples in multiplex. Nat. Methods 5, 235–237. doi: 10.1038/nmeth.1184

PubMed Abstract | CrossRef Full Text | Google Scholar

Herrera, E. R. T., and Heymann, E. W. (2004). Does mom need more protein? Preliminary observations on differences in diet composition in a pair of red titi monkeys (Callicebus cupreus). Folia Primatol. 75, 150–153. doi: 10.1159/000078304

PubMed Abstract | CrossRef Full Text | Google Scholar

Hunt, K. M., Preuss, J., Nissan, C., Davlin, C. A., Williams, J. E., Shafii, B., et al. (2012). Human milk oligosaccharides promote the growth of staphylococci. Appl. Environ. Microbiol. 78, 4763–4770. doi: 10.1128/AEM.00477-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Jari Oksanen, F., Blanchet, G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2019). vegan: Community Ecology. Package. R Package Version 2.5-6.

Google Scholar

Jin, L., Huang, Y., Yang, S., Wu, D., Li, C., Deng, W., et al. (2021). Diet, habitat environment and lifestyle conversion affect the gut microbiomes of giant pandas. Sci. Total Environ. 770:145316. doi: 10.1016/j.scitotenv.2021.145316

PubMed Abstract | CrossRef Full Text | Google Scholar

Kunz, T. (1987). “Post-natal growth and energetics of suckling bats,” in Recent Advances in the Study of Bats. Cambridge: Cambridge University Press, 395–420.

Google Scholar

Kunz, T. H., Whitaker, J., and Wadanoli, M. (1995). Dietary energetics of the insectivorous Mexican free-tailed bat (Tadarida brasiliensis) during pregnancy and lactation. Oecologia 101, 407–415. doi: 10.1007/BF00329419

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S., and Kelleher, S. L. (2016). Biological underpinnings of breastfeeding challenges: the role of genetics, diet, and environment on lactation physiology. Am. J. Physiol. Endocrinol. Metab. 311, E405–E422. doi: 10.1152/ajpendo.00495.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J., Moraes-Vieira, P. M., Castoldi, A., Aryal, P., Yee, E. U., Vickers, C., et al. (2016). Branched fatty acid esters of hydroxy fatty acids (FAHFAs) protect against colitis by regulating gut innate and adaptive immune responses. J. Biol. Chem. 291, 22207–22217. doi: 10.1074/jbc.M115.703835

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Hou, C., Li, N., Zhang, X., Zhang, G., Yang, F., et al. (2019). Microbial and metabolic alterations in gut microbiota of sows during pregnancy and lactation. FASEB J. 33, 4490–4501. doi: 10.1096/fj.201801221RR

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J., Liang, B., Zhang, S., and Metzner, W. (2007). Dietary composition and echolocation call design of three sympatric insectivorous bat species from China. Ecol. Res. 23, 113–119. doi: 10.1007/s11284-007-0344-5

CrossRef Full Text | Google Scholar

Magoč, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507

PubMed Abstract | CrossRef Full Text | Google Scholar

Mallott, E. K., Borries, C., Koenig, A., Amato, K. R., and Lu, A. (2020). Reproductive hormones mediate changes in the gut microbiome during pregnancy and lactation in Phayre's leaf monkeys. Sci. Rep. 10:9961. doi: 10.1038/s41598-020-66865-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Marangoni, F., Cetin, I., Verduci, E., Canzone, G., Giovannini, M., Scollo, P., et al. (2016). Maternal diet and nutrient requirements in pregnancy and breastfeeding. An italian consensus document. Nutrients 8:629. doi: 10.3390/nu8100629

PubMed Abstract | CrossRef Full Text | Google Scholar

Markowiak-Kopec, P., and Slizewska, K. (2020). The effect of probiotics on the production of short-chain fatty acids by human intestinal microbiome. Nutrients 12:1107. doi: 10.3390/nu12041107

PubMed Abstract | CrossRef Full Text | Google Scholar

McCabe, G. M., and Fedigan, L. M. (2007). Effects of reproductive status on energy intake, ingestion rates, and dietary composition of female Cebus capucinus at Santa Rosa, Costa Rica. Int. J. Primatol. 28, 837–851. doi: 10.1007/s10764-007-9159-z

CrossRef Full Text | Google Scholar

McLean, J., and Speakman, J. (1999). Energy budgets of lactating and non-reproductive brown long-eared bats (Plecotus auritus) suggest females use compensation in lactation. Funct. Ecol. 13, 360–372. doi: 10.1046/j.1365-2435.1999.00321.x

CrossRef Full Text | Google Scholar

Muegge, B. D., Kuczynski, J., Knights, D., Clemente, J. C., Gonzalez, A., Fontana, L., et al. (2011). Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. Science 332, 970–974. doi: 10.1126/science.1198719

PubMed Abstract | CrossRef Full Text | Google Scholar

Oftedal, O. T. (1991). The nutritional consequences of foraging in primates: the relationship of nutrient intakes to nutrient requirements. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 334, 161–170. doi: 10.1098/rstb.1991.0105

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillips, C. D., Hanson, J., Wilkinson, J. E., Koenig, L., Rees, E., Webala, P., et al. (2017). Microbiome structural and functional interactions across host dietary niche space. Integr. Comp. Biol. 57, 743–755. doi: 10.1093/icb/icx011

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Racey, P. (2009). “Reproductive assessment of bats,” in Ecological and Behavioural Methods in the Study of Bats. eds. T. H. Kunz and S. Parsons (Baltimore, MD: The Johns Hopkins University Press), 249–264.

Google Scholar

Ramette, A. (2007). Multivariate analyses in microbial ecology. FEMS Microbiol. Ecol. 62, 142–160. doi: 10.1111/j.1574-6941.2007.00375.x

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2019). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna: R Core Team.

Google Scholar

Rowland, I., Gibson, G., Heinken, A., Scott, K., Swann, J., Thiele, I., et al. (2018). Gut microbiota functions: metabolism of nutrients and other food components. Eur. J. Nutr. 57, 1–24. doi: 10.1007/s00394-017-1445-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Rydell, J., and Baagøe, H. J. (1994). Vespertilio murinus. Mamm. Species 146, 1–6. doi: 10.2307/3504265

CrossRef Full Text | Google Scholar

Safi, K., König, B., and Kerth, G. (2007). Sex differences in population genetics, home range size and habitat use of the parti-colored bat (Vespertilio murinus, Linnaeus 1758) in Switzerland and their consequences for conservation. Biol. Conserv. 137, 28–36. doi: 10.1016/j.biocon.2007.01.011

CrossRef Full Text | Google Scholar

Safika, S., Wardinal, W., Ismail, Y. S., Nisa, K., and Sari, W. N. (2019). Weissella, a novel lactic acid bacteria isolated from wild Sumatran orangutans (Pongo abelii). Vet. World 12:1060. doi: 10.14202/vetworld.2019.1060-1065

PubMed Abstract | CrossRef Full Text | Google Scholar

Schoener, T. W. (1971). Theory of feeding strategies. Annu. Rev. Ecol. Syst. 2, 369–404. doi: 10.1146/annurev.es.02.110171.002101

CrossRef Full Text | Google Scholar

Sigsgaard, E. E., Olsen, K., Hansen, M. D. D., Hansen, O. L. P., Hoye, T. T., Svenning, J. C., et al. (2020). Environmental DNA metabarcoding of cow dung reveals taxonomic and functional diversity of invertebrate assemblages. Mol. Ecol. 30, 3374–3389. doi: 10.1111/mec.15734

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, B., Xu, X., Xia, Y., Cheng, Y., Mao, S., Xiang, X., et al. (2020). Variation of gut microbiome in free-ranging female tibetan macaques (Macaca thibetana) across different reproductive states. Animals 11:39. doi: 10.3390/ani11010039

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarnaud, L. (2006). Feeding behavior of lactating brown lemur females (Eulemur fulvus) in Mayotte: influence of infant age and plant phenology. Am. J. Primatol. 68, 966–977. doi: 10.1002/ajp.20288

PubMed Abstract | CrossRef Full Text | Google Scholar

Taylor, S. E., Klein, L. C., Lewis, B. P., Gruenewald, T. L., Gurung, R. A. R., and Updegraff, J. A. (2000). Biobehavioral responses to stress in females: tend-and-befriend, not fight-or-flight. Psychol. Rev. 107, 411–429. doi: 10.1037/0033-295X.107.3.411

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, M. E., and Wrangham, R. W. (2008). Diet and reproductive function in wild female chimpanzees (Pan troglodytes schweinfurthii) at Kibale National Park, Uganda. Am. J. Phys. Anthropol. 135, 171–181. doi: 10.1002/ajpa.20718

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanslambrouck, K., de Kok, B., Toe, L. C., De Cock, N., Ouedraogo, M., Dailey-Chwalibóg, T., et al. (2021). Effect of balanced energy-protein supplementation during pregnancy and lactation on birth outcomes and infant growth in rural Burkina Faso: study protocol for a randomised controlled trial. BMJ Open 11:e038393. doi: 10.1136/bmjopen-2020-038393

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasey, N. (2000). Niche separation in Varecia variegata rubra and Eulemur fulvus albifrons: I. Interspecific patterns. Am. J. Phys. Anthropol. 112, 411–431. doi: 10.1002/1096-8644(200007)112:3<411::AID-AJPA10>3.0.CO;2-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasey, N. (2002). Niche separation in Varecia variegata rubra and Eulemur fulvus albifrons: II. Intraspecific patterns. Am. J. Phys. Anthropol. 118, 169–183. doi: 10.1002/ajpa.10054

PubMed Abstract | CrossRef Full Text | Google Scholar

Vesterinen, E. J., Ruokolainen, L., Wahlberg, N., Pena, C., Roslin, T., Laine, V. N., et al. (2016). What you need is what you eat? Prey selection by the bat Myotis daubentonii. Mol. Ecol. 25, 1581–1594. doi: 10.1111/mec.13564

PubMed Abstract | CrossRef Full Text | Google Scholar

Whitby, M. D., Kieran, T. J., Glenn, T. C., and Allen, C. (2020). Agricultural pests consumed by common bat species in the United States corn belt: the importance of DNA primer choice. Agric. Ecosyst. Environ. 303:107105. doi: 10.1016/j.agee.2020.107105

CrossRef Full Text | Google Scholar

Wichert, B., Schade, L., Gebert, S., Bucher, B., Zottmaier, B., Wenk, C., et al. (2009). Energy and protein needs of cats for maintenance, gestation and lactation. J. Feline Med. Surg. 11, 808–815. doi: 10.1016/j.jfms.2009.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, G. D., Chen, J., Hoffmann, C., Bittinger, K., Chen, Y.-Y., Keilbaugh, S. A., et al. (2011). Linking long-term dietary patterns with gut microbial enterotypes. Science 334, 105–108. doi: 10.1126/science.1208344

PubMed Abstract | CrossRef Full Text | Google Scholar