Skip to main content

Characterization of microbiome and metabolite analyses in patients with metabolic associated fatty liver disease and type II diabetes mellitus



State-of-the-art renewal has indicated the improvement of diagnostics of patients with metabolic associated fatty liver disease (MAFLD) and/or type II diabetes mellitus (T2DM) by dissecting the clinical characteristics as well as genomic analysis. However, the deficiency of the characterization of microbial and metabolite signatures largely impedes the symptomatic treatment.


For the purpose, we retrospectively analyzed the clinical data of 20 patients with MAFLD (short for “M”), 20 cases with MAFLD and T2DM (short for “MD”), together with 19 healthy donors (short for “Ctr”). Microbial and metabolite analyses were further conducted to explore the similarities and differences among the aforementioned populations based on feces and blood samples, respectively.


Compared with those in the Ctr group, patients with M or MD revealed multifaceted similarities (e.g., Age, ALP, LDL, BUN) and distinctions in clinical indicators of liver (e.g., BMI, ALT, PCHE, CAP). With the aid of microbial and metabolite analyses as well as bioinformatic analyses, we found that the characteristics of gut microbiota (e.g., abundance, hierarchical clustering, cladogram, species) and lipid metabolism (e.g., metabolite, correlation coefficient and scatter plot) were distinct among the indicated groups.


The patients with MD revealed multifaceted similarities and distinctions in characteristics of microbiome and metabolites with those in the M and HD groups, and in particular, the significantly expressed microbes (e.g., Elusimicrobiota, Berkelbacteria, Cyanobacteria, Peregrinibacteria) and lipid metabolites (e.g., Lipid-Q-P-0765, Lipid-Q-P-0216, Lipid-Q-P-0034, Lipid-Q-P-0800), which would collectively benefit the clinical diagnosis of MAFLD and T2DM.

Peer Review reports


Metabolic associated fatty liver disease (MAFLD) also known as non-alcoholic fatty liver disease (NAFLD), including nonalcoholic steatohepatitis (NASH) and nonalcoholic fatty liver (NAFL), has been recognized as the leading cause of chronic liver disorder worldwide and the burgeoning public health issue [1,2,3]. Definitely, MAFLD is a spectrum of liver disorders with metabolic dysfunction and the accompanied clinical presence of steatosis in over 5% of hepatocytes, which commonly result in poor hepatic manifestation and prognosis such as cirrhosis, fibrosis and hepatocellular carcinoma (HCC) [2, 4, 5].

Longitudinal studies have suggested the continuous improvement in the diagnostics and management of MAFLD (e.g., liver biopsy, clinical indicators, therapeutic intervention), yet the patients still endure the well-known challenges such as long-suffering relapse, metabolic syndrome and immunological rejection largely attribute to the deficiency of parameter spectrum for dissecting the pathogenesis including gut microbiota and lipid metabolism [6, 7]. In particular, the design of clinical trials and the development of pharmacotherapies are far from satisfaction largely due to the current limitations in understanding the heterogeneity of MAFLD patients with multiple complications including essential hypertension (EH) and T2DM inaccuracies as well as the necessity of reappraising the terminology and definitions of the nomenclature [8,9,10].

Notably, state-of-the-art renewal has revealed diseases linked to dyslipidemia including atherosclerosis, obesity and MAFLD, which are associated with variations in gut microbiota profile and the accompanied host metabolism and physiology [11,12,13,14]. Meanwhile, compelling evidences have indicated that the gut microbiota might impact lipid metabolism in blood and tissue of the body via metabolites (e.g., short-chain fatty acids, trimethylamine, secondary bile acids, indole and its derivatives) and the bacterially proinflammatory factors (e.g., extracellular vesicles, peptidoglycan, DNA, lipopolysaccharide) [15,16,17]. However, the definitive and detailed associations among MAFLD, gut microbiota and lipid metabolism are still woefully inadequate [14, 18].

In this study, we took advantage of the human feces 16S amplicon technology and serum targeted quantification of lipid metabolism in combination with multiple clinical parameters to explore the multifaceted characteristics and candidate diagnostic biomarkers for dissecting patients with MAFLD and/or T2DM from health crowd. Collectively, our study suggested the similarities and distinctions in the clinical indicators, microbiome and metabolites, which would supply new references for the diagnosis and pathogenesis of MAFLD in future.

Methods and materials

Participants and specimen collection

During October and December of 2020, MAFLD and/or T2DM patients in Department of Gastroenterology and healthy volunteers in Physical Examination Center of the First Affiliated Hospital of Kunming Medical University were screened and grouped according to the inclusion and exclusion criteria as well as Declaration of Helsinki and the approval of the Ethics Committee of the First Affiliated Hospital of Kunming Medical University (Approval number: 2020-L-08; Approval date: January 20th, 2020). Informed consent was obtained from all participants (Ctr, MD and M) included in the study.

As to blood sample, 3 ml peripheral blood were collected with anticoagulant tubes containing EDTA as we recently reported. The serum should be isolated from blood samples within 30 min and stored in − 80 °C. As to microflora samples, over 500 mg faeces were collected with sterile centrifugal tubes and stored in − 80 °C as well. The detailed information of the participants was available in Additional file 3: Additional Table S1, and Additional file 4: Additional Table S2.

The inclusion and exclusion criteria

As mentioned above, among the 200 patients (100 M patients, 100 MD patients) and 100 healthy candidates, 20 patients with M or MD and 20 Ctr were enrolled in the study. Generally, the inclusion and exclusion criteria for MAFLD were according to the international expert consensus statement [4]. For instance, the criteria for MAFLD were based on the evidence of hepatic steatosis, and one of the following three criteria, including overweight/obesity, evidence of metabolic dysregulation, or presence of type 2 diabetes mellitus [4]. In details, the inclusion criteria were mainly based on the controlled attenuation parameter (CAP) value and liver stiffness measurement (LSM) value measured by fibroscan (502 model, Echosens, France). The M and MD patients with a CAP value over 238 dB/m, while the Ctr with that lower than 238 dB/m instead. Additionally, the participants should stop using motility drugs and laxatives at least 1 week prior to breath test.

The exclusion criteria were: (1) an alcohol consumption history, over 140 g/week for men and over 70 g/week for women; (2) other diseases that can lead to fatty liver, including viral hepatitis, drug-induced liver disease, total parenteral nutrition, Wilson’s disease and autoimmune liver disease; (3) Age < 45, pacemaker installed, unhealed wounds in right upper abdomen and ascites; (4) malignancy or other terminal diseases including serious liver, kidney, heart, brain diseases and malignant tumors; (5) Subjects who have taken antibiotics in the last 4 weeks, or those have received clean enema and colonoscopy in the last 2 weeks.

Microbiome and metabolite analyses

Faeces and peripheral blood samples were collected from the indicated groups (Ctr, MD, M), and then turned to 16S amplicon technology-based sequencing and gas/liquid chromatography-based serum targeted quantification of lipid metabolism for microbiome and metabolite analyses, respectively. The RNAs were sent to quality test and Novogene (Tianjin, China) for sequencing. The bioinformatic analyses of the data were completed by utilizing the software and platforms including Venn Map diagram, HeatMap, Principal Component Analysis (PCA) as we previously described [19,20,21].

Statistical analysis

All statistical analysis was performed as we recently reported [21,22,23]. The GraphPad Prism 6.0 (GraphPad Software, USA) software was adopted for statistical analysis. The data were shown as mean ± SD (N = 3 independent experiments) and the One-way ANOVA analyses were used for the comparison of the indicated groups. Only when P < 0.05 was considered statistically significant. NS, not significant; *, P < 0.05; **, P < 0.01; ***, P < 0.001.


FBG value revealed preferable correlation with MAFLD and/or T2DM patients

From October to December, 2020, a total number of 40 patients (20 M, 20 MD) and 19 healthy participants (Ctr) were enrolled from the 300 candidates including 100 M, 100 MD and 100 Ctr according to the inclusion and exclusion criteria (Additional file 3: Additional Table S1, Additional file 4: Additional Table S2). As shown by the statistical analyses, patients with M or MD revealed distinctions in multiple clinical parameters compared to those in the Ctr group, including body weight, body mass index (BMI), alanine aminotransferase (ALT), aspartate aminotransferase (AST), triglyceride (TG), uric acid (UA), Girth, pseudocholinesterase (PCHE), Hipline, fasting blood glucose (FBG), while no significant differences between the M and MD groups were observed in the aforementioned indicators except FBG (Fig. 1A-J, Additional file 2: Additional Fig. S1A-S1F, Additional file 1: Additional Information, Additional file 3: Additional Table S1).

Fig. 1
figure 1

The clinicopathological indicators involved in M and MD patients and Ctr. A-M The comparations of clinicopathological parameters among Ctr, M and MD, including Weight (A), BMI (B), ALT (C), AST (D), TG (E), UA (F), Girth (G), PCHE (H), Hipline (I), FBG (J), TP (K), r-GGT (L), ALB (M), CREA (N), and CAP (O). All data are shown as Mean ± SEM. *, P < 0.05; **, P < 0.01; ***, P < 0.001; NS, not significant

Compared to the Ctr group, the values of total protein (TP), γ-Glutamyl transferase (r-GGT) and albumin (ALB) in M patients rather than those in the MD cases exhibited statistical differences (Fig. 1K-M, Additional file 3: Additional Table S1). Furthermore, we also found a certain number of clinical parameters showed distinctions between M and MD patients such as FBG, ALB and creatinine (CREA) (Fig. 1J, M-N). As to other indicators involved in liver disease (e.g., HDL, ALP, BUN) and tumor (AFP, CA199, CEA) diagnosis, minimal differences were observed among the indicated three groups (Additional file 2: Additional Fig. S1B-S1N, Additional file 1: Additional Information). Notably, the controlled attenuation parameter (CAP) value between the Ctr and M, or Ctr and MD groups showed significant differences, which was consistent with our recent reports (Fig. 1O) [2]. Collectively, we verified that only the noninvasive FBG distribution basically satisfied the distinction of patients with MAFLD and/or T2DM from healthy participants.

Patients with MAFLD and T2DM manifested distinguishable profiles of microbiota

Having dissected the clinical relationship of multiple parameters with MAFLD and/or T2DM patients, we next turned to microbiome analysis for further exploring the candidate indicators by utilizing 16S amplicon technology-based sequencing. Generally, we found that the microbiomes in the indicated three groups revealed unique spectrum in spots of distribution and relative abundance as shown by the plots of principal components analysis (PCA) diagrams and species accumulation maps, respectively (Fig. 2A-B). Furthermore, the Venn Map and phylum diagrams intuitively exhibited the distribution of the 4560 kinds of common or distinct microbes among the Ctr, M and MD groups (Fig. 2C-D). For instance, a total of 1801 kind of gut microflora were reflected in the three groups, while only 511, 405 and 345 microflora were specifically enriched in the corresponding Ctr, M and MD groups, respectively (Fig. 2C). Interestingly, from the view of Krona genus evolutionary tree, we noticed that Bacteroideta, Proteobacteria and Firmicutes were the top 3 enriched microflora over those in the rest of phylum (Fig. 2D).

Fig. 2
figure 2

The comparation of microbiota profiles among the indicated groups. A The PCA (principal component analysis) of microbiota among the indicated groups (Ctr, M, MD) based on PC1 and PC2 values. B Histograms of microbiota-associated biological processes in the indicated groups. C Venn Map analysis of the distributions of the distribution of microbiota in stool samples among the indicated groups. D Krona genus evolutionary tree showed the phylum of the enriched microflora. E Hierarchical cluster analysis showed the correlation of the participants with the enriched biofunctions in the indicated groups (Ctr, M and MD). F-G Kyoto encyclopedia of genes and genomes (KEGG) showed the significant differences between Ctr and MD, or between MD and M

Hierarchical cluster analysis further exhibited the detailed correlation of the participants (Ctr, M and MD) with the enriched biofunctions in the indicated groups such as nucleotide metabolism, metabolism of terpenoids, folding and degradation (Fig. 2E). Additionally, within the 95% confidence intervals, a certain number of signaling pathways enriched by Kyoto encyclopedia of genes and genomes (KEGG) revealed significant differences between the enrolled healthy participants and MAFLD and/or T2DM patients (P < 0.05) (Fig. 2F-G).

Patients with MAFLD and T2DM showed multifaceted diversity of microbiota

After evaluating the overview of the profile, we next explored the detailed similarities and differences of microbiota in the aforementioned groups. As shown by the species accumulation boxplot, we intuitively observed the sharp rise of species in the range of 1 and 19, whereas limited increase in species in the range of 28 and 55 (Fig. 3A). According to the rank abundance curve, the species richness and evenness of bacteria community revealed minimal variations among the three groups (Fig. 3B). Simultaneously, with the aid of the LEfSe (LDA effect size) methods, we conducted analysis of different species between groups and identified a series of biomarkers as shown by the cladogram (Fig. 3C). For example, the cladogram of gut microflora with significant differences between the Ctr (Control) and M (MAFLD patients) groups were exhibited such as Prevotellaceae, Chloroplast, Cyanobacteriia, Ruminococcaceae, Oscillospirales and Clostriadia (Fig. 3C).

Fig. 3
figure 3

The multifaceted diversity of microbiota among the indicated groups. A Species accumulation boxplot showed the dynamic changes of species according to the number of samples. B The rank abundance curve showed the species richness and evenness of bacteria community among the indicated groups (Ctr, M, MD). C The cladogram of gut microflora among the indicated groups. D UPGMA (Unweighted Pair-group Method with Arithmetic Mean) analysis revealed the microbial communities and the relative species abundance at phylum level. E Random Forest analysis showed the taxonomy of microflora based on mean decrease accuracy (E) and mean decrease gin (F), respectively. G-H α diversity (G) and β (H) diversity index analyses revealed the between-group variance of species among the indicated groups

Therewith, we took advantage of the weighted unifrac distance matrix for the UPGMA (Unweighted Pair-group Method with Arithmetic Mean) analysis to explore the composition of microbial communities and the relative species abundance of each sample at phylum level in the indicated samples. Interestingly, microbiota within and between the aforementioned groups collectively manifested significant diversity in phylum, which suggested the complexity of bacteria the community and species among the individuals (Fig. 3D). Furthermore, by performing the random forest analysis, we respectively analyzed the top 20 species abundance, and verified that Coprococcus, Subdoligranulum and UCG_005 were the top 3 species based on the analysis of mean decrease accuracy, while Coprococcus, Subdoligranulum and Lachnospira were the top 3 species based on the analysis of mean decrease gin (Fig. 3E-F). Finally, with the aid of α and β diversity index analyses, we found that the between-group variance of species among the indicated groups showed moderate differences (Fig. 3G-H).

The diversity of the characteristics and correlations of 16S-based microbiota and lipid metabolism

To explore the correlation of microbiome and metabolite in the indicated groups, we further turned to 16S-based gut microbiota and lipid metabolism analyses. Compared to the microbiota distribution, all of the aforementioned groups (Ctr, M, MD) revealed distinguishable pattern of metabolite according to the PCA diagrams (Fig. 4A-B). Interestingly, the PC1 and PC2 based PCA analysis further reflected the detailed similarities and differences of lipid metabolism among the individuals of the enrolled participants on a series of levels including phylum, class, order, family, genus and species (Fig. 4C). Intuitively, as shown by the Sperman correlation-based hierarchical clustering analysis, the significantly different microflora among the three groups revealed diversity in the pattern of gut microbiota and lipid metabolism (Fig. 4D-F). For example, the Peregrinibacteria and representative lipid metabolite (e.g., Lipid-Q-P-0862, Lipid-Q-P-0794, Lipid-Q-P-0869) in the MD group showed significantly negative correlation coefficient compared to the Ctr group (P < 0.05), while Firmicutes and Lipid-Q-P-1156 between the Ctr and M groups and Sumerlaeota and Lipid-Q-P-0034 between the MD and M groups revealed positive and negative correlation coefficient, respectively (Fig. 4D-F).

Fig. 4
figure 4

The overview of characteristics and correlations of microbiota and lipid metabolism. A PCA diagrams of metabolite and microbiota in the stool samples based on PC1, PC2 and PC3 values among the indicated groups (Ctr, M, MD). B PCA diagrams of metabolite in the stool samples based on PC1 and PC2 values among the indicated groups. C PCA diagrams of lipid metabolism on a series of levels including phylum, class, order, family, genus and species. D-F Sperman correlation-based hierarchical clustering analysis revealed the diversity in gut microbiota and lipid metabolism between Ctr and MD (D), Ctr and M (E), MD and M (F). *, P < 0.05

Furthermore, the spatial correlation between the species of the indicated gut microbiome and metabolite was intuitive presented by the Circos diagrams among the Ctr, M and MD groups (Fig. 5A). Furthermore, the representative microbiome and metabolite with significant positive or negative correlations (P < 0.05, |r| ≥ 0.3) between the indicated groups were shown by the correlation scatter plots as well (Fig. 5B-D). For instance, Erysipelotrichaceaes had negative correlation with multiple lipid metabolism (e.g., Lipid-Q-P-179, Lipid-Q-P-0225, Lipid-Q-P-0794) between the Ctr and M groups, whereas Lipid-Q-P-0794 showed negative and positive correlations with Cyanobacteria and Cyanobacteria between the Ctr and MD groups, respectively (Fig. 5B-C). As to the M and MD groups, we also observed the diversity of the correlations between the representative microbiome (e.g., Cyanobacteria, Ktedonobacterales, Prevotellaceae) and the corresponding metabolite (e.g., Lipid-Q-P-0216, Lipid-Q-P-0765, Lipid-Q-N-0034) (Fig. 5D).

Fig. 5
figure 5

The detailed correlations of microbiota and lipid metabolism among the indicated groups. A Circos diagrams showed the correlation between the species of gut microbiome and metabolite in the indicated groups (Ctr, M, MD). B-D Correlation scatter plots showed the correlation of representative microbiome and metabolite between Ctr and MD (B), Ctr and M (C), MD and M (D)


The pathophysiology underlying MAFLD is a long-lasting issue, which comprises a multitude of interlinked processes such as insulin resistance, lipotoxicity, infiltration of proinflammatory cells, hepatic stellate cell (HSC) fibrogenesis and over-activation [24, 25]. Despite the aforementioned well-established processes, yet the molecular mechanisms are incompletely uncovered and the approved pharmacological treatments for MAFLD are also unavailable [1, 26, 27]. For the purpose of further verifying the detailed pathogenesis of patients with MAFLD and/or T2DM for developing novel treatment remedies, we took advantage of the 16S-based sequencing and serum targeted quantification of lipid metabolism for gut microbiome and metabolite analyses. Notably, we characterized the similarities and diversities of the microbiome and metabolite among health participants, MAFLD and/or T2DM patients, and verified the correlation of representative microbiome and the corresponding metabolite that held promising prospect as novel biomarkers and therapeutic targets in dissecting and resolving the MAFLD-associated pathophysiology.

Current studies have put forward novel concepts underlying pathophysiology of MAFLD [27, 28]. On the one hand, the excess extracellular matrix (ECM) and proinflammatory cytokines secreted by activated HSCs, proinflammatory cells and damaged hepatocytes have been demonstrated to play a critical role in impairing the physiological structure and biofunctions of the liver [24, 26, 29]. On the other hand, recent advances in microbiome and metabolite have also indicated the involvement of dysfunctional gut microbiome, dysregulated lipid metabolism, intestinal fructose handling and endogenous alcohol production, together with abnormal regulation and alterations of multiple signaling cascades (e.g., Hedgehog, BMP, TGF-β) for modifying the individual susceptibility to MAFLD patients with T2DM [24]. Additionally, the obesity as well as T2DM has been involved with MAFLD and microbiota from pathophysiology to therapeutics as well [30, 31]. Consistently, our study revealed the multifaceted similarities and diversities of MAFLD patients with T2DM from the perspectives of gut microbiome and metabolites as well as clinical indicators.

Emerging evidences have suggested the metabolic disorders with various alterations in the composition of intestinal microbiota and relevant metabolites, which are adequate to translocate from the gut across to liver and thus contribute to metabolic liver disorders (e.g., MAFLD) or T2DM [14, 32, 33]. Patients with T2DM manifested features of insulin resistance including beta cell deterioration, augmented endogenous glucose production and decreased peripheral glucose uptake, which has been considered an integral pathogenesis for MAFLD development [34, 35].

The gut-liver axis is the consequence of bidirectional regulation between the anatomical structure of gut and liver, which thus determines the pivotal modulating effects of gut microbiota and the concomitant microbiome upon liver activity [36, 37]. To date, microbiome has been indicated with the involvement of the pathogenesis of MAFLD or obesity via generating harmful metabolites (e.g., short-chain fatty acids, secondary bile acids, branched amino acids, and indoles) [38, 39]. Very recently, Cerreto and their colleagues proposed the characteristics of the gut-liver axis in liver disorders, and Fianchi et al further confirmed the gut-liver interaction in the pathophysiology of MAFLD and the potential target for personalized treatment [40,41,42]. Furthermore, in this study, we conducted the gut microbiota and lipid metabolite analyses, and suggested the multifaceted characteristics and inherent correlations with pathophysiologic abnormalities among the healthy population and those patients with MAFLD and/or T2DM, which collectively supply new references for further exploring the gut microbiota-mediated pathogenesis as well as gut microbiome- or metabolite-targeted therapies on NAFLD in future.


Gut microbiota and lipid metabolites in MAFLD and/or T2DM hold promising prospective for developing more accurate diagnosis and novel pharmacological treatments. Our findings based on microbiota and metabolite analyses would provide overwhelming new references for understanding the pathogenesis, diagnosis as we as facilitating the development of gut microbiome- or metabolite-based biomarkers and targeted therapies on MAFLD.

Availability of data and materials

All data during this study are included in the published article. Meanwhile, the datasets analyzed during the study are available from the corresponding author upon reasonable request.



Metabolic associated fatty liver disease


Type II diabetes mellitus


Non-alcoholic fatty liver disease


Nonalcoholic steatohepatitis


Nonalcoholic fatty liver


γ-glutamyl transferase


Hepatocellular carcinoma


Controlled attenuation parameter


Liver stiffness measurement


Principal Component Analysis


Body mass index


Alkaline phosphatase


Alanine aminotransferase


Aspartate aminotransferase


Fasting blood glucose




Uric acid


Hepatic stellate cell


Extracellular matrix


Blood urea nitrogen


Total bilirubin


High-density lipoprotein-cholesterol


Low-density lipoprotein-cholesterol




Total cholesterol




Platelet count


  1. Barchetta I, Cimini FA, Cavallo MG. Vitamin D and metabolic dysfunction-associated fatty liver disease (MAFLD): an update. Nutrients. 2020;12(11):3302.

    Article  CAS  PubMed Central  Google Scholar 

  2. Yang H, Li Q, Zhang L, Zhu M, Niu J, Xue F, et al. LncPRYP4-3 serves as a novel diagnostic biomarker for dissecting subtypes of metabolic associated fatty liver disease by targeting RPS4Y2. Clin Exp Med. 2020;20(4):587–600.

    Article  CAS  PubMed  Google Scholar 

  3. Jiang C, Xie C, Li F, Zhang L, Nichols RG, Krausz KW, et al. Intestinal farnesoid X receptor signaling promotes nonalcoholic fatty liver disease. J Clin Invest. 2015;125(1):386–402.

    Article  PubMed  Google Scholar 

  4. Eslam M, Newsome PN, Sarin SK, Anstee QM, Targher G, Romero-Gomez M, et al. A new definition for metabolic dysfunction-associated fatty liver disease: an international expert consensus statement. J Hepatol. 2020;73(1):202–9.

    Article  PubMed  Google Scholar 

  5. Yao J, Chen N, Wang X, Zhang L, Huo J, Chi Y, et al. Human supernumerary teeth-derived apical papillary stem cells possess preferable characteristics and efficacy on hepatic fibrosis in mice. Stem Cells Int. 2020;2020:6489396.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Katsiki N, Mikhailidis DP, Mantzoros CS. Non-alcoholic fatty liver disease and dyslipidemia: an update. Metabolism. 2016;65(8):1109–23.

    Article  CAS  PubMed  Google Scholar 

  7. Mato JM, Alonso C, Noureddin M, Lu SC. Biomarkers and subtypes of deranged lipid metabolism in non-alcoholic fatty liver disease. World J Gastroenterol. 2019;25(24):3009–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Eslam M, Sanyal AJ, George J, International Consensus P. MAFLD: a consensus-driven proposed nomenclature for metabolic associated fatty liver disease. Gastroenterology. 2020;158(7):1999–2014 e1991.

    Article  CAS  PubMed  Google Scholar 

  9. Singh SP, Anirvan P, Reddy KR, Conjeevaram HS, Marchesini G, Rinella ME, et al. Non-alcoholic fatty liver disease: not time for an obituary just yet! J Hepatol. 2021;74(4):972–4.

    Article  PubMed  Google Scholar 

  10. Tilg H, Effenberger M. From NAFLD to MAFLD: when pathophysiology succeeds. Nat Rev Gastroenterol Hepatol. 2020;17(7):387–8.

    Article  PubMed  Google Scholar 

  11. Schoeler M, Caesar R. Dietary lipids, gut microbiota and lipid metabolism. Rev Endocr Metab Disord. 2019;20(4):461–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Jonsson AL, Backhed F. Role of gut microbiota in atherosclerosis. Nat Rev Cardiol. 2017;14(2):79–87.

    Article  CAS  PubMed  Google Scholar 

  13. Xi L, Qin X, Song Y, Han J, Li Z, Zhang J. Gut microbial alterations in diarrheal Baer's Pochards (Aythya baeri). Front Vet Sci. 2021;8:756486.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Canfora EE, Meex RCR, Venema K, Blaak EE. Gut microbial metabolites in obesity, NAFLD and T2DM. Nat Rev Endocrinol. 2019;15(5):261–73.

    Article  CAS  PubMed  Google Scholar 

  15. Chen ML, Yi L, Zhang Y, Zhou X, Ran L, Yang J, et al. Resveratrol attenuates Trimethylamine-N-oxide (TMAO)-induced atherosclerosis by regulating TMAO synthesis and bile acid metabolism via remodeling of the gut microbiota. mBio. 2016;7(2):e02210–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wang F, Xu J, Jakovlic I, Wang WM, Zhao YH. Dietary betaine reduces liver lipid accumulation via improvement of bile acid and trimethylamine-N-oxide metabolism in blunt-snout bream. Food Funct. 2019;10(10):6675–89.

    Article  CAS  PubMed  Google Scholar 

  17. Lee G, You HJ, Bajaj JS, Joo SK, Yu J, Park S, et al. Distinct signatures of gut microbiome and metabolites associated with significant fibrosis in non-obese NAFLD. Nat Commun. 2020;11(1):4982.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Albillos A, de Gottardi A, Rescigno M. The gut-liver axis in liver disease: pathophysiological basis for therapy. J Hepatol. 2020;72(3):558–77.

    Article  CAS  PubMed  Google Scholar 

  19. Huo J, Zhang L, Ren X, Li C, Li X, Dong P, et al. Multifaceted characterization of the signatures and efficacy of mesenchymal stem/stromal cells in acquired aplastic anemia. Stem Cell Res Ther. 2020;11(1):59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Wei Y, Zhang L, Chi Y, Ren X, Gao Y, Song B, et al. High-efficient generation of VCAM-1(+) mesenchymal stem cells with multidimensional superiorities in signatures and efficacy on aplastic anaemia mice. Cell Prolif. 2020;53(8):e12862.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Zhang L, Wei Y, Chi Y, Liu D, Yang S, Han Z, et al. Two-step generation of mesenchymal stem/stromal cells from human pluripotent stem cells with reinforced efficacy upon osteoarthritis rabbits by HA hydrogel. Cell Biosci. 2021;11(1):6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Zhang L, Chi Y, Wei Y, Zhang W, Wang F, Zhang L, et al. Bone marrow-derived mesenchymal stem/stromal cells in patients with acute myeloid leukemia reveal transcriptome alterations and deficiency in cellular vitality. Stem Cell Res Ther. 2021;12(1):365.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Zhang L, Wang H, Liu C, Wu Q, Su P, Wu D, et al. MSX2 initiates and accelerates Mesenchymal stem/stromal cell specification of hPSCs by regulating TWIST1 and PRAME. Stem Cell Rep. 2018;11(2):497–513.

    Article  CAS  Google Scholar 

  24. Kuchay MS, Choudhary NS, Mishra SK. Pathophysiological mechanisms underlying MAFLD. Diabetes Metab Syndr. 2020;14(6):1875–87.

    Article  PubMed  Google Scholar 

  25. Zheng KI, Sun DQ, Jin Y, Zhu PW, Zheng MH. Clinical utility of the MAFLD definition. J Hepatol. 2021;74(4):989–91.

    Article  PubMed  Google Scholar 

  26. Xian YX, Weng JP, Xu F. MAFLD vs. NAFLD: shared features and potential changes in epidemiology, pathophysiology, diagnosis, and pharmacotherapy. Chin Med J. 2020;134(1):8–19.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Rispo A, Imperatore N, Guarino M, Tortora R, Alisi A, Cossiga V, et al. Metabolic-associated fatty liver disease (MAFLD) in coeliac disease. Liver Int. 2021;41(4):788–98.

    Article  CAS  PubMed  Google Scholar 

  28. Eslam M, George J. MAFLD: a holistic view to redefining fatty liver disease. J Hepatol. 2021;74(4):983–5.

    Article  PubMed  Google Scholar 

  29. Huang J, Kumar R, Wang M, Zhu Y, Lin S. MAFLD criteria overlooks a number of patients with severe steatosis: is it clinically relevant? J Hepatol. 2020;73(5):1265–7.

    Article  PubMed  Google Scholar 

  30. Polyzos SA, Kountouras J, Mantzoros CS. Obesity and nonalcoholic fatty liver disease: from pathophysiology to therapeutics. Metabolism. 2019;92:82–97.

    Article  CAS  PubMed  Google Scholar 

  31. Lau LHS, Wong SH. Microbiota, obesity and NAFLD. Adv Exp Med Biol. 2018;1061:111–25.

    Article  CAS  PubMed  Google Scholar 

  32. Tilg H, Zmora N, Adolph TE, Elinav E. The intestinal microbiota fuelling metabolic inflammation. Nat Rev Immunol. 2020;20(1):40–54.

    Article  CAS  PubMed  Google Scholar 

  33. Amar J, Chabo C, Waget A, Klopp P, Vachoux C, Bermudez-Humaran LG, et al. Intestinal mucosal adherence and translocation of commensal bacteria at the early onset of type 2 diabetes: molecular mechanisms and probiotic treatment. EMBO Mol Med. 2011;3(9):559–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Solis-Herrera C, Triplitt C, Cersosimo E, DeFronzo RA. Pathogenesis of type 2 diabetes mellitus. 2021 Sep 27. In: Feingold KR, Anawalt B, Boyce A, Chrousos G, de Herder WW, Dhatariya K, Dungan K, Hershman JM, Hofland J, Kalra S, Kaltsas G, Koch C, Kopp P, Korbonits M, Kovacs CS, Kuohung W, Laferrère B, Levy M, McGee EA, McLachlan R, Morley JE, New M, Purnell J, Sahay R, Singer F, Sperling MA, Stratakis CA, Trence DL, Wilson DP, editors. Endotext [Internet]. South Dartmouth:, Inc.; 2000–. PMID: 25905339.

  35. Mazzini FN, Cook F, Gounarides J, Marciano S, Haddad L, Tamaroff AJ, et al. Plasma and stool metabolomics to identify microbiota derived-biomarkers of metabolic dysfunction-associated fatty liver disease: effect of PNPLA3 genotype. Metabolomics. 2021;17(7):58.

    Article  CAS  PubMed  Google Scholar 

  36. Son G, Kremer M, Hines IN. Contribution of gut bacteria to liver pathobiology. Gastroenterol Res Pract. 2010;2010:453563.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Agus A, Planchais J, Sokol H. Gut microbiota regulation of tryptophan metabolism in health and disease. Cell Host Microbe. 2018;23(6):716–24.

    Article  CAS  PubMed  Google Scholar 

  38. Wu J, Wang K, Wang X, Pang Y, Jiang C. The role of the gut microbiome and its metabolites in metabolic diseases. Protein Cell. 2021;12(5):360–73.

    Article  PubMed  Google Scholar 

  39. Lee CJ, Sears CL, Maruthur N. Gut microbiome and its role in obesity and insulin resistance. Ann N Y Acad Sci. 2020;1461(1):37–52.

    Article  PubMed  Google Scholar 

  40. Fianchi F, Liguori A, Gasbarrini A, Grieco A, Miele L. Nonalcoholic fatty liver disease (NAFLD) as model of gut-liver Axis interaction: from pathophysiology to potential target of treatment for personalized therapy. Int J Mol Sci. 2021;22(12):6485.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Cerreto M, Santopaolo F, Gasbarrini A, Pompili M, Ponziani FR. Bariatric surgery and liver disease: general considerations and role of the gut-liver Axis. Nutrients. 2021;13(8):2649.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Suk KT, Kim DJ. Gut microbiota: novel therapeutic target for nonalcoholic fatty liver disease. Expert Rev Gastroenterol Hepatol. 2019;13(3):193–204.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors thank the nurses and patients as well as healthy participants involved in this study. We also thank the staff in Division of Gastroenterology in The First Affiliated Hospital of Kunming Medical University, Key Laboratory of Molecular Diagnostics and Precision Medicine for Surgical Oncology in Gansu Province of Gansu Provincial Hospital, and Key Laboratory of Radiation Technology and Biophysics in Institute of Biology & Hefei Institute of Physical Science in Chinese Academy of Sciences for their technical support.

Code availability

Not applicable.


The work was supported by Yunnan Science and Technology Project (L-2017018, 2018NS0100), Natural Science Foundation of Jiangxi Province (20212BAB216073), Science and technology projects of Guizhou Province (QKH-J-ZK[2021]-107), The Famous Medical Specialist of “High-level Talent Training Support Program” in Yunnan Province (2020), Natural Science Foundation of Tianjin (19JCQNJC12500), the Non-profit Central Research Institute Fund of Chinese Academy of Medical Sciences (2019PT320005), Yunnan Geriatrics Clinical Research Center (202102AA310069), Jiangxi Provincial Key New Product Incubation Program Funded by Technical Innovation Guidance Program of Shangrao (2020G002), the Reserve Training Project of “Thousand” Project of Health Science and Technology Talents in Kunming (2019-sw-52), Key project funded by Department of Science and Technology of Shangrao City (2020AB002, 2020 K003, 2021F013).

Author information

Authors and Affiliations



Q.Y., L.Z., and Q.L.: collection and assembly of data, manuscript writing; M.G., Q.Q., X.Y., Q.Y., K.G., L.K., M.H., and J.X.: collection and assembly of data; L.Z., and H.Y.: conception and design, data analysis and interpretation, final approval of the manuscript. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Leisheng Zhang or Hongju Yang.

Ethics declarations

Ethics approval and consent to participate

The study was performed according to Declaration of Helsinki, and was approved by Ethics Committee of First Affiliated Hospital of Kunming Medical University (approval number: 2020-L-08). Informed consent was obtained from the healthy donors included in the study.

Consent for publication

Not applicable.

Competing interests

All authors declare there’s no competing interests and consent to publish the data.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Additional Information.

The detailed information accompanied with the main manuscript, including Additional Information (including Additional Figure Legend for Fig. S1).

Additional file 2: Additional Figure S1.

The clinicopathological indicators involved in M and MD patients and Ctr.

Additional file 3: Additional Table S1.

Characteristics of patients with MAFLD, MAFLD/T2DM.

Additional file 4: Additional Table S2.

Clinical Parameters in M and MD patients and Ctr healthy participants.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, Q., Zhang, L., Li, Q. et al. Characterization of microbiome and metabolite analyses in patients with metabolic associated fatty liver disease and type II diabetes mellitus. BMC Microbiol 22, 105 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: