The systems biological effects of gut microbial metabolism are excellently reviewed elsewhere [
35], and therefore, in this review the emphasis will be in relation to NAFLD. It is likely that the severity of NAFLD affects differently the functions of the gut microbiota and vice versa. Therefore, the reviewed studies are divided below in subsections by the disease severity (steatotic without diagnosed NAFLD, NAFL, NASH, fibrosis, cirrhosis, hepatocellular carcinoma). Ultimately, also studies in children with NAFLD are reviewed. The reviewed studies are listed in
Table 1, which is primarily divided by the sections below, and by the order of appearance of the references in the text.
4.1. Fecal Metabolomics and Metagenomics Identifying the Gut Microbial Signatures in Steatotic Adults without Diagnosed NAFLD
The largest metagenomic study conducted by far in fatty liver disease included a sample of 6,269 Finnish participants of which most were overweight, but according to the BMI there were also normal weight individuals included [
38]. It should be noted that the functional metagenomic profiling was qualitative and not quantitative because the sequencing depth did not allow assembling contigs. In addition, instead of measuring liver fat content with imaging the authors used FLI to categorize the participants into liver fat content groups. Nevertheless, FLI is rather widely used and accepted index for NAFLD and its’ stratification [
39]. Ethanol and the SCFA acetate production pathways were found to be positively associate with FLI [
38]. Previously, it has been shown, for instance, that high-alcohol-producing strains of Klebsiella pneumoniae exist in humans with NAFLD [
40]. Thus, it seems that endogenously produced alcohol may play a role in hepatic fat accumulation at least in some populations. Generally, SCFA are considered to exert beneficial functions in the host, such as modulating immune functions [
41] and gastrointestinal permeability [
42]. However, they may also negatively impact inflammatory status of the host [
36].
Table 1.
The reviewed studies presenting metabolomic and/or metagenomic findings on the human gut microbiota and NAFLD.
Table 1.
The reviewed studies presenting metabolomic and/or metagenomic findings on the human gut microbiota and NAFLD.
Reference Number in the Text |
Analysis Method |
Study Patients/Adults |
Ruuskanen et al. [38] |
metagenomics |
fatty liver & all normal weight, overweight or obese |
Hoyles et al. [43] |
metagenomics |
fatty liver, morbidly obese |
Driuchina et al. [47] |
metabolomics |
fatty liver & healthy, obese |
Ge at al. [53] |
metabolomics |
NAFLD & healthy, overweight |
Boursier et al. [57] |
predicted metagenomics |
NAFLD & NAFLD + fibrosis & NASH |
Loomba et al. [58] |
metagenomics |
NAFLD + fibrosis |
Lee et al. [61] |
metabolomics |
NAFLD + obesity + fibrosis & NAFLD + normal weight + fibrosis |
Smirnova et al. [63] |
metabolomics |
NAFL, NASH, fibrosis regardless of BMI |
Sui et al. |
metabolomics |
NASH & healthy, normal weight |
Behary et al. [67] |
metabolomics, metagenomics |
NAFLD-HCC & healthy |
Oh et al. [68] |
metabolomics, metagenomics |
NAFLD-cirrhosis & healthy |
Reference number in the text |
Analysis method |
Study patients/children |
Michail et al. [49] |
metabolomics, metagenomics |
NAFLD + obese, healthy + obese & healthy + normal weight |
Testerman et al. [71] |
metagenomics |
NAFLD & healthy, obese |
Zhao et al. [72] |
metagenomics |
NAFLD & healthy, obese & normal weight |
Kordy et al. [60] |
metabolomics, metagenomics |
NAFL, NASH & healthy, obese & normal weight |
Yu et al. [73] |
metabolomics, metagenomics |
NAFLD & healthy, normal weight |
Schwimmer et al. [74] |
metagenomics |
NAFL, NASH & healthy, overweight |
Fecal metagenomic signatures have been described in non-diabetic morbidly obese women with hepatic steatosis [
43]. By mapping the microbial genes into Kyoto Encyclopedia of Genes and Genomes (KEGG) modules, it was found that steatosis associated positively with carbohydrate, lipid and amino acid metabolism as well as lipopolysaccharide (LPS) and peptidoglycan biosynthesis. LPS are recognized by Toll-like receptor 4 (TLR4), which expression has been shown to be increased in the liver of obese patients with NASH [
44]. Peptidoglycan, in turn, is recognized by multiple pattern-recognition receptors, including nucleotide-binding oligomerization domain-containing proteins (NODs), domain-containing 3 (NLRP3) and Toll-like receptor 2 (TLR2) [
45]. In addition to the above-mentioned pathways, Hoyles et al. also observed that hepatic steatosis was associated with increased number of genes related to the biosynthesis of branched-chain amino acids (BCAA) and aromatic amino acids [
43]. The metagenomic findings of feces were supported further by the elevated concentrations of these particular amino acids in plasma. To show mechanistically that the microbiome contributes to the hepatic steatosis, the authors transplanted feces from steatotic human donors to mice. Indeed, the fecal transplants caused rapid hepatic fat accumulation in mice, which involved elevated circulating BCAA [
43]. These observations are interesting as numerous studies have linked BCAA to obesity and NAFLD (for review [
46]).
We recently compared fecal and plasma metabolomes of humans with low (<5% of fat in liver) and high (>5% of fat in liver) liver fat content without diagnosed NAFLD [
47]. The study groups did not differ from each other in BMI. We found that the fecal histidine metabolism product, N-omega-acetylhistamine was markedly increased in individuals with fatty liver. In addition, another product of histidine degradation, anserine, positively associated with liver fat content [
47]. Previously, plasma levels of histidine have been shown to associate with the grade of hepatic steatosis [
48]. In agreement with a previous metagenome study [
49], we found that the levels of lysine degradation product, saccharopine, were higher in the feces of individuals with high liver fat content [
47]. As an indication of decreased steroid metabolism, the fecal levels of 6-hydroxybetatestosterone were reduced in the steatotic individuals. Contrarily, a previous study reported that low serum testosterone levels were associated with hepatic steatosis in obese male [
50].
To summarize, the two metagenomic studies in individuals with hepatic steatosis show different microbial signatures, which may be due to the different obesity status of the study populations. Neither of the metagenomic studies determined dietary factors between the groups or used them as confounding factors [
38,
43]. In our metabolomic study, there were no major dietary differences between the high and low liver fat content groups, except that higher vitamin E and sucrose intake was observed in individuals with fatty liver [
47]. Interestingly, sucrose is known to contribute to the onset of NAFLD [
51], while vitamin E is considered as a possible treatment for NAFLD [
52]. The main findings of the three reviewed studies in this section are presented in
Figure 2.
4.2. Fecal Metabolomics and Metagenomics Identifying the Gut Microbial Signatures in Patients with NAFL, NASH and Hepatic Fibrosis
Ge and co-workers have studied fecal metabolomics in patients with NAFLD and healthy controls using UHPLC/MS/MS [
53]. All participants were overweight and/or had visceral obesity. Combining the metabolite identification with KEGG pathway analysis, the authors found that the metabolism of nicotinate, nicotinamide and pyrimidine as well as signaling pathways of calcium and oxytocin and pancreatic secretion were altered in NAFLD patients compared to healthy controls. Interestingly, nicotinamide adenine dinucleotide (NAD
+) is being studied in clinical trials as a potential target to treat NAFLD [
54]. However, Ge et al. failed to show differences in the fecal levels of nicotinate between the healthy and NAFLD groups [
53]. Altered pyrimidine metabolism was reflected in lower fecal levels of uracil in the NAFLD patients. In addition, the catabolic byproduct of purine metabolism, xanthine was also lower likely due to the lower enzymatic activity of xanthine oxidase. The authors concluded that xanthine together with the abundance of specific microbial taxa might serve in the diagnostics of NAFLD [
53]. Intriguingly, it has been shown that inhibition of xanthine oxidase can ameliorate hepatic steatosis in mice [
55]. This might be related to decreased production of reactive oxygen species (ROS) [
56].
Based on literature, it seems that microbial functional genes can differentiate less advanced NAFLD from NASH and NAFLD with significant fibrosis [
57] and the stage of fibrosis [
58]. Yet, the authors of one paper did not perform shotgun metagenome sequencing but predicted the functional potential of the gut microbiota from 16S rRNA gene data using PICRUSt [
57]. The name is an abbreviation for Phylogenetic Investigation of Communities by Reconstruction of Unobserved States. PICRUSt is a bioinformatics software package designed to predict metagenome functional content from marker gene (e.g., 16S rRNA gene) surveys and full genomes using a developed algorithm developed by Languille et al. [
59]. However, the enriched functional categories in NASH were mostly related to carbohydrate and lipid metabolism [
57], which is similar to what has been found in simple steatosis [
43] and in NASH [
60] by others. Interestingly, patients with fibrosis could be distinguished from the patients without fibrosis by their microbial functional genes. Fibrosis was associated with enriched functional categories related to carbohydrate and lipid metabolism [
57]. A study that combined real microbial metagenomics and analysis of plasma metabolites in obese NAFLD patients [
58] did not find similar microbial signatures in fibrosis as PICRUSt did [
57]. Of the set of 89 metabolites that are produced both by the host and the microbiota, eleven metabolites could differentiate between mild/moderate NAFLD and advanced fibrosis. These enriched genes/metabolites were mainly related to carbon metabolism in fibrosis grade 2 as well as to nucleotide and steroid degradation in fibrosis grade 1. In addition, higher SCFA butyrate were annotated in grade 1, while the SCFA butyrate and propionate were higher in grade 2. Many findings from plasma metabolites and fecal metagenomic pathways supported each other [
58]. Contrary to these findings, we have found using fecal metabolomics that testosterone metabolism in the gut was lower in obese individuals with fatty liver [
47].
In addition to the stage of hepatic fibrosis, obesity status of the patients appears to be linked to the distinct gut microbial features in NAFLD [
61]. Lee and co-workers reported that while not seen in obese individuals with NAFLD, the levels of SCFA acetate and propionate increased along with the stage of fibrosis in non-obese individuals [
61]. This again indicates that the SCFA are not necessarily always beneficial for the host’s health. Lee et al. further found that the fecal levels of several conjugated and unconjugated bile acids were higher in non-obese individuals with fibrosis, while in obese individuals the levels of total conjugated bile acids inversely associated with the severity of fibrosis [
61]. Supporting the findings of fecal metabolomics, the expression levels of microbial genes encoding for bile salt hydrolase and 7α-hydroxy-3-oxochol-4-en-24-oyl-CoA dehydrogenase were lower in the non-obese patients with fibrosis.
The contribution of the bile acids to the pathophysiology of NAFLD
via gut-liver axis is rather well established, yet the bile acids have been mostly analyzed in plasma and not in feces in NAFLD studies [
62]. However, contrary to the findings of Lee et al. [
61] reviewed above, Smirnova and co-workers repor that the fecal concentrations of several secondary bile acids were lower in obese NASH patients with advanced fibrosis [
63]. Yet, there is a difference between the two studies in the obesity status of the patients. In the study by Smirnova et al. especially taurine conjugated bile acids increased along the stage of fibrosis [
63]. Interestingly, when the authors looked at the NAFLD activity scores (NAS), the findings were opposite. Compared to the healthy controls, the NAFLD patients had higher fecal levels of secondary bile acids and expression of microbial genes involved in biotransformation of bile acids. Further, NASH patients had higher levels of conjugated bile acids than patients with NAFL [
63]. This, in turn, is contrary to another study by Sui et al. showing higher levels of primary bile acids, chenodeoxycholic acid and cholic acid, in the feces of non-diabetic individuals with NASH compared to healthy controls [
64]. Nevertheless, these patients were normal weight, and therefore, the obesity status may be related to the different results between the reviewed studies. Of the other fecal metabolites, increased kynurenine and decreased L-tryptophan levels were good predictors of hepatic steatosis [
64]. The latter finding is interesting as a previous study had shown that tryptophan-derived metabolite kynurenine causes hepatic steatosis in mice by activating Aryl hydrocarbon receptor signaling [
65].
To summarize, while fecal bile acids seem to play an important role in NAFLD and fibrosis, the reviewed studies present different findings, which may, at least partly, depend on the obesity status of the study populations. None of the studies reported whether there were differences in dietary intakes between the study groups, though one study [
64] had collected food consumption data. The main findings on the microbial metagenomes and fecal metabolomics of the studies reviewed in this section are presented in
Figure 3.
4.3. Fecal Metabolomics and Metagenomics Identifying the Gut Microbial Signatures in Patients with NAFLD-Cirrhosis and NAFLD-Hepatocellular Carcinoma
With the global rise in the incidence of obesity and type 2 diabetes, NAFLD-related hepatocellular carcinoma (HCC) and NAFLD-cirrhosis are also becoming most common liver diseases [
66]. By using metagenomics, Behary et al. show that the gut microbiome of the NAFLD-HCC patients was characterized by higher number of SCFA synthesizing genes [
67], which is in line with what was found related to an earlier stage of the disease, namely higher FLI [
38]. Compared to the NAFLD-cirrhosis patients and healthy controls, in the feces of the NAFLD-HCC patients genes related to acetate (phosphate acetyltransferase,
pta), butyrate (phosphate butyryltransferase,
ptb) and propionate (fumarate reductase,
frd and succinate-CoA synthetase,
scs) synthesis were over-expressed [
67]. The gene expression findings were further confirmed using targeted LC/MS/MS and
1H-NMR quantification of the fecal metabolites. The levels of acetate and butyrate as well as oxaloacetate and acetylphosphate that are intermediates of the SCFA metabolism, were higher in the feces of the NAFLD-HCC patients, while propionate concentration did not differ between the study groups [
67]. In the future, it would be interesting to investigate, which other metabolites characterize the NAFLD-HCC patients using a non-targeted metabolomic profiling.
In NAFLD-cirrhosis, the gut microbial signatures have been studied by Oh and co-workers [
68]. Interestingly, similar to what has been shown in steatosis without diagnosed NAFLD [
43], BCAA and aromatic amino acids were predictors of cirrhosis as shown both by metagenomic and metabolomic analyses [
68]. In the NAFLD-cirrhotic patient’s fecal samples, the levels of L-tryptophan were increased due to decreased metabolism of it. This is contrary to another study showing that decreased L-tryptophan levels were good predictors of hepatic steatosis [
64].
To summarize the literature reviewed above, in NAFLD-associated cirrhosis and HCC some similar metabolic signatures can be found as in hepatic steatosis without diagnosed NAFLD, namely higher fecal levels of SCFA, BCAA and aromatic amino acids. Thus, in the future it would be important to explore their role more, and to consider whether these fecal metabolites would be suitable early biomarkers of advanced liver diseases before its’ onset. It should be noted, however, as in the other studies reviewed above that none of the studies in this section considered dietary factors as confounding between the study groups.
4.4. Fecal Metagenomics and Metabolomics Identifying the Gut Microbial Signatures in Children with NAFLD
A multi-omics study by Michail et al. combined metagenomic, metabolomic and proteomic approaches to study the gut microbiome of obese children with and without NAFLD [
49]. On a contrary to one targeted metabolomics study showing higher fecal acetate and propionate levels in NAFLD patients [
69], the
1H-NMR quantification of the fecal metabolites revealed that acetate was lower in pediatric NAFLD patients, while propionate was unaffected [
49]. As reviewed above, supporting again the role of endogenously produced ethanol in NAFLD, the fecal levels of ethanol were ~2-fold higher in obese with NAFLD than in the healthy children and obese children without NAFLD. Based on the microbial metagenomic analysis, the lysine degradation pathway was exclusively identified in children with NAFLD and not in healthy individuals [
49]. This agrees with our metabolomic findings from adults showing higher degradation products of lysine in the feces of individuals with high liver fat content [
47]. The metagenome shotgun sequencing of Michail et al. suggested that the microbial energy metabolism, including fatty acid and carbohydrate biosynthesis is much more enhanced in NAFLD pediatric patients compared to healthy children [
49]. However, this might be solely due to the different obesity status of the patients in the studies as it has been shown that an obese microbiome more efficiently harvests energy from the diet [
70].
Testerman and co-workers studied the fecal metagenomes of children with and without NAFLD [
71]. They observed that multiple pathways for lysine synthesis were increased, and histidine degradation was decreased in NAFLD patients. This is contrary to our recent metabolomics finding in adults with fatty liver showing higher fecal levels of lysine and histidine degradation products [
47]. However, in the NAFLD patient’s metagenomes also microbial metabolic pathways for BCAA, aromatic amino acids and peptidoglycan synthesis were enriched [
71]. Interestingly, similar findings have been reported in morbidly obese adults with hepatic steatosis [
43]. Thus, it could be that these metabolic pathways might help in the detection of the disease at an early age and/or stage of the disease.
Another metagenomic study compared healthy children, obese children, and obese children with NAFLD [
72]. Compared to healthy children, the obese pediatric patients with and without NAFLD had lower abundance of microbial genes related to pathways of replication and reparation, folding, sorting and metabolism of amino acids. The presence of NAFLD differentiated the obese groups showing enrichment of pathways related to digestive system, immune system, and glycan biosynthesis [
72]. Thus, in this case there were microbial signatures in NAFLD that were likely not dependent on the obesity status of the patients.
In agreement with the findings of Hoyles and co-workers in adults [
43], Kordy et al. report that compared to the healthy BMI-matched individuals, the microbiome of the pediatric patients with NASH was characterized by increased carbohydrate, lipid and amino acid metabolism as well as LPS biosynthesis [
60]. While the panels of plasma metabolites could accurately predict NASH, this was not seen with fecal metabolites analyzed from rectal swabs. The author of this current review opinions that the rectal swabs, even if stored at correct temperatures, may not be a reliable way to preserve the fecal metabolites for the subsequent analyses.
A combination of fecal metagenomics and targeted UPLC-MS/MS identified alterations in bile acid metabolism in children with NAFLD [
73]. The number of genes related to the biosynthesis of secondary bile acids were more abundant in NAFLD patients, which agrees with one study in adults [
63]. However, the comparison between the metagenomic and metabolomics results is a bit confusing as the fecal levels many secondary bile acids (a-hyodeoxycholic acid, 7-ketolithocholic acid, 23-nordeoxycholic acid, 7,12-diketolithocholic acid, 3-epideoxycholic acid, and dehydrocholic acid) were reduced in NAFLD, while only chenodeoxycholic acid-3-b-D-glucuronide was increased [
73]. Unfortunately, the authors did not discuss these discrepant results of their publication.
One study in children with NAFL and NASH used metagenomic shotgun sequencing to reveal microbial signatures related to the disease [
74]. Compared to the BMI-matched healthy controls, microbial LPS biosynthesis was significantly enriched in children with NAFL and NASH. This finding agrees with what has been found in non-diabetic obese women with hepatic steatosis [
43]. Several genes related to flagellar assembly were also enriched in children with NASH and in those with moderate-to-severe fibrosis [
74]. We have previously shown in an a mouse model that gut microbial flagellin causes hepatic fat accumulation, which is mediated by Vascular adhesion protein-1 [
75]. Altogether these results suggest that flagellin may promote NAFLD, while other studies have suggested that a knockout of flagellin-recognizing TLR5 protects mice from NAFLD [
36]. However, it should be noted that flagellin can also act via cytosolic nucleotide oligomerization domain (NOD)-like receptors to affect the inflammatory status [
76]. Thus, the effects of flagellin and its’ receptors in relation to the onset of NAFLD should be explored more in the future.
To summarize, studies in children with NAFLD have reported microbial signatures of which some are similar (increased LPS biosynthesis and amino acid metabolism) and some distinct from each other. However, a few of the reviewed findings are like those seen in adults, such as higher fecal levels of ethanol in individuals with hepatic steatosis and increased LPS biosynthesis in NAFLD. None of the reviewed studies in children analyzed whether there were dietary differences between the study groups. The main microbial metagenomic and fecal metabolomic findings from children with NAFLD are summarized in
Figure 4.