Basic Study Open Access
Copyright ©The Author(s) 2018. Published by Baishideng Publishing Group Inc. All rights reserved.
World J Gastroenterol. Jun 21, 2018; 24(23): 2468-2481
Published online Jun 21, 2018. doi: 10.3748/wjg.v24.i23.2468
Dynamic alterations in the gut microbiota and metabolome during the development of methionine-choline-deficient diet-induced nonalcoholic steatohepatitis
Jian-Zhong Ye, Ya-Ting Li, Wen-Rui Wu, Ding Shi, Dai-Qiong Fang, Li-Ya Yang, Xiao-Yuan Bian, Jing-Jing Wu, Qing Wang, Xian-Wan Jiang, Cong-Gao Peng, Lan-Juan Li, State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, School of Medicine, Zhejiang University, Hangzhou 310003, Zhejiang Province, China
Jian-Zhong Ye, Ya-Ting Li, Wen-Rui Wu, Ding Shi, Dai-Qiong Fang, Li-Ya Yang, Xiao-Yuan Bian, Jing-Jing Wu, Qing Wang, Xian-Wan Jiang, Cong-Gao Peng, Lan-Juan Li, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou 310003, Zhejiang Province, China
Wan-Chun Ye, Department of Chemotherapy 2, Wenzhou Central Hospital, Wenzhou 325000, Zhejiang Province, China
Peng-Cheng Xia, Department of Clinical Laboratory, Tai’an Central Hospital, Tai’an 271000, Shandong Province, China
ORCID number: Jian-Zhong Ye (0000-0001-8174-2580); Ya-Ting Li (0000-0002-0761-4967); Wen-Rui Wu (0000-0002-8457-9675); Ding Shi (0000-0003-3988-9321); Dai-Qiong Fang (0000-0002-0758-3988); Li-Ya Yang (0000-0002-7832-983X); Xiao-Yuan Bian (0000-0002-7275-3103); Jing-Jing Wu (0000-0001-9289-0538); Qing Wang (0000-0003-4914-5473); Xian-Wan Jiang (0000-0003-0982-7306); Cong-Gao Peng (0000-0002-0945-0986); Wan-Chun Ye (0000-0002-2732-5150); Peng-Cheng Xia (0000-0001-8990-4900); Lan-Juan Li (0000-0001-6945-0593).
Author contributions: Ye JZ, Li YT and Wu WR have contributed equally to this work. Ye JZ, Li YT, Wu WR and Li LJ conceived and designed the experiments; Ye JZ, Li YT and Wu WR performed the experiments; Ye JZ, Shi D, Fang DQ, Yang LY and Bian XY analyzed the data; Wu JJ, Wang Q and Jiang XW contributed analysis tools; Ye JZ, Li YT and Wu WR wrote the manuscript; Ye JZ and Peng CG conducted a critical revision of the manuscript for important intellectual content; Ye WC, Xia PC and Li LJ reviewed the manuscript; all authors have approved the final version of the manuscript.
Supported by the National Natural Science Foundation of China, No. 81330011, No. 81790631, and No. 81790633; the Science Fund for Creative Research Groups of the National Natural Science Foundation of China, No. 81721091; and the National Basic Research Program of China (973 program), No. 2013CB531401.
Institutional animal care and use committee statement: All procedures involving animals were reviewed and approved by the Institutional Animal Care and Use Committee of Zhejiang University School of Medicine (Permit number: 2017-591).
Conflict-of-interest statement: There is no conflict of interest in this study.
Data sharing statement: No additional data are available.
ARRIVE guidelines statement: The ARRIVE Guidelines have been adopted.
Open-Access: This article is an open-access article which was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Correspondence to: Lan-Juan Li, MD, PhD, Academic Research, Doctor, Professor, State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, School of Medicine, Zhejiang University, No.79 Qingchun Road, Hangzhou 310003, Zhejiang Province, China. ljli@zju.edu.cn
Telephone: +86-571-87236759 Fax: +86-571-87236459
Received: February 28, 2018
Peer-review started: March 1, 2018
First decision: March 15, 2018
Revised: March 27, 2018
Accepted: May 11, 2018
Article in press: May 11, 2018
Published online: June 21, 2018
Processing time: 106 Days and 20.3 Hours

Abstract
AIM

To investigate changes in gut microbiota and metabolism during nonalcoholic steatohepatitis (NASH) development in mice fed a methionine-choline-deficient (MCD) diet.

METHODS

Twenty-four male C57BL/6J mice were equally divided into four groups and fed a methionine-choline-sufficient diet for 2 wk (Control 2w group, n = 6) or 4 wk (Control 4w group, n = 6) or the MCD diet for 2 wk (MCD 2w group, n = 6) or 4 wk (MCD 4w group, n = 6). Liver injury, fibrosis, and intestinal barrier function were evaluated after 2 and 4 wk of feeding. The fecal microbiome and metabolome were studied using 16s rRNA deep sequencing and gas chromatography-mass spectrometry.

RESULTS

The mice fed the MCD diet presented with simple hepatic steatosis and slight intestinal barrier deterioration after 2 wk. After 4 wk of feeding with the MCD diet, however, the mice developed prominent NASH with liver fibrosis, and the intestinal barrier was more impaired. Compared with the control diet, the MCD diet induced gradual gut microbiota dysbiosis, as evidenced by a marked decrease in the abundance of Alistipes and the (Eubacterium) coprostanoligenes group (P < 0.001 and P < 0.05, respectively) and a significant increase in Ruminococcaceae UCG 014 abundance (P < 0.05) after 2 wk. At 4 wk, the MCD diet significantly reduced the promising probiotic Bifidobacterium levels and markedly promoted Bacteroides abundance (P < 0.05, and P < 0.01, respectively). The fecal metabolomic profile was also substantially altered by the MCD diet: At 2 wk, arachidic acid, hexadecane, palmitic acid, and tetracosane were selected as potential biomarkers that were significantly different in the corresponding control group, and at 4 wk, cholic acid, cholesterol, arachidic acid, tetracosane, and stearic acid were selected.

CONCLUSION

The MCD diet induced persistent alterations in the gut microbiota and metabolome.

Key Words: Nonalcoholic steatohepatitis; Methionine-choline deficient diet; Gut microbiota; Metabolome; Nonalcoholic fatty liver disease

Core tip: Nonalcoholic steatohepatitis (NASH) is increasingly prevalent as a remarkable problem worldwide. Increased evidence indicates the critical role of gut microbiota in NASH progression. We aimed to investigate the dynamic alterations in the gut microbiota and the related metabolites during NASH development in mice fed a methionine-choline-deficient (MCD) diet. We for the first time find that the MCD diet may induce persistent gut microbiota and metabolome deterioration.



INTRODUCTION

The prevalence of nonalcoholic fatty liver disease (NAFLD), the most common chronic liver disease worldwide, is increasing, encompassing a pathological spectrum from simple steatosis to nonalcoholic steatohepatitis (NASH)[1]. NASH is a liver characteristic of metabolic syndrome, accounting for lipid accumulation and hepatic inflammation[2] and is projected to be the main indication for liver transplantation in the next decade[3]. NASH is commonly associated with type 2 diabetes, cardiovascular disease, and end-stage kidney disease[4,5], and no approved pharmacological therapies currently exist for NASH[6].

Recently, compelling evidence has indicated the critical role of the gut microbiota in the pathogenesis and progression of NAFLD to NASH[7]. Normally, the commensal gut microbiota and the host maintain a beneficial symbiotic relationship. The liver may function as a vascular firewall, mediating mutualism between the host and the intestinal commensal bacteria[8]. This intimate connection between the gastrointestinal tract and liver defines the term gut-liver axis. However, gut microbiota dysbiosis increases the influx of harmful substances, such as lipopolysaccharide (LPS), ethanol, and bacterial DNA, into the liver through portal vein circulation and promotes NASH development[9,10]. Currently, early and noninvasive diagnosis of NASH with high sensitivity and specificity remains challenging. Metabolomics, along with other “omics” technologies, helps provide a detailed understanding of biochemical events using the systems biology approach and facilitates early diagnosis and exploitation of treatment strategies[11].

Despite the importance of the gut microbiota and metabolome in NASH, detailed information regarding their role during NASH development is limited. The purpose of the present study was to investigate dynamic alterations in the gut microbiota and metabolome during the development of methionine-choline-deficient (MCD) diet-induced NASH in a mouse model.

MATERIALS AND METHODS
Animals and treatments

The animal protocol was designed to minimize the animals’ pain or discomfort and was approved by the Animal Care Committee of Zhejiang University School of Medicine (Permit number: 2017-591). Male C57BL/6J mice (8 wk old) were purchased from SLAC Laboratory (Shanghai, China). After their arrival, the mice were acclimatized for 1 wk in a specific pathogen-free environment (23 °C, 12 h/12 h light/dark, 50% humidity, and ad libitum access to food and water) prior to experimentation. The mice were fed a methionine-choline-sufficient diet (Research Diet, New Brunswick, United States) for 2 wk (the Control 2w group, n = 6) or 4 wk (the Control 4w group, n = 6). Alternatively, the mice were fed an MCD diet (Research Diet) for 2 wk (the MCD 2w group, n = 6) or 4 wk (the MCD 4w group, n = 6). The diet specifications are summarized in Supplementary Table 1. After 2 or 4 wk on the diet, the mice were euthanized by an intraperitoneal injection of 4% chloral hydrate (with 1 mg/100 mL of atropine, to inhibit respiratory secretions) for tissue collection.

Sample collection

Fecal samples were collected from all mice upon defecation and were stored at -80 °C for further procedures. Blood samples were centrifuged (3000 rpm, 15 min) for serum separation, and all serum aliquots were stored in a -80 °C freezer. Liver and colon samples were fixed in either neutral-buffered formalin or optimal cutting temperature (OCT)-embedding media for histological staining or were snap-frozen in liquid nitrogen and kept in a -80 °C freezer for further procedures.

Histopathological evaluation of liver tissue

Paraffin-embedded liver sections were stained with hematoxylin-eosin (HE) or Masson’s trichrome to detect liver injury or fibrosis. The stained sections were scanned using a NanoZoomer Digital Pathology system (Hamamatsu Photonics, KK, Japan), which digitally scanned the sections into a particular image format for further assessment. The HE-stained sections were scored in accordance with the NAFLD activity score (NAS) system (score 0-2: Not NASH, 3-4: Borderline NASH, and 5-8: NASH)[12]. Six fields from each section were selected and analyzed at 200 × magnification. The Masson’s trichrome-stained sections were analyzed to quantify fibrosis using Image-Pro Plus software (version 6.0, Media Cybernetics, Rockville, United States) as previously detailed[13]. For each section, the blue area (collagen) was normalized to the red area (hepatocyte). The fibrosis index (%) was calculated as a percent of the total tissue region and represented the average of six randomly selected fields from each section.

Oil red O staining and liver triglyceride assay

Frozen liver sections fixed in OCT (10 μm) were stained with Oil red O (Sigma-Aldrich, St. Louis, MO, United States). Images were captured using the abovementioned NanoZoomer Digital Pathology system. To quantify intrahepatic lipid accumulation, the mean optical density of the red intensity was assessed using Image-Pro Plus software.

The liver triglyceride (TG) content was determined using a commercial kit from Applygen Technologies Inc. (Beijing, China) according to the manufacturer’s protocols, and the final TG concentrations were normalized to the corresponding protein content.

Immunohistochemistry and immunofluorescence staining

Paraffin-embedded liver sections were stained for F4/80 (anti-active macrophage) (Abcam, Cambridge, United Kingdom) and α-SMA (fibrosis hallmark) with immunohistochemistry (IHC) staining procedures as previously detailed[14]. Briefly, liver sections were incubated with a specific primary antibody, followed by incubation with horseradish peroxidase (HRP)-linked secondary antibody (Dako, Glostrup, Denmark) and 3,3’-diaminobenzidine; the sections were then scanned with the NanoZoomer Digital Pathology system. Image-Pro Plus software was used to count F4/80+ cells and quantitatively analyze the staining intensity of α-SMA as previously described[13]. Six fields of view were randomly selected in each section.

Likewise, paraffin-embedded colon sections were stained for Zonula occludens-1 (ZO-1) (intestinal barrier hallmark) (Proteintech, Rosemont, IL, United States) with standard immunofluorescence staining procedures as previously detailed[15]. Briefly, sections were incubated with the rabbit polyclonal ZO-1 antibody, followed by incubation with Texas Red-conjugated goat anti-rabbit antibody (Jackson ImmunoResearch, West Grove, PA, United States) and 4’,6-diamino-2-phenyl indole (DAPI), and images were captured using a Zeiss LSM T-PMT confocal microscope (Zeiss, Jena, Germany).

Microbial community analysis of fecal samples

Total DNA was extracted from fecal samples with a QIAamp fast DNA stool mini kit (Qiagen, Valencia, CA, United States) following the manufacturer’s handbook. After determining the DNA concentration and integrity, an amplicon sequencing library was constructed based on the PCR-amplified V3-V4 variable regions of 16s rDNA. Then, the qualified libraries were paired-end sequenced on an Illumina MiSeq platform according to manufacturer’s procedures. Raw sequencing data were subjected to filtration using Trimmomatic, FLASH, and QIIME software. Then, clean reads were clustered into operational taxonomic units (OTUs) using UPARSE software with a 97% threshold. The representative read from each OTU was selected using the QIIME package[16]. Representative OTU sequences were annotated and taxonomically classified using Ribosomal Database Project (RDP) Classifier v.2.2, trained on the Silva database version 123[17]. The linear discriminant analysis (LDA) effect size (LEfSe) method (http://huttenhower.sph.harvard.edu/galaxy/) was applied to differentiate taxa with statistical significance and biological relevance[18].

Metabolomic profiling analysis of fecal samples

Metabolomic profiling analysis was performed as previously described[19] with slight modification. Briefly, fecal metabolites were extracted by mixing 15 mg of feces with 800 μL of ice-cold methanol. After homogenization and centrifugation, the supernatant was transferred into an Eppendorf tube containing 20 μL of 1 mg/mL heptadecanoic acid as the internal standard. Then, the sample was dried using a nitrogen stream (Aosheng, Hangzhou, China). The residue was reconstituted in 50 μL of 15 mg/mL methoxylamine hydrochloride in anhydrous pyridine and was incubated at 37 °C for 24 h. Then, 50 μL of N,O-bistrifluoroacetamide (BSTFA) [with 1% trimethylsilyl chloride (TMCS)] (Sigma-Aldrich, St. Louis, MO, United States) was added to the mixture, and the sample was incubated at 70 °C for 120 min. Metabolomic analysis was performed with a gas chromatography-mass spectrometry (GC-MS) on an Agilent 7890A GC system coupled to an Agilent 5975C inert mass selective detector (MSD) system (Agilent Technologies, Santa Clara, CA, United States). For data analysis, ChemStation software (version E.02.02.1431), and ChromaTOF software (version 4.34, LECO, St. Joseph, MI, United States) were used. Metabolites were identified by the National Institute of Standards and Technology (NIST) and Fiehn databases. Principal component analysis (PCA) and orthogonal partial least-squares-discriminant analysis (OPLS-DA) were performed to visualize metabolic differences among the experimental groups. Differential metabolites were selected according to the statistically significant VIP values obtained from the OPLS-DA model, and the P values from two-tailed Student’s t-tests on the normalized peak areas; metabolites with VIP values > 1 and P values < 0.05 were included.

Statistical analysis

The data are shown as the means ± SEM, and the Kolmogorov-Smirnov test was performed to assess data normality. For most data, one-way ANOVA with Tukey’s post hoc test was used to determine the significance between the groups. The Wilcoxon rank sum test was used to evaluate alpha diversity and principal coordinates between the different cohorts in the 16s sequencing analysis. Analysis of similarities (ANOSIM) was performed to test for microbial community clustering according to unweighted UniFrac distance matrices. P values < 0.05 were considered significant. The data were analyzed using SPSS version 17.0 for Windows (SPSS Inc., Chicago, IL, United States).

RESULTS
MCD diet resulted in gradual liver injury and intestinal barrier impairment

After 2 wk, the mice fed the MCD diet developed simple hepatic steatosis: the liver showed predominantly microvesicular fat accumulation without lobular inflammation and ballooning (Figure 1A, Table 1). Hepatic steatosis was significantly increased in the MCD 2w group compared with that in the Control 2w group, as further indicated by Oil red O staining and hepatic TG quantification (Figure 1B and D, Table 1). No significant difference was found in the number of F4/80+ cells between the MCD 2w group and the Control 2w group (Figure 1C and E). Compared with the Control 2w group, the MCD 2w group presented significantly increased serum alanine aminotransferase (ALT) (P < 0.01) and aspartate aminotransferase (AST) (P < 0.01) levels (Table 1). Evidence of liver fibrosis was not found, as shown by Masson’s trichrome and α-SMA immunohistochemical staining (Figure 2). Intestinal barrier destruction is associated with NAFLD progression[10]; therefore, we detected tight junction protein expression in the colon, where the most abundant gut microbiota reside[20]. ZO-1 immunostaining revealed that the colon tissues of the mice in the MCD 2w group exhibited increased disruption and disorganization on the apical surface and in the crypts (Figure 3).

Table 1 Effects of the methionine-choline-deficient diet on nonalcoholic fatty liver disease activity score, alanine aminotransferase, aspartate aminotransferase, and hepatic triglycerides levels.
Control 2wMCD 2wControl 4wMCD 4w
NAS0.17 ± 0.171.17 ± 0.17a0.33 ± 0.214.83 ± 0.31cfj
Steatosis0.00 ± 0.001.17 ± 0.17b0.33 ± 0.212.50 ± 0.22cej
Inflammation0.00 ± 0.000.00 ± 0.000.00 ± 0.001.17 ± 0.17bei
Ballooning0.17 ± 0.170.00 ± 0.000.00 ± 0.001.17 ± 0.31
TG (nmol/mg protein)179.41 ± 21.72592.88 ± 55.84b204.45 ± 6.76e731.60 ± 113.46ah
ALT (U/L)19.33 ± 0.88413.83 ± 81.77a34.83 ± 2.64bd974.17 ± 163.00ah
AST (U/L)105.50 ± 10.55327.33 ± 39.92a109.83 ± 11.91d749.677 ± 92.34bdi
Figure 1
Figure 1 Methionine-choline-deficient diet induced hepatic injury, fat accumulation, and macrophage infiltration. A-C: Representative liver histology assessed by HE staining (A), Oil red O staining (B), and macrophages (F4/80) staining (C). Scale bar: 250 μm. D: Representative staining intensities of Oil red O as designated by mean optical density. E: Percentage of F4/80 positive cells. The data are given as the means ± SEM. n = 6 per group. cP < 0.001 vs Control 2w group, eP < 0.01 and fP < 0.001 vs MCD 2w group, jP < 0.001 vs Control 4w group by post hoc ANOVA one-way statistical analysis. MCD: Methionine-choline-deficient.
Figure 2
Figure 2 Methionine-choline-deficient diet induced liver fibrosis after 4 wk of feeding. A-B: Representative liver histology assessed by Masson’s trichrome staining (A) and α-SMA staining (B). Scale bar: 250 μm. C: Percent of Masson’s trichrome-stained area as indicated by the fibrosis index (%). D: Representative staining intensities of α-SMA as designated by the mean optical density. The data are given as the means ± SEM. n = 6 per group. aP < 0.05 and bP < 0.01 vs Control 2w group, dP < 0.05 and eP < 0.01 vs MCD 2w group, hP < 0.05 and iP < 0.01 vs Control 4w group by post hoc ANOVA one-way statistical analysis. MCD: Methionine-choline-deficient.
Figure 3
Figure 3 Methionine-choline-deficient diet resulted in gradual intestinal barrier impairment. Representative colon histology was assessed by ZO-1 immunofluorescence staining. Scale bar: 100 μm. The data are given as the means ± SEM. MCD: Methionine-choline-deficient.

As expected, the mice in the MCD 4w group developed prominent NASH, as evidenced by major hepatic steatosis with lobular inflammation and ballooning hepatocytes in the liver (Figure 1A, Table 1). Oil red O staining and hepatic TG quantification also revealed that fat accumulation was significantly increased in the livers of the mice fed the MCD diet for 4 wk (the MCD 4w group) compared with that in the mice fed the control diet for 4 wk (the Control 4w group) (Figure 1B and D, Table 1). F4/80+ cell infiltration was significantly increased in the MCD 4w group compared with that in the Control 4w group (Figure 1C and E). In addition, we found that the mice in the MCD 4w group exhibited indications of liver fibrosis, including periportal and interstitial collagen deposition (Figure 2A and C). Immunohistochemical analysis further confirmed our result: α-SMA protein expression was significantly up-regulated in the MCD 4w group compared with that in the Control 4w group (Figure 2B and D). Intestinal barrier function was further impaired in the mice fed the MCD diet for 4 wk compared with that in the mice fed the control diet, as revealed by ZO-1 immunostaining (Figure 3).

MCD diet induced gradual gut microbiota dysbiosis

The MCD diet clearly altered the gut microbiota configuration. No significant differences were observed in alpha diversity between the control and MCD groups after 2 and 4 wk of treatment, as estimated by the Chao1, Shannon, and Simpson indices (data not shown). However, these two groups were clearly separated into different clusters at 2 wk (ANOSIM, P < 0.01, r = 0.6352) and at 4 wk (ANOSIM, P < 0.01, r = 0.9074) in the unweighted UniFrac principal coordinate analysis (PCoA), which was performed to calculate the beta-diversity values (Figure 4A). The most abundant taxa at the phylum, family, and genus levels are shown in Figure 4B-D. At the phylum level, Tenericutes was increasingly more abundant in the fecal microbiota of the MCD group compared with that of the control group at 2 and 4 wk (P < 0.05, and P < 0.01, respectively), while Verrucomicrobia was consistently less abundant in the MCD group at 2 and 4 wk (P < 0.05, and P < 0.05, respectively). Compared with the control group, the MCD group had a significantly higher abundance of Firmicutes and a significantly reduced abundance of Proteobacteria at 2 wk (P < 0.001, and P < 0.05, respectively). At 4 wk, the abundance of Actinobacteria was significantly lower in the MCD group than that in the corresponding control group (P < 0.01). At the family level, the relative abundance levels of Rikenellaceae, Desulfovibrionaceae, and Verrucomicrobiaceae were persistently reduced in the MCD group compared with those in the control group at 2 wk (P < 0.001, P < 0.05, and P < 0.05 respectively) and at 4 wk (P < 0.05, P < 0.05, and P < 0.05 respectively). The relative abundance of the Bacteroidales S24-7 group was significantly lower in the MCD group than that in the control group at 2 wk (P < 0.01), while Ruminococcaceae was significantly higher in the MCD group (P < 0.05). Compared with the control group, the MCD group had significantly higher abundance levels of Bacteroidaceae and Enterobacteriaceae at 4 wk (P < 0.01, and P < 0.05 respectively) and a significantly reduced abundance of Bifidobacteriaceae (P < 0.05). At the genus level, the abundance of the Rikenellaceae RC9 gut group was significantly reduced in the MCD group compared with that in the control group at 2 and 4 wk (P < 0.001, and P < 0.05, respectively). Compared with the control group, the MCD group presented with marked decreases in the abundance levels of Alistipes and (Eubacterium) coprostanoligenes at 2 wk (P < 0.001, and P < 0.05, respectively) and a significant increase in Ruminococcaceae UCG 014 abundance (P < 0.05). However, at 4 wk, the MCD diet significantly reduced the abundance of the promising probiotic Bifidobacterium and markedly promoted Bacteroides abundance (P < 0.05, and P < 0.01, respectively).

Figure 4
Figure 4 Methionine-choline-deficient diet induced gut microbiota dysbiosis. A: PCoA plot of the microbiota based on unweighted UniFrac metric. Each symbol represents one sample (n = 6 per group). B-D: Top most abundant taxa at the phylum (B), family (C), and genus (D) levels. E: LEfSe cladogram represented taxon enriched in the Control 4w group (green) and in the MCD 4w group (red). Rings from the inside out represented taxonomic levels from phylum to genus. Sizes of the circles indicate the relative abundance levels of the taxa. F: Discriminative biomarkers with an LDA score > 4.8. MCD: Methionine-choline-deficient.

To characterize further the distinguishing phylotypes in the gut microbiota of the two groups, LEfSe analysis was performed. No significant differences were found between the MCD 2w group and Control 2w group. At 4 wk, however, we found that the MCD diet increased the abundance levels of Anaerotruncus, Bilophila, SMB53, Clostridium, Anaeroplasma, Turicibacter, Helicobacteraceae, Flexispira, and Bacteroides [LDA score (-log10) > 4.8] and decreased the abundance levels of Allobaculum, S24-7, Bifidobacterium, Adlercreutzia, Lachnospiraceae, Akkermansia, Sutterella, Desulfovibrionaceae, Porphyromonadaceae, Parabacteroides, and Erysipelotrichaceae [LDA score (log10) > 4.8] (Figure 4E and F).

MCD diet altered the fecal metabolomic profile during NASH progression

Using an untargeted strategy, we studied the fecal metabolome associated with functional characteristics of the gut microbiome, and a total of 322 metabolites were ultimately identified and quantified.

The PCA model was established (R2X = 0.526, Q2 = 0.223), corresponding to the four groups (the Control 2w group, MCD 2w group, Control 4w group, and MCD 4w group) (Figure 5A), and the score plot showed a clustering tendency in the first predictive principal component (X axis) and second predictive principal component (Y axis). Then, an OPLS-DA model was constructed. As depicted in Figure 5B, the score plot in the direction of the X axis (the first predictive principal component) and Y axis (the first orthogonal principal component) showed a significant separation in the metabolomics datasets among the four groups. The explained variance, R2, was 0.968. The cross-validated predictive ability Q2 was 0.914, indicating that a random fecal GC-MS spectrum discriminates among the four groups 91.4% of the time. Therefore, these results indicate the distinct clustering of the fecal metabolomic profiles induced by the MCD diet over time. Characteristic metabolites that had been significantly modified by the MCD diet were further identified according to the OPLS-DA model, with VIP values > 1 and P values < 0.05 (Figure 5E and F, and Supplementary Tables 2 and 3). Ultimately, 103 and 93 metabolites were selected at 2 and 4 wk, respectively, most of which were mainly associated with pathways involved in lipid, amino acid, carbohydrate, nucleotide, cofactors, and vitamin metabolism. S-plots were used to identify further potential biomarkers among the metabolites. As shown in Figure 5C and D, arachidic acid, hexadecane, palmitic acid, and tetracosane at 2 wk and cholic acid, cholesterol, arachidic acid, tetracosane, and stearic acid at 4 wk were the farthest from the origin and were selected as potential biomarkers due to their marked contribution to the separation between the control and MCD groups.

Figure 5
Figure 5 Methionine-choline-deficient diet altered the fecal metabolomic profile during nonalcoholic steatohepatitis progression. A: 4-state model of a PCA plot comparing between the Control 2w (cycle), MCD 2w (triangle), Control 4w (square), and MCD 4w groups (diamond). B: OPLS-DA score scatter plot comparing the Control 2w (cycle), MCD 2w (triangle), Control 4w (square), and MCD 4w groups (diamond). Each symbol represents one sample (n = 24 per group, with a quadruplicate technical replicates of 6 samples per group). C-D: S-plot constructed from the OPLS-DA model at 2 wk (C) and 4 wk (D). The selected metabolites (red triangles) are located in the upper far right and lower far left with considerable potential as biomarkers and marked discrimination between the groups. E-F: Euclidean distance hierarchical clustering analysis visualizing the different intensity levels of characteristic metabolites at 2 wk (E) and 4 wk (F). PCA: Principal component analysis; MCD: Methionine-choline-deficient; OPLS-DA: Orthogonal partial least-squares-discriminant analysis.
DISCUSSION

To the best of our knowledge, this is the first study to demonstrate dynamic alterations in the gut microbiota and metabolome in an experimental model of MCD diet-induced steatohepatitis: We sought to determine the key microbiota and metabolites involved in NASH progression over time.

Mice fed the MCD diet developed simple hepatic steatosis at 2 wk and prominent NASH at 4 wk. Among the most abundant taxa, the prevalence of the Rikenellaceae RC9 gut group in the family Rikenellaceae, family Desulfovibrionaceae, and family Verrucomicrobiaceae in the phylum Verrucomicrobia was consistently decreased at 2 and 4 wk in the MCD-supplemented group compared with that in the control group, while the phylum Tenericutes was consistently increased in the MCD-supplemented group. The phylum Tenericutes was also increased in the mice fed a high-fat diet and is correlated with obesity-associated metabolic characteristics[21]. The Rikenellaceae RC9 gut group was previously shown to be decreased in rats with hypertriglyceridemia-related acute necrotizing pancreatitis and in cows receiving oil supplementation[22,23], but the specific function of this taxa remains poorly understood. High Rikenellaceae abundance is reportedly associated with healthy metabolic states[24], and one study reported the ability of some of the bacteria in the family Rikenellaceae to produce butyrate[25]. Desulfovibrionaceae are known as sulfate-reducing bacteria that produce hydrogen sulfide (H2S)[26], which has its pros and cons in relation to gut health. Some studies have demonstrated that H2S functions as a mucosal barrier breaker and therefore leads to inflammation[27,28]. Regarding the advantages of H2S, previous studies have reported that H2S is a crucial mediator in gastrointestinal mucosal defense and repair and in reducing systemic inflammation[26,29]. Verrucomicrobiaceae abundance is positively correlated with gastrointestinal health and negatively correlated with gut inflammation[30]. Therefore, disrupted barrier homeostasis, as evidenced by the ZO-1 immunostaining results in this study, may be associated with a decrease in MCD-sensitive microbiota.

In this study, we found that Bifidobacterium in the family Bifidobacteriaceae, Bacteroides in the family Bacteroidaceae, and the family Enterobacteriaceae respond to the MCD diet in a time-dependent manner. The first of these taxa tended to decrease, while the latter two tended to increase after only 4 wk of feeding. Consistently, Bifidobacterium was significantly decreased in NASH subjects compared with that in healthy and obese subjects[10]. Various members of the genus Bifidobacterium have recently attracted substantial interest due to their multiple beneficial effects on the host, such as protection against enteropathogenic infection through acetate production[31], promotion of antitumor immunity, and facilitation of programmed cell death protein 1 ligand 1 (PD-L1)-specific antibody therapy[32]. Although Bacteroides spp. are overall regarded as beneficial microorganisms[33], they can also cause gut-related bacteremia in conditions of high portal venous pressure[34]. Enterobacteriaceae include many potential pathogens and LPS-producing bacteria and are significantly increased in NASH patients[10]. Furthermore, as evidenced by LEfSe analysis, the MCD diet favored LPS-producing bacteria (Bilophila, and Anaeroplasma[35,36]), pro-inflammatory bacteria (Anaerotruncus, and Turicibacter[37,38]), and other opportunistic pathogens (SMB53, Clostridium, Helicobacteraceae and Flexispira[39-42]) after 4 wk of feeding and inhibited the promising probiotics Sutterella and Akkermansia[43,44]. Therefore, these data indicate an important role for the microbiota in NAFLD progression to NASH.

To link the gut microbiota with their functional states, an untargeted metabolomics analysis was integrated into this study. A core of the MCD diet-responsive signaling pathways were involved in lipid, nucleotide, cofactor, vitamin, carbohydrate, and amino acid metabolism. S-plots further identified potential biomarkers among the metabolites during NASH development. Arachidic acid under-representation and tetracosane over-representation were prominent features in the group with MCD diet-induced steatosis and NASH compared with the corresponding control groups at 2 and 4 wk. However, functional exploration of the association between arachidic acid and tetracosane with non-alcoholic fatty liver disease has not been well documented and warrants further studies. Kuroda et al[45] reported that hexadecane causes nonspecific inflammation. Therefore, we speculated in this study that the increased hexadecane levels after the MCD diet treatment for 2 wk may contribute to disease progression, although we have not yet obtained direct evidence for this effect. Jiao et al[46] found that the cholic acid level is increased in NASH patients. JianHua et al[47] found that cholic acid is significantly decreased in humans with NASH, which is consistent with our study results in the MCD 4w group. Stearic acid is a potent anti-inflammatory lipid and may accelerate hepatic dysfunction recovery in a rat model of liver injury[48]. Our study found that the stearic acid level was significantly decreased in the NASH mice after the MCD diet treatment for 4 wk. Although our study found an association between altered gut microbiota and metabolism and NASH, a causative contribution of the gut microbiota and metabolism to NASH progression has not been sufficiently documented. Establishing a better understanding of the gut microbiota and metabolomics in this disease state will provide beneficial information for the treatment and prevention of NAFLD.

In conclusion, the MCD diet induced gut microbiota and metabolome deterioration. Fundamental observations of these alterations will provide new insight into NASH-associated intestinal disorder and gut-targeted therapies for NASH.

ARTICLE HIGHLIGHTS
Research background

The contributing role of the gut microbiota in the pathogenesis of nonalcoholic steatohepatitis (NASH) has been extensively studied.

Research motivation

Gut microbiota dysbiosis in NASH is mainly depicted as an endpoint, and little is known regarding the microbiota disturbances during NASH progression.

Research objectives

Our goal was to investigate dynamic changes in the gut microbiota and its metabolism during the progression from simple hepatic steatosis to NASH in mice fed a methionine-choline-deficient (MCD) diet.

Research methods

C57BL/6J mice were equally divided into four groups and fed either a methionine-choline-sufficient diet for 2 or 4 wk (the Control 2w group and Control 4w group, respectively) or the MCD diet for 2 or 4 wk (the MCD 2w group and MCD 4w group, respectively) (n = 6 per group). Liver injury, fibrosis, intestinal barrier function, and the fecal microbiome and metabolome were studied.

Research results

The mice fed with the MCD diet for 2 wk developed simple hepatic steatosis, which progressed to prominent NASH with liver fibrosis after 4 wk. Compared with the control diet, the MCD diet induced gradual intestinal barrier impairment and gut microbiota dysbiosis; the fecal metabolomic profile was also substantially altered by the MCD diet.

Research conclusions

The MCD diet induced persistent alteration of the gut microbiota and metabolome.

Research perspectives

We may have for the first time shown that an MCD diet induced persistent gut microbiota and metabolome deterioration. Fundamental observations of these alterations will provide new insight into NASH-associated intestinal disorder and gut-targeted therapies for NASH.

Footnotes

Manuscript source: Unsolicited manuscript

Specialty type: Gastroenterology and hepatology

Country of origin: China

Peer-review report classification

Grade A (Excellent): 0

Grade B (Very good): B, B

Grade C (Good): 0

Grade D (Fair): D

Grade E (Poor): 0

P- Reviewer: Das UN, Shimizu Y, Trovato GMM S- Editor: Wang XJ L- Editor: Filipodia E- Editor: Huang Y

References
1.  Musso G, Cassader M, Gambino R. Non-alcoholic steatohepatitis: emerging molecular targets and therapeutic strategies. Nat Rev Drug Discov. 2016;15:249-274.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 278]  [Cited by in F6Publishing: 323]  [Article Influence: 40.4]  [Reference Citation Analysis (0)]
2.  Ishioka M, Miura K, Minami S, Shimura Y, Ohnishi H. Altered Gut Microbiota Composition and Immune Response in Experimental Steatohepatitis Mouse Models. Dig Dis Sci. 2017;62:396-406.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 34]  [Cited by in F6Publishing: 33]  [Article Influence: 4.7]  [Reference Citation Analysis (0)]
3.  Wong RJ, Aguilar M, Cheung R, Perumpail RB, Harrison SA, Younossi ZM, Ahmed A. Nonalcoholic steatohepatitis is the second leading etiology of liver disease among adults awaiting liver transplantation in the United States. Gastroenterology. 2015;148:547-555.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 1211]  [Cited by in F6Publishing: 1315]  [Article Influence: 146.1]  [Reference Citation Analysis (1)]
4.  Chalasani N, Younossi Z, Lavine JE, Diehl AM, Brunt EM, Cusi K, Charlton M, Sanyal AJ. The diagnosis and management of non-alcoholic fatty liver disease: practice Guideline by the American Association for the Study of Liver Diseases, American College of Gastroenterology, and the American Gastroenterological Association. Hepatology. 2012;55:2005-2023.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 2413]  [Cited by in F6Publishing: 2528]  [Article Influence: 210.7]  [Reference Citation Analysis (1)]
5.  Musso G, Gambino R, Tabibian JH, Ekstedt M, Kechagias S, Hamaguchi M, Hultcrantz R, Hagström H, Yoon SK, Charatcharoenwitthaya P. Association of non-alcoholic fatty liver disease with chronic kidney disease: a systematic review and meta-analysis. PLoS Med. 2014;11:e1001680.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 408]  [Cited by in F6Publishing: 470]  [Article Influence: 47.0]  [Reference Citation Analysis (0)]
6.  Armstrong MJ, Gaunt P, Aithal GP, Barton D, Hull D, Parker R, Hazlehurst JM, Guo K; LEAN trial team, Abouda G, Aldersley MA, Stocken D, Gough SC, Tomlinson JW, Brown RM, Hübscher SG, Newsome PN. Liraglutide safety and efficacy in patients with non-alcoholic steatohepatitis (LEAN): a multicentre, double-blind, randomised, placebo-controlled phase 2 study. Lancet. 2016;387:679-690.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 1100]  [Cited by in F6Publishing: 1304]  [Article Influence: 163.0]  [Reference Citation Analysis (1)]
7.  Kirpich IA, Marsano LS, McClain CJ. Gut-liver axis, nutrition, and non-alcoholic fatty liver disease. Clin Biochem. 2015;48:923-930.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 158]  [Cited by in F6Publishing: 213]  [Article Influence: 23.7]  [Reference Citation Analysis (0)]
8.  Balmer ML, Slack E, de Gottardi A, Lawson MA, Hapfelmeier S, Miele L, Grieco A, Van Vlierberghe H, Fahrner R, Patuto N. The liver may act as a firewall mediating mutualism between the host and its gut commensal microbiota. Sci Transl Med. 2014;6:237ra66.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 283]  [Cited by in F6Publishing: 324]  [Article Influence: 32.4]  [Reference Citation Analysis (0)]
9.  Henao-Mejia J, Elinav E, Jin C, Hao L, Mehal WZ, Strowig T, Thaiss CA, Kau AL, Eisenbarth SC, Jurczak MJ. Inflammasome-mediated dysbiosis regulates progression of NAFLD and obesity. Nature. 2012;482:179-185.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 1620]  [Cited by in F6Publishing: 1800]  [Article Influence: 150.0]  [Reference Citation Analysis (0)]
10.  Zhu L, Baker SS, Gill C, Liu W, Alkhouri R, Baker RD, Gill SR. Characterization of gut microbiomes in nonalcoholic steatohepatitis (NASH) patients: a connection between endogenous alcohol and NASH. Hepatology. 2013;57:601-609.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 1015]  [Cited by in F6Publishing: 1184]  [Article Influence: 107.6]  [Reference Citation Analysis (1)]
11.  Safaei A, Arefi Oskouie A, Mohebbi SR, Rezaei-Tavirani M, Mahboubi M, Peyvandi M, Okhovatian F, Zamanian-Azodi M. Metabolomic analysis of human cirrhosis, hepatocellular carcinoma, non-alcoholic fatty liver disease and non-alcoholic steatohepatitis diseases. Gastroenterol Hepatol Bed Bench. 2016;9:158-173.  [PubMed]  [DOI]  [Cited in This Article: ]
12.  Kleiner DE, Brunt EM, Van Natta M, Behling C, Contos MJ, Cummings OW, Ferrell LD, Liu YC, Torbenson MS, Unalp-Arida A. Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology. 2005;41:1313-1321.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 6807]  [Cited by in F6Publishing: 7788]  [Article Influence: 409.9]  [Reference Citation Analysis (5)]
13.  Julio Junior HR, Costa SF, Costa WS, Barcellos Sampaio FJ, Favorito LA. Structural study of the bladder in fetuses with prune belly syndrome. Neurourol Urodyn. 2018;37:148-152.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 8]  [Cited by in F6Publishing: 9]  [Article Influence: 1.3]  [Reference Citation Analysis (0)]
14.  Wu YS, Tseng JK, Chou CH, Chiu CH, Lin YL, Chen YC. Preventive effects of Ophiocordyceps sinensis mycelium on the liver fibrosis induced by thioacetamide. Environ Toxicol. 2017;32:1792-1800.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 9]  [Cited by in F6Publishing: 12]  [Article Influence: 1.7]  [Reference Citation Analysis (0)]
15.  Chung CY, Alden SL, Funderburg NT, Fu P, Levine AD. Progressive proximal-to-distal reduction in expression of the tight junction complex in colonic epithelium of virally-suppressed HIV+ individuals. PLoS Pathog. 2014;10:e1004198.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 46]  [Cited by in F6Publishing: 50]  [Article Influence: 5.0]  [Reference Citation Analysis (0)]
16.  Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335-336.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 23810]  [Cited by in F6Publishing: 23146]  [Article Influence: 1653.3]  [Reference Citation Analysis (0)]
17.  Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-Syed-Mohideen AS, McGarrell DM, Marsh T, Garrity GM. The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009;37:D141-D145.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 3794]  [Cited by in F6Publishing: 3441]  [Article Influence: 229.4]  [Reference Citation Analysis (0)]
18.  Lu H, Ren Z, Li A, Zhang H, Jiang J, Xu S, Luo Q, Zhou K, Sun X, Zheng S. Deep sequencing reveals microbiota dysbiosis of tongue coat in patients with liver carcinoma. Sci Rep. 2016;6:33142.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 62]  [Cited by in F6Publishing: 73]  [Article Influence: 9.1]  [Reference Citation Analysis (0)]
19.  Liu R, Hong J, Xu X, Feng Q, Zhang D, Gu Y, Shi J, Zhao S, Liu W, Wang X. Gut microbiome and serum metabolome alterations in obesity and after weight-loss intervention. Nat Med. 2017;23:859-868.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 692]  [Cited by in F6Publishing: 948]  [Article Influence: 135.4]  [Reference Citation Analysis (0)]
20.  Wu S, Yi J, Zhang YG, Zhou J, Sun J. Leaky intestine and impaired microbiome in an amyotrophic lateral sclerosis mouse model. Physiol Rep. 2015;3:pii: e12356.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 127]  [Cited by in F6Publishing: 174]  [Article Influence: 19.3]  [Reference Citation Analysis (0)]
21.  Robertson RC, Seira Oriach C, Murphy K, Moloney GM, Cryan JF, Dinan TG, Ross RP, Stanton C. Deficiency of essential dietary n-3 PUFA disrupts the caecal microbiome and metabolome in mice. Br J Nutr. 2017;118:959-970.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 32]  [Cited by in F6Publishing: 35]  [Article Influence: 5.0]  [Reference Citation Analysis (0)]
22.  Huang C, Chen J, Wang J, Zhou H, Lu Y, Lou L, Zheng J, Tian L, Wang X, Cao Z. Dysbiosis of Intestinal Microbiota and Decreased Antimicrobial Peptide Level in Paneth Cells during Hypertriglyceridemia-Related Acute Necrotizing Pancreatitis in Rats. Front Microbiol. 2017;8:776.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 60]  [Cited by in F6Publishing: 78]  [Article Influence: 11.1]  [Reference Citation Analysis (0)]
23.  Zened A, Combes S, Cauquil L, Mariette J, Klopp C, Bouchez O, Troegeler-Meynadier A, Enjalbert F. Microbial ecology of the rumen evaluated by 454 GS FLX pyrosequencing is affected by starch and oil supplementation of diets. FEMS Microbiol Ecol. 2013;83:504-514.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 152]  [Cited by in F6Publishing: 168]  [Article Influence: 14.0]  [Reference Citation Analysis (0)]
24.  Lippert K, Kedenko L, Antonielli L, Kedenko I, Gemeier C, Leitner M, Kautzky-Willer A, Paulweber B, Hackl E. Gut microbiota dysbiosis associated with glucose metabolism disorders and the metabolic syndrome in older adults. Benef Microbes. 2017;8:545-556.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 157]  [Cited by in F6Publishing: 221]  [Article Influence: 31.6]  [Reference Citation Analysis (0)]
25.  Su XL, Tian Q, Zhang J, Yuan XZ, Shi XS, Guo RB, Qiu YL. Acetobacteroides hydrogenigenes gen. nov., sp. nov., an anaerobic hydrogen-producing bacterium in the family Rikenellaceae isolated from a reed swamp. Int J Syst Evol Microbiol. 2014;64:2986-2991.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 124]  [Cited by in F6Publishing: 139]  [Article Influence: 13.9]  [Reference Citation Analysis (0)]
26.  Konturek PC, Koziel J, Dieterich W, Haziri D, Wirtz S, Glowczyk I, Konturek K, Neurath MF, Zopf Y. Successful therapy of Clostridium difficile infection with fecal microbiota transplantation. J Physiol Pharmacol. 2016;67:859-866.  [PubMed]  [DOI]  [Cited in This Article: ]
27.  Sheng L, Jena PK, Liu HX, Kalanetra KM, Gonzalez FJ, French SW, Krishnan VV, Mills DA, Wan YY. Gender Differences in Bile Acids and Microbiota in Relationship with Gender Dissimilarity in Steatosis Induced by Diet and FXR Inactivation. Sci Rep. 2017;7:1748.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 77]  [Cited by in F6Publishing: 81]  [Article Influence: 11.6]  [Reference Citation Analysis (0)]
28.  Figliuolo VR, Dos Santos LM, Abalo A, Nanini H, Santos A, Brittes NM, Bernardazzi C, de Souza HSP, Vieira LQ, Coutinho-Silva R. Sulfate-reducing bacteria stimulate gut immune responses and contribute to inflammation in experimental colitis. Life Sci. 2017;189:29-38.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 49]  [Cited by in F6Publishing: 83]  [Article Influence: 11.9]  [Reference Citation Analysis (0)]
29.  Wallace JL, Motta JP, Buret AG. Hydrogen sulfide: an agent of stability at the microbiome-mucosa interface. Am J Physiol Gastrointest Liver Physiol. 2018;314:G143-G149.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 85]  [Cited by in F6Publishing: 71]  [Article Influence: 11.8]  [Reference Citation Analysis (0)]
30.  Rodriguez C, Taminiau B, Korsak N, Avesani V, Van Broeck J, Brach P, Delmée M, Daube G. Longitudinal survey of Clostridium difficile presence and gut microbiota composition in a Belgian nursing home. BMC Microbiol. 2016;16:229.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 33]  [Cited by in F6Publishing: 35]  [Article Influence: 4.4]  [Reference Citation Analysis (0)]
31.  Fukuda S, Toh H, Hase K, Oshima K, Nakanishi Y, Yoshimura K, Tobe T, Clarke JM, Topping DL, Suzuki T. Bifidobacteria can protect from enteropathogenic infection through production of acetate. Nature. 2011;469:543-547.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 1444]  [Cited by in F6Publishing: 1595]  [Article Influence: 122.7]  [Reference Citation Analysis (0)]
32.  Sivan A, Corrales L, Hubert N, Williams JB, Aquino-Michaels K, Earley ZM, Benyamin FW, Lei YM, Jabri B, Alegre ML. Commensal Bifidobacterium promotes antitumor immunity and facilitates anti-PD-L1 efficacy. Science. 2015;350:1084-1089.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 1979]  [Cited by in F6Publishing: 2594]  [Article Influence: 288.2]  [Reference Citation Analysis (0)]
33.  Pferschy-Wenzig EM, Koskinen K, Moissl-Eichinger C, Bauer R. A Combined LC-MS Metabolomics- and 16S rRNA Sequencing Platform to Assess Interactions between Herbal Medicinal Products and Human Gut Bacteria in Vitro: a Pilot Study on Willow Bark Extract. Front Pharmacol. 2017;8:893.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 21]  [Cited by in F6Publishing: 17]  [Article Influence: 2.4]  [Reference Citation Analysis (0)]
34.  Yao S, Yagi S, Uozumi R, Iida T, Nagao M, Okamura Y, Anazawa T, Okajima H, Kaido T, Uemoto S. A High Portal Venous Pressure Gradient Increases Gut-Related Bacteremia and Consequent Early Mortality After Living Donor Liver Transplantation. Transplantation. 2018;102:623-631.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 7]  [Cited by in F6Publishing: 7]  [Article Influence: 1.2]  [Reference Citation Analysis (0)]
35.  Song JJ, Tian WJ, Kwok LY, Wang YL, Shang YN, Menghe B, Wang JG. Effects of microencapsulated Lactobacillus plantarum LIP-1 on the gut microbiota of hyperlipidaemic rats. Br J Nutr. 2017;118:481-492.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 51]  [Cited by in F6Publishing: 68]  [Article Influence: 9.7]  [Reference Citation Analysis (0)]
36.  Smith PF, Langworthy TA, Mayberry WR. Distribution and composition of lipopolysaccharides from mycoplasmas. J Bacteriol. 1976;125:916-922.  [PubMed]  [DOI]  [Cited in This Article: ]
37.  Zinkernagel MS, Zysset-Burri DC, Keller I, Berger LE, Leichtle AB, Largiadèr CR, Fiedler GM, Wolf S. Association of the Intestinal Microbiome with the Development of Neovascular Age-Related Macular Degeneration. Sci Rep. 2017;7:40826.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 99]  [Cited by in F6Publishing: 118]  [Article Influence: 16.9]  [Reference Citation Analysis (0)]
38.  Suriano F, Bindels LB, Verspreet J, Courtin CM, Verbeke K, Cani PD, Neyrinck AM, Delzenne NM. Fat binding capacity and modulation of the gut microbiota both determine the effect of wheat bran fractions on adiposity. Sci Rep. 2017;7:5621.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 42]  [Cited by in F6Publishing: 37]  [Article Influence: 5.3]  [Reference Citation Analysis (0)]
39.  Horie M, Miura T, Hirakata S, Hosoyama A, Sugino S, Umeno A, Murotomi K, Yoshida Y, Koike T. Comparative analysis of the intestinal flora in type 2 diabetes and nondiabetic mice. Exp Anim. 2017;66:405-416.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 58]  [Cited by in F6Publishing: 76]  [Article Influence: 10.9]  [Reference Citation Analysis (0)]
40.  Collins J, Robinson C, Danhof H, Knetsch CW, van Leeuwen HC, Lawley TD, Auchtung JM, Britton RA. Dietary trehalose enhances virulence of epidemic Clostridium difficile. Nature. 2018;553:291-294.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 278]  [Cited by in F6Publishing: 243]  [Article Influence: 40.5]  [Reference Citation Analysis (0)]
41.  Miller WG, Parker CT, Rubenfield M, Mendz GL, Wösten MM, Ussery DW, Stolz JF, Binnewies TT, Hallin PF, Wang G. The complete genome sequence and analysis of the epsilonproteobacterium Arcobacter butzleri. PLoS One. 2007;2:e1358.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 166]  [Cited by in F6Publishing: 159]  [Article Influence: 9.4]  [Reference Citation Analysis (0)]
42.  Freeman AF, Holland SM. Persistent bacterial infections and primary immune disorders. Curr Opin Microbiol. 2007;10:70-75.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 26]  [Cited by in F6Publishing: 28]  [Article Influence: 1.6]  [Reference Citation Analysis (0)]
43.  Berer K, Gerdes LA, Cekanaviciute E, Jia X, Xiao L, Xia Z, Liu C, Klotz L, Stauffer U, Baranzini SE. Gut microbiota from multiple sclerosis patients enables spontaneous autoimmune encephalomyelitis in mice. Proc Natl Acad Sci. 2017;114:10719-10724.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 470]  [Cited by in F6Publishing: 581]  [Article Influence: 83.0]  [Reference Citation Analysis (0)]
44.  Wu W, Lv L, Shi D, Ye J, Fang D, Guo F, Li Y, He X, Li L. Protective Effect of Akkermansia muciniphila against Immune-Mediated Liver Injury in a Mouse Model. Front Microbiol. 2017;8:1804.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 114]  [Cited by in F6Publishing: 170]  [Article Influence: 24.3]  [Reference Citation Analysis (0)]
45.  Kuroda Y, Ono N, Akaogi J, Nacionales DC, Yamasaki Y, Barker TT, Reeves WH, Satoh M. Induction of lupus-related specific autoantibodies by non-specific inflammation caused by an intraperitoneal injection of n-hexadecane in BALB/c mice. Toxicology. 2006;218:186-196.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 17]  [Cited by in F6Publishing: 17]  [Article Influence: 0.9]  [Reference Citation Analysis (0)]
46.  Jiao N, Baker SS, Chapa-Rodriguez A, Liu W, Nugent CA, Tsompana M, Mastrandrea L, Buck MJ, Baker RD, Genco RJ. Suppressed hepatic bile acid signalling despite elevated production of primary and secondary bile acids in NAFLD. Gut. 2017;pii:gutjnl-2017-314307.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 389]  [Cited by in F6Publishing: 454]  [Article Influence: 75.7]  [Reference Citation Analysis (0)]
47.  Han J, Dzierlenga AL, Lu Z, Billheimer DD, Torabzadeh E, Lake AD, Li H, Novak P, Shipkova P, Aranibar N. Metabolomic profiling distinction of human nonalcoholic fatty liver disease progression from a common rat model. Obesity (Silver Spring). 2017;25:1069-1076.  [PubMed]  [DOI]  [Cited in This Article: ]  [Cited by in Crossref: 36]  [Cited by in F6Publishing: 37]  [Article Influence: 5.3]  [Reference Citation Analysis (0)]
48.  Hashemi Goradel N, Eghbal MA, Darabi M, Roshangar L, Asadi M, Zarghami N, Nouri M. Improvement of Liver Cell Therapy in Rats by Dietary Stearic Acid. Iran Biomed J. 2016;20:217-222.  [PubMed]  [DOI]  [Cited in This Article: ]