1,575
Views
0
CrossRef citations to date
0
Altmetric
Research Paper

Multi-omics reveals deoxycholic acid modulates bile acid metabolism via the gut microbiota to antagonize carbon tetrachloride-induced chronic liver injury

, , , , , , , , , & show all
Article: 2323236 | Received 21 Nov 2023, Accepted 21 Feb 2024, Published online: 28 Feb 2024

Figures & data

Figure 1. Experimental design and procedure HC=Healthy control; CLD=Chronic liver disease; DCA=Deoxycholic acid; Abx=Antibiotic cocktail.

Figure 1. Experimental design and procedure HC=Healthy control; CLD=Chronic liver disease; DCA=Deoxycholic acid; Abx=Antibiotic cocktail.

Figure 2. Characteristics of bile acids profiles in CLD patients (a) heatmap cluster analysis comparing CLD patients and healthy controls. The left cluster line represents the metabolite cluster, and the upper cluster line represents the sample cluster; (b) OPLS-DA score plot illustrating the differentiation between CLD patients and healthy controls; (c) bar chart depicting the differences in metabolite levels between CLD patients and healthy individuals; (d) histograms displaying the levels of 23-DCA, 3β-DCA, and DCA in CLD patients compared to healthy controls. HC: healthy control (n = 20); NASH: nonalcoholic steatohepatitis (n = 10); PBC: primary biliary cholangitis (n = 10); ALD: alcoholic liver disease (n = 10); CHB: chronic hepatitis B; CLD: chronic liver disease (n = 40), *P < .05, **P < .01, ***P < .001.

Figure 2. Characteristics of bile acids profiles in CLD patients (a) heatmap cluster analysis comparing CLD patients and healthy controls. The left cluster line represents the metabolite cluster, and the upper cluster line represents the sample cluster; (b) OPLS-DA score plot illustrating the differentiation between CLD patients and healthy controls; (c) bar chart depicting the differences in metabolite levels between CLD patients and healthy individuals; (d) histograms displaying the levels of 23-DCA, 3β-DCA, and DCA in CLD patients compared to healthy controls. HC: healthy control (n = 20); NASH: nonalcoholic steatohepatitis (n = 10); PBC: primary biliary cholangitis (n = 10); ALD: alcoholic liver disease (n = 10); CHB: chronic hepatitis B; CLD: chronic liver disease (n = 40), *P < .05, **P < .01, ***P < .001.

Figure 3. Effect of DCA on the pathological and plasma biochemical indexes of liver inflammation in mice (a) typical appearance of liver; (b) Representative liver histology by H&E staining (×100 and × 200); (c) immunohistochemical staining of F4/80 in liver tissue (×100 and × 200); (d) immunohistochemical staining of CD86 in liver tissue (×100 and × 200); (e) liver index of mice in different groups; (f) quantitative results of histochemical staining of F4/80 and CD86 in liver tissue. (g) Plasma ALT and AST levels. n = 6,*P < .05,**P < .01,***P < .001,****P < .0001. “ns” indicates no significant difference.

Figure 3. Effect of DCA on the pathological and plasma biochemical indexes of liver inflammation in mice (a) typical appearance of liver; (b) Representative liver histology by H&E staining (×100 and × 200); (c) immunohistochemical staining of F4/80 in liver tissue (×100 and × 200); (d) immunohistochemical staining of CD86 in liver tissue (×100 and × 200); (e) liver index of mice in different groups; (f) quantitative results of histochemical staining of F4/80 and CD86 in liver tissue. (g) Plasma ALT and AST levels. n = 6,*P < .05,**P < .01,***P < .001,****P < .0001. “ns” indicates no significant difference.

Figure 4. Effect of DCA on the hepatic inflammatory cells and plasma inflammatory factors in mice (a) detection of inflammatory cells in liver by flow cytometry [flow cytometry plot (left) and bar graph (right)]. Liver macrophages were identified as CD45+F4/80+CD11b+cells. Liver monocytes were identified as CD45+CD11b+Ly6C hi cells. Liver neutrophils were identified as CD45+Ly6G+CD11b+cells; (b) plasma levels of TNF-α, IL-6 and IL-1β. n = 6, *P < .05, **P < .01, ***P < .001, ****P < .0001. “ns” indicates no significant difference.

Figure 4. Effect of DCA on the hepatic inflammatory cells and plasma inflammatory factors in mice (a) detection of inflammatory cells in liver by flow cytometry [flow cytometry plot (left) and bar graph (right)]. Liver macrophages were identified as CD45+F4/80+CD11b+cells. Liver monocytes were identified as CD45+CD11b+Ly6C hi cells. Liver neutrophils were identified as CD45+Ly6G+CD11b+cells; (b) plasma levels of TNF-α, IL-6 and IL-1β. n = 6, *P < .05, **P < .01, ***P < .001, ****P < .0001. “ns” indicates no significant difference.

Figure 5. Effect of DCA on the pathological and mRNA expression of liver fibrosis in mice (a) Masson staining of liver tissue (×100 and × 200); (b) immunohistochemical staining of α-SMA (×100 and × 200); (c) immunohistochemical staining of collagen I in liver tissue (×100 and × 200); (d) quantitative results of histochemical staining for α-SMA and collagen I; (e) mRNA expression of α-SMA and collagen I in liver tissue. n = 6, *P < .05, **P < .01, ***P < .001, ****P < .0001. “ns” indicates no significant difference.

Figure 5. Effect of DCA on the pathological and mRNA expression of liver fibrosis in mice (a) Masson staining of liver tissue (×100 and × 200); (b) immunohistochemical staining of α-SMA (×100 and × 200); (c) immunohistochemical staining of collagen I in liver tissue (×100 and × 200); (d) quantitative results of histochemical staining for α-SMA and collagen I; (e) mRNA expression of α-SMA and collagen I in liver tissue. n = 6, *P < .05, **P < .01, ***P < .001, ****P < .0001. “ns” indicates no significant difference.

Figure 6. DCA modulates the composition of gut microbiota (a) α-diversity assessed using Chao1 Index, Observed_ species Index, and Shannon Index for each group; (b) Principal coordinate analysis (PCoA) of gut microbiota; (C) the relative abundance of the top 10 bacteria at the phylum level; (D) Representative histograms of gut microbiota at the phylum level; (E) the relative abundance of the top 10 bacteria at the genus level; (f) Representative histograms of gut microbiota at the genus level; (g) the relative abundance of the top 10 bacteria at the species level; (H) Representative histograms of gut microbiota at the species level. Ctrl: Control; CD: CCl4+DCA; CDA: CCl4+DCA+Abx. n = 6, *P < .05, **P < .01, ***P < .001, ****P < .0001. “ns” indicates no significant difference.

Figure 6. DCA modulates the composition of gut microbiota (a) α-diversity assessed using Chao1 Index, Observed_ species Index, and Shannon Index for each group; (b) Principal coordinate analysis (PCoA) of gut microbiota; (C) the relative abundance of the top 10 bacteria at the phylum level; (D) Representative histograms of gut microbiota at the phylum level; (E) the relative abundance of the top 10 bacteria at the genus level; (f) Representative histograms of gut microbiota at the genus level; (g) the relative abundance of the top 10 bacteria at the species level; (H) Representative histograms of gut microbiota at the species level. Ctrl: Control; CD: CCl4+DCA; CDA: CCl4+DCA+Abx. n = 6, *P < .05, **P < .01, ***P < .001, ****P < .0001. “ns” indicates no significant difference.

Figure 7. LEfse analysis of DCA modulating the composition of gut microbiota (a) the classification branch diagram illustrates the hierarchical relationship of main taxa from phylum to genus (from inner circle to outer circle) in the sample community. The node size corresponds to the average relative abundance of the taxon, and nodes of each color represent significant inter-group differences in these taxa. The abundance is higher in the group samples represented by the color; (b) the distribution bar chart of linear discriminant analysis (LDA) values for significantly different species demonstrates the significantly enriched species within each group and their importance. The LDA score reflects significant differences in bacteria among each group, and only bacteria meeting the significant LDA threshold of 2 are shown. The vertical axis represents taxonomic units with significant differences between groups, while the horizontal axis displays the logarithmic scores of LDA analysis for each taxonomic unit. Ctrl: Control; CD: CCl4+DCA; CDA: CCl4+DCA+Abx. n = 6.

Figure 7. LEfse analysis of DCA modulating the composition of gut microbiota (a) the classification branch diagram illustrates the hierarchical relationship of main taxa from phylum to genus (from inner circle to outer circle) in the sample community. The node size corresponds to the average relative abundance of the taxon, and nodes of each color represent significant inter-group differences in these taxa. The abundance is higher in the group samples represented by the color; (b) the distribution bar chart of linear discriminant analysis (LDA) values for significantly different species demonstrates the significantly enriched species within each group and their importance. The LDA score reflects significant differences in bacteria among each group, and only bacteria meeting the significant LDA threshold of 2 are shown. The vertical axis represents taxonomic units with significant differences between groups, while the horizontal axis displays the logarithmic scores of LDA analysis for each taxonomic unit. Ctrl: Control; CD: CCl4+DCA; CDA: CCl4+DCA+Abx. n = 6.

Figure 8. DCA regulates plasma bile acid profiles through gut microbiota (a) unsupervised PCA conducted to evaluate differences among different groups; (b) hierarchical clustering tree of samples from different groups; (c) volcano plot displaying the relative content difference of metabolites in the CCl4 vs. CCl4+DCA group; (d) comparison of plasma BAs levels between CCl4 and CCl4+DCA groups; (e) volcano plot displaying the relative content differences of metabolites in the CCl4+DCA and CCl4+DCA+Abx groups; (f) comparison of plasma BAs levels between CCl4+DCA and CCl4+DCA+Abx groups. Ctrl: Control; CD: CCl4+DCA; CDA: CCl4+DCA+Abx. n = 6, *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001. “ns” indicates no significant difference.

Figure 8. DCA regulates plasma bile acid profiles through gut microbiota (a) unsupervised PCA conducted to evaluate differences among different groups; (b) hierarchical clustering tree of samples from different groups; (c) volcano plot displaying the relative content difference of metabolites in the CCl4 vs. CCl4+DCA group; (d) comparison of plasma BAs levels between CCl4 and CCl4+DCA groups; (e) volcano plot displaying the relative content differences of metabolites in the CCl4+DCA and CCl4+DCA+Abx groups; (f) comparison of plasma BAs levels between CCl4+DCA and CCl4+DCA+Abx groups. Ctrl: Control; CD: CCl4+DCA; CDA: CCl4+DCA+Abx. n = 6, *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001. “ns” indicates no significant difference.

Figure 9. Liver transcriptome sequencing and qPCR analysis (a) cluster heatmap illustrating the gene expression in CCl4 and CCl4+DCA groups, with red indicating high-expression genes and blue indicating low-expression genes; (b) volcano plot displaying the relative content and statistical differences of genes in the CCl4 and CCl4+DCA groups; (c) gene function prediction based on the KEGG database; (d) GSEA analysis of IL-17 and TNF pathways; (e) real-time fluorescence quantitative PCR used to detect the expression of genes related to the IL-17 and TNF pathways.

Figure 9. Liver transcriptome sequencing and qPCR analysis (a) cluster heatmap illustrating the gene expression in CCl4 and CCl4+DCA groups, with red indicating high-expression genes and blue indicating low-expression genes; (b) volcano plot displaying the relative content and statistical differences of genes in the CCl4 and CCl4+DCA groups; (c) gene function prediction based on the KEGG database; (d) GSEA analysis of IL-17 and TNF pathways; (e) real-time fluorescence quantitative PCR used to detect the expression of genes related to the IL-17 and TNF pathways.

Figure 10. Multi-omics analysis (a) correlation heatmap of differential gut microbiota and differential BAs metabolites at the genus and species levels between CCl4 and CCl4+DCA groups; (b) correlation analysis between differential BAs metabolites and differential genes of IL-17 and TNF inflammatory signaling pathway. (c) O2PLS analysis results for metabolome and transcriptome. Bar heights reflect correlation magnitude with respective omics, darker colors signifying stronger correlations; (d) correlation analysis of differentially expressed gut microbiota and BAs metabolites between CCl4 and CCl4+DCA groups, as well as differentially expressed genes involved in IL-17 and TNF inflammatory pathways. The width corresponds to Mantel’s r statistic, representing the correlation strength, and the color represents the statistical significance of Mantel analysis.

Figure 10. Multi-omics analysis (a) correlation heatmap of differential gut microbiota and differential BAs metabolites at the genus and species levels between CCl4 and CCl4+DCA groups; (b) correlation analysis between differential BAs metabolites and differential genes of IL-17 and TNF inflammatory signaling pathway. (c) O2PLS analysis results for metabolome and transcriptome. Bar heights reflect correlation magnitude with respective omics, darker colors signifying stronger correlations; (d) correlation analysis of differentially expressed gut microbiota and BAs metabolites between CCl4 and CCl4+DCA groups, as well as differentially expressed genes involved in IL-17 and TNF inflammatory pathways. The width corresponds to Mantel’s r statistic, representing the correlation strength, and the color represents the statistical significance of Mantel analysis.
Supplemental material

Supplemental Material

Download Zip (2.8 MB)

Data availability statement

The datasets that support the findings of the study are available from the corresponding author upon reasonable request. The 16S rRNA gene sequencing data from this study are available in the NCBI Sequence Read Archive (SRA) repository under the accession number PRJNA1031407.