Jingjing Wang
Tao Lang
Jian Shen
Juanjuan Dai
Ling Tian1*- 1Shanghai Key Laboratory for Pancreatic Diseases, Institute of Translational Medicine, Shanghai General Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- 2Ministry of Education Key Laboratory for Systems Biomedicine, Shanghai Centre for Systems Biomedicine, Shanghai Jiao Tong University, Shanghai, China
Previous studies revealed that there existed great individual variations of gut microbiota in mice, and the gut bacteria of mice were changed with the occurrence and development of diseases. To identify the core gut bacteria in healthy mice and explore their relationships with the host phenotypes would help to understand the underlying mechanisms. In this study, we identified 37 genus-level core bacteria from feces of 101 healthy mice with different ages, sexes, and mouse strains in three previous studies. They collectively represented nearly half of the total sequences, and predominantly included carbohydrate- and amino acids-metabolizing bacteria and immunomodulatory bacteria. Among them, Anaerostipes indwelt the gut of all healthy mice. Co-abundance analysis showed that these core genera were clustered into five groups (Group C1–C5), which were ecologically related. For example, the abundances of Group C2 including probiotics Bifidobacterium and Lactobacillus slightly positively correlated with those of Group C1. Principal component analysis (PCA) and multivariate analysis of variance test revealed that these core gut genera were distinguished with age and sex, and also associated with their health/disease state. Linear discriminant analysis effect size (LEfSe) method showed that bacteria in Group C1 and C2/C3 increased with the age in infancy and early adulthood, and were more abundant in female mice than in male ones. The metabolic syndrome (MS) induced by high fat diet (HFD) and accelerated postnatal growth would decrease Group C2 genera, whereas probiotics intervention would reverse HFD-induced reduction of Group C2. Spearman correlation analysis indicated that the principal components based on the abundance of the 37 core genera were significantly correlated with host characteristic parameters of MS. These results demonstrated that the 37 core genera in five co-abundance groups from healthy mice were related to host phenotypes. It was indicated that these prevalent gut bacterial genera could be representative of the healthy gut microbiome in gnotobiotic animal models, and might also be candidates of probiotics and fecal microbiota transplantation.
Introduction
Gut microbiota is a highly complex ecosystem, with thousands of microbial species and great individual variations. It is known that gut microbiota modulates host immune system development and whole-body metabolism (Kau et al., 2011; Tremaroli and Backhed, 2012). The dysbiosis of gut microbiota will induce immunological and/or metabolic diseases, such as inflammatory bowel disease, colorectal cancer, and type 2 diabetes (Wong et al., 2017; Schirmer et al., 2018; Zhao et al., 2018). Besides the case–control gut microbiota study, the healthy core microbiota study focusing on the stable and permanent members of the community in normal and healthy populations has drawn considerable attention. It was reported that some bacterial phylogroups, including Clostridium XIVa, Faecalibacterium, Ruminococcus, Bacteroides, Alistipes, Parabacteroides, Roseburia, Lachnospiraceae, Sporobacter, Dorea, Clostridium, Eubacterium, Collinsella, Coprococcus, Subdoligranulum, Streptococcus, Holdermania, Butyrivibrio, Anaerotruncus, Enterococcus, Blautia, Bifidobacterium, Anaerostipes, Lactobacillus, Phascolarctobacterium, Prevotella, Odoribacter, Veillonella, Eggerthella, Fusobacterium, Akkermansia, Escherichia, Bilophila, Oscillibacter, Desulfovibrio, Klebsiella, or Shigella, accounted for 50–100% of the gut microbiota of the healthy individuals (Tap et al., 2009; Turnbaugh et al., 2009; Qin et al., 2010; Claesson et al., 2011; Sekelja et al., 2011; Xiao et al., 2015; Zhang J. et al., 2015; Cheng et al., 2016; Maier et al., 2018). No matter how the healthy core bacteria are defined and how different the identification results are in the population studies, it is certain that the prevalent and dominant core bacteria are pivotal to host gut homeostasis and health. So it is important to discover a comprehensive core microbiota profile for defining a healthy gut microbiota and guiding their intervention of host health.
Nevertheless, the ecological relationship of the members of core bacteria is not well understood. In a cohort study with 314 samples of healthy young Chinese, eight in nine core genera, including Blautia, Clostridium, Ruminococcus, Faecalibacterium, Subdoligranulum, Roseburia, Coprococcus, Bacteroides, were significantly positively correlated with each other, yet the core genus Phascolarctobacterium was negatively correlated with the other eight core genera (Zhang J. et al., 2015). In another study of dietary fibers treatment to genetic and simple obesity in children, species from Eubacterium, Faecalibacterium, Roseburia, Clostridium, Alistipes, Subdoligranulum, Bacteroides, Parabacteroides, Flavonifractor were positively correlated with each other, while Bifidobacterium including B. pseudocatenulatum and B. breve, showed negative correlation with other bacteria (Zhang C. et al., 2015). These results indicated that most of the prevalent bacteria were co-abundant and positively correlated. However, a healthy gut microbiota ecosystem should be full of cooperation and competition, which could keep the system stable (Zengler and Zaramela, 2018). Therefore, it is necessary to explore the ecology of core microbiota in healthy state including more bacteria without any selection pressures such as dietary intervention.
Thereinto, another crucial question is whether host phenotype is linked to the core gut bacteria or not. Some previous studies found that the core bacteria showed the ethnicity- and age-associated differences (Zhang J. et al., 2015; Cheng et al., 2016), yet some other ones did not involve the association between gut bacteria and host phenotype. More importantly, healthy gut microbiota are distinct from those of diseased hosts (Round and Mazmanian, 2009; Sekirov et al., 2010), so it should be made clear whether the discrepancy are mainly reflected in the stable core bacteria or the variable ones.
Investigation of the prevalent bacteria in mice feces would help to study human core gut bacteria. Firstly, mouse and human gut microbiota show similarity at the genus level even though the proportions are different (Xiao et al., 2015; Hugenholtz and de Vos, 2018). Secondly, the mouse gut microbiota is functionally similar to its human counterpart (Xiao et al., 2015). Meanwhile, genotype, diets, lifestyles and drugs, which are key factors to affect gut microbiota (Franks, 2011; Claesson et al., 2012; Conlon and Bird, 2014; Maier et al., 2018), can be strictly controlled in the mouse model so as to down out any outside noises and find out how the gut bacteria group work together. Moreover, many conclusions from mice studies are validated on subsequent clinical studies (Sivan et al., 2015; Routy et al., 2018).
Herein, we investigated the full profile of gut core bacteria at the genus level in healthy mouse. By using co-abundance analysis, we divided the gut core bacteria in groups, and further investigated the ecological relationships between groups. With multivariate statistical analysis, we identified the specific core bacteria associated with age, sex and health/disease state, and explored their role in host health. These results provided a comprehensive view on core gut bacteria and host health in mice. These prevalent core gut bacteria could be used as individual or cocktail candidates for probiotics and fecal microbiota transplantation.
Materials and Methods
Datasets
To search for the possible core bacteria in the mouse gut, we analyzed 16S rRNA gene sequences of mouse gut microbiota from three previously published datasets with accession numbers of SRP020353, SRP064646, and SRP076247 (Wang J. et al., 2015; Wang J. et al., 2016; Wang J.J. et al., 2016).
In the first dataset with the ID “PRO” (Wang J. et al., 2015), we collected 301568 fecal bacteria 16S rRNA V3 region sequences from 80 fecal samples of mice with or without metabolic syndrome (MS) and probiotics intervention (Table 1 and Supplementary Table S1). Briefly, 40 male specific pathogen-free (SPF) C57BL/6J mice at the age of 12 weeks were divided into five groups with eight mice a group. One group was fed on normal chow diet (NC), as the control group. The other four groups were fed on high fat diet (HFD) for 12 weeks to induce MS. Among them, three HFD groups were treated by three probiotics strains, respectively. At the end of 12 weeks, we examined the body weight, measured the fasting blood glucose and fasting insulin, carried out the area under the curve (AUC) of oral glucose tolerance tests (OGTT), as well as tested serum adiponectin, serum lipopolysaccharide-binding protein (LBP), and crown-like structures (CLS), cecal short-chain fatty acids (including acetate, propionate, butyrate, isobutyrate, valerate, isovalerate), and recorded the epididymal adipose tissue (eAT) adipocyte size. We also analyzed the mRNA levels of adiponectin, leptin, tumor necrosis factor-α (TNF-α), matrix metalloproteinase-12 (MMP-12), monocyte chemotactic protein-1 (MCP-1), CD11c in eAT, TNF-α in liver and jejunum. These host parameters were also used in our study.
Table 1. The fecal samples of healthy and disease mice from three datasets used in this study.
In the second dataset CUG (Wang J. et al., 2016), we got 432428 fecal bacteria 16S rRNA V3 region sequences from 106 fecal samples of male SPF ICR mice at three time points (4, 12, and 24 weeks of age) (Table 1 and Supplementary Table S1). Three groups were divided in this dataset, and each group contained 12 mice. Besides of the control group, accelerated postnatal growth (APG) groups with normal birth weight (NBW+APG) and low birth weight (LBW+APG) were distinguished. The two APG groups (NBW+APG and LBW+APG) developed MS at the age of 4 weeks. The physiological index data of these grouped mice was also analyzed in this study, including body weight, weight of eAT, perirenal adipose tissue (pAT), subcutaneous adipose tissue (sAT), fasting blood glucose, fasting insulin, and cecal short-chain fatty acids levels.
The third dataset FMT (Wang J.J. et al., 2016) included 698564 fecal bacteria V3–V4 region sequences from 19 fecal samples of gnotobiotic mice (Table 1 and Supplementary Table S1). Ten male and 9 female germ-free C57BL/6J mice at the age of 7 weeks were involved. These mice were inoculated gut microbiota from a healthy man 1 week before fecal sample collection.
Taken together, there were 205 mice fecal samples in the three datasets (Supplementary Table S1). Among them, samples of 48, 34, and 19 healthy mice were in the dataset PRO, CUG, and FMT, respectively (totally 101 healthy mice) (Table 1).
16S rRNA Gene Sequence Processing
Raw data were clustered into operational taxonomic units (OTUs) as before (Wang J. et al., 2015; Wang J. et al., 2016; Wang J.J. et al., 2016), and showed in Supplementary Table S1. For each dataset, the OTUs which were identified as chimeras and singleton, and which contained fewer than 10 reads, and which present in fewer than 1% of samples were removed (Duvallet et al., 2017). The representative sequences of each OTU were subjected to RDP classifier (version 2.11) to determine the taxonomy (RDP 16S rRNA training set 16). The confidence threshold was 50% for the data of 454 pyrosequencing of V3 region, whereas the confidence cut off was 80% for the data of Illumina Miseq sequencing of V3–V4 region. The OTUs to genus was then collapsed, and the relative abundance of each genus in every fecal sample was calculated. The followed statistical analyses were all based on genus-level data to minimize batch effects.
Co-abundance Analysis
The core genera were clustered into the core groups based on their abundance across all the 101 healthy mice samples using spearman correlation analysis. False discovery rate (FDR) was estimated using the Benjamini–Yekutieli method (Benjamini et al., 2001). It was considered significant when the absolute value of the correlation coefficient was more than 0.5 and the FDR value was less than 0.001.
Alignment-Based Phylogenetic Analysis
To identify the phylogenetic relationships of the five core groups, the phylogenetic tree was reconstructed based on the sequences of the most abundant OTUs (named as OTU01 in the figures and tables) belonging to the 37 core gut bacteria in each dataset. Sequence alignment and phylogenetic analysis were performed using Mega 7. Then the neighbor-joining tree was constructed with 1,000 of bootstrap replications.
Statistical Analysis
Principal component analysis (PCA) was performed to globally view the age-, sex-, and disease-associated gut microbiota structural segregation based on the Euclidean distances of the relative abundances of the 37 core genera. The statistical significance of PCA plots was assessed by multivariate analysis of variance (MANOVA) test in MATLAB R2011b (The MathWorks, Inc., Natick, MA, United States). The difference was considered significant when p-value of MANOVA was less than 0.05. Linear discriminant analysis effect size (LEfSe) were used to identify the differential core genera between different ages, sexes, and healthy status based on the relative abundance of 37 core genera. The OTUs were picked out when alpha value of the factorial Kruskal–Wallis test was <0.05 and the logarithmic LDA score was >2.0 (Segata et al., 2011). Spearman correlation analysis was used to evaluate the relationships between core microbiota based on the principle component scores of PCA and the individual host parameters of phenotype, gene expression, and metabolite. The correlations were considered significant when FDR was <0.1.
Results
Thirty-Seven Core Gut Bacteria Were Identified in Healthy Mice
We collected and re-analysis bacterial sequencing data from three previous studies to figure out the prevalent gut bacteria in the healthy mice. These data included 16S rRNA gene V3 or V3–V4 region sequences of fecal bacterium in mice with different ages, genders, strains, and healthy/disease state (Table 1 and Supplementary Table S1). We used OTUs clustered in the previous studies, and removed the OTUs potentially caused by sequencing mistakes (Supplementary Table S1). RDP classifier assigned the OTUs to 72, 72, and 48 genera in the database PRO, CUG, and FMT, respectively, and 101 genera in total (Supplementary Table S1 and Supplementary Figure S1). The 101 genera belonged to the following 10 phyla: Firmicutes (63 genera), Proteobacteria (12 genera), Bacteroidetes (11 genera), Actinobacteria (9 genera), Candidatus Saccharibacteria (1 genus), Deferribacteres (1 genus), Tenericutes (1 genus), Verrucomicrobia (1 genus), Fusobacteria (1 genus), Spirochaetes (1 genus) (Supplementary Figure S1). Among the 101 genera, 26 genera existed in all the three datasets, and 39 genera existed in two of them (Supplementary Figure S2).
We defined healthy core gut bacteria as those genera appeared in at least 50% of healthy mice samples regardless of their age, sex, and strain. In this study, there were 101 healthy mice samples in three datasets, which were from 40 male C57BL/6J mice of 12-week-old and 8 ones of 24-week-old in the database PRO, 11 male ICR mice of 4-week-old, 12 ones of 12-week-old and 11 ones of 24-week-old in the database CUG, 10 male C57BL/6J mice of 8-week-old, and 9 female ones of 8-week-old in the database FMT (Table 1). Across these 101 healthy mice, we identified 37 core genera. Among them, the most prevalence genera were the two, Anaerostipes (100% prevalence) and Parabacteroides (99.0%). The rest genera included Anaerotruncus (94.1%), Oscillibacter (93.1%), Clostridium XlVb (89.1%), Flavonifractor (88.1%), Bacteroides (85.1%), Barnesiella (81.2%), Alistipes (80.2%), Helicobacter (80.2%), Saccharibacteria_genera_incertae_sedis (80.2%), Prevotella (79.2%), Lachnoanaerobaculum (76.2%), Lactobacillus (76.2%), Intestinimonas (75.2%), Roseburia (75.2%), Alloprevotella (73.3%), Rikenella (73.3%), Enterorhabdus (73.3%), Erysipelotrichaceae_incertae_sedis (70.3%), Eggerthella (69.3%), Allobaculum (64.4%), Lachnospiracea_incertae_sedis (64.4%), Pseudoflavonifractor (63.4%), Bifidobacterium (62.4%), Marvinbryantia (60.4%), Mucispirillum (60.4%), Clostridium XIVa (59.4%), Blautia (59.4%), Anaerofilum (58.4%), Parasutterella (57.4%), Odoribacter (55.4%), Olsenella (53.5%), Turicibacter (52.5%), Gordonibacter (52.5%), Ruminococcus (51.5%), and Acetatifactor (50.5%) (Figures 1A,B and Supplementary Table S2).
Figure 1. Thirty-seven core genera were identified in more than 50% of the feces of 101 healthy mice. (A) The proportion of each genus in the total sequences of each dataset. (B) The heatmap of the relative abundance (log2 transformed) of the 37 core genera in 101 healthy mice. Rows corresponded to the 37 core genera in the feces of 101 healthy mice. These genera were presented in order of prevalence (100–50.5%). The taxonomies of the genera were shown on the right, followed by the prevalence, and colored according to phylum. Columns represented 101 healthy mice samples in the three datasets. The 101 mice were grouped based on their sex, age, and strain.
The core genera were also highly abundant, ranking 36.60%, 51.49%, and 60.9% of the total sequences in dataset PRO, CUG, and FMT, respectively, and 46.77% on average (Figure 1A and Supplementary Table S2). Especially, it showed the highest proportion of gut microbiota (60.9%) in the dataset FMT even though the least core genera (22 core genera, 59.46%) were found, which might due to gavage transplantation of human gut microbiota into mice via oral gavage (Figure 1A and Supplementary Figure S3). Besides, there were 10 genera of the average relative abundances were over 1%, including Lactobacillus (10.51%), Parabacteroides (6.63%), Barnesiella (4.85%), Anaerostipes (3.29%), Bacteroides (2.35%), Alistipes (1.52%), Oscillibacter (1.51%), Alloprevotella (1.38%), Turicibacter (1.38%), and Helicobacter (1.21%), respectively (Figure 1B and Supplementary Table S2).
The 37 Genera Were Divided Into Five Potential Core Groups
Spearman correlation analysis was further used to exhibit the relationships between any two co-abundance groups of the 37 core gut genera. The positively correlated co-abundance genera were considered as a core group. Accordingly, the 37 core genera were classified into 5 potential core groups (Figure 2A). Among them, Group C1 included 22 genera, including Anaerostipes, Anaerotruncus, Oscillibacter, Clostridium XlVb, Bacteroides, Barnesiella, Alistipes, Helicobacter, Saccharibacteria_genera_incertae_sedis, Prevotella, Lachnoanaerobaculum, Intestinimonas, Roseburia, Alloprevotella, Rikenella, Allobaculum, Lachnospiracea_incertae_sedis, Pseudoflavonifractor, Marvinbryantia, Mucispirillum, Odoribacter and Acetatifactor. Group C2 contained six genera, including Bifidobacterium, Olsenella, Lactobacillus, Enterorhabdus, Parasutterella, and Turicibacter. Group C3 covered five genera, including Parabacteroides, Flavonifractor, Clostridium XIVa, Blautia, and Anaerofilum. Group C4 comprised three genera, including Erysipelotrichaceae _incertae_sedis, Eggerthella, and Gordonibacter. Group C5 only constituted one genus, Ruminococcus (Figure 2A).
Figure 2. The 37 core genera were classified into five potential core groups by co-abundance analysis. (A) Interaction patterns of the 37 core genera across the 101 samples determined by Spearman correlation analysis. Rows and columns corresponded to the 37 core genera. These genera were with the names shown on the right and top, colored according to phylum, and grouped based on the co-abundance. Red and blue color denoted positive and negative association, respectively. The color intensity represented the degree of association. The black dots in the red/blue cells indicated the associations were significant (0.5 < | r| < 1, and false discovery rate < 0.001). In each core group, the genera were showed in order of prevalence. (B) Phylogenetic reconstruction of sequences to the 37 core genera. The most abundant OTUs (OTU01) of the 37 core genera in each dataset and their nearest neighbors were chosen to construct this neighbor-joining tree using Mega 7. The clades were colored according to the five core groups.
The five core groups were associated with each other (Figure 2A). Members in Group C1 was weakly positively correlated with those in Group C2. Especially, Allobaculum in Group C1 and Bifidobacterium in Group C2 were significantly positively correlated. Similarly, Group C3 showed positive association with Group C4 even if no statistical significance. Group C3 were significantly negatively correlated with Group C2, and Group C4 also showed significantly negative correlation with Group C1. However, Group C5 was the smallest group with only one genus, and was not correlated with any other groups.
To demonstrate the evolutionary relationship of the five core groups, we constructed the phylogenetic tree with the sequences of the most abundant OTUs belonging to the 37 core genera and their nearest neighbors (Figure 2B and Supplementary Tables S3, S4). Each OTU stayed together with its nearest neighbors (shown in Supplementary Table S4), and the OTUs from the three datasets in the same genus got together. Group C1 and Group C3 were widely distributed in the evolutionary tree, and most of them belonged to Firmicutes. Group C2 were clustered into two sub-groups: one belonged to Firmicutes, including Lactobacillus and Turicibacter, and the other included Bifidobacterium, Olsenella, Enterorhabdus (Actinobacteria), and Parasutterella (Proteobacteria), respectively. Group C4 gathered together. Group C5 were clustered with their nearest neighbors and Anaerofilum in Group C3.
The 37 Genera Could Differentiate Gut Microbiota With Age and Sex
We wondered whether these healthy core genera could separate gut microbiota with different ages and sexes. So PCA was performed based on the abundance of the core genera, and MANOVA was used to confirm the statistically significant (Figures 3A,C,G). Further, linear discriminant analysis (LDA) effect size (LEfSe) was applied to identify the specific core bacteria which were separated in the two group mice (Figures 3B,D–F,H).
Figure 3. The core bacteria showed significant age- and sex-associated structural segregation. (A,C,G) PCA score plot based on the core genus abundance of healthy mice in dataset PRO (A), CUG (C), and FMT (G). Each circular point represented the fecal core microbiota of a healthy mouse. P-values were calculated with MANOVA test, and shown on the top of PCA plots. (B,D–F,H) The different abundance of bacterial genus between the two groups were identified by LEfSe. It was significantly different when alpha value of the factorial Kruskal–Wallis test was <0.05 and the logarithmic LDA score was >2.0. The left histogram showed the LDA scores of genera differentially abundant between the two groups. The taxonomy was listed, followed by its core group. The right heatmap showed the relative abundance (log2 transformed) of the differential genera. n = 40 for male C57BL/6J mice (12 weeks), n = 8 for male C57BL/6J mice (24 weeks), n = 11 for male ICR mice (4 weeks), n = 12 for male ICR mice (12 weeks), n = 11 for male ICR mice (24 weeks), n = 10 for male C57BL/6J mice (8 weeks), and n = 9 for female C57BL/6J mice (8 weeks).
Principal component analysis showed that the core gut microbiota of 12-week-old male C57BL/6J mice were significantly different from those of 24-week-old ones (MANOVA, p = 8.97E-10) (Figure 3A). LEfSe LDA score bar and heatmap revealed that these abundance differences between the two group mice were increased in Group C1 (Bacteroides, Oscillibacter, Intestinimonas, Mucispirillum, Anaerotruncus), Group C2 (Olsenella, Bifidobacterium, Enterorhabdus), and Group C4 (Eggerthella), whilst decreased in Group C1 (Barnesiella, Prevotella, Lachnospiracea_incertae_sedis), Group C2 (Parasutterella), and Group C5 (Ruminococcus) (Figure 3B). Among the male ICR mice, PCA plots showed that the core gut bacteria in 4-week-old mice (infancy) were significantly separated with those of 12- and 24-week-old ones (adult) (MANOVA, p = 2.24E-07). As to the C57BL/6J mice, the samples of 12-week-old mice and 24-week-old ones were also significantly separated (MANOVA, p = 0.018) (Figure 3C). LEfSe results showed that the infant mice had more Group C1 and C3 than the adult mice, including Anaerostipes (C1), Prevotella (C1), Acetatifactor (C1), Alloprevotella (C1), Saccharibacteria_genera_incertae_sedis (C1), Alistipes (C1), Rikenella (C1), Oscillibacter (C1), Anaerotruncus (C1), Intestinimonas (C1), Clostridium XIVb (C1), Anaerofilum (C3), and Flavonifractor (C3), whilst less Group C2, including Turicibacter, Bifidobacterium, Enterorhabdus, Olsenella, Parasutterella (Figures 3D,E). In adult ICR mice, Group C1 (Saccharibacteria_genera_incertae_sedis, Prevotella, Alloprevotella, Clostridium XIVb, and Alistipes) and Group C5 (Ruminococcus) were increased, and Group C2 (Olsenella, Lactobacillus) were decreased, when the mice got older (Figure 3F). Taken together, these results indicated that the 37 core genera were age-associated.
The 37 core genera were also related to sex. PCA plots showed the core bacteria were significantly different between female and male C57BL/6J recipient mice transplanted by fecal microbiota from a healthy human (MANOVA, p = 0.049) (Figure 3G). LEfSe showed the female mice had more abundance of Group C1 (Anaerotruncus, Clostridium XIVb), Group C3 (Blautia, Flavonifractor, Anaerofilum) and Group C4 (Erysipelotrichaceae_incertae_sedis) genera than the male mice (Figure 3H).
The 37 Core Genera Could Distinguish Gut Microbiota From Health and Disease State
We also wondered whether these 37 core genera were correlated with health/disease states. PCA was performed based on the abundance of the 37 core genera in healthy and sick mice (metabolic syndrome model, MS for short), and MANOVA was used to test the statistical significance of the differences between them. Spearman correlation analysis between the three principal components (PCs) of PCA was further carried out to determine the relationship of core bacteria and host disease state based on the core bacteria and host disease parameters in the dataset PRO or CUG. Then, different core genera between the two groups were picked out by LEfSe.
In the dataset PRO, PCA plots showed HFD-induced MS group (shown in the figure as HFD) were significantly different from normal chow-fed control (shown in the figure as NC) (MANOVA, p = 4.91E-07), while probiotics treatment (shown in the figure as HFD+probiotics) could significantly reduce the HFD-induced harmful effects (MANOVA, p = 0.044) (Figure 4A). In addition, PC1 and PC3 was positively correlated with the protective parameters, including serum adiponectin/body weight, adiponectin mRNA in eAT, cecal acetate, and cecal butyrate, whilst negatively correlated with obesity and inflammation parameters, including body weight gain, AUC of OGTT, fasting blood glucose, serum LBP, eAT adipose size, eAT CLS, eAT MMP-12, eAT TNF-α mRNA, MCP-1 mRNA, CD11c mRNA, leptin mRNA, liver TNF-α mRNA, and jejunum TNF-α mRNA (Figure 4B). Then, LEfSe analysis indicated that HFD group had more Group C1 than NC (Roseburia, Clostridium XIVb, and Lachnospiracea_incertae_sedis), whilst less Group C2 (Bifidobacterium, Olsenella, and Turicibacter) and Group C1 (Barnesiella and Anaerostipes) (Figure 4C). However, probiotics treatment increased Group C2 including Lactobacillus, Olsenella, Bifidobacterium, and Group C1 including Allobaculum, whilst decreased Group C1 including Roseburia and Clostridium XIVb (Figure 4D). All these results suggested that these core genera were significantly related to HFD-induced MS and probiotics-mediated intervention.
Figure 4. Core bacteria significantly distinguished the gut microbial structure of healthy and disease mice. (A,E) PCA score plot was presented based on the core genus abundance of healthy and disease mice in dataset PRO (A) and CUG (E). Each circular or triangular point represented the fecal core microbiota of a healthy or disease mouse, respectively. P-values were calculated with MANOVA test, and shown on the top of PCA plots. (B,F) Spearman correlation between core microbiota represented by the three principal components (PCs) of PCA and host disease parameters of dataset PRO (B) and CUG (F). Columns corresponded to the three PCs, and rows to the host disease parameters. Red and blue color denoted positive and negative correlations, respectively. The color intensity represented the correlation degree. The black dots in the red/blue cells indicated the correlations were significant (false discovery rate < 0.1). (C,D,G) The different abundance of bacterial genus between two groups was identified by LEfSe. It was significantly different when alpha value of the factorial Kruskal–Wallis test was <0.05 and the logarithmic LDA score was >2.0. The left histogram showed LDA scores of differentially abundant genera between the two groups. The taxonomy was listed, followed by its core group. The right heatmap showed the relative abundance (log2 transformed) of differential genera. n = 8 for NC, HFD, HFD+probiotics 1, HFD+probiotics 2, and HFD+probiotics three groups (n = 24 for HFD+probiotics group in total); n = 12 for Control, NBW+APG, LBW+APG groups (n = 24 for APG group in total).
In dataset CUG, PCA plots exhibited the APG group (shown in the figure as APG) was significantly separated from the control group (shown in the figure as Control) (MANOVA, p = 0.033) (Figure 4E). Spearman correlation heatmap showed PC1 of PCA was significantly positively correlated with cecal butyrate and cecal isovalerate, whilst negatively correlated with MS indices as eAT/body weight and fasting blood glucose. Meanwhile, PC7 showed significantly positive correlation with eAT/body weight, perirenal (p) AT/body weight, subcutaneous (s) AT/body weight and fasting insulin, whilst negative correlation with cecal propionate, isobutyrate, valerate, and isovalerate (Figure 4F). LEfSe results showed that APG raised the abundances of Group C3 (Parabacteroides and Flavonifractor) and Group C5 Ruminococcus, whilst suppressed the abundance of Anaerostipes (C1), Erysipelotrichaceae_incertae_sedis (C4), and Clostridium XIVb (C1) (Figure 4G). Together, the 37 core genera were also critical to the APG-induced MS.
Discussion
Despite of the wide inter-individual structural variations, the composition of core gut bacteria and their roles in host phenotypes have not been completely solved yet. Hereinabove, we had exhibited a more comprehensive profile of core gut bacteria in healthy mouse model, excluding the influence of the diets, lifestyles and drugs. The core gut bacteria were divided into five core groups according to the co-abundant relationship, which might reflect their biological function to some extent. More importantly, the core gut bacteria were closely related to age, sex and health/disease state of the host, suggesting that the core bacteria might play a pivotal role in host gut homeostasis and health.
As mentioned above, the 37 core genera were identified across 101 healthy mice with different ages, sexes, and mouse stains at a prevalence of at least 50% in this study. Different from prior researches (Turnbaugh et al., 2009; Xiao et al., 2015; Zhang J. et al., 2015), we relaxed our restrictions to 50% prevalence to avoid losing the key taxa, which existed in different hosts but performed the same function. The 37 core bacteria presented 46.77% of total sequences on average, and could distinguish gut microbiota with age, sex, and health/disease state. These suggested that the 37 core genera could be the representatives of the whole gut microbiota and the indicators to the host phenotype. Notably, among the 37 core genera, Anaerostipes and Parabacteroides were the two most widespread bacteria, and Anaerostipes resided in the gut of all mice. It was reported that Anaerostipes species, such as A. butyraticus, and A. caccae, were butyrate-producing bacteria (Schwiertz et al., 2002; Eeckhaut et al., 2010), which was an important functional group in the gut of healthy hosts (Sekelja et al., 2011; Zhang J. et al., 2015; Cheng et al., 2016). Parabacteroides was at 99.0% prevalence. Previous studies showed P. distasonis and P. goldsteinii could induce gut or systemic inflammation (Awadel-Kariem et al., 2010; Kverka et al., 2011), and oral administration of P. distasonis antigens attenuated dextran sulfate sodium (DSS)-induced colitis in mice (Kverka et al., 2011). Considering the high prevalence (99%) and relatively high abundance (6.63% of relative abundance on average) of Parabacteroides in healthy mice, it is speculated that moderate Parabacteroides might stimulate the immune system and be essential to host health. Compared with previous studies, among the 37 core genera, 22 genera had been shown popular in the two previous studies, including Alistipes, Bacteroides, Clostridium XIVa, Clostridium XlVb, Helicobacter, Lachnospiracea_incertae_sedis, Lactobacillus, Odoribacter, Parabacteroides, Rikenella, Roseburia, Turicibacter, Prevotella, Marvinbryantia, Ruminococcus, Blautia, Pseudoflavonifractor, Anaerotruncus, Oscillibacter, Eggerthella, Olsenella, and Bifidobacterium (Benson et al., 2010; Xiao et al., 2015). Yet, 17 genera were reported in several human gut microbial compositional studies, such as Anaerostipes, Parabacteroides, Anaerotruncus, Bacteroides, Alistipes, Prevotella, Lactobacillus, Roseburia, Eggerthella, Lachnospiracea_incertae_sedis, Bifidobacterium, Clostridium XIVa, Clostridium XIVb, Blautia, Odoribacter, Ruminococcus and Oscillibacter (Tap et al., 2009; Turnbaugh et al., 2009; Qin et al., 2010; Claesson et al., 2011; Sekelja et al., 2011; Xiao et al., 2015; Zhang J. et al., 2015; Cheng et al., 2016; Maier et al., 2018). Therefore, there were commonalities between human and mouse core gut bacteria.
Now that these genera were indeed stable in the host gut, it is important to identify their function. In healthy mouse model, the natural ecological relationships were more visible, because the effect factors of gut microbiota, such as diets, lifestyles, and drugs, were controllable. Interestingly, these 37 core genera showed ecological networks and were divided into five core groups according to their co-abundant relationships. The core groups might represent the potential functional groups. Group C1, the largest group, included two predominant types of bacteria. One type of bacteria was butyrate-producing bacteria, such as Anaerostipes, Anaerotruncus, Clostridium XIVb, Intestinimonas, Roseburia, Allobaculum, Lachnospiracea_incertae_sedis (the nearest neighbor was Eubacterium ventriosum), Pseudoflavonifractor, and Acetatifactor (Barcenilla et al., 2000; Duncan et al., 2002; Schwiertz et al., 2002; Lau et al., 2006; Eeckhaut et al., 2010; Scott et al., 2011; Pfeiffer et al., 2012; Klaring et al., 2013; Machiels et al., 2014; Wang J. et al., 2015; Petri et al., 2017; Ueki et al., 2017). Butyrate is an essential bacterial metabolite from carbohydrate in the gut, as it is the preferred energy source for colon epithelial cells and maintains gut barrier functions (Peng et al., 2007). Meanwhile, butyrate has immunomodulatory and anti-inflammatory properties, and it attenuates inflammation by promoting differentiation and proliferation of regulatory T cells and inhibiting nuclear factor-κB activation (Furusawa et al., 2013; Kim et al., 2013; Smith et al., 2013). Administration of butyrate could alleviate gut inflammation (Singh et al., 2014). The other type of bacteria was immunity-activing or anti-infection or anti-tumor immunity-promoting bacteria, including Bacteroides, Barnesiella, Alistipes, Helicobacter, Prevotella, and Mucispirillum (Eaton et al., 2011; Ramanan et al., 2014; Bunker et al., 2015; Daillere et al., 2016; Kim and Kim, 2016; Moschen et al., 2016). These bacteria were critical to the immune systems of the host. If there were not these two types of bacteria in the gut, for example in germ-free mice, the immunity would be severely disordered (Hapfelmeier et al., 2010; Alam et al., 2011). Overall, these evidences indicated that Group C1 helped to shape and balance the immune system. Group C2 included carbohydrate-utilizing and lactate and/or acetate-producing bacteria, such as Bifidobacterium, Olsenella, and Lactobacillus (Fukuda et al., 2011; Matsuki et al., 2013; Li et al., 2015). Bifidobacterium and Lactobacillus were often used as probiotics to treat intestinal or systemic inflammation (Mohamadzadeh et al., 2011; Wang J. et al., 2015). Group C3 also contained some butyrate-producing bacteria, such as Clostridium XIVa and Blautia (Van den Abbeele et al., 2013; Takahashi et al., 2016), and some anti-inflammatory bacteria, like Parabacteroides and Flavonifractor (Kverka et al., 2011; Lee et al., 2018). It seemed that Group C3 had similar function with Group C1. This suggested that there were functional redundance in core gut bacteria. Functional redundance would keep system stable when it encountered perturbation. Group C5 only had one genus Ruminococcus, which contain bacteria conversing cellulose to acetate and hydrogen, such as R. flavefaciens and R. albus (Miller and Wolin, 1995; Koike and Kobayashi, 2001). Acetate and hydrogen were also showed to alleviate inflammation (Kajiya et al., 2009; Kim et al., 2013). Group C4 included two Actinobacteria Eggerthella and Gordonibacter, which were most phylogenetically closely to each other. Unlike the other core bacteria, these two bacteria metabolized amino acids rather than sugars (Lau et al., 2004; Woo et al., 2010). Altogether, members of the core microbiota in mice were principally carbohydrate- and amino acids-metabolizing bacteria and immunomodulatory bacteria. This was in accordance with previous functional studies on mice and human gut microbiota (Turnbaugh et al., 2009; Xiao et al., 2015).
In addition to intragroup co-abundance relationship, there were intergroup positive or negative correlation relationships. Group C1 and Group C2 were slightly positively correlated with each other. These interactions relied on nutrient and electron donor exchange (Zengler and Zaramela, 2018). For examples, polysaccharide, such as inulin, fructooligosaccharide or galactooligosaccharide, could be metabolized by Bifidobacterium or Lactobacillus in Group C2 into lactate and/or acetate. Then these products could be converted to butyrate by Anaerostipes, Eubacterium (the nearest neighbor of Lachnospiracea_incertae_sedis), Allobaculum or Roseburia in Group C1 (Schwiertz et al., 2002; Belenguer et al., 2006; Falony et al., 2006). That might be why Bifidobacterium or Lactobacillus could be chosen as probiotics to regulate the balance of gut microbiota. Group C3 was significantly correlated with Group C2 yet weakly negatively with Group C1. It might because Clostridium XIVa, Blautia, Parabacteroides, and Flavonifractor in Group C3 had the similar metabolic and anti-inflammatory functions with Group C1 and C2 (Kverka et al., 2011; Van den Abbeele et al., 2013; Takahashi et al., 2016; Lee et al., 2018). It also had been previously reported that Blautia (one genus in Group C3) were decreased by Lactobacillus (one genus in Group C2) (Wang R. et al., 2015). Bacteria in this group might be probiotics candidate and should be isolated, and their function should be tested in animal experiments. Group C4 was negatively correlated with Group C1, but the underlying mechanisms was unclear. Group C5 Ruminococcus was neutral with the other four groups, indicating that it might have more complete metabolic capacity and was independent to other members in the community. Together, these evidences revealed that there were social networks among the core gut bacteria in the healthy host. The network maintained the balance of host gut microenvironment. It is crucial to understand this complex network for developing strategies to reshape the gut microbial communities, such as making probiotic cocktail or standard fecal microbial transplantation.
We additionally demonstrated that age, sex and health/disease state were associated with the core gut bacteria. Unsupervised statistical method PCA showed that the 37 core genera were significantly separated according to different ages, sexes and health/disease state. Then, supervised multivariate statistical analysis was used to identify the specific genera in response to the corresponding variant. Compared with the infant mice, the adult mice had more abundant of probiotics-containing Group C2, which contrasted our mind that the infant gut resided more probiotics. Low abundant or lack of Group C2 in infancy might cause allergy, and supplement of Group C2 could prevent it (Kirjavainen et al., 2002; Kim et al., 2008). In the adult male C57BL/6J and ICR mice, Group C1 were increased in 24-week-old than in 12-week-old (Bacteroides, Mucispirillum, and Anaerotruncus for C57BL/6J mice; Prevotella, Alistipes, and Clostridium XIVb for ICR mice). Improving the immunoregulatory Group C1 and the probiotics-containing Group C2 with age in infancy and early adulthood would cause the core gut microbiota more and more mature and healthy (Zhang et al., 2012). Meanwhile, the female mice had more butyrate-producing and anti-inflammatory bacteria Anaerotruncus, Clostridium XIVb, Blautia, and Flavonifractor in Groups C1 and C3 than the male ones, suggesting that the females needed more anti-inflammatory bacteria to maintain intestinal health. This was in accordance with previous study that the female recipient mice transferred feces from a healthy man donor had more butyrate/acetate-producing bacteria and opportunistic pathogens (Wang J.J. et al., 2016). Moreover, in HFD-induced MS mice, butyrate-producing and immunity-promoting bacteria Group C1 (Anaerostipes, Barnesiella) and Group C2 (Bifidobacterium, Olsenella, Turicibacter) were decreased; while in APG group, Group C1 (Anaerostipes, Clostridium XIVb) were also reduced. These results indicated that the suffered MS was related to the reduction of Groups C1 and C2. In addition, the correlations between the core microbiota and the MS parameters were significantly. Together, the core microbiota identified in the present study distinguished gut microbiota not only with different ages and sexes in the healthy mice, but also with the health/disease state of mice. Thus, the core gut bacteria should be fully considered during disease treatment.
Although we identified a broad spectrum of mouse core gut bacteria covering 37 genera in this study, some bacteria that were dominant bacteria in human gut and related to the development or treatment of diseases were not in this list. As known, Faecalibacterium showed anti-inflammatory effects (Sokol et al., 2008), Akkermansia controlled obesity and diabetes (Plovier et al., 2017), and both of them were associated with anticancer immune responses in preclinical tumor models and cancer patients (Gopalakrishnan et al., 2018; Routy et al., 2018). However, these two important bacteria were not identified as the core in our present study. That might be due to Faecalibacterium inhabited in only 24 mice among the 101 healthy mice. Notably, 18 mice of these 24 ones were the mice transplanted with the healthy man feces, while only 6 mice were those with the native flora. These results suggested that Faecalibacterium were the dominant bacteria in human gut, but not in mouse gut. Besides, Akkermansia was only at a prevalence of 44.6% (45/101), and it might be identified as core gut bacteria if samples were enough. In addition to the protective bacteria Faecalibacterium and Akkermansia, the opportunistic pathogen Escherichia/Shigella (31/101) and Desulfovibrio (49/101) (Kotlowski et al., 2007; Dzierzewicz et al., 2010) were neither the core gut bacteria in this study. Similarly, Desulfovibrio might also be one of the core gut bacteria if the sample size enlarged. Among the 31 mice with Escherichia/Shigella, 14 mice were the ones with about 4,000 sequencing reads in one fecal sample (totally 19 mice), and 17 mice were with about 30,000 sequencing reads each sample (totally 82 mice). Thus, more core gut bacteria were identified if the sequencing depth was increasing.
In the present study, we clustered the data into genus to reduce batch effects in different sequencing sets (Duvallet et al., 2017). Another reason for the genus-level study was that the mouse gut microbiota was consistent with the humans’ ones at the genus level (Hugenholtz and de Vos, 2018). Even though we lost the sensitivity to detect species or strains, we found that the core gut bacteria, such as Bifidobacterium, Lactobacillus, Clostridium, and Ruminococcus, were in accordance with the metagenomics studies targeting bacterial strains (Qin et al., 2010; Xiao et al., 2015; Zhao et al., 2018). These indicated that the core study in genus level was preliminary but accurate.
Conclusion
In conclusion, our results exhibited a relatively complete picture of core gut bacteria in healthy mice that they were mainly composed of carbohydrate- and amino acids-metabolizing bacteria and immunomodulatory bacteria, as well as their ecological relationships of symbiosis or functional redundancy. It was worth noting that the core bacteria resident in healthy gut was of fundamental importance for not only the healthy physiological characteristics, but also the chronic metabolic disease phenotypes. Our study focused on the gut microbiota of healthy host, and would be an important complement to the case–control study. Although this list of healthy core gut bacteria may need to supplement and correct by more researches, especially human studies, and extend from taxon core ones to functional core ones, it could be a treasury to screen the probiotics candidates or fecal microbial transplantation components for prevention and treatment of the diseases.
Author Contributions
XW and JW designed the research. JS and JW provided the data sets. JW, TL, and JD analyzed the data. LT and JW discussed and interpreted the results. JW, TL, and LT wrote and revised the manuscript.
Funding
This work was supported by grants from National Natural Science Foundation of China (Grant No. 81600409) and Shanghai General Hospital (Grant No. 06N1702003).
Conflict of Interest Statement
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.
Acknowledgments
We thank J. Wang and Z. S. Xue for technical help on data analyzing.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00887/full#supplementary-material
References
Alam, C., Bittoun, E., Bhagwat, D., Valkonen, S., Saari, A., Jaakkola, U., et al. (2011). Effects of a germ-free environment on gut immune regulation and diabetes progression in non-obese diabetic (n.d.) mice. Diabetologia 54, 1398–1406. doi: 10.1007/s00125-011-2097-5
Awadel-Kariem, F. M., Patel, P., Kapoor, J., Brazier, J. S., and Goldstein, E. J. (2010). First report of Parabacteroides goldsteinii bacteraemia in a patient with complicated intra-abdominal infection. Anaerobe 16, 223–225. doi: 10.1016/j.anaerobe.2010.01.001
Barcenilla, A., Pryde, S. E., Martin, J. C., Duncan, S. H., Stewart, C. S., Henderson, C., et al. (2000). Phylogenetic relationships of butyrate-producing bacteria from the human gut. Appl. Environ. Microbiol. 66, 1654–1661. doi: 10.1128/aem.66.4.1654-1661.2000
Belenguer, A., Duncan, S. H., Calder, A. G., Holtrop, G., Louis, P., Lobley, G. E., et al. (2006). Two routes of metabolic cross-feeding between Bifidobacterium adolescentis and butyrate-producing anaerobes from the human gut. Appl. Environ. Microbiol. 72, 3593–3599. doi: 10.1128/AEM.72.5.3593-3599.2006
Benjamini, Y., Drai, D., Elmer, G., Kafkafi, N., and Golani, I. (2001). Controlling the false discovery rate in behavior genetics research. Behav. Brain Res. 125, 279–284. doi: 10.1016/s0166-4328(01)00297-2
Benson, A. K., Kelly, S. A., Legge, R., Ma, F. R., Low, S. J., Kim, J., et al. (2010). Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. Proc. Natl. Acad. Sci. U.S.A. 107, 18933–18938. doi: 10.1073/pnas.1007028107
Bunker, J. J., Flynn, T. M., Koval, J. C., Shaw, D. G., Meisel, M., McDonald, B. D., et al. (2015). Innate and adaptive humoral responses coat distinct commensal bacteria with immunoglobulin A. Immunity 43, 541–553. doi: 10.1016/j.immuni.2015.08.007
Cheng, J., Ringel-Kulka, T., Heikamp-de Jong, I., Ringel, Y., Carroll, I., de Vos, W. M., et al. (2016). Discordant temporal development of bacterial phyla and the emergence of core in the fecal microbiota of young children. ISME J. 10, 1002–1014. doi: 10.1038/ismej.2015.177
Claesson, M. J., Cusack, S., O’Sullivan, O., Greene-Diniz, R., de Weerd, H., Flannery, E., et al. (2011). Composition, variability, and temporal stability of the intestinal microbiota of the elderly. Proc. Natl. Acad. Sci. U.S.A. 108(Suppl. 1), 4586–4591. doi: 10.1073/pnas.1000097107
Claesson, M. J., Jeffery, I. B., Conde, S., Power, S. E., O’Connor, E. M., Cusack, S., et al. (2012). Gut microbiota composition correlates with diet and health in the elderly. Nature 488, 178–184. doi: 10.1038/nature11319
Conlon, M. A., and Bird, A. R. (2014). The impact of diet and lifestyle on gut microbiota and human health. Nutrients 7, 17–44. doi: 10.3390/nu7010017
Daillere, R., Vetizou, M., Waldschmitt, N., Yamazaki, T., Isnard, C., Poirier-Colame, V., et al. (2016). Enterococcus hirae and Barnesiella intestinihominis facilitate cyclophosphamide-induced therapeutic immunomodulatory effects. Immunity 45, 931–943. doi: 10.1016/j.immuni.2016.09.009
Duncan, S. H., Hold, G. L., Barcenilla, A., Stewart, C. S., and Flint, H. J. (2002). Roseburia intestinalis sp. nov., a novel saccharolytic, butyrate-producing bacterium from human faeces. Int. J. Syst. Evol. Microbiol. 52(Pt 5), 1615–1620. doi: 10.1099/00207713-52-5-1615
Duvallet, C., Gibbons, S. M., Gurry, T., Irizarry, R. A., and Alm, E. J. (2017). Meta-analysis of gut microbiome studies identifies disease-specific and shared responses. Nat. Commun. 8:1784. doi: 10.1038/s41467-017-01973-8
Dzierzewicz, Z., Szczerba, J., Lodowska, J., Wolny, D., Gruchlik, A., Orchel, A., et al. (2010). The role of Desulfovibrio desulfuricans lipopolysaccharides in modulation of periodontal inflammation through stimulation of human gingival fibroblasts. Arch. Oral Biol. 55, 515–522. doi: 10.1016/j.archoralbio.2010.05.001
Eaton, K. A., Opp, J. S., Gray, B. M., Bergin, I. L., and Young, V. B. (2011). Ulcerative typhlocolitis associated with Helicobacter mastomyrinus in telomerase-deficient mice. Vet. Pathol. 48, 713–725. doi: 10.1177/0300985810383876
Eeckhaut, V., Van Immerseel, F., Pasmans, F., De Brandt, E., Haesebrouck, F., Ducatelle, R., et al. (2010). Anaerostipes butyraticus sp. nov., an anaerobic, butyrate-producing bacterium from Clostridium cluster XIVa isolated from broiler chicken caecal content, and emended description of the genus Anaerostipes. Int. J. Syst. Evol. Microbiol. 60(Pt 5), 1108–1112. doi: 10.1099/ijs.0.015289-0
Falony, G., Vlachou, A., Verbrugghe, K., and De Vuyst, L. (2006). Cross-feeding between Bifidobacterium longum BB536 and acetate-converting, butyrate-producing colon bacteria during growth on oligofructose. Appl. Environ. Microbiol. 72, 7835–7841. doi: 10.1128/aem.01296-06
Franks, I. (2011). Gut microbiota: FUT2 genotype influences the gut microbiota in patients with Crohn’s disease and healthy individuals. Nat. Rev. Gastroenterol. Hepatol. 9:2. doi: 10.1038/nrgastro.2011.237
Fukuda, S., Toh, H., Hase, K., Oshima, K., Nakanishi, Y., Yoshimura, K., et al. (2011). Bifidobacteria can protect from enteropathogenic infection through production of acetate. Nature 469, 543–547. doi: 10.1038/nature09646
Furusawa, Y., Obata, Y., Fukuda, S., Endo, T. A., Nakato, G., Takahashi, D., et al. (2013). Commensal microbe-derived butyrate induces the differentiation of colonic regulatory T cells. Nature 504, 446–450. doi: 10.1038/nature12721
Gopalakrishnan, V., Spencer, C. N., Nezi, L., Reuben, A., Andrews, M. C., Karpinets, T. V., et al. (2018). Gut microbiome modulates response to anti-PD-1 immunotherapy in melanoma patients. Science 359, 97–103. doi: 10.1126/science.aan4236
Hapfelmeier, S., Lawson, M. A., Slack, E., Kirundi, J. K., Stoel, M., Heikenwalder, M., et al. (2010). Reversible microbial colonization of germ-free mice reveals the dynamics of IgA immune responses. Science 328, 1705–1709. doi: 10.1126/science.1188454
Hugenholtz, F., and de Vos, W. M. (2018). Mouse models for human intestinal microbiota research: a critical evaluation. Cell. Mol. Life Sci. 75, 149–160. doi: 10.1007/s00018-017-2693-8
Kajiya, M., Silva, M. J., Sato, K., Ouhara, K., and Kawai, T. (2009). Hydrogen mediates suppression of colon inflammation induced by dextran sodium sulfate. Biochem. Biophys. Res. Commun. 386, 11–15. doi: 10.1016/j.bbrc.2009.05.117
Kau, A. L., Ahern, P. P., Griffin, N. W., Goodman, A. L., and Gordon, J. I. (2011). Human nutrition, the gut microbiome and the immune system. Nature 474, 327–336. doi: 10.1038/nature10213
Kim, D., and Kim, W. U. (2016). Editorial: can prevotella copri be a causative pathobiont in rheumatoid arthritis? Arthritis Rheumatol. 68, 2565–2567. doi: 10.1002/art.39807
Kim, J. Y., Choi, Y. O., and Ji, G. E. (2008). Effect of oral probiotics (Bifidobacterium lactis AD011 and Lactobacillus acidophilus AD031) administration on ovalbumin-induced food allergy mouse model. J. Microbiol. Biotechnol. 18, 1393–1400.
Kim, M. H., Kang, S. G., Park, J. H., Yanagisawa, M., and Kim, C. H. (2013). Short-chain fatty acids activate GPR41 and GPR43 on intestinal epithelial cells to promote inflammatory responses in mice. Gastroenterology 145, 396–406. e1–10. doi: 10.1053/j.gastro.2013.04.056
Kirjavainen, P. V., Arvola, T., Salminen, S. J., and Isolauri, E. (2002). Aberrant composition of gut microbiota of allergic infants: a target of bifidobacterial therapy at weaning? Gut 51, 51–55. doi: 10.1136/gut.51.1.51
Klaring, K., Hanske, L., Bui, N., Charrier, C., Blaut, M., Haller, D., et al. (2013). Intestinimonas butyriciproducens gen. nov., sp. nov., a butyrate-producing bacterium from the mouse intestine. Int. J. Syst. Evol. Microbiol. 63(Pt 12), 4606–4612. doi: 10.1099/ijs.0.051441-0
Koike, S., and Kobayashi, Y. (2001). Development and use of competitive PCR assays for the rumen cellulolytic bacteria: fibrobacter succinogenes, Ruminococcus albus and Ruminococcus flavefaciens. FEMS Microbiol. Lett. 204, 361–366. doi: 10.1016/s0378-1097(01)00428-1
Kotlowski, R., Bernstein, C. N., Sepehri, S., and Krause, D. O. (2007). High prevalence of Escherichia coli belonging to the B2+D phylogenetic group in inflammatory bowel disease. Gut 56, 669–675. doi: 10.1136/gut.2006.099796
Kverka, M., Zakostelska, Z., Klimesova, K., Sokol, D., Hudcovic, T., Hrncir, T., et al. (2011). Oral administration of Parabacteroides distasonis antigens attenuates experimental murine colitis through modulation of immunity and microbiota composition. Clin. Exp. Immunol. 163, 250–259. doi: 10.1111/j.1365-2249.2010.04286.x
Lau, S. K., Woo, P. C., Woo, G. K., Fung, A. M., Ngan, A. H., Song, Y., et al. (2006). Bacteraemia caused by Anaerotruncus colihominis and emended description of the species. J. Clin. Pathol. 59, 748–752. doi: 10.1136/jcp.2005.031773
Lau, S. K., Woo, P. C., Woo, G. K., Fung, A. M., Wong, M. K., Chan, K. M., et al. (2004). Eggerthella hongkongensis sp. nov. and Eggerthella sinensis sp. nov., two novel Eggerthella species, account for half of the cases of Eggerthella bacteremia. Diagn. Microbiol. Infect. Dis. 49, 255–263. doi: 10.1016/j.diagmicrobio.2004.04.012
Lee, T. E., Shin, K. C., and Oh, D. K. (2018). Biotransformation of fructose to allose by a one-pot reaction using flavonifractor plautiiD-allulose 3-epimerase and Clostridium thermocellum ribose 5-phosphate isomerase. J. Microbiol. Biotechnol. 28, 418–424. doi: 10.4014/jmb.1709.09044
Li, X., Jensen, R. L., Hojberg, O., Canibe, N., and Jensen, B. B. (2015). Olsenella scatoligenes sp. nov., a 3-methylindole- (skatole) and 4-methylphenol- (p-cresol) producing bacterium isolated from pig faeces. Int. J. Syst. Evol. Microbiol. 65(Pt 4), 1227–1233. doi: 10.1099/ijs.0.000083
Machiels, K., Joossens, M., Sabino, J., De Preter, V., Arijs, I., Eeckhaut, V., et al. (2014). A decrease of the butyrate-producing species Roseburia hominis and Faecalibacterium prausnitzii defines dysbiosis in patients with ulcerative colitis. Gut 63, 1275–1283. doi: 10.1136/gutjnl-2013-304833
Maier, L., Pruteanu, M., Kuhn, M., Zeller, G., Telzerow, A., Anderson, E. E., et al. (2018). Extensive impact of non-antibiotic drugs on human gut bacteria. Nature 555, 623–628. doi: 10.1038/nature25979
Matsuki, T., Pedron, T., Regnault, B., Mulet, C., Hara, T., and Sansonetti, P. J. (2013). Epithelial cell proliferation arrest induced by lactate and acetate from Lactobacillus casei and Bifidobacterium breve. PLoS One 8:e63053. doi: 10.1371/journal.pone.0063053
Miller, T. L., and Wolin, M. J. (1995). Bioconversion of cellulose to acetate with pure cultures of Ruminococcus albus and a hydrogen-using acetogen. Appl. Environ. Microbiol. 61, 3832–3835.
Mohamadzadeh, M., Pfeiler, E. A., Brown, J. B., Zadeh, M., Gramarossa, M., Managlia, E., et al. (2011). Regulation of induced colonic inflammation by Lactobacillus acidophilus deficient in lipoteichoic acid. Proc. Natl. Acad. Sci. U.S.A. 108(Suppl. 1), 4623–4630. doi: 10.1073/pnas.1005066107
Moschen, A. R., Gerner, R. R., Wang, J., Klepsch, V., Adolph, T. E., Reider, S. J., et al. (2016). Lipocalin 2 protects from inflammation and tumorigenesis associated with gut microbiota alterations. Cell Host Microbe 19, 455–469. doi: 10.1016/j.chom.2016.03.007
Peng, L., He, Z., Chen, W., Holzman, I. R., and Lin, J. (2007). Effects of butyrate on intestinal barrier function in a Caco-2 cell monolayer model of intestinal barrier. Pediatr. Res. 61, 37–41. doi: 10.1203/01.pdr.0000250014.92242.f3
Petri, R. M., Pourazad, P., Khiaosa-Ard, R., Klevenhusen, F., Metzler-Zebeli, B. U., and Zebeli, Q. (2017). Temporal dynamics of in-situ fiber-adherent bacterial community under ruminal acidotic conditions determined by 16S rRNA gene profiling. PLoS One 12:e0182271. doi: 10.1371/journal.pone.0182271
Pfeiffer, N., Desmarchelier, C., Blaut, M., Daniel, H., Haller, D., and Clavel, T. (2012). Acetatifactor muris gen. nov., sp. nov., a novel bacterium isolated from the intestine of an obese mouse. Arch. Microbiol. 194, 901–907. doi: 10.1007/s00203-012-0822-1
Plovier, H., Everard, A., Druart, C., Depommier, C., Van Hul, M., Geurts, L., et al. (2017). A purified membrane protein from Akkermansia muciniphila or the pasteurized bacterium improves metabolism in obese and diabetic mice. Nat. Med. 23, 107–113. doi: 10.1038/nm.4236
Qin, J., Li, R., Raes, J., Arumugam, M., Burgdorf, K. S., Manichanh, C., et al. (2010). A human gut microbial gene catalogue established by metagenomic sequencing. Nature 464, 59–65. doi: 10.1038/nature08821
Ramanan, D., Tang, M. S., Bowcutt, R., Loke, P., and Cadwell, K. (2014). Bacterial sensor Nod2 prevents inflammation of the small intestine by restricting the expansion of the commensal Bacteroides vulgatus. Immunity 41, 311–324. doi: 10.1016/j.immuni.2014.06.015
Round, J. L., and Mazmanian, S. K. (2009). The gut microbiota shapes intestinal immune responses during health and disease. Nat. Rev. Immunol. 9, 313–323. doi: 10.1038/nri2515
Routy, B., Le Chatelier, E., Derosa, L., Duong, C. P. M., Alou, M. T., Daillere, R., et al. (2018). Gut microbiome influences efficacy of PD-1-based immunotherapy against epithelial tumors. Science 359, 91–97. doi: 10.1126/science.aan3706
Schirmer, M., Franzosa, E. A., Lloyd-Price, J., McIver, L. J., Schwager, R., Poon, T. W., et al. (2018). Dynamics of metatranscription in the inflammatory bowel disease gut microbiome. Nat. Microbiol. 3, 337–346. doi: 10.1038/s41564-017-0089-z
Schwiertz, A., Hold, G. L., Duncan, S. H., Gruhl, B., Collins, M. D., Lawson, P. A., et al. (2002). Anaerostipes caccae gen. nov., sp. nov., a new saccharolytic, acetate-utilising, butyrate-producing bacterium from human faeces. Syst. Appl. Microbiol. 25, 46–51. doi: 10.1078/0723-2020-00096
Scott, K. P., Martin, J. C., Chassard, C., Clerget, M., Potrykus, J., Campbell, G., et al. (2011). Substrate-driven gene expression in Roseburia inulinivorans: importance of inducible enzymes in the utilization of inulin and starch. Proc. Natl. Acad. Sci. U.S.A. 108(Suppl. 1), 4672–4679. doi: 10.1073/pnas.1000091107
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Sekelja, M., Berget, I., Naes, T., and Rudi, K. (2011). Unveiling an abundant core microbiota in the human adult colon by a phylogroup-independent searching approach. ISME J. 5, 519–531. doi: 10.1038/ismej.2010.129
Sekirov, I., Russell, S. L., Antunes, L. C., and Finlay, B. B. (2010). Gut microbiota in health and disease. Physiol. Rev. 90, 859–904. doi: 10.1152/physrev.00045.2009
Singh, N., Gurav, A., Sivaprakasam, S., Brady, E., Padia, R., Shi, H., et al. (2014). Activation of Gpr109a, receptor for niacin and the commensal metabolite butyrate, suppresses colonic inflammation and carcinogenesis. Immunity 40, 128–139. doi: 10.1016/j.immuni.2013.12.007
Sivan, A., Corrales, L., Hubert, N., Williams, J. B., Aquino-Michaels, K., Earley, Z. M., et al. (2015). Commensal Bifidobacterium promotes antitumor immunity and facilitates anti-PD-L1 efficacy. Science 350, 1084–1089. doi: 10.1126/science.aac4255
Smith, P. M., Howitt, M. R., Panikov, N., Michaud, M., Gallini, C. A., Bohlooly, Y. M., et al. (2013). The microbial metabolites, short-chain fatty acids, regulate colonic Treg cell homeostasis. Science 341, 569–573. doi: 10.1126/science.1241165
Sokol, H., Pigneur, B., Watterlot, L., Lakhdari, O., Bermudez-Humaran, L. G., Gratadoux, J. J., et al. (2008). Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc. Natl. Acad. Sci. U.S.A. 105, 16731–16736. doi: 10.1073/pnas.0804812105
Takahashi, K., Nishida, A., Fujimoto, T., Fujii, M., Shioya, M., Imaeda, H., et al. (2016). Reduced abundance of butyrate-producing bacteria species in the fecal microbial community in Crohn’s disease. Digestion 93, 59–65. doi: 10.1159/000441768
Tap, J., Mondot, S., Levenez, F., Pelletier, E., Caron, C., Furet, J. P., et al. (2009). Towards the human intestinal microbiota phylogenetic core. Environ. Microbiol. 11, 2574–2584. doi: 10.1111/j.1462-2920.2009.01982.x
Tremaroli, V., and Backhed, F. (2012). Functional interactions between the gut microbiota and host metabolism. Nature 489, 242–249. doi: 10.1038/nature11552
Turnbaugh, P. J., Hamady, M., Yatsunenko, T., Cantarel, B. L., Duncan, A., Ley, R. E., et al. (2009). A core gut microbiome in obese and lean twins. Nature 457, 480–484. doi: 10.1038/nature07540
Ueki, A., Goto, K., Ohtaki, Y., Kaku, N., and Ueki, K. (2017). Description of Anaerotignum aminivorans gen. nov., sp. nov., a strictly anaerobic, amino-acid-decomposing bacterium isolated from a methanogenic reactor, and reclassification of Clostridium propionicum, Clostridium neopropionicum and Clostridium lactatifermentans as species of the genus Anaerotignum. Int. J. Syst. Evol. Microbiol. 67, 4146–4153. doi: 10.1099/ijsem.0.002268
Van den Abbeele, P., Belzer, C., Goossens, M., Kleerebezem, M., De Vos, W. M., Thas, O., et al. (2013). Butyrate-producing Clostridium cluster XIVa species specifically colonize mucins in an in vitro gut model. ISME J. 7, 949–961. doi: 10.1038/ismej.2012.158
Wang, J., Tang, H., Wang, X., Zhang, X., Zhang, C., Zhang, M., et al. (2016). The structural alteration of gut microbiota in low-birth-weight mice undergoing accelerated postnatal growth. Sci. Rep. 6:27780. doi: 10.1038/srep27780
Wang, J. J., Wang, J., Pang, X. Y., Zhao, L. P., Tian, L., and Wang, X. P. (2016). Sex differences in colonization of gut microbiota from a man with short-term vegetarian and inulin-supplemented diet in germ-free mice. Sci. Rep. 6:36137. doi: 10.1038/srep36137
Wang, J., Tang, H., Zhang, C., Zhao, Y., Derrien, M., Rocher, E., et al. (2015). Modulation of gut microbiota during probiotic-mediated attenuation of metabolic syndrome in high fat diet-fed mice. ISME J. 9, 1–15. doi: 10.1038/ismej.2014.99
Wang, R., Chen, S., Jin, J., Ren, F., Li, Y., Qiao, Z., et al. (2015). Survival of Lactobacillus casei strain Shirota in the intestines of healthy Chinese adults. Microbiol. Immunol. 59, 268–276. doi: 10.1111/1348-0421.12249
Wong, S. H., Zhao, L., Zhang, X., Nakatsu, G., Han, J., Xu, W., et al. (2017). Gavage of fecal samples from patients with colorectal cancer promotes intestinal carcinogenesis in germ-free and conventional mice. Gastroenterology 153, 1621.e6–1633.e6. doi: 10.1053/j.gastro.2017.08.022
Woo, P. C., Teng, J. L., Lam, K. K., Tse, C. W., Leung, K. W., Leung, A. W., et al. (2010). First report of Gordonibacter pamelaeae bacteremia. J. Clin. Microbiol. 48, 319–322.
Xiao, L., Feng, Q., Liang, S., Sonne, S. B., Xia, Z., Qiu, X., et al. (2015). A catalog of the mouse gut metagenome. Nat. Biotechnol. 33, 1103–1108. doi: 10.1038/nbt.3353
Zengler, K., and Zaramela, L. S. (2018). The social network of microorganisms – How auxotrophies shape complex communities. Nat. Rev. Microbiol. 16, 383–390. doi: 10.1038/s41579-018-0004-5
Zhang, C., Yin, A., Li, H., Wang, R., Wu, G., Shen, J., et al. (2015). Dietary modulation of gut microbiota contributes to alleviation of both genetic and simple obesity in children. EBioMedicine 2, 968–984. doi: 10.1016/j.ebiom.2015.07.007
Zhang, J., Guo, Z., Xue, Z., Sun, Z., Zhang, M., Wang, L., et al. (2015). A phylo-functional core of gut microbiota in healthy young Chinese cohorts across lifestyles, geography and ethnicities. ISME J. 9, 1979–1990. doi: 10.1038/ismej.2015.11
Zhang, C., Zhang, M., Pang, X., Zhao, Y., Wang, L., and Zhao, L. (2012). Structural resilience of the gut microbiota in adult mice under high-fat dietary perturbations. ISME J. 6, 1848–1857. doi: 10.1038/ismej.2012.27
Keywords: core gut bacteria, ecological relationship, carbohydrate- and amino acid-metabolizing bacteria, immunomodulatory bacteria, probiotics, fecal microbiota transplantation
Citation: Wang J, Lang T, Shen J, Dai J, Tian L and Wang X (2019) Core Gut Bacteria Analysis of Healthy Mice. Front. Microbiol. 10:887. doi: 10.3389/fmicb.2019.00887
Received: 05 September 2018; Accepted: 05 April 2019;
Published: 24 April 2019.
Edited by:
George Tsiamis, University of Patras, GreeceReviewed by:
Jinshui Zheng, Huazhong Agricultural University, ChinaAggeliki Saridaki, Technical University of Crete, Greece
Copyright © 2019 Wang, Lang, Shen, Dai, Tian and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ling Tian, tl09168@hotmail.com Xingpeng Wang, richardwangxp@163.com
†These authors have contributed equally to this work