Next Article in Journal
Differential and Synergistic Effects of Low Birth Weight and Western Diet on Skeletal Muscle Vasculature, Mitochondrial Lipid Metabolism and Insulin Signaling in Male Guinea Pigs
Next Article in Special Issue
Microbiome and Human Aging: Probiotic and Prebiotic Potentials in Longevity, Skin Health and Cellular Senescence
Previous Article in Journal
Graptopetalum paraguayense Extract Ameliorates Proteotoxicity in Aging and Age-Related Diseases in Model Systems
Previous Article in Special Issue
Changes in Circulating Metabolites during Weight Loss and Weight Loss Maintenance in Relation to Cardiometabolic Risk
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 13
Issue 12
10.3390/nu13124318
nutrients-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

Examining the Interaction of the Gut Microbiome with Host Metabolism and Cardiometabolic Health in Metabolic Syndrome

by
Serena Galié
1,2,
Christopher Papandreou
2,3,*,
Pierre Arcelin
2,4,
David Garcia
5,
Antoni Palau-Galindo
1,4,
Laia Gutiérrez-Tordera
1,2,
Àlex Folch
1,2 and
Mònica Bulló
1,2,3,*
1
Department of Biochemistry and Biotechnology, Faculty of Medicine and Health Sciences, University RoviraiVirgili (URV), 43201 Reus, Spain
2
Institute of Health Pere Virgili—IISPV, University Hospital Sant Joan, 43202 Reus, Spain
3
CIBER Fisiopatología de la Obesidad y Nutrición (CIBEROBN), Instituto de Salud Carlos III, 28029 Madrid, Spain
4
Atención Basica de Salut (ABS) Reus V. Centre d’Assistència Primària Marià Fortuny, SAGESSA, 43204 Reus, Spain
5
ABS Alt Camp Oest, Centre d’Atenció Primària, 43460 Alcover, Spain
*
Authors to whom correspondence should be addressed.
Nutrients 2021, 13(12), 4318; https://doi.org/10.3390/nu13124318
Submission received: 25 October 2021 / Revised: 20 November 2021 / Accepted: 28 November 2021 / Published: 29 November 2021
(This article belongs to the Special Issue Nutrition, Metabolites, and Human Health)
Browse Figures
Versions Notes

Abstract

:
(1) Background: The microbiota-host cross-talk has been previously investigated, while its role in health is not yet clear. This study aimed to unravel the network of microbial-host interactions and correlate it with cardiometabolic risk factors. (2) Methods: A total of 47 adults with overweight/obesity and metabolic syndrome from the METADIET study were included in this cross-sectional analysis. Microbiota composition (151 genera) was assessed by 16S rRNA sequencing, fecal (m = 203) and plasma (m = 373) metabolites were profiled. An unsupervised sparse generalized canonical correlation analysis was used to construct a network of microbiota-metabolite interactions. A multi-omics score was derived for each cluster of the network and associated with cardiometabolic risk factors. (3) Results: Five multi-omics clusters were identified. Thirty-one fecal metabolites formed these clusters and were correlated with plasma sphingomyelins, lysophospholipids and medium to long-chain acylcarnitines. Seven genera from Ruminococcaceae and a member from the Desulfovibrionaceae family were correlated with fecal and plasma metabolites. Positive correlations were found between the multi-omics scores from two clusters with cholesterol and triglycerides levels. (4) Conclusions: We identified a correlated network between specific microbial genera and fecal/plasma metabolites in an adult population with metabolic syndrome, suggesting an interplay between gut microbiota and host lipid metabolism on cardiometabolic health.

1. Introduction

The metabolic variability derives from a complex and dynamic interaction between endogenous metabolism, environmental factors and the gut microbial ecosystem. Although gut microbiota is considered quite stable in adulthood, several modifiable factors can modulate its composition and activity. Whether these changes could compromise health is not yet well understood [1]. Therefore, a better understanding of how microbial communities affect or modify the complex environment of the human organism is a promising approach for the design of preventive and therapeutic strategies in different conditions [2]. In this sense, previous studies have suggested a role of some gut microbiota-derived metabolites, like bile acids and short-chain fatty acids (SCFAs), branched-chain amino acids (BCAAs), trimethylamine N-oxide, tryptophan and indole derivatives in the host-microbial cross-talk [3,4]. Furthermore, several fecal and circulating metabolites have been associated with clinical features associated with cardiometabolic risk [5], but it remains difficult to determine whether they are fully microbiota-derived or if other sources, including diet or the host itself, are also involved. An integration of metagenomics and metabolomics information may advance our knowledge on microbiota-host interactions [6]; however, to date, limited multi-omics analyses have been conducted. Two recent studies suggested a microbiota-host cross-talk analyzing the correlations of fecal and blood metabolites with gut microbiota composition by using both 16S rRNA and whole metagenomic shotgun sequencing in a large sample of UK adults (TwinsUK) [7,8]. However, whether a multi-omics profile characterizing this cross-talk could be associated with cardiometabolic health is unknown. Metabolic syndrome (MetS), a cluster of cardiometabolic conditions, is often accompanied by an imbalance of the gut microbiota [9] and alterations in metabolic pathways [5]. Identifying a host-microbial cross-talk in a population with MetS and its relationship with MetS features could further advance the understanding of biochemical processes preceding the development of cardiometabolic diseases.
Therefore, the aim of this study was to decipher the network of correlations between microbial genera, fecal and plasma metabolites using a multi-omics integrative approach in adults with overweight/obesity and MetS. Furthermore, we examined whether the identified multi-omics profiles were associated with cardiometabolic risk factors.

2. Materials and Methods

This is a cross-sectional analysis nested within the METADIET study, a randomized, controlled, crossover, dietary-intervention trial conducted in 50 adults with overweight/obesity and MetS [10]. Community-dwelling adults aged 30–65 years, with a body mass index (BMI) of 25–34.9kg/m2 who met at least three of the five diagnostic criteria of MetS and who regularly consumed a non-MedDiet were included in the study. Subjects were excluded if they suffered from type 2 diabetes (T2D), chronic diseases, had secondary obesity or related pathologies, non-controlled hypertension, LDL-cholesterol > 160 mg/dL, triglycerides > 400 mg/dL, followed specific pharmacological treatments (anti-inflammatory, corticoids, hormones or antibiotics), were alcohol or drug abusers and consumed prebiotics, probiotics or laxatives. Written informed consent was obtained from all study participants. The Institutional Review Board approved the study protocol, which accomplishes the ethical standards of the Declaration of Helsinki.
Weight, height and waist circumference were determined with calibrated scales and a wall-fixed stadiometer, and BMI was calculated. Blood pressure was measured in duplicate using a validated semiautomatic oscillometer (Omron Electronics Iberia S.A.U., Madrid, Spain). Blood and fecal samples were collected in fasting conditions before any intervention. Glucose and lipid profiles were measured using standard enzymatic automated methods. LDL-cholesterol was estimated using the Friedewald formula in subjects with triglycerides < 400 mg/dL. Circulating insulin levels were measured by commercial ELISA (Deltaclon SL, Madrid, Spain). The homeostatic model assessment of insulin resistance (HOMA-IR) was estimated [11]. Participants were instructed to collect stool samples in hermetic sterile flasks and freeze them immediately at −20 °C. Frozen samples were delivered to the laboratory within 1–2 days after collection and stored at −80 °C.
The fecal metabolomics profiling included 226 metabolites derived from a dual analytical approach (Supplementary Table S1 [12]). Ninety-four metabolites were quantified by a targeted analysis, using nuclear magnetic resonance (NMR) and liquid chromatography coupled to triple quadrupole mass spectrometry (LC-qTOF), while 132 metabolites were analyzed and semi-quantified with an untargeted approach. Particularly for the targeted analysis, NMR-based metabolic profiling of 37 metabolites included SCFAs, alcohols and organic acids, while LC-qTOF was used for the determination of 16 bile acids and 41 amino acids. The analytical procedures are specified in the Supplemental methods of [12]. Plasma metabolites (m = 378) were analyzed by different analytical platforms (Supplementary Table S2 [12]). LC-QqQ was used to analyze TMAO and derivatives, acylcarnitines, amino acids and serotonin. LC-qTOF was used for the lipidomic analysis. Measurements of total fatty acids, together with organic acids and sugar metabolites, mainly belonging to the tricarboxylic acid cycle, were obtained by a GC-qTOF analytical platform. A further description of the analytical procedures can be found in the Supplemental methods of [12].
Fecal DNA extraction was performed using a QIAmpPowerFecal DNA kit (Qiagen, Germantown, TN, USA) with a previous 5-minute lysis step (FastPrep-24-5G Homogenizer, MP Biomedicals). The 16S rRNA gene was amplified (Ion Metagenomics kitTM (Life Technology, Carlsbad, CA, USA), performing two separated PCR reactions with two primer sets (to amplify different hypervariable regions of 16S rRNA: V2, V4, V8 and V3, V6–7, V9). Amplicons were processed to obtain DNA libraries (Ion Plus Fragment Library kit and Ion Xpress Barcodes Adapters, 1-64 (Life Technology, Carlsbad, CA, USA) and adapter-ligated and nick-repaired libraries were purified (CleanNGSkit, CleanNA, Waddinxveen, The Netherlands). The libraries were further amplified (Ion Plus Fragment Library kit, Life Technology, Carlsbad, CA, USA) and quantified with Bioanalyzer (Agilent DNA 7500 Reagents, Agilent Technologies, Santa Clara, CA, USA). Equimolar amounts of all the libraries (60 µM) were sequenced in 4 different runs with Ion 520 and Ion 530 Kit-Chef (Life Technologies, Carlsbad, CA, USA) in an S5 sequencer from the Ion Torrent platform. Fastq data from sequencing were pre-processed with an adapted in-house script [13] in order to split only forward reads of each sample data into 6 subsets of 6 hypervariable regions. Forward reads from the V4 region were used for this study. Quality control, length filtering at 280 bp, and denoising of sequences with DADA2 pipeline other than taxonomy assignment were performed in QIIME2 software package using the latest version of Silva 132 as 16S rRNA gene classifier database. Finally, a further filtering step of ASV (Amplicon Sequence Variant) table at 10% prevalence cut-off at the taxonomic level of genera was achieved in R (using phyloseq package functions).
Twenty-three fecal metabolites were removed from the analysis due to the high number of missing values (>20%), and a total of 203 metabolites were included. In those metabolites with less than 20%, missing values were imputed using the random forest imputation approach (“missForest” function of “randomForest” R package version 4.6-14). The concentrations of metabolites were normalized and scaled to multiples of 1 SD with the rank-based inverse normal transformation.
Five out of 378 plasma metabolites were removed because of the high number of missing values (>20%), and the remaining missing values were imputed using the same approach as above. The rank-based inverse normal transformation was used to normalize their concentrations.
Baseline values of absolute abundances of the ASV table at the taxonomic level of 151 genera were center-normalized with clr function from the “composition” package on R (version 1.4-40). Due to the high dimensionality and collinear nature of the data, a Sparse Generalized Canonical Correlation Analysis (SGCCA) in an unsupervised mode was conducted to select the most relevant variables from three omics datasets (fecal metabolites, plasma metabolites, 16S rRNA). The model was implemented by using the canonical mode of the “wrapper.sgcca” algorithm [14], available in the mixOmics package in R (version 6.14.0) (Available online: http://mixomics.org/ (accessed on 20 November 2021)). The tuning procedure to identify the optimal number of components for each omic dataset was performed by a separated performance analysis on both tri and dual-omics datasets by using the “perform” function, which is available for both “block.splsda” and spls models. We then validated the choice of the optimal parameters by evaluating the AVE values of our obtained model of spls models (perf function in mixOmics package in R, version 6.14.0). A network analysis was implemented by the network function with a cut-off correlation value of 0.6 [15]. The network was further analyzed in Cytoscape software (version 3.8.2) (Available online: https://cytoscape.org/ (accessed on 20 November 2021)) in order to better visualize the presence of connected components and the relevant associations between fecal and plasma metabolites and microbial genera. A prefuse force-directed layout with a color grouping visualization based on the nature of components was selected to better feature the clusters of the network. A prefuse force-directed layout in Cytoscape is based on the force-directed-layout algorithm, which uses repulsive forces between nodes and attractive forces between adjacent nodes. The multi-omics score was calculated based on the weighted sum of the selected components in each cluster. The weights were obtained from the scaling with the eigengene centrality scores of each component in each module measured with the Page Ranking algorithm in Cytoscape. Linear regression models were fitted to examine the association between the derived scores and cardiometabolic risk factors (glucose, insulin, HOMA-IR, total cholesterol, HDLc, LDLc, VLDLc, triglyceride levels, systolic and diastolic blood pressure) adjusting for age, sex and BMI values. Furthermore, a Pearson partial correlation analysis adjusting for age, sex and BMI, was implemented between individual components of each multi-omic score and the cardiometabolic parameters significantly correlated with the scores. Values of the selected metabolites in feces and plasma, as well as the selected genera, were correlated with the cardiometabolic parameters, using the “associate” function from the “microbiome” package in R (version 1.12.0). A Benjamini-Hochberg false discovery rate (FDR) approach [16] was used to correct p-values for multiple testing.
All analyses were performed using R, version 3.6.2. All tests were two-sided, and significance was defined as p < 0.05.

3. Results

Of the 50 participants initially included in the METADIET study, three were excluded because of the unavailability of either 16S rRNA sequencing data or fecal metabolomics, resulting in a final total number of 47 participants. The general characteristics of participants are shown in Table 1. The mean age of participants was 50.6 ± 7.13 years, and the mean BMI was 30.5 ± 2.28 Kg/m2.

3.1. Multi-Omics Network of Correlations between Gut Microbiota, Fecal and Plasma Metabolites

Figure 1 shows the network deriving from the correlation network analysis with a cut-off point value of 0.6. Ninety-four nodes and 263 edges characterized the resulting network, which displays the most relevant correlations deriving from the SGCCA between the different multi-omics data. The network shows five distinct multi-omics clusters, including a total number of 9 genera, 31 fecal metabolites and 41 plasma metabolites. A detailed description of the network and cluster characteristics is given in Supplementary Table S3 and Figure S1.
In the tri-omics cluster one, we observed a hub of connected components, mainly constituted by a fecal primary bile acid [chenodeoxycholic acid (CDCA)] and its derivative secondary bile acid, ursodeoxycholic acid (UDCA), fecal propionate, fecal derivatives of arachidonic acid, as well as, fecal LPC 16:0 and LPC 20:4, that were negatively correlated with three uncultured genera from the Ruminococcaceae family and Christensenellaceae R7 group. These microbial genera were also positively correlated with several plasma sphingomyelins (SMs) with one and two double bonds. Cluster two mainly consisted of fecal amino acids positively correlated with plasma species of LPCs (16:0, 17:1, 18:2, 18:1, 18:3, 20:5, 22:6) and LPIs (18:1, 18:2, 20:4). Plasma levels of LPI 22:6 were also negatively correlated with Subdoligranum. The multi-omics cluster three was constituted by plasma levels of medium/long-chain mainly saturated acylcarnitines (C8:0, C10:0, C12:0, C12:1, C14:1, C14:2) found to be positively correlated with fecal levels of spermidine and cadaverine. Plasma levels of C10:0, C12:0, C14:1, C12:0-OH were also negatively correlated with two genera from the Ruminococcaceae family (Ruminococcaceae UCG09 group and Acetanaereobacterium). The multi-omics cluster four was formed by fecal fatty acids such as myristic acid, methyladipic acid, capric acid and dodecanoic acid that were positively correlated with plasma levels of different androsterones, including dehydroepiandrosterone sulfate (DHEAS), and with a member of Firmicutes (Anaerotruncus). The multi-omics cluster five was a dual-omics hub of plasma levels of phosphocholines (PCs 32:2, 32:1, 34:4, 30:0), triglycerides (TGs 46:1, 46:2, 48:1, 48:2, 48:3), lysophosphatidylcholine (LPC) 14:0, lysophosphatidylethanolamine (LPE) 14:0 and lysophosphatidylinositol (LPI) 14:0 negatively correlated with a genus from the Desulfovibrionaceae family. Better visualization of the direction of each correlation is shown in Supplementary Figures S2–S4 [12].

3.2. Associations between Multi-Omics Scores and Cardiometabolic Risk Factors

The linear regression analysis of associations between the five multi-omics scores (one for each cluster) and cardiometabolic parameters are shown in Table 2.
Significant positive associations of the multi-omics cluster one score with total and LDL cholesterol levels were observed. Significant positive associations were also found between the multi-omics cluster five score and levels of cholesterol (total, LDL and VLDL) and triglycerides.
Figure 2 and Figure 3 show the heatmaps of Pearson correlations between individual components of the multi-omics clusters one and five, and cardiometabolic parameters.
Figure 2 shows significant positive correlations between specific plasma levels of certain SMs (32:1, 34:1, 38:1, 39:1, 40.1, 41:1, 43:1) from the multi-omics cluster one and cholesterol (total and LDL) levels as well as a positive correlation of SM32:1 and SM34:1 with HDLc. Figure 3 shows significant positive correlations of PCs (30:0, 32:1, 32:2, 34:4), LPC 14:0, LPI 14:0, LPE 14:0 and different TGs species (46:2, 48:1, 48:2, 48:3) from cluster five with serum levels of cholesterol (total, LDL, VLDL) and triglycerides.

4. Discussion

In this secondary analysis conducted in the METADIET study, we identified a network of correlations between microbial genera and specific fecal and circulating metabolites that constituted five different multi-omics clusters. Two of the identified clusters, and especially plasma lipid species, were associated with cardiometabolic risk factors highlighting their importance on cardiometabolic health [17]. Despite the absence of a direct correlation between the identified microbial genera and cardiometabolic health features, their co-occurrence with their microbial functional readout represented by fecal metabolites and known secondary biomarkers of cardiometabolic health in plasma offers new insights for further mechanistic studies.
Cluster one demonstrated a cross-talk between circulating SM species and fecal bile acids and arachidonic acid derivatives in which members from the Ruminococcaceae and Christensenellaceae families were involved. The correlation of these interrelated metabolic components with total and LDL cholesterol adds evidence on the role of SMs on cholesterol metabolism. Indeed, SMs are the most abundant sphingolipids in lipoproteins, such as LDL, and their levels in plasma have been previously associated with cardiovascular disease risk [18]. SMs also play a role in cholesterol homeostasis [19] by affecting LDL’s ability to bind surface receptors and subsequent internalization [19]. A growing body of evidence regarding their favorable interaction with sterols indicates SMs as a potential key regulator of cholesterol distribution within cellular membranes. Interestingly, the Christensenellaceae R7 group and Ruminococcaceae UCG002, which we found to be positively correlated with these SM species, have both been negatively associated with VLDL diameter in an observational study conducted in The Netherlands [20]. Similarly, other studies have also found negative associations of the Christensenellaceae family with total and LDL cholesterol and apolipoprotein particles [21,22], suggesting its role in overall cholesterol homeostasis. Whether these associations could be partially explained by SMs is a hypothesis that needs further investigation. At the same time, we observed negative correlations between the Christensenellaceae R7 group and Ruminococcaceae UCG002 with fecal bile acids CDCA and its microbial derivative UDCA [23], supporting the mediation of the microbial activity on fecal bile acids in cluster two. Reinforcing our findings of the potential involvement of this cluster in cholesterol metabolism, a recent study conducted in an Italian elderly cohort was able to characterize a healthier metabolic profile by Christensenellaceae-enterotype, which was also associated with an improved visceral lipid composition, as well as with a trend towards lower serum levels of CDCA [24], whose harmful impact on cholesterol metabolism has been previously demonstrated [25]. On the contrary, the therapeutical benefits of secondary bile acid UDCA in reducing cholesterol solubilization in blood has been demonstrated [23]. In this cluster, circulating levels of SMs were also positively correlated with fecal 13-methylmyristic acid, also known as 13-methytetradecanoic acid (13-MTD). 13-MTD is a saturated iso-fatty acid that derives from microbial fermentation and has been extensively studied for its apoptotic properties, especially in certain cancer cells [26], other than being proposed as a marker of adipose tissue turnover [27]. In this same cluster, fecal levels of the SCFA propionate were also related to plasma levels of different SM species (38:2, 41:2) and Ruminococcaceae UCG-002. Despite the fact that no specific genera from Ruminococcaceae have been previously addressed to propionate production pathways [28], we could speculate its indirect involvement in propionate metabolism.
In cluster two, fecal amino acids, including microbial-derived BCAAs, aromatic amino acids and proline, were correlated with plasma LPC and LPI species, and LPI 22:6 was negatively correlated with Subdoligranulum. The genus Subdoligranulum, from the Ruminococcaceae family, is a Gram-negative, strictly anaerobe and butyrate-producer [29], which has been recently proposed as beneficial bacteria for metabolic health [30]. Strong evidence supports a relationship between circulating BCAAs and cardiometabolic diseases [31,32], and a similar pattern seems to be reflected in feces regarding the association of fecal BCAAs with insulin resistance [33]. However, in our study, fecal BCAAs were not correlated with cardiometabolic risk parameters. The positive correlation we found between fecal levels of BCAAs and LPCs and LPI species in plasma is not supported by previous findings; therefore, the potential intermediation of this metabolic interaction in cardiometabolic health needs further mechanistic studies.
In cluster three, plasma levels of long-chain acylcarnitines were negatively correlated with fecal levels of the polyamines cadaverine and spermidine, as well as negatively correlated with members of the Ruminococcaceae family (Acetanaerobacterium and Ruminococcaceae UCG-009). Plasma levels of long (C12, C14) and medium (C8, C10) chain acylcarnitines are associated with β oxidation of fatty acids [34]. A similar pattern of correlations between medium and long-chain acylcarnitines and both spermidine and cadaverine were associated with increased BMI in the Northern Finland Birth Cohort, while some acylcarnitines were correlated with an unknown genus from the Firmicutes phylum [35]. Acetanaerobacterium and Ruminococcaceae UCG-009 are both butyrate-producers belonging to the Ruminococcaceae family and are known for their potential protective role in the metabolism of high-fat diet-induced obese mice [36]. Different members of the Ruminococcaceae family have been previously correlated with lower levels of medium-chain acylcarnitines in a randomized crossover clinical trial comparing vegetarians versus meat-eaters [37]. Other studies have also demonstrated the ability of some bacteria to ferment cadaverine for butyrate production [38], partially explaining our observations about the co-appearance of the Ruminococcaceae members with fecal cadaverine. Therefore, there might be a potential involvement of this multi-omics profile of butyrate-producers in energy homeostasis through the regulation of fatty acids oxidation.
In cluster four, Anaerotruncus from the family of Clostridiaceae was correlated with fecal levels of dodecanoic acid and methyladipic acid, as well as with plasma levels of different steroid species of androsterone sulfate, including DHEAS. Previous studies have observed the involvement of Anaerotruncus in the biosynthesis of both steroids and terpenoids in colorectal cancer [39]. Although a recent in vivo study found a relationship between lower plasma concentrations of DHEAS and impaired glucose tolerance [40], we did not find any significant correlation between plasma levels of steroids and cardiometabolic parameters.
Finally, a dual-omics profile was identified in our study. To our knowledge, this is the first study finding negative correlations of plasma phospholipid species and TGs with 46 or 48 carbon atoms and <3 double bonds with bacteria from the Desulfovibrioaceae family. Previous evidence suggests that the Desulfovibrionaceae family members might have an adverse effect on dyslipidemia, even if the biological mechanisms behind this association remain still unknown [41]. Therefore, our results in relation to the correlations of lipid species with serum cholesterol (total, LDL, VLDL) and triglycerides support the involvement of this microbial-host cross-talk in lipid metabolism and transport.
The present study has some strengths that deserve to be mentioned, like being one of the few studies analyzing the microbial-host cross-talk by using the combination of three different omics, the use of a multi-platform metabolomics analysis, which allowed to combine targeted and untargeted analytical approaches, in order to cover a wide range of metabolites. Furthermore, the availability of a pretty uniform population characterized by MetS allowed us to relate the microbial-host cross-talk with cardiometabolic risk factors. Regarding limitations, this is a cross-sectional study in a small-sized population, and no causal relationship between the microbiome and the metabolome or between the metabolome and the panel of cardiometabolic risk factors can be inferred. Due to the fact that the implemented model is right below the limits to be considered as a robust model, some of our findings could be the result of spurious correlations. Furthermore, the use of 16S rRNA sequencing does not allow unraveling of the taxonomical composition of gut microbiota up to the strain level. A further limitation of the study is the lack of a “healthy” control group. This may have partially accounted for the lack of significant associations between the multi-omics scores and cardiometabolic measures such as insulin and HOMA-IR due to their low variation among the study population. As well, we evaluated a sample of individuals with MetS that could limit the generalizability of our results to other populations. Other larger studies would be necessary to confirm our novel findings, ideally, prospectively.

5. Conclusions

In conclusion, this study demonstrated a network of five host-microbial cross-talk profiles in adults with MetS. A total of 41 plasma and 31 fecal metabolites, as well as 9 microbial genera, were identified in 5 clusters. Plasma levels of different lipid species (SMs, LPCs, LPIs and long-chain acylcarnitines) were mainly correlated with members from the Ruminococcaceae and Christenseneellaceae families, as well as with fecal bile acids, BCAAs, AAs and polyamines. The plasma lipid correlations with serum total LDL and VLDL cholesterol and triglycerides levels suggest a link between the microbial-host cross-talk and cardiometabolic health.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/nu13124318/s1, Figure S1. Description of the connected components of the network; Figure S2. Heatmap of correlations from the network similarity matrix between selected plasma and fecal metabolites; Figure S3. Heatmap of correlations from the network similarity matrix between selected plasma metabolites and microbial genera; Figure S4. Heatmap of correlations from the network similarity matrix between selected fecal metabolites and microbial genera; Supplementary Table S1. List of fecal metabolites analyzed; Supplementary Table S2. List of plasma metabolites analyzed; Supplementary Table S3. Summary statistics of the network analysis.

Author Contributions

M.B. designed the research and was the coordinator of the study; M.B. and S.G. recruited participants and conducted the research; S.G. conducted the bioinformatics analysis; S.G. and C.P. conducted the statistical analyses; S.G., C.P. and M.B., drafted the paper; P.A., D.G., A.P.-G., L.G.-T. and À.F. and all other authors revised the manuscript for important intellectual content, and read and approved the final version. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Fondo de Investigación Sanitaria (PI16/00516) and cofounded by the FondoEuropeo de Desarrollo Regional (FEDER). SG received a pre-doctoral fellowship from the Generalitat de Catalunya’s Department of Universities (2018FI_B_00444). CP is a recipient of the Instituto de Salud Carlos III Miguel Servet fellowship (grant CP 19/00189).

Institutional Review Board Statement

The trial was registered in the ISRCTN with identifier ISRCTN88780852, https://doi.org/10.1186/ISRCTN88780852 (accessed on 20 November 2021).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Further data will be provided under request to the authors.

Acknowledgments

Metabolomics analyses and 16S rRNA gene sequencing were conducted in the facilities of the Centre for Omics Sciences (COS) Joint Unit of Universitat Rovira i Virgili-Eurecat.

Conflicts of Interest

The authors declare no conflict of interest in the study.

References

  1. Xie, H.; Guo, R.; Zhong, H.; Feng, Q.; Lan, Z.; Qin, B.; Ward, K.J.; Jackson, M.; Xia, Y.; Chen, X.; et al. Shotgun Metagenomics of 250 Adult Twins Reveals Genetic and Environmental Impacts on the Gut Microbiome. Cell Syst. 2016, 3, 572–584.e3. [Google Scholar] [CrossRef] [Green Version]
  2. Turnbaugh, P.; Gordon, J.I. An Invitation to the Marriage of Metagenomics and Metabolomics. Cell 2008, 134, 708–713. [Google Scholar] [CrossRef] [Green Version]
  3. Rooks, M.G.; Garrett, W.S. Gut microbiota, metabolites and host immunity. Nat. Rev. Immunol. 2016, 16, 341–352. [Google Scholar] [CrossRef]
  4. Sonnenburg, J.L.; Bäckhed, F. Diet–microbiota interactions as moderators of human metabolism. Nat. Cell Biol. 2016, 535, 56–64. [Google Scholar] [CrossRef] [PubMed]
  5. Monnerie, S.; Comte, B.; Ziegler, D.; Morais, J.A.; Pujos-Guillot, E.; Gaudreau, P. Metabolomic and Lipidomic Signatures of Metabolic Syndrome and its Physiological Components in Adults: A Systematic Review. Sci. Rep. 2020, 10, 669. [Google Scholar] [CrossRef] [Green Version]
  6. Chen, M.X.; Wang, S.-Y.; Kuo, C.-H.; Tsai, I.-L. Metabolome analysis for investigating host-gut microbiota interactions. J. Formos. Med. Assoc. 2019, 118, S10–S22. [Google Scholar] [CrossRef]
  7. Zierer, J.; Jackson, M.A.; Kastenmüller, G.; Mangino, M.; Long, T.; Telenti, A.; Mohney, R.P.; Small, K.S.; Bell, J.; Steves, C.; et al. The fecal metabolome as a functional readout of the gut microbiome. Nat. Genet. 2018, 50, 790–795. [Google Scholar] [CrossRef] [PubMed]
  8. 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] [PubMed] [Green Version]
  9. Festi, D.; Schiumerini, R.; Eusebi, L.H.; Marasco, G.; Taddia, M.; Colecchia, A. Gut microbiota and metabolic syndrome. World J. Gastroenterol. 2014, 20, 16079–16094. [Google Scholar] [CrossRef] [PubMed]
  10. Galié, S.; García-Gavilán, J.; Papandreou, C.; Camacho-Barcía, L.; Arcelin, P.; Palau-Galindo, A.; Rabassa, A.; Bulló, M. Effects of Mediterranean Diet on plasma metabolites and their relationship with insulin resistance and gut microbiota composition in a crossover randomized clinical trial. Clin. Nutr. 2021, 40, 3798–3806. [Google Scholar] [CrossRef]
  11. Antuna-Puente, B.; Disse, E.; Rabasa-Lhoret, R.; Laville, M.; Capeau, J.; Bastard, J.-P. How can we measure insulin sensitivity/resistance? Diabetes Metab. 2011, 37, 179–188. [Google Scholar] [CrossRef] [PubMed]
  12. Galiè, S. Supplementary Data, Mendeley Data, V3, 2021. Available online: https://data.mendeley.com/datasets/nhh5cpsz4d/3 (accessed on 20 November 2021). [CrossRef]
  13. Mas-Lloret, J.; Obón-Santacana, M.; Ibáñez-Sanz, G.; Guinó, E.; Pato, M.L.; Rodriguez-Moranta, F.; Mata, A.; García-Rodríguez, A.; Moreno, V.; Pimenoff, V.N. Gut microbiome diversity detected by high-coverage 16S and shotgun sequencing of paired stool and colon sample. Sci. Data 2020, 7, 90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Tenenhaus, A.; Philippe, C.; Guillemot, V.; Le Cao, K.-A.; Grill, J.; Frouin, V. Variable selection for generalized canonical correlation analysis. Biostatistics 2014, 15, 569–583. [Google Scholar] [CrossRef] [Green Version]
  15. Benedetti, E.; Pučić-Baković, M.; Keser, T.; Gerstner, N.; Büyüközkan, M.; Štambuk, T.; Selman, M.H.J.; Rudan, I.; Polašek, O.; Hayward, C.; et al. A strategy to incorporate prior knowledge into correlation network cutoff selection. Nat. Commun. 2020, 11, 5153. [Google Scholar] [CrossRef] [PubMed]
  16. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  17. Amar, J. Microbiota–Host Crosstalk: A Bridge between Cardiovascular Risk Factors, Diet, and Cardiovascular Disease. Am. J. Hypertens. 2018, 31, 941–944. [Google Scholar] [CrossRef]
  18. Jiang, X.-C.; Paultre, F.; Pearson, T.A.; Reed, R.G.; Francis, C.K.; Lin, M.; Berglund, L.; Tall, A.R. Plasma Sphingomyelin Level as a Risk Factor for Coronary Artery Disease. Arter. Thromb. Vasc. Biol. 2000, 20, 2614–2618. [Google Scholar] [CrossRef] [Green Version]
  19. Zalloua, P.; Kadar, H.; Hariri, E.; Farraj, L.A.; Brial, F.; Hedjazi, L.; Le Lay, A.; Colleu, A.; Dubus, J.; Touboul, D.; et al. Untargeted Mass Spectrometry Lipidomics identifies correlation between serum sphingomyelins and plasma cholesterol. Lipids Health Dis. 2019, 18, 38. [Google Scholar] [CrossRef]
  20. Vojinovic, D.; Radjabzadeh, D.; Kurilshikov, A.; Amin, N.; Wijmenga, C.; Franke, L.; Ikram, M.A.; Uitterlinden, A.G.; Zhernakova, A.; Fu, J.; et al. Relationship between Gut Microbiota and Circulating Metabolites in Population-Based Cohorts. Nat. Commun. 2019, 10, 5813. [Google Scholar] [CrossRef] [Green Version]
  21. López-Contreras, B.E.; Morán-Ramos, S.; Villarruel-Vázquez, R.; Macías-Kauffer, L.; Villamil-Ramírez, H.; León-Mimila, P.; Vega-Badillo, J.; Sánchez-Muñoz, F.; Llanos-Moreno, L.E.; Canizalez-Román, A.; et al. Composition of gut microbiota in obese and normal-weight Mexican school-age children and its association with metabolic traits. Pediatr. Obes. 2018, 13, 381–388. [Google Scholar] [CrossRef]
  22. Russell, D.W. The Enzymes, Regulation, and Genetics of Bile Acid Synthesis. Annu. Rev. Biochem. 2003, 72, 137–174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Tavella, T.; Rampelli, S.; Guidarelli, G.; Bazzocchi, A.; Gasperini, C.; Pujos-Guillot, E.; Comte, B.; Barone, M.; Biagi, E.; Candela, M.; et al. Elevated gut microbiome abundance of Christensenellaceae, Porphyromonadaceae and Rikenellaceae is associated with reduced visceral adipose tissue and healthier metabolic profile in Italian elderly. Gut Microbes 2021, 13, 1–19. [Google Scholar] [CrossRef]
  24. Porez, G.; Prawitt, J.; Gross, B.; Staels, B. Bile acid receptors as targets for the treatment of dyslipidemia and cardiovascular disease. J. Lipid Res. 2012, 53, 1723–1737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Simental-Mendía, L.E.; Simental-Mendía, M.; Sánchez-García, A.; Banach, M.; Serban, M.-C.; Cicero, A.F.G.; Sahebkar, A. Impact of ursodeoxycholic acid on circulating lipid concentrations: A systematic review and meta-analysis of randomized placebo-controlled trials. Lipids Health Dis. 2019, 18, 1–13. [Google Scholar] [CrossRef] [Green Version]
  26. Klein, R.A.; Halliday, D.; Pittet, P.G. The use of 13-methyltetradecanoic acid as an indicator of adipose tissue turnover. Lipids 1980, 15, 572–579. [Google Scholar] [CrossRef]
  27. Louis, P.; Flint, H.J. Formation of propionate and butyrate by the human colonic microbiota. Environ. Microbiol. 2017, 19, 29–41. [Google Scholar] [CrossRef] [Green Version]
  28. Waters, J.L.; Ley, R.E. The human gut bacteria Christensenellaceae are widespread, heritable, and associated with health. BMC Biol. 2019, 17, 83. [Google Scholar] [CrossRef]
  29. Van Hul, M.; Le Roy, T.; Prifti, E.; Dao, M.C.; Paquot, A.; Zucker, J.-D.; Delzenne, N.M.; Muccioli, G.G.; Clément, K.; Cani, P.D. From correlation to causality: The case of Subdoligranulum. Gut Microbes 2020, 12, 1–13. [Google Scholar] [CrossRef]
  30. Ruiz-Canela, M.; Toledo, E.; Clish, C.; Hruby, A.; Liang, L.; Salas-Salvadó, J.; Razquin, C.; Corella, D.; Estruch, R.; Ros, E.; et al. Plasma Branched-Chain Amino Acids and Incident Cardiovascular Disease in the PREDIMED Trial. Clin. Chem. 2016, 62, 582–592. [Google Scholar] [CrossRef] [Green Version]
  31. Ruiz-Canela, M.; Guasch-Ferré, M.; Toledo, E.; Clish, C.B.; Razquin, C.; Liang, L.; Wang, D.D.; Corella, D.; Estruch, R.; Hernáez, Á.; et al. Plasma branched chain/aromatic amino acids, enriched Mediterranean diet and risk of type 2 diabetes: Case-cohort study within the PREDIMED Trial. Diabetologia 2018, 61, 1560–1571. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Palomino-Schätzlein, M.; Mayneris-Perxachs, J.; Caballano-Infantes, E.; Rodríguez, M.A.; Palomo-Buitrago, M.-E.; Xiao, X.; Mares, R.; Ricart, W.; Simó, R.; Herance, J.R.; et al. Combining metabolic profiling of plasma and faeces as a fingerprint of insulin resistance in obesity. Clin. Nutr. 2020, 39, 2292–2300. [Google Scholar] [CrossRef]
  33. Lu, Y.; Wang, Y.; Ong, C.-N.; Subramaniam, T.; Choi, H.W.; Yuan, J.-M.; Koh, W.-P.; Pan, A. Metabolic signatures and risk of type 2 diabetes in a Chinese population: An untargeted metabolomics study using both LC-MS and GC-MS. Diabetologia 2016, 59, 2349–2359. [Google Scholar] [CrossRef] [Green Version]
  34. Gannon, N.P.; Schnuck, J.K.; Vaughan, R.A. BCAA Metabolism and Insulin Sensitivity-Dysregulated by Metabolic Status? Mol. Nutr. Food Res. 2018, 62, e1700756. [Google Scholar] [CrossRef]
  35. Loftfield, E.; Herzig, K.-H.; Caporaso, J.G.; Derkach, A.; Wan, Y.; Byrd, D.A.; Vogtmann, E.; Männikkö, M.; Karhunen, V.; Knight, R.; et al. Association of Body Mass Index with Fecal Microbial Diversity and Metabolites in the Northern Finland Birth Cohort. Cancer Epidemiol. Biomark. Prev. 2020, 29, 2289–2299. [Google Scholar] [CrossRef] [PubMed]
  36. Hildebrandt, M.A.; Hoffmann, C.; Sherrill–Mix, S.A.; Keilbaugh, S.A.; Hamady, M.; Chen, Y.-Y.; Knight, R.; Ahima, R.S.; Bushman, F.; Wu, G.D. High-Fat Diet Determines the Composition of the Murine Gut Microbiome Independently of Obesity. Gastroenterology 2009, 137, 1716–1724.e2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Djekic, D.; Shi, L.; Brolin, H.; Carlsson, F.; Särnqvist, C.; Savolainen, O.; Cao, Y.; Bäckhed, F.; Tremaroli, V.; Landberg, R.; et al. Effects of a Vegetarian Diet on Cardiometabolic Risk Factors, Gut Microbiota, and Plasma Metabolome in Subjects with Ischemic Heart Disease: A Randomized, Crossover Study. J. Am. Hear. Assoc. 2020, 9, e016518. [Google Scholar] [CrossRef]
  38. Roeder, J.; Schink, B. Syntrophic Degradation of Cadaverine by a Defined Methanogenic Coculture. Appl. Environ. Microbiol. 2009, 75, 4821–4828. [Google Scholar] [CrossRef] [Green Version]
  39. Loke, M.F.; Chua, E.G.; Gan, H.M.; Thulasi, K.; Wanyiri, J.W.; Thevambiga, I.; Goh, K.L.; Wong, W.F.; Vadivelu, J. Metabolomics and 16S rRNA sequencing of human colorectal cancers and adjacent mucosa. PLoS ONE 2018, 13, e0208584. [Google Scholar] [CrossRef]
  40. Ma, J.; Yue, J.; Huang, R.; Liao, Y.; Li, S.; Liu, W. Reversion of aging-related DHEAS decline in mouse plasma alleviates aging-related glucose tolerance impairment by potentiation of glucose-stimulated insulin secretion of acute phase. Biochem. Biophys. Res. Commun. 2018, 500, 671–675. [Google Scholar] [CrossRef] [PubMed]
  41. Qin, J.; Li, Y.; Cai, Z.; Li, S.; Zhu, J.; Zhang, F.; Liang, S.; Zhang, W.; Guan, Y.; Shen, D.; et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 2012, 490, 55–60. [Google Scholar] [CrossRef]
Nutrients 13 04318 g001 550
Figure 1. Network analysis of correlation for the multi-omics clusters.
Figure 1. Network analysis of correlation for the multi-omics clusters.
Nutrients 13 04318 g001
Nutrients 13 04318 g002 550
Figure 2. Heatmap of Pearson partial correlations between individual components of the multi-omics cluster 1 and cardiometabolic risk factors.
Figure 2. Heatmap of Pearson partial correlations between individual components of the multi-omics cluster 1 and cardiometabolic risk factors.
Nutrients 13 04318 g002
Nutrients 13 04318 g003 550
Figure 3. Heatmap of Pearson partial correlations between individual components of the multi-omics cluster 5 and selected cardiometabolic risk factor.
Figure 3. Heatmap of Pearson partial correlations between individual components of the multi-omics cluster 5 and selected cardiometabolic risk factor.
Nutrients 13 04318 g003
Table
Table 1. General characteristics of study participants.
Table 1. General characteristics of study participants.
CharacteristicsBaseline
n = 47
Age (years)50.6 (48.6, 52.6)
Women n (%)30 (63.8%)
BMI (Kg/m2)30.5 (29.9, 31.2)
Waist Circumference (cm)102.1 (99.3, 104.8)
SBP (mmHg)135.1 (131.7, 138.6)
DBP (mmHg)85.0 (82.3, 87.7)
Total Cholesterol (mg/dL)215.3 (206.1, 224.5)
LDLc (mg/dL)135.8 (128.1, 143.6)
HDLc (mg/dL)50.5 (47.6, 53.5)
VLDLc (mg/dL)28.1 (24.7, 31.5)
Triglycerides (mg/dL)147.7 (126.3, 169.2)
Glucose (mg/dL)100.2 (96.5, 103.9)
Insulin (mcUI/mL)13.4 (11.3, 15.5)
HOMA-IR3.3 (2.8, 3.9)
All values are given as means (95% CI). Abbreviations: BMI; body mass index, SBP; systolic blood pressure, DBP; diastolic blood pressure, LDLc; low-density lipoprotein cholesterol, HDLc; high-density lipoprotein cholesterol, VLDLc; very-low-density lipoprotein cholesterol, HOMA-IR; homeostatic model assessment of insulin resistance.
Table
Table 2. Linear regression analysis examining the associations of 1-SD increase of multi-omics scores with cardiometabolic risk factors.
Table 2. Linear regression analysis examining the associations of 1-SD increase of multi-omics scores with cardiometabolic risk factors.
Multi-Omic Score 1Multi-Omic Score 2Multi-Omic Score 3Multi-Omic Score 4Multi-Omic Score 5
FactorMean ± SEp 1Mean ± SEp 1Mean ± SEp 1Mean ± SEp 1Mean ± SEp 1
Cholesterol (mg/dL)1.666 ± 0.4350.0010.249 ± 0.0870.0700.099 ± 0.2150.9830.121 ± 0.3110.7430.621 ± 0.1480.0003
LDLc (mg/dL)1.269 ± 0.3360.0050.161 ± 0.06920.0830.062 ± 0.1660.9830.079 ± 0.2390.7430.327 ± 0.1260.032
HDLc (mg/dL)0.374 ± 0.1490.0530.0452 ± 0.0290.330−0.001 ± 0.0680.9830.089 ± 0.0980.7430.0326 ± 0.0550.803
VLDLc (mg/dL)−0.068 ± 0.1790.7550.036 ± 0.0330.4950.030 ± 0.0760.9830.050 ± 0.1100.7430.224 ± 0.0520.0003
Triglycerides (mg/dL)0.373 ± 1.1860.7550.234 ± 0.2210.4950.222 ± 0.5050.983−0.498 ± 0.7280.7431.420 ± 0.3510.0003
Glucose (mg/dL)0.358 ± 0.1910.1700.087 ± 0.0350.0830.132 ± 0.08240.9830.116 ± 0.1210.7430.043 ± 0.0690.803
Insulin (mcUI/mL)0.088 ± 0.1140.7380.007 ± 0.0210.838−0.028 ± 0.04860.983−0.062 ± 0.0700.743−0.010 ± 0.0400.803
HOMA IR0.034 ± 0.0290.4720.004 ± 0.0050.656−0.0005 ± 0.0120.983−0.008 ± 0.0180.743−0.003 ± 0.0100.803
SBP (mmHg)0.077 ± 0.1910.7550.019 ± 0.0360.7490.105 ± 0.0800.983−0.117 ± 0.1170.7430.019 ± 0.0670.803
DBP (mmHg)0.065 ± 0.1480.7550.002 ± 0.0280.9510.006 ± 0.0630.9830.045 ± 0.0910.7430.028 ± 0.0520.803
Values are given as means ± standard error. Model was adjusted for sex, age and BMI. Abbreviations, DBP; diastolic blood pressure, LDLc; low-density lipoprotein cholesterol, HDLc; high-density lipoprotein cholesterol, VLDLc; very-low-density lipoprotein cholesterol, HOMA-IR; homeostatic model assessment of insulin resistance, SBP; systolic blood pressure, DBP; diastolic blood pressure. 1 Adjusted with the Benjamini-Hochberg False Discovery Rate method.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Galié, S.; Papandreou, C.; Arcelin, P.; Garcia, D.; Palau-Galindo, A.; Gutiérrez-Tordera, L.; Folch, À.; Bulló, M. Examining the Interaction of the Gut Microbiome with Host Metabolism and Cardiometabolic Health in Metabolic Syndrome. Nutrients 2021, 13, 4318. https://doi.org/10.3390/nu13124318

AMA Style

Galié S, Papandreou C, Arcelin P, Garcia D, Palau-Galindo A, Gutiérrez-Tordera L, Folch À, Bulló M. Examining the Interaction of the Gut Microbiome with Host Metabolism and Cardiometabolic Health in Metabolic Syndrome. Nutrients. 2021; 13(12):4318. https://doi.org/10.3390/nu13124318

Chicago/Turabian Style

Galié, Serena, Christopher Papandreou, Pierre Arcelin, David Garcia, Antoni Palau-Galindo, Laia Gutiérrez-Tordera, Àlex Folch, and Mònica Bulló. 2021. "Examining the Interaction of the Gut Microbiome with Host Metabolism and Cardiometabolic Health in Metabolic Syndrome" Nutrients 13, no. 12: 4318. https://doi.org/10.3390/nu13124318

APA Style

Galié, S., Papandreou, C., Arcelin, P., Garcia, D., Palau-Galindo, A., Gutiérrez-Tordera, L., Folch, À., & Bulló, M. (2021). Examining the Interaction of the Gut Microbiome with Host Metabolism and Cardiometabolic Health in Metabolic Syndrome. Nutrients, 13(12), 4318. https://doi.org/10.3390/nu13124318

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