Maternal microbiota-derived metabolic profile in fetal murine intestine, brain and placenta
BMC Microbiology volume 22, Article number: 46 (2022)
The maternal microbiota affects the development of the offspring by microbial metabolites translocating to the fetus. To reveal the spectrum of these molecular mediators of the earliest host-microbe interactions, we compared placenta, fetal intestine and brain from germ-free (GF) and specific pathogen free (SPF) mouse dams by non-targeted metabolic profiling.
One hundred one annotated metabolites and altogether 3680 molecular features were present in significantly different amounts in the placenta and/or fetal organs of GF and SPF mice. More than half of these were more abundant in the SPF organs, suggesting their microbial origin or a metabolic response of the host to the presence of microbes. The clearest separation was observed in the placenta, but most of the molecular features showed significantly different levels also in the fetal intestine and/or brain. Metabolites that were detected in lower amounts in the GF fetal organs included 5-aminovaleric acid betaine, trimethylamine N-oxide, catechol-O-sulphate, hippuric and pipecolic acid. Derivatives of the amino acid tryptophan, such as kynurenine, 3-indolepropionic acid and hydroxyindoleacetic acid, were also less abundant in the absence of microbiota. Ninety-nine molecular features were detected only in the SPF mice. We also observed several molecular features which were more abundant in the GF mice, possibly representing precursors of microbial metabolites or indicators of a metabolic response to the absence of microbiota.
The maternal microbiota has a profound impact on the fetal metabolome. Our observations suggest the existence of a multitude of yet unidentified microbially modified metabolites which pass through the placenta into the fetus and potentially influence fetal development.
The intestinal microbiota has a great impact on the life and wellbeing of the host. Microbes residing in the gut participate in digestion and metabolic modification of nutrients, producing substances that are absorbed by the host . The cell numbers of gut microbiota are estimated to at least equal and its gene pool exceed that of its host [2, 3]. The potential of the microbiota to influence host metabolism is illustrated by a comparison of the serum metabolome of conventionally colonized and germ-free (GF) mice: 3.5% of the > 4000 molecular features detected were unique for conventional mice and 10% of the shared molecular features had significantly different levels between the groups . All the organ systems of the host are affected to a varying degree .
While a majority of the compounds originating from microbial metabolism detected in mammalian tissues still remain uncharacterized, some of these substances and their effects on the host are well documented. These include short chain fatty acids (SCFAs) which the host utilizes as an essential part of its metabolism . The SCFAs produced from complex carbohydrates by microbes residing in the alimentary tract are an important source of energy for the host . The SCFAs have also been shown to contribute to the maintenance of the gut epithelium and regulation of the immune responses by facilitating regulatory T cell generation in the colonic mucosa [8,9,10,11,12]. Gut-residing microbes are also known to modify endogenous primary bile acids creating molecular species such as deoxycholate that stimulates colonic enteroendocrine cells and thus affects the regulation of the intestinal function of the host . Other microbial metabolites or their host-produced derivatives, such as trimethylamine N-oxide (TMAO) and 5-aminovaleric acid betaine (5-AVAB), are known to modify specific host reactions of lipid metabolism [14,15,16]. Microbiota also affects the levels of intestinal and absorbed nutrients, particularly amino acids [4, 17, 18]. Microbial metabolites of amino acids, such as the indole derivatives of tryptophan, are signaling molecules with both local and systemic effects .
The metabolic coexistence between the animal and the bacteria begins already before birth. While it is still unclear whether small numbers of live microbes exist in the healthy fetus, hundreds of microbial metabolites originating from the dam pass through the placenta [20, 21]. Very little is known of their properties and physiological effects across developing organs of the fetus . Microbe-derived aryl hydrocarbon receptor (AhR) ligands and microbially regulated retinoids are essential for fetal development of the immune system [20, 23, 24]. Maternal SCFAs are also readily transmitted to the fetus, programming the fetal metabolic and neural systems . Other maternally derived microbial metabolites have been primarily studied in the context of toxicology . These observations suggest that whole bacteria are not necessarily required to inflict inflammatory immune responses by the host cells .
To examine the extent of the cross-placental transfer of microbial metabolites during pregnancy, we analyzed placenta, fetal intestine and brain samples from germ-free (GF) and specific pathogen free (SPF) murine dams using a broad non-targeted metabolomics approach (Fig. 1). Ultra-high performance liquid chromatography (UHPLC) coupled with quadrupole time-of-flight (QTOF) mass spectrometry allowed the detection of thousands of differentially abundant molecular features in the organ samples.
Differences in all observed molecular features
The non-targeted metabolomics data consisted of a total of 12,166 molecular features from four analytical modes after data cleanup. The metabolic profiles were clearly different in all studied organs (Fig. 2 and Table S1). The GF and SPF mice clustered separately in t-distributed stochastic neighbor embedding (TSNE) analysis, especially when each organ was analyzed individually (Fig. 2). The clearest separation between GF and SPF samples was observed in the placenta. The gender of the fetus did not influence the separation (not shown).
The clustering by organ is also evident in the heatmap of all observed molecular features (Fig. 3). At this level, the difference in signal abundance related to germ-free status can be observed from a few relatively small clusters of molecular features in each analyzed organ.
The concentrations of 3680 molecular features differed between germ-free (GF) and specific pathogen free (SPF) mice in at least one of the organs investigated (FDR-adjusted p < 0.05 and Cohen’s d > 0.8). There were 2200 features which were more abundant in SPF mice in at least one organ (Fig. 4a). These were most numerous in the fetal intestine. One hundred sixty-eight features were more abundant in SPF mice in all three organs investigated. Similarly, 1533 features were more abundant in GF mice in at least one organ, most commonly in placenta (Fig. 4b). Eighty-eight features were more abundant in GF mice in all organs.
A total of 99 features were only observed in SPF mice (Fig. 5). These were most commonly detected in all three SPF organs (n = 37), or in both placenta and fetal intestine (n = 36). None were detected only in both fetal organs, or only in the fetal brain. In contrast, only 6 features were observed exclusively in GF mice (not shown).
Among the differentially abundant molecular features, 101 metabolites were annotated. Out of these, 54 were identified (MSI level 1), 35 putatively annotated (level 2), 10 putatively characterized for compound class (level 3), and 2 unknowns (level 4) given a molecular formula (Table S1). A heatmap of significantly differential annotated metabolites is shown in Fig. 6. Sixty-one of these metabolites were more abundant in SPF mice in at least one organ, most of these in all three organs or in intestine and/or placenta (FDR-adjusted p < 0.05, Cohen’s d > 0.8; Figs. 6 and 8, Table S1).
Volcano plots of FDR-adjusted p values and Cohen’s d values for each organ are shown in Fig. 7. Twenty-five metabolites were more abundant in all three SPF organs. Of these, 5-AVAB [also known as δ-valerobetaine (δVB), N,N,N-trimethyl-5-aminovalerate (TMAV) and N,N,N-trimethyl-5-aminovaleric acid (TMAVA)] was the most affected metabolite in all three organs. Other metabolites significantly affected included TMAO, alanine / β-alanine betaine, solanidine, catechol-O-sulphate, hippuric and pipecolic acid, amino acids and their derivatives (such as kynurenine, 3-indolepropionic acid and aminoisobutyric acid) and some small peptides. Five of the annotated compounds were observed exclusively in SPF mice: benzamide, 4-hydroxybenzenesulfonic acid, two unidentified alkaloids and a triterpenoid. Forty annotated metabolites were more abundant in GF mice in at least one organ, primarily in placenta and/or brain (Fig. 8). These included several acylcarnitines, phosphatidylcholine, amino acids (such as L-threonine, L-arginine, ergothioneine and L-ornithine) and several small peptides.
In order to predict metabolic pathways affected by the lack of microbiota in germ-free mice, we performed the MS Peaks to Pathways analysis in the MetaboAnalyst pipeline . This module utilizes the mummichog algorithm and gene set enrichment analysis (GSEA), which fit the mass spectrometry peak data into known metabolic pathways without pre-existing compound annotations .
In terms of murine KEGG (Kyoto Encyclopedia of Genes and Genomes) and BioCyc Genome Database Collection pathways, the metabolism of several essential amino acids, as well as tRNA charging / aminoacyl-tRNA biosynthesis were significantly different in both fetal organs (Supplementary Table 2) [28, 29]. In the fetal intestine, phosphonate & phosphinate metabolism, glycolysis, caffeine metabolism, degradation of putrescine and nicotine, methionine salvage, and lipoate biosynthesis and incorporation were affected. In placenta, the affected pathways included lactose degradation, methionine salvage, nicotine degradation and folate biosynthesis.
This is the first study probing the effects of a complete maternal microbiota on the metabolite composition in mammalian placental and fetal tissues. We compared fetuses of germ-free (GF) and specific pathogen free (SPF) murine dams utilizing non-targeted metabolomics. The lack of maternal microbiota affected the metabolite profile of the placenta, fetal intestine and brain. A total of 2200 detected molecular features were more abundant in SPF mice in at least one organ, while more than 1500 showed higher levels in the GF mice. Approximately 100 molecular features could be detected only in the SPF organ samples. The numbers of compounds depleted in GF mice were largest in the fetal intestine and/or placenta. These observations indicate that the maternal microbiota strongly affects the host metabolism in placenta and in the fetus not only by directly producing metabolites but also pervasively impacting host physiology.
One hundred one of the differentially abundant metabolites could be annotated based on current databases. The limited number of identified metabolites, which is typical for non-targeted metabolomics study, is explained by several factors. Our metabolite identification was focused to the molecular features that showed significant difference between the experimental groups, i.e., to those metabolites that are associated with mouse gut microbiota. There are very limited reference data available for mouse gut microbiota associated metabolites when compared to, e.g., with metabolites found in human plasma or serum, reducing the number of possible identifications. Furthermore, we report each metabolite only once as an identified metabolite in the final results although a single metabolite can produce several molecular features in the data (the metabolite is measured in several modes of analysis and produces adducts and/or fragments).
Several betaines, amino acids, and their derivatives, such as pipecolic acid and small peptides, certain alkaloids, catechol-O-sulphate and hippuric acid were more abundant in SPF mice. Several acylcarnitines, some amino acids, several small peptides and phosphatidylcholine were more abundant in GF mice.
Trimethylated compounds 5-AVAB, (β-)alanine betaine and TMAO
Considerably lower levels of the trimethylated compounds, 5-AVAB, (β-)alanine betaine and TMAO, were detected in all the studied organs of the GF mice. These compounds are zwitterions containing a positively charged trimethylammonium group and a negatively charged oxygen. They can act as osmoprotective substances and methyl donors in various cellular metabolic processes [30,31,32]. 5-AVAB is a bacterial metabolite [30, 33] with reported effects on fatty acid metabolism. It inhibits oxygen consumption due to β-oxidation of fatty acids in cultured mouse cardiomyocytes . Thus, it may protect the heart tissue in ischemic conditions. The cord plasma of pre-eclamptic infants has been shown to contain increased levels of 5-AVAB [34, 35]. 5-AVAB has also been detected in the fetal mouse brain where its quantity has been shown to be strongly dependent on the presence of maternal gut microbiota . Lower levels of (β-) alanine betaine have been detected in the intestine of GF mice than in conventional mice suggesting that it can also be produced by gut microbiota . Its role in mammalian physiology is poorly known.
TMAO is the mammalian end product of dietary phosphatidylcholine, choline, and carnitine, which are first metabolized by gut microbes into trimethylamine (TMA) and further into its N-oxide form in the liver [37, 38]. Whilst red meat is the main dietary source for the TMAO precursors choline and carnitine, TMAO as such is abundant in seafood . TMAO inhibits reverse cholesterol transport by affecting bile acid synthesis on multiple levels and increases deposition of cholesterol to arterial walls . Elevated levels are associated with increased risk of cardiovascular diseases, such as atherosclerosis and thrombosis [40, 41], but a causative role of TMAO has not been definitely established . A recent study indicates that TMAO and possibly also 5-AVAB have a role in development of the fetal brain . More than twofold reductions of TMAO and 5-AVAB were detected in the maternal blood and fetal brain of both GF and antibiotic-treated pregnant mice relative to SPF controls. Colonization of GF dams with ingenious spore-forming bacteria increased the levels of these compounds in the fetal brain. The maternal microbiota was shown to promote fetal thalamocortical axonogenesis and the effect could be reproduced by supplementing microbiota-depleted dams with TMAO .
Amino acids and related metabolites
Concentrations of several amino acids and their derivatives were significantly different in the organs of the test groups. The pathway analysis based on unannotated MS peak data further substantiated the broad impacts of maternal microbiota on amino acid metabolism and tRNA charging especially regarding essential amino acids. Earlier studies have shown that the GF status is associated with differences in the levels of several amino acids in plasma and intestine [4, 17, 18]. Notably, microbial metabolism can diversify the fates of the amino acid tryptophan. On the one hand, microbes can produce an array of microbial indole metabolites, such as indole-3-propionic acid (IPA), in the gut lumen. On the other, they can also act indirectly by modifying the host metabolism of tryptophan to its neuromodulatory metabolite serotonin by enterochromaffin cells and to kynurenine by immune and epithelial cells [13, 42,43,44]. Adult GF mice have been reported to exhibit lower levels of serotonin in plasma and brain but higher levels of tryptophan in plasma than in conventional mice [4, 42].
We observed lower levels of tryptophan and its metabolites serotonin, hydroxyindoleacetic acid (HIAA), kynurenine and 3-indolepropionic acid (3-IPP) in intestine or brain or both organs of the GF fetuses. The levels of serotonin and 3-IPP were also reduced in placenta. Kynurenine and HIAA can act as ligands for AhR receptors expressed on gut epithelia and many types of immune cells and may have a significant role in modifying the host mucosal immune system to promote the survival of commensal microbiota and provide protection against pathogens . Modulation of the levels of tryptophan, serotonin and kynurenine by microbiota may have effects on both the enteral and central nervous systems and the immune system. Disturbances in the intestinal microbiota can affect the levels of these metabolites, which may contribute to the pathogenesis of a multitude of disorders, such as inflammatory bowel diseases, metabolic syndrome and obesity, and various neuropsychiatric disorders . 3-indolepropionic acid is a strong antioxidant [45, 46]. It has been shown to regulate the intestinal barrier function by acting as a ligand to the pregnane X receptor (PXR) . It also has neuroprotective functions  and was recently shown to promote thalamocortical axonogenesis in fetal mouse brain . The changes in the levels of tryptophan metabolites in fetal intestine and brain reported here suggest that the microbial influences may begin already during fetal development.
Methylimidazoleacetic acid is the main metabolite of histamine. In this study, it was not detected in brain of the GF fetuses and its levels were also considerably lower in the other organs compared to SPF mice. Methylimidazoleacetic acid may be associated with miscarriage, potentially by the dysregulation of cytokine networks caused by imbalance in gut bacteria . Another main metabolic pathway of histamine leads to the formation of imidazoleacetic acid. Its levels were decreased in the placenta and fetal intestine of the GF mice. The 4 to 6-fold lower levels in GF mice fetal organs of pipecolic acid, a degradation product of L-lysine produced by intestinal bacteria, may reflect lack of the contribution of gut microbiota to the pipecolic acid pathway of lysine catabolism in GF mice . Other amino acid-related metabolites, including phenylacetylglycine, 1-methylhistidine, 3-methylhistidine, 3-hydroxybutanoic acid, Glu-Gln and Glu-Tyr, also had lower levels in the GF organs, indicating that the presence of microbiota may increase their production. Phenylacetylglycine has been identified as a gut microbial metabolite of phenylalanine with associations to health and disease in humans . It has also been detected in forebrain of mice from neonatal period into adulthood . Certain amino acids and small peptides, such as L-threonine and two tripeptides with unknown structure, had higher levels in the GF organs compared to SPF, suggesting that they were accumulated in the organs due to lack of microbial metabolism.
Plant-derived metabolites catechol-O-sulfate, hippuric acid and solanidine
The mice obtain plant-derived metabolites, such as polyphenols and polyalkaloids, by ingestion of plant-based feed. These molecules can be modified in the microbial metabolism and subsequently absorbed by the host. Catechol-O-sulfate and hippuric acid are metabolites of dietary polyphenols, such as flavonoids and phenolic acids, via degradation by gut microbiota and subsequent sulfation or glycine conjugation in the liver [52, 53]. Hippuric acid is also produced from breakdown products of aromatic amino acids in liver, intestine and kidney for excretion into urine ; gut microbiota is involved in this process . Levels of hippuric acid were significantly reduced in the placenta, fetal intestine and brain. Decreased levels of hippuric acid have been reported in serum of GF mice and in urine of antibiotic-treated rats [4, 54] and recently also in the brain of fetuses of GF or antibiotic-treated mice . Hippuric acid is one of the metabolites reported to be present in neonatal mouse brain with amounts decreasing by age . Apart from participating in removal of benzoate, a metabolite potentially toxic to mitochondria , hippuric acid has been suggested a role in regulation of blood glucose levels, insulin secretion by β-cells and glucose utilization in skeletal muscle [52, 56]. Solanidine is a steroidal glycoalkaloid, which is obtained via dehydroxylation from other glycoalkaloids, such as α-chaconine and α-solanine, which are abundant in potato, a component of the RM3-A-P breeding diet. Solanidine has been detected as the main metabolite in rats after oral ingestion of α-chaconine, and the current findings support the hypothesis that solanidine is a gut microbial metabolite of dietary glycoalkaloids in mice .
Intestinal microbiota is known to affect the energy metabolism of the host. For example, SCFAs produced from dietary fibers by gut microbiota can directly contribute the host lipogenesis and gluconeogenesis. On the other hand, SCFAs and microbiota-produced secondary bile acids also act as regulators increasing insulin sensitivity and satiety thus decreasing energy consumption . SCFAs may also have roles in the programming of fetal metabolism during pregnancy .
There were few lipids among the differentially abundant metabolites in the placenta, fetal intestine, and brain even though lipids were widely represented in the LC-MS data. Moderately higher abundances of some long-chain polyunsaturated fatty acids were measured in SPF mice and higher abundances of palmitic acid in GF mice. It is possible that more consistent and extensive effects of the microbiota on the lipid metabolite profiles would be detected in fetal tissues central in lipid metabolism, such as the liver, muscle, and adipose tissue.
We found increased levels of carnitine and various acylcarnitines in placenta and fetal organs of GF mice. These are precursors of TMAO , which suggests that they were not metabolized in the gut of the dam due to the lack of gut microbiota. Increased levels of carnitine could also be due to decreased levels of 5-AVAB in the GF mice, since 5-AVAB inhibits the cell membrane carnitine transporter responsible for cellular intake of carnitine . Carnitine and acyl carnitines are involved in beta oxidation of fatty acids. They also have a multifaceted role in the developing brain .
Unannotated molecular features undetected in GF mice
In total 3680 molecular features were differentially expressed in SPF and GF organs and 99 were not detected in GF mice. Most of these features could not yet be annotated. Unidentified metabolites strongly affected by the presence of maternal microbiota may contain hitherto unknown molecular species with potential bioactivity in the mammalian host, as was shown recently by identification of novel amino acid conjugations of bile acids by extensive bioinformatic analysis of metabolomics data from GF and SPF mice .
Limitations of the study
The LC-MS method used to detect the metabolites in this study does not allow the detection of short-chain fatty acids (formic to valeric acid), many volatile organic compounds (VOCs) and molecules above 1500 Da (including lipopolysaccharides, large peptides, proteins, and nucleic acids). The SCFAs have been extensively studied already previously [6, 10,11,12, 43]. SCFAs are known to cross the placenta to fetal tissues and have an impact on the fetal development [61, 62].
The current databases available for metabolite identification contain the most well-known mammalian endogenous metabolites. Reference data on exogenous (e.g., plant-derived) and microbially produced metabolites is limited. Therefore, a significant proportion of molecular features in the dataset could not be annotated. A tentative characterization of the chemical class was acquired for 10 of the unknowns with the prediction of molecular formula and comparison to in silico generated MS/MS fragmentation patterns.
The germ-free status of the dam strongly affects the metabolic landscape of the placenta and the developing fetus. Several known metabolites originating in microbial metabolism or combined metabolism of microbiota and the host were among the differential compounds, further highlighting the impact of gut microbiota on the host metabolism. Intriguingly, the most significantly differing metabolites included 5-AVAB and TMAO with reported effects on neural development in mice. Among these were also the tryptophan metabolite serotonin, an important regulator of neural and enteroendocrine functions, and kynurenine and HIAA that can act as effectors of the immune system through ArH receptors. Additionally, hundreds of unannotated molecular features were significantly more abundant in SPF mice or missing from GF mice, indicating that a majority of microbially processed metabolites in the fetus and placenta are still unknown. These represent a pool of compounds with potentially important bioactivities in the host and effects on the fetal development of the brain and other organ systems.
Fetal and placental mouse organ samples
Fetal and placental mouse organ samples from six pregnant germ-free (GF) and six specific pathogen free (SPF) C57BL/6 J dams were obtained from the EMMA Axenic Service at Instituto Gulbenkian de Ciência, Portugal. The GF dams were 3–4 months old and the SPF dams 4–5 months old. Obtaining exactly matching age groups was not possible without wasting pregnant animals, which would violate the 3R principles of animal experimentation. Significant physiological differences between groups are not likely, as both groups represent adult dams at optimal reproductive age, litter sizes were similar (in average, 8.5 pups per SPF dam and 8.33 pups per GF dam), and the sizes of the fetuses were not significantly different when controlled for litter size. The GF and SPF statuses were regularly monitored by culture and 16S qPCR (for details, see Supplementary Information: Mouse health surveillance program; the SPF dams came from the Production Area and were housed in E0 until sampling). All dams were fed identical RM3-A-P (Rat and Mouse No.3 Breeding Autoclavable) breeding diets (SDS Special Diet Services, Essex, UK), autoclaved at 121 °C. The SPF feed was autoclaved for 20 min and the GF feed for 30 min due to logistical reasons.
The dams were euthanized 18.5 days post coitum by cervical dislocation. Whole brain, intestine and placenta were collected from two fetuses per dam, to a total of 12 fetuses per experimental group. The fetal organ samples were immediately frozen in liquid nitrogen, stored at − 80 °C and transferred on dry ice to the research laboratory. In both groups, there were 7 male fetuses and 5 female fetuses.
The number of experimental animals was estimated to be sufficient based on our previous metabolomics research. No animals were excluded from the analyses. The dams were selected from the colony by the presence of copulatory plugs and pregnancy. Confounders were not specifically controlled; no obvious confounders to this study were observed.
Metabolomics analyses were performed at Afekta Technologies Ltd. Frozen organ samples were thawed at 8 °C for 2 h and then weighed (approx. 100 mg) in homogenizer tubes. For the metabolite extraction, cold methanol (80% v/v) was added in a ratio of 500 μl per 100 mg of sample. The samples were homogenized (TissueLyser II bead mill, Qiagen, Hilden, Germany) using metal beads at 6 m/s for 30 s. The samples were then shaken for 5 min in room temperature and centrifuged at 14,000 rpm at 4 °C for 10 min. After the centrifugation, the samples were kept on ice for 5 to 10 min, after which the supernatant was filtered (Acrodisc 0.2 μm PTFE membrane, Pall, Port Washington, NY, USA) into HPLC vials for analysis. The pooled quality control (QC) sample was prepared by collecting 20 μl from each sample vial and combining the material to two vials.
The laboratory staff was blinded to the identity of the samples and was not aware of the experimental setup or grouping. The samples were randomized to analysis batches.
The samples were analyzed by liquid chromatography–mass spectrometry (LC-MS), consisting of a 1290 Infinity Binary UPLC coupled with a 6540 UHD Accurate-Mass Q-TOF (Agilent Technologies Inc., Santa Clara, CA, USA), as described previously . In brief, a Zorbax Eclipse XDB-C18 column (2.1 × 100 mm, 1.8 μm; Agilent Technologies) was used for the reversed-phase (RP) separation and an Acquity UPLC BEH amide column (Waters Corporation, Milford, MA, USA) for the HILIC separation. After each chromatographic run, the ionization was carried out using jet stream electrospray ionization (ESI) in the positive and negative mode, yielding four data files per sample. The collision energies for the MS/MS (tandem mass spectrometry) analysis were selected as 10, 20 and 40 V, for compatibility with spectral databases.
Peak detection and alignment were performed in MS-DIAL ver. 4.00 . For the peak collection, m/z values between 50 and 1500 and all retention times were considered. The amplitude of minimum peak height was set at 2000. The peaks were detected using the linear weighted moving average algorithm. For the alignment of the peaks across samples, the retention time tolerance was 0.05 min, and the m/z tolerance was 0.015 Da. The heatmaps were produced with Multiple Experiment Viewer (MeV) version 4.9.0. MetaboAnalyst 4.0 was used for the pathway analysis of the annotated metabolites . Data clean-up (for each mode separately) and statistics (for all signals remaining after clean-up) were performed in R version 3.5.1. Molecular features were only kept if they met all the following quality metrics criteria: low number of missing values, present in more than 70% of the QC samples, present in at least 60% of samples in at least one study group, RSD* (the non-parametric version of relative standard deviation) below 20%, D-ratio* (non-parametric measure of the spread of the QC samples compared to the biological samples) below 10%. In addition, if either RSD* or D-ratio* was above the threshold, the features were still kept if their classic RSD, RSD* and basic D-ratio were all below 10%. Low-quality features were flagged and discarded from statistical analyses. Drift correction was applied to the data.
The cleaned data matrices of the four modes were combined before imputation. Features were then imputed using random forest imputation with an OOB error of 0.009. QC samples were removed prior to imputation to prevent them from biasing the procedure. Differential features between the treatment (GF) and control (SPF) were determined using a simple linear model (Student’s t-test) fit separately for each feature. The results were adjusted for multiple comparisons using Benjamini–Hochberg false discovery rate (FDR). FDR-adjusted p values (q values) below 0.05 were considered significant.
For the MS Peaks to Pathways analysis in MetaboAnalyst 4.0 , the data was first normalized by medians, cube root transformed, automatically scaled, and parametric statistical significances calculated with equal variances and adjusted p value (FDR) cutoff 0.05. In Peaks to Pathways, the molecular weight tolerance was set to 10 ppm, primary ions enforced, and adducts set based on the experimental data. Gene Set Enrichment analysis  and Mummichog version 1.0.10  were used.
The chromatographic and mass spectrometric characteristics (retention time, exact mass, and MS/MS spectra) of the significantly differential molecular features were compared with entries in an in-house standard library and publicly available databases, such as METLIN and HMDB, as well as with published literature. The annotation of each metabolite, and the level of identification was given based on the recommendations published by the Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI): level 1 refers to confirmed identifications based on reference standards analyzed with the same instrument with identical conditions; level 2 means putative annotations with matching m/z and MS/MS fragmentation spectra with publicly available databases; level 3 signifies a putative characterization of compound class based on the observed physicochemical characteristics of the molecular feature; and level 4 covers all the remaining (unknown) signals . Level 1–3 annotations included comparison to reference MS/MS spectra. The small peptides were annotated based on exact mass and comparison of the MS/MS spectra with publicly available databases, such as METLIN. The order of the amino acid residues could be deduced from the MS/MS spectra. Level 2 annotation was also given if MS/MS data was not available but the retention time and calculated molecular formula matched with that of a reference standard.
Availability of data and materials
The metabolomics dataset for this study can be found at EUDAT: https://doi.org/10.23728/b2share.4be0ea9f87b84a06be960d6a1c4b0b42.
- 5-AVAB, also known as δ-valerobetaine (δVB), N,N,N-trimethyl 5-aminovalerate (TMAV) and N,N,N-trimethyl-5-aminovaleric acid (TMAVA):
5-aminovaleric acid betaine
Aryl hydrocarbon receptor
Chemical Analysis Working Group
Gene set enrichment analysis
Jet stream electrospray ionization
Kyoto Encyclopedia of Genes and Genomes
Liquid chromatography–mass spectrometry
Metabolomics Standards Initiative
Pregnane X receptor
Rat and Mouse No.3 Breeding Autoclavable
Specific pathogen free
Short chain fatty acids
Tandem mass spectrometry
t-Distributed stochastic neighbor embedding
Ultra-high performance liquid chromatography
Zhang LS, Davies SS. Microbial metabolism of dietary components to bioactive metabolites: opportunities for new therapeutic interventions. Genome Med. 2016;8(1):46. https://doi.org/10.1186/s13073-016-0296-x.
Lloyd-Price J, Mahurkar A, Rahnavard G, Crabtree J, Orvis J, Hall AB. Strains, functions and dynamics in the expanded human microbiome project. Nature. 2017;550:61–6.
Sender R, Fuchs S, Milo R. Are we really vastly outnumbered? Revisiting the ratio of bacterial to host cells in humans. Cell. 2016;164(3):337–40. https://doi.org/10.1016/j.cell.2016.01.013.
Wikoff WR, Anfora AT, Liu J, Schultz PG, Lesley SA, Peters EC, et al. Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites. Proc Natl Acad Sci U S A. 2009;106(10):3698–703. https://doi.org/10.1073/pnas.0812874106.
Quinn RA, Melnik AV, Vrbanac A, Fu T, Patras KA, Christy MP, et al. Global chemical effects of the microbiome include new bile-acid conjugations. Nature. 2020;579(7797):123–9. https://doi.org/10.1038/s41586-020-2047-9.
LeBlanc JG, Chain F, Martin R, Bermudez-Humaran LG, Courau S, Langella P. Beneficial effects on host energy metabolism of short-chain fatty acids and vitamins produced by commensal and probiotic bacteria. Microb Cell Factories. 2017;16(1):79. https://doi.org/10.1186/s12934-017-0691-z.
Bergman EN. Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol Rev. 1990;70:567–90.
Arpaia N, Campbell C, Fan X, Dikiy S, van der Veeken J, deRoos P, et al. Metabolites produced by commensal bacteria promote peripheral regulatory T-cell generation. Nature. 2013;504(7480):451–5. https://doi.org/10.1038/nature12726.
Furusawa Y, Obata Y, Fukuda S, Endo TA, Nakato G, Takahashi D, et al. Commensal microbe-derived butyrate induces the differentiation of colonic regulatory T cells. Nature. 2013;504(7480):446–50. https://doi.org/10.1038/nature12721.
Inan MS, Rasoulpour RJ, Yin L, Hubbard AK, Rosenberg DW, Giardina C. The luminal short-chain fatty acid butyrate modulates NF-kappaB activity in a human colonic epithelial cell line. Gastroenterology. 2000;118(4):724–34. https://doi.org/10.1016/s0016-5085(00)70142-9.
Roediger WEW, Moore A. Effect of short-chain fatty acid on sodium absorption in isolated human colon perfused through the vascular bed. Dig Dis Sci. 1981;26:100–6.
Smith PM, Howitt MR, Panikov N, Michaud M, Gallini CA, Bohlooly YM, et al. The microbial metabolites, short-chain fatty acids, regulate colonic Treg cell homeostasis. Science. 2013;341(6145):569–73. https://doi.org/10.1126/science.1241165.
Yano JM, Yu K, Donaldson GP, Shastri GG, Ann P, Ma L, et al. Indigenous bacteria from the gut microbiota regulate host serotonin biosynthesis. Cell. 2015;161(2):264–76. https://doi.org/10.1016/j.cell.2015.02.047.
Koeth RA, Wang Z, Levison BS, Buffa JA, Org E, Sheehy BT, et al. Intestinal microbiota metabolism of L-carnitine, a nutrient in red meat, promotes atherosclerosis. Nat Med. 2013;19(5):576–85. https://doi.org/10.1038/nm.3145.
Kärkkäinen O, Tuomainen T, Koistinen V, Tuomainen M, Leppanen J, Laitinen T, et al. Whole grain intake associated molecule 5-aminovaleric acid betaine decreases beta-oxidation of fatty acids in mouse cardiomyocytes. Sci Rep. 2018;8(1):13036. https://doi.org/10.1038/s41598-018-31484-5.
Wang Z, Klipfell E, Bennett BJ, Koeth R, Levison BS, Dugar B, et al. Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature. 2011;472(7341):57–63. https://doi.org/10.1038/nature09922.
Mardinoglu A, Shoaie S, Bergentall M, Ghaffari P, Zhang C, Larsson E, et al. The gut microbiota modulates host amino acid and glutathione metabolism in mice. Mol Syst Biol. 2015;11(10):834. https://doi.org/10.15252/msb.20156487.
Yamamoto Y, Nakanishi Y, Murakami S, Aw W, Tsukimi T, Nozu R, et al. A metabolomic-based evaluation of the role of commensal microbiota throughout the gastrointestinal tract in mice. Microorganisms. 2018;6(4):101. https://doi.org/10.3390/microorganisms6040101.
Agus A, Planchais J, Sokol H. Gut microbiota regulation of tryptophan metabolism in health and disease. Cell Host Microbe. 2018;23(6):716–24. https://doi.org/10.1016/j.chom.2018.05.003.
Gomez de Aguero M, Ganal-Vonarburg SC, Fuhrer T, Rupp S, Uchimura Y, Li H, et al. The maternal microbiota drives early postnatal innate immune development. Science. 2016;351(6279):1296–302. https://doi.org/10.1126/science.aad2571.
Walker RW, Clemente JC, Peter I, Loos RJF. The prenatal gut microbiome: are we colonized with bacteria in utero? Pediatr Obes. 2017;12(Suppl 1):3–17. https://doi.org/10.1111/ijpo.12217.
Ganal-Vonarburg SC, Hornef MW, Macpherson AJ. Microbial-host molecular exchange and its functional consequences in early mammalian life. Science. 2020;368(6491):604–7. https://doi.org/10.1126/science.aba0478.
van de Pavert SA, Ferreira M, Domingues RG, Ribeiro H, Molenaar R, Moreira-Santos L, et al. Maternal retinoids control type 3 innate lymphoid cells and set the offspring immunity. Nature. 2014;508(7494):123–7. https://doi.org/10.1038/nature13158.
Grizotte-Lake M, Zhong G, Duncan K, Kirkwood J, Iyer N, Smolenski I, et al. Commensals suppress intestinal epithelial cell retinoic acid synthesis to regulate Interleukin-22 activity and prevent microbial Dysbiosis. Immunity. 2018;49(6):1103–15.e6. https://doi.org/10.1016/j.immuni.2018.11.018.
Kimura I, Miyamoto J, Ohue-Kitano R, Watanabe K, Yamada T, Onuki M, et al. Maternal gut microbiota in pregnancy influences offspring metabolic phenotype in mice. Science. 2020;367(6481). https://doi.org/10.1126/science.aaw8429.
Horn DL, Morrison DC, Opal SM, Silverstein R, Visvanathan K, Zabriskie JB. What are the microbial components implicated in the pathogenesis of sepsis? Report on a symposium. Clin Infect Dis. 2000;31:851–8.
Chong J, Soufan O, Li C, Caraus I, Li S, Bourque G, et al. MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018;46(W1):W486–W94. https://doi.org/10.1093/nar/gky310.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14. https://doi.org/10.1093/nar/gkr988.
Karp PD, Billington R, Caspi R, Fulcher CA, Latendresse M, Kothari A, et al. The BioCyc collection of microbial genomes and metabolic pathways. Brief Bioinform. 2019;20(4):1085–93. https://doi.org/10.1093/bib/bbx085.
Koistinen VM, Kärkkäinen O, Borewicz K, Zarei I, Jokkala J, Micard V, et al. Contribution of gut microbiota to metabolism of dietary glycine betaine in mice and in vitro colonic fermentation. Microbiome. 2019;7(1):103. https://doi.org/10.1186/s40168-019-0718-2.
Raman SB, Rathinasabapathi B. beta-alanine N-methyltransferase of Limonium latifolium. cDNA cloning and functional expression of a novel N-methyltransferase implicated in the synthesis of the osmoprotectant beta-alanine betaine. Plant Physiol. 2003;132(3):1642–51. https://doi.org/10.1104/pp.103.020453.
Velasquez MT, Ramezani A, Manal A, Raj DS. Trimethylamine N-oxide: the good, the bad and the unknown. Toxins (Basel). 2016;8(11):326. https://doi.org/10.3390/toxins8110326.
Servillo L, D'Onofrio N, Giovane A, Casale R, Cautela D, Castaldo D, et al. Ruminant meat and milk contain delta-valerobetaine, another precursor of trimethylamine N-oxide (TMAO) like gamma-butyrobetaine. Food Chem. 2018;260:193–9. https://doi.org/10.1016/j.foodchem.2018.03.114.
Jääskeläinen T, Kärkkäinen O, Jokkala J, Klavus A, Heinonen S, Auriola S, et al. A non-targeted LC-MS metabolic profiling of pregnancy: longitudinal evidence from healthy and pre-eclamptic pregnancies. Metabolomics. 2021;17(2):20. https://doi.org/10.1007/s11306-020-01752-5.
Jääskeläinen T, Kärkkäinen O, Jokkala J, Litonius K, Heinonen S, Auriola S, et al. A non-targeted LC-MS profiling reveals elevated levels of carnitine precursors and Trimethylated compounds in the cord plasma of pre-Eclamptic infants. Sci Rep. 2018;8(1):14616. https://doi.org/10.1038/s41598-018-32804-5.
Vuong HE, Pronovost GN, Williams DW, Coley EJL, Siegler EL, Qiu A, et al. The maternal microbiome modulates fetal neurodevelopment in mice. Nature. 2020;586(7828):281–6. https://doi.org/10.1038/s41586-020-2745-3.
al-Waiz M, Mikov M, Mitchell SC, Smith RL. The exogenous origin of trimethylamine in the mouse. Metabolism. 1992;41(2):135–6. https://doi.org/10.1016/0026-0495(92)90140-6.
Lang DH, Yeung CK, Peter RM, Ibarra C, Gasser R, Itagaki K, et al. Isoform specificity of trimethylamine N-oxygenation by human flavin-containing monooxygenase (FMO) and P450 enzymes: selective catalysis by FMO3. Biochem Pharmacol. 1998;56(8):1005–12. https://doi.org/10.1016/s0006-2952(98)00218-4.
Landfald B, Valeur J, Berstad A, Raa J. Microbial trimethylamine-N-oxide as a disease marker: something fishy? Microb Ecol Health Dis. 2017;28(1):1327309. https://doi.org/10.1080/16512235.2017.1327309.
Tang WH, Wang Z, Levison BS, Koeth RA, Britt EB, Fu X, et al. Intestinal microbial metabolism of phosphatidylcholine and cardiovascular risk. N Engl J Med. 2013;368(17):1575–84. https://doi.org/10.1056/NEJMoa1109400.
Zhu W, Gregory JC, Org E, Buffa JA, Gupta N, Wang Z, et al. Gut microbial metabolite TMAO enhances platelet hyperreactivity and thrombosis risk. Cell. 2016;165(1):111–24. https://doi.org/10.1016/j.cell.2016.02.011.
Clarke G, McKernan DP, Gaszner G, Quigley EM, Cryan JF, Dinan TG. A distinct profile of tryptophan metabolism along the kynurenine pathway downstream of toll-like receptor activation in irritable bowel syndrome. Front Pharmacol. 2012;3:90. https://doi.org/10.3389/fphar.2012.00090.
Reigstad CS, Salmonson CE, Rainey JF 3rd, Szurszewski JH, Linden DR, Sonnenburg JL, et al. Gut microbes promote colonic serotonin production through an effect of short-chain fatty acids on enterochromaffin cells. FASEB J. 2015;29(4):1395–403. https://doi.org/10.1096/fj.14-259598.
Zelante T, Iannitti RG, Cunha C, De Luca A, Giovannini G, Pieraccini G, et al. Tryptophan catabolites from microbiota engage aryl hydrocarbon receptor and balance mucosal reactivity via interleukin-22. Immunity. 2013;39(2):372–85. https://doi.org/10.1016/j.immuni.2013.08.003.
Chyan YJ, Poeggeler B, Omar RA, Chain DG, Frangione B, Ghiso J, et al. Potent neuroprotective properties against the Alzheimer beta-amyloid by an endogenous melatonin-related indole structure, indole-3-propionic acid. J Biol Chem. 1999;274(31):21937–42. https://doi.org/10.1074/jbc.274.31.21937.
Hardeland R, Zsizsik BK, Poeggeler B, Fuhrberg B, Holst S, Coto-Montes A. Indole-3-pyruvic and -propionic acids, kynurenic acid, and related metabolites as luminophores and free-radical scavengers. In: Huether G, Kochen W, Simat TJ, Steinhart H, editors. Tryptophan, serotonin, and melatonin. Boston: Springer US; 1999. p. 389–95.
Venkatesh M, Mukherjee S, Wang H, Li H, Sun K, Benechet AP, et al. Symbiotic bacterial metabolites regulate gastrointestinal barrier function via the xenobiotic sensor PXR and toll-like receptor 4. Immunity. 2014;41(2):296–310. https://doi.org/10.1016/j.immuni.2014.06.014.
Liu Y, Chen H, Chen D, Feng L, Zhang J. Gut Dysbacteriosis is associated with an imbalanced cytokines network in women with unexplained miscarriage. Review. 2020. https://doi.org/10.21203/rs.3.rs-20415/v1.
He M. Pipecolic acid in microbes: biosynthetic routes and enzymes. J Ind Microbiol Biotechnol. 2006;33(6):401–7. https://doi.org/10.1007/s10295-006-0078-3.
Nemet I, Saha PP, Gupta N, Zhu W, Romano KA, Skye SM, et al. A cardiovascular disease-linked gut microbial metabolite acts via adrenergic receptors. Cell. 2020;180(5):862–77.e22. https://doi.org/10.1016/j.cell.2020.02.016.
Swann JR, Spitzer SO, Diaz Heijtz R. Developmental signatures of microbiota-derived metabolites in the mouse brain. Metabolites. 2020;10(5):172. https://doi.org/10.3390/metabo10050172.
de Mello VD, Lankinen MA, Lindstrom J, Puupponen-Pimia R, Laaksonen DE, Pihlajamaki J, et al. Fasting serum hippuric acid is elevated after bilberry (Vaccinium myrtillus) consumption and associates with improvement of fasting glucose levels and insulin secretion in persons at high risk of developing type 2 diabetes. Mol Nutr Food Res. 2017;61(9):1700019. https://doi.org/10.1002/mnfr.201700019.
Feliciano RP, Boeres A, Massacessi L, Istas G, Ventura MR, Nunes Dos Santos C, et al. Identification and quantification of novel cranberry-derived plasma and urinary (poly)phenols. Arch Biochem Biophys. 2016;599:31–41. https://doi.org/10.1016/j.abb.2016.01.014.
Lee SH, An JH, Park HM, Jung BH. Investigation of endogenous metabolic changes in the urine of pseudo germ-free rats using a metabolomic approach. J Chromatogr B Analyt Technol Biomed Life Sci. 2012;887-888:8–18. https://doi.org/10.1016/j.jchromb.2011.12.030.
Badenhorst CP, Erasmus E, van der Sluis R, Nortje C, van Dijk AA. A new perspective on the importance of glycine conjugation in the metabolism of aromatic acids. Drug Metab Rev. 2014;46(3):343–61. https://doi.org/10.3109/03602532.2014.908903.
Bitner BF, Ray JD, Kener KB, Herring JA, Tueller JA, Johnson DK, et al. Common gut microbial metabolites of dietary flavonoids exert potent protective activities in beta-cells and skeletal muscle cells. J Nutr Biochem. 2018;62:95–107. https://doi.org/10.1016/j.jnutbio.2018.09.004.
Norred WP, Nishie K, Osman SF. Excretion, distribution and metabolic fate of 3H-alpha-chaconine. Res Commun Chem Pathol Pharmacol. 1976;13(2):161–71.
Yu Y, Raka F, Adeli K. The role of the gut microbiota in lipid and lipoprotein metabolism. J Clin Med. 2019;8(12). https://doi.org/10.3390/jcm8122227.
Zietek M, Celewicz Z, Szczuko M. Short-chain fatty acids, maternal microbiota and metabolism in pregnancy. Nutrients. 2021;13(4):1244. https://doi.org/10.3390/nu13041244.
Ferreira GC, McKenna MC. L-carnitine and acetyl-L-carnitine roles and neuroprotection in developing brain. Neurochem Res. 2017;42(6):1661–75. https://doi.org/10.1007/s11064-017-2288-7.
Priyadarshini M, Thomas A, Reisetter AC, Scholtens DM, Wolever TM, Josefson JL, et al. Maternal short-chain fatty acids are associated with metabolic parameters in mothers and newborns. Transl Res. 2014;164(2):153–7. https://doi.org/10.1016/j.trsl.2014.01.012.
Shekhawat P, Bennett MJ, Sadovsky Y, Nelson DM, Rakheja D, Strauss AW. Human placenta metabolizes fatty acids: implications for fetal fatty acid oxidation disorders and maternal liver diseases. Am J Physiol Endocrinol Metab. 2003;284(6):E1098–105. https://doi.org/10.1152/ajpendo.00481.2002.
Klåvus A, Kokla M, Noerman S, Koistinen VM, Tuomainen M, Zarei I, et al. “notame”: workflow for non-targeted LC-MS metabolic profiling. Metabolites. 2020;10(4):135. https://doi.org/10.3390/metabo10040135.
Tsugawa H, Cajka T, Kind T, Ma Y, Higgins B, Ikeda K, et al. MS-DIAL: data-independent MS/MS deconvolution for comprehensive metabolome analysis. Nat Methods. 2015;12(6):523–6. https://doi.org/10.1038/nmeth.3393.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50. https://doi.org/10.1073/pnas.0506580102.
Li S, Park Y, Duraisingham S, Strobel FH, Khan N, Soltow QA, et al. Predicting network activity from high throughput metabolomics. PLoS Comput Biol. 2013;9(7):e1003123. https://doi.org/10.1371/journal.pcbi.1003123.
Sumner LW, Amberg A, Barrett D, Beale MH, Beger R, Daykin CA, et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics. 2007;3(3):211–21. https://doi.org/10.1007/s11306-007-0082-2.
We thank Kirsi Lahti and Santeri Suokas for expert technical assistance and Anton Klåvus for running the statistical analyses.
This study was funded by the Finnish Cultural Foundation. The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript. Open access funded by Helsinki University Library.
Ethics approval and consent to participate
The research was carried out in compliance with European Directive 2010/63/EU and Portuguese Law. The EMMA Axenic Service has an approval for each research project and for the animal euthanization method (cervical dislocation without anesthesia) by the Ethics Committee of Instituto Gulbenkian de Ciência, the Animal Welfare Body, and the National Competent Authority (DGAV), and follows the recommendations of the Federation of Laboratory Animal Science Associations (FELASA). The Animal Facility is licensed by Direcção Geral de Veterinária (N° of identification and registration: PT 05 052 OICB). The procedures in this study (breeding of non-harmful phenotypes and organ collection upon euthanasia) are not considered experimental animal procedures (2010/63/EU and Portuguese Law).
Consent for publication
The authors declare that they have no financial or non-financial competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Pessa-Morikawa, T., Husso, A., Kärkkäinen, O. et al. Maternal microbiota-derived metabolic profile in fetal murine intestine, brain and placenta. BMC Microbiol 22, 46 (2022). https://doi.org/10.1186/s12866-022-02457-6