Skip to main content

Calorie restriction improves metabolic state independently of gut microbiome composition: a randomized dietary intervention trial

Abstract

Background

The gut microbiota has been suggested to play a significant role in the development of overweight and obesity. However, the effects of calorie restriction on gut microbiota of overweight and obese adults, especially over longer durations, are largely unexplored.

Methods

Here, we longitudinally analyzed the effects of intermittent calorie restriction (ICR) operationalized as the 5:2 diet versus continuous calorie restriction (CCR) on fecal microbiota of 147 overweight or obese adults in a 50-week parallel-arm randomized controlled trial, the HELENA Trial. The primary outcome of the trial was the differential effects of ICR versus CCR on gene expression in subcutaneous adipose tissue. Changes in the gut microbiome, which are the focus of this publication, were defined as exploratory endpoint of the trial. The trial comprised a 12-week intervention period, a 12-week maintenance period, and a final follow-up period of 26 weeks.

Results

Both diets resulted in ~5% weight loss. However, except for Lactobacillales being enriched after ICR, post-intervention microbiome composition did not significantly differ between groups. Overall weight loss was associated with significant metabolic improvements, but not with changes in the gut microbiome. Nonetheless, the abundance of the Dorea genus at baseline was moderately predictive of subsequent weight loss (AUROC of 0.74 for distinguishing the highest versus lowest weight loss quartiles). Despite the lack of consistent intervention effects on microbiome composition, significant study group-independent co-variation between gut bacterial families and metabolic biomarkers, anthropometric measures, and dietary composition was detectable. Our analysis in particular revealed associations between insulin sensitivity (HOMA-IR) and Akkermansiaceae, Christensenellaceae, and Tanerellaceae. It also suggests the possibility of a beneficial modulation of the latter two intestinal taxa by a diet high in vegetables and fiber, and low in processed meat.

Conclusions

Overall, our results suggest that the gut microbiome remains stable and highly individual-specific under dietary calorie restriction.

Trial registration

The trial, including the present microbiome component, was prospectively registered at ClinicalTrials.govNCT02449148 on May 20, 2015.

Background

Obesity is an important risk factor for chronic diseases such as type 2 diabetes (T2D), cardiovascular diseases (CVD), and several types of cancer [1]. Aside from well-established mechanisms such as obesity-induced inflammation, alterations in sugar and lipid metabolism, and steroid hormone signaling [2,3,4], imbalances in the composition of the gut microbiome have also been linked to the progression of obesity and its cardio-metabolic sequelae [5]. Obesity has been associated with a lower abundance of gut bacteria in the Bacteroidetes phylum, in contrast to a higher abundance of Firmicutes [5]. Furthermore, gut bacteria-dependent turnover of metabolites, e.g., short-chain fatty acids (SCFAs) [6], secondary bile acids [7], trimethylamine-N-Oxide (TMAO) [8] or amino acids [9], acting locally or systemically, may mediate the links between the gut microbiome, obesity, and obesity-associated chronic diseases.

Dietary weight loss may lead to improvements in host metabolism and alleviate the risk of possible obesity-associated comorbidities. Despite evidence to suggest that these benefits are due to modulations of the gut microbiome and related metabolites upon calorie restriction [10,11,12], findings on the effects of dietary weight loss on the composition and function of the microbiome from randomized trials have been equivocal [13]. This may be attributable to small sample sizes, short trial durations, and diverse microbiome assessment techniques in previous studies, but also to the drastic and restrictive dietary methods used to achieve weight loss in many intervention studies [13, 14]. Moreover, distinguishing transient variability in gut microbiome composition [15, 16] from intervention-associated changes has been challenging since the majority of earlier studies lacked repeated gut microbiome assessments. In addition, it is not known whether some reported initial changes in gut microbiome composition after weight loss persist in the long term, and to what extent this relates to variations in body composition, anthropometric, and clinical biomarkers as well as circulating metabolites.

While mixed results have been reported from the abovementioned studies on different types of continuous calorie restriction (CCR) and the microbiome, intermittent calorie restriction (ICR) has recently been shown to improve several health outcomes in animal models, possibly via effects on the microbiome [17,18,19,20,21]. However, despite these promising findings from animal models and first smaller studies among humans [13], there is a lack of well-powered trials on ICR in relation to the microbiome. Thus, in the present study, we investigated whether ICR (operationalized as the 5:2 diet) or CCR induced alterations in the gut microbiome and to which extent these were associated with overall weight loss irrespective of the dietary intervention in overweight or obese adults. To this end, we analyzed repeated fecal samples over 1 year, i.e., at baseline, 12, 24, and 50 weeks with 16S rRNA gene sequencing and also performed repeated targeted profiling of gut microbiome-related metabolites in the circulation. This study was conducted using data and samples of the HELENA Trial (NCT02449148), a randomized controlled trial undertaken to compare the metabolic effects of ICR vs. CCR [22].

Methods

Study design and participants

The protocol of the HELENA Trial, a parallel-arm randomized controlled trial which included a 12-week controlled intervention phase with additional dietary counselling, followed by a 12-week maintenance phase (week 24) and a final observational follow-up phase of 26 weeks (week 50) (Fig. 1) has previously been published [23]. Briefly, the trial recruited 150 overweight or obese (BMI > 25 and < 40 kg/m2), non-smoking adults (age range 35–65 years, 50% women), free of any gastrointestinal tract diseases as well as major chronic diseases such as CVD or diabetes in Heidelberg, Germany, in 2015. Recruitment was carried out from 6th May, 2015, to 17th May, 2016. At the time of enrolment, none of the participants were taking antibiotics. Participants were randomly assigned to one of three groups, i.e., an intermittent calorie restriction (ICR) (n = 49), a continuous calorie restriction (CCR) (n = 49), or a control group (CTR) (n = 52) over a 50-week period in a 1:1:1 ratio. Overall, fecal samples were obtained from 147 participants at baseline (Fig. 1). In the ICR group, participants reduced calorie intake by ~75% on two non-consecutive days of the week (so-called “5:2 diet”) whereas those in the CCR group reduced daily calorie intake by 20% [22]. All participants (ICR, CCR, and CTR) were advised to adhere to the dietary recommendations of the German Nutrition Society (DGE e.V.), i.e., 55% energy from carbohydrates, 15% from protein, and 30% from fat [24]. The results on the effects of ICR vs. CCR vs. CTR with respect to the pre-specified primary endpoints, i.e., changes in the expression of 82 pre-selected genes in subcutaneous adipose tissue (SAT), as well as secondary endpoints, i.e., anthropometric, body composition, and routine metabolic blood biomarkers, have already been published [22]. There were no differential effects of ICR and CCR [22]. The present analyses on changes in clinical parameters (secondary endpoints) as well as microbiome and metabolites (exploratory endpoints) across the initial trial arms and by overall weight loss were pre-specified in the study protocol of the HELENA Trial [23]. However, all exploratory analyses of the gut microbiome beyond the comparisons across trial arms and weight-loss quartiles had a post hoc character. These analyses include (i) association tests of gut microbial families with body composition, clinical markers, and dietary intake, (ii) the prediction of subsequent weight loss from baseline microbiome features, as well as (iii) comparisons of microbiome data from the Helena Trial with that from published studies. The ethics committee of the Heidelberg University Hospital (Heidelberg, Germany, Reference: vote S299-2014) as well as the EMBL Bioethics Advisory Committee (Reference: 2020-022) approved the trial, which was designed and conducted in agreement with the principles of the World Medical Association’s Declaration of Helsinki [25]. All participants signed informed consent forms before enrolment. The trial, including the present microbiome component, was registered at clinicaltrials.gov (NCT02449148).

Fig. 1
figure 1

Modified CONSORT diagram showing the flow of participants in the HELENA Trial. ICR – intermittent calorie restriction; CCR – continuous calorie restriction; CTR – control group; CONSORT – consolidated standards of reporting trials

Dietary assessment

Dietary intake was assessed using food records for 7 consecutive days. Information on the time and location of food intake, portion sizes, and details on the characteristics of the meals/foods consumed were obtained with the 7-day food records [22]. Pictured portion size books were provided to all study participants in order to aid with portion size estimation. Energy and nutrient intake were analyzed based on the German Nutrient Database (Bundeslebensmittelschlüssel, version 3.02) using the software PRODI 6.8 (Nutri-Science GmbH, Hausach, Germany). Daily energy and nutrient intake were estimated as the averages of all reported days [22].

Blood sample collection and biochemical analyses

Measurements of routine metabolic biomarkers were carried out at the central laboratory of the Heidelberg University Hospital by clinical standard assays right after fasting blood sampling at baseline, and after 12, 24, and 50 weeks after the start of the trial. Fasting blood samples were subsequently aliquoted and immediately stored at − 80 °C at the Biobank of the National Center for Tumor Diseases (NCT), Heidelberg, Germany, until further biochemical analyses. The concentrations of routine blood-based biomarkers such as LDL cholesterol, HDL cholesterol, total cholesterol, triglycerides, glucose, GGT, ALT, AST, and of serum biomarkers (adiponectin, CRP, IL-8, IL-6, IFN-γ, TNFα, resistin, leptin, insulin, IGF-1) were measured as described in detail previously [22]. The collection of all biospecimens, and subsequent biochemical analyses in the HELENA Trial was conducted in line with standard operating procedures [22].

Plasma SCFA concentrations were determined by high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) at the Institute of Agricultural and Nutritional Sciences, Martin Luther University Halle-Wittenberg using a modified form of the methods by Chan et al. [26] and Zeng and Cao [27]. Further details of the measurements of the SCFAs have been described previously in a related publication [28].

Targeted metabolomics analysis of plasma acylcarnitines and amino acids as their butyl esters was performed on dry blood spot samples by tandem mass spectrometry (MS/MS), using a triple quadrupole tandem mass spectrometer with [2H3] methionine and [2H5] phenylalanine as internal deuterated standards as previously described [29]. Plasma acylcarnitine and amino acid analysis were carried out at the Stoffwechsellabor (new-born screening laboratory), University Hospital Heidelberg, Germany.

The concentrations of bile acids were measured at the Analytical Unit of the Health Research Institute Hospital La Fe (Valencia, Spain) by means of a validated ultra-performance liquid chromatography/multiple reaction monitoring/mass spectrometry (UPLC-MRM-MS) method, details of which have already been published [30]. Liquid chromatography-high-resolution mass spectrometry was used to measure TMAO, betaine, and choline levels in plasma at the Institute of Clinical Chemistry, University Hospital Zurich (Zurich, Switzerland) as described previously [31].

Fecal sample collection and DNA extraction

Fecal samples were collected from participants at four different time points, i.e., baseline and weeks 12, 24, and 50 into sterile RNAlater™-filled tubes (Thermo Fisher Scientific, Germany) and stored at − 80°C at the Biobank of the National Center for Tumor Diseases, Heidelberg, Germany, until DNA extraction. Bacterial genomic DNA was extracted from fecal samples according to the manufacturer’s protocol using the AllPrep® PowerFaecal® DNA/RNA Kit (QIAGEN, Hilden, Germany). DNA from all fecal samples belonging to the same participant were extracted in the same run. In brief, approximately 200 mg of feces were added to lysis tubes containing glass beads and a lysis buffer. The tubes were then put into a Precellys® 24 Homogenizer (Bertin instruments, France) for mechanical bacterial cell lysis, which was performed twice—45 s each with an intermediate incubation period of 5 min. Genomic DNA was bound to silica membranes in spin columns, washed, and finally eluted into 1.5-ml tubes after a final incubation period of 1 min. The Qubit® 2.0 Fluorometer with broad range assay (Thermo Fisher Scientific, Germany) was used to quantify DNA concentrations (ng/μl) according to the manufacturer’s instructions. Genomic DNA was immediately stored at − 20 °C until further analyses.

Library preparation and 16S rRNA gene sequencing

In the first stage of a two-stage PCR protocol, the V4 region of the 16S rRNA bacterial gene was targeted for selective PCR amplification using V4 region-specific 16S rRNA primers 515F (GTGCCAGCMGCCGCGGTAA) and 806R (GGACTACHVGGGTWTCTAAT) [32]. For each DNA sample, a volume of 3 μl (~ 10 ng) was used as a template in a total PCR volume of 12 μl in the first-stage PCR with conditions of an initial denaturation step at 98 °C for 2 min, followed by 20 cycles of a three-step process comprising (i) denaturation at 98 °C for 10 s (ii) annealing at 65 °C for 20 s (iii) extension 72 °C for 20 s. The first-stage PCR ended with a final extension step at 72 °C for 2 min.

In the second PCR amplification, a volume of 3 μl PCR product from the first PCR step was used as a template in a ~ 12 μl total PCR volume. The PCR conditions and steps were similar to the first, but with 15 rather than 20 cycles. In order to allow for pooling/multiplexing of samples for the 16S rRNA gene sequencing after the second PCR step, each sample was amplified with a unique 12-base barcode primer (Nextera® barcodes, Illumina Inc., USA). Repeat DNA samples from the same participant were amplified in the same PCR run.

The success of amplification of each sample was assessed by visualization of bands on a 2.5% agarose gel against a 100-bp ladder. Samples that did not amplify were repeated, after which all samples were pooled. A clean-up of the pool was performed with SPRIselect® magnetic beads (Beckman Coulter Inc., USA), using 0.8 left side size selection according to the manufacturer’s instructions. The cleaned pool was again quantified using the Qubit® 2.0 Fluorometer with high sensitivity assay. The Agilent 2100 Bioanalyzer® system (Agilent Technologies, USA) was used to assess the average size distribution, the quality, and again the concentration of the pooled sample. Afterwards, a 2 × 250 paired-end sequencing of the V4 region of the 16S rRNA gene was carried out on the Illumina Miseq platform (Illumina, Inc., USA) according to the manufacturer’s instructions as previously outlined [33]. Sequencing was performed in two batches, but all samples from the same participant were sequenced in the same batch. The library preparation and 16S rRNA gene sequencing for a selection of the samples were performed twice for the technical validation of both steps (Fig. 2d). Library preparation and 16S rRNA gene sequencing were carried out at the Genomics Core Facility of the European Molecular Biology Laboratory (EMBL), Heidelberg.

Fig. 2
figure 2

Overview of the study and intervention effects on the gut microbiome composition. a Graphical overview of the HELENA Trial design showing details of all data collected over the course of the trial. In total, 150 overweight/obese adults were randomly assigned to the ICR (5 days eucaloric diet, 2 days of ~25% energy requirement), CCR (daily reduction of ~20% of caloric intake) or CTR group, 147 of them provided stool samples at baseline. Anthropometric measurements were performed for all participants. Blood samples were collected at all study timepoints for the assessment of plasma/serum concentrations of routine biomarkers. Dietary information was collected at baseline, week 12, and week 50. The amount of weight loss was significantly higher among ICR and CCR participants compared to CTR participants at all study timepoints, as previously published by Schübel et al. [22]. Differences in weight change across intervention groups at each timepoint was assessed using linear mixed models adjusted for age and sex, with participant identifiers as random effects in the model. A significant difference, i.e., p < 0.05 is indicated by (*). ICR – intermittent calorie restriction; CCR – continuous calorie restriction; CTR – controls. b The gut microbiota was predominantly composed of bacteria belonging to Clostridia followed by Bacteroidia as shown by stacked bar charts of taxonomic classes across timepoints and subdivided by intervention group. c Non-metric multidimensional scaling (NMDS) plots of Bray-Curtis dissimilarity distances based on ASV relative abundances between samples after 12 weeks of the intervention. Each point represents the microbial community of a sample, and the colour indicates the intervention group to which it belongs. NMDS plots do not show clustering by intervention groups and PERMANOVA does not indicate a statistically significant difference between groups. d Intra-individual variability (same subject compared across time points) was significantly lower compared with inter-individual variability (different subjects, same group and different subjects, different groups). Technical replicates in this study showed a very low variation in Bray-Curtis dissimilarity distances based on ASV relative abundances, attesting to a high technical reproducibility. Statistical significance as indicated above boxplots was assessed by Wilcoxon signed-rank test. e There was no significant intervention effect on overall gut microbiota alpha diversity (Shannon Index) across all intervention timepoints. All boxplots show the interquartile ranges (IQRs) as boxes, with the median as a black horizontal line and the whiskers extending up to the most extreme points within 1.5-fold IQR. Additional file 3: Source data 2

Analysis of 16S rRNA amplicon sequences

After sequencing, the raw sequence reads were demultiplexed based on their unique primer barcode identifiers using MiSeq instrument software and prepared for further downstream processing and analysis. The demultiplexed sequence reads were processed with DADA2 (Divisive Amplicon Denoising Algorithm), a tool able to identify the original sequences that generated the observed 16S sequences, known as amplicon sequence variants (ASVs) [34]. Using DADA2, raw sequence reads were quality filtered and denoised, after which paired-end reads were merged, and chimeric sequences discarded. A total of 13,767,541 reads were generated after Miseq paired-end sequencing, and 11,546,617 reads remained after filtering, denoising and chimeric sequences removal. The reads were grouped into amplicon sequence variants (ASVs) based on sequence similarity, and subsequently rarefied/subsampled to an equal depth of 5000 sequences per sample. This was performed to account for variability in sequencing depth. Only three samples were excluded after rarefaction. Standard quality control checks identified a probable case of sample swap involving 2 samples, which was resolved by relabeling based on Euclidean distance clustering of the samples (Additional File 1: Figure S1).

Alpha diversity (within-sample species diversity) at baseline and all other time points was estimated based on the Shannon index. For the assessment of gut microbiome compositional differences between the intervention groups, i.e., beta diversity, non-metric multidimensional scaling (NMDS) plots based on the Bray-Curtis dissimilarity matrices were generated on ASV relative abundances of samples for visual inspection of clustering and multivariate testing.

Taxonomic annotation at phylum, class, order, family, and genus levels was performed by matching pre-processed sequence reads against the SILVA rRNA gene database [35]. For each taxonomic rank at each specific time point(s), only bacteria clades that were present in at least 5% of the samples were retained for further analyses. The read counts were further converted to relative abundances. An unbiased assessment of the technical reproducibility using technical replicates of samples showed a low variability in the Bray-Curtis dissimilarity distances between technical replicates in this study (Fig. 2d; Technical replicates), which points to high reproducibility.

Prevotella-to-Bacteroides ratio (P/B ratio) was calculated as the summed relative abundance of the ASVs that map to Prevotella (“Prevotella”, “Prevotella_2,” etc.), divided by the cumulative relative abundance of the ASVs that map to the Bacteroides genus in Silva. In cases where the sum of Prevotella was zero, the P/B ratio was set to 1/10,000. The P/B ratio was set to 10,000 when the sum of Bacteroides relative abundances was zero. Lastly, when the sum of both Prevotella and Bacteroides relative abundances was zero, the P/B ratio was set to 1. We subsequently applied a log10 transformation to improve visualization and compared the P/B ratio across intervention groups and also assessed the association of P/B ratio with weight loss.

Statistical analysis

Data for continuous variables, including changes, are presented as means ± SEM (standard error of the mean), and as frequencies and/or percentages for categorical variables. Relative changes in body composition and in the concentrations of clinical biomarkers and circulating metabolites from baseline to each post-intervention phase, i.e., weeks 12, 24, and 50 were calculated across the initial intervention groups (ICR, CCR, and CTR). Differential changes in metabolite concentrations between the intervention groups at each time point were assessed using linear mixed effect (LME) models. The LME model included age, sex, time, intervention group, and a time-by-intervention interaction term, with participants’ identifier or subject ID set as the random effect. Differences between the intervention groups with respect to changes in the outcomes of interest were considered significant at p < 0.05 for the time-by-intervention interaction effect. When the overall time-by-intervention effect was significant, pairwise comparisons were performed, again using LME models as described.

For 16S rRNA data, bacteria relative abundances were log-transformed, after zeros had been replaced with a small number (pseudocount of 1E−04) to avoid taking the log of zero. Log-fold changes in bacteria relative abundances at each post-intervention time point were computed relative to baseline. Non-parametric Kruskal-Wallis tests were performed to assess the differences in the relative abundance of each individual gut microbial taxon between the intervention groups at baseline.

The differential effects of interventions on individual gut microbiome relative abundance and diversity at weeks 12, 24, and 50 from baseline were evaluated using LME models as already described for the intervention effects on routine clinical markers and circulating metabolites. When the difference in microbiome relative abundance across all intervention groups was significant, pairwise comparisons were performed, again with LME models with FDR corrections using the Benjamin-Hochberg FDR correction method [36]. The associations of changes in the abundances of individual bacteria which were found to be differentially abundant at specific timepoints with changes in clinical and dietary variables were also assessed using Spearman’s correlations. An association was considered to be statistically significant at an FDR-corrected p < 0.05.

In order to examine the unique contribution of baseline bacterial abundance on post-intervention abundances, i.e., the proportion of variance in post-intervention abundance explained by baseline abundances, we fitted linear regression models with age, sex, intervention group, and baseline bacterial relative abundance as independent variables, and each post-intervention abundance as dependent variable. In each model, we examined the coefficient of partial determination or partial R-squared of each independent variable and their significance in the model.

To further analyze differences in the overall gut microbiome composition between the intervention groups post-intervention, permutational analysis of variance (PERMANOVA) (permutations = 999) was performed using the adonis function in the vegan R Package [37] on Bray-Curtis dissimilarity matrices generated between samples at each timepoint.

In a further analyses step, we examined the effect of overall weight loss irrespective of the dietary regimen, i.e., ICR, CCR, or CTR on gut microbiome composition. The changes in gut microbiome abundances, alpha diversity, and overall composition across the time points were then profiled with respect to overall weight loss using LME models. The changes in metabolic parameters and circulating metabolites were also evaluated according to overall weight loss using LME models. Trends for linear associations between weight loss and changes in dependent variables, i.e., bacterial relative abundances or concentrations of circulating metabolites, were modelled with weight loss as a continuous parameter.

In order to examine associations between core families (prevalence > 50%) within the gut microbiome and body composition, clinical markers, and also dietary intake, we used LME models with the R package lmerTest [38]. LME models include both fixed and random effects and can therefore uncover associations in data which are not independent. In our case, we used LME modelling to find common associations between bacterial families and clinical markers (as fixed effects), including the participant ID as random effect to account for the non-independence of repeated measurements. The formula for the LME model was family ~ marker + (1 | Participant_ID), allowing for a different intercept for each individual. This way, inter-individual differences at baseline are accounted for, as a means of dealing with the well-established strong inter-individual variation of gut microbiome composition [39]. In a second analysis, we also included the weight loss quartile as a modifying effect by adjusting the formula to: family ~ marker + (1 + weight_loss_quartile | marker) + (1 | Participant_ID). For visualization of the results, we extracted the participant-specific intercepts and used them to adjust the log-transformed family abundances. P values from LME models were corrected with the Benjamini-Hochberg method [36].

To examine whether gut microbiota profiles at baseline were predictive of weight loss outcomes, we pooled participants across all intervention groups to identify associations between microbiome features and weight loss at later time points. At each time point, participants were grouped into weight loss quartiles based on the percentage of weight lost compared to the initial time point. We then used linear models as well as ROC analysis to find associations between the top and bottom weight lost quartile and different microbial features: we investigated (i) the baseline abundance of bacterial genera, (ii) the baseline Prevotella-to-Bacteroides ratio, (ii) baseline gut microbiota richness, and (iv) microbiome plasticity (defined as Bray-Curtis dissimilarity between the first and second time point). To compare our findings to those reported in the study by Grembi et al. [40], which investigated associations between initial microbiome composition and weight development over 1 year, we downloaded the raw reads of all technically comparable samples (paired-end sequencing, corresponding to the validation cohort in Grembi et al., n = 56, SRA project ID=PRJNA542910) and conducted the same profiling procedure described above (DADA2 plus annotation to the Silva database).

For all analyses, visual inspection of histograms and normal probability plots, as well as Kolmogorov-Smirnov tests, were performed when relevant, for the assessment of normality assumptions. All statistical analyses were two-tailed.

All statistical analyses were carried out using R statistical software [41].

Results

Baseline characteristics of participants

Overall, 147 participants (74 men and 73 women), out of the 150 enrolled, provided fecal samples at baseline (ICR, n = 49, CCR, n = 46, CTR, n = 52) (Fig. 1, see also Fig. 2), indicative of high compliance at baseline, which was also observed at subsequent timepoints. Their mean age at baseline was 50.3 ± 7.9 years, with an average BMI of 31.4 ± 3.7 kg/m2 across all groups (Table 1, see also Additional file 2: Table S1 for further baseline characteristics).

Table 1 Baseline characteristics of study participants according to intervention groupa

As expected, there were no differences in core baseline characteristics by group. When categorizing the sample by quartiles of overall weight lost during the trial, there were no significant differences at baseline either (Table 2, see also Additional file 2: Table S2 for further characteristics).

Table 2 Baseline characteristics of participants according to weight loss quartilesa

Consistent with randomization, there was no significant difference in the relative abundance of bacteria between the intervention groups (ICR vs CCR vs CTR) at baseline (Kruskal-Wallis test; FDR-corrected p-value, i.e., q value > 0.05 for individual bacteria at all taxonomic levels, for bacterial classes see Fig. 2b). Expectedly, overall gut microbiome composition was dominated mainly by two bacterial classes: Clostridia followed by Bacteroidia (Fig. 2b), which belong to the two most abundant phyla of gut microbiota, i.e., Firmicutes and Bacteroidetes, respectively [42].

Dietary intervention did not cause major changes in gut microbiome composition and diversity

Both ICR and CCR interventions, compared to CTR, resulted in significantly greater reductions in weight over the course of the intervention (Fig. 2a). In contrast, gut microbiome composition remained relatively stable and ICR and CCR did not induce major changes (Fig. 2b). Beta diversity analysis (quantifying between-sample differences) using Bray-Curtis dissimilarity on amplicon sequence variants (ASVs) did not reveal any clustering by intervention groups after 12 weeks (Fig. 2c), as is apparent from visualizations based on non-metric multidimensional scaling (NMDS) (Fig. 2c, PERMANOVA p = 0.12). Consistent with this finding, beta diversity at later time points (week 24 and week 50) was not significantly affected by the intervention either (Additional file 1: Figure S2). Likewise, we did not observe an intervention effect on the Prevotella-to-Bacteroides ratio (P/B ratio) (Additional file 1: Figure S3).

Instead, we found microbiome composition to primarily vary between individuals irrespective of intervention group. Consistent with previous observational microbiome studies [40, 43], inter-individual differences were substantially higher than intra-individual differences between samples taken from the same subject (Fig. 2d). We moreover observed a slow, but significant, decrease in intra-individual similarity in gut microbiome composition over time (Fig. 2d). However, inter-individual differences between versus within study groups were insignificant (Fig. 2d). Similar to our findings on beta diversity, the relative changes in gut microbiome alpha diversity, assessed by the Shannon index, were not significantly different between the intervention groups at any of the time points either (Fig. 2e, see also Additional file 2: Table S3).

Intervention affected intestinal Lactobacillales abundance

When assessing changes in the relative abundances of individual bacterial taxa over time, we only found the order Lactobacillales to be significantly increased in the ICR group, in contrast to marginal decreases in both the CCR and CTR (Additional file 2: Table S4) right after the intervention period (Fig. 3a). This increase in the ICR group was still detectable for the class of Bacilli, to which it belongs (Additional file 2: Table S5). At lower taxonomic ranks within Lactobacillales, the family Streptococcaceae showed the strongest change after 12 weeks (q = 0.12), with median relative changes (25th–75th percentile) of 0.62 (− 0.12,0.89) vs. 0.00 (− 0.56,0.46) vs. 0.00 (-0.51,0.42) for ICR vs. CCR vs. CTR respectively (see also Additional file 2: Table S14).

Fig. 3
figure 3

Changes in Lactobacillales abundance and processed meat intake, and associations between changes in Lactobacillales abundance and processed meat intake. a The relative abundance of the order Lactobacillales increased in the ICR group after the 12-week intervention but returned to levels similar to baseline at the end of follow-up as shown by boxplots (defined as in Fig. 3); statistical significance was assessed by Wilcoxon signed-rank test. b Participants in the ICR group significantly reduced their intake of processed meat during the 12-week intervention period, although intakes increased again at follow-up. In contrast, there was no significant change in the consumption of processed meat in either the CCR or CTR groups as shown by boxplots. c The changes in the abundance of Lactobacillales after 12 weeks were significantly associated with the changes in processed meat intake (Spearman correlation, q = 0.02). Scatter plot showing Spearman correlation between changes in the relative abundance of Lactobacillales and changes in processed meat intake. d The genus Streptococcus, which belongs to the order Lactobacillales, also significantly co-varied with the changes in processed meat intake after the 12-week intervention (q = 0.02). e Abundance changes in Lactobacillales and its member genus Streptococcus were significantly associated with changes in processed meat intake after 12 weeks. Bar length indicates q value (FDR-corrected p-value) and effect size is color-coded. Additional file 4: Source data

Given that weight loss did not significantly differ between ICR and CCR (Fig. 2a), we further analyzed whether the differential increases in Lactobacillales and Bacilli with ICR were attributable to changes in dietary composition. The only dietary item that showed a differential relative change with ICR was processed meat consumption (its decrease during intervention was likely due to higher baseline levels despite randomization, Fig. 3b). In fact, the changes in Lactobacillales relative abundance was negatively associated with changes in processed meat intake after 12 weeks (rho = − 0.19, q = 0.02, Fig. 3c. See also Fig. 3e). Within the Lactobacillales order, a slightly stronger negative association was apparent for the genus Streptococcus after 12 weeks (rho = − 0.23, q = 0.02, Fig. 3d). Details of changes in dietary intake according to intervention groups is presented in Additional file 1: Figure S4.

With respect to circulating metabolites, we did not observe significant effects of the interventions on the concentrations of plasma bile acids, most amino acids, and acylcarnitines as well as TMAO and its precursors across timepoints, consistent with the lack of significant effects of the interventions on gut microbiome composition in our study.

Weight loss was associated with changes in metabolic biomarkers but not gut microbiome composition or microbial metabolites

To assess whether the potential metabolic improvements that may accompany dietary weight loss were linked to the gut microbiome, we investigated the effect of overall weight loss, irrespective of the dietary regimen, on plasma concentrations of routine biomarkers as well as gut microbiome composition and concentrations of microbial-related metabolites in circulation. Weight loss was strongly associated with variations in body composition parameters such as visceral adipose tissue (VAT), subcutaneous adipose tissue (SAT), and liver and pancreatic fat after 12 and 50 weeks (Fig. 4a). Moreover, we found highly significant metabolic improvements associated with weight loss: established biomarkers of lipid metabolism (cholesterol, low-density lipoprotein (LDL) and leptin), liver, and kidney function (ALT and GGT) as well as cardiovascular health (thrombocytes and E-selectin as well as blood pressure) were all positively linearly associated with weight loss (Fig. 4a, d–i and Additional file 2: Table S6).

Fig. 4
figure 4

Association of weight loss with anthropometric, clinical markers, metabolic biomarkers and the gut microbiome. a Weight loss was significantly associated with decreases in body fat compartments (e.g., VAT and SAT), and routine biomarkers (e.g., cholesterol and LDL), but not with changes in gut microbial and small-molecule metabolites, with the exception of choline. Only variables significantly associated with weight loss are labelled (red dots). Measured metabolites (unlabelled due to a lack of association with weight loss) aside routine biomarkers such as those of glucose metabolism, e.g., insulin, comprised bile acids, SCFAs, TMAO, and its precursors, amino acids, and acylcarnitines, but also inflammatory cytokines (Additional file 5: Source data 4). Variables shown by dark red spots were significantly associated with weight loss also at follow-up, i.e., week 50. b The changes in weight after 12 weeks were not significantly associated with any individual bacterial taxon across all taxonomic levels, after FDR correction (all taxonomic levels included in the plot) c There was no differential association of weight loss with gut microbial alpha diversity across all timepoints. Boxplots show changes in alpha diversity over time according to quartiles of weight loss (Q1—quartile 1, Q2—quartile 2, Q3—quartile 3, and Q4—quartile 4, see Fig. 2 for definition of boxplots). Noticeably, metabolic improvements based on the changes in anthropometric and body composition parameters as well as blood concentrations of some biomarkers were greater among participants in Q4 who achieved the highest amount of weight loss after 12 weeks (See Additional file 2: Table S8). As also indicated by the volcano plot in a, there was a significant linear association between weight loss and changes in d visceral adipose tissue, e liver fat, f leptin, g gamma-glutamyl transferase (GGT), h choline, and i cholesterol after the 12-week intervention period. j Baseline abundances, rather than intervention group or weight loss accounted for the largest variation in post-intervention relative abundance of bacteria (here shown at phylum level). Bar plot showing the variations in phyla relative abundance (see colour key) after 12 weeks of intervention accounted for by age at recruitment, sex, baseline phyla abundance, and weight loss (%). Models were generated separately for weight loss and intervention, while partial R-squares of age, sex, and baseline abundances are from the model including weight loss. Additional file 5: Source data 4

However, despite the clear metabolic improvements mentioned above, there were no significant associations between overall weight loss and changes in bacterial relative abundances at any taxonomic level after the 12-week intervention (Fig. 4b). Likewise, bacterial alpha diversity was not significantly associated with weight loss at any time post-intervention (Fig. 4c, see also Table S7), nor did we observe any association with the P/B ratio (pre- or post-treatment) or changes therein (Additional file 1: Figure S3), in contrast to what was reported previously [44, 45]. The notion of metabolic improvements occurring independently of gut microbial changes was further supported by the lack of weight loss effects on circulating microbial metabolites, such as plasma SCFAs, bile acids, TMAO, and betaine after 12 and 50 weeks (Fig. 4a, see also Tables S8, S9, S10). However, choline concentrations were observed to decrease over time (Fig. 4h).

When jointly analyzing factors influencing gut microbiome composition, we found baseline microbial abundance to account for the largest proportion of variance in post-intervention abundances by far, with only minimal contributions from intervention group, weight loss, and other sources investigated (as illustrated for bacteria at the phylum level in Fig. 4j). The absence of observable shifts over time due to both intervention and weight loss was further confirmed by ordination analysis revealing predominant clustering by individuals (Additional file 1: Figure S5).

Baseline gut microbiota features predictive of weight loss

When exploring in more detail if baseline microbial features would be predictive of post-intervention weight loss, we observed the abundance of Dorea at baseline, as the only one among all genera tested, to significantly differ among participants in the lowest compared to the highest quartile of weight loss after 12 weeks (FDR-corrected p = 0.013, Fig. 5a). Low Dorea abundance at baseline was moderately predictive of weight loss after 12 weeks when evaluated as a predictive biomarker for classifying participants into the lowest versus the highest weight loss quartiles with an area under the receiver operating characteristic curve (AUROC) of 0.74 (Fig. 5b). Predictive accuracy remained similar when evaluated against weight loss at later time points (Figure S6). We further compared Dorea abundance to other biomarkers proposed earlier as predictors of weight loss, such as the Prevotella-to-Bacteroides (P/B) ratio, or gut microbiota richness [46]. However, neither the P/B ratio (AUROC: 0.56) nor microbiota richness (AUROC: 0.54) were predictive of weight loss in the HELENA Trial (Fig. 5b). We further investigated whether lower Dorea abundance at baseline would also be associated with higher weight loss in an independent dataset (Fig. 5c) collected by Grembi et al. [40] in a dietary intervention study with participant characteristics similar to ours. Also in this external data set, we observed a similar association (Pearson’s r = 0.21 compared to 0.27 in the Helena Trial), which however did not reach statistical significance (P = 0.13) (Fig. 5c). We finally also investigated whether gut microbiota plasticity, as proposed by Grembi et al. and defined as the variability in the gut microbiota composition estimated with a β-diversity metric [40], was associated with weight loss. As multiple pre-intervention microbiome samples were not available in the Helena Trial, we assessed Bray-Curtis dissimilarity-based plasticity between baseline and week 12, i.e., T0–T1 as a predictor of long-term weight loss (i.e., at week 50). In contrast to Dorea baseline abundance, microbiota plasticity was not predictive of long-term weight loss in the Helena Trial (AUROC: 0.56, Figure S6).

Fig. 5
figure 5

Association of baseline gut microbiota composition with weight loss. A higher abundance of the genus Dorea at baseline may be associated with a difficulty in losing excess weight. a Volcano plot showing the difference in the relative abundance of Dorea between participants in the highest quartile of weight loss and those in the lowest at baseline. Baseline Dorea abundances were significantly higher among participants in the lowest quartile of weight loss. b Baseline Dorea abundances rather than Prevotella-to-Bacteroides ratio or gut microbiota richness was predictive of post-intervention weight loss. AUROC showing how the aforementioned parameters were predictive of weight loss. c The relative abundance of Dorea at baseline was positively associated with weight loss, i.e., a higher abundance corresponds to a lower amount of weight loss. A scatter plot showing the association between the relative abundance of Dorea at baseline and weight loss after 12 weeks. Additional file 6: Source data 5

Gut microbial families associated with clinical biomarkers, body composition, and dietary intake

Even though we did not observe significant longitudinal changes in gut microbial abundances between intervention or weight loss groups, our study was uniquely powered to identify changes in gut microbiota associated with anthropometric, dietary, clinical, and metabolic biomarkers of obesity and early-stage metabolic syndrome when assessed irrespectively of intervention type. To leverage both the longitudinal sampling and the large number of participants, we subjected gut bacterial families to analysis with linear mixed effect (LME) models [38], which can reveal linear co-variation across timepoints that is consistent across study participants (Fig. 6a, see also Additional file 1: Figure S7 and Methods). Additionally, we also constructed LME models with the weight loss quartiles as a modifying effect, which revealed a subset of associations between bacterial families and clinical biomarkers to be of different strength across weight loss groups (Fig. 6a,b). This analysis showed Akkermansiaceae, the only member of the Verrucomicrobia phylum found in the human gut, to be inversely correlated with established biomarkers of metabolic dysfunction, namely insulin, HOMA-IR, and triglyceride concentrations (Fig. 6a,c–d). Our data also showed Christensenellaceae to be inversely associated with weight and metabolic indicators of (pre-) diabetes and cardiovascular conditions. Specifically, we found this bacterial family to be anti-correlated with BMI, VAT, liver fat, HOMA-IR, and blood concentrations of triglycerides, and cholesterol (Fig. 6a, e–g). Marinifilaceae and Rikenellaceae, which positively correlated with Christensenellaceae, were also inversely associated with VAT. In contrast to Akkermansiaceae and Christensenellaceae, we observed the abundance of Tannerellaceae to be positively associated with HOMA-IR (Fig. 6a, h). Beyond biomarkers of lipid and sugar metabolism, our analyses also uncovered interesting covariation between gut microbial families (especially Peptostreptococcaceae, but also Christensenellaceae and Marinifilaceae) and metabolic indicators of liver damage, i.e., ALT and AST, which are biomarkers of non-alcoholic fatty liver disease (NAFLD).

Fig. 6
figure 6

Association of the microbiome with anthropometric measurements, metabolic biomarkers, and dietary intake across timepoints. Several core families within the gut microbiome significantly co-vary with some of the clinical biomarkers, body composition measures, and dietary intake assessed in the trial. a Heatmap showing the strength of association calculated irrespective of weight loss quartile or outcome. Phylum affiliation for each family is indicated as a colour strip on the left. Significance for each association after FDR correction is indicated by asterisks (q < 0.1 (*), q < 0.05 (**), and q < 0.01 (***)). Asterisks colored yellow were no longer significant after adjustment for weight loss in the LME models. b Many of the associations shown in a were consistent, even after adjusting for weight loss in the LME models. The yellow points represent the associations for which significance was lost after adjustment for weight loss. c–i Scatterplots for selected significant associations based on the LME models in a. Log-transformed relative abundances of gut bacterial families were corrected for participant-specific offsets using the regression intercepts. The regression slopes for each weight loss quartile has been shown for associations that were significantly influenced by the degree of weight loss (refer to legend) (raw LME model plots without individual-specific abundance correction are shown in Figure S7). j Correlation network between bacterial families and anthropometric measurements, clinical markers and food and energy intake based on the data from a. The network includes bacterial families with at least one significant association to any of the aforementioned parameters and other families that show significant correlation across LME model coefficients (rows of the heatmap in a) with any of those families (Spearman correlation q < 0.05, edge thickness proportional to rho). Edges between families and the parameters are included if the absolute effect size (estimated by the LME) exceeded 0.075 (significant associations are indicated by stronger edges). Additional file 7: Source data 6.

We also included (self-reported) dietary intake into the LME model analysis. Both Tannerellaceae and Coriobacteriaceae were significantly inversely correlated with intake of fiber and vegetables (Fig. 6a, h–i). On the other hand, we detected a negative association between Christensenellaceae (and to a lesser extent, Akkermansiaceae) and processed meat and energy intake (Fig. 6a, g).

To investigate the similarity in association profiles between bacterial families and metabolic biomarkers, anthropometric measures, and food groups, we constructed a correlation network across bacterial families based on the LME model coefficients (Fig. 6j). Focusing on families with at least one significant association, we detected three distinct clusters, two of which were anti-correlated: One cluster, comprising both Akkermansiaceae and Christensenellaceae among other gut bacterial families, showed inverse correlations with biomarkers of metabolic dysfunction and processed meat and total energy intake. Their co-variation profiles were strongly anti-correlated to those of another small cluster, comprising Coriobacteriaceae and Tannerellacae, which were associated with lower vegetable and fiber intake, and higher levels of HOMA-IR and other biomarkers of metabolic dysfunction. The two anti-correlated clusters in this network thus revealed a broader grouping of gut microbes by their association profiles with metabolic dysfunction, suggesting that broader shifts in community composition might be associated with more pronounced changes in metabolic dysfunction or its reversal by targeted dietary interventions.

The majority of the described associations between bacteria, metabolic biomarkers, and questionnaire-derived dietary items were not differential depending on weight loss over time (Fig. 6b). Nevertheless, a few associations showed heterogeneity across weight loss quartiles in LME models, with stronger correlations between Christensenellaceae and both HOMA-IR and energy intake, and between Tannerellaceae and fiber intake, among individuals in the highest weight loss quartile (Fig. 6f–g, i, see also Additional file 1: Figure S7 and “Methods”).

Discussion

Many studies showing intriguing plasticity in gut microbiome composition in response to dietary changes [47,48,49] have raised the question of whether intestinal microbiota mediate some of the physiological effects of such interventions. In particular, understanding whether the health benefits of weight loss achieved through calorie restriction are linked to changes in the gut microbiome and its metabolites would open up possibilities for optimizing dietary interventions for intestinal microbiome modulation. Thus, we conducted the first study to longitudinally evaluate gut microbiome composition and metabolite changes following calorie restriction among 147 overweight or obese adults (non-smokers free of any gastrointestinal or chronic diseases) over 1 year in a randomized controlled dietary intervention trial. Notwithstanding significant improvements in body composition and several metabolic parameters, the type of calorie restriction, or the amount of weight lost were not accompanied by substantial and consistent shifts in gut microbiome composition or the abundance of individual bacterial taxa.

We did not observe significant changes in circulating microbial-related metabolites, such as bile acids or SCFAs either, with the exception of choline, which could be of dietary or as it has recently been shown, of microbial origin [50]. The levels of choline decreased with weight loss, consistent with previous reports associating this with improvements in insulin sensitivity [51] and a more favorable cardio-metabolic risk factor pattern in general [52]. Our data overall support the interpretation that the decrease in choline and other observed metabolic improvements upon weight loss were not dependent on the gut microbiome. By contrast, the only strong predictors of post-intervention microbial abundances in our study were their baseline abundances suggesting that changes in gut microbiome composition over time are strongly constrained by an individual’s unique microbiome configuration [53]. Taken together, these results suggest that weight loss (i.e., ~7.5–20% in the highest quartile), achieved by calorie restriction is neither associated with substantial short-term nor persistent long-term changes in gut microbiome composition or microbial metabolites in circulation and that metabolic improvements were mainly a consequence of the favorable changes in body composition (decreases in visceral fat and liver fat), rather than in the gut microbiome. This result on gut microbiome stability is in line with two recent randomized trials that have also reported only minimal effects on gut microbiome composition despite significant weight loss [54, 55].

While dietary intervention was not accompanied by significant changes in the gut microbiome over time, we found microbiome composition at baseline to be moderately predictive of weight loss suggesting that higher Dorea abundance may be associated with difficulty in losing excess weight through calorie restriction. This finding adds to previous reports that baseline gut microbiota composition may influence response to dietary weight loss interventions [13, 56]. Dorea belongs to the Clostridia class within the Firmicutes phylum and a higher abundance of this phylum has also been associated with greater energy harvest and obesity [57]. Moreover, most bacteria previously suggested to be predictive of weight loss belong to the class Clostridia [56]. However, failed validation of previously proposed baseline microbiome biomarkers in our data and the fact that the Dorea association with weight loss did not reach statistical significance in an external data set, despite a similar effect size, highlights the need for additional large-scale dietary intervention studies to resolve these discrepancies.

While our study was likely sufficiently powered (N = 147) and entailed a longitudinal component which helped to partially adjust for individual-specific microbiome composition, there are a few caveats to our findings. We cannot completely rule out potential local effects in the intestine (undetectable in plasma or fecal samples) or changes in gut microbiome function or activity that remain undetected in our study. Functional and strain-level analyses would require shotgun metagenomic or metatranscriptomic approaches, and the lack of such data is a limitation of our study. This also precluded an analysis of bacterial toxin genes, such as tcdA and tcdB encoded in the Clostridioides difficile genome, the products of which have recently been reported to affect the response to very low caloric diets among morbidly obese patients [58]. Future studies will have to resolve whether relevant alterations in gut microbial functions can occur without significant compositional changes, complementing findings of gut microbiome plasticity being often higher at the taxonomic than the functional level [39]. Finally, we acknowledge that fasting concentrations of metabolites, as used in our study, may not always best reflect metabolic dysfunction, and particularly systemic bile acid concentrations may be more informative in the postprandial state [59].

While our results are not necessarily in agreement with previous reports [13, 60], it should be noted that many previous weight loss interventions involved much more drastic reductions in calorie intake or stronger manipulations of fiber and macronutrient proportions [13, 61]. These studies mostly found abundance changes in specific gut bacterial taxa in association with specific food groups or macronutrients, particularly fiber [10, 62]. In contrast to these studies, throughout our intervention, participants consumed an “everyday diet” based on healthy eating guidelines. Fiber consumption remained relatively stable at a moderate level, while the most significant change of macronutrient proportions occurred through reduction of calories from fat (Figure S8). These modest changes in dietary composition within the limits of healthy dietary intake guidelines nevertheless resulted in substantial weight loss (~5% on average and ~10.7% in the highest quartile), similar to what other trials have achieved [10, 13], and may have contributed to the high retention rates observed throughout the intervention suggesting much better applicability for disease prevention than some previously reported extreme interventions [47, 58]. At the same time, our findings suggest that a balanced dietary composition during weight loss may help to avoid potential detrimental effects to the microbiome that have been observed as a consequence of restrictive approaches for energy reduction [13, 58].

In contrast to several recent animal experiments and first smaller human trials [63], our study does not point to differential effects of ICR vs. CCR on the microbiome, despite a differential increase in Lactobacillales with ICR after 12 weeks. Previous studies reported the genus Lactobacillus (order Lactobacillales) to increase with ICR in mice [18, 19] and among obese humans after a low-fat dietary intervention [64] or successful weight loss after Roux-en-Y Gastric Bypass [65]. However, considering that the increase in Lactobacillales in our study was associated with concurrent decreases in processed meat intake during ICR, we cannot rule out that a reduction in red meat intake rather than ICR may have led to this increase. Such an association would also be consistent with previous studies investigating dietary protein in animal models [66,67,68]. While adherence to ICR was high in our study, especially during the 12-week intervention phase [22], we acknowledge that the 5:2 diet is a mild form of ICR and that further studies on other types of ICR in relation to the microbiome are needed. However, recent meta-analyses of randomized studies do not suggest stronger effects of ICR vs. CCR diets on established metabolic biomarkers, which may speak against differential effects of ICR on the microbiome [69,70,71].

In a pooled analysis of all study groups, we found that many gut microbial families co-varied significantly and consistently across participants with host metabolic parameters, forming a network with emerging modules of intercorrelated taxa. In particular, our findings from robust LME models leveraging repeated measurements over time on significant associations of Christensenellaceae, Akkermansiaceae, Marinifilaceae, and Tannerellaceae with HOMA-IR (and triglycerides for the former two) underline the potential importance of these bacteria for host metabolic status. Consistent with our observation, previous cross-sectional studies have associated a healthy lipid profile and normal glucose homeostasis with increased abundance of the genus Akkermansia [72, 73]. Intriguingly, oral supplementation of Akkermansia muciniphila was shown to alleviate insulinemia and to lower the plasma concentrations of cholesterol, inflammatory biomarkers, and markers of liver dysfunction among overweight/obese insulin-resistant adults [74]. Also the Christensenellaceae family was previously associated with an overall healthy metabolic status [75, 76] and had been shown to be more abundant among individuals with low BMI in cross-sectional twin studies [77]. Our results complement the earlier studies, as they show that Christensenellaceae fluctuates over time around individual-specific levels and inversely co-varies with anthropometric and metabolic biomarkers of obesity and diabetes in a consistent manner across individuals. The significant inverse association of Peptostreptococcaceae with both ALT and AST, and to a lesser extent with GGT and CRP, is interesting in light of the discrepant reports on the putative role of this intestinal bacterial family in liver disease (sometimes found to increase, sometimes to decrease, with liver dysfunction) [78, 79]. Our results rather support a potentially disease-protective role. Finally, Marinifilaceae have recently been shown to be enriched in adipose tissue samples of obese people without diabetes compared to obese people with type 2 diabetes [80], which is somewhat in line with the inverse association between Marinifilaceae and VAT in our study.

As opposed to the above-described network of bacteria related to a favorable pattern of metabolic biomarkers, we also observed an anti-correlated profile characterized by co-variation of Tannerellaceae and Coriobacteriaceae. To our knowledge, Tannerellaceae, which were negatively associated with vegetable and fiber intake but positively with HOMA-IR, have not previously been associated with obesity or diabetes in human studies. In mice, however, the abundance of Tannerellaceae was found to increase upon a high-fat diet intervention with concomitant elevations in plasma glucose and lipid concentrations [81]. Tannerellaceae were positively correlated with Coriobacteriaceae that also showed inverse associations with vegetable and fiber intake. Collinsella, the only intestinal genus in the latter bacterial family, has been reported to be enriched among individuals with low fiber intake and to be inversely associated with a Mediterranean dietary pattern [82, 83].

We found a few cases of heterogeneity across quartiles of weight loss in the associations between bacteria, metabolic biomarkers, and food intakes in integrated LME analyses of repeated measurements. However, while our analysis of intervention outcomes suggests that they cannot be easily manipulated by calorie restriction alone, their associations with food components, such as vegetables or processed meat, highlight the promise of future dietary intervention studies; based on our post hoc observations, such studies could attempt to specifically modulate these bacterial families while monitoring concomitant metabolic changes. Together with randomized controlled trials of bacterial supplementation [74], these could help to untangle causality and to achieve weight loss and its associated metabolic benefits not only through direct effects on the host, but also through synergistic modulation of the gut microbiota and its metabolic products.

Conclusions

In conclusion, our findings suggest that moderate ICR or CCR interventions as well as an overall moderate weight loss induced by calorie restriction (irrespective of which form) may not be associated with significant changes in the gut microbiome of overweight and obese adults, notwithstanding observed metabolic improvements. Nevertheless, baseline microbiome composition, in particular Dorea, was moderately predictive of intervention outcome in our study, but further work is needed to independently validate this in comparison to other proposed predictors of weight loss. Despite the lack of consistent intervention effects on microbiome composition, LME model analysis, which leveraged data across all time points and participants, while also accounting for weight loss, revealed several significant and plausible associations of gut bacteria with anthropometric measurements, body composition indices, and biomarkers of metabolic health.

Availability of data and materials

The 16S rRNA amplicon sequencing data from this study are available from the European Nucleotide Archive (ENA) database: accession number PRJEB40697 (https://www.ebi.ac.uk/ena/browser/view/PRJEB40697) (see also Additional file 2: Table S21). Changes in the microbiome relative abundance at the phylum, family, and genus taxonomic levels according to intervention groups are also available on Zenodo at (https://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.4287277) [84]. Information on the concentrations of bile acids, acylcarnitines, amino acids, and choline derivatives at baseline, as well as the changes in the concentrations of these metabolites across timepoints according to both intervention group and weight loss quartiles, is publicly available on Zenodo (https://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.4287277) [84]. Source data for all figures shown in this paper have been provided. All other data is available from the corresponding authors on request, on the basis of a data transfer agreement, which is required according to the informed consent statement used for the study.

The 16S rRNA amplicon sequencing data from this study are available from the European Nucleotide Archive (ENA) database: accession number PRJEB40697 (https://www.ebi.ac.uk/ena/browser/view/PRJEB40697) [85] (see also Additional file 2: Table S21). All other data is available from the corresponding authors on request, on the basis of a data transfer agreement, which is required according to the informed consent statement used for the study.

Abbreviations

ALT:

Alanine aminotransferase

AST:

Aspartate aminotransferase

AUROC:

Area under the receiver operating characteristic curve

BMI:

Body mass index

BP:

Blood pressure

CCR:

Continuous calorie restriction

CRP:

C-reactive protein

CTR:

Control group

CVD:

Cardiovascular disease

DNA:

Deoxyribonucleic acid

GGT:

Gamma-glutamyl transferase

HDL:

High-density lipoprotein

ICR:

Intermittent calorie restriction

IFN-γ:

Interferon gamma

IGF-1:

Insulin-like growth factor 1

IL-6:

Interleukin 6

IL-8:

Interleukin 8

LC:

Liquid chromatography

LME:

Linear mixed effect

MS:

Mass spectrometry

NMDS:

Non-metric multidimensional scaling

PERMANOVA:

Permutational analysis of variance

RCT:

Randomized controlled trial

rRNA:

Ribosomal ribonucleic acid

SAT:

Subcutaneous adipose tissue

TNF-α:

Tumor necrosis factor-α

VAT:

Visceral adipose tissue

References

  1. Calle EE, Rodriguez C, Walker-Thurmond K, Thun MJ. Overweight, obesity, and mortality from cancer in a prospectively studied cohort of U.S. adults. N Engl J Med. 2003;348(17):1625–38.

    PubMed  Google Scholar 

  2. Pi-Sunyer FX. The obesity epidemic: pathophysiology and consequences of obesity. Obes Res. 2002;10 Suppl 2:97s–104s.

    PubMed  Google Scholar 

  3. Smith SR, Lovejoy JC, Greenway F, Ryan D, de Jonge L, de la Bretonne J, et al. Contributions of total body fat, abdominal subcutaneous adipose tissue compartments, and visceral adipose tissue to the metabolic complications of obesity. Metabolism. 2001;50(4):425–35.

    CAS  PubMed  Google Scholar 

  4. Shoelson SE, Herrero L, Naaz A. Obesity, inflammation, and insulin resistance. Gastroenterology. 2007;132(6):2169–80.

    CAS  PubMed  Google Scholar 

  5. Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006;444(7122):1027–31.

    PubMed  Google Scholar 

  6. Chambers ES, Preston T, Frost G, Morrison DJ. Role of gut microbiota-generated short-chain fatty acids in metabolic and cardiovascular health. Curr Nutr Rep. 2018;7(4):198–206.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Yoshimoto S, Loo TM, Atarashi K, Kanda H, Sato S, Oyadomari S, et al. Obesity-induced gut microbial metabolite promotes liver cancer through senescence secretome. Nature. 2013;499(7456):97–101.

    CAS  PubMed  Google Scholar 

  8. Wang Z, Klipfell E, Bennett BJ, Koeth R, Levison BS, Dugar B, et al. Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature. 2011;472(7341):57–63.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Liu R, Hong J, Xu X, Feng Q, Zhang D, Gu Y, et al. Gut microbiome and serum metabolome alterations in obesity and after weight-loss intervention. Nat Med. 2017;23(7):859–68.

    CAS  PubMed  Google Scholar 

  10. Cotillard A, Kennedy SP, Kong LC, Prifti E, Pons N, Le Chatelier E, et al. Dietary intervention impact on gut microbial gene richness. Nature. 2013;500(7464):585–8.

    CAS  PubMed  Google Scholar 

  11. Frost F, Storck LJ, Kacprowski T, Gärtner S, Rühlemann M, Bang C, et al. A structured weight loss program increases gut microbiota phylogenetic diversity and reduces levels of Collinsella in obese type 2 diabetics: A pilot study. Plos One. 2019;14(7):e0219489.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Heinsen FA, Fangmann D, Müller N, Schulte DM, Rühlemann MC, Türk K, et al. Beneficial effects of a dietary weight loss intervention on human gut microbiome diversity and metabolism are not sustained during weight maintenance. Obes Facts. 2016;9(6):379–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Seganfredo FB, Blume CA, Moehlecke M, Giongo A, Casagrande DS, Spolidoro JVN, et al. Weight-loss interventions and gut microbiota changes in overweight and obese patients: a systematic review. Obes Rev. 2017;18(8):832–51.

    CAS  PubMed  Google Scholar 

  14. Leeming ER, Johnson AJ, Spector TD, Le Roy CI. Effect of diet on the gut microbiota: rethinking intervention duration. Nutrients. 2019;11(12):2862.

  15. Gerber GK. The dynamic microbiome. FEBS Lett. 2014;588(22):4131–9.

    CAS  PubMed  Google Scholar 

  16. David LA, Materna AC, Friedman J, Campos-Baptista MI, Blackburn MC, Perrotta A, et al. Host lifestyle affects human microbiota on daily timescales. Genome Biolo. 2014;15(7):R89.

    Google Scholar 

  17. Cignarella F, Cantoni C, Ghezzi L, Salter A, Dorsett Y, Chen L, et al. Intermittent fasting confers protection in CNS autoimmunity by altering the gut microbiota. Cell Metab. 2018;27(6):1222–35.e6.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Liu Z, Dai X, Zhang H, Shi R, Hui Y, Jin X, et al. Gut microbiota mediates intermittent-fasting alleviation of diabetes-induced cognitive impairment. Nat Commun. 2020;11(1):855.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Beli E, Yan Y, Moldovan L, Vieira CP, Gao R, Duan Y, et al. Restructuring of the gut microbiome by intermittent fasting prevents retinopathy and prolongs survival in db/db Mice. Diabetes. 2018;67(9):1867–79.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Li G, Xie C, Lu S, Nichols RG, Tian Y, Li L, et al. Intermittent fasting promotes white adipose browning and decreases obesity by shaping the gut microbiota. Cell Metab. 2017;26(5):801.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Muñoz-Hernández L, Márquez-López Z, Mehta R, Aguilar-Salinas CA. Intermittent fasting as part of the management for T2DM: from animal models to human clinical studies. Curr Diab Rep. 2020;20(4):13.

    PubMed  Google Scholar 

  22. Schübel R, Nattenmuller J, Sookthai D, Nonnenmacher T, Graf ME, Riedl L, et al. Effects of intermittent and continuous calorie restriction on body weight and metabolism over 50 wk: a randomized controlled trial. Am J Clin Nutr. 2018;108(5):933–45.

    PubMed  PubMed Central  Google Scholar 

  23. Schubel R, Graf ME, Nattenmuller J, Nabers D, Sookthai D, Gruner LF, et al. The effects of intermittent calorie restriction on metabolic health: rationale and study design of the HELENA Trial. Contemp Clin Trials. 2016;51:28–33.

    PubMed  Google Scholar 

  24. Oberritter H, Schäbethal K, von Ruesten A, Boeing H. The DGE Nutrition Circle – presentation and basis of the food-related recommendations from the German Nutrition Society (DGE). Ernährungs Umschau. 2013;2:24–9.

    Google Scholar 

  25. World Medical Association. World Medical Association Declaration of Helsinki: ethical principles for medical research involving human subjects. JAMA. 2013;310(20):2191–4.

    Google Scholar 

  26. Chan JC, Kioh DY, Yap GC, Lee BW, Chan EC. A novel LCMSMS method for quantitative measurement of short-chain fatty acids in human stool derivatized with (12) C- and (13) C-labelled aniline. J Pharm Biomed Anal. 2017;138:43–53.

    CAS  PubMed  Google Scholar 

  27. Zeng M, Cao H. Fast quantification of short chain fatty acids and ketone bodies by liquid chromatography-tandem mass spectrometry after facile derivatization coupled with liquid-liquid extraction. J Chromatog B Analyt Technol Biomed Life Sci. 2018;1083:137–45.

    CAS  Google Scholar 

  28. Sowah SA, Hirche F, Milanese A, Johnson TS, Grafetstatter M, Schubel R, et al. Changes in plasma short-chain fatty acid levels after dietary weight loss among overweight and obese adults over 50 weeks. Nutrients. 2020;12(2):452.

  29. Lindner M, Abdoh G, Fang-Hoffmann J, Shabeck N, Al-Sayrafi M, Al-Janahi M, et al. Implementation of extended neonatal screening and a metabolic unit in the State of Qatar: developing and optimizing strategies in cooperation with the Neonatal Screening Center in Heidelberg. J Inherit Metab Dis. 2007;30(4):522–9.

    CAS  PubMed  Google Scholar 

  30. Garcia-Canaveras JC, Donato MT, Castell JV, Lahoz A. Targeted profiling of circulating and hepatic bile acids in human, mouse, and rat using a UPLC-MRM-MS-validated method. J Lipid Res. 2012;53(10):2231–41.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Mueller DM, Allenspach M, Othman A, Saely CH, Muendlein A, Vonbank A, et al. Plasma levels of trimethylamine-N-oxide are confounded by impaired kidney function and poor metabolic control. Atherosclerosis. 2015;243(2):638–44.

    CAS  PubMed  Google Scholar 

  32. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci USA. 2011;108 Suppl 1(Suppl 1):4516–22.

    PubMed  Google Scholar 

  33. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6(8):1621–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(Database issue):D590–6.

    CAS  PubMed  Google Scholar 

  36. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B (Methodological). 1995;57(1):289–300.

    Google Scholar 

  37. Jari O, Guillaume B, Michael Friendly, Roeland Kindt, Pierre Legendre, Dan McGlinn, et al. Vegan: Community Ecology Package. R package version 2.5-6. 2019.

    Google Scholar 

  38. Kunzetsova A, Brockhoff P, Christensen R. lmerTest package: tests in linear mixed effect models. J Stat Softw. 2017;82:1–26.

    Google Scholar 

  39. Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486(7402):207–14.

    Google Scholar 

  40. Grembi JA, Nguyen LH, Haggerty TD, Gardner CD, Holmes SP, Parsonnet J. Gut microbiota plasticity is correlated with sustained weight loss on a low-carb or low-fat dietary intervention. Sci Rep. 2020;10(1):1405.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2019.

    Google Scholar 

  42. Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, et al. Diversity of the human intestinal microbial flora. Science. 2005;308(5728):1635–8.

    PubMed  PubMed Central  Google Scholar 

  43. Caporaso JG, Lauber CL, Costello EK, Berg-Lyons D, Gonzalez A, Stombaugh J, et al. Moving pictures of the human microbiome. Genome Biol. 2011;12(5):R50.

    PubMed  PubMed Central  Google Scholar 

  44. Hjorth MF, Christensen L, Larsen TM, Roager HM, Krych L, Kot W, et al. Pretreatment Prevotella-to-Bacteroides ratio and salivary amylase gene copy number as prognostic markers for dietary weight loss. Am J Clin Nutr. 2020;111(5):1079–86.

    PubMed  Google Scholar 

  45. Hjorth MF, Blaedel T, Bendtsen LQ, Lorenzen JK, Holm JB, Kiilerich P, et al. Prevotella-to-Bacteroides ratio predicts body weight and fat loss success on 24-week diets varying in macronutrient composition and dietary fiber: results from a post-hoc analysis. Int J Obes. 2019;43(1):149-57.

  46. Biesiekierski JR, Jalanka J, Staudacher HM. Can gut microbiota composition predict response to dietary treatments? Nutrients. 2019;11(5):1134.

  47. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505(7484):559–63.

    CAS  PubMed  Google Scholar 

  48. Wu GD, Chen J, Hoffmann C, Bittinger K, Chen YY, Keilbaugh SA, et al. Linking long-term dietary patterns with gut microbial enterotypes. Science. 2011;334(6052):105–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  49. Walker AW, Ince J, Duncan SH, Webster LM, Holtrop G, Ze X, et al. Dominant and diet-responsive groups of bacteria within the human colonic microbiota. ISME J. 2011;5(2):220–30.

    CAS  PubMed  Google Scholar 

  50. Chittim CL, Martínez Del Campo A, Balskus EP. Gut bacterial phospholipase Ds support disease-associated metabolism by generating choline. Nat Microbiol. 2019;4(1):155–63.

    CAS  PubMed  Google Scholar 

  51. Heianza Y, Sun D, Li X, DiDonato JA, Bray GA, Sacks FM, et al. Gut microbiota metabolites, amino acid metabolites and improvements in insulin sensitivity and glucose metabolism: the POUNDS Lost trial. Gut. 2019;68(2):263–70.

    CAS  PubMed  Google Scholar 

  52. Lent-Schochet D, McLaughlin M, Ramakrishnan N, Jialal I. Exploratory metabolomics of metabolic syndrome: a status report. World J Diabetes. 2019;10(1):23–36.

    PubMed  PubMed Central  Google Scholar 

  53. Costea PI, Hildebrand F, Arumugam M, Bäckhed F, Blaser MJ, Bushman FD, et al. Enterotypes in the landscape of gut microbial community composition. Nat Microbiol. 2018;3(1):8–16.

    CAS  PubMed  Google Scholar 

  54. Roager HM, Vogt JK, Kristensen M, Hansen LBS, Ibrugger S, Maerkedahl RB, et al. Whole grain-rich diet reduces body weight and systemic low-grade inflammation without inducing major changes of the gut microbiome: a randomised cross-over trial. Gut. 2019;68(1):83-93.

  55. Fragiadakis GK, Wastyk HC, Robinson JL, Sonnenburg ED, Sonnenburg JL, Gardner CD. Long-term dietary intervention reveals resilience of the gut microbiota despite changes in diet and weight. Am J Clin Nutr. 2020;111(6):1127–36.

    PubMed  PubMed Central  Google Scholar 

  56. Korpela K, Flint HJ, Johnstone AM, Lappi J, Poutanen K, Dewulf E, et al. Gut microbiota signatures predict host and microbiota responses to dietary interventions in obese individuals. Plos One. 2014;9(6):e90702.

    PubMed  PubMed Central  Google Scholar 

  57. Ley RE, Backhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI. Obesity alters gut microbial ecology. Proc Natl Acad Sci USA. 2005;102(31):11070–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. von Schwartzenberg RJ, Bisanz JE, Lyalina S, Spanogiannopoulos P, Ang QY, Cai J, et al. Caloric restriction disrupts the microbiota and colonization resistance. Nature. 2021;595(7866):272–7.

    Google Scholar 

  59. Kuipers F, Bloks VW, Groen AK. Beyond intestinal soap--bile acids in metabolic control. Nat Rev Endocrinol. 2014;10(8):488–98.

    CAS  PubMed  Google Scholar 

  60. Nadal I, Santacruz A, Marcos A, Warnberg J, Garagorri JM, Moreno LA, et al. Shifts in clostridia, bacteroides and immunoglobulin-coating fecal bacteria associated with weight loss in obese adolescents. Int J Obes. 2009;33(7):758-67.

    CAS  Google Scholar 

  61. Damms-Machado A, Mitra S, Schollenberger AE, Kramer KM, Meile T, Konigsrainer A, et al. Effects of surgical and dietary weight loss therapy for obesity on gut microbiota composition and nutrient absorption. BioMed Res Int. 2015;2015:806248.

    PubMed  PubMed Central  Google Scholar 

  62. Duncan SH, Belenguer A, Holtrop G, Johnstone AM, Flint HJ, Lobley GE. Reduced dietary intake of carbohydrates by obese subjects results in decreased concentrations of butyrate and butyrate-producing bacteria in feces. Appl Environ Microbiol. 2007;73(4):1073–8.

    CAS  PubMed  Google Scholar 

  63. Daas MC, de Roos NM. Intermittent fasting contributes to aligned circadian rhythms through interactions with the gut microbiome. Benef Microbes. 2021;12(2):147-161.

  64. Remely M, Tesar I, Hippe B, Gnauer S, Rust P, Haslberger AG. Gut microbiota composition correlates with changes in body fat content due to weight loss. Benef Microbes. 2015;6(4):431–9.

    CAS  PubMed  Google Scholar 

  65. Fouladi F, Brooks AE, Fodor AA, Carroll IM, Bulik-Sullivan EC, Tsilimigras MCB, et al. The role of the gut microbiota in sustained weight loss following Roux-en-Y gastric bypass surgery. Obes Surg. 2019;29(4):1259–67.

    PubMed  Google Scholar 

  66. Schutkowski A, König B, Kluge H, Hirche F, Henze A, Schwerdtle T, et al. Metabolic footprint and intestinal microbial changes in response to dietary proteins in a pig model. J Nutr Biochem. 2019;67:149–60.

    CAS  PubMed  Google Scholar 

  67. Zhu Y, Lin X, Zhao F, Shi X, Li H, Li Y, et al. Meat, dairy and plant proteins alter bacterial composition of rat gut bacteria. Sci Rep. 2015;5(1):15220.

    CAS  PubMed  PubMed Central  Google Scholar 

  68. Zhu Y, Lin X, Li H, Li Y, Shi X, Zhao F, et al. Intake of meat proteins substantially increased the relative abundance of Genus Lactobacillus in rat feces. Plos one. 2016;11(4):e0152678.

    PubMed  PubMed Central  Google Scholar 

  69. Allaf M, Elghazaly H, Mohamed OG, Fareen MFK, Zaman S, Salmasi AM, et al. Intermittent fasting for the prevention of cardiovascular disease. Cochrane Database Syst Rev. 2021;1:Cd013496.

    PubMed  Google Scholar 

  70. Schwingshackl L, Zähringer J, Nitschke K, Torbahn G, Lohner S, Kühn T, et al. Impact of intermittent energy restriction on anthropometric outcomes and intermediate disease markers in patients with overweight and obesity: systematic review and meta-analyses. Crit Rev Food Sci Nutr. 2021;61(8):1293–1304.

  71. Vitale R, Kim Y. The effects of intermittent fasting on glycemic control and body composition in adults with obesity and type 2 diabetes: a systematic review. Metab Syndr RelatDisord. 2020;18(10):450–61.

    CAS  Google Scholar 

  72. Schneeberger M, Everard A, Gómez-Valadés AG, Matamoros S, Ramírez S, Delzenne NM, et al. Akkermansia muciniphila inversely correlates with the onset of inflammation, altered adipose tissue metabolism and metabolic disorders during obesity in mice. Sci Rep. 2015;5:16643.

    CAS  PubMed  PubMed Central  Google Scholar 

  73. Zhang X, Shen D, Fang Z, Jie Z, Qiu X, Zhang C, et al. Human gut microbiota changes reveal the progression of glucose intolerance. Plos One. 2013;8(8):e71108.

    CAS  PubMed  PubMed Central  Google Scholar 

  74. Depommier C, Everard A, Druart C, Plovier H, Van Hul M, Vieira-Silva S, et al. Supplementation with Akkermansia muciniphila in overweight and obese human volunteers: a proof-of-concept exploratory study. nat Med. 2019;25(7):1096–103.

    CAS  PubMed  PubMed Central  Google Scholar 

  75. Waters JL, Ley RE. The human gut bacteria Christensenellaceae are widespread, heritable, and associated with health. BMC Biol. 2019;17(1):83.

    PubMed  PubMed Central  Google Scholar 

  76. Tavella T, Rampelli S, Guidarelli G, Bazzocchi A, Gasperini C, Pujos-Guillot E, et al. Elevated gut microbiome abundance of Christensenellaceae, Porphyromonadaceae and Rikenellaceae is associated with reduced visceral adipose tissue and healthier metabolic profile in Italian elderly. Gut Microbes. 2021;13(1):1–19.

    PubMed  Google Scholar 

  77. Goodrich JK, Waters JL, Poole AC, Sutter JL, Koren O, Blekhman R, et al. Human genetics shape the gut microbiome. Cell. 2014;159(4):789–99.

    CAS  PubMed  PubMed Central  Google Scholar 

  78. Jiang W, Wu N, Wang X, Chi Y, Zhang Y, Qiu X, et al. Dysbiosis gut microbiota associated with inflammation and impaired mucosal immune function in intestine of humans with non-alcoholic fatty liver disease. Sci Rep. 2015;5:8096.

    CAS  PubMed  PubMed Central  Google Scholar 

  79. Caussy C, Tripathi A, Humphrey G, Bassirian S, Singh S, Faulkner C, et al. A gut microbiome signature for cirrhosis due to nonalcoholic fatty liver disease. Nat Commun. 2019;10(1):1406.

    PubMed  PubMed Central  Google Scholar 

  80. Anhê FF, Jensen BAH, Varin TV, Servant F, Van Blerk S, Richard D, et al. Type 2 diabetes influences bacterial tissue compartmentalisation in human obesity. Nat Metab. 2020;2(3):233–42.

    PubMed  Google Scholar 

  81. Wang J-H, Shin NR, Lim S-K, Im U, Song E-J, Nam Y-D, et al. Diet control more intensively disturbs gut microbiota than genetic background in wild type and ob/ob mice. Front Microbiol. 2019;10:1292.

  82. Gomez-Arango LF, Barrett HL, Wilkinson SA, Callaway LK, McIntyre HD, Morrison M, et al. Low dietary fiber intake increases Collinsella abundance in the gut microbiota of overweight and obese pregnant women. Gut Microbes. 2018;9(3):189–201.

    PubMed  PubMed Central  Google Scholar 

  83. Ghosh TS, Rampelli S, Jeffery IB, Santoro A, Neto M, Capri M, et al. Mediterranean diet intervention alters the gut microbiome in older people reducing frailty and improving health status: the NU-AGE 1-year dietary intervention across five European countries. Gut. 2020;69(7):1218–28.

    CAS  PubMed  Google Scholar 

  84. Sowah SA, Milanese A, Schübel R, Wirbel J, Kartal E, Johnson TS, et al. Dataset - Calorie restriction improves metabolic state independently of gut microbiome composition Figshare; 2020. https://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.4287277.

    Book  Google Scholar 

  85. Sowah SA, Milanese A, Schübel R, Wirbel J, Kartal E, Johnson TS, et al. Calorie restriction improves metabolic state independently of gut microbiome composition: 16S rRNA amplicon sequencing data; ENA Accession Number PRJEB40697, European Nucleotide Archive (ENA) database. 2020. https://www.ebi.ac.uk/ena/browser/view/PRJEB40697.

    Google Scholar 

Download references

Acknowledgements

The authors thank the study participants of the HELENA Trial. We would like to acknowledge the support of Anja Telzerow and Vladimir Benes from the EMBL Genomics Core Facility during the library preparation and 16S rRNA sequencing. The authors would like to thank Marita Wenzel, Biobank of the National Center for Tumor Diseases (NCT) Heidelberg, Heidelberg for her assistance with sample storage and retrieval, and Muhabbet Celik, Bettina Ehret and Christine Niesik from the Division of Cancer Epidemiology, DKFZ for their support in the lab.

Funding

Open Access funding enabled and organized by Projekt DEAL. The HELENA Trial was sponsored by the Helmholtz Association of German Research Centers (Cross Program Topic: Ageing and Metabolic Dysfunction, AMPro). Bile acid measurements in the present study were funded by Helmholtz Association of German Research Centers via the project Aging and Metabolic Programming (AMPro; DKFZ). MRI examinations were performed at and financed by the Department of Diagnostic and Interventional Radiology, University Hospital Heidelberg. MRI financing was also supported by the Stiftung zur Förderung der Erforschung der Zivilisationserkrankungen, Baden-Baden, Germany. Data analysis was supported by the Federal Ministry of Education and Research (BMBF; grant no. 031L0181A to G.Z.). C.M.U was supported by the Huntsman Cancer Foundation, Salt Lake City, UT, and by NIH grant U01 CA 206110. S.A.S is supported by the German Academic Exchange Service (DAAD) via the Graduate School Scholarship Programme at the Helmholtz International Graduate School for Cancer Research in Heidelberg, Germany. The funders had no role in the design and conduct of this study.

Author information

Authors and Affiliations

Authors

Contributions

S.A.S, A.M., T. K, and G. Z conceived and designed the study. R. S, R.K, T. K, C.M.U, and G. Z acquired funding. R. S, M. G, L. G, T. N, J. N, T.S.J, R. K, and T. K conducted the study. S.A.S performed the microbiome experiments with support from E.K. A. M performed the bioinformatics analysis with supervision from G.Z. S.A.S, A. M, and J.W. performed the statistical analysis with supervision from G. Z and T.K. S.A.S wrote the paper with contributions from A.M. T. K and G. Z revised the first draft of the manuscript, provided critical comments, and supervised all aspects of the present study. All authors provided critical comments on the manuscript, and all authors read and approved the final manuscript. S.A.S and A. M contributed equally to this paper, T. K and G. Z equally supervised this work and are joint last authors. All author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Tilman Kühn or Georg Zeller.

Ethics declarations

Ethics approval and consent to participate

The ethics committee of the Heidelberg University Hospital (Heidelberg, Germany , Reference: vote S299-2014) and the EMBL Bioethics Advisory Committee approved the trial (Reference: 2020-022), which was designed and conducted in agreement with the principles of the World Medical Association’s Declaration of Helsinki. All participants signed informed consent forms before enrolment.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Tilman Kühn and Georg Zeller are joint last authors.

Supplementary Information

Additional file 1: Figure S1.

Identification of possibly swapped samples. In the figure, samples are displayed as subject id and time point. The affected samples are 700119_3 in the top panel, 700053_T3 in the middle panel, 700085_T0 and 700085_T2 in the bottom panel. Figure S2. NMDS plots of Bray-Curtis dissimilarity matrices did not reveal distinct clustering patterns according to intervention groups at weeks (a) 24 and (b) 50. Figure S3. Prevotella and Bacteroides relative abundance of all 16S amplicon samples (a), and histogram of the distribution of the Prevotella-to-Bacteroides ratio (b). Bar Plot showing the distribution of the P/B ratio per intervention group (c), time point (d) and combination of intervention group and time point (e). Correlation between weight loss and P/B ratio pre-intervention (f), post-intervention (g) and changes in P/B ratio (post- minus pre-intervention) (h). There is a nominal association in panel c which disappears with multiple testing correction. Figure S4. Body weight and energy, fiber, macronutrients, fruits, vegetables and processed meat consumption across timepoints according to intervention groups. Figure S5. Non-metric multidimensional scaling (NMDS) plots of Bray-Curtis dissimilarity distances between samples across all timepoints according to (a) intervention groups and (b) weight loss quartiles. Figure S6. Association of baseline Dorea abundance and gut microbiota plasticity (baseline to week 12) with long-term weight loss i.e. week 24 and week 50. Figure S7. LME plots for all significant associations between bacterial families and anthropometric measures (see Figure 6a in the main text). Figure S8. Body weight and intakes energy, fiber, macronutrients, fruits, vegetables and processed meat across timepoints according to weight loss quartiles.

Additional file 2: Table S1.

Further baseline characteristics of participants according to intervention group. Table S2. Further baseline characteristics of participants according to weight loss quartiles. Table S3. Diversity indices across timepoints according to intervention groups. Table S4. Changes in bacteria (order) abundance across time points according to intervention group. Table S5. Changes in bacteria (class) abundance across time points according to intervention group. Table S6. Changes in anthropometric measurements, body composition, clinical markers and blood-based biomarkers by weight loss quartiles after 24 and 50 weeks. Table S7. Diversity indices across timepoints according to weight loss quartiles. Table S8. Intervention effects on plasma bile acids concentrations (nmol/L) after 12 weeks according to weight loss quartiles. Table S9. Changes in the concentrations of acylcarnitines and amino acids according to weight loss quartiles after 12, 24 and 50 weeks. Table S10. Intervention effects on plasma TMAO, betaine and choline concentrations (μmol/L) across timepoints according to weight lossquartiles. Table S11. Comparison of bacteria (phylum) relative abundances across timepoints according to intervention groups. Table S12. Comparison of bacteria (class) relative abundances across timepoints according to intervention groups. Table S13. Comparison of bacteria (order) relative abundances across timepoints according to intervention groups. Table S14. Comparison of bacteria (family) relative abundances across timepoints according to intervention groups. Table S15. Comparison of bacteria (genus) relative abundances across timepoints according to intervention groups. Table S16. Comparison of within-group bacteria (phylum) relative abundances at each timepoint relative to baseline. Table S17. Comparison of within-group bacteria (class) relative abundances at each timepoint relative to baseline. Table S18. Comparison of within-group bacteria (order) relative abundances at each timepoint relative to baseline. Table S19. Comparison of within-group bacteria (family) relative abundances at each timepoint relative to baseline. Table S20. Comparison of within-group bacteria (genus) relative abundances at each timepoint relative to baseline. Table S21. Study, sample, experiment and run accession numbers of all 16S data generated in the HELENA Trial.

Additional file 3. Source data for Figure 2.

Additional file 4. Source data for Figure 3.

Additional file 5. Source data for Figure 4.

Additional file 6. Source data for Figure 5.

Additional file 7. Source data for Figure 6.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sowah, S.A., Milanese, A., Schübel, R. et al. Calorie restriction improves metabolic state independently of gut microbiome composition: a randomized dietary intervention trial. Genome Med 14, 30 (2022). https://0-doi-org.brum.beds.ac.uk/10.1186/s13073-022-01030-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13073-022-01030-0

Keywords