What this is
- The study investigates the impact of the Green-Mediterranean (Green-MED) diet on changes and .
- Participants with abdominal obesity/dyslipidemia were randomized into three groups: standard dietary guidelines, Mediterranean diet (MED), and Green-MED.
- The Green-MED diet emphasizes increased plant-based foods and reduced meat intake, supplemented with walnuts and green tea.
Essence
- The Green-MED diet significantly alters composition and improves markers, with the microbiome mediating some of these benefits.
Key takeaways
- The Green-MED diet led to greater changes in composition compared to standard dietary guidelines and the MED diet.
- Specific microbial changes, such as increases in Prevotella and decreases in Bifidobacterium, were linked to the Green-MED diet.
- Improvements in body weight and cardiometabolic markers were associated with adherence to the Green-MED diet, with changes mediating some effects.
Caveats
- The study had a high proportion of male participants (88%), which may limit generalizability to females.
- Dietary intake assessments relied on self-reports, which can introduce bias.
- Only a subset of fecal samples underwent detailed sequencing, potentially limiting the detection of functional differences.
Definitions
- Green-Mediterranean diet: A modified Mediterranean diet emphasizing plant-based foods, reduced red and processed meat, and inclusion of polyphenol-rich foods.
- gut microbiome: The community of microorganisms residing in the gastrointestinal tract, influencing health and disease.
- cardiometabolic health: A term encompassing the interrelated conditions of cardiovascular disease, metabolic syndrome, and diabetes.
Simplified
Background
The Mediterranean (MED) diet, high in nuts, vegetables, and legumes and low in red meat intake, is recommended for the prevention of cardiometabolic diseases [1–3].
It has been reported that adherence to MED dietary patterns is associated with a distinct gut microbiome profile [4]. An 8-week study (n = 82) of an isocaloric MED diet intervention revealed changes in multiple microbial features in the gut, including an increased abundance of major dietary fiber metabolizers such as Bacteroides cellulosilyticus [5], concurrent with beneficial changes in cardiometabolic biomarkers [6]. However, to the best of our knowledge, no study has tested the degree to which the MED diet exerts its beneficial effects on cardiometabolic health through modulation of the gut microbiome.
While randomized controlled trials showed that dietary interventions targeting overall eating patterns lead to a substantial reduction in cardiometabolic disease risk [2, 7], optimizing dietary regimens known to promote health could have vast implications. Recently, we tested the effects of a modified MED diet, the Green-MED diet, which further increases the plant-sourced components while minimizing red and processed meat intake [8, 9]. The two key components of the Green-MED diet are daily intake of green tea adding large quantities of polyphenol compounds [10] and Mankai, a novel Wolffia globosa strain, providing additional polyphenols and plant-based protein.
We recently reported the beneficial effects of this modified MED diet on the reduction of liver fat, a marker of insulin resistance and increased cardiovascular risk [11, 12], beyond the established benefits of the traditional MED diet [8]. Furthermore, we suggested that the Green-MED diet might modify the gut microbiome in a manner that optimizes the efficacy of autologous fecal microbiota transplantation for the prevention of weight regain [13].
Only a minor proportion of ingested polyphenols are absorbed in the upper gastrointestinal tract, with estimates ranging between 5 and 10%. The remaining fraction of polyphenols are metabolized by the gut microbiome, with microbial-derived phenolic metabolites subsequently absorbed by the portal vein, ultimately reaching the systemic circulation as bioactive compounds [14, 15]. It was previously shown that polyphenols are not merely subject to microbial metabolism, but that certain doses of selected polyphenols may modify the gut microbial composition, thus forming a bidirectional relationship [14, 15]. For example, flavonoids, the dominant class of polyphenol compounds found in Mankai [16], are metabolized and transformed by the gut microbiome, funneling this diverse range of compounds to a reduced number of metabolites, while altering the taxonomic structure of the gut microbial community [17].
Another group of metabolites shown to be affected by the gut microbiome is circulating branched-chain amino acids (BCAA). It was previously established that dysregulation of BCAA metabolism is a hallmark of obesity and insulin resistance, with obese/insulin-resistant individuals being characterized by elevated levels of serum BCAA [18]. It was later suggested that variation in circulating BCAA levels, and subsequently the host’s metabolic state, are driven by microbial biosynthesis and uptake in the gut microbiome [19].
The previously reported DIRECT-PLUS dietary trial was a three-armed randomized trial including 294 subjects, with primary outcomes that included abdominal fat deposition and obesity [8, 9]. In the current analysis, predefined as a secondary outcome, we evaluated the effect of MED-based dietary interventions on the gut microbiome composition and function. We sought to identify specific microbial genera and metabolic pathways modified by the plant-enriched Green-MED diet while shedding light on the important role of the rare, “non-core” taxa as drivers of change in composition. Lastly, we aimed to determine the mediating role of the gut microbiome in the association between dietary adherence and the notable cardiometabolic risk reduction introduced by the Green-MED diet.
Methods
Recruitment to the study
The DIRECT-PLUS (clinicaltrials.gov ID: NCT03020186↗) was a prospectively registered clinical trial carried out in the Nuclear Research Center Negev, Dimona, Israel, a facility with an in-house dining room and medical clinic. The first subject was enrolled on January 28, 2017, and the last subject was enrolled on April 30, 2017. The trial was initiated and conducted in a single phase between May 2017 and November 2018.
Inclusion criteria were age > 30 years with abdominal obesity [waist circumference (WC): men > 102 cm, women > 88 cm] or dyslipidemia [triglycerides (TG) > 150mg/dL and high-density lipoprotein cholesterol (HDLc) ≤ 40 mg/dL for men, ≤ 50 mg/dL for women]. Exclusion criteria are detailed in Additional file: Methods. The study was approved by the Soroka University Medical Centre Institutional Review Board. Participants provided written informed consent and received no compensation. 1
Study design, randomization, and intervention
Design and intervention of the DIRECT-PLUS, as well as primary outcomes of the trial (defined as abdominal fat tissues, body weight, and waist circumference), were previously reported in detail [8, 9]. This report is focused on the analysis of the trials’ secondary outcomes, including gut microbiome profile, lipid profile, glycemic control, inflammatory state, and cardiometabolic risk.
All eligible participants were randomized in a 1:1:1 ratio, stratified by gender and work site (to ensure equal workplace-related lifestyle features between groups), into one of the three intervention groups: healthy dietary guidelines (HDG), MED, Green-MED, all combined with physical activity accommodation (Additional file: Table S1). 2
The trial was conducted simultaneously, with all participants randomized and followed in a single phase. Participants were aware of their assigned dietary group (open-label protocol). All participants were provided with a free gym membership and were guided to engage in moderate-intensity physical activity (Additional file: Methods). 1
Dietary guidance was provided as follows, with sessions identical in time and guidance across intervention groups to achieve equivalent intervention intensity: (1) HDG group: In addition to physical activity, participants received standard nutritional counseling based on the Harvard T. H. Chan School of Public Health’s “The Nutrition Source,” a website sharing informed nutrition knowledge with the public, providing science-based guidance for healthy living. (2) MED group: In addition to physical activity, participants were instructed to adopt a calorie-restricted Mediterranean diet [2, 20]. The MED diet assigned was rich in vegetables, with poultry and fish replacing beef and lamb intake. The diet also included 28g/day of walnuts.
(3) Green-MED group: In addition to physical activity and the provision of 28g/day walnuts, the Green-MED diet was further restricted in processed and red meat and was richer in plants and polyphenols. The participants were provided with 3–4 cups/day of green tea and 100g/day of frozen Wolffia globosa (Mankai strain) plant frozen cubes, as a green shake for dinner. Both the MED and green-MED diets were equally calorie-restricted (1500–1800kcal/day for men and 1200–1400kcal/day for women). Two hundred eighty-six subjects, all of which provided a fecal sample at baseline (97%), were included in the current analyses.
Details regarding the lifestyle interventions and motivation techniques are provided in Additional file: Methods. 1
Anthropometric measurements and blood and stool sample collection
Measurements and sample collection were obtained at baseline and 6 months. Participants were weighed without shoes to the nearest 0.1 kg. Waist circumference was measured halfway between the last rib and the iliac crest to the nearest millimeter.
Two blood pressure (BP) measurements and pulse were recorded after resting, using an automatic BP monitor (Accutorr-4; Datascope). Blood pressure was calculated as the mean of the two measurements.
Blood samples were obtained after a 12-h fast, centrifuged, and stored at –80°C pending analyses; further information regarding the assays can be found in Additional file 1: Methods. The 10-year Framingham risk score (FRS) was calculated based on gender, age, total cholesterol and HDLc levels, systolic BP and antihypertensive drug dependence, smoking status, and diabetes [21].
Fecal samples were collected at the study site, immediately frozen to –20°C for 1–3 days, then transferred to –80°C pending DNA extraction. Participants prescribed antibiotic therapy 2 months prior to the delivery of baseline fecal samples were excluded from all microbiome analyses. All blood biochemical assays were performed at the University of Leipzig, Germany. 16S rRNA gene sequencing was performed at Fondazione Edmund Mach, Italy. Fecal metagenomic sequencing was performed at CoreBiome, MN, USA.
16S rRNA sequencing and processing
To characterize the microbiome composition and to assess taxon abundance, fecal DNA from 546 fecal samples was extracted and 16S rRNA gene amplicon sequencing was performed by targeting the V3–V4 variable region of approximately 460-bp length using the bacterial primer pair 341F (5′ CCTACGGGNGGCWGCAG 3′) and 806R (5′ GACTACNVGGGTWTCTAATCC 3′). Paired-end sequencing reads were quality filtered, trimmed, de-noised, and merged using the DADA2 software version 1.14.0 [22]. To reject low-quality bases, forward reads were truncated to 280 bases and 25 bases were cropped from start. Reverse reads were truncated to 220 bases and 14 bases were cropped from start. Reads with a number of expected errors higher than two were discarded. Taxonomy was assigned to all identified amplicon sequence variants (ASVs) based on the SILVA rRNA reference database version 132 [23].
Shotgun metagenomics sequencing and processing
To characterize the microbiome’s functional capacity, a subset of 173 samples (31% of samples) underwent shotgun metagenomics. Samples were collected from 90 participants as part of the pre-transplantation sampling of the DIRECT-PLUS autologous fecal microbiota transplantation study [13], with 31 samples from the HDG group, 69 from the MED group, and 73 from the Green-MED group. DNA was extracted, sequenced, and normalized with an average depth of 15.4±2.6 million reads per sample (mean±standard deviation). DNA sequences were aligned using an accelerated version of the Needleman-Wunsch algorithm to a curated database containing all representative genomes in RefSeq v86 [24]. Each input sequence was assigned the lowest common ancestor that was consistent across at least 80% of all reference sequences tied for best hit. Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology groups were observed directly using alignment against a gene database derived from the strain database used above. Further details regarding DNA extraction and bioinformatics pipelines can be found in Additional file 1: Methods.
Dietary and adherence assessment
Assessment of nutritional intake and lifestyle habits was performed using self-reported food frequency questionnaires administered through a computer at baseline and after 6 months [25, 26]. Frequencies and portions of each individual food item were converted to average daily intake for each participant at each time point. Average daily energy intake was calculated by multiplying the frequency of consumption of each item by its caloric content and summing across all foods.
We then applied a Green-MED adherence score to measure the degree of adherence to the Green-MED pattern. The adherence score was created based on the MED index by Willett et al. [27] and Trichopoulou et al. [28] and modified to suit the Green-MED intervention. The adherence score was based on the intake of 9 items: walnuts, vegetables, processed meat, red meat, legumes, fruits, fish, green tea, and Mankai. Each component’s daily intake was normalized to the average daily intake (component intake in grams/total calories). For beneficial components (vegetables, legumes, fruit, walnuts, fish, green tea, and Mankai), individuals whose consumption was below the median were assigned a value of 0, and those whose consumption was at or above the median were assigned a value of 1. For red/processed meat intake, participants whose consumption was below the median were assigned a value of 1, and those at or above the median were assigned a value of 0. The final Green-MED score ranged from 0 (minimal adherence) to 9 (perfect adherence).
Statistical analysis
All statistical analyses were performed in R v3.5.0. To assess the effect of dietary intervention on the composition of microbial communities, we applied permutational multivariate analysis of variance (PERMANOVA) on dissimilarity matrices of weighted UniFrac, using the adonis function from the R package vegan with 999 permutations [29]. To further evaluate the intervention effect on the microbiome composition over time, a PERMANOVA model was applied with the interaction term of lifestyle intervention and group as independent variables (group X time), both for all groups together and for all three combinations of group pairs in a leave-one-out manner.
To quantify the percentage of variance explained by change in biomarkers and dietary variables on the gut microbiome’s composition, vectors of person-specific ASV changes were calculated for each ASV per individual as follows: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Log}_2\frac{ASV_{t6}}{ASV_{t0}}$$\end{document}Log2ASVt6ASVt0. Euclidian distances between the vectors were calculated, forming a distance matrix between all individuals. Next, a PERMANOVA model was applied with the distance matrix as the dependent variable and percent change of biomarkers or dietary intake score as independent variables. All the P values from the PERMANOVA were corrected for multiple comparisons using the Benjamini-Hochberg procedure. For all analyses stratified by core/non-core taxonomic features, the cutoff was defined as >50% prevalence [30]. To assess to what extent the intervention effect on the gut microbiome was convergent by nature, we generated a distance matrix between all subjects at 6 months, measured in weighted UniFrac and based on the 16S rRNA sequences. We next compared the observed differences between the mean within-group and between-group dissimilarities, to a null model generated by 106 permutations of random sample labeling. For per-feature analysis, we used the R package MaAsLin2 [31], with microbial features, taxonomic (16S rRNA) or pathways (shotgun sequencing) filtered, leaving features with a relative abundance of 10−5 in at least 5% of samples. To detect differences in changes of microbial features between groups over time, we built linear mixed models that include group, time, and their product term, with participant identifier as random. A significant P value for the product term indicates that over time changes in microbial features differ between groups. For changes in biomarkers and specific microbial features, we determined the gradual dietary effect using Kendall’s tau-b correlation on the percent change in each biomarker across dietary groups and used generalized linear models to assess between-group differences. Dietary groups were considered a ranked variable, taking into account the increase in polyphenols and decrease in red and processed meat across the groups. To assess the correlation between specific taxa and biomarkers, the Spearman correlation test was used.
Assessing the mediatory effect of the microbiome, we first employed the meditation analysis suggested by Imai et al. [32] by the “mediation” package in R [33], with the Green-MED score being the independent variable, PCo1 loading of the log2 change matrix accounting for mediation by microbiome composition change, and biomarkers being the dependent variables.
We next used the high-dimensional mediation analysis suggested by Zhang et al. [34], using the HIMA package in R [34], using the same dependent/independent variables, with the change of addressing the high-dimensional nature of the microbiome data, and identifying changes in specific taxa as mediators. In the per-feature tests, unless otherwise noted, all the P values were corrected for multiple comparisons using the Benjamini-Hochberg procedure. All multiple hypothesis testing was controlled for false discovery rate with a target of q < 0.25.
Results
Baseline characteristics of participants and adherence rates in the DIRECT-PLUS trial
Two hundred ninety-four subjects were recruited in the DIRECT-PLUS weight-loss trial [8, 9] and randomized to one of the three lifestyle intervention groups.
Changes in gut microbiome composition following 6-month dietary interventions.Study design, adherence rate, and total sample count by category.PCoA of all samples by group and time. The Green-MED group had undergone the most prominent compositional shift, with a significant interaction between group and time.values denote the significance level of time-group interaction terms with all three groups (bold) and group-pair combinations by leave-one-out analysis (PERMANOVA).,Heatmap of changes in genus-level bacteria () and KEGG metabolizing modules () by lifestyle intervention group. For each cell, colors indicate the within-group change coefficient by MaAsLin2 between baseline and 6 months. The three dietary arms lead to distinct changes in patterns across genera and microbial metabolic pathways A B C D C D P
The Green-MED diet modifies the composition and function of the gut microbiome
Following 6 months of dieting, all three groups had undergone a significant change in microbiome composition, as assessed by 16S rRNA sequencing, though it appeared that the Green-MED dieters exhibited the most prominent change, with a significant interaction between time (baseline to six months) and group (P of interaction = 0.01) as revealed by principal coordinates analysis (PCoA) based on weighted UniFrac distances. To compare the compositional change between groups in a subject-specific manner, we created a log2 change ASV matrix between the two timepoints of each subject. The Green-MED group differed significantly in its change compared to the other two groups (PGreen-MED vs HDG = 0.001, PGreen-MED vs MED = 0.04), whereas no significant difference was observed between the MED and HDG groups (PMED vs HDG = 0.19; Fig. 1B and Additional file 4: Fig. S3). We next evaluated the effect of diet on change at the genus level of bacteria and KEGG pathway modules. At month 6, different dietary interventions led to distinct changes across different phyla and functional groups (Fig. 1C, D).
The non-core microbiome as a driver of general composition change
We next asked the question whether the dietary regimens modified the composition of the microbiome in a converging manner, i.e., whether changes in the microbiome compositions of participants resulted in the subject’s gut microbiomes being more similar to each other during the intervention, and to what extent the converging effect was affected by the dietary regimen. Comparing within-group and between-group dissimilarities to a null model, we observed that the Green-MED dieters’ microbiome composition became more similar following the diet, as compared to the HDG and MED groups. Furthermore, performing the same comparison only on the non-core taxonomic features, it could be observed that convergence was more dominant in the non-core microbiome (Fig. 2C).
The non-core microbiome as a driver of general composition change.All ASVs with prevalence above 3%, ordered by prevalence across all samples. “core” taxa were defined as > 50% prevalence, with the remaining taxa defined as “non-core.”PCoA stratified by the non-core (<50% prevalence, top) and core microbiome (>50% prevalence, bottom). The compositional shift of the Green-MED group originated from the non-core (rare) taxa with a significant interaction between time and group in PERMANOVA.Differences in within-group dissimilarity and between-group dissimilarity, by non-core (green points) and all taxa (gray points) microbiome composition. The Green-MED dieters’ microbiome composition assimilated each other at the 6-month point. This effect was even more dominant in the non-core fraction of the microbiome. Violin plots describe the distribution of 10random permutations of the same measure, shuffling sample labels at each iteration A B C 6
The Green-MED diet increases abundance ofand BCAA degrading pathways and decreasesand BCAA biosynthesis pathways in the gut microbiome Prevotella Bifidobacterium
Comparing changes in microbial metabolic pathways (KEGG modules) that utilize dietary components (i.e., carbohydrate, lipid, amino acid metabolism), we found a stepwise increase in BCAA degrading pathways (isoleucine degradation, Kendall tau; P = 0.02) and a decrease in pathways of BCAA biosynthesis (valine and isoleucine biosynthesis, Kendall tau; P < 0.05) comparing MED diet groups to the HDG group. MED diet interventions also led to an enrichment in lipid metabolizing pathways, possibly due to the increased lipid intake and decreased carbohydrate intake commensurate with MED diets (Fig. 3C, D).
MED diets’ effects on microbial genera and metabolic pathways.Heatmap of changes in genus-level bacteria by lifestyle intervention group, for taxa with significant time*group interaction. For each cell, colors indicate the within-group change coefficient over time and asterisks denote significance. Black-White annotations on the left denote the significance of between-group change difference (by MaAsLin2 time*group interaction).Green-MED was distinguished by several changes, including an increase inabundance and a decrease inHeatmap of changes in metabolism-related KEGG modules by lifestyle intervention group, for pathways with significant time*group interaction. For each cell, colors indicate the within-group change coefficient over time and asterisks denote significance. Black-white annotations on the left denote the significance of between-group change difference (by MaAsLin2 time*group interaction). Colors on the left denote the metabolism category.The Green-MED diet was associated with increased BCAA degradation pathways and decreased BCAA biosynthesis****< 0.001, ***< 0.01, **< 0.05, *< 0.25 A B C D Prevotella Bifidobacterium. . FDR FDR FDR FDR
The beneficial effects of MED diets on cardiometabolic risk factors were attributable to plant-sourced components
Evaluating the association between changes in biomarkers and shifts in microbiome composition (assessed by 16S rRNA sequencing), body-weight change explained the highest proportion of variation (0.9%). Significant associations with compositional changes were observed in the change in several biomarkers, including plasma leptin (0.7%), waist circumference (0.66%), systolic blood pressure (0.5%), HOMA-IR (0.5%), and Framingham risk score (0.65%, Fig. 4).
Six-month changes in cardiometabolic markers and their associations with changes in global microbiome composition. Left: Forest plot of percent change in cardiometabolic health markers following 6 months of dietary intervention across groups. Significant changes in biomarkers are shown (FDR-correctedKendell tau < 0.25). Right: Associations between changes in the corresponding biomarker on the left and composition change (log2 fold change of all ASVs) by performing PERMANOVA test on a Euclidean distance matrix. The Green-MED diet improved cardiometabolic health, with changes in a subset of biomarkers, including body weight, waist circumference, blood pressure, and glycemic profile, being associated with microbiome compositional shift. ****< 0.001, ***< 0.01, **< 0.05, *< 0.25. FRS, Framingham Risk Score; GGT, gamma glutamyl transferase; SBP/DBP, systolic/diastolic blood pressure; MAP, mean arterial pressure; HDG, healthy dietary guidance; HDLc, high-density lipoprotein cholesterol; HOMA-IR, homoeostatic model assessment of insulin resistance; LDLc, low-density lipoprotein cholesterol; MED, Mediterranean; TC, total cholesterol; TG, triglycerides P FDR FDR FDR FDR
Changes in gut microbiome mediate the protective effect of adherence to the Green-MED diet on cardiometabolic risk factors
We next identified 21 specific genus-level features (by 16S rRNA sequencing) associated with the Green-MED diet, i.e., features with 6-month changes in abundance significantly associated with the Green-MED adherence score or significantly different across the groups (FDR-corrected q < 0.25). Among Green-MED-associated taxonomic features, we found that change in Prevotella abundance was positively associated with the Green-MED adherence score, while Bifidobacterium was inversely associated with the score (Fig. 5B).
Following the identification of Green-MED-associated bacteria, we examined the association of the specific taxonomic features with cardiometabolic risk markers and found that a reduction in Bifidobacterium abundance was significantly associated with weight loss (Spearman’s rho = 0.31, FDR-corrected q < 0.001). We also found that increased abundance of several genera from the Ruminococcaceae family was associated with weight loss, namely Ruminococcaceae UCG10 (Spearman’s rho = −0.28, FDR-corrected q < 0.001), Ruminococcaceae NK4A214 group (rho = −0.21, FDR-corrected q = 0.005), and Ruminococcaceae UCG3 (rho = −0.14, FDR-corrected q = 0.08). In addition, a reduction in Bifidobacterium and Dorea and an increase in genera from the Ruminococcaceae and Lachnospiraceae families were associated with a reduction in Framingham risk score (Fig. 5B).
In a similar analysis performed on KEGG modules (assessed by shotgun sequencing), we identified significant associations of the Green-MED adherence score with 10 microbial metabolic pathways (FDR-corrected q < 0.25). The leucine degradation pathway was positively correlated with Green-MED adherence (Spearman’s rho = −0.21, FDR-corrected q = 0.165). Reduction in the levels of the cysteine biosynthesis module was associated with reductions of body weight (Spearman’s rho = 0.33, FDR-corrected q = 0.02) and Framingham risk score (Spearman’s rho = 0.25, P = 0.03), an observation that is consistent with a previous report demonstrating that high abundance of cysteine biosynthesis pathways is associated with unsuccessful weight loss [41] (Fig. 5C).
As we hypothesized that the dietary interventions, through modulation of the gut microbiome, exert their effects on cardiometabolic risk, we next sought to quantify the mediation effect. To this aim, we estimated the causal mediation effect as proposed by Imai et al. [32], using the first principal component of variation of the ASV change matrix as a summary measure of the gut microbiome composition. The association between the Green-MED diet adherence and weight loss (12% mediation by taxonomic change), WC loss (10%), and FRS reduction (18%) was found to be mediated by the change in the gut microbial composition. When assessing the mediatory effect specific to the change of Green-MED-associated genera, Chemerin change was also found to be significantly mediated (Fig. 5B, D). In order to assess the mediatory effect of specific bacteria, we applied a high-dimensional mediation analysis [34]. Among the significant mediatory taxa, Bifidobacterium change was found to mediate the Green-MED effect on weight loss (4%) and Ruminococcaceae UCG10 mediated the effect on FRS reduction (1.3%). Changes in the cysteine biosynthesis pathway were a prominent mediator of weight loss (14.2%) and FRS reduction (29.3%).
Green-MED adherence, cardiometabolic markers, and the association with specific bacteria and microbial metabolic pathways.Left: Forest plot of Green-MED adherence score and its components following 6 months of dietary intervention across groups. Adherence was assessed using a 9-dimensional index, ranging from 0 (non-adherence) to 9 (full adherence). Right: Associations between changes in the corresponding dietary component on the left and compositional change (log2 fold change of all ASVs) by performing PERMANOVA test on the Euclidean distance matrix. Green-MED dieters adhered to their assigned intervention, with intake of specific components (tea and walnuts) being associated with microbiome compositional shift.,Heatmap showing the association patterns of changes in the Green-MED adherence score and change in genus taxa () and KEGG modules (). Features are arranged from top to bottom by their correlations with the adherence score. For each cell, colors indicate the Spearman rho values. Several links between change in microbial features and change in biomarkers were observed. Specifically, changes in specific Green-MED-associated bacteria and microbial metabolic, e.g.,and Leucin (BCAA) degradation, were associated with weight change following the intervention.Mediation analysis: assessing the proportional mediatory effect of change in all genus-level taxa (calculated as #1 principal component of their change matrix) in the association between lifestyle intervention and weight change. Weight loss and Framingham risk score were found to be partially mediated by changes in the gut microbiome. ****< 0.001, ***< 0.01, **< 0.05, *< 0.25,value < 0.05. FRS, Framingham risk score; MAP, mean arterial pressure; HDG, healthy dietary guidance; HOMA-IR, homoeostatic model assessment of insulin resistance; LDLc, low-density lipoprotein cholesterol; MED, Mediterranean A B C B C D Bifidobacterium FDR FDR FDR FDR P #
Discussion
In the current study, we demonstrate the effect of MED and Green-MED diets on the gut microbiome and determine the mediatory role of the microbiome in the associations between adherence to the Green-MED diet and subsequent improvement in markers of cardiometabolic risk. We showed that the Green-MED diet induced a prominent change in the gut microbiome composition, driven by the low-prevalent “non-core” fraction of the gut microbiome. This change was convergent by nature, with the taxonomic composition of the microbiome of individuals in the Green-MED becoming more similar to each other. The different diets were also characterized by distinct taxonomic and functional changes, inducing a gradual increase in the abundance of Prevotella and BCAA degradation pathway concomitant with the increment of plant components in MED diets, and reductions in the abundance of Bifidobacterium and BCAA biosynthesis pathways when comparing MED diets vs. the control diet. The MED and Green-MED diets improved cardiometabolic markers, including the Framingham risk score, body weight, blood pressure, HOMA-IR, and fasting plasma leptin. These beneficial changes in levels of cardiometabolic biomarkers were associated with a concurrent shift in the gut microbiome composition, estimated to mediate 12% of body-weight reduction and 18% of Framingham risk score reduction.
Our study has several limitations. The high proportion of male participants in the trial (88%) may limit the generalizability of our results to females. However, a previous report showed that, in western populations, gut microbiome composition is not associated with gender [42]. We addressed this limitation by employing gender stratification at the trial’s randomization. Furthermore, assessments of dietary intake and adherence to intervention regimes were based on self-reports. This limitation was addressed by the measurement of serum folic acid, a biomarker reflecting green leaf consumption [38], that increased across groups, by the use of previously validated dietary questionnaires [26, 43] and by focusing the analysis on individual-level dietary changes. As only a subset of approximately one-third of the fecal samples were sequenced by shotgun metagenomics, with a relatively shallow coverage of 15 million reads per sample on average, this study was possibly underpowered to detect subtle functional differences between the gut microbiome of the three intervention groups.
The strengths of the study include the intense dietary instruction and multidisciplinary supervision, the simultaneous design, and the high retention rate.
Although the function and composition of the gut microbiome have dynamic features, with some changes observed on a daily scale, most of this variation is person-specific [44]. Moreover, while short-term dietary changes induce rapid and brief changes in microbiome composition, they tend to regress to the baseline personal microbiome structure [45]. Long-term dietary habits however are considered a major determinant of gut microbiome composition, possibly due to consistent selective pressure induced by the diet and specific nutrients [46, 47]. A recent 1-year dietary trial reported that the MED diet showed a notable impact on the composition of the gut microbiome [48]. In the current analysis, we showed that long-term adherence to the Green-MED diet prompted a greater composition change than the shift induced by the traditional MED diet. Furthermore, we demonstrated that the Green-MED composition change was convergent by nature (i.e., the microbiome of the Green-MED dieters became more similar to each other), corroborating the hypothesis that change was driven by a selective pressure [49]. Interestingly enough, convergence was likely driven by the rare, non-core fraction of the microbiome, presumably utilizing its pool of genetic resources given a change in environmental conditions [50].
The major determinants differentiating the Green-MED diet from the MED diet are the supplementation of Mankai and green tea intake while reducing the consumption of red and processed meat.
Mankai was previously established as a high-quality substitute source for animal protein in a controlled test meal, with blood concentrations of essential amino acids being similar after Mankai intake to that of well-established animal and plant protein sources. As an exception, Mankai intake induced only a minor increase in circulating BCAA levels [51]. This observation is in accord with our current analysis, showing that the Green-MED diet was associated with a reduction in the gut microbiomes’ capacity for BCAA biosynthesis and an increase in BCAA degradation pathways.
Furthermore, both Mankai and green tea were shown to be rich sources of polyphenols [16, 52, 53], with Mankai containing ~200 different phenolic metabolites. We previously reported that in a Mankai-supplemented artificial gut bioreactor, the functional profile of the gut microbiome shifted towards polyphenol degradation, e.g., degradation of the phenolic compounds protocatechuate and catechol [52]. Of note, in this study, we were not able to identify a differential effect of the Green-MED diet on polyphenol-related pathways, possibly due to our power limitation mentioned above.
While the Mankai plant is considered fiber-rich, with a fiber content of ~25% [54], the Green-MED diet lead to an increase in Prevotella abundance, a genus previously associated with high fiber intake [55] and vegetarian diets [47, 56]. Moreover, low abundance of Prevotella was shown to be a hallmark for the transition to a western lifestyle [56]. Previous studies have reached contradictory conclusions regarding the role of the dominant Prevotella species, Prevotella copri, in host metabolism. It was recently suggested that the abundance of Prevotella copri, as an indication of the preexisting, diet-related biosynthetic capacity of the gut microbiome, potentially attenuates the protective effect of the Mediterranean diet on cardiometabolic disease risk [5]. Furthermore, De Vadder et al. demonstrated that elevated Prevotella copri abundance is associated with improvements in glucose metabolism and insulin sensitivity by the production of succinate [57], while a study by Pedersen et al. suggested that Prevotella copri abundance is associated with insulin resistance, presumably due to increased BCAA biosynthesis [19]. Here, we report that across MED diets with a gradual increment in plant components, intervention intensity was associated with a reduction in BCAA biosynthesis, an increase in BCAA degradation pathways, and increased insulin sensitivity, measured as HOMA-IR. These findings support the hypothesis that Prevotella may not be beneficial or deleterious per se, but should be examined in light of other environmental and host-specific factors, requiring further investigations [5, 58]. The association of Bifidobacterium with cardiometabolic disease risk has also been inconsistent in previous studies, with an increase in Bifidobacterium abundance linked to both successful [59] and unsuccessful [36] weight loss. Our study demonstrated that in the setting of a plant-enriched (green) MED diet, weight loss is associated with a reduction in Bifidobacterium.
The associations between dietary patterns, the gut microbiome, and cardiometabolic health have been extensively evaluated in recent years [37, 60, 61], with several studies exploring microbial modifications related to beneficial diets, including the MED diet [4, 5, 48]. Our study links MED diets with several taxa, including the depletion of Blautia [48] and enrichment of Bacteroides [48], adding to the body of evidence by examining the effect of expanding the plant fraction of the diet while minimizing meat intake. We showed that overall, the gut microbiome contributes a small but significant mediatory effect in the association between MED diets and weight and cardiometabolic risk reduction.
Conclusions
MED diet enriched in plant components and with decreased red and processed meat consumption (Green-MED diet) has extensive effects on the composition and function of the host gut microbiome, with the latter partially mediating the beneficial effects of the diet on cardiometabolic health. In future studies, the diet-microbiome-host interplay should be further explored and could guide novel beneficial modifications in known dietary patterns.
Supplementary Information
Additional file 1. Methods – Exclusion criteria, Physical activity protocol, Lifestyle sessions and motivation techniques, Blood sample analysis, 16s rRNA sequencing pipeline, Metagenomics analysis BoosterShot pipeline.Additional file 2: Table S1. Outline of the lifestyle interventions.Additional file 3: Table S2. Characteristics of the study population.Additional file 4: Figure S1. Flow Diagram of the DIRECT-PLUS trial.Additional file 5. Contains Figures S2 and S3.