Composition of the gut microbiota in the three groups
The rarefaction curve of each sample was close to saturation, the amount of sequencing data was reasonable, and the sequencing depth was sufficient (Figure 1A). Sample valid sequences with 97% identity were clustered into OTUs. The gut microbiota was annotated with representative OTUs, with a total of 2,126 OTUs. Moreover, 660, 616, and 850 OTUs were identified in the SC, CLP, and GM groups, respectively, of which 495 OTUs were shared by three groups (Figure 1B). In the phylum-level gut microbiota, the predominant phyla in the SC and GM groups were Firmicutes, Proteobacteria, and Verrucomicrobiota in the CLP group (Figures 1C). The Firmicutes–Proteobacteria ratio was significantly decreased in the CLP group (Figures 1D). The predominant genera in the SC and GM groups were Lactobacillus, Ruminococcus, and Blautia. The predominant genera in the CLP group were Escherichia–Shigella, Akkermansia, and Enterococcus (Figure 1E).
α-Diversity analysis
The Simpson index represented the diversity and evenness of species distribution. The Wilcoxon rank-sum test showed that the diversity and evenness (Simpson index) of species distribution in the SC and GM groups were higher than those in the CLP group. The Shannon index, which measures the total number of microbial communities and their proportions, was also higher in the SC and GM groups than in the CLP group. The Chao1 index was used to evaluate the total number of species. The Chao1 index was higher in the SC group than in the CLP and GM groups. ACE was used to estimate the number of OTUs in the community. The ACE was higher in the SC group than in the CLP and GM groups. Good’s coverage index represents the sequencing depth. Good’s coverage rates of the three groups were all greater than 99.8%, indicating that the sequencing depth of the gut microbiota was ideal. The PD whole tree index represents the genetic relationship of the species in the community. The PD whole tree index in the SC and GM groups was higher than in CLP group (Table 1).
β-Diversity analysis
To assess the degree of similarity between microbial communities, conducted sampling comparative analysis (β-diversity analysis) for different groups. The PCoA results showed that samples from the SC, CLP, and GM groups formed three distinct clusters, among which samples from the SC and GM groups tended to cluster together, indicating that the two groups had similar species composition. The CLP group was far away, indicating that it had a larger community difference compared with the SC and GM groups (Figure 1F).
Differential microbiota analysis in the SC, CLP, and GM groups
To identify the specific microbiota associated with the SC, CLP, and GM groups, we used LEfSe analysis. The LSD score was used to evaluate the influence of the relative abundance of species in different groups on the difference effect (LDA > 4 and P < 0.05). There are three differential genera, namely, Lactobacillus, Romboutsia, and Escherichia–Shigella, between the CLP and SC groups (Figure 2A). Similarly, there were five differential genera between the CLP and GM groups, namely, Lactobacillus, Pygmaiobacter, Erysipelotrichaceae_UCG_003, Escherichia–Shigella, and Akkermansia (Figure 2B). The only differential genera between the SC group and GM group was Parasutterella (Figure 2C).
Correlation analysis of genera between the CLP and GM groups
The gut microbiota is interconnected to maintain homeostasis, and studying the interactions between different genera contributes to finding the important role of microbiota of the CLP and GM groups (Spearman correlation coefficients <−0.80 or >0.80, P < 0.05). Notably, the GM-enriched genera had higher associations than the CLP group-enriched genera. Akkermansia was strongly positively correlated with Desulfovibrio. Lachnospiraceae UCG-006 was strongly negatively correlated with Lactobacillus in the CLP group. Moreover, Lactobacillus and Desulfovibrio were strongly negatively correlated (Supplementary Figure 1). In the GM group, Ruminococcus was strongly negatively correlated with Escherichia–Shigella and Enterobacteria, and Candidatus Saccharimonas was strongly negatively correlated with Escherichia–Shigella (Supplementary Figure 2). The enriched Muribaculum in the CLP group was strongly negatively correlated with the enriched Lactobacillus in the GM group and strongly positively associated with the enriched Lachnospiraceae UCG-006 in the CLP group. The enriched Akkermansia in the CLP group was strongly positively correlated with Terrisporobacter, Proteus, and Odoribacter. Candidatus Saccharimonas enriched in the GM group was strongly negatively correlated with Escherichia–Shigella enriched in the CLP group (Supplementary Figure 3).
Functional alterations of GM on gut microbiota in the CLP group
To characterize the functional alterations of gut microbiota following GM application in rats with sepsis, we used Tax4Fun to analyze 16S rRNA and predicted functional composition profiles. Eighteen pathways were differentially enriched (P < 0.05) between the CLP and GM groups (Supplementary Figure 4). Then, we investigated the correlation between different species (n = 27) and KEGG pathways (n = 44) using Spearman correlation analysis (Student t-test, Spearman correlation coefficient <−0.80 or >0.80, P < 0.05) (Supplementary Figure 5). CLP-enriched Akkermansia was negatively correlated with phospholipase D signaling pathway and carbohydrate metabolism signaling pathway. GM-enriched Lactobacillus is associated with several pathways, such as cell cycle, tight junction, phosphotransferase system, isoflavone biosynthesis, caffeine metabolism, thyroid hormone synthesis, mineral absorption, and cellular motility and secretion, and positively correlated with transport and cancer pathways. GM-enriched Ruminococcus was positively correlated with amino sugar and nucleotide sugar metabolism, while Lactobacillus was negatively associated with apoptosis, lipopolysaccharide biosynthesis, steroid hormone biosynthesis, transport, sphingolipid metabolism, and pathways in cancer, fluid shear stress, and atherosclerosis. The dominant Phascolarctobacterium in the CLP group was positively correlated with G protein coupled receptors. Moreover, the dominant Colidextribacter in the CLP group was positively correlated with polyketide sugar unit biosynthesis, terpomycin biosynthesis, and N glycan biosynthesis (Supplementary Figure 6).
Metabolite analysis of the SC, CLP, and GM groups
Six samples and a QC sample were run during the running process. The aggregation of the QC samples was good, indicating that the method was stable and experimental data had high quality (Supplementary Figure 7). The total ion chromatograms of the SC, CLP, and GM groups were observed (Supplementary Figure 8). In PCA, the SC group was closer to the GM group and was significantly separated from the CLP group (Figures 3A and 3D). Using OPLS-DA, the SC group and CLP group and CLP group and GM group can be divided according to their metabolic differences, and the distinction is significant (Figures 3B-C, 3E-F).
Differential metabolite identification
Multivariate statistical analysis was used to screen differential metabolites (VIP > 1.0 and P < 0.05). In the positive ion mode, the SC group had 28 differential metabolites compared with the CLP group, eight were significantly upregulated in the SC group, and 20 were significantly downregulated in the SC group (Figure 4A). The GM group had 26 differential metabolites compared with the CLP group: 13 were significantly upregulated and 13 were significantly downregulated (Figure 4B). In the negative ion mode, the SC group had 17 differential metabolites compared with the CLP group: 12 were significantly upregulated and 5 were significantly downregulated (Figure 4C). The GM group had 20 differential metabolites compared with the CLP group: 16 were significantly upregulated and 4 were significantly downregulated (Figure 4D). There were 15 significantly changed metabolites in the SC, CLP, and GM groups. The detailed information of each metabolite is shown in Table 2. As shown in Figure 4E, 15 metabolites are mainly divided into 11 fatty acids and 4 other classified metabolites, namely, 10E, 12Z-octadecadienoic acid, 9,10-epoxyoctadecenoic acid, 15-KETE, and alpha-linolenic acid. These four fatty acids are grouped together, and the relative abundance is lower in the CLP group but higher in SC and GM groups.
Metabolic pathway analysis of key differential metabolites
Using MetaboAnalyst 3.0 and the KEGG database, the different metabolite pathways between the CLP and GM groups were analyzed. It was found that GM mainly regulates sphingolipid metabolism, histidine metabolism, steroid biosynthesis, glycerophospholipid metabolism, primary bile acid biosynthesis, and alpha-linolenic acid metabolism (Supplementary Figure 9). They are consistent with the predicted pathway of Tax4Fun function, and the GM group was negatively related to apoptosis, lipopolysaccharide biosynthesis, steroid hormone biosynthesis, and sphingolipid metabolism. Therefore, changes in gut metabolites after GM application may be related to gut microbiota dysbiosis.
Links between different genera and metabolites
To further investigate the microbiota–intestinal metabolite relationship after GM application, Spearman’s correlation analysis was used to determine the connections between 27 different intestinal genera and differential metabolites between the CLP and GM groups. Moreover, we used a heatmap to illustrate their correlations (Spearman’s correlation coefficient <−0.80 or >0.80, P < 0.05). Some CLP-rich genera (e.g., Akkermansia) are positively associated with chenodeoxycholic acid 3-sulfate, cholesterol, lysophosphatidylethanolamine (LysoPE) (16:1(9Z)/0:0), LysoPE (16:0/0:0), LysoPE (15:0/0:0), and PC-M6 and negatively correlated with 15-KETE, sphinganine, and Choline. CLP-rich genera (e.g., Escherichia–Shigella) was positively correlated with chenodeoxycholic acid 3-sulfate, calystegin A3, and [10]-gingerdione and negatively correlated with 15-KETE. CLP-rich genera (e.g., Proteus) was positively correlated with chenodeoxycholic acid 3-sulfate, cholesterol, LysoPE (16:0/0:0), LysoPE (15:0/0:0), and [10]-gingerdione and negatively correlated with 15-KETE. Moreover, the GM-rich genus (Lactobacillus) was negatively correlated with LysoPE (15:0/0:0). Among them, LysoPE (15:0/0:0) was positively correlated with CLP-rich genera (Akkermansia and Proteus) and had a negative correlation with GM-rich genera (Lactobacillus) (Figure 5A).
As shown in Figure 5B, some CLP-rich genera (e.g., Akkermansia) were positively correlated with quinolinic acid and negatively correlated with leukotriene B4, 9,10-DHOME, deoxycholic acid, L-cysteine, phenylacetylglycine, N-acetylvaline, and alpha-linolenic acid. Escherichia–Shigella enriched in the CLP group was positively correlated with quinolinic acid and 5a-tetrahydrocorticosterone and negatively correlated with leukotriene B4, 9,10-DHOME, and m-coumaric acid. CLP-rich genera (Proteus) were positively correlated with quinolinic acid, 5a-tetrahydrocorticosterone, and cortisol and negatively correlated with leukotriene B4, 9,10-DHOME, deoxycholic acid, N-alpha-acetyllysine, oleoyl glycine, m-coumaric acid, and alpha-linolenic acid. Furthermore, GM-rich genus (Lactobacillus) was positively correlated with leukotriene B4, 9,10-DHOME, N-alpha-acetyllysine, N-acetylvaline, phenylacetylglycine, and oleoyl glycine and negatively correlated with quinolinic acid. Among them, leukotriene B4 and 9,10-DHOME were negatively correlated with CLP-rich genera (Akkermansia, Escherichia–Shigella, and Proteus) and positively correlated with the GM-rich genus (Lactobacillus). N-Alpha-acetyllysine and oleoyl glycine were negatively correlated with the CLP-rich genus (Proteus) and positively correlated with the GM-rich genus (Lactobacillus). N-Acetylvaline and phenylacetylglycine were negatively correlated with CLP-rich Akkermansia and positively correlated with GM-rich Lactobacillus. Quinolinic acid was positively correlated with CLP-rich genera (Akkermansia, Escherichia–Shigella, and Proteus) and negatively correlated with GM-rich genera (Lactobacillus).
Effects of GM on lung damage and survival in rats with sepsis
The mortality rate in the CLP group (10/20) at 24 h after operation was higher than in the SC group (0/20) (P = 9.60335E-05) and GM group (4/20) (P = 0.048). Inflammatory infiltration, edema, and structural damage of the lung tissue in the CLP group were significantly worse than those in SC and GM groups (Figure 6).