Warning: fopen(/home/virtual/pediatrics/journal/upload/ip_log/ip_log_2025-10.txt) [function.fopen]: failed to open stream: Permission denied in /home/virtual/pediatrics/journal/ip_info/view_data.php on line 93

Warning: fwrite(): supplied argument is not a valid stream resource in /home/virtual/pediatrics/journal/ip_info/view_data.php on line 94
Dual-strain probiotics <i xmlns="">Bifidobacterium bifidum</i> and <i xmlns="">Lactobacillus acidophilus</i> reverse gut dysbiosis in preterm neonates: a randomized controlled trial

Dual-strain probiotics Bifidobacterium bifidum and Lactobacillus acidophilus reverse gut dysbiosis in preterm neonates: a randomized controlled trial

Article information

Clin Exp Pediatr. 2025;68(10):763-771
Publication date (electronic) : 2025 August 6
doi : https://doi.org/10.3345/cep.2025.00374
1Department of Pediatrics, Faculty of Medicine Ramathibodi Hospital, Mahidol University, Bangkok, Thailand
2Chakri Naruebodindra Medical Institute, Faculty of Medicine Ramathibodi Hospital, Mahidol University, Bangkok, Thailand
Corresponding author: Pornthep Tanpowpong, MD, MPH. Department of Pediatrics, Faculty of Medicine Ramathibodi Hospital, Mahidol University, Bangkok, Thailand Email: pornthep.tan@mahidol.ac.th
Received 2025 February 12; Revised 2025 June 6; Accepted 2025 June 7.

Abstract

Background

Preterm neonates exhibit gut dysbiosis, characterized by increased numbers of pathogenic bacteria and decreased Bifidobacterium and Lactobacillus levels. Supplementation with the probiotic Bifidobacterium bifidum/Lactobacillus acidophilus (BB/LA) may reverse gut dysbiosis.

Purpose

To study the effects of BB/LA on the gut microbiota of preterm neonates.

Methods

We enrolled neonates born between July 2022 and September 2023 with a gestational age of <33 weeks or birth weight of <1,500 g. After randomization into probiotic (PG) and control (CG) groups, stool samples were collected at 3 time points: birth (V1), 35 weeks' postmenstrual age (PMA) (V2), and 4 months of age (V3). BB/LA was administered to the PG until PMA 35 weeks. All neonates received a feeding protocol similar to that of predominant breast milk. Stool samples were stored at -80°C, a DNA extraction performed, and 16S rRNA gene sequencing used to define alpha and beta diversities and the relative abundances of the bacteria. Baseline characteristics and clinical outcomes were collected.

Results

We analyzed 68 neonates (33 in the PG, 35 in the CG). The alpha diversities did not differ significantly between the groups at any time point. At V1, beta diversity was not significantly different between the 2 groups. After BB/LA supplementation (V2), beta diversity was significantly greater in the PG versus CG (P=0.004). The relative abundances of Bifidobacterium and Lactobacillus were higher in the PG (both P<0.001), whereas that of Clostridium senso stricto 1 was higher in the CG (P=0.017). Growth parameters, necrotizing enterocolitis, and mortality rate did not differ between groups. No adverse events were observed.

Conclusion

BB/LA led to healthier gut microbiota in preterm neonates as demonstrated by a reversal of gut dysbiosis characterized by increased beta diversity, increased the relative abundances of Bifidobacterium and Lactobacillus, and decreased the relative abundance of Clostridium senso stricto 1.

Key message

Question: Can probiotics BB/LA reverse gut dysbiosis in preterm neonates?

Finding: BB/LA supplementation induced more diverse beta diversity and increased relative abundances of Bifidobacterium, Lactobacillus and decreased relative abundance Clostridium.

Meaning: Early BB/LA supplementation could reverse gut dysbiosis in preterm neonates.

Graphical abstract

Introduction

Preterm neonates often exhibit gut dysbiosis among other gastrointestinal problems. Gut dysbiosis in preterm neonates is caused by many disadvantageous factors such as Cesarean delivery, nonbreast feeding, exposure to broad spectrum antibiotics and prolonged hospital stay with multiple harsh environmental exposures in the neonatal intensive care unit (NICU) [1-3]. Consequently, gut microbiota of preterm neonates typically has aberrant microbial colonization (i.e., reduced diversity) with notable numbers of pathogenic bacteria such as the Enterobacteriaceae which are Gram-negative facultative anaerobes (i.e., Escherichia coli, Klebsiella pneumoniae) and a lower abundance of commensal obligate anaerobes such as Bifidobacterium, Lactobacillus, and Bacteroides [4], the aforementioned microbiota condition can be defined as gut dysbiosis [5]. Intestinal dysbiosis in preterm neonates has been linked to necrotizing enterocolitis (NEC) [6-8]. Thus, preterm neonates are one of the most vulnerable populations to encounter detrimental gastrointestinal problems, most notably NEC [1]. Fecal microbiota analysis of preterm neonates with NEC demonstrated an increased abundance of Enterobacteriaceae and Clostridiaceae and a decrease of Bifidobacteriaceae, Firmicutes, and Bacteroidetes [1,3].. Besides gut dysbiosis, impaired maturation of gastrointestinal tract as well as the underdevelopment of the gut immune system make it difficult for the preterm neonates to tolerate optimal enteral feeding and undergo optimal and sustainable growth [2,9].

In contrast to gut dysbiosis, healthy gut microbial colonization in neonates also known as gut eubiosis, has been shown as a pivotal role in maintaining balanced intestinal immunity, nutrient digestion and absorption, intestinal motility, and energy utilization [1]. Intestinal eubiosis is typically observed in term, vaginally delivered, exclusively breastfed infants [3,5,10]. The proposed healthy gut microbiota is characterized by an abundance in Gram-positive Firmicutes such as Staphylococcus, Bifidobacterium, and Lactobacillus with the ability to produce glycosidases necessary to digest human milk oligosaccharides (HMOs) [3]. Probiotics supplementation in the early life of preterm neonates might reverse or modify the gut dysbiosis to resemble those of term neonates with gut eubiosis.

Probiotics are live microorganisms that, when administered in adequate amounts, confer a healthy benefit to the host [11]. Early supplementation of probiotics, either multistrain or single-strain, to preterm neonates aims to alter gut microbiota environment by colonizing the immature gut with commensal bacteria which are naturally present in the presumably eubiotic environment of exclusively breastfed term newborns [9,12]. When the gastrointestinal tract is colonized by these bacteria derived from the probiotics, studies have shown that the preterm neonates are less likely to develop late-onset sepsis and NEC [1,3,13]. However, the beneficial mechanisms of probiotics remain largely unknown but are thought to involve changes in intestinal permeability, enhanced mucosal IgA responses to foreign antigens and increased production of anti-inflammatory cytokines [14]. Promising beneficial effects of probiotics have been shown in various clinical trials in the prevention of NEC [14-19], late-onset sepsis [15,16], and reduction in mortality [14,15,17,19]. However, many other previous studies in this area demonstrated inconsistent results, in which some studies did not show clinical benefits of probiotics supplementation in preterm neonates. Two meta-analyses concluded that only specific strains of probiotics such as various strains of Bifidobacterium spp. yielded promising clinical benefits while others did not [16,18]. The 2020 European guideline for using probiotics in preterm neonates only provides a conditional recommendation on administering specific strains of probiotics in the prevention of NEC [20].

A dual-strain probiotics consisting of Bifidobacterium bifidum/Lactobacillus acidophillus (BB/LA) has been shown in neonates to increase the abundance of Bifidobacterium spp. and Lactobacillus spp. with a reduction of Enterobacteriaceae, compared with placebo/the control group [12,21]. The use of probiotics BB/LA has increased among preterm newborns in many NICUs across Europe with good safety profiles [22,23] and has been shown to significantly reduce the risk of NEC, nosocomial bacteremia, sepsis, and death [24-26]. It is well known that breast milk feeding reduces NEC risk in preterm neonates [27,28]. Probiotics has a synergistic effect among the breastfed neonates in altering the gut microbiota [21,29] and preventing NEC [18,23]. A reasonable explanation would be that breast milk contains prebiotic HMOs which act as an important metabolic source of growth for the Bifidobacterium [11,25]. Nowadays many NICUs encourage breast milk as a primary source of enteral feeding in the preterm neonates, we believe that BB/LA would yield advantageous gut microbiota profile for the predominantly breastfed preterm neonates.

Probiotics BB/LA administration in the neonatal period has been shown to maintain the colonization of Bifidobacterium after discontinuing it during the infancy period [21], indicating the sustainable beneficial effects. According to previously mentioned studies, we hypothesized that preterm neonates generally exhibit gut dysbiosis earlier in life, which is characterized by high levels of Enterobacteriaceae and Clostridiacieae and relatively low levels of Bifidobacterium and Lactobacillus. With the aid of probiotics BB/LA, the gut dysbiosis might be reversed or attenuated. Since there are only a handful of randomized controlled studies on the prevention of gut dysbiosis by BB/LA, therefore we aim to compare gut microbial profiling of the preterm infants who receive in-hospital probiotics BB/LA and those who did not.

Methods

We enrolled preterm newborns who were born between July 2022 and September 2023. This clinical trial was reviewed and approved by the Institutional Review Board of Faculty of Medicine Ramathibodi Hospital (MURA 2022/153) and registered in the Thai Clinical Trials Registry (TCTR20220829005). This study was carried out with strict adherence to the World Medical Association’s Declaration of Helsinki. Written informed consent was obtained from the parents or guardians of all participants. The inclusion criteria were as follows: preterm newborns with a gestational age (GA) <33 weeks or birth weight <1,500 grams (i.e., very low birth weight). The exclusion criteria were cyanotic congenital heart disease, chromosomal or gastrointestinal tract anomalies, or enteral feeding was withheld for more than 14 consecutive days. The enrolled infants were randomized into the probiotics group (PG) and control group (CG).

Preterm newborns in PG received BB/LA (Infloran, Laboratorio Farmaceutico SIT Mede, Italy). One capsule of BB/LA (250 mg) contains 1 billion colony forming unit (CFU) lyophilized Bifidobacterium bifidum and 1 billion CFU lyophilized Lactobacillus acidophilus. After having been fed breast milk for one day without complications, one capsule of BB/LA was administered once daily via dissolving in 1–2 mL of sterile water prepared by neonatal registered nurses, then tube-fed before feeding milk. BB/LA preparations were freshly prepared right before administering to the infants. Throughout the trial, BB/LA was kept in a 4°C temperature-controlled refrigerator, separated from the milk preparation station to avoid cross-contamination with CG. BB/LA was administered until postmenstrual age (PMA) of 35 weeks and continued until the infants were discharged. In CG, besides a standard neonatal care similar to the PG, preterm newborns did not receive BB/LA during the study. The feeding protocol strongly supported the use of breast milk as the first choice in all studied infants. However, when there was insufficient breast milk in some infants, pasteurized donor human milk (DHM) would be given after written informed consent was obtained from the parents. Preterm newborns in both groups received similar standardized feeding protocol and medical care without interference from the researchers.

Perinatal and maternal data was obtained by interviewing the mother according to the case record form and from electronic medical records. All preterm newborns would receive empirical antibiotics if there was clinical suspicion of neonatal sepsis. Then, if blood cultures were negative at 48 to 72 hours and the newborns had no signs of sepsis, the antibiotics were discontinued. Enteral feeding was initiated for studied newborns within the first few days of life with maternal breast milk (MBM) or DHM in case of unavailable mother’s breast milk. MBM or DHM was continuously given at least 2 weeks after birth but not more than 4 weeks. After 2 weeks, if there was no MBM, premature formula milk was gradually introduced, alternating with MBM/DHM. At any point, if MBM became available, the infants would receive MBM as the first choice. Full enteral feeding was defined as a total volume of 150 mL/kg/day.

The focused clinical outcomes associated with gut dysbiosis were NEC with modified Bell’s stage ≥2A, late-onset neonatal sepsis, and death. Time to full enteral feeding, duration of hospitalization, body weight at discharge, 2 months and 4 months, respiratory distress syndrome and bronchopulmonary dysplasia were also recorded. Factors altering gut microbiota of preterm newborns such as duration of antibiotics administration, route of delivery and types of feeding were also recorded.

1. Collection of stool samples

Nurses in the NICU collected fresh stool samples from the diapers of the infants into a sterile container at 3 different time points. Stool samples were collected within 72 hours after birth as meconium (V1) and after PMA 35 weeks (V2). After being discharged from the hospital, infants in both groups did not further receive BB/LA. At chronological age of 4 months, the infants were scheduled for follow up and vaccination at our outpatient clinic. At the clinic, stool samples were collected as the third sample (V3) for postdischarge analysis.

2. DNA extraction

All stabilized samples were then transferred to the Center for Medical Genomics, Faculty of Medicine Ramathibodi Hospital, where they were stored at -80°C until processing. Metagenomic DNA was extracted using QIAamp PowerFecal DNA Kit (QIAGEN, Germany) according to the manufacturer's protocol. DNA quality and quantity were measured using NanoDrop One spectrophotometer and Qubit 3.0 fluorometer (Thermo Fisher Scientific, USA), with acceptable quality thresholds of A260/280 ratio between 1.8-2.0 and A260/230 ratio above 1.4.

3. Polymerase chain reaction amplification and DNA sequencing

The V3–V4 hypervariable regions of bacterial 16s ribosomal RNA genes were amplified using a 2-step polymerase chain reaction (PCR) approach. The first PCR amplification was performed using primers 341F (5′-CCTAYGGGRBGCASCAG-3′) and 806R (5′-GGACTACHVGGGTWTCT AAT-3′). Each PCR reaction mixture (total volume 15 μL) contained Phusion High-Fidelity PCR Master Mix (New England Biolabs, USA), 0.2 μM of each primer, and approximately 10 ng of template DNA. The PCR conditions were as follows: Initial denaturation at 98°C for 1 minute; 30 cycles of denaturation at 98°C for 10 seconds, annealing at 50°C for 30 seconds, and extension at 72°C for 30 seconds; followed by a final extension at 72°C for 5 minutes. Index sequences were subsequently added in the second PCR step using NEBNext Ultra II FS DNA PCR-free Library Prep Kit (New England Biolabs). The prepared libraries underwent paired-end sequencing (2×250 bp) using NovaSeq 6000 SP Reagent Kit v1.5 (500 cycles) on the Illumina NovaSeq 6000 platform.

4. Bioinformatics analysis of gut microbiota composition

Raw sequencing data underwent quality control and preprocessing using FLASH software before being analyzed with Quantitative Insights into Microbial Ecology 2 (QIIME 2). The sequences were denoised using DADA2 to generate amplicon sequence variants (ASVs), which were then taxonomically classified using the SILVA 138 ribosomal RNA database. ASV count tables were generated and grouped at various taxonomic levels.

Microbial community diversity was assessed through both alpha and beta diversity analyses. Alpha diversity, which reflects the complexity of species diversity within samples, was evaluated using 4 indices: abundance-based coverage estimator, Chao1 (species richness estimation), Shannon index (diversity considering both richness and evenness), and Simpson index (community dominance). For beta diversity analysis, principal coordinate analysis (PCoA) with Bray-Curtis dissimilarity matrices was performed to compare microbial community compositions between the 2 study groups, as well as across different visiting timepoints. Differential abundance analysis was conducted using linear discriminant analysis effect size (LEfSe) when differences in taxonomic abundance were statistically different and then visualized through taxonomic bar charts. All diversity metrics and LEfSe analyses were performed using MicrobiomeAnalyst 2.0 (https://www.microbiomeanalyst.ca/). The data were then subjected to comprehensive statistical analyses to identify significant associations and differences between groups.

QIIME generates operational taxonomic unit tables containing information regarding the microbiome composition and relative abundances of each bacterial taxon (from phylum to genus) in each sample. QIIME was used to assign taxonomy and determine alpha and beta diversities. To measure alpha diversity, 4 parameters were calculated; Chao 1, Shannon, and Simpson indices. To assess beta diversity, PCoA of Bray-Curtis dissimilarity was performed to depict the dissimilarity of bacterial populations between the PG and CG. If differences in taxonomic abundance were statistically different between the groups, LEfSe was performed.

5. Biostatistical analyses

We conducted statistical analysis using STATA/SE program ver. 18.0 (StataCorp, USA). For comparison of categorical data, we used Fisher exact test and chi-square test. The normal distribution of continuous data was evaluated using the Shapiro-Wilk test. For comparison of continuous data with normal distribution, we used Student t test. For continuous data with nonnormal distribution, we used quantile regression.

For microbiome statistical analyses, alpha diversity comparisons between groups were performed using univariate and multivariate linear regression analyses. For nonparametric comparisons, Wilcoxon rank-sum tests were applied with false discovery rate control using the Benjamini-Hochberg method. P values for beta diversity were calculated using permutational multivariate analysis of variance. Statistical significance was set at P value <0.05.

Results

1. Perinatal and maternal characteristics

In our study, 93 preterm newborns were initially enrolled and 25 of them were excluded; 2 had cyanotic congenital heart disease, 2 had gastrointestinal tract anomalies, 9 died before randomization, 2 referred out, and 10 of their parents refused to participate. We had 68 patients remained in the trial (33 in PG and 35 in CG). Upon the process of fecal DNA extraction and PCR amplification, some samples were unqualified for sequencing, thus were excluded from analysis (Fig. 1). Baseline data including maternal age, mode of delivery, sex, GA, birth weight, Apgar scores and perinatal antibiotics exposure were comparable in both groups (Table 1). All infants were given mixed feeding of 50% fortified breast milk and 50% premature formula when the second stool samples (after PMA 35 weeks) were collected.

Fig. 1.

Flow diagram of the study. GA, gestational age; BW, body weight; PMA, postmenstrual age.

Maternal and perinatal characteristics of the preterm newborns

2. Gut microbiota profile

The alpha diversity of meconium samples (V1) between PG and CG were not statistically different (Supplementary Fig. 1). Also, the beta diversities of meconium samples represented by Bray-Curtis index were comparable (P=0.323) (Supplementary Fig. 2). Taxonomy abundance of meconium samples in both groups revealed predominance of Staphylococcus as shown in dark green color in Fig. 2A. The relative abundance of Staphylococcus in infants delivered via Cesarean section was significantly higher than those who were vaginally delivered (P=0.04) (Supplementary Fig. 3).

Fig. 2.

Taxonomy abundances of stool samples across 3 time points (V1–V3). (A) At V1 (birth), an overall predominance of Staphylococcus was seen in both groups. (B) At V2 (35 weeks' postmenstrual age), several peaks of Bifidobacterium were seen in the probiotics group along with a predominance of Klebsiella in the control group. (C) At V3 (4 months of age), an overall predominance of Bifidobacterium was noted in both groups.

The alpha diversity of V2 was also not statistically different between the 2 groups (Supplementary Fig. 4). Whereas, the beta diversity of V2 in the PG was significantly diverse when compared with the CG group (P=0.004) (Fig. 3). The taxonomy abundance of V2 showed a more predominance of Bifidobacterium in the PG as shown in dark pink color in Fig. 2B. LefSe analysis depicted several differences in the relative abundance of certain bacteria including higher abundances of Bifidobacterium spp. (P<0.001) (Fig. 4A) and Lactobacillus spp. in the PG (P<0.001) (Fig. 4B), and a higher abundance of Clostridium senso stricto 1 in the CG (P=0.017) (Fig. 4C). The relative abundances of Klebsiella, Escherichia-Shigella, Enterococcus, Enterobacter, Staphylococcus, Streptococcus, and Serratia were not statistically different between the 2 groups.

Fig. 3.

Beta diversity at time point V2 (35 weeks' postmenstrual age). PERMANOVA, permutational multivariate analysis of variance; CG, control group; PG, probiotics group.

Fig. 4.

Relative abundances of Bifidobacterium (A), Lactobacillus (B), and Clostridium senso stricto 1 (C) at V2 (35 weeks' postmenstrual age).

The alpha and beta diversities at V3 were not statistically different between the PG and CG (Supplementary Figs. 5 and 6, respectively). Taxonomy abundance of V3 showed a predominance of Bifidobacterium spp. in both groups as shown in dark green color (Fig. 2C).

Supplementary Fig. 7 demonstrates a higher Simpson index at V2 in both groups but the alpha diversity defined by this index seemed to decrease more in the PG vs. CG group in V3 without statistical significance (P=0.12).

3. Clinical outcomes

Most clinical outcomes including early-onset and late-onset neonatal sepsis, NEC stage ≥2A and lung complications were comparable as shown in Table 2. The number of days to full enteral feeding was approximately 10 days in both groups with comparable length of hospitalization of 44 days. There were 2 deaths in the PG; one died without known cause at home and the other died from severe congenital brain anomaly, which was diagnosed after randomization. One infant in the CG died due to spontaneous intestinal perforation. Also, the body weight at 2 months and 4 months between the PG and CG were not statistically different.

Clinical outcomes

Discussion

The gut microbiota profiles of meconium samples in preterm newborns during the first few days of life in our study were predominantly colonized with Staphylococcus, which is consistent with previous studies [30-32]. This can be explained by the fact that most of the preterm infants in this study were delivered via Cesarean delivery, thus making the gut microbiome of their early life resembled those of the skin surface of their mothers [3,33,34]. Therefore, the overall profile of gut microbiome in the first days of life was therefore not diverse. In accordance with our hypothesis, we found significant difference in the gut microbiota profile after we administered probiotics BB/LA to the newborns. The gut microbiota demonstrated a more eubiotic environment [3,5,10], with significantly higher abundances of Bifidobacterium and Lactobacillus [31]. Whereas in the CG, the relative abundance of Clostridium senso stricto 1 was significantly higher. In previous studies, Clostridium senso stricto 1 was more abundant preceding the development of NEC in the preterm newborns [35]. Another signal of gut dysbiosis was that Bifidobacterium and Lactobacillus were virtually non-existent in CG at PMA 35 weeks. Lower numbers of Bifidobacteriaceae and higher numbers of Clostridiaceae in the preterm newborns could be a signal of gut dysbiosis [3], which was also associated with NEC [36]. Although not statistically different, probably due to a relatively small study population, preterm newborns in the CG had a slightly higher incidence of NEC stage ≥2A (12% vs. 3% in the PG, Table 2). Our research findings are consistent with the previous studies that showed the ability of probiotics to modulate gut microbiota of preterm newborns into a healthier environment [12,21,29,30,37].

We did not find any differences in alpha diversities between the 2 groups at all 3 time points, which are consistent with a previous study by Chang et al. [30]. Alpha diversity measures the richness of bacteria in each sample, hence at birth, the considerably low alpha diversities in both groups were not unexpected because of similar perinatal characteristics including route of delivery (mostly Cesarean delivery), one of the crucial factors contributing to gut dysbiosis early after birth [1,2]. Despite not statistically different after PMA 35 weeks, the alpha diversity in both groups increased overtime which is consistent with the study by Korpela et al. [32], where the diversity of gut microbiota in preterm newborns increased along with postnatal age. The cause of increased alpha diversity might be different between the neonates in both groups. In the PG, their gut microbiota profile was largely modulated by the BB/LA resulting in the higher abundance of Bifidobacterium and Lactobacillus. Infants in the CG were also exposed to the harsh environment such as NICU setting, antibiotics administration, and formula feeding, thus potentially resulting in a higher number of Clostridium senso stricto 1 [35]. Beta diversity at PMA 35 weeks demonstrated a clear dissimilarity of gut microbiota between the 2 groups, in which PG had a more diverse profile consistent with the study by Chang et al. [30].

From discharged until 4 months of age, infants in both groups relocated to their home environment for several weeks, which likely altered their gut microbiota profiles. We found no difference of alpha and beta diversities between the 2 groups. This might be explained by the fact that they were not exposed to hospital environment for several weeks and most of them received similar type of feeding. Probiotics BB/LA could not maintain gut microbiota diversity after its discontinuation. Surprisingly, gut microbiota in both groups showed a predominance of Bifidobacterium despite discontinuing BB/LA after discharged. In addition, probiotics BB/LA might pertain to the colonization of Bifidobacterium in the gut of former preterm infants [21] until 4 months of age after discontinuation. In our research, most patients after discharged receive half of their intake as breast milk and the other half as formula milk. This might explain that the predominance of Bifidobacterium was enhanced by breast milk that contained HMOs which serve as a prebiotic for Bifidobacterium [11,25]. Unfortunately, we were not able to perform the LefSe analysis to depict a statistical difference of Bifidobacterium between the 2 groups because of the limited number of stool samples.

There were several advantages in this clinical trial. We compared the gut microbiota profiles of preterm newborns in 3 different time points; at birth, at PMA 35 weeks and at age 4 months between the infants who received probiotics BB/LA and infants who did not. This comparison elucidated the characteristics of gut microbiota of preterm newborns which exhibited gut dysbiosis. Our NICU has a unified feeding protocol where all preterm newborns mostly receive breast milk or DHM in early weeks after birth. This research also showed that continuing breast milk after discharge from the hospital helped sustaining the abundance of Bifidobacterium. This result again stresses the importance of breast milk in preterm newborns in improving the chance of maintaining gut eubiosis. However, there were some limitations in this study. The number of infants and stool samples for gut microbiota analysis were relatively small. Therefore, we were not able to demonstrate significant differences of gut microbiome at 4 months of follow-up and did not find clinical significances on the use of probiotics BB/LA such as the prevention of NEC and late-onset neonatal sepsis. Our technique of 16s rRNA sequencing only detected the bacteria at the genus level, so we could not tell whether the Bifidobacterium and Lactobacillus found in the stool samples belonged to the probiotic’s species (BB/LA) or not.

In conclusion, this research showed that probiotics BB/LA modulated gut microbiota of preterm newborns into a more eubiotic environment with diverse commensal bacteria such as Bifidobacterium and Lactobacillus and low pathogenic bacteria such as Clostridium senso stricto 1. Studying the metabolomic profile of preterm gut microbiota affected by probiotics may shed the lights on mechanism of action of probiotics and further explain their benefits on clinical outcome.

Supplementary materials

Supplementary Figs. 1-7 are available at https://doi.org/10.3345/cep.2025.00374.

Supplementary Fig. 1.

Alpha diversity at V1 (birth). (A) Abundance-based coverage estimator (ACE) index. (B) Chao1 index. (C) Shannon index. (D) Simpson index. No significant intergroup differences were noted in any parameters.

cep-2025-00374-Supplementary-Fig-1.pdf
Supplementary Fig. 2.

Beta diversity at V1 (birth). PERMANOVA, permutational multivariate analysis of variance; CG, control group; PG, probiotics group.

cep-2025-00374-Supplementary-Fig-2.pdf
Supplementary Fig. 3.

Relative abundance of Staphylococcus after vaginal versus cesarean delivery. At V1 (birth), cesarean-delivered infants showed a significantly higher relative abundance of Staphylococcus than those who were vaginally delivered.

cep-2025-00374-Supplementary-Fig-3.pdf
Supplementary Fig. 4.

Alpha diversity at V2 (35 weeks' postmenstrual age). (A) Abundancebased coverage estimator (ACE) index. (B) Chao1 index. (C) Shannon index. (D) Simpson index. No significant intergroup differences were noted in any parameters.

cep-2025-00374-Supplementary-Fig-4.pdf
Supplementary Fig. 5.

Alpha diversity at V3 (4 months of age). (A) Abundance-based coverage estimator (ACE) index. (B) Chao1 index. (C) Shannon index. (D) Simpson index. No significant intergroup differences were noted in any parameters.

cep-2025-00374-Supplementary-Fig-5.pdf
Supplementary Fig. 6.

Beta diversity at V3 (4 months of age). PERMANOVA, permutational multivariate analysis of variance; CG, control group; PG, probiotics group.

cep-2025-00374-Supplementary-Fig-6.pdf
Supplementary Fig. 7.

Simpson index alpha diversity at V1–V3 by study group. The Simpson index was low at birth (V1) but increased after 35 weeks' postmenstrual age (V2) in both groups. At 4 months of age (V3), the Simpson index was decreased more in the probiotics versus control group, but the difference was insignificant. The dashed line represents the mean value.

cep-2025-00374-Supplementary-Fig-7.pdf

Notes

Conflicts of interest

No potential conflict of interest relevant to this article was reported.

Funding

This study received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Acknowledgments

This study was financially supported by a grant from MU’s Strategic Research Fund from Mahidol University. PT received a mid-career Career Development Grant from Faculty of Medicine Ramathibodi Hospital.

Author Contribution

Conceptualization: SS, PT, PP, PN; Data curation: SS; Formal analysis: SS; Funding acquisition: PT ; Methodology: SS, PT, PP, PN; Project administration: SS; Visualization: SS; Writing-original draft: SS; Writing-review & editing: PT, PP, PN

References

1. Baldassarre ME, Di Mauro A, Capozza M, Rizzo V, Schettini F, Panza R, et al. Dysbiosis and prematurity: is there a role for probiotics? Nutrients 2019;11:1273.
2. Xu L, Wang Y, Wang Y, Fu J, Sun M, Mao Z, et al. A double-blinded randomized trial on growth and feeding tolerance with Saccharomyces boulardii CNCM I-745 in formula-fed preterm infants. J Pediatr (Rio J) 2016;92:296–301.
3. Underwood MA, Mukhopadhyay S, Lakshminrusimha S, Bevins CL. Neonatal intestinal dysbiosis. J Perinatol 2020;40:1597–608.
4. Navarro-Tapia E, Sebastiani G, Sailer S, Almeida Toledano L, Serra-Delgado M, García-Algar Ó, et al. Probiotic supplementation during the perinatal and infant period: effects on gut dysbiosis and disease. Nutrients 2020;12:2243.
5. Henderickx JG, Zwittink RD, van Lingen RA, Knol J, Belzer C. The preterm gut microbiota: an inconspicuous challenge in nutritional neonatal care. Front Cell Infect Microbiol 2019;9:85.
6. Zhao W, Ho HE, Bunyavanich S. The gut microbiome in food allergy. Ann Allergy Asthma Immunol 2019;122:276–82.
7. Lee KH, Song Y, Wu W, Yu K, Zhang G. The gut microbiota, environmental factors, and links to the development of food allergy. Clin Mol Allergy 2020;18:5.
8. Matsui S, Kataoka H, Tanaka JI, Kikuchi M, Fukamachi H, Morisaki H, et al. Dysregulation of intestinal microbiota elicited by food allergy induces IgA-mediated oral dysbiosis. Infect Immun 2019;88:e00741-19.
9. Indrio F, Riezzo G, Tafuri S, Ficarella M, Carlucci B, Bisceglia M, et al. Probiotic supplementation in preterm: feeding intolerance and hospital cost. Nutrients 2017;9:965.
10. Arboleya S, Sanchez B, Milani C, Duranti S, Solis G, Fernandez N, et al. Intestinal microbiota development in preterm neonates and effect of perinatal antibiotics. J Pediatr 2015;166:538–44.
11. Poindexter B, Committee On F. Use of probiotics in preterm infants. Pediatrics 2021;147e2021051485.
12. Esaiassen E, Hjerde E, Cavanagh JP, Pedersen T, Andresen JH, Rettedal SI, et al. Effects of probiotic supplementation on the gut microbiota and antibiotic resistome development in preterm infants. Front Pediatr 2018;6:347.
13. Meyer MP, Chow SS, Alsweiler J, Bourchier D, Broadbent R, Knight D, et al. Probiotics for prevention of severe necrotizing enterocolitis: experience of New Zealand neonatal intensive care units. Front Pediatr 2020;8:119.
14. Olsen R, Greisen G, Schroder M, Brok J. Prophylactic probiotics for preterm infants: a systematic review and meta-analysis of observational studies. Neonatology 2016;109:105–12.
15. Sun J, Marwah G, Westgarth M, Buys N, Ellwood D, Gray PH. Effects of probiotics on necrotizing enterocolitis, sepsis, intraventricular hemorrhage, mortality, length of hospital stay, and weight gain in very preterm infants: a meta-analysis. Adv Nutr 2017;8:749–63.
16. Aceti A, Gori D, Barone G, Callegari ML, Di Mauro A, Fantini MP, et al. Probiotics for prevention of necrotizing enterocolitis in preterm infants: systematic review and meta-analysis. Ital J Pediatr 2015;41:89.
17. AlFaleh K, Anabrees J. Probiotics for prevention of necrotizing enterocolitis in preterm infants. Cochrane Database Syst Rev 2014;(4):CD005496.
18. Thomas JP, Raine T, Reddy S, Belteki G. Probiotics for the prevention of necrotising enterocolitis in very low-birth-weight infants: a meta-analysis and systematic review. Acta Paediatr 2017;106:1729–41.
19. Chang HY, Chen JH, Chang JH, Lin HC, Lin CY, Peng CC. Multiple strains probiotics appear to be the most effective probiotics in the prevention of necrotizing enterocolitis and mortality: an updated meta-analysis. PLoS One 2017;12e0171579.
20. van den Akker CH, van Goudoever JB, Shamir R, Domellof M, Embleton ND, Hojsak I, et al. Probiotics and preterm infants: a position paper by the European Society for Paediatric Gastroenterology Hepatology and Nutrition Committee on Nutrition and the European Society for Paediatric Gastroenterology Hepatology and Nutrition Working Group for Probiotics and Prebiotics. J Pediatr Gastroenterol Nutr 2020;70:664–80.
21. Abdulkadir B, Nelson A, Skeath T, Marrs EC, Perry JD, Cummings SP, et al. Routine use of probiotics in preterm infants: longitudinal impact on the microbiome and metabolome. Neonatology 2016;109:239–47.
22. Repa A, Thanhaeuser M, Endress D, Weber M, Kreissl A, Binder C, et al. Probiotics (Lactobacillus acidophilus and Bifidobacterium bifidum) prevent NEC in VLBW infants fed breast milk but not formula. Pediatr Res 2015;77:381–8.
23. Samuels N, Van De Graaf R, Been JV, De Jonge RC, Hanff LM, Wijnen RM, et al. Necrotising enterocolitis and mortality in preterm infants after introduction of probiotics: a quasiexperimental study. Sci Rep 2016;6:31643.
24. Denkel LA, Schwab F, Garten L, Geffers C, Gastmeier P, Piening B. Protective effect of dual-strain probiotics in preterm infants: a multi-center time series analysis. PLoS One 2016;11e0158136.
25. Fortmann I, Marißen J, Siller B, Spiegler J, Humberg A, Hanke K, et al. Lactobacillus acidophilus/Bifidobacterium infantis probiotics are beneficial to extremely low gestational age infants fed human milk. Nutrients 2020;12:850.
26. Robertson C, Savva GM, Clapuci R, Jones J, Maimouni H, Brown E, et al. Incidence of necrotising enterocolitis before and after introducing routine prophylactic Lactobacillus and Bifidobacterium probiotics. Arch Dis Child Fetal Neonatal Ed 2020;105:380–6.
27. Alshaikh B, Kostecky L, Blachly N, Yee W. Effect of a quality improvement project to use exclusive mother's own milk on rate of necrotizing enterocolitis in preterm infants. Breastfeed Med 2015;10:355–61.
28. Sisk PM, Lovelady CA, Dillard RG, Gruber KJ, O'Shea TM. Early human milk feeding is associated with a lower risk of necrotizing enterocolitis in very low birth weight infants. J Perinatol 2007;27:428–33.
29. Westaway JA, Huerlimann R, Kandasamy Y, Miller CM, Norton R, Staunton KM, et al. The bacterial gut microbiome of probiotic-treated very-preterm infants: changes from admission to discharge. Pediatr Res 2022;92:142–50.
30. Chang HY, Lin CY, Chiang Chiau JS, Chang JH, Hsu CH, Ko MH, et al. Probiotic supplementation modifies the gut microbiota profile of very low birth weight preterm infants during hospitalization. Pediatr Neonatol 2024;65:55–63.
31. Bäckhed F, Roswall J, Peng Y, Feng Q, Jia H, Kovatcheva-Datchary P, et al. Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe 2015;17:852.
32. Korpela K, Blakstad EW, Moltu SJ, Strømmen K, Nakstad B, Rønnestad AE, et al. Intestinal microbiota development and gestational age in preterm neonates. Sci Rep 2018;8:2453.
33. Zhang C, Li L, Jin B, Xu X, Zuo X, Li Y, et al. The effects of delivery mode on the gut microbiota and health: state of art. Front Microbiol 2021;12:724449.
34. Rios-Covian D, Langella P, Martin R. From short- to long-term effects of c-section delivery on microbiome establishment and host health. Microorganisms 2021;9:2122.
35. Zhou Y, Shan G, Sodergren E, Weinstock G, Walker WA, Gregory KE. Longitudinal analysis of the premature infant intestinal microbiome prior to necrotizing enterocolitis: a case-control study. PLoS One 2015;10e0118632.
36. Chae H, Kim SY, Kang HM, Im SA, Youn YA. Dysbiosis of the initial stool microbiota increases the risk of developing necrotizing enterocolitis or feeding intolerance in newborns. Sci Rep 2024;14:24416.
37. Alcon-Giner C, Dalby MJ, Caim S, Ketskemety J, Shaw A, Sim K, et al. Microbiota supplementation with Bifidobacterium and Lactobacillus modifies the preterm infant gut microbiota and metabolome: an observational study. Cell Rep Med 2020;1:100077.

Article information Continued

Fig. 1.

Flow diagram of the study. GA, gestational age; BW, body weight; PMA, postmenstrual age.

Fig. 2.

Taxonomy abundances of stool samples across 3 time points (V1–V3). (A) At V1 (birth), an overall predominance of Staphylococcus was seen in both groups. (B) At V2 (35 weeks' postmenstrual age), several peaks of Bifidobacterium were seen in the probiotics group along with a predominance of Klebsiella in the control group. (C) At V3 (4 months of age), an overall predominance of Bifidobacterium was noted in both groups.

Fig. 3.

Beta diversity at time point V2 (35 weeks' postmenstrual age). PERMANOVA, permutational multivariate analysis of variance; CG, control group; PG, probiotics group.

Fig. 4.

Relative abundances of Bifidobacterium (A), Lactobacillus (B), and Clostridium senso stricto 1 (C) at V2 (35 weeks' postmenstrual age).

Table 1.

Maternal and perinatal characteristics of the preterm newborns

Characteristic Probiotics group (n=33) Control group (n=35) P value
Maternal age (yr) 33.0±4.6 31.0±6.8 0.11
Prenatal steroid use 26 (72) 23 (64) 0.45
Cesarean delivery 25 (76) 23 (77) 0.46
Male sex 20 (61) 20 (57) 0.77
Gestational age (wk) 29.0±2.5 29.0±2.5 0.97
Birthweight (g) 1,124±418 1,167±354 0.65
Birth weight <1,000 g 15 (42) 12 (33) 0.47
Small for gestational age 6 (18) 4 (11) 0.43
Apgar score < 7 at 5 min 9 (25) 7 (19) 0.57
Apgar score at 5 min 8.0±1.5 8.0±2.1 0.63
Perinatal antibiotics exposure 6 (17.1) 10 (28.6) 0.26

Values are presented as mean±standard deviation or number (%).

Table 2.

Clinical outcomes

Clinical outcome Probiotics group (n=33) Control group (n=35) P value
Early-onset neonatal sepsis 19 (58) 22 (65) 0.55
Late-onset neonatal sepsis 6 (18) 3 (9) 0.26
Necrotizing enterocolitis ≥stage 2A 1 (3) 4 (12) 0.17
Days of antibiotics 2 (0–5) 3 (1–6) 0.44
Respiratory distress syndrome 14 (41) 12 (36) 0.69
Bronchopulmonary dysplasia 10 (29) 12 (35) 0.55
Days to full enteral feeding 10 (8–14) 10 (7–14) >0.99
Length of stay (day) 44 (35–69) 44 (36–80) >0.99
Death 2* (5.7) 1 (2.9) >0.99
Body weight at discharge 2,547±686 2,560±543 0.92
Body weight at 2 mo 2,853±798 (n=26) 2,760±583 (n=24) 0.64
Body weight at 4 mo 4,343±804 (n=20) 4,595±912 (n=21) 0.35

Values are presented as number (%), median (interquartile range), or mean±s tandard deviation.

*

Both died after randomization