Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 26 November 2025

Sec. Microorganisms in Vertebrate Digestive Systems

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1668148

This article is part of the Research TopicNew and advanced mechanistic insights into the influences of the infant gut microbiota on human health and disease, Volume IIView all 20 articles

Integrated multi-omics reveals distinct maternal and neonatal gut microbial and metabolic signatures associated with small for gestational age status

  • 1State Key Laboratory for Quality and Safety of Agro-products & Food Sciences Institute, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
  • 2Obstetrical Department of Hangzhou Linping District Maternal & Child Health Care Hospital, Hangzhou, Zhejiang, China

Introduction: Children born small for gestational age (SGA) have an elevated risk of developing metabolic disorders in later life. However, the underlying gut microbiota and metabolomic alterations in SGA mother-infant dyads remain poorly understood.

Methods: We performed an integrated analysis of fecal metagenomics, metabolomics, and short-chain fatty acids (SCFAs) profiling in 10 SGA and 10 appropriate for gestational age (AGA) mother-infant dyads at term. Taxonomic composition, microbial functional pathways, carbohydrate-active enzyme (CAZyme) profiles, differential metabolites, and metabolite pathway enrichment were systematically evaluated.

Results and discussion: SGA neonates exhibited reduced microbial richness (Chao1 index), distinct beta-diversity, and differential abundance of key bacterial species including increased Enterococcus faecalis and Escherichia coli. Functionally, SGA maternal subjects showed divergent profiles in CAZyme genes, with lower abundance of glycoside hydrolase family 13 subfamily 16, glycosyl transferase family 66, and carbohydrate-binding module family 6, and altered structural polysaccharide degradation capacity. Metabolomic profiling revealed significant perturbations in tryptophan metabolism pathways, notably enriched in kynurenine and taurine derivatives in SGA mother and neonates. Notably, SCFA profiles were disrupted, with increased butyrate in SGA mother and reduced propionate and isobutyrate in SGA neonates. Microbe-metabolite correlation networks revealed strong associations between SGA-specific bacterial taxa and fecal metabolites. In conclusion, our analysis identifies distinct features of the early fecal microbiome and metabolome within 48 h of birth in SGA neonates compared with AGA peers, reflecting differences in initial colonization and metabolism that warrant longitudinal follow-up.

Introduction

Low birth weight has been defined as weight at birth of <2,500 g (5.5 pounds) by the World Health Organization. Small for gestational age (SGA), newborns whose birth weight below the 10th percentile or 2 standard deviations of the many other babies of the same gestational age (Alkalay et al., 1998; Usher and McLean, 1969). SGA infants have become a significant public health concern due to their higher risk of complications compared to appropriate for gestational age (AGA) infants. United Nations International Children's Emergency Fund (2023) reports that 14.7% of newborns globally in 2020 had low birth weight, and its incidence was 17.0% in Asia. The pathogenesis of SGA infants involves multiple factors, including genetics, placental dysfunction, and maternal nutritional status. SGA neonates face higher risks of both short-term and long-term health complications, such as neonatal distress, feeding intolerance, necrotizing enterocolitis, insulin resistance, and an increased likelihood of metabolic disorders in adulthood. Understanding the underlying mechanisms of SGA development is crucial for early interventions and improving long-term health outcomes.

The gut microbiota plays a crucial role in infant health, and its composition, diversity, and abundance vary among neonates with different physiological conditions. Chen et al. (2022) analyzed the changes in gut microbiota of SGA infants from birth to 6 months and found that certain pathogenic and conditional pathogenic bacteria, such as Escherichia_Shigella, Ralstonia and Clostridium increased. SGA neonates also exhibit lower microbial diversity and a reduced proportion of beneficial bacteria, including Klebsiella and Enterobacter (Chang et al., 2022). Dietary changes throughout infancy are consistent with increasing demands for carbohydrate-active enzyme (CAZyme). The complex structures of human milk oligosaccharides (HMOs) necessitate a specific repertoire of glycoside hydrolases (GHs), notably those belonging to the GH13 and GH26 families, which are abundant in Bacteroidetes and Bifidobacterium species (Shoemaker et al., 2001; Xu et al., 2024). CAZyme profiling in SGA infant feces provides functional insight into gut microbial adaptation to nutritional challenges, potentially revealing enzymatic deficiencies or metabolic adaptations that influence growth and health outcomes. Additionally, changes in microbiome were accompanied by changes in fecal metabolites. SGA neonates exhibit increased levels of hexanoylcarnitine alongside reduced abundances of 3-hydroxybutyrate and 3-iminopropanoate, indicating disruptions in taurine metabolism and ketone body pathways (Wang et al., 2022). The perinatal period is a critical window for the establishment of the infant gut microbiota. During the perinatal period, microorganisms in the maternal gut can share genes with those in the infant gut, and this process begins even before birth and continues for several weeks postpartum (Vatanen et al., 2022). This horizontal gene transfer allows maternal microbes to influence the functional properties of the infant's microbiome, even though the microbial strains themselves are not continuously transmitted. However, there is currently no direct evidence confirming whether differences in initial microbial colonization in SGA infants are associated with insufficient or abnormal transmission of the maternal microbiota.

This study aims to investigate the metagenomic and untargeted metabolomic characteristics of the feces in maternal-infant dyads. In addition, to clarify the role of the gut microbiota and metabolites in the development of SGA, we performed gut microbiota functional prediction and correlation analysis between differential microbiota and fecal metabolites. The findings of this study will provide a theoretical basis for the early prediction and effective intervention of SGA infants.

Subjects and methods

Subjects

A case-control study design was employed, with women delivering SGA infants matched 1:1 with those delivering AGA infants. Matching criteria included full-term birth (defined as a gestational age between 37 and 42 weeks), body mass index (BMI, kg/m2) within the normal range, maternal age within ±5 years, and vaginal delivery. SGA and AGA neonates and their mother were recruited from the pediatric neonatal ward of Hangzhou City Linping District Maternal and Child Care Hospital between July 2021 and February 2022. A total of 20 eligible pregnant women were enrolled, including 10 mothers of SGA infants (SGA_M group) and 10 mother of AGA infants (AGA_M group). Correspondingly, 20 newborns were included, with 10 SGA infants assigned to the SGA_S group and 10 AGA infants to the AGA_S group. Based on intervention data, measurement robustness was quantified as the sample size required to achieve 80% power at α = 0.05 and effect size = 0.60 (two-tailed), calculated a priori using G*Power 3.19.7. All infants in this study were exclusively breastfed, thereby standardizing feeding mode and eliminating formula-related confounding factors.

Inclusion criteria

The following inclusion criteria were used: (a) neonates in the SGA group were required to meet the diagnostic criteria of SGA infant: newborns whose birth weight was less than the 10th percentile of the birth weight for the same gestational age and sex; (b) gestational age was defined as ≥37 and < 42 weeks; (c) newborns without asphyxia, neonatal hypoxic-ischemic encephalopathy, intrauterine infection, severe intracranial hemorrhage, neonatal infectious diseases, chromosomal aneuploidies, and genetic metabolic diseases; and (d) informed consent was obtained from the parent or legal guardian of the infant.

Exclusion criteria

The following exclusion criteria were used: (a) the pregnant women were exposed to antibiotics and probiotics 2 weeks before sample collection; (b) a self-reported special medication history, smoking, or a history of drug or alcohol abuse; (c) presence of severe neonatal conditions such as birth asphyxia, intrauterine infection, hypoxic-ischemic encephalopathy, or neonatal infectious diseases; (d) congenital malformations or genetic disorders in the mother or infant.

Methods

Data collection

Clinical data, including sex, gestational age, birth weight, and mode of delivery were collected. Meconium from neonates and stool samples from pregnant women were collected within 48 h after delivery. The specimens were stored in sterile freezing tubes and stored at −80°C for analysis.

Targeted metabolites measurement of short chain fatty acids (SCFAs)

Quantitative determination of SCFAs was performed as described previously (Li et al., 2025). Briefly, the stool samples were mixed with 50 μL of 15% phosphoric acid, 100 μL of internal standard (isohexanoic acid) solution (125 μg/mL), and 400 μL of ether. The mixture was homogenized for 1 min and then centrifuged at 12,000 rpm and 4°C for 10 min. The supernatant was used for SCFA determination (Trace 1310 GC coupled to a Thermo ISQ mass spectrometer, Thermo Scientific, USA). Chromatographic separation was performed on an Agilent HP-INNOWAX capillary column (30 m × 0.25 mm × 0.25 μm). The concentrations of SCFAs were quantified via a single-point internal standard method. Peak identification and internal response factors were determined using a calibration cocktail containing acetate, propionate, butyrate, isobutyrate, valerate, isovalerate. Concentrations were quantified via an internal standard method. Normalization was performed using a standard curve from QC samples within each batch to eliminate inter-batch differences.

Stool metagenomic sequencing and analysis

Genomic DNA from 50 mg stool samples were extracted using the QIAamp DNA Stool Mini Kit (Qiagen, Dusseldorf, Germany), following the manufacturer's instructions. After fragmenting the extracted stool DNA, whole-genome shotgun metagenomic sequencing was used to sequence total metagenomic DNA on the Illumina TruSeq Nano DNA LT Library Prep at Suzhou PANOMIX Biomedical Tech (Suzhou, China), using a paired-end 150 bp shotgun sequencing strategy and an insert size of 400 bp. DNA libraries were constructed using a fragmentation-based approach. The library was multiplexed and sequenced using Illumina Hiseq X-ten platform. Libraries were quantified by Qubit TM 3.0 Fluometer (ThermoFisher).

Raw sequencing reads were quality-checked using FastQC. Low-quality reads, adapter sequences, and host-derived contamination were removed using Trimmomatic. High-quality reads were retained for downstream analysis. Cleaned reads were assembled into contigs using MEGAHIT, and gene prediction was performed using Prodigal. Taxonomic classification was assigned using Kraken2 against a curated microbial reference database, while functional annotation was performed using the KEGG and EggNOG databases. Species and functional gene abundance tables were normalized using the centered log-ratio transformation. Sequencing batch was included as a covariate in the generalized linear model to control for its effects.

Alpha and beta diversity indices were calculated using QIIME2 to compare microbial composition between groups. Differentially abundant taxa and functional pathways were identified using linear discriminant analysis effect size (LEfSe) and DESeq2, with a significance threshold of p < 0.05. CAZyme family abundances were quantified using Transcripts Per Kilobase Million (TPM) to account for gene length and sequencing depth. Between-sample normalization was performed using DESeq2′s size factors for metagenomic data, and BH-FDR was applied to control multiple testing.

Metabolomic profiling of stool samples

Chemical compounds in maternal stool samples were assessed by untargeted liquid chromatography-mass spectrometry (LC-MS) at Suzhou PANOMIX Biomedical Tech (Suzhou, China). 50 mg stool sample was added to 600 μL L-2-Chlorophenylalanine in methanol and centrifuged at 12,000 rpm at 4°C for 10 min. The supernatant was filtered through a 0.22 μm membrane. Chromatography was carried out with an ACQUITY UPLC® HSS T3 (150 × 2.1 mm, 1.8 μm) (Waters, Milford, MA, USA), with the column maintained at 40°C. The flow rate and injection volume were set at 0.3 mL/min and 2 μL, respectively. The mobile phases for the gradient conditions were performed as described by Eva et al. (Zelena et al., 2009). Mass spectrometric detection of metabolites was performed on Orbitrap Exploris 120 (Thermo Fisher Scientific, USA). Electrospray ionization mass spectrometric (ESI-MS) experiments were performed with spray voltages of 3.5 and −2.5 kV in positive and negative modes. Sheath and auxiliary gas and capillary temperature were set at 30 and 10 arbitrary units and 325°C. The scanning quality range of the analyzer was m/z 81–1,000, and the mass resolution is 60,000. Raw LC-MS data were converted to mzXML format and processed using XCMS for peak detection, retention time alignment, and normalization. Noise filtering and batch effect correction were applied to improve data reliability. Metabolites were annotated by searching the Human Metabolome (HMDB) Database (www.hmdb.ca).

Raw peak area data were log-transformed. Instrumental drift was corrected based on a regression model established from QC samples, and the ComBat algorithm was employed to adjust for batch effects arising from sample preparation. Multivariate statistical analyses, including principal component analysis (PCA) and orthogonal partial least squares-discriminant analysis (OPLS-DA), were conducted to identify metabolic differences between groups. Differential metabolites were selected based on a variable importance in projection (VIP) score >1.0 and p < 0.05. Pathway enrichment analysis was performed using MetaboAnalyst to identify key metabolic pathways associated with the observed changes.

Ethical approval

The study was approved by the Ethics Committee of the Obstetrical Department of Hangzhou Linping District Maternal &Child Health Care Hospital. The legal guardians of each participant provided written informed consent (Ethics ID: LLSC-KYKT-2022-0048-A).

Statistical analyses

Demographic characteristics were calculated and presented as median (upper or lower quartile) for continuous variable and proportions for categorical variables. Categorical variables were compared by means of chi-square test. Continuous variables were compared using Students' t-test. All analyses were performed with IBM SPSS Statistics (version 25.0) and STATA (version 16.0). Two-tailed p-values of < 0.05 were considered statistically significant.

Results

Demographics and clinical characteristics of the participants

A total of 10 SGA neonates were enrolled in the SGA_S group, including 7 males. The AGA_S group included 10 AGA neonates, with 6 males. All neonates were born at a gestational age of 37–42 weeks and were fed a diet of breast milk. The clinical characteristics of the enrolled participants are presented in Table 1. Compared to the AGA_S group, neonates in the SGA_S group had significantly higher gestational age and lower birth weight. However, no significant differences were found between the SGA_M and AGA_M groups in terms of maternal BMI, geographic region, maternal education level, presence of siblings. Furthermore, none of the enrolled neonates were infected with the novel coronavirus, and their mother and family members did not exhibit any pandemic-related psychological or behavioral abnormalities.

Table 1
www.frontiersin.org

Table 1. Demographics and clinical characteristics of the participants.

Compositional profiles of the fecal microbiome and taxonomic differences between the AGA and SGA participants

To assess the diversity of gut microbiota in maternal and newborn subjects, the alpha diversity indices (Shannon, Simpson, and Chao1 index) were calculated. There are no significant differences in the Shannon and Simpson indices between AGA_M and SGA_M groups, as well as AGA_S and SGA_S groups (Figures 1a, b). As shown in Figure 1c, Chao1 index was significantly reduced in SGA_S group as compared to the AGA_S group; there was no significant difference between AGA_M and SGA_M groups. The distribution of maternal and neonatal microorganisms included bacteria, viruses, archaea, and eukaryotes (Figure 1d). Principal Coordinates Analysis (PCoA) by bray-curtis distance showed no significant difference between the bacterial profiles of AGA_M and SGA_M groups, as well as AGA_S and SGA_S groups (Figure 1e). Partial least squares discriminant analysis (PLS-DA) revealed distinct differences between AGA_M and SGA_M groups, as well as AGA_S and SGA_S groups, at the bacterial species level (Figure 1f). Therefore, only bacterial taxa were considered for downstream comparative analysis between the AGA and SGA groups in maternal and infant.

Figure 1
Several panels display statistical graphs comparing microbial diversity indices between groups. Panel (a) shows violin plots for the Shannon Index for AGA_M vs. SGA_M and AGA_S vs. SGA_S. Panel (b) shows similar plots for the Simpson Index. Panel (c) presents violin plots for the Chao1 Index, highlighting significant differences in AGA_S vs. SGA_S. Panel (d) features pie charts showing microbial composition distinguished by categories like Bacteria and Eukaryota for each group. Panel (e) displays PCA plots comparing microbial compositions between groups. Panel (f) contains PLS-DA plots for separating groups by AGA_M vs. SGA_M and AGA_S vs. SGA_S. Graphs indicate clustering patterns and significant differences.

Figure 1. Alpha and beta diversity analysis of microbial communities across different groups. (a) Changes in the Shannon index. (b) Changes in the Simpson index. (c) Changes in the Chao1 index. (d) Alternations of microbial domains between AGA and SGA subjects. (e) PCoA plot of gut microbiota of subjects according to the Bray-Curtis distance. (f) PLS-DA models discriminated microbiome between the AGA and SGA groups. The p-value was calculated using the two-tailed unpaired Students' t-test. *p < 0.05. The p-value interaction between groups was corrected by Bonferroni-correction.

At the phylum level, the relative abundance of Tenericutes increased significantly in SGA_M group when compared to AGA_M group (p < 0.05, by a Wilcoxon rank sum test, Figure 2a, left). Compared to the AGA_S group, the relative abundance of Synergistetes was significantly decreased in the SGA_S group (p < 0.05, by a Wilcoxon rank sum test, Figure 2a, right). At the species level, compared to the AGA_S group, the relative abundance of Desulfovibrio fairfieldensis, Actinomyces sp.HMT 175, Candidatus Saccharibacteria bacterium, and Bacteroides zhangwenhongii were significantly increased, while Pediococcus pentosaceus, Uncultured Tenericutes bacterium, Bifidobacterium angulatum, Prevotella jejuni, Streptococcus sp. NCTC 11567, Citrobacter freundii, and Desulfovibrio desulfuricans were significantly decreased in the SGA_M (p < 0.05, by a Wilcoxon rank sum test, Table 2). Compared to the AGA_S group, the relative abundance of Bacteroides caecimuris, Muribaculaceae bacterium MF13079, Firmicutes bacterium ASF500, Muribaculum gordoncarteri, Duncaniella dubosii, Faecalibacterium duncaniae, Bacillus pacificus, Roseburia intestinalis, Bacteroides nordii, Romboutsia ilealis, Flavonifractor plautii, Roseburia hominis, Parabacteroides goldsteinii, and Granulicatella adiacens were significantly increased, while Clostridioides difficile, Acinetobacter sp. MYb10, Clostridium perfringens, Acinetobacter johnsonii, Streptococcus ilei, and Bifidobacterium breve were significantly decreased in the SGA_S (p < 0.05, by a Wilcoxon rank sum test, Table 3). Additionally, the relative abundance of Enterococcus faecalis and Escherichia coli increased significantly in SGA_S group when compared to AGA_S group (p < 0.05, by a Students' t-test, Figures 2c, d).

Figure 2
Five-part image showing microbiome analysis. Panels (a) and (b) are stacked bar graphs displaying the relative abundance of microbial phyla and species in AGA_M and SGA_M groups. Panels (c) and (d) show bar graphs of relative abundance for *Enterococcus faecalis* and *Escherichia coli*, illustrating significant differences between AGA_S and SGA_S groups, marked with asterisks. Panel (e) shows two LEfSe analysis plots with LDA scores, indicating different microbial groups in AGA_M and SGA_M, as well as AGA_S and SGA_S, with colored bars for group differentiation.

Figure 2. Changes in the taxonomic distribution of AGA and SGA subjects. (a) Composition of the microbial community at the phylum. (b) Composition of the microbial community at the species level. (c) Relative abundance of Enterococcus faecalis. (d) Relative abundance of Escherichia coli. (e) LDA score plot based on LEfSe analysis. The p-value was calculated using the Kruskal-Walli's test. *p < 0.05. The p-value interaction between groups was corrected by Bonferroni-correction. **p < 0.01

Table 2
www.frontiersin.org

Table 2. Significantly differential microbes at species level in the metagenomic sequencing datasets between the AGA_M and SGA_M group.

Table 3
www.frontiersin.org

Table 3. Significantly differential microbes at species level in the metagenomic sequencing datasets between the AGA_S and SGA_S group.

The LEfSe algorithm was used to assess differences in microbial abundance between AGA_M and SGA_M groups, as well as AGA_S and SGA_S groups (Figure 2e). The differential species score chart revealed that Desulfovibrio fairfieldensis was relatively abundant in the AGA_M group, whereas Eubacterium callanderi was more abundant in the SGA_M group. In the AGA_S group, the genus Bacteroides and Bacteroides caecimuris were enriched, while the genus Streptococcus was more prevalent in the SGA_S group (LDA score > 2). Take together, these findings suggest altered microbial ecosystems associated with SGA.

Microbial function differences between the AGA and SGA participants

Functional annotation and KEGG pathway enrichment analyses were performed to characterize the biological functions of the annotated genes. As shown in Figure 3a, most genes were primarily associated with metabolism-related pathways, including amino acid metabolism, carbohydrate metabolism, and energy metabolism. Other enriched categories included cellular processes, environmental information processing, and human diseases, although to a lesser extent. Figure 3b further illustrates the relative abundance of functional categories across different groups. Metabolic pathways were dominant in all groups, with notable contributions from carbohydrate metabolism, amino acid metabolism, and lipid metabolism.

Figure 3
Figure “a” depicts bar charts showing counts of KEGG pathways in AGA_M vs SGA_M and AGA_S vs SGA_S, categorized into cellular processes, metabolism, and more. Figure “b” shows stacked bar charts of relative abundance percentages for pathways and processes in AGA_M vs SGA_M and AGA_S vs SGA_S. Figure “c” presents bar charts of relative abundance for certain groups with key provided.

Figure 3. Functional annotation and KEGG pathway enrichment analysis. (a) KEGG classification of annotated genes across samples, showing predominant enrichment in metabolic pathways, particularly amino acid metabolism, carbohydrate metabolism, and energy metabolism; (b) Relative abundance of functional categories based on KEGG annotations across different groups, with metabolism-related pathways dominating; (c) Proportions of major KEGG categories in each group, highlighting metabolism as the primary functional category.

To further investigate the microbial potential for carbohydrate degradation in the AGA and SGA subjects, we screened the assembled contigs for the presence of CAZymes. The distribution of different CAZyme classes, including GH, glycosyl transferases (GT), carbohydrate-binding modules (CBM), carbohydrate esterases (CE), auxiliary activities (AA), and polysaccharide lyases (PL), is presented in Figure 3c. The most abundant enzyme class in was the GHs followed by GTs. Less abundant were CBMs, CEs, enzymes with AAs, and PLs.

GHs are responsible for the degradation or hydrolysis of glycosidic bonds in carbohydrates. Oligosaccharide-degrading enzymes are major parts of GH families. Compared with the AGA_M group, we found that the abundance of GH15_2 was significantly increased (p < 0.01), and GH20 and GH13_16 were significantly decreased in the SGA_M group (Table 4, p < 0.05). Debranching enzymes were identified belonging to GH5_9, GH95, GH163, GH166, GH117, and GH110 families. GH5_9 and GH163 were more abundant (p < 0.05) in the SGA_S group than in the AGA_S group, and the abundance of GH95, GH166, GH117, and GH110 in the SGA_S group were lower than in the AGA_S group (Table 5, p < 0.05).

Table 4
www.frontiersin.org

Table 4. Significantly differential genes encoding GHs between the AGA_M and SGA_M group.

Table 5
www.frontiersin.org

Table 5. Significantly differential genes encoding GHs between the AGA_S and SGA_S group.

PLs employ unique mechanisms to degrade polysaccharides containing uronic acids. Furthermore, the abundance of PL1 in the SGA_S group were higher than in the AGA_S group (Table 5, p < 0.05).

GTs are responsible for the addition of saccharides onto other biomolecules. The abundance of GT66 in the SGA_M group was lower than in the AGA_M group (Table 4, p < 0.05). Compared with the AGA_S group, we found that the abundance of GT75, GT100, and GT33 genes were significantly decreased in SGA_S group (Table 5).

The CBM domain facilitates the binding of CAZymes to carbohydrate substrates, thereby enhancing enzymatic activity. Compared with the AGA_M group, we found that the abundance of CBM77 gene was significantly increased in SGA_M group (Table 4). The abundance of CBM6 in the SGA_S group was lower than in the AGA_S group (Table 5, p < 0.05). These data indicate divergent microbial capacities for carbohydrate metabolism and structural polysaccharide degradation in the SGA subjects.

Metabolomic variations in AGA and SGA neonatal meconium and maternal feces

To characterize metabolic alterations associated with SGA status, orthogonal partial least squares discriminant analysis (OPLS-DA) was performed in both negative and positive ion modes. Clear group separation was observed between SGA and AGA subjects (Figures 4a, b), indicating distinct metabolic signatures. The OPLS-DA score plots revealed robust model performance between AGA_M and SGA_M groups (R2 = 0.92 and 0.93 for negative and positive ion modes, respectively), while the neonatal samples (AGA_S vs SGA_S) showed moderate separation (R2 = 0.99 and 0.97, respectively), although with relatively low predictive power (Q2 values < 0.05). Although OPLS-DA revealed group separation, the low Q2 value suggests limited predictive performance and a risk of overfitting.

Figure 4
Comparison plots in both negative and positive ion modes for AGA and SGA groups. Panels a and b show OPLS-DA score plots, similarity plots, and loading plots for AGA_M vs SGA_M and AGA_S vs SGA_S. Panel c presents volcano plots for differential expression with highlighted metabolites, indicating significant changes in the AGA vs SGA comparisons.

Figure 4. Changes in the fecal metabolome of AGA and SGA subjects. OPLS-DA score plots, permutation test, S-plots, and loading plots comparing metabolic profiles between AGA and SGA groups under negative ion mode (a) and positive ion mode (b). Symbols representing samples from the two groups are shown in different colors. (c) Significant differential metabolites between the AGA and SGA groups in maternal (left) and newborn (right) subjects, detected by LC-MS. The x-axis shows fold change (log2 ratio scale), and the y-axis shows the negative log10 of p-values (a higher value represents a stronger statistical significance). Significant p-values were corrected using the Benjamini–Hochberg method.

Volcano plot analysis identified multiple significantly dysregulated metabolites between AGA and SGA participants (Figure 4c). One hundred and fifty two differential metabolites were detected between AGA_M and SGA_M groups. A total of 382 differential metabolites were identified between AGA_S and SGA_S groups. Compared to AGA_M group, quantification levels of 3-hydroxykynurenamine, xanthurenic acid, citrulline, and oxoglutaric acid were significantly increased, while acetyl-N-formyl-5-methoxykynurenamine was decreased in SGA_M group (p < 0.05, Table 6). Whereas, in infants, SGA_S subjects exhibited significantly higher levels of xanthurenic acid, indole-3-lactic acid, and citraconic acid and lower level of acetyl-N-formyl-5-methoxykynurenamine and 3-hydroxykynurenamine relative to AGA_S counterparts (p < 0.05, Table 7), Collectively, these results demonstrate that metabolic profiles are significantly altered in SGA maternal and neonatal samples, reflecting systemic metabolic perturbations potentially linked to fetal growth restriction.

Table 6
www.frontiersin.org

Table 6. Significantly differential metabolites between the AGA_M and SGA_M group.

Table 7
www.frontiersin.org

Table 7. Significantly differential metabolites between the AGA_S and SGA_S group.

Altered metabolic pathway enrichment distinguishes SGA and AGA subjects in meconium and fecal samples

To gain mechanistic insight into the metabolic disturbances associated with SGA maternal and neonatal subjects, pathway enrichment and topology analyses were performed using differential metabolites identified between SGA and AGA groups. In maternal participants (AGA_M vs. SGA_M), significant enrichment was observed in taurine and hypotaurine metabolism, tryptophan metabolism, and histidine metabolism pathways (Figure 5a, left), with high pathway impact and strong statistical significance (FDR-adjusted p < 0.05). Additional alterations were noted in central carbon metabolism, citrate cycle (TCA), and arginine biosynthesis, suggesting disrupted energy and amino acid metabolism in SGA neonates. In neonatal subjects (AGA_S vs. SGA_S), tryptophan metabolism again emerged as a key altered pathway, alongside beta-alanine metabolism, arginine and proline metabolism, and pathways related to branched-chain amino acid biosynthesis (valine, leucine, and isoleucine) (Figure 5a, right). These findings suggest systemic perturbations in tryptophan metabolism in the SGA group.

Figure 5
Three panels of charts comparing metabolic pathway impacts and p-values. Panel (a) features two dot plots for “AGA_M vs SGA_M” and “AGA_S vs SGA_S” with pathway names on the y-axes and impact on the x-axes. Colors represent p-values, and dot sizes indicate count. Panel (b) shows bubble charts with “AGA_M vs SGA_M” and “AGA_S vs SGA_S”. It plots pathway impacts on the x-axes and -log10(p-values) on the y-axes. Bubble sizes and colors represent impact and -log10(p-values). Panel (c) includes two scatterplots for “AGA_M vs SGA_M” and “AGA_S vs SGA_S”. It highlights compounds with x-axes showing log2 fold changes and y-axes showing negative log10 p-values. Colors and sizes of circles convey p-value significance. Dotted borders separate panels.

Figure 5. Metabolic pathway enrichment and functional analysis reveal altered metabolic networks in SGA and AGA subjects. (a) KEGG pathway enrichment of differential metabolites between AGA and SGA groups in maternal (left) and newborn subjects (right); (b) Bubble plots for key enriched metabolic pathways in AGA and SGA groups; (c) Functional enrichment maps based on differentially abundant metabolites.

Bubble plot visualization further highlighted the significance and metabolic relevance of these pathways. In maternal subjects (Figure 5b, left), taurine/hypotaurine and tryptophan metabolism were characterized by both high impact and low p-value, while in neonatal participants (Figure 5b, right), tryptophan metabolism, along with alanine, aspartate, glutamate metabolism, and arginine and proline metabolism were prominently enriched.

Functional enrichment analysis of key differential metabolites (Figure 5c) reinforced these findings. In maternal subjects, metabolites involved in kynurenine pathway, pyrimidine metabolism, and oxidative stress response were enriched in the SGA group, including kynurenic acid, uracil, and 4-hydroxyphenylacetic acid. In neonatal participants, dysregulation was evident in galactose metabolism, bile acid biosynthesis, and linoleic acid metabolism. Notably, common enrichment across matrices in tryptophan and taurine-related pathways points to consistent metabolic stress and adaptive shifts in SGA groups. Together, these data reveal that SGA neonates exhibit broad metabolic pathway alterations spanning amino acid metabolism, energy production, and redox balance, with both shared and matrix-specific signatures across maternal and neonatal samples.

Distinct microbiota-metabolite correlation networks in AGA and SGA dyads reveal potential microbial contributions to growth programming

To investigate the potential interactions between gut microbiota and metabolic profiles, we performed Spearman correlation analyses between differential metabolites and microbial taxa in both neonates and maternal samples. In maternal subjects (Figure 6a), strong correlations were observed between specific microbial species and amino acid or bile acid-related metabolites. Notably, 1,4,5,6-Tetrahydro-6-oxonicotinic acid showed a significantly positive correlation with Desulfovibrio desulfuricans, Bifidobacterium angulatum, and uncultured Tenericutes bacterium while urocanic acid was inversely correlated with Prevotella_jejuni. Several other host-derived metabolites (e.g., cystic acid, 3-hydroxykynurenamine) displayed significant associations with Candidatus Saccharibacteria bacterium and Actinomyces sp.HMT 175, suggesting that microbial activity may contribute to the modulation of tryptophan and histidine metabolism in AGA and SGA groups. In neonatal participants (Figure 6b), broader metabolic-microbiota associations were identified. Bacteroidota such as Bacteroides caecimuris, Duncaniella dubossi, Muribaculum gordoncarteri exhibited strong positive correlations with acetyl-N-formyl-5-methoxykynurenamine, gamma-glutamylcysteine, and palmitic acid, while negative correlations with alpha-hydroxyisobutyric acid, citraconic acid, indole-3-lactic acid, and xanthurenic acid. Firmicutes such as Enterococcus faecalis was positively correlated with indole-3-lactic acid and xanthurenic acid, while negatively correlated with uric acid and quinolinic acid. Proteobacteria such as Escherichia coli showed consistent positive relationships with hydroxykynurenine. These results reveal distinct microbiota-metabolite correlation networks in neonates and their mother, supporting the hypothesis that altered maternal and infant gut microbiomes may shape metabolic phenotypes related to intrauterine growth.

Figure 6
Two clustered heatmaps display correlations between chemical compounds and microbial species. Heatmap a shows correlations using shades of red for positive and blue for negative with compounds like urocanic acid. Heatmap b uses shades of orange and green for the same, featuring compounds like acetyl-N-formyl-5-methoxykynurenamine. Each heatmap includes hierarchical clustering along the axes.

Figure 6. Spearman correlation analysis between differential metabolites and microbial species in neonates and their mother. (a) Correlation matrix in the AGA_M and SGA_M groups; (b) Correlation matrix in the AGA_S and SGA_S groups. The p-value was calculated using the two-tailed unpaired Students' t-test. *p < 0.05. Significant p-values were corrected using the Benjamini–Hochberg method. **p < 0.01, ***p < 0.001

Altered SCFA profiles reflect microbial functional shifts in mother and infants

To explore functional shifts in microbial fermentation, we quantified fecal SCFAs in both mother and infants. Quantification of fecal SCFAs revealed significantly elevated levels of butyrate in SGA_M group compared to AGA_M group (p < 0.05, Figure 7a), whereas in infants, SGA_S samples exhibited significantly lower levels of propionate and isobutyrate relative to AGA_S counterparts (p < 0.05, Figure 7b), suggesting that the observed differences in SCFA profiles may represent an early metabolic signature in SGA infants.

Figure 7
Bar charts comparing short-chain fatty acids (SCFAs) contents between groups. Chart (a) shows AGA_M and SGA_M; chart (b) shows AGA_S and SGA_S. Acetate has the highest concentration in both charts. Butyrate and isobutyrate levels show significant differences noted by asterisks.

Figure 7. Altered SCFA profiles in maternal and infant fecal samples between AGA and SGA groups. (a) Comparison of SCFA concentrations in maternal feces between AGA_M and SGA_M groups; (b) SCFA concentrations in infant feces (AGA_S vs. SGA_S). Data are presented as mean ± SD. The p-value was calculated using the two-tailed unpaired Students' t-test. *p < 0.05. The p-value interaction between groups was corrected by Bonferroni-correction.

Discussion

This study provides a multi-omics characterization of the metagenome and metabolome in SGA neonates and their mother, revealing distinct microbial, metabolic, and functional signatures compared to AGA counterparts. Our findings support the notion that SGA is associated with early-life alterations in the maternal-infant gut axis, potentially influencing neonatal development and metabolic programming.

Although alpha diversity did not differ significantly in Shannon and Simpson indices, a reduction in Chao1 richness in SGA neonates suggests a contraction in microbial species richness, consistent with previous observations linking low birth weight to reduced microbial complexity early in life (Chang et al., 2022; Chen et al., 2022). Notably, beta diversity analyses via PCoA failed to reveal significant compositional divergence, while PLS-DA demonstrated clear separability, indicating subtle but biologically relevant shifts in microbial taxa. Taxonomic profiling revealed an increased abundance of Tenericutes in SGA mother and reduced Synergistetes in SGA neonates, echoing earlier reports implicating shifts in low-abundance phyla in adverse perinatal outcomes (Li et al., 2019). Tenericutes are important constituents of the gut microbiota, whose abundance and diversity may be modulated by hormonal and immune changes during pregnancy. Alterations in Tenericutes have been associated with various pregnancy complications, including preterm birth and fetal growth abnormalities (Saraf et al., 2022).

At the species level, SGA infants exhibited enrichment in opportunistic or conditionally pathogenic taxa such as Escherichia coli (facultative anaerobe) and Enterococcus faecalis (microaerophile). The gastrointestinal tract of neonates is colonized immediately by bacteria after birth. Escherichia coli is among the first bacterial species to colonize the neonatal gut within hours after birth, establishing a mutualistic relationship with the host (Villa et al., 2013). Enterococcus faecalis is also an early colonizer of the infant gut and contributes to immune modulation (Kanasugi et al., 1996; Satonaka et al., 1996). However, the role of Enterococcus faecalis in human health remains controversial. In immune-compromised individuals, Enterococcus faecalis can act as an opportunistic pathogen (Pinholt et al., 2014). Whereas, strains isolated from healthy adults have been shown to exhibit probiotic properties (Nueno-Palop and Narbad, 2011). During early colonization, facultative anaerobes such as Escherichia coli consume residual oxygen in the intestinal lumen, thereby reducing the local redox potential and facilitating the transition to an anaerobic environment (Dominguez-Bello et al., 2011; Byndloss et al., 2018). Previous studies have demonstrated that SGA infants exhibit higher oxidative stress at birth, the severity of SGA was significantly correlated with oxidative stress levels in SGA neonates (Ashina et al., 2021). Reactive oxidative metabolites can influence the redox potential of the gut environment, thereby potentially favoring taxa with specific adaptive mechanisms. Studies have reported elevated levels of lactic acid accumulation in SGA infants at delivery and in the first few postnatal hours (Hawdon and Ward Platt, 1993). Lactic acid serves as energy source for the aerobic growth of Enterococcus faecalis (Pritchard and Wimpenny, 1978; Clarke and Knowles, 1980; London, 1968), which may partly explain its high abundance in the meconium of SGA infants. Notably, metagenomic sequencing revealed no substantial overlap between differentially abundant microbial species identified in SGA infants and those in their mother. This finding suggests a limited contribution of maternal gut microbiota to the initial microbial colonization of SGA neonates, potentially arguing against a dominant role of vertical microbial transmission in this context. In parallel with alterations in microbial taxonomic profiles, we observed significant changes in the abundance of specific CAZymes. The reduced abundance of key glycoside hydrolases (GH20, GH13_16) and glycosyltransferases (GT75, GT100, GT33) in SGA maternal subjects suggests diminished capacity for complex carbohydrate degradation and host-gut mutualism. The upregulation of GH5_9 and PL1 alongside the marked downregulation of mucin- and fucose-targeting hydrolases (GH95, GH110), glycosyltransferases (GT5, GT33), and carbohydrate-binding modules (CBM6) in SGA newborn subjects, suggests a microbiota-level metabolic shift from host-glycan-dependent colonization toward plant polysaccharide-based metabolism, potentially driven by impaired mucosal development and altered postnatal nutrient exposures. In breastfed neonates, HMOs can modulate the expression and activity of specific GH families, which in turn influences bacterial composition and HMO utilization pathways, indicating a dynamic HMO-GH-microbiota interplay (Samarra et al., 2025). This bidirectional HMO-GH-microbiota axis may help contextualize the compositional and functional differences observed in early life. Taken together, our findings suggest that impaired microbial maturation in SGA infants may reduce debranching CAZyme activity, limiting complex carbohydrate utilization and exacerbating growth restriction.

The vertical transfer of metabolome from mother to infant has been best characterized in human subjects (Lovrić et al., 2024). Tryptophan metabolism plays a critical role in regulating infant growth and development. Tryptophan is an essential aromatic amino acid primarily metabolized through the kynurenine pathway, serotonin/melatonin pathway, and indole pathway by gut microbiota. Recent studies have demonstrated that maternal-fetal tryptophan transfer remains intact in fetuses with SGA, but its conversion to quinolinic acid is impaired (Sano et al., 2016). In the present study, 3-hydroxykynurenine and xanthurenic acid were increased in SGA mother, while 3-hydroxykynurenine and quinolinic acid were decreased and xanthurenic acid was increased in SGA infants. Quinolinic acid is a key downstream metabolite of the kynurenine pathway and responsible for nicotinamide adenine dinucleotide (NAD) synthesis. Kynurenine 3-monooxygenase, an essential enzyme of the kynurenine pathway, which converts kynurenine into 3-hydroxykynurenine (Badawy, 2017). 3-Hydroxykynurenine has been shown to possess anti-inflammatory and antioxidant properties by modulating STAT1 and NF-κB signaling (Clement et al., 2021). 3-Hydroxykynurenine is degraded by kynureninase into 3-hydroxyanthranilic acid and is further degraded by kynurenine aminotransferase into xanthurenic acid. Xanthurenic acid was found to activate type-2 and−3 metabotropic glutamate receptors (mGlu2 and−3 receptors) (Fazio et al., 2017). Previous study has identified a protein in Escherichia coli with kynurenine aminotransferase activity (Han et al., 2001). This may explain the elevated levels of xanthurenic acid in the meconium of SGA infants. Acetyl-N-formyl-5-methoxykynurenamine, a metabolite of melatonin, was known to be produced or modified by the gut microbiome, and has potent antioxidant and anti-inflammatory abilities (Mayo et al., 2005; Tan et al., 2007). The level of acetyl-N-formyl-5-methoxykynurenamine was decreased in SGA mother and infants. Apart from kynurenine and melatonin metabolites, we also observed a significantly increased level of indole-3-lactic acid, another tryptophan catabolite, in SGA infants. Increased excretion of indole-3-lactic acid may reflect a shift in tryptophan metabolism toward microbial indole pathways. Previous studies have shown that Enterococcus and Escherichia can metabolize tryptophan to produce indole derivatives (include indole, indoleacetic acid, and indole-3-lactic acid) (Roager and Licht, 2018; Lee and Lee, 2010). Indole-3-lactic acid, a tryptophan-derived microbial metabolite, serves as a crucial mediator in host-microbiota crosstalk (Gao et al., 2018; Aragozzini et al., 1979; Narayanan and Rao, 1976). This increase may represent a compensatory response to impaired host kynurenine metabolism or an early microbial adaptation that supports gut immune homeostasis in the context of SGA infants. Collectively, these results indicate that SGA neonates may induce metabolic disorder by disturbing tryptophan catabolism.

In the mammalian gut, obligate anaerobe bacterium produce SCFAs such as butyrate, a key metabolite central to maintaining host-microbiota symbiosis (Jha et al., 2014). Butyrate levels increase significantly throughout pregnancy despite overall SCFA stability, potentially counteracting heightened late-gestation inflammation due to its well-established anti-inflammatory properties (Martín-Grau et al., 2022; Liu et al., 2018). Studies have shown that fecal levels of SCFAs (acetate, propionate, butyrate, isobutyrate, valerate, isovalerate, and caproate) are significantly reduced in SGA rats (An et al., 2021). In the present study, fecal SCFA profiles further corroborated microbial functional divergence, with increased butyrate levels in SGA mother and reduced propionate and isobutyrate in SGA neonates. While elevated butyrate in mother may represent a compensatory microbial adaptation to the metabolic demands of pregnancy under growth-restricted conditions, the diminished SCFA levels observed in infants may reflect delayed colonization, and differences in initial bacterial species. These findings indicate potential interruptions in the early vertical transmission of microbial functions related to SCFA production.

In summary, our results demonstrate that SGA is associated with altered gut microbiota colonization and metabolic disturbances in newborns. These findings lay the groundwork for future mechanistic studies and microbiome-targeted interventions aimed at improving health outcomes in growth-restricted infants. Several limitations warrant consideration. Despite the use of a priori power calculations, the limited sample size due to the rarity of SGA cases may still constrain our ability to detect small to moderate effect sizes, particularly given the high-dimensional nature of the omics data and the multiple comparisons performed. As a result, some biologically meaningful associations may have been overlooked. Future studies with larger cohorts or independent validation sets are needed to confirm and extend our findings. Longitudinal data are needed to evaluate the persistence and developmental consequences of these microbiota and metabolic alterations. Additionally, maternal diet, antibiotic exposure, and perinatal stress-though controlled for-may still exert residual confounding effects.

Data availability statement

The sequences of raw data generated in the present study have been deposited at the NCBI SRA database under the accession number PRJNA1277160. Untargeted metabolomics raw data and metadata have been uploaded to National Genomics Data Center under the accession number OMIX011963. Standard curve, precision, limit of quantitation, quality control frequency, and absolute concentration data for SCFAs have been provided as Supplementary Tables 1, 2.

Ethics statement

The study was approved by the Ethics Committee of the Obstetrical Department of Hangzhou Linping District Maternal & Child Health Care Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants' legal guardians/next of kin. Written informed consent was obtained from the individual(s), and minor(s)' legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.

Author contributions

XB: Conceptualization, Formal analysis, Writing – original draft. HX: Investigation, Methodology, Software, Writing – original draft. JiaL: Investigation, Visualization, Writing – original draft. JK: Resources, Supervision, Writing – original draft. FS: Software, Validation, Writing – original draft. XL: Data curation, Visualization, Writing – original draft. JinL: Conceptualization, Funding acquisition, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Medical and Health Science and Technology Plan of Department of Health of Zhejiang Province (2023XY135).

Acknowledgments

We are grateful to Suzhou PANOMIX Biomedical Tech Co., LTD. for experimental support.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1668148/full#supplementary-material

References

Alkalay, A. L., Graham, J. M. Jr., and Pomerance, J. J. (1998). Evaluation of neonates born with intrauterine growth retardation: review and practice guidelines. J. Perinatol. 18, 142–151.

PubMed Abstract | Google Scholar

An, J., Wang, J., Guo, L., Xiao, Y., Lu, W., Li, L., et al. (2021). The impact of gut microbiome on metabolic disorders during catch-up growth in small-for-gestational-age. Front. Endocrinol. 12:630526. doi: 10.3389/fendo.2021.630526

PubMed Abstract | Crossref Full Text | Google Scholar

Aragozzini, F., Ferrari, A., Pacini, N., and GualaNDRIS, R. (1979). Indole-3-lactic acid as a tryptophan metabolite produced by Bifidobacterium spp. Appl. Environ. Microbiol. 38, 544–546. doi: 10.1128/aem.38.3.544-546.1979

PubMed Abstract | Crossref Full Text | Google Scholar

Ashina, M., Kido, T., Kyono, Y., Yoshida, A., Suga, S., Nakasone, R., et al. (2021). Correlation between severity of fetal growth restriction and oxidative stress in severe small-for-gestational-age infants. Int. J. Environ. Res. Public Health 18:10726. doi: 10.3390/ijerph182010726

PubMed Abstract | Crossref Full Text | Google Scholar

Badawy, A. A. (2017). Kynurenine pathway of tryptophan metabolism: regulatory and functional aspects. Int. J. Tryptophan Res. 10:1178646917691938. doi: 10.1177/1178646917691938

PubMed Abstract | Crossref Full Text | Google Scholar

Byndloss, M. X., Pernitzsch, S. R., and Bäumler, A. J. (2018). Healthy hosts rule within: ecological forces shaping the gut microbiota. Mucosal Immunol. 11, 1299–1305. doi: 10.1038/s41385-018-0010-y

PubMed Abstract | Crossref Full Text | Google Scholar

Chang, H.-Y., Chiang Chiau, J.-S., Chang, J.-H., Hsu, C.-H., Lin, C.-Y., Ko, M. H.-J., et al. (2022). Characteristics of gut microbiota in small for gestational age infants with very low birth weight. Nutrients 14:5158. doi: 10.3390/nu14235158

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, X., Yan, Z., Liu, L., Zhang, R., Zhang, X., Peng, C., et al. (2022). Characteristics of gut microbiota of term small gestational age infants within 1 week and their relationship with neurodevelopment at 6 months. Front. Microbiol. 13:912968. doi: 10.3389/fmicb.2022.912968

PubMed Abstract | Crossref Full Text | Google Scholar

Clarke, D. J., and Knowles, C. J. (1980). The effect of haematin and catalase on Streptococcus faecalis var. zymogenes growing on glycerol. J. Gen. Microbiol. 121, 339–347. doi: 10.1099/00221287-121-2-339

PubMed Abstract | Crossref Full Text | Google Scholar

Clement, C. C., D'alessandro, A., Thangaswamy, S., Chalmers, S., Furtado, R., Spada, S., et al. (2021). 3-hydroxy-L-kynurenamine is an immunomodulatory biogenic amine. Nat. Commun. 12:4447. doi: 10.1038/s41467-021-24785-3

PubMed Abstract | Crossref Full Text | Google Scholar

Dominguez-Bello, M. G., Blaser, M. J., Ley, R. E., and Knight, R. (2011). Development of the human gastrointestinal microbiota and insights from high-throughput sequencing. Gastroenterology 140, 1713–1719. doi: 10.1053/j.gastro.2011.02.011

PubMed Abstract | Crossref Full Text | Google Scholar

Fazio, F., Lionetto, L., Curto, M., Iacovelli, L., Copeland, C. S., Neale, S. A., et al. (2017). Cinnabarinic acid and xanthurenic acid: two kynurenine metabolites that interact with metabotropic glutamate receptors. Neuropharmacology 112, 365–372. doi: 10.1016/j.neuropharm.2016.06.020

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, J., Xu, K., Liu, H., Liu, G., Bai, M., Peng, C., et al. (2018). Impact of the gut microbiota on intestinal immunity mediated by tryptophan metabolism. Front. Cell. Infect. Microbiol. 8:13. doi: 10.3389/fcimb.2018.00013

PubMed Abstract | Crossref Full Text | Google Scholar

Han, Q., Fang, J., and Li, J. (2001). Kynurenine aminotransferase and glutamine transaminase K of Escherichia coli: identity with aspartate aminotransferase. Biochem. J. 360, 617–623. doi: 10.1042/bj3600617

PubMed Abstract | Crossref Full Text | Google Scholar

Hawdon, J. M., and Ward Platt, M. P. (1993). Metabolic adaptation in small for gestational age infants. Arch. Dis. Child. 68:262. doi: 10.1136/adc.68.3_Spec_No.262

PubMed Abstract | Crossref Full Text | Google Scholar

Jha, A. K., Li Jianzheng, L. J., Yuan Yixing, Y. Y., Baral, N., and Ai Binling, A. B. (2014). A Review on Bio-butyric Acid Production and Its Optimization.

Google Scholar

Kanasugi, H., Hasegawa, T., Yamamoto, T., Abe, S., and Yamaguchi, H. (1996). Optimal dose of enterococcal preparation (FK-23) supplemented perorally for stimulation of leukocyte reconstitution in dogs treated with cyclophosphamide. J. Vet. Med. Sci. 58, 563–565. doi: 10.1292/jvms.58.563

PubMed Abstract | Crossref Full Text | Google Scholar

Lee, J.-H., and Lee, J. (2010). Indole as an intercellular signal in microbial communities. FEMS Microbiol. Rev. 34, 426–444. doi: 10.1111/j.1574-6976.2009.00204.x

PubMed Abstract | Crossref Full Text | Google Scholar

Li, H., He, Z., Gao, D., Lv, Y., Zhou, Q., Xiao, B., et al. (2019). Characteristics of the intestinal microbiota in very low birth weight infants with extrauterine growth restriction. Front. Pediatr. 7:99. doi: 10.3389/fped.2019.00099

PubMed Abstract | Crossref Full Text | Google Scholar

Li, J., Zhao, Z., Deng, Y., Li, X., Zhu, L., Wang, X., et al. (2025). Regulatory roles of quercetin in alleviating fructose-induced hepatic steatosis: targeting gut microbiota and inflammatory metabolites. Food Sci. Nutr. 13:e4612. doi: 10.1002/fsn3.4612

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, H., Wang, J., He, T., Becker, S., Zhang, G., Li, D., et al. (2018). Butyrate: a double-edged sword for health? Adv. Nutr. 9, 21–29. doi: 10.1093/advances/nmx009

PubMed Abstract | Crossref Full Text | Google Scholar

London, J. (1968). Regulation and function of lactate oxidation in Streptococcus faecium. J. Bacteriol. 95, 1380–1387. doi: 10.1128/jb.95.4.1380-1387.1968

PubMed Abstract | Crossref Full Text | Google Scholar

Lovrić, M., Horner, D., Chen, L., Brustad, N., Malby Schoos, A. M., Lasky-Su, J., et al. (2024). Vertical metabolome transfer from mother to child: an explainable machine learning method for detecting metabolomic heritability. Metabolites 14:136. doi: 10.3390/metabo14030136

PubMed Abstract | Crossref Full Text | Google Scholar

Martín-Grau, C., Díaz-López, A., Aparicio, E., and Arija, V. (2022). Short-chain fatty acid reference ranges in pregnant women from a mediterranean region of Northern Spain: ECLIPSES study. Nutrients 14:798. doi: 10.3390/nu14183798

PubMed Abstract | Crossref Full Text | Google Scholar

Mayo, J. C., Sainz, R. M., Tan, D. X., Hardeland, R., Leon, J., Rodriguez, C., et al. (2005). Anti-inflammatory actions of melatonin and its metabolites, N1-acetyl-N2-formyl-5-methoxykynuramine (AFMK) and N1-acetyl-5-methoxykynuramine (AMK), in macrophages. J. Neuroimmunol. 165, 139–149. doi: 10.1016/j.jneuroim.2005.05.002

PubMed Abstract | Crossref Full Text | Google Scholar

Narayanan, T. K., and Rao, G. R. (1976). Beta-indoleethanol and beta-indolelactic acid production by Candida species: their antibacterial and autoantibiotic action. Antimicrob. Agents Chemother. 9, 375–380. doi: 10.1128/AAC.9.3.375

PubMed Abstract | Crossref Full Text | Google Scholar

Nueno-Palop, C., and Narbad, A. (2011). Probiotic assessment of Enterococcus faecalis CP58 isolated from human gut. Int. J. Food Microbiol. 145, 390–394. doi: 10.1016/j.ijfoodmicro.2010.12.029

PubMed Abstract | Crossref Full Text | Google Scholar

Pinholt, M., Ostergaard, C., Arpi, M., Bruun, N. E., Schønheyder, H. C., Gradel, K. O., et al. (2014). Incidence, clinical characteristics and 30-day mortality of enterococcal bacteraemia in Denmark 2006-2009: a population-based cohort study. Clin. Microbiol. Infect. 20, 145–151. doi: 10.1111/1469-0691.12236

PubMed Abstract | Crossref Full Text | Google Scholar

Pritchard, G. G., and Wimpenny, J. W. (1978). Cytochrome formation, oxygen-induced proton extrusion and respiratory activity in Streptococcus faecalis var. zymogenes grown in the presence of haematin. J. Gen. Microbiol. 104, 15–22. doi: 10.1099/00221287-104-1-15

PubMed Abstract | Crossref Full Text | Google Scholar

Roager, H. M., and Licht, T. R. (2018). Microbial tryptophan catabolites in health and disease. Nat. Commun. 9:3294. doi: 10.1038/s41467-018-05470-4

PubMed Abstract | Crossref Full Text | Google Scholar

Samarra, A., Renwick, S., Arzamasov, A. A., Rodionov, D. A., Spann, K., Cabrera-Rubio, R., et al. (2025). Human milk oligosaccharide metabolism and antibiotic resistance in early gut colonizers: insights from bifidobacteria and lactobacilli in the maternal-infant microbiome. Gut Microbes 17:2501192. doi: 10.1080/19490976.2025.2501192

PubMed Abstract | Crossref Full Text | Google Scholar

Sano, M., Ferchaud-Roucher, V., Kaeffer, B., Poupeau, G., Castellano, B., and DarmaUN, D. (2016). Maternal and fetal tryptophan metabolism in gestating rats: effects of intrauterine growth restriction. Amino Acids 48, 281–290. doi: 10.1007/s00726-015-2072-4

PubMed Abstract | Crossref Full Text | Google Scholar

Saraf, V. S., Ali Zaidi, S. S., Zhu, J., Gillevet, P., Sikaroodi, M., Bokhari, H., et al. (2022). Pathogenic microbe detection in placental tissues supports placental pathobiome association with preterm birth risk in Pakistani women: a brief snapshot. Placenta 126, 160–163. doi: 10.1016/j.placenta.2022.07.002

PubMed Abstract | Crossref Full Text | Google Scholar

Satonaka, K., Ohashi, K., Nohmi, T., Yamamoto, T., Abe, S., Uchida, K., et al. (1996). Prophylactic effect of Enterococcus faecalis FK-23 preparation on experimental candidiasis in mice. Microbiol. Immunol. 40, 217–222. doi: 10.1111/j.1348-0421.1996.tb03337.x

PubMed Abstract | Crossref Full Text | Google Scholar

Shoemaker, N. B., Vlamakis, H., Hayes, K., and Salyers, A. A. (2001). Evidence for extensive resistance gene transfer among Bacteroides spp. and among Bacteroides and other genera in the human colon. Appl. Environ. Microbiol. 67, 561–568. doi: 10.1128/AEM.67.2.561-568.2001

PubMed Abstract | Crossref Full Text | Google Scholar

Tan, D. X., Manchester, L. C., Terron, M. P., Flores, L. J., and Reiter, R. J. (2007). One molecule, many derivatives: a never-ending interaction of melatonin with reactive oxygen and nitrogen species? J. Pineal Res. 42, 28–42. doi: 10.1111/j.1600-079X.2006.00407.x

PubMed Abstract | Crossref Full Text | Google Scholar

United Nations International Children's Emergency Fund (2023). Low Birthweight. Available online at: https://data.unicef.org/topic/nutrition/low-birthweight/ (Accessed January 20, 2025).

Google Scholar

Usher, R., and McLean, F. (1969). Intrauterine growth of live-born Caucasian infants at sea level: standards obtained from measurements in 7 dimensions of infants born between 25 and 44 weeks of gestation. J. Pediatr. 74, 901–910. doi: 10.1016/S0022-3476(69)80224-6

PubMed Abstract | Crossref Full Text | Google Scholar

Vatanen, T., Jabbar, K. S., Ruohtula, T., Honkanen, J., Avila-Pacheco, J., Siljander, H., et al. (2022). Mobile genetic elements from the maternal microbiome shape infant gut microbial assembly and metabolism. Cell 185 4921-4936.e15. doi: 10.1016/j.cell.2022.11.023

PubMed Abstract | Crossref Full Text | Google Scholar

Villa, R., Martorana, A. M., Okuda, S., Gourlay, L. J., Nardini, M., Sperandeo, P., et al. (2013). The Escherichia coli Lpt transenvelope protein complex for lipopolysaccharide export is assembled via conserved structurally homologous domains. J. Bacteriol. 195, 1100–1108. doi: 10.1128/JB.02057-12

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, X., Liu, J., Hui, X., and Song, Y. (2022). Metabolomics applied to cord serum in preeclampsia newborns: implications for neonatal outcomes. Front. Pediatr. 10:869381. doi: 10.3389/fped.2022.869381

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, X., Feng, Q., Zhang, T., Gao, Y., Cheng, Q., Zhang, W., et al. (2024). Infant age inversely correlates with gut carriage of resistance genes, reflecting modifications in microbial carbohydrate metabolism during early life. iMeta 3:e169. doi: 10.1002/imt2.169

PubMed Abstract | Crossref Full Text | Google Scholar

Zelena, E., Dunn, W. B., Broadhurst, D., Francis-Mcintyre, S., Carroll, K. M., Begley, P., et al. (2009). Development of a robust and repeatable UPLC-MS method for the long-term metabolomic study of human serum. Anal. Chem. 81, 1357–1364. doi: 10.1021/ac8019366

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: small for gestational age infants, metagenome, metabolome, tryptophanmetabolism, short-chain fatty acid

Citation: Bian X, Xu H, Li J, Kuang J, Shi F, Li X and Li J (2025) Integrated multi-omics reveals distinct maternal and neonatal gut microbial and metabolic signatures associated with small for gestational age status. Front. Microbiol. 16:1668148. doi: 10.3389/fmicb.2025.1668148

Received: 18 July 2025; Accepted: 31 October 2025;
Published: 26 November 2025.

Edited by:

Mingsong Kang, Canadian Food Inspection Agency (CFIA), Canada

Reviewed by:

Anne Christine Aust, University of Giessen, Germany
Bedour Abdelghany, Port Said University, Egypt
Kristy Thomas, Medical University of South Carolina, United States

Copyright © 2025 Bian, Xu, Li, Kuang, Shi, Li and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jinjun Li, bGlqaW5qdW5AemFhcy5hYy5jbg==

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.