Serum metabolite profiling yields insights into health promoting effect of A. muciniphila in human volunteers with a metabolic syndrome.

Reduction of A. muciniphila relative abundance in the gut microbiota is a widely accepted signature associated with obesity-related metabolic disorders. Using untargeted metabolomics profiling of fasting plasma, our study aimed at identifying metabolic signatures associated with beneficial properties of alive and pasteurized A. muciniphila when administrated to a cohort of insulin-resistant individuals with metabolic syndrome. Our data highlighted either shared or specific alterations in the metabolome according to the form of A. muciniphila administered with respect to a control group. Common responses encompassed modulation of amino acid metabolism, characterized by reduced levels of arginine and alanine, alongside several intermediates of tyrosine, phenylalanine, tryptophan, and glutathione metabolism. The global increase in levels of acylcarnitines together with specific modulation of acetoacetate also suggested induction of ketogenesis through enhanced β-oxidation. Moreover, our data pinpointed some metabolites of interest considering their emergence as substantial compounds pertaining to health and diseases in the more recent literature.

Introduction

The past two decades of research in the field of obesity-associated comorbidities were marked by the recognition of the gut microbiota as a key factor in disease etiology.[1] Efforts to characterize the human gut microbiota in health and disease, encompassing extensive translational research in rodent models, have led to the identification of the Gram-negative bacterium Akkermansia muciniphila (A. muciniphila) as a next-generation beneficial microbe.[2] We demonstrated that the bacterium exerts anti-obesogenic and health-promoting effects in mice,[3] which was subsequently widely confirmed by numerous publications and extended to other metabolic disorders (for the latest reviews see[4-7]). Following the unexpected discovery of pasteurization as a way to amplify its benefits in a similar murine model,[8] we implemented the first human pilot intervention, called Microbes4U©, consisting of daily ingestion of either alive or pasteurized A. muciniphila for 12 weeks. The initial primary objectives of the Microbes4U© study were to evaluate the tolerance, safety, and feasibility of A. muciniphila administration in a population with excess body weight suffering from prediabetes and metabolic syndrome. Not only we confirmed these aims but our data also showed several proofs of efficacy of the bacterium with stronger effects of the pasteurized form.[9] Beneficial impacts were shown on the levels of blood lipids, glycemic parameters such as insulin resistance, hepatic status with reduced blood levels of liver enzymes, and endotoxemia. We also observed some interesting trends regarding anthropometric measurements associated with obesity. Further in-depth toxicological analysis in rats confirmed the safety of the administration of the pasteurized A. muciniphila.[10] Although many metabolic outcomes and potential causal factors were identified in preclinical studies, some pieces of machinery remain to be uncovered. The collection of biological samples during the study allowed us to further explore the mechanisms underlying host–A. muciniphila crosstalk in the context of the Microbes4U© study. We notably highlighted that beneficial effects of A. muciniphila were independent of any major shift of the circulating endocannabinoidome. However, we identified a potential positive influence on the plasma levels of peroxisome proliferator-activated receptor alpha (PPARα) agonist 2-palmitoyl-glycerol.[1]

As A. muciniphila shows such a large effect on a variety of parameters, it is likely to influence a large series of other bioactive molecules. Hence, we were interested in carrying out a metabolomic study on blood samples to identify potential pathways affected by the bacterium as well as novel metabolite signatures associated with the metabolic effects of A. muciniphila. In this placebo-controlled interventional study conducted in treatment-naive subjects with prediabetes or metabolic syndrome, we have shown that A. muciniphila administration promotes differential and common alterations of the metabolome in individuals according to the form of the administered bacterium, i.e. the pasteurized versus the alive form.

Results

Pasteurized and alive A. muciniphila effects are linked to distinct and shared modulations of the serum metabolome

Because this work was initiated in the continuity of what we previously published, the approach adopted for the analysis of the metabolome is similar to what has been described previously.[1,9] Hence, for each metabolite, we measured the interventional effect per group, by calculating the mean differential value between the two main time-points, that is before and after 3 months of supplementation with the bacterium or the placebo (i.e., T0 and T3). The numbers obtained for both treatment groups were then subtracted from the ones of the placebo group to get the mean difference from placebo for each metabolite. We then assessed metabolites that were specifically and significantly modulated by the intervention through a comparison of differences between treatment groups and the placebo group (Figure 1). The two volcano plots displayed in Figure 1 provide an overview of the differential response between the two treatments at the level of the metabolome.

Figure 1.

Differential and common metabolomic effects of both treatments relative to the control group. (a) Study design of the Microbes4U© intervention. 52 obese or overweight individuals diagnosed with a pre-diabetic state and a metabolic syndrome were recruited and started the intervention. After three months of supplementation, blood samples were collected from 32 individuals. (b-c) Volcano plots depicting the metabolites which most notably changed following 3 months of supplementation with pasteurized (b) or alive A. muciniphila (c), under comparison with the placebo interventional effect. Labeled points represent metabolites for which the results of the Mann–Whitney U test comparing deltas of the treatment groups versus the placebo group were significant. Metabolites marked in red, exhibiting positive value for ‘mean difference from placebo’ are compounds that increased significantly following treatment with A. muciniphila. Metabolites marked in blue, exhibiting negative value for ‘mean difference from placebo’ are compounds that decreased significantly following treatment with A. muciniphila. The compounds which did not satisfy the condition of p < .05 for the univariate Mann–Whitney U test appear in yellow, and are not labeled. For readability, metabolites whose global delta exceeded the absolute value of 4, and for which the – result of the Mann–Whitney U test was not significant, were omitted from the plot. Metabolites for which the – log10Pvalue was below 0.5 were also discarded from the plot. The two dash lines correspond to p-value cutoffs of 0.05 and 0.01, respectively. See also Table S1 for the detailed list of the metabolites, with corresponding value and pathway

Compared to the placebo group, a total of 54 and 41 metabolites were identified as specifically and significantly modified following 3 months of supplementation with alive or pasteurized A. muciniphila, respectively. Briefly, 9 metabolites increased and 32 decreased in the pasteurized group, while 11 increased and 43 decreased in the alive group. The list of metabolites, with details of the corresponding pathways, is provided in Supplementary Tables 1 and 2. Of note, several metabolites related to the oxidative state were modulated. Cysteinylglycine (cys-gly) and oxidized cys-gly, two glutathione cycle intermediates, were significantly decreased by the alive from compared to the placebo group. Moreover, plasma levels of cysteine glutathione disulfide were significantly reduced in the alive group compared to baseline value (matched paired t-test, pv = 0.04). Global screening reveals differential modulation of the circulating metabolome between both treatment groups, with distinct signatures according to the form administrated. Nevertheless, it is worth underlining that alanine, arginine, caproate, 1-carboxyethylphenylalanine, and 3-hydroxyhexanoylcarnitine are among the metabolites that were altered in the same way regardless of the form administered. In conclusion, untargeted metabolomic profiling reveals both shared and distinct effects of alive and pasteurized forms on the serum metabolome in individuals with excess body weight and insulin resistance.

Table 1.

Influence of A. muciniphila on the metabolome is closely linked to ketone body-related metabolism

Top metabolites correlating most to acetoacetate at baseline (n = 52 biological samples). Metabolites for which the percentage of detection was below 100% percent are shown in red. A green background was used to highlight metabolites belonging to the lipid super pathway. Metabolites of the glycolytic pathway are shown on a blue background. For each metabolite, listed values for the corresponding Spearman’s coefficient of correlation, the fold change ratio, and the mean difference from placebo according to the treatment group are indicated. Adjustment for multiple testing was performed using the Benjamini-Hochberg method. Only the metabolites that significantly correlated with acetoacetate are presented (Adjusted p value >0.05). “T0” refers to the median-scaled baseline value, while “T3” refers to median scale final value. Abbreviations: A, Alive; N, Non-treated (Placebo); P, Pasteurized.

A. muciniphila negatively modulates circulating levels of several amino acid-derived metabolites potentially associated with hepatic function

Comparison of both volcano plots led to the observation that a substantial number of metabolites belonging to the amino-acid super pathway were significantly affected by the intervention; irrespective of the form, suggesting a common path induced by both forms of the bacterium. More specifically, several intermediates of tyrosine, phenylalanine, and tryptophan metabolism arose in the negative frame of both volcano plots. Based on this observation, we laid out pathway mapping for the quantified amino acids-related intermediates, adding annotations for the effect of the intervention (Figure 2). Briefly, this mapping aimed at summarizing all outputs from statistical analysis seeking for differences between groups (placebo versus treated) and within-groups (T0 versus T3) in serum metabolome. Exploration of the tyrosine and phenylalanine pathways showed a significant increase of 4-hydroxyphenylpyruvate, 4-hydroxyphenylacetate, 1-carboxyethylphenylalanine, N-acetyl-phenylalanine, and N-acetyl tyrosine in the placebo group (Figure 2a). Conversely, most intermediates, if not significantly decreased, tended to be reduced in individuals treated with A. muciniphila compared to baseline value and/or to the placebo group (Figure 2a). Similarly, although the effects were mostly dominated by trends when the markers were considered individually, the analyses of the tryptophan pathways showed an overall downregulation of the serum levels of several intermediates following both interventions compared with the placebo effect (Figure 2b). Collectively, administration of A. muciniphila seemed to alleviate the elevation observed in the serum levels of diverse intermediates of tyrosine, phenylalanine, and tryptophan metabolism occurring over time in treatment-naïve overweight subjects with prediabetes. Because numerous elements from the literature suggest a link between elevated amounts of tyrosine and phenylalanine-derived compounds and severity of liver dysfunction,[11-14] bivariate correlational analysis was performed on all subjects at baseline to assess relationships between serum concentration of hepatic enzymes, including transaminases, and quantified metabolites related to tyrosine and phenylalanine pathway in our cohort. The obtained Spearman’s correlation matrix indicated that both amino acids and several derived compounds were significantly and positively correlated with the serum levels of hepatic enzymes (Figure 2c). Conversely, three metabolites showed an inverse correlation with some enzymes. More specifically, phenol sulfate correlated negatively with aspartate aminotransferase (AST), phenylacetate correlated negatively with gamma-glutamyl transferase (GGT) and 2-hydroxyphenylacetate correlated negatively with alkaline phosphatase (AlKP). Interestingly, those metabolites were increased in response to the intervention while being decreased in the placebo group (Figure 2a).

Figure 2.

A. muciniphila negatively modulated circulating levels of several amino acids-derived metabolites that correlated with hepatic dysfunction Qualitative changes in phenylalanine, and tyrosine (a), and tryptophan (b) metabolites in plasma following the intervention. Intermediate metabolites which were not assayed are shown in gray and written in a smaller font. The routes including several steps and intermediaries not shown are represented by a dashed line. Enzymes are omitted. Compounds shown in red had a percentage of detection below 100% in our cohort. The trio of boxes below each metabolite represents the three groups, from left to right, given placebo, pasteurized, and alive A. muciniphila. Matched-pairs t-tests were performed on median scaled, log-transformed data, to verify changes from baseline (intragroup changes). Boxes with filled background mark metabolites that showed significant changes (p ≤ 0.05). The color represents the mean change from baseline value per group: red and green indicate a significant mean increase or decrease following intervention for the corresponding metabolite, respectively. Mann–Whitney U tests were performed to compare the differential values of both treated groups versus the placebo group (intergroup changes). Light red – and light green-shaded cells indicate 0.05 < p-value<0.10. The box was left emptied when the test was not significant. The direction of the arrow below both intervention groups indicates whether the metabolite had globally increased or decreased compared to the placebo effect. When the Mann–Whitney U test was significant (p ≤ 0.05), the arrow was colored according to the direction of the global change (red: increase, green, decrease). Light red – and light green-shaded cells indicate 0.05< pvalue<0.10. (c) Baseline Spearman’s correlation matrix between plasmatic hepatic enzymes and metabolites of the tyrosine and phenylalanine metabolism in the Microbes4U© cohort (n = 52). Negative correlations are colored in shades of blue and positive correlations in red. *: p < .05; **: p < .01; *** : p < .001

A. muciniphila-induced metabolome changes are closely linked to ketone bodies-related metabolism, including acylcarnitines

Impacts of A. muciniphila daily administration on the serum metabolite profiles encompassed a significant increase of acetoacetate (Figure 1). Given the numerous similarities existing between the beneficial effects driven by A. muciniphila and ketone bodies, we further explored the relationship existing between ketone bodies and other metabolites in our cohort by performing multiple correlation analyses. Table 1 lists the top metabolites that correlated significantly with acetoacetate at baseline. The table also intends to show how the metabolome associated with ketone bodies evolved during the intervention. Perturbations observed in the listed metabolites were quite similar between both A. muciniphila treatments compared to placebo. Moreover, although modest for most of them, modifications were often closely proportional to the value of the correlation coefficient. In other words, when a metabolite correlated positively with acetoacetate at baseline, then this metabolite increased following treatment, regardless of the form of administration, and inversely. For instance, pyruvate and lactate, two endproducts of the glycolytic pathway, not only correlated negatively with acetoacetate at baseline but also decreased in the serum as a result of both interventions. Noteworthily, a large majority of the top metabolites belong to the lipid super pathway, and among those, many are acylcarnitines. This prompted us to look at the profile of all acylcarnitines assayed in our study. We observed that most of them increased following both treatments (Figure 3a). Higher serum levels of acylcarnitines could be a downstream consequence of enhanced fatty acid β-oxidation. Consistent with our previous correlational analysis, when performing principal component analysis (PCA), we observed that acetoacetate correlated strongly and significantly with the dimension one of the PCA that best resumes variance among acylcarnitines (Figure 3 B-C).

Figure 3.

A. muciniphila administration was associated with a positive modulation of acylcarnitine-circulating levels. (a) Bars represent the mean difference from placebo of the relative plasma concentrations of measured acylcarnitines for the groups given pasteurized and alive A. muciniphila. The color indicates the direction of the global delta. (b) Principal component analysis of variables for acylcarnitines and direction of vector for acetoacetate if applied to the acylcarnitine variables plot. (c) Scatterplot showing the linear relation between first principal components and acetoacetate

A. muciniphila cross-talks with ketogenesis, lipid metabolism, and glycolysis-related metabolites

To strengthen the association between ketone bodies, glycolysis, and acylcarnitines in the context of our intervention, we generated a Spearman correlation matrix, calculated between the differential value of relevant metabolites and A. muciniphila (Figure 4a). Since acylcarnitines represent markers of fatty acid efflux destined for β-oxidation, we also considered the plasma concentration of non-esterified fatty acids (NEFA). Figure 4b is the direct representation of the matrix in a network format. Ketone bodies correlated positively with acylcarnitine and NEFA, suggesting that ketogenesis was enhanced alongside an increased transport of lipids. This cluster of metabolites correlated negatively with lactate and pyruvate, suggesting that ketogenesis preferentially used as substrate the acetyl-CoA produced by β-oxidation rather than glycolysis. Finally, although non-significant, except for NEFA, A. muciniphila correlated positively with elements relative to ketogenesis and fatty acids transport, and negatively with elements belonging to the glycolysis.

Figure 4.

A. muciniphila’s cross-talk with ketogenesis, lipid metabolism and glycolysis related metabolites. (a) Spearman’s rank correlation matrix of A. muciniphila and plasma metabolites belonging to ketone bodies metabolism, TCA cycle, and lipid metabolism, expressed as differential value (delta) (*P < .05). (b) Correlation Network map illustrating the Spearman’s correlation matrix. Metabolites that were highly correlated are brought together. The positioning of the metabolites was calculated by multidimensional scaling of the absolute values of the correlations. Gradient color, distance, and thickness of the lines were applied to metabolites nodes depending on coefficients of correlation. Negative and positives correlations are colored in shades of red and blue, respectively. Abbreviations: Akk, A. muciniphila; BHB, 3-hydroxybutyrate; NEFA, non-esterified fatty acids; X3.HHC, 3-hydroxhexanoylcarnitine

Discussion

Observations accumulated from both correlational studies in humans and preclinical studies in mice strongly positioned A. muciniphila as a “next-generation beneficial microbe”.[2] While many questions remained unanswered, previous and ongoing studies sought to shed light on the mechanisms behind the promotion of healthy status by A. muciniphila. In the framework of the Microbes4U© study, we conducted a large prospective serum metabolomic analysis to identify novel metabolite signatures connected with the beneficial effects of 3 months daily supplementation with either alive or pasteurized A. muciniphila to humans in the context of metabolic syndrome.

Following the characterization of metabolomic profiles of our cohort, we conducted an intensive search through the literature for the key metabolite discriminating the placebo effect from treatment effect, to uncover its connection with human health. Myo-inositol and phosphate were both significantly increased in response to pasteurized A. muciniphila exposure in comparison with the placebo. Administration of myo-inositol is associated with positive impacts in term of cardiovascular risk factors and diabetes incidence, notably through improvement of insulin sensitivity, in various pathological conditions such as childhood obesity, gestational diabetes, or polycystic ovary syndrome.[15-17] Moreover, A. muciniphila was strongly positively correlated with myo-inositol in a depressive rat model exhibiting reduced intestinal amount of the bacterium.[18] Notably, two bacteria (Dysosmobacter welbionis and Anaerostipes hadrus) recently identified as potentially beneficial against obesity and diabetes and also linked to higher thermogenesis, were shown to use myo-inositol to produce butyrate or propionate.[19,20] Regarding phosphate, two studies have described positive associations between a low phosphorus status and a higher risk for incident obesity.[21,22] Our results showed decreased levels of α-hydroxycaproate (Sup. Table 1) and 1-palmitoyl-2-stearoyl-GPC in the group given the pasteurized A. muciniphila. Interestingly, α-hydroxycaproate predicted cardiovascular mortality in atherosclerotic patients,[23] while 1-palmitoyl-2-stearoyl-GPC was found to be elevated in adipose tissues from obese individuals, in the plasma of dyslipidemic mice, and to be linked with obesity-cardiac lipotoxicity.[24-26] Collectively, these findings echo the widely accepted anti-atherogenic and anti-obesogenic properties of the bacterium (reviewed in[27]).

Major changes in the metabolome of the participants treated with the living form included down-modulation of caproate and 2-methylbutyrylglycine. Medium-chain fatty acid caproate is well documented in the literature as positively associated with type 2 diabetes (T2D), inflammatory bowel disease (IBD), rheumatoid arthritis, or Hepatitis B virus-induced cirrhosis.[28-31] Interestingly, one common feature for all of those conditions is a poor representation of A. muciniphila.[5,32,33] In the same line, A. muciniphila has been suggested to exert neuroprotective effects[2,5] and it is worth underlining that 2-methylbutyrylglycine is a metabolite described as a neurotoxin and that A. muciniphila treated subjects had lower levels of this metabolite.[34] Effects of living A. muciniphila were also characterized by an elevation of circulating levels of gluconate and 3-hydroxyhexanoylcarnitine. Although largely described in microbiology, the function of gluconate in human physiology remains unclear.[35] However, recent metabolomic studies highlighted a potential role in redox homeostasis through promotion of GSH recycling, in accordance with its metabolomic connection to the pentose phosphate pathway.[36] Finally, consistent with A. muciniphila being a signature of longevity and being underrepresented in age-related diseases,[37,38] 3-hydroxyhexanoylcarnitine was described as a potent predictor of healthy aging.[39]

Our analytic approach suggested some divergences in the mode of actions between living and pasteurized A. muciniphila. However, the metabolomics response triggered by both forms shared some striking similarities, such as a down-regulation of tyrosine, phenylalanine, and tryptophan metabolism. Based on these findings, a recent metagenome-wide association study showed that A. muciniphila was negatively associated with circulating levels of phenylalanine in a human cohort comprising lean and obese individuals.[40] In the same study, a large number of microbial species overrepresented in obese individuals were positively correlated with tyrosine. Noteworthily, as indicated by several human studies, tryptophan, tyrosine, and phenylalanine metabolisms are connected positively with body mass index, metabolic syndrome score, and T2D risk.[41-44] Some reports even underscored strong links between an altered tyrosine or phenylalanine metabolism and premalignant liver diseases, acute or chronic liver failure, and related-mortality.[11-14] More importantly, these observations were linked to the gut microbiota. Along with this, we further speculate that alterations in the aforementioned amino acids metabolism might be linked to the hepatoprotective effect of the bacterium. Spearman correlation analysis confirmed the existence of a global link between biomarkers of hepatic dysfunction and relative elevations of intermediates of tyrosine and phenylalanine metabolism. These findings reflect our previous observations of significantly reduced levels of transaminases and gamma glutamyl transferase in participants from the same cohort treated with the pasteurized A. muciniphila.[9] Of note, bariatric surgery, which leads to sustained weight loss, improved glycemic control and a bloom in gut A. muciniphila,[40,45-47] also results in down-modulation of tryptophan and intermediates of the pathway, especially the kynurenine pathway, which in our study was attenuated.[48] This being said, one exception relates to indolepropionate, a gut-derived product of tryptophan metabolism and a largely recognized health-associated metabolite, the reduced blood level of which constitutes a hallmark of various metabolic diseases.[49,50] In our studies, the serum of participants treated with pasteurized A. muciniphila showed reduced levels of indolepropionate. Whether this feature is unbeneficial is a relevant question, but given that all the other metabolites are pointing in another direction, this finding warrants further investigation.

Tyrosine, tryptophan, and phenylalanine are amino acid required in gut protein putrefaction, a process that is associated with disruption of gut homeostasis,[51] which is affected by bariatric surgery.[52,53] Thus, our result could potentially reflect a decrease of protein putrefaction. Except for hydroxyphenylacetate, typical putrefaction-derived metabolites were not significantly modulated by the intervention (i.e. putrescine, p-cresol, phenylacetylglutamine, spermidine).[54] However, that our results might related on modulation of gut protein or putrefaction remains suppositional without fecal analysis, since blood metabolome does not always reflect colon luminal metabolome.

Our data indicated significant reduction of alanine and arginine irrespective of the form of the administered A. muciniphila. Both amino acids were shown to be highly correlated with insulin resistance in a cohort of obese children.[55] Connections of these amino acids to metabolic unwellness were further confirmed in other human studies, including a recent systemic review of human studies addressing metabolomic profiling of obesity.[43,44,56,57]

Despite being a heterogeneous group of compounds, acylcarnitines were markedly increased in the serum of both treatment groups. Levels of most acylcarnitines naturally rise in the peripheral blood upon fasting condition or following exercise, as a reflect of increased transport of fatty acids, and subsequent increase in fatty acid β-oxidation rate.[58,59] Thereby, our observations are likely to be due to increased mitochondrial efficiency for fatty acid β-oxidation. In line with this, a recent in vivo study in mice confirmed that A. muciniphila inoculation not only promoted efflux of fatty acids through positive modulation of genes involved in hepatic lipid transport, but also resulted in a higher number of mitochondria.[60] Our hypothesis is further substantiated by recent research, including studies performed in our lab, showing that direct inoculation or indirect diet-mediated A. muciniphila enrichment in rodents activates a network of genes involved in lipid oxidation in adipose and liver tissue.[3,60-67] These include gene encoding for PPARα, peroxisome proliferator-activated receptor-gamma coactivator-1 alpha (PGC1α), and carnitine palmitoyl‐transferase I (CPT‐I). PPARα is a key sensor of fatty acid flux whose activation, together with PGC1α, stimulated the expression of CPT-I isoform α.[68] This enzyme catalyzes the formation of acylcarnitines by condensation of activated fatty acids (acyl‐CoAs) to carnitine, thus regulating the entry of acyl‐CoAs into the mitochondrial matrix.[69] Of note, we have recently shown in the Microbes4U© cohorts that A. muciniphila supplementation induced an increase blood level of 2-palmitoyl glycerol, a newly identify PPARα agonist.[1]

A metabolic switch toward the use of fatty acid as the major energy substrate was suggested to potentially result in higher oxygen consumption and energy expenditure.[70] Of note, we have previously described that administration of pasteurized A. muciniphila in a diet-induced obesity model led to increased energy expenditure, oxygen consumption, and carbon dioxide production in a thermogenesis-independent manner.[71]

If not oxidized in the TCA cycle, as suggested by reduced amount of fumarate and malate, an alternative route for β-oxidation-derived Acetyl-coA is ketogenesis. Our results indicated a significant increase of acetoacetate in the group given live A. muciniphila. We also noted a similar trend for 3-hydroxybutyrate. Of note, both ketone bodies tend to increase in the group given pasteurized A. muciniphila. Finally, we demonstrated that the metabolites closely associated with ketone bodies shifted toward the direction of that correlation upon both interventions. Recent literature tends to emphasize the key role of ketone bodies as signaling molecules (for review see[72]). Ketogenesis can be induced by the consumption of high-fat, low-carbohydrate diets.[73] The so-called ketogenic diet is currently getting popularity as a new efficient approach for obesity management since it promotes weight loss, glycemic control, reduced systemic inflammation among other effects.[74] Positive effects were also well described in the context of epilepsy,[75] a condition associated with an altered gut microbiota.[76] Moreover, colonization with A. muciniphila exerts anti-seizure effects in a mice model.[77] Adding to the physiological importance of the pathway, impaired ketogenesis was shown to exert a strong impact on steatosis development and to limit hepatic gluconeogenesis in humans.[41] As for our result regarding amino-acid metabolism, this is also in line with numerous murine studies showing that A. muciniphila lowers hepatic inflammation and steatosis.[60,78-80]

This proposed mechanism of enhanced ketogenesis is supported by several lines of evidence. Indeed, Crohn’s disease, ulcerative colitis, multiple sclerosis, Alzheimer’s disease, and autism are among pathological conditions for which ketogenic diet-based therapies were reported as beneficial.[81] Interestingly, those are conditions in which intestinal A. muciniphila is underrepresented, but also those in which A. muciniphila administration in murine models improves the phenotype,[82] as reviewed in.[5,8,79]

On top of this, ketone bodies and A. muciniphila share common features in the context of healthy-aging and cancer prevention. This fact was notably highlighted in a recent metabolomic study carried out in mice receiving alive or pasteurized A. muciniphila. The authors described an elevation of intestinal and circulating levels of the ketone body 2β-hydroxybutyrate, and linked this to the anti-aging and anti-cancer effect of A. muciniphila.[83] Along the same line, colonization of germ-free mice with A. muciniphila led to significant modulation of genes related to ketone body metabolism in the ileum.[84] Finally, one study also showed increased plasma acetoacetate in chow-fed mice treated with live A. muciniphila.[85]

Consistent with their implications in the maintenance of cellular redox homeostasis, ketone bodies, together with the gluconate metabolite, might be partially responsible for the observed lower levels of indirect markers of oxidative stress (i.e. cys-gly, oxidized cys-gly, cysteine glutathione disulfide). Favorable redox state, in turn, preserves the function of the carnitine/acylcarnitine carrier involved in entry of fatty acids in the mitochondria.[86]

Our assumptions need further validation in larger prospective human studies. However, when all the observations described in this paper are grouped, a metabolic pattern emerges. A proposed metabolic mechanism for the plasmatic metabolome associated with A. muciniphila exposure is displayed in Figure 5.

Figure 5.

Overview of the mechanism proposed for the A. muciniphila-mediated metabolomic switch toward β-oxidation and ketogenesis. In physiological conditions, overnight fast is associated with high lipolytic rates to increase the availability of substrates for β – oxidation. This results in a substantial efflux of non-esterified fatty acids and esterified carnitine in the plasma. Entry of acyl-CoA into the mitochondria is facilitated by the acylcarnitine transport system, while part of the generated acylcarnitines are released in the plasma. Fatty acids are then broken down into Acetyl-coA through β-oxidation in the mitochondria. Administration of A. muciniphila in insulin-resistant overweight individual turndowned glycolysis and amplified the aforementioned pathway and triggered a switch toward the use of Acetyl-CoA for ketogenesis rather than toward the TCA cycle. In turn, elevated levels of ketones bodies, alongside gluconate, may contribute to reduce the redox state, favoring the activity of CACT. The thicker the arrow the stronger the corresponding pathway was triggered. Elements and arrows in pink reflect known facts from literature relative to A. muciniphila action. Increase and decrease metabolites are shown in blue and green, respectively. CACT, Carnitine-acylcarnitine translocase; CPT1α, carnitine palmitoyl‐transferase I; G6P, Glucose-6-phosphatase; PPARα, peroxisome proliferator-activated receptor alpha; TCA, tricarboxylic acid cycle

Briefly, our data raise the hypothesis that administration of A. muciniphila in prediabetic subjects with metabolic syndrome enhances hepatic fatty acids delivery and subsequent β-oxidation, but attenuates acetyl-CoA oxidation in the TCA cycle, leading preferentially to ketogenesis. Downstream metabolomic effects of A. muciniphila also lead to downregulation of the metabolism of several amino acids, whose implications in various metabolic disorders are convincingly demonstrated in the scientific literature. We notably identified a global downregulation of the tyrosine and phenylalanine metabolism as a potential mechanism contributing to the hepato-protective effects of A. muciniphila. Our analysis also revealed specific modulations of metabolites with emerging implications in health and diseases, including disorders for which A. muciniphila was shown to be negatively associated with. To close the discussion, it is noteworthy to raise the question of whether or not our observations are the cause or the consequence of the improvement of the metabolic state, constituting thus a limitation of our study. Although we acknowledge the modest scope of each observation taken individually, the strength of our study lies in the coherence between the global scheme that we propose, born from the junction between all observations, and the existing literature in relation to A. muciniphila. Those new insights bridge the recognized links between obesity-associated metabolomic hallmarks (including oxidative stress, alteration of amino acids metabolism, markers of hepatic disorders) and isolated observations connecting A. muciniphila to its health-promoting status.

Materials and methods

Participants and study design

The Microbes4U© study used a randomized, double-blinded, parallel, placebo-controlled design. The cohort consisted of overweight or obese individuals (body mass index (BMI)>25 kg/m2), aged between 18 and 70 years. Subjects were eligible for the study if they present with a newly diagnosed pre-diabetic state and a metabolic syndrome. Pre-diabetes was evaluated using HOMA-modeling of insulin sensitivity (HOMA Calculator, University of Oxford), with insulin resistance considered if participant’s insulin sensitivity was below 75%. The diagnosis of a metabolic syndrome was established according to the National Cholesterol Education Program Adult Treatment Panel III definition, that is, at least three of the five following criteria: fasting glycemia >100 mg dL-1; blood pressure >130/85 mmHg or antihypertensive treatment; fasting triglyceridemia >150 mg dL-1; high-density lipoprotein (HDL) cholesterol <40 mg dL-1 for men, <50 mg dL-1 for women; and/or waist circumference >102 cm for men, >88 cm for women. Exclusion criteria and details regarding study design were previously extensively described.[9] Briefly, between 2015 and 2018, 52 volunteers were enrolled at the Cliniques Universitaires Saint-Luc in Belgium and started the intervention. Following screening, participants were randomized into four parallels groups, receiving for 12 weeks either a sterile liquid solution of PBS-containing glycerol as a placebo, or an equivalent volume of the same solution plus 1010 alive A. muciniphila, 109 alive A. muciniphila, or 1010 pasteurized A. muciniphila. Pasteurization consisted of heat treatment at 70°C for 30 min of fresh A. muciniphila. Participants and physicians were both blinded to treatment allocation. Blood samples were collected in the fasting state at the beginning of the intervention prior to any administration and on the last day of the trial. The study protocol was approved by the local ethical committee on July 2015 (comité d’éthique hospitalo-facultaire UCLouvain, Cliniques Universitaires Saint-Luc) under the number 2015/02JUL/369, and the study was registered at clinicaltrial.gov under the number NCT02637115.

This study aimed at further extending our preliminary analysis, beyond safety, tolerability, and feasibility assessment.[8,9] Aiming will to continue the effort undertaken, on the one hand, in the characterization of the effects and, on the other hand, in the exploration of the underlying biological mechanisms, our metabolomic analysis focused on the comparison of the control group with the groups treated with the highest dose tested, that is 1010 bacteria per day. Accordingly, 32 subjects completed the trial after 3 months of daily supplementation. Baseline characteristics of that cohort were previously described.[9]

Blood Serum biochemical analysis

Plasma samples were collected after overnight fasting (8 hours minimum) in lithium-heparin coated tubes. One set of tubes was sent directly to the hospital laboratory for several blood analyses including liver enzymes activities. The remaining samples were transported on ice from the sampling center to the research laboratory. Plasma was immediately isolated from whole blood by centrifugation at 4200 g for 10 min at 4°C and stored at −80°C. For metabolomics profiling, 100 μl of plasma was aliquoted and transported on dry ice to Metabolon Inc.

Alkaline phosphatase (AlkP) and gamma-glutamyl transferase (GGT) were assayed by enzymatic colorimetric method (Cobas 8000 – Roche Diagnostics), while aspartate aminotransferase (AST) and alanine aminotransferase (ALT) were assayed by enzymatic methods (International Federation of Clinical Chemistry and Laboratory Medicine) without activation by pyridoxal phosphate (Cobas 8000 – Roche Diagnostics). Plasma non-esterified fatty acids (NEFAs) were measured using kits coupling an enzymatic reaction with spectrophotometric detection of the reaction endproducts (Diasys Diagnostic and Systems, Holzheim, Germany) according to the manufacturer’s instructions.

Untargeted metabolomics assays

Untargeted metabolomics profiling using high-resolution mass spectrometry with hydrophilic interaction chromatography was applied to all samples via Metabolon’s HD4 multi-platform techniques. The non-targeted mass spectrometry analysis of Metabolon Inc. (North Carolina, USA) was previously described.[87] Briefly, it comprised ultra-performance liquid chromatography/mass spectrometry with a heated electrospray ionization source and mass analyzer. Following proper handling, samples were first prepared using the automated MicroLab STAR® system from Hamilton Company. The resulting extract was divided into several fractions, analyzed by four ultra-high-performance liquid chromatography-tandem mass spectrometry according to the Metabolon pipeline. Biochemical identification of metabolites contained in one sample was then performed by comparison to a reference library of purified standards consisting of more than 33000 metabolites. The comparison was based on retention time/index, mass-to-charge ratio (m/z), and chromatographic data (MS/MS spectral data) using software developed at Metabolon. Further details regarding quality controls, data extraction, curation, quantification, and bioinformatics were previously described.[88]

Statistical analysis

Metabolomic profiling detected a total of 1169 compounds, 947 of which were of known identity, while 222 were of unknown identity (X-number). Because of a very weak detection rate, metabolites related to drugs (i.e. analgesic, neurological, respiratory, antibiotic, psychoactive) were not considered in the analysis. Prior to formal statistical analysis, each compound was scaled to the median value in order for the median value for each metabolite to equal one. Missing data were imputed with the minimum observed value for that compound. Evaluation of the interventional effect within each group was assessed using Matched paired T-test on log-transformed, median-scaled data (intragroup change). The ratio or fold change was calculated for each metabolite by dividing the final value by the baseline value, to evaluate the mean average individual changes between both timings, according to the following formula:

where ‘ind’ referred to each participant and ‘group’ to the group.

Similarly, the mean change from baseline to end of treatment was then calculated for each metabolite, within each group, by subtracting the value obtained at Time 0 from the value obtained at Time 3 months for each participant, according to the following formula:

By subtracting the mean difference, thus calculated for each treatment group, to the one calculated for the placebo group, we obtained the “mean difference from placebo” for each metabolite, according to the following formula:

This value gave a global measurement of the interventional effect when compared to the evolution of the placebo group. Mann–Whitney U tests were then performed to compare whether the mean difference obtained for the two treatment groups was significantly different from the one calculated for the placebo group. Intergroup statistical analyses were conducted using SPSS v.27.0 (IBM Corporation). All tests were two-tailed and significance was set at p < .05.

Results from those extended univariate analyses were collected and presented using two volcano plots, one per treatment group. X-axis and Y-axis represent the mean difference from the placebo, and statistical significance as the negative log of Mann–Whitney P-values, respectively. The significantly increased or decreased metabolites are labeled in red or blue, respectively, while non-significant genes are shown as color dots without labeling.

Spearman’s correlation of acetoacetate against all others metabolites of the dataset was computed on Rstudio using the package ‘Hmisc’ (version 4.5–0) with multiple testing correction via FDR estimation according to the Benjamini and Hochberg procedure. The metabolites that significantly correlated with acetoacetate and with the highest absolute coefficient of correlation were extracted and listed in a table format.

Principal component analyses were performed on Rstudio (R version 3.6.3, Rstudio Team, Boston, MA, USA) using factoextra (version 1.0.7) and factoMiner (version 2.3) packages. Linear regression was then built between the metabolite acetoacetate and the first principal component summarizing the highest variance in acylcarnitine, as previously described.[89]

The Spearman’s correlation matrix related to ketone bodies was drawn with Prism software v.9.0 (GraphPad Software). All the other plots (i.e. Volcano plots, Spearman’s correlation matrix, Barplots, the correlative network plot, and the Scatterplot) were constructed on Rstudio (R version 3.6.3, Rstudio Team, Boston, MA, USA) using the R packages “ggplot2” (version 3.3.2), “tidyverse” (version 1.3.0) and “corrr”(version 0.4.3).

Click here for additional data file.