Next Article in Journal
Untargeted Metabolomics Reveals Metabolic Stress Alleviation by Prepartum Exercise in Transition Dairy Cows
Next Article in Special Issue
Evaluating Cefoperazone-Induced Gut Metabolic Functional Changes in MR1-Deficient Mice
Previous Article in Journal
Amplified Drought and Seasonal Cycle Modulate Quercus pubescens Leaf Metabolome
Previous Article in Special Issue
High-Throughput UHPLC-MS to Screen Metabolites in Feces for Gut Metabolic Health
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 12
Issue 4
10.3390/metabo12040308
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

NMR Metabolomics Reveal Urine Markers of Microbiome Diversity and Identify Benzoate Metabolism as a Mediator between High Microbial Alpha Diversity and Metabolic Health

by
Johannes Hertel
1,*,
Daniel Fässler
1,
Almut Heinken
2,
Frank U. Weiß
3,
Malte Rühlemann
4,
Corinna Bang
4,
Andre Franke
4,
Kathrin Budde
5,
Ann-Kristin Henning
5,
Astrid Petersmann
5,6,
Uwe Völker
7,
Henry Völzke
8,
Ines Thiele
2,9,10,11,
Hans-Jörgen Grabe
1,12,
Markus M. Lerch
3,13,
Matthias Nauck
5,14,
Nele Friedrich
5,14,† and
Fabian Frost
3,†
1
Department of Psychiatry and Psychotherapy, University Medicine Greifswald, D-17475 Greifswald, Germany
2
School of Medicine, National University of Ireland, H91 CF50 Galway, Ireland
3
Department of Internal Medicine A, University Medicine Greifswald, D-17475 Greifswald, Germany
4
Institute of Clinical Molecular Biology, Kiel University, D-24105 Kiel, Germany
5
Institute of Clinical Chemistry and Laboratory Medicine, University Medicine Greifswald, D-17475 Greifswald, Germany
6
Institute of Clinical Chemistry and Laboratory Medicine, University Medicine Oldenburg, D-26129 Oldenburg, Germany
7
Department of Functional Genomics, Interfaculty Institute for Genetics and Functional Genomics, University Medicine Greifswald, D-17475 Greifswald, Germany
8
Institute for Community Medicine, University of Greifswald, D-17475 Greifswald, Germany
9
Discipline of Microbiology, National University of Galway, H91 CF50 Galway, Ireland
10
APC Microbiome Ireland, University College Cork, T12 CY82 Cork, Ireland
11
Ryan Institute, National University of Galway, H91 CF50 Galway, Ireland
12
German Center for Neurodegenerative Diseases (DZNE), Partner Site Rostock/Greifswald, D-17475 Greifswald, Germany
13
Faculty of Medicine, Ludwig-Maximilian University Munich, D-80539 Munich, Germany
14
German Centre for Cardiovascular Research (DZHK), Partner Site, D-17475 Greifswald, Germany
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Metabolites 2022, 12(4), 308; https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12040308
Submission received: 25 February 2022 / Revised: 22 March 2022 / Accepted: 25 March 2022 / Published: 31 March 2022
(This article belongs to the Special Issue Biomarker Discovery through Microbiome Metabolism Analysis)
Browse Figures
Versions Notes

Abstract

:
Microbial metabolites measured using NMR may serve as markers for physiological or pathological host–microbe interactions and possibly mediate the beneficial effects of microbiome diversity. Yet, comprehensive analyses of gut microbiome data and the urine NMR metabolome from large general population cohorts are missing. Here, we report the associations between gut microbiota abundances or metrics of alpha diversity, quantified from stool samples using 16S rRNA gene sequencing, with targeted urine NMR metabolites measures from 951 participants of the Study of Health in Pomerania (SHIP). We detected significant genus–metabolite associations for hippurate, succinate, indoxyl sulfate, and formate. Moreover, while replicating the previously reported association between hippurate and measures of alpha diversity, we identified formate and 4-hydroxyphenylacetate as novel markers of gut microbiome alpha diversity. Next, we predicted the urinary concentrations of each metabolite using genus abundances via an elastic net regression methodology. We found profound associations of the microbiome-based hippurate prediction score with markers of liver injury, inflammation, and metabolic health. Moreover, the microbiome-based prediction score for hippurate completely mediated the clinical association pattern of microbial diversity, hinting at a role of benzoate metabolism underlying the positive associations between high alpha diversity and healthy states. In conclusion, large-scale NMR urine metabolomics delivered novel insights into metabolic host–microbiome interactions, identifying pathways of benzoate metabolism as relevant candidates mediating the beneficial health effects of high microbial alpha diversity.

Graphical Abstract

1. Introduction

An incredible wealth of research has highlighted the importance of the gut microbiome’s metabolic activity to health and disease [1,2,3,4,5]. Microbial metabolites, such as short-chain fatty acids or vitamin derivatives, are beneficial for health [6,7,8], while others, such as hydrogen sulphide [9,10], secondary bile acids [11,12,13], or trimethylamine-N-oxide [14,15], are involved in pathophysiological processes. Yet, our understanding of the complex interplay between the microbiome and host metabolism is incomplete, in particular because there is no bijection between gut microbiome composition and its metabolic activity. As metabolic functions widely overlap across the main bacterial phyla in the gut, two microbial communities may diverge strongly in their composition but be largely equivalent in their metabolic functions [16]. Moreover, since the gut microbiome acts within the unique ecosystem shaped by host factors, such as behaviour [17], genetics [18,19,20,21], and physiology [22,23], two gut microbiomes that are largely equivalent in metabolic function may lead to different health outcomes. Thus, the question arises regarding how alterations in composition translate into differences in metabolic host functions.
The mechanisms conferring the health-promoting effects of high alpha diversity are of particular interest [24]. Multiple studies have revealed associations of reduced alpha diversity with a wide range of disease phenotypes, including gastrointestinal disorders [25,26], obesity [27], or diabetes [28,29,30,31,32]. While the hypothesis exists that high diversity implicates a higher functional redundancy and thus more stable host–microbiome interactions [33], the concrete metabolic traits underlying the relation between microbial diversity and human health remain elusive.
A frequently utilised paradigm to shine a light on the complex interplay between gut microbiome composition and metabolic function in the host is the ‘ome-wide’ association study. In this paradigm, the statistical association pattern between abundances of taxonomic units and metabolite concentrations in the host is derived via sequential multivariable regressions. ‘Ome-wide’ association studies can serve as powerful tools for discovery and hypotheses generation. For example, we previously showed that reduced microbial diversity is associated with long-term microbiome instability, which is a phenotype characterised by an increase in the microbial biosynthesis capability for proinflammatory lipopolysaccharides over time [33]. Moreover, Wilmanski et al. [34] found that blood metabolites were highly predictive of the gut microbiome alpha diversity with prominent roles for human–microbial co-metabolites, such as hippurate, cinnamoylglycine, and p-cresol sulfate. However, microbiome–metabolome association studies on large cohorts are still scarce such that putative associations are often in need of validation.
Following these lines of research, the overarching goal of this study was to broaden the understanding of metabolite–microbiome association patterns by using data from the TREND cohort of the Study of Health in Pomerania [35] (SHIP), which is a population-based cohort study from North-Eastern Germany, where subsets were characterised by targeted 1H nuclear magnetic resonance (NMR) urine metabolomics and faecal 16S rRNA gene sequencing. Importantly, rich phenotype data for the SHIP-TREND participants were also available, allowing for comprehensive covariate adjustments. NMR urine metabolomics is especially valuable in the context of microbiome research since a wide range of topical microbiome-related metabolites, such as indoxyl-3-sulphate, trimethylamine-N-oxide, formate, or acetate, can be readily quantified with high accuracy in the human urine through NMR technology. To the best of our knowledge, no large cohort study has been published using NMR urine metabolomics and 16S rRNA gene-sequencing data in conjunction to explore association patterns between the microbiome and host metabolome.
We herein report (1) the urine metabolite–genus association pattern, (2) the urine metabolome association pattern with measures of microbiome alpha diversity, (3) the results of machine learning modelling predicting the urine metabolites from microbiome abundances, and (4) the clinical association pattern of the microbiome-based hippurate prediction score in relation to the Shannon diversity. Finally, we explore the relationship between microbial benzoate metabolism and alpha diversity via functional annotations of an independent metagenomics dataset. In conjunction, our results highlight the usefulness of NMR metabolomics in exploring host–microbiome association patterns, while indicating that benzoate metabolism is one of the mediating metabolic functions underlying the association between alpha diversity and human health.

2. Results

In SHIP-TREND, gut microbiome quantifications from stool samples using 16S sequencing were available for 3637 individuals, from which, 950 were additionally characterised using NMR urine metabolomics. The sample characteristics are given in Table 1. As known diabetes cases were excluded from NMR characterisation for reasons unrelated to the presented study, diabetes-related comorbidities and physiological traits were under-represented in the sample with NMR characterisation. Thus, the NMR analyses were performed on a predominantly healthy population. For more information on sample selection, see the Methods section.

2.1. NMR Metabolomics Revealed Genus Metabolite Associations

The utilised panel of metabolite quantifications contained 60 metabolites, from which, 42 were included in analyses after excluding all metabolites with more than 50% missing quantifications. All metabolite concentrations were log-transformed after regression-based normalisation [36], and outliers, as defined using the 4-SD rule, were excluded from analyses. The final panel of metabolites included a wide range of metabolites known to be implicated in host–microbiome interactions. Descriptive statistics for all metabolites measured can be found in Supplementary Table S1.
The metabolite concentrations were regressed in multivariable mixed linear regressions on the genus abundances, including known influence factors of the urine metabolome and the microbiome as covariates to reduce residual confounding and using batch as the random effects variable (see the Methods section for details). We ran two sets of metabolite genus associations: (1) we included the genus presence (binary: present vs. absent) as the predictor of interest, (2) we utilised the genus abundance (dimensional, %) as the predictor of interest. Those two sets of analyses are not equivalent, as it is known that sometimes the presence of a certain species or genus may be more informative than its abundance, especially in the case of keystone species contributing important metabolic capabilities to the community [37,38].
After correction for multiple testing, we found three genus presence metabolite associations (Figure 1A,B), with urinary hippurate being positively associated with the presence of Catabacter and Barnesialla and urinary succinate being positively associated with Eisenbergiella. In the domain of genus abundances, we detected ten associations after correction for multiple testing (Figure 1C,D). Hippurate was associated with nine different genera, while formate was additionally negatively associated with the Clostridium XIVa cluster and indoxyl sulfate was positively associated with Escherichia/Shigella. Both sets of associations remained stable upon full covariate adjustments (Figure 1B,D). In the sensitivity analyses, by utilising non-parametric bootstrapping for p-value calculation, all abundance metabolite associations stayed significant. However, the p-value of the indoxyl sulfate Escherichia/Shigella association dropped to nominal significance, indicating that this specific result may need further validation. The summary statistics of the corresponding association analyses can be found in Supplementary Tables S2 and S3.

2.2. NMR Metabolomics Revealed Markers of Microbial Alpha Diversity

Next, we screened the urine metabolome for markers of alpha diversity to investigate species richness and the Shannon entropy following the same regression methodology as above.
After correction for multiple testing, we retrieved two associations regarding Shannon entropy and three associations regarding species richness. The top hit of the analyses was hippurate (Figure 2A,B). Additionally, we found two additional markers of microbiome diversity: formate was positively associated with both species richness and Shannon entropy, while 4-hydroxyphenylacetate was negatively associated with species richness. Complete results on metabolite diversity associations can be found in Supplementary Tables S4 and S5.

2.3. Microbiome-Based Predictions Scores for Urinary Hippurate Mediated the Associations of SHANNON Diversity to Markers of Metabolic Health

Then, we generated prediction scores for each metabolite using the genus abundance data as input for elastic net regressions (see the Methods section for details), assessing the model validity via internal cross-validations. However, the machine learning approach only led to substantial R-squared values for hippurate (Figure 2C), reaching 12.7% of the explained variance, whereas, for all other metabolites, the explained variance was below 3%, as determined in cross-validations (Supplementary Table S6). Therefore, we focused on the microbiome-based hippurate prediction score in further analyses.
Interestingly, the microbiome-based hippurate prediction score was substantially correlated (r = 0.53; 95% CI: 0.49, 0.58; p = 1.07 × 10−71) with measures of alpha diversity (Figure 2C). Furthermore, in fully adjusted models, the microbiome-based hippurate prediction score showed robust negative associations with markers of human diseases (Table 2), including inflammation markers (high-sensitive C-reactive protein), triglyceride levels, and, most prominently, markers of liver injury (plasma concentrations of gamma-glutamyl transferase (GGT), alanine aminotransferase (ALAT), and aspartateaminotransferase (ASAT)). Notably, the association pattern was more pronounced than for both the Shannon entropy and urinary hippurate levels themselves.
Given the strong association of the hippurate prediction score with the Shannon entropy, we tested whether the hippurate prediction score could mediate the associations of the Shannon entropy with liver markers (GGT and ALAT) and triglycerides. Indeed, in all three cases, the microbiome-based hippurate prediction score mediated the associations nearly completely (GGT: 77.6% of the effect mediated, 95% CI: 56.8%, 119.7%; ALAT: 116.7% of the effect mediated, 95% CI: 64.6%, 432.1%; triglycerides: 77.9% of the effect mediated, 95% CI: 54.3%, 132.5%). Note that mediation effects of more than 100% indicate that the direct effect had a different sign than the total effect.
Thus, the microbiome-based prediction score for hippurate contained substantial information on markers of human health and mediated the corresponding effects of microbial alpha diversity.

2.4. Functional Annotations of an Independent Published Dataset Indicated a Direct Relationship between Microbial Diversity and Benzoate Metabolism

Above, we show the results of the analysed association patterns between urinary metabolites and microbiome diversity via statistical screening approaches. To strengthen our findings, we used an independent public dataset with metagenomic gut microbiome characterisations from Yachida et al. (n = 616) [39] (see Table S7 for sample characteristics) to analyse the relationship between microbial benzoate metabolism and microbial diversity. Note that the study of Yachida et al. [39] was a case-control study investigating colorectal cancer (CRC).
For functional annotations, we mapped the species abundances as reported in the supplement of Yachida et al. [39] onto AGORA2 [40]. AGORA2 is a resource of over 7000 genome-scale reconstructions of gut microbes that were semi-automatically generated and manually curated to match the experimental behaviour of the microbes. We then calculated reaction abundances for all microbial reactions that either produced or degraded benzoate (Figure 3A). Benzoate is the precursor of the human metabolite hippurate, which is predominantly produced in the liver through glycine conjugation (Figure 3A). Note that reaction abundances are not equivalent to gene abundances, as several genes may facilitate the same reactions and one gene may catalyse different reactions.
Then, we calculated the association between microbial diversity measured through the Shannon entropy and the reaction abundances of reactions either degrading, producing, or transporting benzoate (Figure 4A,B). Surprisingly, in all cases, we found a negative association between reaction abundances and Shannon entropy, which was most pronounced for the benzoate exchange reaction. In the sensitivity analyses, we tested whether these findings held up in both study groups of the Yachida et al. study (healthy controls and CRC patients), finding that this is a robust result regarding health and at least one disease known to alter the gut microbiome profoundly.
Furthermore, we utilised AGORA2 to retrieve the three biomarkers (benzoate, formate, and 4-hydroxyphenylacetate) of microbial diversity from the list of strains included in AGORA2 with the capacity to produce them (Supplementary Table S8). The results per phyla for benzoate are given in Figure 3B. Notably, benzoate production capabilities are distributed across several of the main phyla in the gut, including Actinobacteria, Bacteroidetes, Proteobacteria, and Firmicutes.
In conclusion, our functional analyses of the published data from Yachida et al. indicated a direct relationship between microbial diversity and benzoate metabolisation. However, the indicated relationship was negative with highly diverse communities showing lower reaction abundances of reactions transporting, degrading, or producing benzoate.

3. Discussion

Although immense research efforts have been undertaken, we are still at the beginning of understanding the multitude of effects the microbiome has on human health and disease. One key to understanding the relationship between the host and microbiome is to unravel the metabolic interplay between human metabolism and microbiome metabolism. This study contributed to the latter by investigating the urine metabolome via NMR technology [41] in conjunction with 16S rRNA gene sequencing of stool samples stemming from a large population-based study.
To date, a wide range of studies associating the urine metabolome with the gut microbiome have been performed, revealing a rich association pattern between the two types of omics data in various setups [42,43,44,45,46,47,48]. However, only a few studies have so far utilised NMR technology, and those studies were small studies relying on small-to-medium sample sizes [49,50,51,52]. Our present work was, to the best of our knowledge, the largest study so far that made use of NMR metabolomics to reveal the association patterns between the urine metabolome and the gut microbiome.
As NMR technology allows for quantifications of a wide range of microbiome-related metabolites, such as dimethylamine, formate, methanol, and acetate, which are normally not included in mass-spectrometric analyses, we could deliver novel insights while replicating earlier results, in particular the associations of urinary hippurate with measures of alpha diversity [34,51,52].
With formate and 4-hydroxyphenylacetate, we identified two novel urinary markers of microbial alpha diversity. Formate is a known microbial fermentation product of a wide range of microbes across various phyla [53,54] (Table S8). However, formate is also a product of human metabolism, e.g., as a by-product of multiple human pathways, such as tryptophan degradation, sterol metabolism, and one-carbon metabolism [53]. Thus, while it is plausible that urinary formate is a marker of microbial diversity, the relatively small effect size indicates that most of the variance in urinary formate concentration was due to host factors. Yet, urinary formate may be an interesting candidate for further research into host–microbiome metabolic interactions.
4-Hydroxyphenylacetate is an intermediate of tyrosine metabolism, both in humans and in microbes. However, the production of 4-hydroxyphenylacetate is a rather rare capacity among the microbes in AGORA2 (333 strains noted; in comparison, 6010 strains can produce formate and 1552 can produce benzoate; for complete lists, see Supplementary Table S8). Moreover, several important opportunistic pathogens are among the 4-hydroxyphenylacetate-producing species, such as Burkholderia cepacia [55], Acinetobacter baumannii [56], Clostridioides difficile [57], and Pseudomonas aeruginosa [58]. Notably, the association between microbiome diversity and urinary 4-hydroxyphenylacetate was negative, allowing for the speculation that 4-hydroxyphenylacetate is a metabolic marker of dysbiosis, in particular given the prominence of pathogens in the list of microbes capable of producing 4-hydroxyphenylacetate. Nevertheless, urinary 4-hydroxyphenylacetate may be largely determined by host factors, as we failed to predict a substantial amount of variance in urinary 4-hydroxyphenylacetate levels using genus abundances. However, it should be noted that small intestinal bacterial overgrowth in children has been associated with a marked increase in 4-hydroxyphenylacetate in urine [59].
A further result was the association between Escherichia/Shigella and indoxyl-3-sulfate, although sensitivity analyses indicated that this specific result needs careful interpretation. Escherichia/Shigella is a Gram-negative opportunistic pathogen [60], which is implicated in a variety of intraabdominal infections and is linked to gut microbiome instability [33]. It is also an important producer of pro-inflammatory lipopolysaccharides [61]. Indoxyl-3-sulfate is a metabolite produced by the liver through sulfation of the microbial fermentation product indole, which is produced during tryptophan degradation by various microbes, including members of the genera Escherichia and Shigella. Thus, while not reported previously, the association is plausible given the known metabolic traits of the associated genera. Importantly, indoxyl-3-sulfate is a known uremic toxin that accumulates in chronic kidney disease and is also said to be pro-inflammatory in general [62]. Thus, our results corroborated the role of Escherichia/Shigella, one of the most important opportunistic pathogens in the gut microbiome, in the production of indoxyl-3-sulfate with potential adverse clinical implications.
The main results of our study, however, were with respect to the metabolite hippurate. This metabolite was found to have a strong positive association with microbial alpha diversity, as well as a rich association pattern with various genera. These results provide an additional validation of earlier studies [34,51,52]. Thus, cumulative evidence points towards hippurate being a marker of alpha diversity and gut microbiome metabolism. Moreover, hippurate was also shown to be associated with metabolic health [63,64], and certain experimental evidence exists showing beneficial effects in mice models of diabetes [52]. Hippurate is produced via glycine conjugation of benzoate (Figure 3B), predominantly in the liver, although the pathway was also reported for the kidney [65]. Benzoate can either be directly ingested through the consumption of benzoate-containing food, such as berries, seafood, and dairy products, or produced by the microbiome by metabolising other phenols [66]. However, benzoate is also degraded by microbes in various pathways, including pathways that feed benzoate into central carbon metabolism [61]. Notably, the pathways for benzoate production and benzoate degradation co-occur in microbes [38]. Previous studies have shown that the gut microbiome has a clear effect on human benzoate pools [61]. Thus, the question arises whether the microbiome is a net consumer or producer of benzoate. In this respect, the analyses of the metagenomic data from Yachida et al. [39], which we functionally annotated using the AGORA2 resource [38], pointed towards the gut microbiome being a net consumer of benzoate. As our analyses demonstrated, highly diverse communities had lower reaction abundances of benzoate-transporting reactions. Thus, it follows plausibly that the higher the diversity, the lower the capacity to metabolise benzoate. Less benzoate consumption by the gut microbiome could then result in more benzoate being transported through the portal vein to the liver for hippurate formation. An important observation in this context is that benzoate exhibits antimicrobial activity such that high benzoate diets may lead to microbiomes with high-benzoate-metabolising capacities [66]. Although this hypothesis needs further validation, it would have consequences for interventional studies aiming at promoting human health through the microbiome. Our analyses would indicate that reducing microbial benzoate metabolism could be beneficial through the effects of hippurate.
Finally, our results indicated that the microbiome-related metabolic traits associated with urinary hippurate were of interest for human health, as they were robustly associated with markers of inflammation, liver injury, and metabolic health. As such, they fully mediated the effects of microbial diversity on the same markers, indicating that microbial benzoate metabolism is one of the metabolic functions underlying the rich associations of microbial diversity to human health.
Although our study successfully replicated known findings while expanding the scope of knowledge regarding metabolome–microbiome relations, certain strengths and limitations have to be kept in mind. First, while our study was of reasonable size, the statistical power of finding associations was low for rare genera and metabolites with missing quantifications. Moreover, since the human microbiome shows strong regional differences, generalisations of our results towards other regions are not guaranteed, and an independent validation investigating the herein-reported biomarkers is needed [67,68]. Moreover, due to the usage of 16S rRNA gene sequencing, the taxonomic resolution was limited to the genus level in most cases. Therefore, it was not possible to investigate species or strain-specific microbe–metabolite associations. Importantly, our analyses were cross-sectional in design, and while we could adjust for a wide range of covariates relevant to metabolomics and microbiome data, there was still the risk of unmeasured confounding. Moreover, our analyses of metabolite–microbiome relations were performed on a predominantly healthy population, which limits the ability to detect associations with manifest diseases. As our study was a fasting study, certain aspects of diet–microbiome-host interactions could not be investigated. In relation to microbial benzoate metabolism, it can be conceived that our current knowledge is still incomplete [61]. Thus, the hypothesis regarding the relationship between benzoate metabolism and alpha diversity based on functional annotations representing our current knowledge base needs further validation, preferably in an experimental set-up.
In conclusion, our study demonstrated the usefulness of NMR urine metabolomics in the assessment of host–microbiome interactions in the domain of metabolism. As NMR is inherently robust, reproducible, non-destructive, and involves minimal sample preparation, we believe NMR metabolomics will have a future in the emerging field of microbiome–host interaction in health and disease.

4. Materials and Methods

4.1. Study Population

Data of SHIP-TREND (recruited 2008–2012, n = 4420) was used, consisting of individuals from the region of Pomerania, North-Eastern Germany. The SHIP cohorts were designed to investigate the prevalence and incidence of clinical and subclinical phenotypes and their risk factors. For detailed information on the design, biomaterials, and available data, see [35].
In the presented study, we focused on the sub-cohort of SHIP-TREND with available 16S rRNA gene-sequencing data of stool samples (n = 3637) and the sub-cohort with NMR urine metabolomics (n = 996). Both 16S rRNA gene sequencing data and NMR metabolomics were available for 951 individuals. Due to the design of the SHIP-TREND study, individuals with known diabetes were not characterised using NMR metabolomics. This was the only exclusion criteria. However, as a consequence, the SHIP-TREND NMR sample was healthier than the average population, as diabetes-related comorbidities and traits were underrepresented.
The institutional review board of the University of Greifswald approved the survey and methods of the SHIP studies and all analyses followed the Declaration of Helsinki. Written informed consent was provided by all participants.

4.2. Assays and Phenotypes

The medical history and sociodemographic factors of each participant were evaluated using a computer-assisted face-to-face interview. Next, participants underwent extensive medical examinations, including measurements of the waist circumference, body height weight, and blood pressure. Waist circumference was used as an indicator of abdominal fat in all analyses. Participants were asked to bring their medication prescription sheets or packing containers of all medication they have been taken within the last seven days. By counting the number of alcoholic beverages that had been consumed on average per day over the last 30 days, alcohol consumption was estimated.
Blood samples were drawn between 07:00 a.m. and 12:30 p.m. Fasting time was assessed by asking for the last time the participants ate or drank beverages other than water. In SHIP-TREND, participants were explicitly asked not to eat or drink before blood sampling. Fasting time was an average of 13:28 h, standard deviation (SD) = 1:32 h. Blood samples were taken from the cubital vein and analysed directly or stored at −80 °C in the Integrated Research Biobank of the University Medicine Greifswald [69]. Spontaneous urine specimens were collected and immediately stored at −80 °C. White blood cell count (WBC), red blood cell count (RBC), and thrombocytes count (PLT) were measured on the XE 5000 from Sysmex (Sysmex Deutschland GmbH, Norderstedt, Germany). Glycated hemoglobin (HbA1c) concentrations were determined using high-performance liquid chromatography (Bio-Rad Diamat, Munich, Germany). Highly sensitive CRP, triglycerides, total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), GGT, ASAT, ALAT, and creatinine (Jaffé) were determined on the Dimension VISTA 1500 according to the recommendations of the manufacturer (Siemens Healthcare Diagnostics GmbH, Eschborn, Germany). From the serum creatinine, the estimated glomerular filtration rate (eGFR) was calculated using the Modification of Diet in Renal Disease equation, as described before [35]. All requirements of the corresponding guideline according to quality specifications were at least fulfilled. Plasma fibrinogen concentrations were assayed according to Clauss (BCS, Siemens Healthcare Diagnostics GmbH; Eschborn). Exocrine pancreatic function was determined using a pancreatic elastase ELISA (BIOSERV Diagnostics, Greifswald, Germany) based on faecal samples, as described before.

4.3. 16S rRNA Gene Sequencing and Taxonomic Assignments

Stool sample sequencing followed established procedures, as outlined before in detail [70]. Briefly, study participants collected faecal samples at home and stored them in a tube containing a stabilising DNA buffer. Next, faecal samples were transported to our laboratory. Then, DNA from the faecal samples was extracted (PSP Spin Stool DNA Kit; Stratec Biomedical AG, Birkenfeld, Germany) and stored at −20 °C until analysis using 16S rRNA gene sequencing of the V1–V2 region utilising a MiSeq platform (Illumina, San Diego, CA, USA). For taxonomic assignments, DADA2 16 (V.1.10) was employed for amplicon data processing, enabling single-nucleotide resolutions of amplicons. All samples were normalised to 10,000 16S rRNA gene read counts.

4.4. NMR Measurements in SHIP-TREND

After thawing, urine specimens were centrifuged for 5 min at 3000× g, and the supernatant was used for spectroscopic analysis. To this purpose, we mixed 450 µL of urine with 50 µL of phosphate buffer to stabilise the urinary pH at 7.0 (±0.35). The buffer was prepared with D2O and contained sodium 3-trimethylsilyl-(2,2,3,3-D4)-1-propionate (TSP) as a reference. Spectra were recorded on a Bruker DRX-400 NMR spectrometer (Bruker BioSpin GmbH, Ettlingen, Germany) at a 1H frequency of 400.13 MHz. The instrument was equipped with a 4 mm selective inverse flow probe (FISEI, 120 µL active volume) with a z-gradient. Specimens were automatically delivered to the spectrometer via flow injection. The acquisition temperature was set to 300 K. A standard one-dimensional 1H-NMR pulse sequence with suppression of the water peak (NOESYPREAST) was used: RD − P(90°) − 4 µsec − P(90°) − tm − P(90°) − acquisition of the free induction decay (FID). The non-selective 90° hard pulse P(90°) was adjusted to 9.4 µsec. The relaxation delay (RD), mixing time (tm), and acquisition time were set to 4 s, 100 msec, and 3.96 s, respectively, resulting in a total recycle time of ~8.0 s. Low-power continuous-wave irradiation on the water resonance at a field strength of ~25 Hz was applied during RD and tm for pre-saturation. After the application of 4 dummy scans, 64 FIDs were collected into 65,536 (64 K) complex data points using a spectral width of 20.689 parts per million (ppm). FIDs were multiplied with an exponential function corresponding to a line broadening of 0.3 Hz before Fourier transformation. Spectra were processed within TOP-SPIN 1.3 (Bruker BioSpin, Billerica, MA, USA).
The NMR spectra were analysed by Chenomx Inc. (Edmonton, AB, Canada) and manually annotated by spectral pattern matching using the Chenomx Worksuite 7.0 (Chenomx Inc., Edmonton, Canada) to deduce absolute urinary metabolite concentrations. The method of Chenomx’s patented NMR-based metabolomics platform is based on software and an extensive reference database that contains a compound list of over 300 metabolite models and was used for the analysis. Manual phasing, baseline correction, and targeted metabolic profiling using the Chenomx library were performed. This resulted in a list of 60 identified metabolites (including creatinine, which was used for normalisation) and their concentrations in millimoles per liter (mM). The spectrum regions of water (δ = 4.6–5.0) and the regions below δ = 0.0 and above δ = 10.0 were removed from the analysis for all groups in order to prevent variation in each sample. Each NMR variable was normalised to the total area in order to allow for a spectrum-to-spectrum comparison.

4.5. Data Normalisation and Outlier Detection

For statistical analyses, regression-based normalisation was performed on the log-transformed urinary concentrations based on probabilistic quotient normalisation using restricted cubic splines, as described in [36]. Regression-based normalisation leads to slightly higher statistical power, and it accounts for potential metabolite-specific dilution–concentration dependencies [36]. Details on regression-based normalisation can be found in Hertel et al. [36]. After the regression-based normalisation, outliers based on the 4-standard-deviation rule were excluded.

4.6. Statistical Analyses in SHIP-TREND

For descriptive statistics, metric variables were described using means and standard deviations, while categorical and ordinal variables were described using proportions. All reported p-values were two-tailed, and multiple testing correction was performed via applying the false discovery rate (FDR). All major routes of analysis were performed through multivariable regressions from the class of generalised linear models. All regression analyses included a basic set of covariates if not specified otherwise, consisting of age, sex, age–sex interaction terms, and waist circumference. Furthermore, we included smoking, hypertonia, HbA1c, years of education, the eGFR, urinary pH, and alcohol intake as covariates for the full adjustment. Age, waist circumference, and eGFR were introduced as restricted cubic splines using four knots at the 5% percentile, the 33% percentile, the 66% percentile, and the 95% percentile of the respective distributions. Nonlinear modelling for these variables was chosen, as previous analyses indicated nonlinearity in respect to basic covariates [71].
In the first set of regression models, we screened urinary metabolite concentrations on associations with genus abundances. We included only metabolites and genera with at least 50% non-zero measurements, resulting in 1681 metabolite–genus combinations tested. This rather strict criterion was applied to filter out metabolite–genus combinations with low case numbers and thus low statistical power. The regressions then included normalised metabolite concentration as a response variable, while the genus abundance was utilised as a predictor of interest adjusting for the set of basic covariates. Lastly, we used the 16S rRNA gene-sequencing batch variable as a random effect variable to account for batch effects. The herein identified metabolite–genus pairs with FDR < 0.05 were then subject to sensitivity analyses using the full list of covariates. Furthermore, we utilised non-parametric bootstrapping (2000 replications) for recalculating p-values, ensuring that our results are not biased by certain distributional features of the microbial abundances. In the main text, we report the fully adjusted effects, whereas Supplementary Tables S2 and S3 contain the summary statistics for all regressions performed.
In the second set of analogous regression models, we exchanged the genus abundance with a binary variable indicating genus presence or absence. Here, we included all genera that were found in at 20% of the samples and maximally in 80% of the samples. This led to 1848 genus–metabolite pairs. As before, these combinations were screened on significant metabolite–genus associations using the basic set of covariates, and significant associations were then re-tested using the full adjustment.
Third, to test the associations between the urine metabolome and metrics of gut microbiome alpha diversity, we ran mixed-effect regression models as before. However, now the predictor of interest was either the Shannon entropy or the species richness, as calculated via the Vegan R package [72]. Multivariable linear mixed regressions were performed for the 42 metabolites having more than 50% non-zero measurements testing on associations between metabolites and alpha diversity measures, both for the basic adjustment and the full list of covariates. After graphical examination of the results indicating the presence of nonlinear diversity, we furthermore checked on nonlinear associations via restricted cubic splines as before.
Fourth, we explored the predictive value of the microbiome for each of the 42 metabolites that was measured in more than 50% of the cases. To this end, we used an elastic net regression methodology, which is a machine learning approach from the class of penalised regressions based on the objective function
Q = 1 2 n ∑ i = 1 n ( y i − β 0 − x i β ′ ) 2 + λ ∑ j = 1 p ( 1 − α 2 β j   2 + α | β j | )
where n is the sample size, y is the response variable (the normalised metabolite concentration), x is the vector of predictors (genus abundances), p is the number of predictors, β is the regression coefficients, and λ and α are the penalty parameters. Then, a grid search over various λ and α values was performed using internal cross-validation to find the parameters with the lowest prediction error. Note that if α = 1, the elastic net is equivalent to the Lasso, while with α = 0, it becomes equivalent to the Ridge regression. We performed a grid search on α = {0, 0.1, …, 1} and using the default parameters for λ in the STATA 17/mp elasticnet function. We utilised 10-fold cross-validation and utilised the out-of-sample R-squared value as an indicator of model fit. Full results for all 42 metabolites can be found in Supplementary Table S6.
As only the microbiome-based prediction score for urinary hippurate explained the substantial proportions of variance, we explored the information encoded in this prediction score further. To this end, we utilised linear mixed-effect regressions with either the normalised urinary hippurate concentrations, the Shannon entropy, or the microbiome-based hippurate prediction score as a response variable, the full set of covariates, and one additional clinical parameter as a predictor of interest. In this way, we screened inflammation markers (high sensitivity CRP, fibrinogen, white blood cell counts), markers of metabolic health (triglycerides, ratio of total cholesterol and HDL cholesterol, HbA1c), and markers of liver injury (GGT, ASAT, ALAT) on the association with urinary hippurate concentration, the Shannon entropy, and the microbiome-based hippurate prediction score. Finally, we performed mediation analyses following the procedure of Hicks et al. [73] to test whether the microbiome-based hippurate prediction score mediates the effects of microbiome diversity on markers of human health.

4.7. Statistical Analyses on Yachida et al.’s Metagenome Data

To enlighten the relation between microbial benzoate metabolism and microbial alpha diversity further, we went on to analyse a public metagenomics dataset (n = 616, see Yachida et al. [39] for details) with functional annotations based on the AGORA2 platform [40]. In brief, species abundances were retrieved via the MetaPhlaN2 pipeline and then mapped onto AGORA2. AGORA2 is a collection of semi-automatically derived and manually curated genome-scale reconstructions of microbes that was curated against known experimental data reflecting the true metabolic capacities of the microbes [40]. We utilised the genome-scale reconstructions to analyse which microbes were able to secrete benzoate, formate, and 4-hydroxyphenylacetate. Moreover, we calculated reaction abundances for all reactions involving benzoate based on the AGORA2 functional annotations using established pipelines [74] for all metagenomes measured in Yachida et al. [39]. Then, we regressed the reaction abundances via fractional regressions on the Shannon entropy including age, sex, body mass index, and the case–control status as covariates.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/metabo12040308/s1, Table S1: Descriptive statistics on urinary metabolites that were included and excluded from the analyses; Table S2: Association analyses between genus presences and urinary metabolite concentrations; Table S3: Association analyses between genus presences and urinary metabolite concentrations; Table S4: Summary statistics for metabolite Shannon diversity associations; Table S5: Summary statistics for metabolite species richness associations; Table S6: Overview on performances on predicting metabolite concentrations from genus abundances via elastic net regression methodologies; Table S7: Sample characteristics of the utilised metagenomic data from Yachida et al. (2019); Table S8: Lists of strains from AGORA2 producing benzoate, formate, and 4-hydroxyphenylacetate.

Author Contributions

Conceptualisation, J.H. and F.F.; methodology, J.H., F.F. and A.H.; formal analysis, J.H., F.F. and A.H.; resources, I.T., H.V., N.F., M.N., F.U.W., H.-J.G., M.M.L. and U.V.; data curation, N.F., K.B., A.-K.H., A.P., F.U.W., F.F., A.F., C.B. and M.R.; writing—original draft preparation, J.H., F.F. and D.F.; writing—review and editing, all authors; visualisation, J.H.; supervision, N.F. All authors have read and agreed to the published version of the manuscript.

Funding

SHIP is part of the Research Network Community Medicine of the University Medicine Greifswald, which is supported by the German Federal State of Mecklenburg-West Pomerania. This study was, in part, funded by a grant from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 757922) to I.T.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of the University of Greifswald.

Informed Consent Statement

Written informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are not publicly available due to the data use policy of SHIP. All SHIP data can be requested free of charge from the SHIP transfer office (https://www.fvcm.med.uni-greifswald.de/cm_antrag/index.php accessed on 29 March 2022).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. H.J.G. has received travel grants and speakers honoraria from Fresenius Medical Care, Neuraxpharm, Servier, and Janssen Cilag. A.P. has received lecturer honoraria from Technopath Clinical Diagnostics, Tosoh Bioscience, Roche Diagnostics, Beckmann Coulter GmbH, and Radiometer. M.N. has received travel grants by German Medical Association, German Centre for Cardiovascular Research, German Society for Clinical Chemistry and Laboratory Medicine, German National Cohort, German Research Foundation, Deutsche Akkreditierungsstelle, Sysmex, MDI Limbach, medpoint GmbH and Diasys; speaker’s honoraria from Novartis Pharma, Radiometer, AstraZeneca, Technopath Clinical Diagnostics, Sysmex, MDI Limbach, and medpoint Medizinkommunikations GmbH; and research funding from Aerocom GmbH and Profil Institut für Stoffwechselforschung GmbH.

References

  1. Nicholson, J.K.; Holmes, E.; Kinross, J.; Burcelin, R.; Gibson, G.; Jia, W.; Pettersson, S. Host-Gut Microbiota Metabolic Interactions. Science 2012, 336, 1262–1267. [Google Scholar] [CrossRef] [Green Version]
  2. Cho, I.; Blaser, M.J. The Human Microbiome: At the Interface of Health and Disease. Nat. Rev. Genet. 2012, 13, 260–270. [Google Scholar] [CrossRef] [Green Version]
  3. Lynch, S.V.; Pedersen, O. The Human Intestinal Microbiome in Health and Disease. N Engl. J. Med. 2016, 375, 2369–2379. [Google Scholar] [CrossRef] [Green Version]
  4. Hollister, E.B.; Gao, C.; Versalovic, J. Compositional and Functional Features of the Gastrointestinal Microbiome and Their Effects on Human Health. Gastroenterology 2014, 146, 1449–1458. [Google Scholar] [CrossRef] [Green Version]
  5. Visconti, A.; Le Roy, C.I.; Rosa, F.; Rossi, N.; Martin, T.C.; Mohney, R.P.; Li, W.; de Rinaldis, E.; Bell, J.T.; Venter, J.C.; et al. Interplay between the Human Gut Microbiome and Host Metabolism. Nat. Commun. 2019, 10, 4505. [Google Scholar] [CrossRef] [Green Version]
  6. LeBlanc, J.G.; Chain, F.; Martín, R.; Bermúdez-Humarán, L.G.; Courau, S.; Langella, P. Beneficial Effects on Host Energy Metabolism of Short-Chain Fatty Acids and Vitamins Produced by Commensal and Probiotic Bacteria. Microb. Cell Fact. 2017, 16, 79. [Google Scholar] [CrossRef] [Green Version]
  7. Puertollano, E.; Kolida, S.; Yaqoob, P. Biological Significance of Short-Chain Fatty Acid Metabolism by the Intestinal Microbiome. Curr. Opin. Clin. Nutr. Metab. Care 2014, 17, 139–144. [Google Scholar] [CrossRef]
  8. Tan, J.; McKenzie, C.; Potamitis, M.; Thorburn, A.N.; Mackay, C.R.; Macia, L. The Role of Short-Chain Fatty Acids in Health and Disease. In Advances in Immunology; Elsevier: Amsterdam, The Netherlands, 2014; Volume 121, pp. 91–119. ISBN 978-0-12-800100-4. [Google Scholar]
  9. Buret, A.G.; Allain, T.; Motta, J.-P.; Wallace, J.L. Effects of Hydrogen Sulfide on the Microbiome: From Toxicity to Therapy. Antioxid. Redox Signal. 2021, 36, 211–219. [Google Scholar] [CrossRef]
  10. Blachier, F.; Andriamihaja, M.; Larraufie, P.; Ahn, E.; Lan, A.; Kim, E. Production of Hydrogen Sulfide by the Intestinal Microbiota and Epithelial Cells and Consequences for the Colonic and Rectal Mucosa. Am. J. Physiol. Gastrointest. Liver Physiol. 2021, 320, G125–G135. [Google Scholar] [CrossRef]
  11. Ridlon, J.M.; Kang, D.J.; Hylemon, P.B.; Bajaj, J.S. Bile Acids and the Gut Microbiome. Curr. Opin. Gastroenterol. 2014, 30, 332–338. [Google Scholar] [CrossRef] [Green Version]
  12. Chen, J.; Thomsen, M.; Vitetta, L. Interaction of Gut Microbiota with Dysregulation of Bile Acids in the Pathogenesis of Nonalcoholic Fatty Liver Disease and Potential Therapeutic Implications of Probiotics. J. Cell Biochem. 2019, 120, 2713–2720. [Google Scholar] [CrossRef]
  13. Wahlström, A.; Sayin, S.I.; Marschall, H.-U.; Bäckhed, F. Intestinal Crosstalk between Bile Acids and Microbiota and Its Impact on Host Metabolism. Cell Metab. 2016, 24, 41–50. [Google Scholar] [CrossRef] [Green Version]
  14. Tang, W.H.W.; Hazen, S.L. Microbiome, Trimethylamine N-Oxide, and Cardiometabolic Disease. Transl. Res. 2017, 179, 108–115. [Google Scholar] [CrossRef] [Green Version]
  15. Zeisel, S.H.; Warrier, M. Trimethylamine N-Oxide, the Microbiome, and Heart and Kidney Disease. Annu. Rev. Nutr. 2017, 37, 157–181. [Google Scholar] [CrossRef]
  16. The Human Microbiome Project Consortium Structure, Function and Diversity of the Healthy Human Microbiome. Nature 2012, 486, 207–214. [CrossRef] [Green Version]
  17. Sandrini, S.; Aldriwesh, M.; Alruways, M.; Freestone, P. Microbial Endocrinology: Host–Bacteria Communication within the Gut Microbiome. J. Endocrinol. 2015, 225, R21–R34. [Google Scholar] [CrossRef] [Green Version]
  18. Bonder, M.J.; Kurilshikov, A.; Tigchelaar, E.F.; Mujagic, Z.; Imhann, F.; Vila, A.V.; Deelen, P.; Vatanen, T.; Schirmer, M.; Smeekens, S.P.; et al. The Effect of Host Genetics on the Gut Microbiome. Nat. Genet. 2016, 48, 1407–1412. [Google Scholar] [CrossRef]
  19. Goodrich, J.K.; Waters, J.L.; Poole, A.C.; Sutter, J.L.; Koren, O.; Blekhman, R.; Beaumont, M.; Van Treuren, W.; Knight, R.; Bell, J.T.; et al. Human Genetics Shape the Gut Microbiome. Cell 2014, 159, 789–799. [Google Scholar] [CrossRef] [Green Version]
  20. Rühlemann, M.C.; Hermes, B.M.; Bang, C.; Doms, S.; Moitinho-Silva, L.; Thingholm, L.B.; Frost, F.; Degenhardt, F.; Wittig, M.; Kässens, J.; et al. Genome-Wide Association Study in 8956 German Individuals Identifies Influence of ABO Histo-Blood Groups on Gut Microbiome. Nat. Genet. 2021, 53, 147–155. [Google Scholar] [CrossRef]
  21. Kurilshikov, A.; Medina-Gomez, C.; Bacigalupe, R.; Radjabzadeh, D.; Wang, J.; Demirkan, A.; Le Roy, C.I.; Raygoza Garay, J.A.; Finnicum, C.T.; Liu, X.; et al. Large-Scale Association Analyses Identify Host Factors Influencing Human Gut Microbiome Composition. Nat. Genet. 2021, 53, 156–165. [Google Scholar] [CrossRef]
  22. Hasan, N.; Yang, H. Factors Affecting the Composition of the Gut Microbiota, and Its Modulation. PeerJ 2019, 7, e7502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Aleman, F.D.D.; Valenzano, D.R. Microbiome Evolution during Host Aging. PLoS Pathog. 2019, 15, e1007727. [Google Scholar] [CrossRef] [PubMed]
  24. Prehn-Kristensen, A.; Zimmermann, A.; Tittmann, L.; Lieb, W.; Schreiber, S.; Baving, L.; Fischer, A. Reduced Microbiome Alpha Diversity in Young Patients with ADHD. PLoS ONE 2018, 13, e0200728. [Google Scholar] [CrossRef] [PubMed]
  25. Schirmer, M.; Franzosa, E.A.; Lloyd-Price, J.; McIver, L.J.; Schwager, R.; Poon, T.W.; Ananthakrishnan, A.N.; Andrews, E.; Barron, G.; Lake, K.; et al. Dynamics of Metatranscription in the Inflammatory Bowel Disease Gut Microbiome. Nat. Microbiol. 2018, 3, 337–346. [Google Scholar] [CrossRef]
  26. Hirano, A.; Umeno, J.; Okamoto, Y.; Shibata, H.; Ogura, Y.; Moriyama, T.; Torisu, T.; Fujioka, S.; Fuyuno, Y.; Kawarabayasi, Y.; et al. Comparison of the Microbial Community Structure between Inflamed and Non-Inflamed Sites in Patients with Ulcerative Colitis: Dysbiosis in Ulcerative Colitis. J. Gastroenterol. Hepatol. 2018, 33, 1590–1597. [Google Scholar] [CrossRef] [PubMed]
  27. Stanislawski, M.A.; Dabelea, D.; Lange, L.A.; Wagner, B.D.; Lozupone, C.A. Gut Microbiota Phenotypes of Obesity. npj Biofilms Microbiom. 2019, 5, 18. [Google Scholar] [CrossRef] [PubMed]
  28. Menni, C.; Zhu, J.; Le Roy, C.I.; Mompeo, O.; Young, K.; Rebholz, C.M.; Selvin, E.; North, K.E.; Mohney, R.P.; Bell, J.T.; et al. Serum Metabolites Reflecting Gut Microbiome Alpha Diversity Predict Type 2 Diabetes. Gut Microb. 2020, 11, 1632–1642. [Google Scholar] [CrossRef]
  29. Frost, F.; Weiss, F.U.; Sendler, M.; Kacprowski, T.; Rühlemann, M.; Bang, C.; Franke, A.; Völker, U.; Völzke, H.; Lamprecht, G.; et al. The Gut Microbiome in Patients With Chronic Pancreatitis Is Characterized by Significant Dysbiosis and Overgrowth by Opportunistic Pathogens. Clin. Trans. Gastroenterol. 2020, 11, e00232. [Google Scholar] [CrossRef]
  30. Pietzner, M.; Budde, K.; Rühlemann, M.; Völzke, H.; Homuth, G.; Weiss, F.U.; Lerch, M.M.; Frost, F. Exocrine Pancreatic Function Modulates Plasma Metabolites Through Changes in Gut Microbiota Composition. J. Clin. Endocrinol. Metab. 2021, 106, e2290–e2298. [Google Scholar] [CrossRef]
  31. Frost, F.; Kacprowski, T.; Rühlemann, M.; Weiss, S.; Bang, C.; Franke, A.; Pietzner, M.; Aghdassi, A.A.; Sendler, M.; Völker, U.; et al. Carrying Asymptomatic Gallstones Is Not Associated with Changes in Intestinal Microbiota Composition and Diversity but Cholecystectomy with Significant Dysbiosis. Sci. Rep. 2021, 11, 6677. [Google Scholar] [CrossRef]
  32. Frost, F.; Storck, L.J.; Kacprowski, T.; Gärtner, S.; Rühlemann, M.; Bang, C.; Franke, A.; Völker, U.; Aghdassi, A.A.; Steveling, A.; 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, e0219489. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Frost, F.; Kacprowski, T.; Rühlemann, M.; Pietzner, M.; Bang, C.; Franke, A.; Nauck, M.; Völker, U.; Völzke, H.; Dörr, M.; et al. Long-Term Instability of the Intestinal Microbiome Is Associated with Metabolic Liver Disease, Low Microbiota Diversity, Diabetes Mellitus and Impaired Exocrine Pancreatic Function. Gut 2021, 70, 522–530. [Google Scholar] [CrossRef] [PubMed]
  34. Wilmanski, T.; Rappaport, N.; Earls, J.C.; Magis, A.T.; Manor, O.; Lovejoy, J.; Omenn, G.S.; Hood, L.; Gibbons, S.M.; Price, N.D. Blood Metabolome Predicts Gut Microbiome α-Diversity in Humans. Nat. Biotechnol. 2019, 37, 1217–1228. [Google Scholar] [CrossRef]
  35. Volzke, H.; Alte, D.; Schmidt, C.O.; Radke, D.; Lorbeer, R.; Friedrich, N.; Aumann, N.; Lau, K.; Piontek, M.; Born, G.; et al. Cohort Profile: The Study of Health in Pomerania. Int. J. Epidemiol. 2011, 40, 294–307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Hertel, J.; Rotter, M.; Frenzel, S.; Zacharias, H.U.; Krumsiek, J.; Rathkolb, B.; Hrabe de Angelis, M.; Rabstein, S.; Pallapies, D.; Brüning, T.; et al. Dilution Correction for Dynamically Influenced Urinary Analyte Data. Anal. Chim. Acta 2018, 1032, 18–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Hertel, J.; Heinken, A.; Martinelli, F.; Thiele, I. Integration of Constraint-Based Modeling with Fecal Metabolomics Reveals Large Deleterious Effects of Fusobacterium Spp. on Community Butyrate Production. Gut Microb. 2021, 13, 1915673. [Google Scholar] [CrossRef]
  38. Hertel, J.; Harms, A.C.; Heinken, A.; Baldini, F.; Thinnes, C.C.; Glaab, E.; Vasco, D.A.; Pietzner, M.; Stewart, I.D.; Wareham, N.J.; et al. Integrated Analyses of Microbiome and Longitudinal Metabolome Data Reveal Microbial-Host Interactions on Sulfur Metabolism in Parkinson’s Disease. Cell Rep. 2019, 29, 1767–1777.e8. [Google Scholar] [CrossRef] [Green Version]
  39. Yachida, S.; Mizutani, S.; Shiroma, H.; Shiba, S.; Nakajima, T.; Sakamoto, T.; Watanabe, H.; Masuda, K.; Nishimoto, Y.; Kubo, M.; et al. Metagenomic and Metabolomic Analyses Reveal Distinct Stage-Specific Phenotypes of the Gut Microbiota in Colorectal Cancer. Nat. Med. 2019, 25, 968–976. [Google Scholar] [CrossRef]
  40. Heinken, A.; Acharya, G.; Ravcheev, D.A.; Hertel, J.; Nyga, M.; Okpala, O.E.; Hogan, M.; Magnúsdóttir, S.; Martinelli, F.; Preciat, G.; et al. AGORA2: Large Scale Reconstruction of the Microbiome Highlights Wide-Spread Drug-Metabolising Capacities. Systems Biology. 2020. Available online: https://www.biorxiv.org/content/10.1101/2020.11.09.375451v1 (accessed on 28 January 2022).
  41. Wishart, D.S. Quantitative Metabolomics Using NMR. TrAC Trends Anal. Chem. 2008, 27, 228–237. [Google Scholar] [CrossRef]
  42. Wu, S.; Wang, M.; Zhang, M.; He, J.-Q. Metabolomics and Microbiomes for Discovering Biomarkers of Antituberculosis Drugs-Induced Hepatotoxicity. Arch. Biochem. Biophys. 2022, 716, 109118. [Google Scholar] [CrossRef]
  43. Xu, H.; Tamrat, N.E.; Gao, J.; Xu, J.; Zhou, Y.; Zhang, S.; Chen, Z.; Shao, Y.; Ding, L.; Shen, B.; et al. Combined Signature of the Urinary Microbiome and Metabolome in Patients With Interstitial Cystitis. Front. Cell Infect. Microbiol. 2021, 11, 711746. [Google Scholar] [CrossRef] [PubMed]
  44. Yaskolka Meir, A.; Tuohy, K.; von Bergen, M.; Krajmalnik-Brown, R.; Heinig, U.; Zelicha, H.; Tsaban, G.; Rinott, E.; Kaplan, A.; Aharoni, A.; et al. The Metabolomic-Gut-Clinical Axis of Mankai Plant-Derived Dietary Polyphenols. Nutrients 2021, 13, 1866. [Google Scholar] [CrossRef]
  45. Laforest-Lapointe, I.; Becker, A.B.; Mandhane, P.J.; Turvey, S.E.; Moraes, T.J.; Sears, M.R.; Subbarao, P.; Sycuro, L.K.; Azad, M.B.; Arrieta, M.-C. Maternal Consumption of Artificially Sweetened Beverages during Pregnancy Is Associated with Infant Gut Microbiota and Metabolic Modifications and Increased Infant Body Mass Index. Gut Microb. 2021, 13, 1857513. [Google Scholar] [CrossRef] [PubMed]
  46. Wang, X.; Liu, H.; Li, Y.; Huang, S.; Zhang, L.; Cao, C.; Baker, P.N.; Tong, C.; Zheng, P.; Qi, H. Altered Gut Bacterial and Metabolic Signatures and Their Interaction in Gestational Diabetes Mellitus. Gut Microb. 2020, 12, 1840765. [Google Scholar] [CrossRef] [PubMed]
  47. Liu, T.; Gu, X.; Li, L.-X.; Li, M.; Li, B.; Cui, X.; Zuo, X.-L. Microbial and Metabolomic Profiles in Correlation with Depression and Anxiety Co-Morbidities in Diarrhoea-Predominant IBS Patients. BMC Microbiol. 2020, 20, 168. [Google Scholar] [CrossRef]
  48. Jeffery, I.B.; Das, A.; O’Herlihy, E.; Coughlan, S.; Cisek, K.; Moore, M.; Bradley, F.; Carty, T.; Pradhan, M.; Dwibedi, C.; et al. Differences in Fecal Microbiomes and Metabolomes of People With vs Without Irritable Bowel Syndrome and Bile Acid Malabsorption. Gastroenterology 2020, 158, 1016–1028.e8. [Google Scholar] [CrossRef]
  49. Dewulf, E.M.; Cani, P.D.; Claus, S.P.; Fuentes, S.; Puylaert, P.G.; Neyrinck, A.M.; Bindels, L.B.; de Vos, W.M.; Gibson, G.R.; Thissen, J.-P.; et al. Insight into the Prebiotic Concept: Lessons from an Exploratory, Double Blind Intervention Study with Inulin-Type Fructans in Obese Women. Gut 2013, 62, 1112–1121. [Google Scholar] [CrossRef]
  50. Prochazkova, M.; Budinska, E.; Kuzma, M.; Pelantova, H.; Hradecky, J.; Heczkova, M.; Daskova, N.; Bratova, M.; Modos, I.; Videnska, P.; et al. Vegan Diet Is Associated With Favorable Effects on the Metabolic Performance of Intestinal Microbiota: A Cross-Sectional Multi-Omics Study. Front. Nutr. 2021, 8, 783302. [Google Scholar] [CrossRef]
  51. West, K.A.; Kanu, C.; Maric, T.; McDonald, J.A.K.; Nicholson, J.K.; Li, J.V.; Johnson, M.R.; Holmes, E.; Savvidou, M.D. Longitudinal Metabolic and Gut Bacterial Profiling of Pregnant Women with Previous Bariatric Surgery. Gut 2020, 69, 1452–1459. [Google Scholar] [CrossRef]
  52. Wijeyesekera, A.; Wagner, J.; De Goffau, M.; Thurston, S.; Rodrigues Sabino, A.; Zaher, S.; White, D.; Ridout, J.; Peters, M.J.; Ramnarayan, P.; et al. Multi-Compartment Profiling of Bacterial and Host Metabolites Identifies Intestinal Dysbiosis and Its Functional Consequences in the Critically Ill Child. Crit. Care Med. 2019, 47, e727–e734. [Google Scholar] [CrossRef]
  53. Pietzke, M.; Meiser, J.; Vazquez, A. Formate Metabolism in Health and Disease. Mol. Metab. 2020, 33, 23–37. [Google Scholar] [CrossRef] [PubMed]
  54. Hughes, E.R.; Winter, M.G.; Duerkop, B.A.; Spiga, L.; Furtado de Carvalho, T.; Zhu, W.; Gillis, C.C.; Büttner, L.; Smoot, M.P.; Behrendt, C.L.; et al. Microbial Respiration and Formate Oxidation as Metabolic Signatures of Inflammation-Associated Dysbiosis. Cell Host Microb. 2017, 21, 208–219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. McClean, S.; Callaghan, M. Burkholderia Cepacia Complex: Epithelial Cell–Pathogen Confrontations and Potential for Therapeutic Intervention. J. Med. Microbiol. 2009, 58, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Howard, A.; O’Donoghue, M.; Feeney, A.; Sleator, R.D. Acinetobacter Baumannii: An Emerging Opportunistic Pathogen. Virulence 2012, 3, 243–250. [Google Scholar] [CrossRef] [PubMed]
  57. Neumann-Schaal, M.; Jahn, D.; Schmidt-Hohagen, K. Metabolism the Difficile Way: The Key to the Success of the Pathogen Clostridioides Difficile. Front. Microbiol. 2019, 10, 219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Stover, C.K.; Pham, X.Q.; Erwin, A.L.; Mizoguchi, S.D.; Warrener, P.; Hickey, M.J.; Brinkman, F.S.L.; Hufnagle, W.O.; Kowalik, D.J.; Lagrou, M.; et al. Complete Genome Sequence of Pseudomonas Aeruginosa PAO1, an Opportunistic Pathogen. Nature 2000, 406, 959–964. [Google Scholar] [CrossRef]
  59. McGrath, K.H.; Pitt, J.; Bines, J.E. Small Intestinal Bacterial Overgrowth in Children with Intestinal Failure on Home Parenteral Nutrition. JGH Open 2019, 3, 394–399. [Google Scholar] [CrossRef]
  60. Sims, G.E.; Kim, S.-H. Whole-Genome Phylogeny of Escherichia Coli/Shigella Group by Feature Frequency Profiles (FFPs). Proc. Natl. Acad. Sci. USA 2011, 108, 8329–8334. [Google Scholar] [CrossRef] [Green Version]
  61. Cani, P.D.; Amar, J.; Iglesias, M.A.; Poggi, M.; Knauf, C.; Bastelica, D.; Neyrinck, A.M.; Fava, F.; Tuohy, K.M.; Chabo, C.; et al. Metabolic Endotoxemia Initiates Obesity and Insulin Resistance. Diabetes 2007, 56, 1761–1772. [Google Scholar] [CrossRef] [Green Version]
  62. Cheng, T.-H.; Ma, M.-C.; Liao, M.-T.; Zheng, C.-M.; Lu, K.-C.; Liao, C.-H.; Hou, Y.-C.; Liu, W.-C.; Lu, C.-L. Indoxyl Sulfate, a Tubular Toxin, Contributes to the Development of Chronic Kidney Disease. Toxins 2020, 12, 684. [Google Scholar] [CrossRef]
  63. Brial, F.; Chilloux, J.; Nielsen, T.; Vieira-Silva, S.; Falony, G.; Andrikopoulos, P.; Olanipekun, M.; Hoyles, L.; Djouadi, F.; Neves, A.L.; et al. Human and Preclinical Studies of the Host–Gut Microbiome Co-Metabolite Hippurate as a Marker and Mediator of Metabolic Health. Gut 2021, 70, 2105–2114. [Google Scholar] [CrossRef] [PubMed]
  64. Williams, H.R.; Cox, I.J.; Walker, D.G.; Cobbold, J.F.; Taylor-Robinson, S.D.; Marshall, S.E.; Orchard, T.R. Differences in Gut Microbial Metabolism Are Responsible for Reduced Hippurate Synthesis in Crohn’s Disease. BMC Gastroenterol. 2010, 10, 108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Pallister, T.; Jackson, M.A.; Martin, T.C.; Zierer, J.; Jennings, A.; Mohney, R.P.; MacGregor, A.; Steves, C.J.; Cassidy, A.; Spector, T.D.; et al. Hippurate as a Metabolomic Marker of Gut Microbiome Diversity: Modulation by Diet and Relationship to Metabolic Syndrome. Sci. Rep. 2017, 7, 13670. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Yadav, M.; Lomash, A.; Kapoor, S.; Pandey, R.; Chauhan, N.S. Mapping of the Benzoate Metabolism by Human Gut Microbiome Indicates Food-Derived Metagenome Evolution. Sci. Rep. 2021, 11, 5561. [Google Scholar] [CrossRef]
  67. Gupta, V.K.; Paul, S.; Dutta, C. Geography, Ethnicity or Subsistence-Specific Variations in Human Microbiome Composition and Diversity. Front. Microbiol. 2017, 8, 1162. [Google Scholar] [CrossRef] [Green Version]
  68. Nishijima, S.; Suda, W.; Oshima, K.; Kim, S.-W.; Hirose, Y.; Morita, H.; Hattori, M. The Gut Microbiome of Healthy Japanese and Its Microbial and Functional Uniqueness. DNA Res. 2016, 23, 125–133. [Google Scholar] [CrossRef] [Green Version]
  69. Winter, T.; Friedrich, N.; Lamp, S.; Schäfer, C.; Schattschneider, M.; Bollmann, S.; Brümmer, D.; Riemann, K.; Petersmann, A.; Nauck, M. The Integrated Research Biobank of the University Medicine Greifswald. Open J. Bioresour. 2020, 7, 2. [Google Scholar] [CrossRef] [Green Version]
  70. Frost, F.; Kacprowski, T.; Rühlemann, M.; Bülow, R.; Kühn, J.-P.; Franke, A.; Heinsen, F.-A.; Pietzner, M.; Nauck, M.; Völker, U.; et al. Impaired Exocrine Pancreatic Function Associates With Changes in Intestinal Microbiota Composition and Diversity. Gastroenterology 2019, 156, 1010–1015. [Google Scholar] [CrossRef]
  71. Hertel, J.; Friedrich, N.; Wittfeld, K.; Pietzner, M.; Budde, K.; Van der Auwera, S.; Lohmann, T.; Teumer, A.; Völzke, H.; Nauck, M.; et al. Measuring Biological Age via Metabonomics: The Metabolic Age Score. J. Proteome Res. 2016, 15, 400–410. [Google Scholar] [CrossRef]
  72. Dixon, P. VEGAN, a Package of R Functions for Community Ecology. J. Veg. Sci. 2003, 14, 927–930. [Google Scholar] [CrossRef]
  73. Hicks, R.; Tingley, D. Causal Mediation Analysis. Stata J. 2011, 11, 605–619. [Google Scholar] [CrossRef] [Green Version]
  74. Heinken, A.; Magnúsdóttir, S.; Fleming, R.M.T.; Thiele, I. DEMETER: Efficient Simultaneous Curation of Genome-Scale Reconstructions Guided by Experimental Data and Refined Gene Annotations. Bioinformatics 2021, 37, 3974–3975. [Google Scholar] [CrossRef] [PubMed]
Metabolites 12 00308 g001 550
Figure 1. Overview of the significant metabolite–genus associations. (A) Boxplots for genus presence metabolite association with a false discovery rate < 0.05. The Y-axis denotes the log-transformed and regression-normalised urinary concentrations. (B) Table displaying the genus presence metabolite associations with FDR < 0.05. (C) Scatter plots of selected genus abundances against urinary metabolite levels after log-transformation and regression-based normalisation. The red line indicates the linear regression line, where the dashed red lines display the 95% confidence interval. (D) Table displaying the genus abundance metabolite associations with FDR < 0.05. FDR—false discovery rate; b—unstandardised regression coefficients; 95% CI—95% confidence intervals.
Figure 1. Overview of the significant metabolite–genus associations. (A) Boxplots for genus presence metabolite association with a false discovery rate < 0.05. The Y-axis denotes the log-transformed and regression-normalised urinary concentrations. (B) Table displaying the genus presence metabolite associations with FDR < 0.05. (C) Scatter plots of selected genus abundances against urinary metabolite levels after log-transformation and regression-based normalisation. The red line indicates the linear regression line, where the dashed red lines display the 95% confidence interval. (D) Table displaying the genus abundance metabolite associations with FDR < 0.05. FDR—false discovery rate; b—unstandardised regression coefficients; 95% CI—95% confidence intervals.
Metabolites 12 00308 g001
Metabolites 12 00308 g002 550
Figure 2. Overview of microbiome alpha diversity–metabolite associations. (A) Scatter plots of Shannon diversity or species richness against urinary metabolite levels after log-transformation and regression-based normalisation. A red line indicates the linear regression line, where the dashed red lines display the 95% confidence interval. All shown associations have an FDR < 0.05. (B) Table displaying the alpha diversity metabolite associations with an FDR < 0.05. (C) Scatter plots of Shannon diversity or urinary hippurate concentrations (log-transformed and regression-based normalised) against the microbiome-based hippurate prediction score from elastic net regressions. A red line indicates the linear regression line, where the dashed red lines display the 95% confidence interval. FDR—false discovery rate; b—unstandardised regression coefficients; 95% CI—95% confidence intervals.
Figure 2. Overview of microbiome alpha diversity–metabolite associations. (A) Scatter plots of Shannon diversity or species richness against urinary metabolite levels after log-transformation and regression-based normalisation. A red line indicates the linear regression line, where the dashed red lines display the 95% confidence interval. All shown associations have an FDR < 0.05. (B) Table displaying the alpha diversity metabolite associations with an FDR < 0.05. (C) Scatter plots of Shannon diversity or urinary hippurate concentrations (log-transformed and regression-based normalised) against the microbiome-based hippurate prediction score from elastic net regressions. A red line indicates the linear regression line, where the dashed red lines display the 95% confidence interval. FDR—false discovery rate; b—unstandardised regression coefficients; 95% CI—95% confidence intervals.
Metabolites 12 00308 g002
Metabolites 12 00308 g003 550
Figure 3. Microbe–host interactions regarding benzoate metabolism. (A) Microbe–host benzoate co-metabolism, as noted in AGORA2 and the Virtual Metabolic Human database (https://www.vmh.life, accessed on 28 January 2022). The reaction naming follows the AGORA2 nomenclature. (B) Benzoate production capabilities across phyla, as noted in AGORA2 as a share within phyla (left panel) and absolute number (right panel).
Figure 3. Microbe–host interactions regarding benzoate metabolism. (A) Microbe–host benzoate co-metabolism, as noted in AGORA2 and the Virtual Metabolic Human database (https://www.vmh.life, accessed on 28 January 2022). The reaction naming follows the AGORA2 nomenclature. (B) Benzoate production capabilities across phyla, as noted in AGORA2 as a share within phyla (left panel) and absolute number (right panel).
Metabolites 12 00308 g003
Metabolites 12 00308 g004 550
Figure 4. Overview of the results regarding reaction abundances after functional annotation of the Yachida et al. dataset using AGORA2. (A) Scatter plots of reaction abundances against the Shannon entropy of metagenomes after mapping onto AGORA2. A red line indicates the linear regression line, where the dashed red lines display the 95% confidence interval. All displayed associations had a false discovery rate < 0.05. (B) Table displaying the reaction abundance Shannon entropy associations for benzoate-producing, -degrading, or -transporting reactions noted in AGORA2. OR—odds ratio; 95% CI—95% confidence interval.
Figure 4. Overview of the results regarding reaction abundances after functional annotation of the Yachida et al. dataset using AGORA2. (A) Scatter plots of reaction abundances against the Shannon entropy of metagenomes after mapping onto AGORA2. A red line indicates the linear regression line, where the dashed red lines display the 95% confidence interval. All displayed associations had a false discovery rate < 0.05. (B) Table displaying the reaction abundance Shannon entropy associations for benzoate-producing, -degrading, or -transporting reactions noted in AGORA2. OR—odds ratio; 95% CI—95% confidence interval.
Metabolites 12 00308 g004
Table
Table 1. Descriptive statistics for the utilised SHIP-TREND-0 sub-cohorts.
Table 1. Descriptive statistics for the utilised SHIP-TREND-0 sub-cohorts.
SHIP-TREND with Faecal Samples (n = 3637)SHIP-TREND with Faecal Samples and NMR Metabolite Measurements (n = 951)
VariableMissing Values, %Mean (SD) or Share, %Missing Values, %Mean (SD) or Share, %
Age, years0.0051.33 (14.94)0.0050.21 (13.63)
Female, %0.0051.69%0.0056.68%
Body mass index, kg/m20.1628.02 (5.15)0.0027.37 (4.57)
Waist circumference, cm0.2790.79 (14.35)0.0088.08 (12.82)
Current smoking, %0.2526.82%0.1121.79%
Average alcohol consumption over the last 30 d, g/d0.918.83 (13.79)0.638.56 (13.31)
Diabetes a0.1611.54%0.002.73%
Hypertonia b0.3346.43%0.1139.58%
HbA1c, %0.195.34 (0.83)0.115.19 (0.56)
eGFR, mL/min0.1689.73 (18.81)0.0092.12 (17.12)
White blood cell count, Gpt/L1.956.08 (2.70)0.215.73 (1.48)
Triglycerides, mmol/L0.161.64 (1.24)0.001.42 (0.86)
Ratio of TC/HDL-C0.164.03 (1.26)0.003.93 (1.14)
Fibrinogen, g/L2.643.07 (0.74)0.953.02 (0.73)
CRP (high sensitive), mg/L4.672.52 (3.93)3.362.29 (3.67)
GGT, μkat/L0.190.70 (0.80)0.000.65 (0.63)
ALAT, μkat/L0.220.45 (0.30)0.110.44 (0.29)
ASAT, μkat/L0.300.33 (0.19)0.210.32 (0.17)
a Diabetes is defined by either an HbA1c > 6.5% or intake of antidiabetic medication. b Hypertonia is defined by the intake of antihypertensive medication or blood pressure higher than 140/90 mmHg; SHIP—Study of Health in Pomerania, SD—standard deviation, CRP—C-reactive protein, TC—total cholesterol, HDL-C—high-density lipoprotein cholesterol, HbA1c—glycated hemoglobin, GGT—gamma-glutamyl-transferase, ALAT—alanine-amino-transferase, ASAT—aspartate-amino-transferase.
Table
Table 2. Biomarker associations from multivariable regressions for urinary hippurate, Shannon entropy, and the microbiome-based hippurate prediction score.
Table 2. Biomarker associations from multivariable regressions for urinary hippurate, Shannon entropy, and the microbiome-based hippurate prediction score.
Urinary Hippurate (n = 951)Shannon Entropy (n = 3637) bMicrobiome-Based Hippurate Prediction Score (n = 3637) b
Markerb (95% CI) *p-Value *b (95% CI) *p-Value *b (95% CI) *p-Value *
Log hs-CRP−0.09 (−0.17, −0.01)2.40 × 10−2−0.01 (−0.03, 0.00)7.55 × 10−2−0.01 (−0.02, −0.00)1.65 × 10−2
Fibrinogen−0.02 (−0.09, 0.04)4.59 × 10−1−0.00 (−0.02, 0.01)6.27 × 10−10.01 (−0.00, 0.02)2.05 × 10−1
White blood cell count−0.15 (−0.28, −0.03)1.68 × 10−2−0.00 (−0.01, 0.00)1.76 × 10−1−0.00 (−0.01, 0.00)3.76 × 10−1
Triglycerides−0.02 (−0.09, 0.05)5.24 × 10−1−0.03 (−0.04, −0.02)3.83 × 10−6−0.03 (−0.04, −0.02)3.35 × 1014
Ratio of TC/HDL-C0.04 (−0.05, 0.12)3.95 × 10−1−0.01 (−0.02, 0.00)1.58 × 10−1−0.00 (−0.01, 0.01)8.08 × 10−1
Baseline glucose−0.03 (−0.09, 0.02)1.87 × 10−10.00 (−0.01, 0.01)9.32 × 10−1−0.00 (−0.01, 0.00)1.67 × 10−1
HbA1c0.004 (−0.04, 0.05)8.45 × 10−10.01 (−0.01, 0.03)2.13 × 10−10.01 (−0.00, 0.02)1.69 × 10−1
Log GGT−0.05 (−0.09, −0.01)2.05 × 10−2−0.07 (−0.10, 0.04)1.64 × 10−7−0.10 (−0.12, 0.08)4.13 × 1022
Log ALAT−0.03 (−0.06, 0.01)1.07 × 10−1−0.04 (−0.07, 0.01)9.85 × 10−3−0.07 (−0.09, 0.04)4.77 × 10−9
Log ASAT−0.03 (−0.06, 0.00)6.29 × 10−2−0.02 (−0.06, 0.01)1.62 × 10−1−0.06 (−0.08, 0.03)4.33 × 10−6
* Estimates from multivariable (mixed-effect b) linear regressions including age (nonlinear), sex, age–sex interaction terms, waist circumference (nonlinear), smoking, hypertonia, years of education, kidney function (nonlinear), urinary pH, and alcohol intake. b—unstandardised regression coefficient, 95% CI—95% confidence interval, hs-CRP—high-sensitivity C-reactive protein, TC—total cholesterol, HDL-C—high-density lipoprotein cholesterol, HbA1c—glycated hemoglobin, GGT—gamma-glutamyl transferase, ALAT—alanine aminotransferase, ASAT—aspartate aminotransferase. Bold indicates statistical significance.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hertel, J.; Fässler, D.; Heinken, A.; Weiß, F.U.; Rühlemann, M.; Bang, C.; Franke, A.; Budde, K.; Henning, A.-K.; Petersmann, A.; et al. NMR Metabolomics Reveal Urine Markers of Microbiome Diversity and Identify Benzoate Metabolism as a Mediator between High Microbial Alpha Diversity and Metabolic Health. Metabolites 2022, 12, 308. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12040308

AMA Style

Hertel J, Fässler D, Heinken A, Weiß FU, Rühlemann M, Bang C, Franke A, Budde K, Henning A-K, Petersmann A, et al. NMR Metabolomics Reveal Urine Markers of Microbiome Diversity and Identify Benzoate Metabolism as a Mediator between High Microbial Alpha Diversity and Metabolic Health. Metabolites. 2022; 12(4):308. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12040308

Chicago/Turabian Style

Hertel, Johannes, Daniel Fässler, Almut Heinken, Frank U. Weiß, Malte Rühlemann, Corinna Bang, Andre Franke, Kathrin Budde, Ann-Kristin Henning, Astrid Petersmann, and et al. 2022. "NMR Metabolomics Reveal Urine Markers of Microbiome Diversity and Identify Benzoate Metabolism as a Mediator between High Microbial Alpha Diversity and Metabolic Health" Metabolites 12, no. 4: 308. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12040308

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop