Article Text

Original article
Effects of dietary fat on gut microbiota and faecal metabolites, and their relationship with cardiometabolic risk factors: a 6-month randomised controlled-feeding trial
Free
  1. Yi Wan1,
  2. Fenglei Wang1,2,
  3. Jihong Yuan3,
  4. Jie Li3,
  5. Dandan Jiang3,
  6. Jingjing Zhang4,
  7. Hao Li1,
  8. Ruoyi Wang1,2,
  9. Jun Tang1,
  10. Tao Huang5,
  11. Jusheng Zheng6,
  12. Andrew J Sinclair7,
  13. Jim Mann8,
  14. Duo Li1,9
  1. 1 Department of Food Science and Nutrition, Zhejiang University, Hangzhou, China
  2. 2 Department of Nutrition, Harvard T. H. Chan School of Public Health, Boston, USA
  3. 3 No. 1 Department of Nutrition, Chinese People’s Liberation Army General Hospital, Beijing, China
  4. 4 Department of Gastroenterology, School of Medicine, Zhejiang University, Hangzhou, China
  5. 5 Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China
  6. 6 Institute of Basic Medical Science, Westlake University, Hangzhou, China
  7. 7 Department of Nutrition, Dietetics and Food, Monash University, Melbourne, Australia
  8. 8 Department of Human Nutrition and Medicine, University of Otago, Otago, New Zealand
  9. 9 Institute of Nutrition and Health, Qingdao University, Qingdao, China
  1. Correspondence to Professor Duo Li, Institute of Nutrition and Health, Qingdao University, Qingdao 266071, China; duoli{at}qdu.edu.cn

Abstract

Objective To investigate whether diets differing in fat content alter the gut microbiota and faecal metabolomic profiles, and to determine their relationship with cardiometabolic risk factors in healthy adults whose diet is in a transition from a traditional low-fat diet to a diet high in fat and reduced in carbohydrate.

Methods In a 6-month randomised controlled-feeding trial, 217 healthy young adults (aged 18–35 years; body mass index <28 kg/m2; 52% women) who completed the whole trial were included. All the foods were provided during the intervention period. The three isocaloric diets were: a lower-fat diet (fat 20% energy), a moderate-fat diet (fat 30% energy) and a higher-fat diet (fat 40% energy). The effects of the dietary interventions on the gut microbiota, faecal metabolomics and plasma inflammatory factors were investigated.

Results The lower-fat diet was associated with increased α-diversity assessed by the Shannon index (p=0.03), increased abundance of Blautia (p=0.007) and Faecalibacterium (p=0.04), whereas the higher-fat diet was associated with increased Alistipes (p=0.04), Bacteroides (p<0.001) and decreased Faecalibacterium (p=0.04). The concentration of total short-chain fatty acids was significantly decreased in the higher-fat diet group in comparison with the other groups (p<0.001). The cometabolites p-cresol and indole, known to be associated with host metabolic disorders, were decreased in the lower-fat diet group. In addition, the higher-fat diet was associated with faecal enrichment in arachidonic acid and the lipopolysaccharide biosynthesis pathway as well as elevated plasma proinflammatory factors after the intervention.

Conclusion Higher-fat consumption by healthy young adults whose diet is in a state of nutrition transition appeared to be associated with unfavourable changes in gut microbiota, faecal metabolomic profiles and plasma proinflammatory factors, which might confer adverse consequences for long-term health outcomes.

Trial registration number NCT02355795; Results.

  • diet
  • nutrition
  • intestinal microbiology
  • gut inflammation
  • clinical trials

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Significance of this study

What is already known on this subject?

  • Rates of obesity and other cardiometabolic disorders in most developing countries have increased rapidly in parallel with a transition from the traditional lower-fat diet to a diet relatively high in fat but lower in carbohydrate.

  • Gut microbiota dysbiosis has been shown to be associated with a high risk of obesity, type 2 diabetes and many other cardiometabolic diseases.

  • Studies using rodent models suggest that a high-fat diet unbalances the gut microbiota and impairs the gut barrier, resulting in cardiometabolic diseases.

  • Observational studies conducted in globally distinct human populations suggest that diet has a strong effect on the composition of gut microbiota and related faecal metabolomics, while evidence from dietary intervention studies suggests that the effect of diet on gut microbiota might be more modest.

Significance of this study

What are the new findings?

  • At the phylum level, the moderate-fat and higher-fat diets decreased the ratio of Firmicutes to Bacteroidetes after the intervention.

  • At the genus level, the higher-fat diet decreased the abundance of Faecalibacterium, increased the abundance of Alistipes and Bacteroides, while the lower-fat diet increased the Faecalibacterium and Blautia abundance after the intervention.

  • Concentration of faecal butyric acids was increased after the lower-fat diet intervention and decreased after the higher-fat diet intervention. The cometabolites p-cresol and indole were decreased after consumption of the lower-fat diet, while palmitic acid (C16:0), stearic acid (C18:0) and arachidonic acid were increased after the higher-fat diet consumption.

  • Change in relative abundance of Blautia was negatively associated with the changes in serum total cholesterol, low-density lipoprotein cholesterol and non-high-density lipoprotein cholesterol, whereas the change in Bacteroides abundance was positively correlated with the changes in these three blood lipid markers.

  • For circulating proinflammatory factors, the higher-fat diet was associated with increased plasma high sensitivity C-reactive protein and thromboxane B2 relative to the lower-fat diet, while leukotriene B4 and prostaglandin E2 in the lower-fat diet were decreased the most among the three diet groups.

How might it impact on clinical practice in the foreseeable future?

  • Our findings demonstrate that among healthy young adults whose diet is in a state of transition, a higher-fat diet had a notably unfavourable impact on gut microbial taxa, faecal metabolomic profiles and plasma proinflammatory biomarkers, whereas the lower-fat diet was associated with a more favourable profile of these biomarkers. These findings provide confirmatory evidence that nutritional guidelines in countries in a state of nutrition transition should advise against increasing intakes of dietary fat. The results might also have relevance in developed countries in which fat intake is already high.

Introduction

The parallel increases in dietary carbohydrate and the prevalence of obesity and type 2 diabetes (T2D) in the United States and some other Western countries have led to the suggestion that this dietary change, resulting from an increase in carbohydrate consumption and a reduction in total fat intake, might be a determinant of cardiometabolic disorders.1–3 By contrast in China, the nutritional transition from the traditional low-fat, high-carbohydrate diet to a diet relatively higher in fat and lower in carbohydrate has been associated with a dramatic increase in the risk of obesity, T2D and cardiovascular diseases in the past 30 years.4 5 We have previously shown in a randomised controlled-feeding trial that among healthy young adults, a lower-fat, higher-carbohydrate diet is likely to be associated with a lower risk of excessive weight gain and increase in waist circumference and a more favourable lipid profile than  a higher-fat, lower-carbohydrate diet.6

There has been interest in the role of gut microbiota in the development of obesity and cardiometabolic diseases.7 Rodent studies have suggested that the high-fat Western-type diet has a strong effect on the genetic composition and metabolic activity of gut microbiota.8 9 Evidence has shown that in humans gut microbiota diversity and richness are reduced when comparing such high-fat diets with more traditional diets with relatively higher proportions of carbohydrate.10 11 Such diet-induced ‘dysbiosis’ in gut-associated microbial communities has been postulated as a major trigger of metabolic impairments associated with obesity.12 Much of the relevant existing literature on humans is based on observational studies and the relatively few short-term dietary intervention trials suggest that the effect of dietary change might be relatively modest.13–16 No information is available relating to healthy young people or to populations whose diet is in transition from that traditionally consumed to one associated with a high risk of cardiometabolic disease.

Thus in a 6-month randomised controlled-feeding trial among healthy young adults, we compared a number of attributes of gut microbiota and faecal metabolomic profiles as well as markers of inflammation for diets of differing proportions of fat and carbohydrate.

Methods

Study design and population

This study investigated the gut microbiota, faecal metabolomic profiles response and their relationship with cardiometabolic risk factors in 217 healthy adults who had provided faecal samples in the Optimal Dietary Macronutrient Distribution in China trial.6 The trial was conducted at People’s Liberation Army General Hospital in north China and Zhejiang University in south China and was approved by the ethics committees. Each participant provided written informed consent and received no financial compensation or gifts. The trial was registered at ClinicalTrials.gov, number NCT02355795.

The primary aim of the trial was to determine whether the traditional lower-fat diet was more effective than a higher-fat diet at  weight control among a healthy young population whose diet was in a state of nutrition transition. The study design, method and primary results have been described elsewhere.6 17 In brief, the trial involved screening 1145 potentially eligible adults in two centres from January 2015. The main inclusion criteria included body mass index <28 kg/m2, the lower cut-off point according to the Chinese obesity criteria. Key exclusion criteria were blood pressure ≥140 mm Hg systolic or ≥90 mm Hg diastolic, total cholesterol (TC) ≥6.19 mmol/L, low-density lipoprotein cholesterol (LDL-C) ≥4.12 mmol/L, fasting triglyceride ≥2.25 mmol/L and fasting glucose ≥6.11 mmol/L. Finally, 307 eligible participants were randomised (1:1:1) to one of three diets with differing proportions of dietary fat, stratified by study centre, age, sex and body mass index using a computer-generated random number list. After 6 months of controlled-feeding intervention, 245 participants (79.8%) completed the whole trial. Among them, a total of 217 participants who provided faecal samples at baseline and the end of the trial were included in the 16S rRNA sequencing. Additionally, 120 participants (n=40 in each diet group) with enough faecal samples after 16S rRNA sequencing were further included in the faecal metabolomic analysis.

Study diets and intervention

The three diets were isocaloric, the primary distinguishing feature being their fat and carbohydrate content. We achieved the required macronutrient distribution in the three diet groups by replacing a proportion of energy derived from carbohydrate (mainly white rice and wheat flour, the most consumed carbohydrate sources in Asia) with fats (mainly soybean oil, the most widely used edible oil in Asia). Otherwise, the menus for the three diet groups were almost identical. The intervention diets were a lower-fat diet (fat 20% and carbohydrate 66% energy, corresponding to the macronutrient distribution 30 years ago during which obesity was rare in China); a moderate-fat diet (fat 30% and carbohydrate 56% energy, comparable with current macronutrient intake in China and also the upper limit of fat intake recommended by the Chinese Nutrition Society) and a higher-fat diet (fat 40% and carbohydrate 46% energy, approximating the current consumption of Chinese residents of some megacities). Protein provided 14% energy in all three diets. In order to avoid an additional dietary variable, intake of dietary fibre on all the three diets was maintained at the baseline level of consumption, around 14 g/day. Details of the study diets are provided in the online supplementary table 1. Participants were provided with all of their food and most beverages throughout the intervention. During the controlled-feeding period, they were asked to complete a daily diary in which they recorded whether they had eaten all the study foods and provided a list of non-study foods eaten.

Supplementary file 1

Blood specimen collection and assessment of proinflammatory factors

Blood samples were collected by venepuncture after a 12-hour fast. Plasma and serum were centrifuged from the blood samples and immediately stored at −80°C for further analysis. Serum concentrations of clinical biomarkers, including lipid profiles, glucose and insulin, were measured as described previously.6 Plasma proinflammatory factor concentrations, including interleukin (IL)−1β, IL-6, IL-8, tumour necrosis factor-α (TNF-α), prostaglandin E2 (PGE2), thromboxane B2 (TXB2) and leukotriene B4 (LTB4), were measured using ELISA kits (Cayman Chemical Co, Ann Arbour, Michigan, USA and MultiSciences Biotech Co, Hangzhou, China). In addition, high-sensitive C-reactive protein (hs-CRP) was measured by a latex-enhanced turbidimetric immunoassay method. All measurements were performed at the end of the study to minimise variability.

16S rRNA sequencing and measurement of faecal metabolomics

Faecal samples were collected at baseline and the end of the 6-month intervention. Each faecal sample was snap frozen in liquid nitrogen within minutes of donation and then kept at −80°C. DNA was extracted from stool samples (200 mg) using QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions, with additional homogenised steps in a bead beater (FastPrep, Thermo Electron Co, Boston, Massachusetts, USA). The V3-V4 hypervariable regions of the bacteria 16S rRNA gene were amplified with barcode-indexed primers 338F (5’-ACTCCTACGGGAGGCAGCAG-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’). The purified amplicons were pooled in equimolar concentration and further paired-end sequencing was performed using an Illumina Miseq instrument (Illumina, San Diego, California, USA). Raw sequencing data have been deposited in the National Center for Biotechnology Information Sequence Read Archive database with accession number PRJNA480547. Targeted identification and measurement of faecal metabolomics were performed based on a validated method.18 Details of the 16S rRNA sequence data preparation and analysis as well as the quantitative measurement of faecal metabolomics are available in online supplementary materials.

Statistical analyses

Descriptive statistics included crude comparison of demographic factors, clinical parameters and dietary intake among the three diet groups. Differences among diets for 6-month changes from baseline of proinflammatory markers were tested by analysis of variance or Kruskal-Wallis test based on whether the data were normally distributed. The α-diversity indices evaluating gut microbial community richness (the Ace and Chao1 estimators) and community diversity (the Shannon estimator) were calculated using Mothur.19 Principal coordinate analysis (PCoA) based on Bray-Curtis distance and permutational multivariate analysis of variance (PERMANOVA) were performed to compare the global microbiota composition before and after intervention in each group at phylum, genus and operational taxonomic unit (OTU) levels, respectively. Enterotypes were determined using the abundance of genus and OTUs separately using the methods suggested by Arumugam et al.20 In univariate analysis of gut microbiota and predicted Kyoto Encyclopaedia of Genes and Genomes (KEGG) biochemical pathways in each group, a paired t-test or a Wilcoxon matched-pairs test was adopted and p values were adjusted for multiple comparison using the Benjamini-Hochberg false discovery rate (FDR). Correlations between changes in cardiometabolic risk factors, including body weight, waist and lipid profiles selected based on the previous results,6 and changes in genus relative abundance, were calculated using Spearman’s rank or Pearson’s correlation depending on the distribution of the data.

Orthogonal projection to latent structure-discriminant analysis (OPLS-DA) with unit variance scaling was performed to examine the overall microbial metabolites distribution before and after intervention in each group. The qualities of all OPLS-DA models were assessed with R2X (the total variation being explained by the model) and Q2 (denoting the predictability of the model). The significance of the models was further validated by cross-validation analysis of variance (CV-ANOVA). Univariate analysis was also conducted with p values adjusted for multiple testing by the Benjamini-Hochberg approach. Correlations between changes in faecal metabolites and changes in genus relative abundance, as well as proinflammatory markers were calculated by Spearman’s rank test or Pearson’s correlation test, when appropriate, in the full cohort. P values for the correlations between changes in faecal metabolites and changes in genus relative abundance were corrected for multiple comparisons.

Analyses were also conducted to examine the difference in gut microbiota and faecal metabolomic profiles across groups. PCoA based on Bray-Curtis distance and PERMANOVA were performed to compare the global microbiota composition among the three groups at phylum, genus and OTU levels, at baseline and after intervention, respectively. OPLS-DA with unit variance scaling and CV-ANOVA were performed to examine the overall microbial metabolites distribution across the three groups at baseline and after intervention. For the significantly changed bacteria and faecal metabolites after FDR correction in each group, Wilcoxon test was further performed to compare the changes in values between the dietary groups.

OPLS-DA and CV-ANOVA were performed with the software SIMCA-P +version 14.0, and all other statistical analyses were performed in R version 3.4. The main R packages used were ‘ape’ and ‘vegan’ for PCoA and PERMANOVA analyses and ‘ade4’ and ‘clusterSim’ for enterotype analysis.

Results

Characteristics of study participants

No significant differences in anthropometric measurements or clinical parameters were observed among the three groups at baseline (n=217, table 1). Mean carbohydrate intake of the three groups was 55% energy at baseline and during the intervention period 67% energy in the lower-fat diet group, 57% energy in the moderate-fat diet group and 48% energy in the higher-fat diet group (online supplementary figure 1). Mean dietary fat consumption of the three groups was 32% at baseline, which changed to around 20%, 30% and 40% energy in the lower-fat diet group, moderate-fat diet group and higher-fat diet group, respectively, during the intervention period. Protein intake remained stable, at around 13–14% energy. After the 6-month controlled-feeding intervention, all groups lost weight, the reduction in the lower-fat diet group being significantly greater than in the higher-fat groups. Similarly, reduction in waist circumference, TC, high-density lipoprotein cholesterol (HDL-C), LDL-C and non-HDL-C for the lower-fat diet group were significantly greater than those observed in the higher-fat diet group. However, the change in ratio of TC to HDL-C among the groups did not differ significantly throughout the intervention period. Changes in anthropometric data and clinical parameters during the 6-month intervention among participants included in the gut microbiota study (n=217) are presented in online supplementary figure 2.

Table 1

Baseline characteristics of the participants included in the gut microbiota study (n=217)

Effects of dietary interventions on gut microbiota

The microbial community richness indicated by Ace and Chao1 estimators showed no significant changes among the three groups, whereas the community diversity estimated by the Shannon index was significantly increased in the lower-fat diet group relative to the higher-fat diet group (p=0.03) (figure 1A). No significant shift was observed on the three diets in the overall composition of the gut microbiota at phylum, genus and OTU levels as indicated by PERMANOVA analysis (figure 1B-D). Similarly, the enterotypes did not change significantly in any of the three groups (figure 2 and online supplementary figure 3).

We further performed PCoA and PERMANOVA analyses among groups to investigate the potential differences at baseline and month 6. At phylum level, we observed no statistically significant differences among the three diet groups at baseline, but marginally significant differences at month 6 (p=0.08). However, at genus and OTU level, there were no significant differences among groups at baseline or month 6 (online supplementary figure 4). Univariate analysis at phylum level showed that the relative abundance of Bacteroidetes was significantly increased in the moderate-fat diet group after intervention (p=0.005 after FDR correction) (figure 3A). In the higher-fat diet group, the relative abundance of Firmicutes was decreased, while Bacteroidetes was increased after intervention (p<0.001 after FDR correction) (figure 3A). The ratio of Firmicutes to Bacteroidetes was significantly decreased in the moderate-fat and higher-fat diet group after intervention (p=0.004 for moderate fat diet group and p<0.001 for higher-fat diet group). Compared with the lower-fat diet, the higher-fat diet significantly increased the abundance of Bacteroidetes and decreased the abundance of the Firmicutes (all p<0.01, online supplementary figure 5). At genus level, the lower-fat diet resulted in a significant increase in relative abundance of Blautia (p=0.007) and Faecalibacterium (p=0.04) after the intervention; while the higher-fat diet led to a reduction in Faecalibacterium (p=0.04) and an increase in Alistipes (p=0.04) and Bacteroides (p<0.001), after FDR correction for multiple testing (figure 3B and online supplementary table 2). Compared with lower-fat diet, the higher-fat diet led to a significant decrease in Blautia and Faecalibacterium abundance, while there was a significant increase in the abundance of Alistipes and Bacteroides (all p<0.01, online supplementary figure 5). The change in relative abundance of Blautia was negatively associated with the changes in serum TC (r=−0.27, p<0.001), LDL-C (r=−0.26, p<0.001) and non-HDL-C (r=−0.24, p<0.001), whereas the change in Bacteroides was positively correlated with the changes in these three blood lipid markers (r=0.44, p<0.001 for TC; r=0.34, p<0.001 for LDL-C and r=0.40, p<0.001 for non-HDL-C, respectively) (figure 3C).

No significant changes in KEGG database biochemical pathways were seen in the lower-fat and moderate-fat groups after intervention (online supplementary table 3). Among the nine most responding pathways (FDR corrected p<0.10) in the lower-fat group, lipopolysaccharide biosynthesis (ko00540) decreased to a notable degree (median fold change after intervention=0.75). In the higher-fat group, 29 pathways were identified as significantly responding to the higher-fat intervention after FDR correction (FDR corrected p<0.05). Among these, pathways with notable changes (median fold change after intervention ≤0.80 or ≥1.20) were steroid hormone biosynthesis (ko00140), lysosome pathway (ko04142), arachidonic acid metabolism (ko00590) and lipopolysaccharide biosynthesis (ko00540). All four pathways were significantly increased in the higher-fat group.

Figure 1

Changes of global gut microbiota after intervention in each group. (A) α-Diversity at genus level estimated by Ace, Chao 1 and Shannon estimator. *p<0.05, **p<0.01. (B) Principal coordinate analysis (PCoA) score plots based on Bray-Curtis distance at phylum level. (C) PCoA score plots based on Bray-Curtis distance at genus level. (D) PCoA score plots based on Bray-Cutis distance at operational taxonomic unit (OTU) level.

Figure 2

Clustering of gut microbiota into enterotypes at genus-level using Jensen-Shannon distance. (A) Two clusters exist most naturally in the dataset by the PAM method. The x-axis shows the cluster number, the y-axis shows the Calinski-Harabasz (CH) index. (B) Clustering of the first two principal components. (C) Relative abundance of bacterial taxa characteristic of each enterotype. Boxes represent the IQR and the line within represents the median. Whiskers denote the lowest and highest values within 1.5×IQR. (D) Proportions of enterotypes in each group before and after the intervention. No statistically significant changes were seen in each group before and after the intervention by χ2 tests.

Figure 3

Analyses based on individual genera. (A) Changes of relative abundance of individual phylum before and after diet intervention. (B) Heatmap of the median fold changes in relative abundance of individual genera before and after the intervention. (C) Significant associations between changes in lipid profiles and changes in genus abundance as measured by the Spearman’s correlations. **p<0.01. FDR, false discovery rate.

Effects of dietary interventions on faecal metabolomic profiles

Significant shifts in the composition of faecal metabolites were observed after the lower-fat and higher-fat diets intervention (p=0.02 for lower-fat and p=0.01 for higher-fat diet) (figure 4A). Furthermore, no statistically significant differences were seen among groups at baseline (online supplementary figure 6). However, we observed significant differences among groups at month 6 (p=0.002, online supplementary figure 6). Further stratified analysis by metabolite categories showed that amino acid metabolites were significantly altered in the lower-fat group (p=0.001) and long-chain saturated fatty acids (with 14 or more carbons) were significantly altered in the higher-fat group (p<0.001) (online supplementary figure 7). In addition, total short-chain fatty acid (SCFAs) concentrations were significantly decreased in the higher-fat group, compared with other diet groups (p<0.001) (figure 4B). Individual metabolite analysis identified four significantly changed faecal metabolites, including three amino acid metabolites (p-cresol, indole and 3-indolepropionic acid) and one SCFA (butyric acid), in the lower-fat group after FDR correction (figure 4C and online supplementary table 4). Among these four metabolites, two were decreased (p-cresol and indole) and two were increased (3-indolepropionic acid and butyric acid). While in the higher-fat group, eight faecal metabolites including long-chain saturated fatty acids, SCFAs and amino acid metabolites were found to be significantly changed. Among them, three were decreased (butyric acid, valeric acid and ethylmethylacetic acid) and five were increased (indole, palmitic acid, stearic acid, arachidonic acid and indoleacetic acid). Compared with the lower-fat group, the higher-fat diet group significantly increased the faecal concentration of palmitic acid, stearic acid, arachidonic acid, indoleacetic acid, indole and p-cresol (all p<0.01), while significantly decreasing the faecal concentration of butyric acid, valeric acid and 3-indolepropionic acid (all p<0.05, online supplementary figure 8). Since arachidonic acid is the precursor of eicosanoids and other lipid mediators involved in the inflammation process, we specifically tested the correlation between changes in faecal arachidonic acid concentration and changes in plasma proinflammatory markers concentrations. The change in faecal concentration of arachidonic acid was found to be positively associated with the changes in plasma concentrations of PGE2 (r=0.35, p<0.001), TXB2 (r=0.28, p=0.002) and hs-CRP (r=0.20, p=0.03) (figure 4D). Additionally, Spearman correlation analysis showed that the change in Bacteroides was positively associated with the changes in palmitic acid (r=0.35, p=0.004), indole (r=0.39, p<0.001) and p-cresol (r=0.64, p<0.001); the change in Faecalibacterium was positively associated with the change in butyric acid (r=0.35, p=0.004); whereas the change in Ruminococcus was negatively associated with the change in p-cresol (r=0.32, p=0.01, figure 5).

Figure 4

Faecal metabolomics changed after the intervention. (A) OPLS-DA score plots comparing the faecal metabolites before and after the intervention in each group. P values are evaluated by CV-ANOVA. (B) Faecal total short-chain fatty acids concentration before and after the intervention. *p<0.05 and **p<0.01. (C) Heatmap of the median fold changes in faecal metabolite concentrations before and after the intervention. (D) Associations between change in faecal arachidonic acid concentration and changes in plasma proinflammatory marker levels. CV-ANOVA, cross-validation analysis of variance; FDR, false discovery rate; hs-CRP, high sensitivity C-reactive protein; IL, interleukin; LTB4, leukotriene B4; OPLS-DA, orthogonal projection to latent structure-discriminant analysis; PGE2, prostaglandin E2; TNF-α, tumour necrosis factor α; TXB2, thromboxane B2.

Figure 5

Correlation between changes in faecal metabolites and changes in genus abundance. The heat-map shows the Spearman correlation coefficient between changes in faecal metabolite concentrations and changes in relative abundance of individual genera. The intensity of the colours represents the degree of association between changes in the concentration of faecal metabolites and changes in relative abundance of individual genera as measured by the Spearman’s correlations. P values were corrected for multiple testing using the Benjamini-Hochberg false discovery rate. *p<0.05, **p<0.01. TSCFA, total short-chain fatty acids.

Effects of dietary interventions on inflammatory markers

Plasma concentration of hs-CRP was significantly increased during the higher-fat diet intervention, compared with the moderate-fat (p for higher-fat vs moderate-fat <0.001) and lower-fat (p for higher-fat vs lower-fat <0.001) diet groups (table 2). Likewise, compared with the lower-fat diet, plasma concentration of the inflammatory mediators TXB2 was increased on the higher-fat diet (p for higher-fat vs lower-fat=0.02). Additionally, LTB4 and PGE2 were significantly decreased in the lower-fat diet group, when compared with moderate-fat diet group (p for lower-fat vs moderate-fat=0.006, 0.003, respectively) and higher-fat diet group (p for lower-fat vs higher-fat=0.02, <0.001, respectively). However, the changes in proinflammatory cytokines including IL-1β, IL-6, IL-8 and TNF-α showed no significant differences among the three groups after the 6-month diet intervention.

Table 2

Effects of dietary interventions on plasma proinflammatory markers (n=217)

Discussion

We believe that the present randomised, controlled-feeding trial, is the first study to investigate the effect of different proportions of dietary fat intake on gut microbiota, faecal metabolomic profiles and plasma inflammatory factors in healthy young adults drawn from a population whose diet is in a state of transition. The findings of our study suggest that although the three diets did not induce major changes of the global gut microbiota composition, the higher-fat diet had overall unfavourable effects on several important biomarkers, while the lower-fat diet with a macronutrient distribution corresponding to those consumed in China until some 30 years ago, appeared to have beneficial effects.

Numerous rodent studies have reported that a high-fat diet might appreciably unbalance the gut microbiota composition.21 However, the effect of diet on the human gut is complex given the considerable individual variation in the response of the gut microbiota to dietary intervention.22 23 Our findings are in keeping with previous studies, which indicate that the overall structure of gut microbiota is principally influenced by long-term dietary intake.10 13–15 23 24 Short-term effects have mainly been reported in studies with small sample size,25 26 and characterised by extreme and unrealistic dietary intakes25 or when the subjects studied were at high cardiometabolic risk.27–29 A 10-day clinical trial with nine participants reported that the composition and function of the microbiota were rapidly changed when carbohydrates were eliminated from the diet, a dietary intervention which is neither sustainable nor desirable in the long term.25 In much more moderate dietary interventions that can be sustained in humans, high interindividual variations among subjects might have masked minor changes in distinct taxonomic groups.30 It has been reported that diet might have a greater influence on the metabolic activity of the gut microbiota than on its taxonomic composition.16 In accordance with this proposition, in our study diets differing in dietary fat content provoked a marked shift of faecal concentrations of metabolites produced by the gut microbiota.

We found Blautia and Faecalibacterium, genera known to contain butyrate-producing bacteria, were increased in the lower-fat diet group. Blautia is a group of bacteria containing various acetate and butyrate producers and reported to have a lower relative abundance in patients with T2D than in healthy controls.31 32 Faecalibacterium contains anti-inflammatory and functionally important bacteria.33 Higher abundance of F. prausnitzii was associated with a reduction of systematic low-grade inflammation among obese subjects receiving bariatric surgery.34 Many previous studies have reported the benefits of dietary fibre on gut microbiota.35 In our study, the amount of carbohydrate was highest in the lower-fat diet group, mainly from white rice and wheat flour (bread). The proportion of resistant starch has been reported to be 3% in white rice and 2% in bread.36 37 Given that dietary fibre intake did not differ among the three groups, it is possible that the favourable effects of the lower-fat diet in our study might be due to the increased amount of resistant starch which can, like dietary fibre, be fermented by the gut microbiota, with associated health benefits.

In contrast, the abundance of Faecalibacterium was decreased and Bacteroides and Alistipes were increased after the higher-fat diet intervention. Among the Chinese population, Bacteroides and Alistipes were found to be more abundant in patients with T2D than in subjects with normal glucose metabolism.38 These changes at genus level indicated that dietary fat content had a selective effect on the human gut microbiota, which might have clinical significance among healthy young adults. The higher-fat diet induced a reduction in faecal butyrate acid and SCFAs as opposed to the lower-fat diet group. This is probably due to the reduction in dietary carbohydrate intake.39 It has been suggested that the benefits of a ‘healthy diet’ are mediated via the anti-inflammatory effects of SCFAs and other bioactive compounds produced by the gut microbiota.40 Several observations in our study pointed towards an effect of gut microbiota-mediated proinflammatory factors. Among the individual genera which responded to the dietary intervention, the abundance of Faecalibacterium was directly related to the concentration of butyric acid, as suggested by the positive association between changes in butyrate concentration and the abundance of Faecalibacterium. A study using faecal microbiota transplantation from lean donors to insulin-resistant patients with metabolic syndrome demonstrated that faeces from lean subjects were associated with enhanced numbers of butyrate-producing bacteria and increased insulin sensitivity.41 In contrast, a significant reduction of dietary carbohydrates in omnivores led to much lower levels of faecal SCFAs and decreased bacterial numbers.39 The protective role of SCFAs against different types of disease is well documented; in particular, butyric acid can serve as the energy substrate for epithelial cells of the gut and confer anti-inflammatory effects.42 The reduction in SCFAs induced by the higher-fat diet was presumably detrimental for metabolic health.

In addition to the reduction of SCFAs, the higher-fat diet induced significant alterations in long-chain fatty acids metabolism. The faecal concentrations of palmitic acid (C16:0) and stearic acid (C18:0) were significantly increased in the higher-fat group. Palmitic acid and stearic acid, the main saturated fatty acids, in food and tissues, have been considered to stimulate inflammatory signalling in macrophages, adipocytes, myocytes, hepatocytes, etc.43 Epidemiological studies showed that plasma phospholipid palmitic acid and stearic acid were positively associated with the incidence of T2D and cardiovascular diseases.44 Notably, the predicted lipopolysaccharide biosynthesis and arachidonic acid metabolism pathways were also increased in response to the higher-fat diet. Lipopolysaccharide is known to induce the release of arachidonic acid and its inflammation-involving metabolites, such as prostaglandins, thromboxane and leukotrienes.45 It should be noted that the intake of polyunsaturated fatty acids (PUFAs) was relatively high in the higher-fat diet group (24% of total energy) owing to exclusive use of soybean oil, which is rich in n-6 PUFA. A higher intake of n-6 PUFA has been reported to have proinflammatory effects.46 Consistent with the predicted pathway analysis, we found the faecal concentration of arachidonic acid (the precursor of proinflammatory compounds like PGE2, TXB2 and LTB4) was indeed increased in the higher-fat group. The positive association between changes in faecal arachidonic acid concentration and the changes in plasma levels of inflammatory makers—namely, PGE2 and TXB2, suggested a potential mechanism through which the amount of dietary fat consumption affects inflammation via gut microbiota.

Interestingly, although the same amount and quality of protein was provided in each diet group, we found that the lower-fat diet group induced a significantly alteration in the amino acids metabolism, indicating enhanced protein degradation by the gut microbiota. The results reinforced the view that luminal substrate availability is also an important variable for fixing bacterial metabolite concentration.16 The lower-fat diet group induced reduction of faecal concentrations of indole and p-cresol, the host microbiota cometabolites derived from amino acids, while the higher-fat diet led to an increase in indole concentration. Indole is the precursor of indoxyl sulfate, which has been linked with hypertension and cardiovascular disease in people with chronic kidney disease.47 Likewise, p-cresol, the precursor of p-cresol sulfate, has been shown to act as an inhibitor of colonocyte oxygen consumption and as a genotoxic agent in vitro.48 Production rates of these potentially toxic compounds are found to be markedly lower in vegetarians than in people consuming an unrestricted diet.49 It has also been reported that bacterial diversity is higher among populations adopting a lower-fat diet like a vegetarian diet and can be modified through dietary intervention in obese subjects.50 It should be noted that the microbial α-diversity (indicated by the Shannon index) was increased in the lower-fat group compared with the higher-fat diet group. Increased diversity has been associated with improved health in the elderly, and reduction of gut microbial diversity has been linked to increased risk of gastrointestinal diseases and proinflammatory characteristics.14 51 Thus, in addition to increased SCFAs, our data suggest that the lower-fat diet can decrease the luminal concentration of these deleterious compounds while increasing the gut microbial α-diversity.

Correlation analysis allowed us to identify several new bacterial genera potentially implicated in the host metabolic health. We observed negative associations of Blautia abundance with serum concentration of TC, LDL-C and non-HDL-C, whereas the opposite was seen for Bacteroides, a genus that we found to be increased after higher-fat consumption. A previous clinical trial reported that T2D and hyperlipidaemia patients were ameliorated by metformin and a herbal formula with increased abundance of Bautia, which was correlated significantly with the improvement in circulating glucose and lipids.52 Also, in line with a clinical trial in which Blautia was found to be negatively associated with TC and non-HDL-C, we observed similar negative associations between Blautia and blood lipid profiles.53 The enrichment of Blautia in the lower-fat group suggested that the diet-induced change in gut microbiota might be related to the host lipid homoeostasis in healthy young adults and possibly, be a future target for the management of subjects with cardiometabolic diseases.

A limitation of this study was that faecal sampling was carried out only at baseline and at the end of the trial. A more complete picture of changes in gut microbiota would have been possible with more frequent sampling. In addition, all three groups lost weight and weight loss differed among groups after dietary intervention. Whether weight loss led to the changes in gut microbiota and faecal metabolomics or vice versa needs to be studied further. Another limitation relates to the generalisability of our study since our subjects were healthy young non-obese adults; the results might not apply to people at high risk of cardiometabolic diseases. The strengths of this study included the controlled-feeding design, the large sample size and the relatively long intervention period. In contrast to previous studies, we applied 16S rRNA amplicon profiling and quantitatively targeted bacterial metabolomics, which allowed us to understand both gut microbiota response and bacterial metabolites to gain more information about host–gut microbiota metabolic interactions in response to dietary fat content.

Conclusion

In summary, compared with a lower-fat diet, long-term consumption of a higher-fat diet appears to be undesirable owing to changes in gut microbiota, faecal metabolomic profiles and proinflammatory factors for healthy young adults whose diet is in transition from the traditionally consumed lower fat, higher carbohydrate diet to one characterised by an appreciably higher fat content. These findings might also have relevance in developed countries in which fat intake is already high.

Acknowledgments

We thank all the participants, the investigators, on-site staff, kitchen staff and all members of the Duo Li Laboratory for participating in scientific discussions about the manuscript.

References

Supplementary materials

Footnotes

  • YW and FW contributed equally.

  • Contributors Conception and design of the study: DL, YW, and FW. Collection of data: JL, DJ, JY, JiZ, HL and RW. Analysis and interpretation of data: YW, FW, JT, TH and JuZ. Drafting of the manuscript: YW and FW. Critical revision of the manuscript for important intellectual content: AJS, JM, and DL. Administrative support and study supervision: DL and JY. All authors read, revised and approved the final draft.

  • Funding This study was funded by the National Basic Research Program of China (2015CB553604) and China Postdoctoral Science Foundation (2018M642466). The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the manuscript.

  • Competing interests None declared.

  • Patient consent Obtained.

  • Ethics approval Institutional review board.

  • Provenance and peer review Not commissioned; externally peer reviewed.

Linked Articles