Gut microbiota link dietary fiber intake and short-chain fatty acid metabolism with eating behavior

The gut microbiome has been speculated to modulate feeding behavior through multiple factors, including short-chain fatty acids (SCFA). Evidence on this relationship in humans is however lacking. We aimed to explore if specific bacterial genera relate to eating behavior, diet, and SCFA in adults. Moreover, we tested whether eating-related microbiota relate to treatment success in patients after Roux-en-Y gastric bypass (RYGB). Anthropometrics, dietary fiber intake, eating behavior, 16S-rRNA-derived microbiota, and fecal and serum SCFA were correlated in young overweight adults (n = 27 (9 F), 21–36 years, BMI 25–31 kg/m2). Correlated genera were compared in RYGB (n = 23 (16 F), 41–70 years, BMI 25–62 kg/m2) and control patients (n = 17 (11 F), 26–69 years, BMI 25–48 kg/m2). In young adults, 7 bacteria genera, i.e., Alistipes, Blautia, Clostridiales cluster XVIII, Gemmiger, Roseburia, Ruminococcus, and Streptococcus, correlated with healthier eating behavior, while 5 genera, i.e., Clostridiales cluster IV and XIVb, Collinsella, Fusicatenibacter, and Parabacteroides, correlated with unhealthier eating (all | r | > 0.4, FDR-corrected p < 0.05). Some of these genera including Parabacteroides related to fiber intake and SCFA, and to weight status and treatment response in overweight/obese patients. In this exploratory analysis, specific bacterial genera, particularly Parabacteroides, were associated with weight status and eating behavior in two small, independent and well-characterized cross-sectional samples. These preliminary findings suggest two groups of presumably beneficial and unfavorable genera that relate to eating behavior and weight status, and indicate that dietary fiber and SCFA metabolism may modify these relationships. Larger interventional studies are needed to distinguish correlation from causation.

However, nutrition-or body weight-related microbial changes and their functional relevance are still relatively unclear. In mice, gastric bypass-related differences in the microbiota profile, such as a higher abundance of the genera Escherichia (phylum Proteobacteria) and Akkermansia (phylum Verrucomicrobia), induced weight loss when transferred to germ-free animals [4]. In humans, bariatric surgery similarly led to higher overall microbiota diversity and to higher abundance of the species Escherichia coli and in some studies to further abundance changes within the phylum Bacteroidetes, such as a higher post-surgery ratio of the genera Bacteroides to Prevotella [5] and less Firmicutes (phylum level) or to more Gammaproteobacteria (class level) [6]. The ratio of Bacteroides to Prevotella at baseline predicted dietary weight loss success after 24 weeks in an intervention study in 80 overweight individuals [7]. Further, a oneweek dietary intervention trial in 20 individuals found that microbial composition predicted glycemic response [8].
Human-to-mouse fecal transplant experiments further underline the causal role of specific microbiota to facilitate weight loss [9], and human-to-human fecal microbiota transplantation (FMT) experiments increased insulin sensitivity according to [10]. In a recent human study, accompanied by mouse model data, an individual's microbiota profile, extracted from fecal samples during periods of dietary weight loss, prevented weight regain when transferred back to the same individuum orally, known as autologuos FMT [11].
Mechanistic insights into how specific gut bacteria modulate human eating behavior and weight status are still limited. The gut microbiota is supposed to affect the host's metabolism by altering energy extraction from food, and by modulating dietary or hostderived compounds that modify the metabolic pathways of the host [12]. For example, short-chain fatty acids (SCFA) are excreted by certain gut bacteria as a result of carbohydrate fermentation, and SCFA stimulate the secretion of anorexigenic hormones, such as peptide YY (peptide tyrosine tyrosine or PYY) and glucagonlike-peptide-1 (GLP-1) in the colon, which further signal to hypothalamic nuclei as one mechanism of homeostatic regulation [13]. SCFA can also cross the blood-brain barrier and act as signaling molecules in the brain to directly modulate appetite and food-decision making [1]. First interventional studies showed that intake of butyrate (one type of SCFA) or the butyrate-producing bacteria Akkermansia spp. exert beneficial effects on body weight depending on treatment intention in humans [14] and on brain functions in mice [15], including reduced food intake [16]. Notably, specific pre-biotic nutrients, such as dietary fibers, are known to nourish SCFA-producing bacteria in the gut, rendering diet a potent modifier of gut-brain signaling [17].
In sum, the gut microbiome may influence feeding behavior, e.g., by modulating reward and homeostatic signaling [18,19] and by stimulating the vagal nerve [20], in particular in dysregulated biological systems, such as in food addiction [21] or eating disorders [2]. Yet, direct knowledge if specific genera are linked to eating behavior via dietary intake and SCFA in humans is lacking. Here, we asked whether gut microbial diversity and genera abundance relate to eating behavior, and to SCFA metabolites in the colon (feces) and in the periphery (blood) in a homogenous sample of young overweight adults. In addition, we tested whether the abundance of microbiota that related to eating behavior in that overweight sample correlate with weight status, eating behavior, and treatment success (i.e., achieved weight loss) in another sample, i.e., patients at two years after bariatric surgery and control overweight/obese patients.

Samples characteristics and data collection
We included all participants with available microbiota datasets measured at a cross-sectional timepoint from two studies. Sample 1 comprised 27 healthy young overweight adults (9 F, 21-36 years, BMI 25-31 kg/m 2 ) drawn from a randomized clinical trial (Clinical Trials registration NCT03829189), where baseline data was available from ongoing data collection until January 2021. All participants were included if following a typical Western omnivorous diet and thoroughly screened for habitual dietary patterns (exclusion criteria were assessed via an interview at pre-screening included any sort of restrictive diet (incl. vegan, vegetarian, gluten-free, lactose-free, food allergies), regular excessive caffeine intake (more than 6 cups a day), regular alcohol intake (>1.25 L beer/day or equivalent) or smoking >10 cigarettes/day). Estimated nutrient intake represents a Western style omnivorous diet (10.4 ± 3.6 g/day/ 1000 kcal), with lower than recommended fiber intake (intake recommendations by WHO and EU nutritional agencies state >25 g or 25-35 g of dietary fiber per day are required to meet healthy intake levels). For further information on dietary and coffee intake see SI (see Dietary Intake).
All participants donated feces (see SI) for microbiota analysis (Shannon effective [22] and relative abundance of microbiota genera), underwent anthropometric measurements, and filled in questionnaires to quantify eating behavior traits. Also, data on dietary fiber intake, hunger ratings after a standardized meal, and SCFA in blood and feces were available in sample 1 (see below).

Microbiota assessment
To assess microbiota community structure we used 16 S rRNA gene profiling of the fecal samples. Therefore, DNA was extracted and V3-V4 variable regions of the 16 S rRNA genes were amplified by PCR and a library was constructed, followed by paired-end 2x250bp Illumina sequencing. These analyses were done by GENEWIZ Germany GmbH, Leipzig. Next, the inhouse Galaxy server using a pipeline implemented with the DADA2 R package processed raw data in fastq format. For each sample, paired-end reads were joined, low-quality reads were removed, reads were corrected, chimeras removed, and Amplicon Sequence Variants (ASVs) were obtained. Taxonomy was annotated to the ASVs using the RDP database [23]. The read counts per ASV with taxonomic annotation were normalized and relative abundances of each ASV and taxa were calculated using the R scripts Rhea. Visualization of all library-indexed genera was done as in [24] by inhouse written R-tools using ggplot2.

Eating behavior
To characterize eating behavior traits, questionnaires based on self-report were used: the Three-Factor Eating Questionnaire (TFE-Q) (German version, [25]) and the Eating Disorder Examination Questionnaire (EDE-Q) (German version, [26]) as available for sample 1 and sample 2, respectively. The TFE-Q assesses three domains of eating behavior (cognitive restraint, disinhibition, hunger), and the EDE-Q covers the subscales dietary restraint, eating concern, weight concern, and shape concern. Scoring was performed according to the respective manuals.

Additional analyses in sample 1
From all measurements available in sample 1 in the context of the RCT (see above), we additionally considered all available hunger ratings after a standardized meal (3 out of 4 measures, 1 with missing data) and all available dietary fiber intake data (from a quantitative food frequency questionnaire, fiber in g/day and fiber per 1000 kcal). We further considered anthropometric assessments to be of interest in this study and limited those to two major health indicators, i.e., systolic blood pressure (mean of three consecutive measurements) and relative body fat (%) obtained from bioelectrical impedance analysis (see SI). Blood was obtained in fasting state (12 ± 3 h fasted) and samples were centrifuged at 3500 revolutions per minute at 7°C for 6 min. Serum was aliquoted within 1 h of obtainment. Processed aliquots were stored at −80°C until data analysis. For SCFA in blood and stool, analyzed according to [27], we focused on three major and most abundant SCFAs out of eight measured, i.e., acetate, butyrate, and propionate (see SI). All other measures were not considered of interest to the current analyses.

Obesity surgery in sample 2
For sample 2, RYGB (see SI) patients were selected for microbiota analysis based on their response to the surgical treatment. Specifically, RYGB patients were identified from the database of the University of Leipzig if their surgery dated back at least 2 years and all those were further divided in percentiles according to pre-defined relative excessive weight loss (EWL) thresholds defined more conservatively than previous literature (most common <50% EWL at 18 months, according to [28]. This resulted in 23 RYGB patients good responders: sustained EWL > 70%, mean 93% ± 4 SD, range 86-98%, n = 12; bad responders: sustained EWL < 40%, mean 20% ± 13 SD, range 3-35%, n = 11). Next, obese and overweight control patients were selected from the database based on age, sex, and BMI to match those two groups of RYGB patients. Afterwards, RYGB patients only filled in a series of questionnaires, performed cognitive tests, and donated blood for another study purpose; and fecal samples of all patients were analyzed. From this dataset, we considered of interest to the current analysis the following variables: weight loss after surgery (in kg and in BMI), all available eating questionnaire data (four EDEQ scales, see above), and microbiota genera abundances based on 16 S rRNA sequencing.

Statistical analysis
Correlational analysis. Relative taxa abundance (%) on the genera level was used as primary variables of interest. Non-normally distributed variables were log-or Tukey-transformed, so that skewness of <|1| was reached (for details Supplementary Fig. 1). No observations were eliminated, instead all cases with microbiota data were complete and included. For the main analysis, 20 out of 121 genera were included as they appeared in at least 80% of individuals [29] and fed into a correlation matrix with all variables of interest in sample 1 (37 variables in total, see above), i.e., Shannon index, 3 TFEQ traits, 3 hunger ratings, body fat, systolic blood pressure, dietary fiber intake (g/day and g/1000 kcal), and 3 SCFA each in feces and blood, respectively. All values were FDR-corrected E. Medawar et al. and statistical significance was set to p < 0.05. Those genera that were significantly associated with eating behavior (TFEQ traits and/or hunger ratings, p-FDR < 0.05) were then correlated with weight status and RYGB treatment success in sample 2. Group differences across overweight, obese, good and bad RYGB responders were tested with non-parametric Kruskal-Wallis tests. Further correlations were tested with Pearson's correlation coefficient r for normally distributed variables or with Spearman's rho for non-normally distributed variables. Explorative analysis considerations were addressed according to [30] (see Additional SI).
To further investigate, if the interplay of correlated genera-rather as a holobiont than individually-is determinative of the observed relations, the relation between correlated to non-correlated genera was computed by three composite scores (1)- (3).
Mediation analysis. Using simple mediation analysis using medmod (https://cran.r-project.org/web/packages/medmod/index.html) in RStudio version 3.6.1, we checked for statistical mediation in sample 1 for variables showing bivariate correlations in the following paths: i. fiber -> correlated genera or sumscores -> eating behavior (TFEQ, hunger ratings) ii. eating behavior (TFEQ, hunger ratings) -> fiber -> correlated genera or sumscores iii. correlated genera or sumscores -> SCFA -> eating behavior (TFEQ, hunger ratings) Significance was set to p < 0.05, and the main analysis for sample 1 was corrected for multiple testing using the false-detection rate (FDR)correction. All analyses were performed in RStudio version 3.6.1.

RESULTS
Characteristics of sample 1 and 2 are listed below (see Tables 1-2). Data from post-RYGB patients was on average collected 4.7 ± 1.4 years after surgery. Eating behavior traits varied across both samples, and in sample 2, restrained eating and shape/weight concerns differed between those that achieved long-term excessive weight loss after bariatric surgery compared to those that did not (good vs. bad responders, all W > 58.5, p < 0.001, Table 2, Supplementary Fig. 2).
Overall microbiota diversity at the phylum level was relatively comparable across participants of samples 1 and 2 except higher ratio of Firmicutes to Bacteroidetes in sample 1, and Prevotellaceae and Fusobacteriaceae families were more abundant in patients after RYGB surgery (Fig. 1, Supplementary Figs. 3-4, see SI for details). Additionally, we tested for sex/gender-specific differences in alpha diversity in sample 1 and found none (sample 1: male (n = 18) 111 ± 15, female (n = 7) 110 ± 13, t(13) = −0.08, p < 0.94). Due to limited sample size we refrained from further sexsegregated analyses, yet we encourage future meta-analyses to include our datasets (see open data).
Relation to dietary fiber intake and SCFA Out of the 12 genera that were significantly associated with eating behavior (from now on called "(inversely) health-related" genera), three were associated with lower (Collinsella and Parabacteroides) or higher (Clostridium XVIII) dietary fiber intake (all 0.73 < |r | > 0.49, p-FDR < 0.05, Fig. 3A-C). Moreover, higher dietary fiber intake per se was significantly associated with lower disinhibited eating (r = − 0.58, p-FDR < 0.01) and lower body fat mass (r = −0.75, p-FDR < 0.0001, Fig. 3D-E). SCFA concentrations in feces were highly variable and up tõ 1000 times higher compared to serum for all three measured SCFA (all t(24) > 11.6, p < 0.001). Serum acetate was 2.5 times higher compared to butyrate and propionate in serum (Supplementary Table 1).
We observed that higher abundance of some of the inversely health-related genera correlated with higher levels of different SCFA in feces and serum (all r > 0.50, p-FDR < 0.01). In addition, most health-related genera correlated with some feces and serum SCFA markers, however revealing both positive and negative associations (only those associated with eating behavior were considered, all 0.65 < |r | > 0.44, p-FDR < 0.05, Supplementary Fig.  6). Note, that some genera showed differential correlations within the different SCFA, e.g., higher Alistipes correlated with higher acetate in both feces and serum and with fecal butyrate, but with lower fecal propionate. Moreover, considering the inversely health-related genera, Fusicatenibacter and Parabacteroides correlated significantly with higher fecal concentrations of propionate and acetate, respectively.

Genera sumscore and mediation analyses
The negative sumscore of the five inversely health-related genera abundances resulted in significant correlations for two of the eating traits (cognitive restraint r = 0.59, p-uncorr = 0.001; disinhibition r = 0.65, p-uncorr <0.001, Fig. 5A-B). The positive sumscore of the seven health-related genera abundances showed no significant associations (all p-uncorr <0.95, Supplementary Fig. 8). Neither sumscore correlated with fecal or serum SCFA levels.
Exploratory mediation path analyses of the proposed models did not show statistically significant mediating paths for differences in diet, eating behavior or hunger ratings through differences in Parabacteroides or positive/negative sumscores (Supplementary Tables 2-3, Supplementary Fig. 9). Considering SCFA, similar results emerged, except for acetate: here, while the direct effect c' did not reach significance (ß = −0.3, p = 0.13), higher Parabacteroides abundance was linked with higher postprandial hunger ratings through higher fecal acetate levels (indirect effect, a*b, ß = 0.36, 95% CI [0.05 0.66], p = 0.02, Supplementary Table 4).
Exploratory analysis based on reviewer suggestions showed that, when adjusting the correlational analysis for body fat mass, associations with inversely health-related bacterial genera remained largely significant (TFEQ and Clostridum XIVb, Collinsella, Fusicatenibacter, Parabacteroides, all p < 0.05), yet positively healthrelated correlations do not (only for hunger ratings with Clostridium XVIII and Roseburia, all p < 0.05) (relating to Fig. 2, for details see SI Table 5). For fiber intake associations, only the negative association with Collinsella abundance and TFEQ disinhibition scores remain significant (relating to Fig. 3, for details see SI Table 5). The association between propionate levels in feces and TFEQ cognitive restraint when adjusted for body fat mass is no longer significant (relating to Fig. 4).
Potential confounders of the gut microbiome Besides body fat mass, several confounding factors have been proposed to influence gut microbiota, such as time of day of stool collection [31], seasonality [32], coffee consumption [33], and others. Note that statistical tests showed no significant associations of the above mentioned confounders on alpha diversity in our analysis ( Table 3). Influences of medication was not tested further, since medical products varied largely in sample 1 and was quite similar across patient groups in sample 2.

DISCUSSION
Combining data from two human cross-sectional datasets, this exploratory analysis finds two groups of microbiota genera that were either positively or inversely associated with both healthier eating behavior and anthropometrics (1) in a deeply phenotyped sample of young overweight adults and (2) when comparing microbiota observed in (1) in patients showing a good or bad response two years after bariatric surgery with matched controls, respectively. More specifically, in young overweight adults, 7 bacterial genera, i.e., Alistipes, Blautia, Clostridium XVIII, Gemmiger, Roseburia, Ruminococcus, and Streptococcus, correlated with healthier eating behavior traits and lower subjective hunger ratings, indicating potential benefits for the host metabolism, while 5 bacterial genera, i.e., Clostridum IV, Clostridium XIVb, Collinsella, Fusicatenibacter, and Parabacteroides, correlated with unhealthier eating traits and higher subjective hunger ratings. Collinsella was further related to higher body fat mass and Streptococcus to lower systolic blood pressure. The health-related bacterial genera were also more abundant in the overweight good responder controls, compared to the obese bad responder controls and RYGB-operated patients, while the inversely healthrelated genera showed a less clear distribution across groups, with Parabacteroides being significantly less abundant in good vs. bad RYGB-operated patients. Moreover, relative abundance of Parabacteroides as well as a composite score of all inversely correlated genera, were associated with higher eating restraint and with lower post-operative weight loss across both RYGB groups. Considering diet and SCFA-related pathways, we observed that higher dietary fiber intake in overweight adults correlated with more abundant Clostridium XVIII, and less abundant Collinsella and Parabacteroides, as well as with healthier eating behavior and anthropometrics. While SCFA showed a rather mixed pattern of correlations with the different markers, Fusicatenibacter and Parabacteroides abundance correlated with higher fecal propionate and acetate, respectively, that again correlated with elevated hunger. Contrastingly, higher acetate and butyrate in serum correlated with lower fat mass, indicating a possible inverse association of acetate in feces and serum with respect to health indicators. Together, these results indicate that presumably beneficial and unfavorable microbiota genera relate to eating behavior and weight status, and that dietary fiber intake and SCFA metabolism may modify these relationships.

Bacterial genera
Due to the lack of associations with alpha diversity, except for subjective hunger, it remains difficult to draw strong conclusions on relations of eating behavior and microbial diversity, measured with Shannon index, based on the present BMI-defined overweight sample. The health-related microbiota group is comprised of bacterial genera that have been described as beneficial for the host in previous literature. For example, Alistipes and Blautia were found to produce SCFA [34,35]. Similarily, Gemmiger, Roseburia, and Ruminococcus belong to the families of Ruminococcaceae or Lachnospiraceae, which share a common role as active plant degraders [36]. These positive metabolic effects on the host could eventually contribute to improved adiposity control, as e.g., higher Blautia was correlated to lower body fat [37], Roseburia was linked to lower blood glucose and Ruminococcus to higher weight loss in mice after vertical sleeve gastrectomy via regulation of nuclear Fig. 2 Pearson's correlations between eating behavior traits (TFEQ and hunger ratings) or health indicators and bacterial genera in overweight adults (all | r | > 0.42, all p-FDR < 0.05; sample 1). A inversely health-related genera (blue), B health-related genera (yellow), C Collinsella and body fat mass (black), and D Streptococcus and mean systolic blood pressure (black).
receptor binding of bile acids [38]. A microbial transfer study from human to mice showed obesity-promoting effects of the species C. ramosum, which is part of Clostridium XVIII [39]. However, studies on the genera Clostridium XVIII and Streptococcus in relation to health are scarce. Clostridia are known to be key commensals for gut homeostasis [40], but classification of the genus Clostridium remains challenging due to the high heterogeneity of the listed species [41]. Also, there are currently 50 species identified in the genus Streptococcus alone, rendering different functionality in these genera likely. Yet, we found that Clostridium XVIII abundance related to higher dietary fiber intake, and Streptococcus abundance to lower blood pressure. Indeed, fiber intake related to healthy eating behavior and lower body fat mass in overweight young adults in the present analyses may point towards rather beneficial fiber-correlating Clostridium XVIII and Streptococcus genera species that underly those associations. Moreover, these results underline the potential impact of a fiberrich diet for health indicators. Due to the exclusive occurrence of fiber in plants, fiber-rich diets are oftentimes attributed to plantbased (vegetarian or vegan) diets, and plant-based diets have been shown extensively to be beneficial for weight status, gut, and overall health [42,43].
Considering the inversely health-related group of microbiota, some genera were described to include pathogens, e.g., in Clostridium XIVb the species C. piliforme, the causative agent of Tyzzer's disease [41] and Parabacteroides as an opportunistic pathogen in infectious diseases [44]. Of note, in the Parabacteroides genera, also beneficial species, e.g., P. distasonis, have been described [45]. The anaerobic Collinsella colonizes mucosal surfaces and has recently been reported to degrade potentially toxic food contaminants found in processed foods [46]. While this could be beneficial for the host, unhealthier eating behaviors (such as intake of processed food) and higher body weight could then likely be related to higher abundances of Collinsella. Likewise, studies showed that Collinsella linked to less dietary fiber intake, which is in line with our results in overweight adults, and higher weight loss in cross-sectional [47] and dietary intervention studies [48]. Fusicatenibacter, including the species F. saccharivorans, are strictly anaerobic sugar fermenters, again linking to unhealthier eating behavior and obesity [49]. The genus Clostridum IV however has rather been reported as beneficial SCFA producers, e.g., the species Faecalibacterium prausnitzii (F. prausnitzii), which play a noticeable role in intestinal homeostasis [50]. Yet again, those genera comprise many different species and it can also be speculated that some bacteria species or genera underlying the observed correlations could have likely been taxonomically misplaced [41]. Taken together, the negatively correlated microbiota genera seem to consist on the one hand of pathogens, indicative of a rather pro-inflammatory milieu in participants with higher weight status, which is well in line with our findings showing that higher Parabacteroides correlated with unhealthier eating traits and poorer weight loss maintenance in RYGB patients. On the other hand, those negatively correlated genera are comprised of those bacteria that metabolize processed food and sugars, again indicative of higher weight and unhealthy eating behavior. Future studies now need to integrate microbiota data at the species level and randomized interventional trials are required to eventually understand cause and effect of these eating behavior-microbiota-diet interrelations.

SCFA metabolism
We could not establish reliable links between serum and fecal concentrations of those metabolites. The overall weak relationship might be explained by rapid metabolization of SCFAs, as for example butyrate is rapidly absorbed by the gut mucosa and reaches blood circulation [51], therefore, fecal levels of butyrate may not directly relate to butyrate-producing bacteria abundance nor to serum levels of butyrate. In addition, biosamples of serum and feces were not collected in a time-locked way, therefore a time difference of hours to days might have blurred potential (inverse) correlations. Indeed, it has been shown, that fecal SCFA levels decrease throughout the day due to metabolization and that overnight-fast duration influenced these results [52].
Still, we found that higher fecal SCFA levels (i.e., acetate, butyrate, and propionate) linked to higher subjective hunger ratings and also to higher cognitive restraint (i.e., propionate), whereas lower acetate and butyrate in serum correlated with higher fat mass. Statistical path analyses proposed that higher Parabacteroides abundance link to higher hunger through higher fecal acetate. Bearing in mind that higher fecal SCFA levels may indicate less efficient absorption in the gut, leading to lower SCFA availability in serum [53], these findings are somewhat in line with studies showing reduced appetite and less weight gain after acetate intake [1,54]. Note however, that we did not adjust mediation statistics for multiple testing, rendering false positives likely. In addition, it has been discussed that only a minimal fraction of the colon-derived SCFA directly reaches the brain. Instead, more downstream targets of SCFA signaling might be more important for gut-brain communication, such as SCFAinduced release of GLP-1 and PYY at the gut epithelium, modulation of liver metabolism or indirect signaling via the vagus nerve [1]. Future studies could help to further disentangle the different mechanisms at play by assessing further blood-, tissueor imaging-based biomarkers of these pathways.
In an exploratory analysis, we found that body fat mass explains some of the variance in the observed relationships, especially in those with health-related commensals, and less with inversely health-related ones. Although BMI spanned within a very homogenous overweight status group (25-30 kg/m 2 ), body fat mass was quite variable (7-40%) in sample 1 and showed a significant influence of microbiome-behavior associations. This may hint to body fat as an important determinant of gut-behavior relations [55] that should receive further attention when designing dietary interventions targeting the gut microbiome.
Besides body fat, exploratory analyses showed no effects of sex/ gender or common lifestyle factors on alpha diversity, yet these findings remain speculative because of small size and the crosssectional nature of our analysis, therefore we cannot rule out that these or other factors such as medication might have confounded our analyses. Indeed, some studies reported on gut-modulating effects of nutrient supplementation such as in vitamin D [56] or vitamin B12 [57]. We recommend to document and report potential confounders in all microbiome analyses and encourage future data pooling and meta-analyses including our datasets (see open data).

Limitations
Firstly, all analyses are based on cross-sectional data, therefore no conclusions about causal relationships can be drawn. We performed exploratory analyses centered around core hypotheses with the aim to gain more specific testable hypotheses for upcoming intervention trials. In addition, both samples are limited by size, especially with regard to the larger number of variables of interest. Due to these constraints, more elaborate statistical analyses (such as structural equation modeling) could not be performed. A major strength of this study is the inclusion of two independent samples integrating next-generation sequencing and SCFA metabolomics with psychological markers in well-  Sumscore of inversely health-related genera is positively associated with eating traits (TFEQ). Pearson's correlations shown for microbial genera abundance sumscore for inversely health-related genera with respect to eating behavior outcomes from TFEQ and hunger ratings shown for A) cognitive restraint (r = 0.59, p-uncorr = 0.001) B) disinhibition (r = 0.65, p-uncorr < 0.001) C) hunger score D) 10 minpostprandial hunger E) 40 min-postprandial hunger and F) 65 min-postprandial hunger. Data from sample 1, all p-uncorrected.
characterized adults at risk for future weight gain that yielded similar associations of eating behavior with gut microbiota at the genera level.

CONCLUSION
The combination of data from cross-sectional samples of overweight, obese, and post-bariatric surgery individuals showed multivariate associations between specific bacterial gut genera, particularly beneficial SCFA-producing genera and presumably unfavorable pathogens or sugar-/processed-food digesting bacteria, with anthropometrics, eating traits, and dietary fiber intake. While speculative concerning causality, our results propose key microbiota candidates for diet-gut-brain-behavior interactions in humans and may help to develop novel hypotheses how to prevent and treat unhealthy food craving through microbiotal modulation of the gut-brain axis. Longitudinal and interventional studies integrating metagenomic approaches and functional pathway analysis are needed to disentangle correlation from causality and to further characterize eating behavior-relevant microbiota genera at the species level.

DATA AVAILABILITY
All data generated during this study are included in this published article and its supplementary information files. Raw data cannot be shared due to data protection regulations. Fig. 6 Differences in relative bacterial genera abundance across groups of overweight, obese and post-RYGB adults in sample 2. For A inversely health-related (blue), B positively related genera (yellow), C negative (blue), and D positive sumscores of all related genera (yellow) detected in sample 1. E Correlation of negative genera sumscore and/or Parabacteroides abundance with eating behavior and/or weight loss after RYGB surgery (green: good responders; red: bad responders). RYGB, Roux-en-Y gastric bypass. Mean + SD.  Table 5 Most of these additional exploratory analyses could be performed in sample 1 only, due to limited information in sample 2.
E. Medawar et al.