Probiotic modulation of symbiotic gut microbial-host metabolic interactions in a humanized microbiome mouse model.

The transgenomic metabolic effects of exposure to either Lactobacillus paracasei or Lactobacillus rhamnosus probiotics have been measured and mapped in humanized extended genome mice (germ-free mice colonized with human baby flora). Statistical analysis of the compartmental fluctuations in diverse metabolic compartments, including biofluids, tissue and cecal short-chain fatty acids (SCFAs) in relation to microbial population modulation generated a novel top-down systems biology view of the host response to probiotic intervention. Probiotic exposure exerted microbiome modification and resulted in altered hepatic lipid metabolism coupled with lowered plasma lipoprotein levels and apparent stimulated glycolysis. Probiotic treatments also altered a diverse range of pathways outcomes, including amino-acid metabolism, methylamines and SCFAs. The novel application of hierarchical-principal component analysis allowed visualization of multicompartmental transgenomic metabolic interactions that could also be resolved at the compartment and pathway level. These integrated system investigations demonstrate the potential of metabolic profiling as a top-down systems biology driver for investigating the mechanistic basis of probiotic action and the therapeutic surveillance of the gut microbial activity related to dietary supplementation of probiotics.

Introduction

The pan class="Species">gut microbiome–pan class="Species">mammalian ‘Superorganism' (Lederberg, 2000) represents a level of biological evolutionary development in which there is extensive ‘transgenomic' modulation of metabolism and physiology that is a characteristic of true symbiosis. By definition, superorganisms contain multiple cell types, and the coevolved interacting genomes can only be effectively studied as an in vivo unit in situ using top-down systems biology approaches (Nicholson, 2006; Martin ). Interest in the impact of gut microbial activity on human health is expanding rapidly and many mammalian–microbial associations, both positive and negative, have been reported (Dunne, 2001; Verdu ; Nicholson ; Gill ; Ley ). Mammalian–microbial symbiosis can play a strong role in the metabolism of endogenous and exogenous compounds and can also be influential in the etiology and development of several diseases, for example insulin resistance (Dumas ), Crohn's disease (Gupta ; Marchesi ), irritable bowel syndrome (Sartor, 2004; Martin ), food allergies (Bjorksten ), gastritis and peptic ulcers (Warren, 2000; Marshall, 2003), obesity (Ley ; Turnbaugh ), cardiovascular disease (Pereira and Gibson, 2002) and gastrointestinal cancers (Dunne, 2001). Activities of the diverse gut microbiota can be highly specific and it has been reported that the establishment of Bifidobacteria is important for the development of the immune system and for maintaining gut function (Blum and Schiffrin, 2003; Salminen ; Ouwehand, 2007). In particular, elevated counts in Bifidobacterium with reduced Escherichia coli, streptococci, Bacteroides and clostridia counts in breast-fed babies compared to formula-fed neonates may result in the lower incidence of infections, morbidity and mortality in breast-fed infants (Dai ; Kunz ). As the microbiome interacts strongly with the host to determine the metabolic phenotype (Holmes and Nicholson, 2005; Gavaghan McKee ) and metabolic phenotype influences outcomes of drug interventions (Nicholson ; Clayton ), there is clearly an important role of understanding these interactions as part of personalized healthcare solutions (Nicholson, 2006).

One of the current approaches used to modulate the bpan class="Chemical">alance of intestinal microflora is based on oral administration of probiotics. A probiotic is generally defined as a ‘live microbial feed supplement which beneficially affects the host animal by improving its intestinal microbial balance' (Fuller, 2004). The gastrointestinal system is populated by potentially pathogenic bacteria that are capable of degrading proteins (putrefaction), releasing ammonia, amines and indoles, which in high concentrations can be toxic to humans (Cummings and Bingham, 1987). Probiotic supplementation aims at replacing or reducing the number of potentially harmful E. coli and Clostridia in the intestine by enriching the populations of gut microbiota that ferment carbohydrates and that have little proteolytic activity. Probiotics, most commonly Lactobacillus and Bifidobacteria, can be used to modulate the balance of the intestinal microflora in a beneficial way (Collins and Gibson, 1999). Although Lactobacilli do not predominate among the intestinal microflora, their resistance to acid conditions and bile salts toxicity results in their ubiquitous presence throughout the gut (Corcoran ), hence they can exert metabolic effects at many levels. Fermented dairy products containing Lactobacillus have traditionally been used to modulate the microbial ecology (Dunne, 2001). In particular, L. paracasei was shown to modulate the intestinal physiology, to prevent infection of pathogenic bacteria (Sarker ), to stimulate the immune system (Ibnou-Zekri ), and to normalize gastrointestinal disorders (Martin ). L. rhamnosus is also a significant probiotic strain with proven health benefits and therapeutic applications in the treatment of diarrhea (Szynanski ), irritable bowel syndrome (Kajander ), atopic eczema (Corcoran ) and the prevention of urinary tract infections (Reid and Bruce, 2006). However, the functional effects of probiotic interventions cannot be fully assessed without probing the biochemistry of the host at multiple compartmental levels, and we propose that top-down systems biology provides an ideal approach to further understanding in this field. The microbiota observed in human baby flora (HBF) mice have a number of similarities with that found in formula-fed neonates (Mackie ), which makes it to be a well-adapted and simplified model to assess probiotics impact on gut microbial functional ecosystems (in particular on metabolism of Bifidobacteria and potential pathogens) and subsequent effects on host metabolism.

Metabolic profiling using high-density data genepan class="Species">rating spectroscopic techniques, in combination with multivariate mathematical modelling is a tool which is well suited to generate metabolic profiles that encapsulate the top-down system response of an organism to a stressor or intervention (Nicholson and Wilson, 2003). Multivariate metabolic profiling offers a practical approach to measuring the metabolic endpoints that link directly to whole system activity and which are determined by both host genetic and environmental factors (Nicholson ). Recently, metabolic profiling strategies have been successfully applied to characterizing the metabolic consequences of nutritional intervention (Rezzi ; Wang ) the effects of the gut microflora on mammalian metabolism (Martin , 2007a, 2007b) and mechanisms of insulin-resistance (Dumas ). In the current study, 1H nuclear magnetic resonance (NMR) spectroscopy and targeted ultra performance liquid chromatography-mass spectrometry (UPLC-MS) analysis have been applied to characterize the global metabolic responses of humanized microbiome mice subsequently exposed to placebo, Lactobacillus paracasei or Lactobacillus rhamnosus supplementation. Correlation of the response across multiple biofluids and tissue, using plasma, urine, fecal extracts, liver tissues and ileal flushes as the biological matrices for the detection of dietary intervention, generates a top-down systems biology view of the response to probiotics intervention.

Results

Gut bacterial composition

Microbiological analyses were performed on fecal samples to assess the growth of the HBF in germ-free pan class="Species">mice and to ascertain the effects of probiotics on the development of gut bacteria. The measured terminal composition of the fecal microbiota is detailed in Table I, where the statistically significant differences between the various groups were calculated using a two-tailed pan class="Species">Mann–Whitney test. The bacterial populations of Bifidobacteria longum and Staphylococcus aureus were reduced after introduction of both probiotics. Additionally, unique effects of L. rhamnosus supplementation caused decreased populations of Bifidobacterium breve, Staphylococcus epidermidis and Clostridium perfringens but an increase of E. coli.

Table 1

Microbial species counts in mouse feces at the end of the experiment

Groups/log10 CFU	HBF (n=10)	HBF+L. paracasei (n=9)	HBF+L. rhamnosus (n=9)
L. paracasei	—	8.5±0.2	—
L. rhamnosus	—	—	7.8±0.2
E. coli	9.2±0.3	9.4±0.3	9.8±0.5**
B. breve	9.1±0.2	7.78±2.13	8.7±0.3*
B. longum	8.2±0.6	5.6±1.9***	6.3±0.5***
S. aureus	7.4±0.3	6.3±0.3***	6.6±0.5***
S. epidermidis	4.8±0.4	4.9±1.2	4.0±0.5**
C. perfringens	7.2±0.3	7.0±0.5	5.7±1.0***
Bacteroides	10.3±0.2	10.4±0.2	10.1±0.4

log10 CFU (colony-forming unit) given per gram of wet weight of feces. Data are presented as mean±s.d. Absence of specific bacterial strains in the gut microflora is indicated by “—”. The values for the HBF mice supplemented with probiotics were compared to HBF control mice, *, ** and *** indicate a significant difference at 95, 99 and 99.9% confidence levels, respectively.

Gut levels of short-chain fatty acids

pan class="Chemical">Short-chain fatty acids (pan class="Chemical">SCFAs), namely acetate, propionate, isobutyrate, n-butyrate and isovalerate, were identified and quantified from the cecal content using GC-FID. The results, presented in Table II, are given in μmol per gram of dry fecal material and as mean±s.d. for each group of mice. The production of some of the SCFAs, that is, acetate and butyrate, by the HBF mice supplemented with both of the probiotics was reduced. In addition, increases of the concentrations in isobutyrate and isovalerate were observed in the mice fed with L. paracasei.

Table 2

Short-chain fatty acid content in the cecum from the different groups

Amounts of SCFAs given in μmol per gram of dry feces for each group	Acetate	Propionate	Isobutyrate	Butyrate	Isovalerate
HBF (n=10)	77.6±17.6	22.3±4.3	0.9±0.2	3±0.6	2.1±0.6
HBF+L. paracasei (n=9)	52.3±23.6***	22.2±10.8	1.2±0.5***	1.5±0.8***	2.7±1.2**
HBF+L. rhamnosus (n=9)	40.6±8***	20.3±2.8	0.8±0.2	2.1±0.4***	2.1±0.5

Data are presented in μmol per gram of dry feces and are presented as means±s.d. The amounts of SCFAs for the HBF mice supplemented with probiotics were compared to HBF control mice, ** and *** indicate a significant difference at 99 and 99.9% confidence levels, respectively.

Analysis of1H NMR spectroscopic data on plasma, urine, liver and fecal extracts

A series of pairwise O-PLS-DA models of 1H NMR spectra were performed to extract information on the metabolic effects of probiotic modulation. A statistically significant metabolic phenotype sepapan class="Species">ration between untreated pan class="Species">mice and probiotic supplemented animals was observed as reflected by the high value of QY2 for each model (Cloarec ; Table III). The corresponding coefficients describing the most important metabolites in plasma, liver, urine and fecal extracts that contributed to group separation are also listed in Supplementary Table 1. The area normalized intensities (101 a.u.) of representative metabolite signals are given as means±s.d. in Table III. The O-PLS-DA coefficients plots are presented in Figure 1 using a back-scaling transformation and projection to aid biomarker visualization (Cloarec ). The direction of the signals in the plots relative to zero indicates positive or negative covariance with the probiotic-treated class. Each variable is plotted with a color code that indicates its discriminating power as calculated from the correlation matrix thus highlighting biomarker-rich spectral regions.

Table 3

Summary of influential metabolites for discriminating NMR spectra of liver, plasma, fecal extracts and urine

Metabolites	Chemical shift and multiplicity	HBF controls	HBF+L. paracsei	HBF+L. rhamnosus
Liver			QY2=21%, RX2=44%	QY2=41%, RX2=37%
Leu	0.92(t)	2.4±0.6	1.7±0.3***	1.9±0.5*
Ileu	0.94 (t)	0.8±0.1	0.6±0.05***	0.7±0.2
Lactate	1.32(d)	38.4±5.8	46.2±8.2a	39.2±9.7
Succinate	2.41(s)	0.2±0.1	1.0±0.6**	0.7±0.3*
MA	2.61(s)	0.1±0.06	0.04±0.002**	0.08±0.05
TMA	2.91 (s)	0.2±0.04	0.07±0.03***	0.2±0.09
TMAO	3.27(s)	10.3±2.2	13.1±3.7	18.5±8.0**
Gln	2.44(m)	0.4±0.1	0.3±0.1*	0.3±0.1
Glycogen	5.38–5.45	3.4±1.9	1.5±0.6*	3.2±1.9
Plasma			QY2=44%, RX2=50%	QY2=51%, RX2=32%
Lipoproteins	0.84 (m)	13.7±1.8	10.1±0.9***	9.8±4.2**
Citrate	2.65(d)	1.4±0.3	0.9±0.4**	1.1±0.2**
Choline	3.2(s)	11.6±2.6	16.2±5.7*	20.5±3.8***
GPC	3.22(s)	44.1±4.6	57.3±12.5**	68.1±11.2***
Glyceryls	3.91(m)	2.0±0.3	2.5±0.4**	2.7±0.4**
Feces			QY2=90%, RX2=48%	QY2=89%, RX2=49%
Caprylate	1.27(m)	2.5±0.1	3.5±0.2***	2.4±0.1
Lys	3.00(m)	3.2±0.8	5.0±0.3***	4.9±1.2**
Osides	5.42(m)	0.9±0.09	1.2±0.1***	1.4±0.1***
Bile acids	0.72(s)	3.1±0.9	1.8±0.7**	2.0±0.6*
Ethanol	1.18(t)	2.5±0.1	2.0±0.08***	1.9±0.09***
Choline	3.20(s)	48.0±19.5	11.3±4.1***	20.3±10.9***
NAM	2.06(m)	7.1±1.0	5.4±0.3***	5.5±0.3***
Acetate	1.91(s)	58.7±34.2	27.0±9.2**	32.9±12.9*
U1	3.71(s)	9.2±0.5	7.3±0.3***	8.2±0.5***
Urine			QY2=91%, RX2=47%	QY2=59%, RX2=46%
IAG	7.55(d)	0.1±0.03	0.6±0.2***	0.4±0.2**
PAG	7.37(m)	0.8±0.1	1.5±0.3***	1.2±0.4*
Tryptamine	7.70(d)	0.1±0.04	0.4±0.1***	0.2±0.1**
UGLp	1.27(m)	1.7±0.1	2.7±0.4***	1.7±0.2
Glycero-metabolites	4.04 (m)	1.7±0.1	2.2±0.2***	1.9±0.2*
NAG	2.04(s)	3.5±0.2	4.3±0.2***	3.8±0.5
Butyrate	0.90(t)	6.9±0.8	5.2±0.7**	5.2±0.9**
α-keto-isocaproate	0.94(d)	13.8±4.6	6.1±2.3**	7.9±2.1**
Propionate	1.05(t)	0.9±0.2	0.8±0.04*	0.8±0.1
3-hydroxy-isovalerate	1.24(s)	3.0±0.4	2.1±0.4**	2.6±0.2
Citrate	2.55(d)	10.8±7.6	1.6±0.6**	2.2±0.9*
Creatine	3.92(s)	5.7±2.1	4.5±1.5	3.5±0.3**
Citrulline	1.88(m)	3.8±0.5	3.3±0.4	3.0±0.4**

O-PLS models were generated for comparing probiotics treated to HBF control mice using one predictive and two orthogonal components, RX value shows how much of the variation is explained, QY value represents the predictability of the models, and relates to its statistical validity. Data are presented as area normalized intensities (101 a.u.) of representative metabolite signals as means±s.d. The values for the HBF mice supplemented with probiotics were compared to HBF control mice, a, *, ** and *** indicate a significant difference at 90, 95, 99 and 99.9% confidence levels, respectively.

s, singlet; d, doublet; t, triplet; q, quartet; m, multiplet; dd, doublet of doublets. For metabolite abbreviations, refer to the key in Figure 1, MA: methylamine.

Figure 1

O-PLS-DA coefficient plots derived from 1H MAS NMR CPMG spectra of liver (A, D), 1H NMR CPMG spectra of plasma (B, E), 1H NMR standard spectra of fecal extracts (C, F) and urine (G, H), indicating discrimination between HBF mice fed with probiotics (positive) and HBF control mice (negative). The color code corresponds to the correlation coefficients of the variables with the classes. BAs, Bile acids; DMA, dimethylamine; Glc, glucose; Gln, glutamine; GPC, glycerophosphorylcholine; IAG, indoleacetylglycine; Ileu, isoleucine; Leu, leucine; Lys, lysine; NAG, N-acetylated glycoproteins; NAM, N-acetylated metabolites; Osides, glycosides; PAG, phenylacetylglycine; TBAs, taurine conjugated to bile acids; TMA, trimethylamine; TMAO, trimethylamine-N-oxide; UGLp, unidentified glycolipids.

Liver metabolic profiles

Livers of pan class="Species">mice fed with pan class="Species">L. paracasei showed relative decreases in dimethylamine (DMA), trimethylamine (TMA), leucine, isoleucine, glutamine, and glycogen and increased levels of succinate and lactate (Figure 1A). Mice supplemented with L. rhamnosus showed relative decreases in leucine and isoleucine and relative increases in succinate, TMA and trimethylamine-N-oxide (TMAO) in the liver compared to controls (Figure 1D).

Plasma metabolic profiles

Plasma samples showed relative decreases in the levels of lipoproteins and increases in the concentpan class="Species">rations of pan class="Chemical">glycerophosphorylcholine (GPC) and triglycerides in mice fed with both probiotics compared to controls (Figure 1B and E). Elevated choline levels were observed in plasma of mice fed with L. rhamnosus and reduced plasma citrate levels were observed in mice fed with L. paracasei compared to controls.

Fecal extract metabolic profiles

Marked changes were observed in the metabolic profiles of fecal extracts from all supplemented pan class="Species">mice, for example relative decreased concentrations of choline, acetate, ethanol, a range of putative N-acetylated metabolites (NAMs), unconjugated bile acids (BAs) and tauro-conjugated bile acids (Figure 1C and F). Furthermore, relative higher levels of glucose, lysine and polysaccharides were detected in the feces from mice fed with probiotics. A relative increased level of n-caproate (chemical shifts δ at 0.89(t), 1.27(m), 1.63(q), 2.34(t)) appeared to be associated with mice supplemented with L. paracasei.

Urine metabolic profiles

Urine samples of pan class="Species">mice supplemented with both probiotics showed relative increased concentrations of indoleacetylglycine (IAG), phenylacetylglycine (PAG), tryptamine and a relative decrease in the levels of α-keto-isocaproate and citrate (Figure 1G and H). Relative increased concentrations of a mixture of putative glycolipids (UGLp, chemical shifts of multiplets at δ 0.89, 1.27, 1.56, 1.68, 2.15, 2.25, 3.10, 3.55, 3.60), N-acetyl-glycoproteins (NAGs) and a reduction in 3-hydroxy-isovalerate were also observed in mice supplemented with L. paracasei compared to controls. Urine of mice fed with L. rhamnosus showed a reduction in levels of creatine and citrulline.

UPLC-MS analysis of bile acids in ileal flushes

The proportion of pan class="Chemical">bile acids in ileal flushes from the different groups are given in Table IV and are shown as mean±s.d. of the percentage of the total bile acid content. O-PLS-DA of the data set revealed that the relative concentrations of bile acids obtained from unsupplemented HBF mice are separated from those treated with probiotics, the correlation observed with L. paracasei being more significant than with L. rhamnosus as noted by the values of the cross-validated model parameter QY2 (Figure 2). For example, HBF mice supplemented with L. paracasei showed strong correlations with higher amounts of GCA, CDCA and UDCA and lower levels of α-MCA in the ileal flushes when compared to controls. HBF mice fed with L. rhamnosus also showed higher levels of GCA associated with lower levels of TUDCA and TCDCA in the ileal flushes when compared to untreated HBF mice.

Table 4

Bile acids composition in gut flushes for the different microbiota

Microbiota/ Bile acids	HBF	HBF+L. paracasei	HBF+L. rhamnosus
CDCA	ND	0.04±0.07	0.01±0.02
UDCA	ND	0.1±0.1	ND
CA	0.4±0.5	0.3±0.7	0.3±0.2
ωMCA	ND	0.02±0.07	ND
αMCA	0.3±0.2	0.3±0.6	0.2±0.2
βMCA	0.9±0.7	1.3±2.7	0.8±0.7
GCA	ND	0.1±0.1	0.1±.08
TCDCA	3.3±1	4±2.1	2.5±1.3
TUDC A	6.6±1.4	6.6±5.2	5.1±1.02
TβMCA	50±4.8	49.2±7.5	50.5±5.2
TCA	38±3.5	38±6.2	40.4±5

Relative composition in bile acids given in percentage of total bile acid content. Species not detected with UPLC-MS experiment are shown as ND. The key is given in UPLC-MS material and methods.

Figure 2

O-PLS-DA coefficient plots derived from the bile acid composition obtained by UPLC-MS analysis of ileal flushes, which indicate discrimination between HBF control mice (negative) and HBF mice treated with probiotics (positive), (A) L. paracasei and (B) L. rhamnosus. The color code corresponds to the correlation coefficients of the variables. One predictive and one orthogonal component were calculated; the respective QY2 and RX2 are (76.4, 52.2%) and (50.3, 51.2%).

Integration of multicompartment metabolic data using hierarchical-principal component analysis

A principal component analysis (PCA) model was initially constructed sepapan class="Species">rately for the metabolic data from each individual biological matrix (plasma, urine, liver, fecal extracts and pan class="Chemical">bile acid composition in ileal flush; Figures 3 and 4). Three principal components were calculated for each cross-validated PCA model, except for the plasma where four principal components were retained to maximize the explained variance R2X and the cross-validation parameter Q2 following the standard sevenfold cross-validation method (Cloarec ). These PCA models descriptors (R2X/Q2) were 0.70/0.36 (plasma), 0.41/0.11 (urine), 0.80/0.71 (liver), 0.94/0.67 (bile acid) and 0.61/0.42 (feces). The score vectors tb from each model were then assigned as new X-variables (Figure 3). Thus, the top level X-matrix contained 16 descriptors, denoted Pi (plasma PCs 1–4), Li (liver PCs 1–3), Ui (urine PCs 1–3), Bi (bile acid PCs 1–3) and Fi (fecal PCs 1–3), which comprise only the systematic variation from each of the blocks/compartments. The two first principal components (p1 and p2) calculated for the hierarchical-principal component analysis (H-PCA; Westerhuis ) model (R2X=0.60) accounted for 37 and 23% of the total variance in the combined multi compartment data respectively. The cross validation for the H-PCA model failed due to the high degree of orthogonality within the X-matrix, that is, within each of the blocks all variables are orthogonal to each other, while each of the biological matrices could be cross-validated at the individual level.

Figure 3

Schematic overview of H-PCA modelling: (Gunnarsson ). In the sublevel, each block of data XB is modelled locally by a PCA model. Each block is summarized by one or more loading vectors pb and score vectors tp (‘super variables'), which can be combined to form a new data matrix than can then be modelled using PCA, which generates the ‘super scores' tT and the ‘super loadings' pT. All conventional PCA statistics and diagnostics are retained.

Figure 4

H-PCA scores (A) and loadings (B) plots for the two first components derived from scores of separate PCA constructed separately for the metabolic data from each individual biological matrix from HBF mice (▪), HBF-L. paracasei mice (⧫) and HBF-L. rhamnosus mice (○). These PCA models explained 94% (bile acid, C), 70% (plasma, D), 80% (liver, E), 41% (urine, F) and 61% (feces, data not shown) of the total variation in the data, respectively. The systematic variation from each of the block/compartment is summarized by its score vectors denoted Pi (plasma PCs 1–4), Li (liver PCs 1–3), Ui (urine PCs 1–3), Bi (bile acid PCs 1–3) and Fi (fecal PCs 1–3), which can be combined to form a new data matrix than can then be modelled using PCA. The individual PCA loadings were color coded according to their contribution to the H-PCA model in red (principal component 1) and in blue (principal component 2). The model has been calculated from Pareto scaled data using two cross-validated PCs, R2X=60%. Ala, alanine; see Figure 1.

The H-PCA scores plot illustpan class="Species">rated a degree of clustering with respect to the groups of HBF pan class="Species">mice (Figure 4A). The corresponding H-PCA loadings plot indicated the contribution of the 16 descriptors to the differences observed between the samples in the H-PCA scores plot. Probiotic-supplemented HBF animals are separated from the controls along the first principal component, and this arises from the main variations modelled at the base level PCA of the individual plasma (P1, P2), liver (L2, L3), ileal flush (B2) and urine (U1, U2) data sets (Figure 4B). HBF mice fed with L. paracasei were separated from those fed with L. rhamnosus along the second principal component, which was mainly due to the variations modelled at the base level PCA of plasma (P2, P3), liver (L1) and urine (U1, U3). Interestingly, the metabolic variations in the fecal samples have no weight in discriminating the bacterial supplementation from the corresponding controls in the global model.

To uncover variables contributing to the H-PCA super scores, the loadings at the base level PCA model were interrogated (Figure 4C–F). HBF pan class="Species">mice supplemented with probiotics showed higher concentrations of glucose, choline, GPC, glutamine, glutamate and lysine in the plasma profiles associated with elevated concentrations of glucose in the liver and higher levels of TUDCA and TCDCA in the ileal flushes. Controls showed higher levels of lipoproteins in the plasma, elevated concentrations of lipids, glycogen, glutamine, glutamate, alanine, TMAO and lactate in the liver, associated with higher levels of TCA and TβMCA in the ileal content. Unsupplemented HBF mice also showed elevated urinary excretions of creatine, citrate, citrulline, lysine, UGLp, NAG and α-keto-isocaproate compared to animals fed with probiotics. Moreover, H-PCA also revealed that HBF mice treated with L. rhamnosus had higher levels of hepatic lipids and plasma lipoproteins but lower concentrations of lactate and amino acids in plasma and lower urinary excretion of PAG, IAG, tryptamine and taurine than HBF mice treated with L. paracasei.

Integration of correlations between bile acids and fecal flora

We further investigated the connections between fecal flora and intestinal pan class="Chemical">bile acids using a correlation analysis based bipan>rtite graphical modelling approach (see Materials and methods; Figure 5) used previously to investigate the effects of pan class="Species">gut microbiome humanization in germfree mice (Martin ). In the current study, we are working with a superior model where all the major bacteria strains are identified, which was not possible when considering conventional microflora. Positive and negative correlations show the multicolinearity between bile acids and gut bacteria, whose concentrations are interdependent such as in the case of substrate–product biochemical reactions. Additional pixel maps of the correlation matrices are given to help interpretation in Supplementary Figure 1.

Figure 5

Integration of bile acid and fecal flora correlations. The bipartite graphs were derived from correlations between fecal flora and bile acids in each group: HBF mice (A), HBF mice supplemented with L. paracasei (B) or L. rhamnosus (C). The cut-off value of 0.5 was applied to the absolute value of the coefficient ∣r∣ for displaying the correlations between fecal flora and bile acids. Bile acids and fecal bacteria correspond to blue ellipse nodes and green rectangle nodes, respectively. Edges are coded according to correlation value: positive and negative correlations are respectively displayed in blue and in red. aMA, α-muricholic acid; Ba, Bacteroides; Bb, B. breve; Bl, B. longum; bMA, β-muricholic acid; CA, cholic acid; CDCA, chenodeoxycholic acid; Cp, C. perfringens; Ec, E. coli; GCA, glycocholic acid; La, Lactobacillus probiotics; Sa, S. aureus; Se, S. epidermidis; TCA, taurocholic acid; TCDCA, taurochenodeoxycholic acid; TMCA, tauro-β-muricholic acid; TUDCA, tauroursocholic acid; UDCA, ursocholic acid. * and *** indicate a statistically significant correlation at 95 and 99.9% confidence levels, respectively.

Control HBF pan class="Species">mice and HBF pan class="Species">mice supplemented with probiotics show remarkably different bile acid/fecal flora correlation networks (Figure 5A–C), indicating that small modulations in the species composition of the microbiome can result in major functional ecological consequences. Network statistics reveal that microbiome/metabolome bipartite graphs from HBF mice supplemented with Lactobacilli show a totally different nodal structure (given for the cut-off value of 0.5). In particular, we observed that in the network obtained from control HBF mice, the most connected bacteria are Bifidobacteria, Staphylococci, and Clostridia, for which the variations are intrinsically correlated with the balance in tauro-conjugated bile acids (TCDCA, TβMCA, TCA, TUDCA) and unconjugated bile acids (CA, αMCA, βMCA). In particular, the potentially harmful opportunist C. perfringens shows functional correlation of opposite sign for TβMCA, TCDCA, αMCA and βMCA when compared to Bifidobacteria and S. aureus.

Network analysis for HBF pan class="Species">mice supplemented with pan class="Species">L. paracasei reveals that Lactobacilli supplementation resulted in decreasing the functional links between Bifidobacteria and bile acids, while new significant correlations were observed between bile acids, and Bacteroides, S. aureus, S. epidermidis and L. paracasei. Moreover, E. coli has several connections with UDCA, αMCA and TCDCA. Interestingly, Bacteroides shows functional correlations of opposite signs for TβMCA, TUDCA and TCA when compared to Lactobacilli and Staphylococci.

When HBF pan class="Species">mice received pan class="Species">L. rhamnosus probiotic, the microbiome/metabolome network shows a significant lower level of complexity (given for the cut-off value of 0.5). The balance within B. breve, S. aureus and S. epidermidis appears highly correlated to the composition in tauro-conjugated bile acids (TβMCA, TCA) and unconjugated bile acids (βMCA, CA).

Discussion

The cometabolic processes regulating pan class="Species">mammalian systems and their coexisting gut microbiota are an essential evolutionary driver towards providing more refined control mechanisms on the host physiology (Pereira and Gibson, 2002; Pereira ; Backhed ; Holmes and Nicholson, 2005; Nicholson ; Dumas ; Martin , 2007a). In the present study, we demonstpan class="Species">rate a significant association between the probiotic modulation of the gut microbiome and the metabolism of SCFAs, amino acids and methylamines, bile acids and plasma lipoproteins, and also an association with stimulated glycolysis, showing the remarkable diversity of symbiotic cometabolic connections.

Gut microbiotal links to host energy metabolism

We have recently described that HBF pan class="Species">mice supplemented with pan class="Species">L. paracasei were characterized by a high gut content of tauro-conjugated bile acids due to the inability of gut flora to deconjugate the bile acids, which resulted in high intestinal absorption of dietary lipids, accumulation of lipids in the liver and reduction of lipoprotein levels in plasma when compared to conventional animals (Martin ). Here, we show that probiotics supplementation of HBF mice resulted in a decrease in the plasma concentrations of VLDL and low-density lipoproteins (LDL) and increased triglyceride and GPC concentrations in plasma when compared to HBF controls (Figure 1B). In addition, the Lactobacillus supplementation resulted in decreased fecal excretion of bile acids (Figure 1C), that may be caused by accumulation of bile acids in Lactobacillus probiotics (Kurdi ). These observation indicated probiotic-induced changes in the enterohepatic recirculation of bile acids, which were shown to lower cholesterol and systemic levels of blood lipids (Pereira and Gibson, 2002). Moreover, probiotic-specific modulation of the ileal concentrations of UDCA and CDCA (Figure 2; Table IV) may also contribute to modulation of the synthesis and secretion of VLDL into the blood (Lin ; Watanabe ).

Moreover, the main source of dietary pan class="Chemical">lipids in animal chow is pan class="Species">soybean oil, which is composed at 65% of long-chain polyunsaturated fatty acids. It is well known that Lactobacillus hydrolyzes soy oil to conjugated linoleic acid efficiently (Xu ), which results in a reduction of plasma lipoprotein concentrations and hepatic cholesterol (Fukushima , 1997; Al-Othman, 2000) and in the inhibition of S. aureus growth (Das, 2002), as observed in the current study. Furthermore, probiotics supplementation was associated with significant reduction of acetate in the cecal content (Table II) and in a reduced hepatic acetate:propionate ratio, for which a serum lipids lowering effect has previously been described (Wong ). Our results illustrate the fine relationship between a specific gut microbial population modulation and the host's lipid metabolism and that a probiotic intervention can provide refined control mechanisms on the host's physiology.

Furthermore, the molecular foundations of beneficial symbiotic host–bacteria relationships lie in the critical involvement of the microbiome in calorie recovery through further processing of dietary nutrients and indigestible fibers. Levels of pan class="Chemical">leucine and pan class="Chemical">isoleucine were reduced together with their keto-acid derivative (α-keto-isocaproate) in L. paracasei-supplemented mice. These observations suggest higher catabolism of branched-chain amino acids to produce acetyl-CoA and glucose via gluconeogenesis. Decreased levels of citrate in urine and plasma, but increased liver succinate levels may also indicate the shunt of the tricarboxylic acid cycle towards production of phosphoenolpyruvate for gluconeogenesis in L. paracasei-supplemented mice. Moreover, reduction of liver glycogen observed with L. paracasei supplementation is consistent with our other observations of generalized host mobilization of other metabolic fuels.

Probiotics induce specific microbiome–host transgenomic metabolic interactions

We investigated the relationship between probiotic-induced changes in gut microbes and pan class="Chemical">bile acid cometabolism using bipan>rtite graphs to display correlation patterns between fecal flora and pan class="Chemical">bile acids (Figure 5). Correlation analysis derived from bile acid and fecal flora profiles offers a unique approach to capture subtle variations in bile acid composition that may be directly related to changes in gut microbial levels, and that may be induced by accumulation of bile acids in Lactobacillus probiotics for instance. Control HBF mice and HBF mice supplemented with probiotics show remarkably different bile acid/fecal flora correlation networks. The different bacterial strains of Bacteroides, Clostridia, Streptococci and Lactobacilli share similar abilities to deconjugate the hepatic tauro-conjugated bile acids (Midtvedt and Norman, 1967; Floch, 2002). In that regard, the overwhelming contrast between the balance of these bacteria on one hand, and conjugated bile acids (TCDCA, TβMCA, TCA, TUDCA) and unconjugated bile acids (CA, αMCA, βMCA) on the other hand, highlights the metabolic flexibility of the gut microbiota in response to probiotics supplementation. These different correlative patterns further characterize the microbial–mammalian transgenomic metabolic interactions, whereby probiotics-induced modulation of the gut microbial functional ecosystem results in different bile acid composition (Figure 2) and enterohepatic recirculation.

The relationship between specific gut microbial strains and pan class="Chemical">bile acid cometabolism is well illustpan class="Species">rated with the contrast between Lactobacillus, which shows resistance to bile salt toxicity (Corcoran ), and C. perfringens, which is sensitive to the strong growth inhibitory effects of unconjugated bile acids and TCDCA (Kishinaka ; Floch, 2002). In the absence of Lactobacillus supplementation, C. perfringens has anticorrelated connections with TCDCA, αMCA and βMCA, which highlights the strong inhibitory effects of these bile acids on C. perfringens growth. Interestingly, the probiotic supplementation was either associated with maintenance or decrease of the Clostridial population (Table I), while no functional correlations between C. perfringens and bile acids were observed. In that regard, these observations may indicate different nutritional competition leading to modulation of C. perfringens population and maintainance of the intestinal ecology. Consequently, interbacterial cooperation to transform bile acids is an important factor that needs to be considered not only for the fine tuning of microbial balance but also to modulate dietary fat emulsification and absorption.

Gut-bacterial production of methylamines via choline metabolism

The elevation of pan class="Chemical">methylamines (pan class="Chemical">TMA, TMAO) in liver, choline and GPC in the plasma and decreased choline in feces from HBF mice supplemented with L. rhamnosus (Figure 1) are additional illustrations of the complex ‘microbial–mammalian metabolic axis,' as the microbiota are involved in the synthesis and metabolism of these methylamines (al-Waiz ). The first reaction of the methylamine pathway involves conversion of dietary choline into TMA by gut microbiota (Zeisel ), which is then detoxified to TMAO in the liver via the flavine monooxygenase system (Smith ; (Figure 6). L. rhamnosus supplementation contributes to higher absorption of free choline and may induce elevated production of methylamines by Bacteroides and C. perfringens (Allison and Macfarlane, 1989) through nutritional competition.

Figure 6

Gut-microbiota–mammalian cometabolism of methylamines and aromatic amino acids. DMG, dimethylglycine; IA, indoleacetate; MA, methylamine; PA, phenylacetate (see Figure 1).

Interestingly, pan class="Species">L. paracasei consumption may favor a different metabolic fate for pan class="Chemical">choline through different bacterial reprocessing. Decreased fecal choline was associated with reduced concentrations of TMA and DMA in the liver, an increase in plasma GPC, but with no changes in liver TMAO and plasma choline after L. paracasei supplementation. These animals also showed a greater reduction in plasma lipoproteins when compared to other groups. Here, these observations may result from elevated bacterial consumption for cholesterol assimilation (Rasic, 1992) and phospholipid metabolism (Jenkins and Courtney, 2003; Taranto ; Kankaanpaa ). Thus, the reduced availability of choline to other bacterial strains may have led to lower production of methylamines and absorption of free choline into host metabolism.

Probiotic modulation of amino-acid metabolism

Investigation of the urine metabolic profiles showed significant increases in the concentpan class="Species">rations of microbial cometabolites pan class="Chemical">PAG, IAG and tryptamine in probiotic supplemented mice (Goodwin ; Smith and Macfarlane, 1996). These metabolites are produced from amino acids, which after depolymerization of dietary proteins (casein) by pancreatic endopeptidases and bacterial proteases and peptidases, become available for fermentation by the gut microflora (Smith and Macfarlane, 1996), as outlined in Figure 6. Metabolism of the aromatic amino acids phenylalanine, tyrosine and tryptophan generates phenylacetate (PA), p-cresol, indoleacetate (IA) and tryptamine, respectively (Donaldson, 1962; Smith and Macfarlane, 1996). PA and IA may be detoxified in the gut mucosa and the liver by glycine conjugation forming PAG and IAG prior to excretion via the urine (Donaldson, 1962; Smith and Macfarlane, 1996). The production of PA and IA has been restricted to a certain taxonomic group of gut bacteria, including Bacteroides, Clostridia and E. coli, which count among the dominant species in HBF colonized mice (Smith and Macfarlane, 1996; Xu ). The increased urinary excretion of phenolic and indolic compounds reflects variations in gut microflora composition in relation to nutritional competition (Smith and Macfarlane, 1996; Nowak and Libudzisz, 2006). For instance, IA has been reported to inhibit the growth and survival of Lactobacilli, and specifically L. paracasei (Nowak and Libudzisz, 2006).

In addition, increased levels of pan class="Chemical">lysine in feces from pan class="Species">mice supplemented with probiotics is another example of bacterial contribution to the mammalian amino-acid homeostasis (Metges, 2000) and hepatic protein synthesis (Metges ; Figure 1). Higher production of bacterial isobutyrate and isovalerate (Table II) suggests increased bacterial fermentation of leucine and valine that can also influence host energy metabolism (Macfarlane ). Moreover, L. paracasei supplementation specifically induced higher urinary excretion of NAG and decreased N-acetylated metabolites in fecal extracts. A relationship between high casein diet and urinary excretion of NAG has been described previously (Hallson ), which also suggested elevated bacterial proteolysis.

Altogether, our data suggest that the probiotic-induced increased proteolytic activities may reflect the basal metabolism of these pan class="Species">Lactobacillus strains, in particular pan class="Species">L. paracasei for which the proteolytic activities on casein medium are known to be elevated (Sasaki ; Ikram and Mukhtar, 2006). Moreover, increased urinary NAG levels have been reported as a biomarker of increased tubular activity and tubular cell toxicity. Changes in bacterial fermentation of carbohydrates can lead to different ion absorption from the gut (Scholz-Ahrens ), which may contribute to altered kidney metabolism and tubular activity. Therefore, the gut microbial contribution to the modulation of the NAG biomarker is certainly of potential toxicological assessment significance.

H-PCA modelled multicompartmental matrices related to lipid metabolism

H-PCA has been explored for the first time as a top-down systems approach to model and integpan class="Species">rate metabolic profiles from diverse biological compan>rtments. H-PCA modelling summarized clearly the intercorrelated changes induced by probiotic-treatment in plasma, urine and liver matrices and composition in pan class="Chemical">bile acids (Figure 4). One benefit of the hierarchical approach lies in the much simplified simultaneous visualization of global system biochemical changes in multibiological matrices and improves interpretability. For instance, H-PCA resulted in a separation between the treated groups, which was not observed in separate PCA models. In addition, the H-PCA loadings plot gives the relative importance of the different blocks (biological matrices) in carrying diet-induced discriminant information. Moreover, the relationships between the descriptors in the simplified H-PCA space (H-PCA loadings plot) indicate correlations and anticorrelations between block variables, which may give insight into correlations between the metabolic variations in different biological matrices as exemplified here by metabolic changes in liver, plasma and urine and related to lipid metabolism. The H-PCA model also efficiently summarized the intercorrelated changes related to higher systemic glycolysis in plasma, urine and liver matrices, that is, reduced ketone body formation, anaerobic glycolysis, tricarboxylic cycle perturbation and amino-acid catabolism (Figure 4). Such observations might lead to better description of multiorgan metabolic perturbations (Figure 4B). This multicompartmental top-down approach offers a way forward to study the systemic biochemical profiles and regulation of function in whole organisms by analyzing simultaneously several metabolite pools from different biofluids and tissues. This approach also provides a new strategy for the quantitative and qualitative evaluation of different probiotic (or indeed any functional food or nutriceutical) interventions in relation to host biochemistry.

Conclusion

Significant associations between host metabolic phenotypes and a nutritionally modified gut-microbiota strongly supports the idea that changes across a whole range of metabolic pathways are the product of extended genome perturbations that can be oriented using probiotic supplementation, and which may play a role in host metabolic health. Bipartite network analysis highlights the metabolic flexibility of the gut microbiota, whereby bacterial strains communicate with each other to metabolize differently pan class="Chemical">bile acids in a gut microbial ecosystem modulated with probiotics. In this case, probiotic consumption exerted a modification over the microbiome resulting in different hepan>tic influx and efflux of pan class="Chemical">fatty acids in the liver, as observed with increased enterohepatic recycling of bile acids and dietary fats, lowered plasma lipoprotein levels and stimulated glycolysis. Probiotics also induced a different microbial proteolytic activity as well as modulation of bacterial metabolism of amino acids, methylamines and SCFAs. We showed the novel application of H-PCA as a means to study perturbation of metabolic profiles triggered by symbiotic microbiota at a ‘global system' level by analyzing several metabolite pools simultaneously from different biofluids and tissues. These integrated system investigations demonstrate the potential of metabolic profiling as a top-down systems biology driver for investigating the mechanistic basis of probiotic action and the therapeutic surveillance of the gut microbial activity related to dietary supplementation of probiotics and their health consequences.

Materials and methods

Animal handling procedure

All animal studies were carried out under appropriate national guidelines at the Nestlé Research Center (Lausanne, Switzerland). The HBF is constituted of a total of seven bacterial strains, isolated from stool of a 20-day-old female baby that was naturally delivered and breast-fed, namely pan class="Species">E. coli, pan class="Species">B. breve and B. longum, S. epidermidis and S. aureus, C. perfringens and Bacteroides distasonis. Bacterial cell mixtures contain approximately 1010 cells/ml for each strain and were kept in frozen aliquots until use. Lactobacillus paracasei NCC2461 and L. rhamnosus NCC4007 probiotics were obtained from the Nestlé Culture Collection (Lausanne, Switzerland).

A total of 28 female germ-free pan class="Species">mice (C3H strain), aged 6 weeks, received a single dose of HBF bacteria mixture and will be called HBF pan class="Species">mice in the current manuscript. The experimental design is detailed in Supplementary Figure 2. The animals were fed with a standard pathogen-free rodent diet constituted of 50% cornstarch, 20% casein, 10% sucrose, 7% soybean oil, 5% cellulose, 0.25% choline bitartrate, 0.3% cystine and vitamin and mineral mixtures (Reeves ) for 2 weeks. A control group of HBF mice (n=10) received a saline drink ad libitum containing Man, Rogosa and Sharpe (MRS) culture medium and was fed with a basal mix diet containing in composition 2.5% of a glucose–lactose mixture (1.25% each) for 2 additional weeks. Two groups of HBF mice were given a daily probiotic supplement, either L. paracasei (group A, n=9) or L. rhamnosus (group D, n=9), containing around 108 probiotic bacteria in MRS per day mixed with the saline solution ad libitum and were also fed with the basal mix diet.

Fecal pellets were first collected in the morning for microbiological analysis and frozen at −80°C for NMR spectroscopic analysis. Urine was then collected from animals prior to euthanasia and frozen at −80°C. Urine samples were not obtained for every animal, as some pan class="Species">mice had an empty bladder at the time of killing (group A, n=8; group B, n=6; group C, n=8). Blood (400 μl) was collected into Li-pan class="Chemical">heparin tubes and plasma was obtained after centrifugation and frozen at −80°C. The ileum and liver were dissected and snap-frozen. Ileal flush samples were obtained by rinsing the ileal lumen using a 1 ml sterile syringe containing a phosphate buffer solution (0.2 M Na2HPO4/0.04 M NaH2PO4, pH 7.4) for UPLC-MS analysis. Cecal content was collected upon animal autopsy, snap-frozen immediately and maintained at −80°C prior to analysis.

Microbial profiling of fecal contents

Immediately after collection, fecal pellets were homogenized in 0.5 ml Ringer solution (Oxoid, UK) supplemented with 0.05% (w/v) pan class="Chemical">L-Cystein (pan class="Chemical">HCl). The enumeration of specific microorganisms was performed after plating and incubation of different dilutions of the bacterial solution on selective and semiselective media, for example, Bifidobacteria on Eugom Tomato medium, Lactobacillus on MRS+antibiotic (phosphomycine, sulfamethoxazole, trimethoprim) medium, C. perfringens on NN-agar medium, Enterobacteriaceae on Drigalski medium and Bacteroides on Shaedler Neo Vanco medium. Enterobacteriaceae cultures were incubated at 37°C under aerobic conditions for 24 h, and other cultures were incubated under anaerobic conditions over a 48 h period.

Gas-chromatography on cecal content

An aliquot of cecal content was extracted with 4 ml buffer (0.1% (w/v) pan class="Chemical">HgCl2 and 1% (v/v) pan class="Chemical">H3PO4 supplemented with 0.045 mg/ml 2,2-dimethylbutyric acid (as internal standard) per gram fresh weight. The resulting slurry was centrifuged for 30 min at 5000 g at 4°C. Fecal SCFAs were analyzed using a gas-chromatograph (HP 6890) equipped with a flame ionization detector and a DB-FFAP column (J&W Scientific, MSP Friedli & Co, Switzerland) of 30 m length, 530 μm diameter and 1 μm film thickness. The system was run with helium gas at an inlet constant pressure of 10 psi at 180°C. Each sample run was preceded with a cleaning injection of 1.2% formic acid. Samples were run at an initial temperature of 80°C for 1.2 min followed by heating to 145°C in 6.5 min, heating to 200°C in 0.55 min and an additional 0.5 min at 200°C. SCFAs were identified using external standards (acetate, propionate, iso-butyrate, n-butyrate, iso-valerate, n-valerate) and the concentration was calculated using the internal standard.

1H NMR spectroscopic analysis

A volume of 100 μl of blood plasma was added to 450 μl of saline solution containing 10% pan class="Chemical">D2O, which was used as a spectrometer field frequency lock, into 5 mm NMR tubes. Urine samples were prepan>red by mixing 20 μl of samples with 30 μl of a phosphate buffer solution containing 90% D2O and 0.25 mM 3-trimethylsilyl-1-[2,2,3,3-2H4] propionate (TSP), which was used as chemical shift reference (δ 0.0), into 1.7 mm NMR tubes. Fecal pellets were homogenized in 650 μl of a phosphate buffer solution containing 90% D2O and 0.25 mM TSP. The fecal samples were sonicated for 30 min at 25°C and then centrifuged at 13 000 r.p.m. for 20 min to remove particulates. The supernatants were removed and centrifuged at 13 000 r.p.m. for 10 min. A 580 μl aliquot of the fecal supernatant was then pipetted into a 5 mm NMR tube for spectroscopic analysis. Portions of intact liver samples (∼15 mg) were bathed in ice-cold 0.9% saline D2O solution and packed into a zirconium oxide 4 mm outer diameter rotor.

1H NMR spectra were acquired for each sample using a Bruker DRX 600 NMR spectrometer (Rheinstetten, Gerpan class="Species">many) opepan class="Species">rating at 600.11 MHz for 1H. The instrument was equipped with a Bruker 5 mm TXI triple resonance probe maintained at 298 K for liquid samples and a standard Bruker high-resolution MAS probe under magic-angle-spinning conditions at a spin rate of 5000 Hz for intact tissues (Waters ). Tissue samples were regulated at 283 K to minimize biochemical degradation.

One-dimensional (1D) 1H NMR spectra were obtained from each sample using a standard solvent suppression pulse sequence (RD-90°-t1-90°-tm-90°-acquire FID) with tm fixed at 100 ms and t1 at 3 μs (Wang ). Additional spin echo Carr-Purcell-Meiboom-Gill (CPMG) spectra were acquired for plasma and liver samples using the pulse sequence (RD-90°-(t-180°-t)n-acquire FID), with a spin–spin relaxation delay, 2nτ, of 160 ms for plasma and 200 ms for tissue (Meiboom and Gill, 1958). The 90° pulse length was 9.0–12 μs. A total of 128 transients were collected into 32K data points with a recycle delay (RD) of 2 s. The assignment of the 1H NMR spectral peaks to specific metabolites was achieved based on the litepan class="Species">rature (Nicholson ; Fan, 1996), and confirmed by 2D COrrelation SpectroscopY (COSY) (Hurd, 1990) and TOtal Correlation SpectroscopY (TOCSY) (Bax and Davis, 1985). 2D NMR spectra were acquired on selected samples. Further assignment of the metabolites was also accomplished with the use of Statistical TOtal Correlation SpectroscopY (STOCSY) on 1D spectra (Cloarec ).

UPLC-MS methods

The Ultra Perforpan class="Species">mance™ liquid chromatography of ileal flushes was performed on a ACQUITY UPLC system (Waters, Milford, MA, USA) equipped with a Tof™ LCT-Premier (Waters MS Technologies, Manchester, UK) for mass detection using the method we described previously (Martin ). The same conditions were applied for analysis of the following bile acid standards: cholic (CA), taurocholic (TCA), glycocholic (GCA), deoxycholic (DCA), taurodeoxycholic (TDCA), glycodeoxycholic (GDCA), chenodeoxycholic (CDCA), taurochenodeoxycholic (TCDCA), glycochenodeoxycholic (GCDCA), lithocholic (LCA), taurolithocholic (TLCA), glycolithocholic (GLCA), ursodeoxycholic (UDCA), tauroursodeoxycholic (TUDCA), glycoursodeoxycholic (GUDCA), hyocholic (HCA), α-muricholic (αMCA), β-muricholic (βMCA), ω-muricholic (ωMCA), tauro-β-muricholic (TβMCA) and tauro-α-muricholic (TαMCA) acids. We have previously given the molecular structure, the retention time and the mass to charge ratio (and not the molecular weight as published previously) of the observed ions. Unconjugated bile acids formed a formate adduct.

Data analysis

Microbial counts and pan class="Chemical">SCFAs composition in the cecum were analyzed using a two-tailed pan class="Species">Mann–Whitney test.

1H NMR spectra were corrected pan class="Species">manually for phase and baseline distortion and referenced to the chemical shift of the CH3 resonance of pan class="Chemical">alanine at δ 1.466 for plasma and liver samples, to the TSP resonance of alanine at δ 0.0 for urine and fecal samples using XwinNMR 3.5 (Bruker Biospin, Rheinstetten, Germany). The spectra were converted into 22K data points over the range of δ 0.2–10.0 using an in-house developed MATLAB routine. The regions containing the water resonance (δ 4.5–5.19), and for urine urea resonance (δ 4.5–6.2), were removed. Chemical shift intensities were normalized to the sum of all intensities within the specified range before chemometric analysis.

UPLC-MS data were processed using the Micromass MarkerLynx™ applications mapan class="Gene">nager Version 4.0 (Waters Corp, Milford, USA). The peaks of pan class="Chemical">bile acids were identified by comparing the m/z ratio and retention time to the set of standard bile acids measured under the same conditions. Data were noise-reduced in both of the UPLC and MS domains using MarkerLynx standard routines. Integration of the UPLC-MS bile acid peaks was performed using ApexTrack2™. Each peak integral was expressed as a ratio to the sum of integrals of the 21 measured bile acids.

The multivariate pattern recognition techniques used in this study were based on PCA (Wold ), H-PCA (Westerhuis ) and the orthogonal-projection to latent structure (O-PLS) (Trygg and Wold, 2003). PCA was carried out using the SIMCA-P 11 software (Umetrics, Umeå, Sweden) in order to detect the presence of inherent similarities between metabolic profiles. Both NMR and LC-MS variables were subjected to pan class="Chemical">Pareto scaling, by dividing each variable by the square root of its standard deviation. Data were visualized by means of principan>l component scores and loadings plots. Each point in the scores plot represents an individual biochemical profile of a sample, whereas on the loadings plot each coordinate represents a single NMR spectral region or LC-MS retention time and mass to charge pan class="Species">ratio (m/z). Because the scores and loadings plots are complementary, biochemical components responsible for the differences between samples detected in the scores plot can be extracted from the corresponding loadings plot.

O-PLS-DA was also carried out using the method developed by Trygg et al (Trygg and Wold, 2003) and implemented for NMR spectral data by Cloarec to exclusively focus on the effects of probiotic supplementation. All O-PLS-DA models were constructed using one predictive and two orthogonal components using data scaled to unit variance (i.e. by dividing each variable by its standard deviation). Here, the test for the significance of the Pearson product–moment correlation coefficient was used to calculate the cut-off value of the correlation coefficients at the level of P<0.05. To test the validity of the model against over-fitting, the cross-validation parameter Q2 was computed and the standard seven-fold cross validation method was used (Cloarec ). The interpretation of the model was achieved from correlation coefficient plots that incorpopan class="Species">rated a back-transformation method such that the coefficients resembled the original structure of the NMR spectral data (Cloarec ).

H-PCA methods have been proposed in the recent litepan class="Species">rature to improve the interpretability of multiple models (Westerhuis ) generated from several blocks of descriptor variables measured on the same objects (Bergman ; Lundstedt ; Janne ). The data are divided into well-defined logical blocks X1, …, XB (where the index B refers to each biosample type), according to an event, instrumentation or biological compartment. A consensus direction among all the blocks is sought. Here, logical blocks were defined according to the biological nature of the samples obtained from the same animals (Figure 3). Hence, five independent blocks of variables were comprised of NMR data from plasma, liver, urine and fecal extracts as well as bile acids composition obtained by UPLC-MS. PCA was performed on each data matrix individually, and block loadings pb and scores tb were generated for each block (Figure 3). All block scores tb were then combined into a super block T. The PCA (H-PCA) was then applied to the super block T to generate the super scores tT and the super loadings pT. The super scores tT gives the relationships between the observations and all the sources of variations. The super loadings pT gives the relative importance of the different blocks XB for each principal component as well as information on the relationship/covariance between the different variables in the different compartments.

Bipartite graph representation of bile acid and fecal flora profiles

The bipartite graph (Rgraphsviz) package from R was used to display the correlation matrix derived from pan class="Chemical">bile acid and fecal flora profiles to assess the probiotic-induced changes to the microbial-pan class="Species">mammalian transgenomic interactions (Figure 5). Pearson's correlation coefficients were computed between bile acid variables and fecal flora variables from the same mice and a cut-off value of 0.5 was applied to the absolute value of the coefficient ∣r∣ so that the bipartite graph only represents the correlations between the two types of nodes (fecal flora and bile acids) above the cut-off (Martin ). The sign of the initial correlation was then color coded (red: negative, blue: positive), and the correlation value displayed on the bipartite graph. In that context, presence of edges between two specific nodes (one of each type) reveals a functional correlation (above the cut-off) between these entities.

Supplementary Table 1

Supplementary Figure 1

Supplementary Figure 2