
==== Front
Heliyon
Heliyon
Heliyon
2405-8440
Elsevier

S2405-8440(24)13188-X
10.1016/j.heliyon.2024.e37157
e37157
Research Article
The effect of oral probiotics in the last trimester on the human milk and infant gut microbiotas at six months postpartum: A randomized controlled trial
Ma Guangyu a1
Li Yimi c1
Tye Kian Deng d
Huang Ting a
Tang Xiaomei a
Luo Huijuan a
Wang Dongju a
Zhou Juan a
Li Zhe mcal@foxmail.com
b⁎⁎
Xiao Xiaomin jnuxiaoxiaomin@163.com
a⁎
a Department of Obstetrics and Gynecology, The First Affiliated Hospital of Jinan University, Guangzhou, China
b Department of Obstetrics and Gynecology, The Third Affiliated Hospital of Sun Yat-Sen University, Guangzhou, China
c Department of Obstetrics and Gynecology, Dangyang People's Hospital, Dangyang, China
d Department of Obstetrics and Gynecology, The Fifth Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
⁎ Corresponding author. jnuxiaoxiaomin@163.com
⁎⁎ Corresponding author. mcal@foxmail.com
1 These authors contributed equally to this work.

30 8 2024
15 9 2024
30 8 2024
10 17 e371576 2 2024
30 7 2024
28 8 2024
© 2024 The Authors
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
Objective

The main aim of this study was to evaluate the effect of oral probiotics on the human milk microbiota and determine whether that influenced infant microbiota development.

Methods

A total of 27 pregnant women were recruited; 14 were assigned to the probiotic group, and the rest were assigned to the control group. Their infants were likewise assigned to the probiotic group or the control group. Pregnant women in the probiotic group received probiotic supplementation from 32 weeks of gestation until delivery. Human milk samples and infant fecal samples were collected at 6 months after delivery, and 16S rRNA sequencing was used to analyze the composition of the human milk and infant gut microbiota (NCT06241222).

Results

In the control group, bacterial microbiota were detected in 8 out of 13 milk samples, whereas in the probiotic group, only 6 out of 14 milk samples contained bacterial microbiota. We examined the composition of the human milk and infant gut microbiota in both the control and probiotic groups. Spearman correlation analysis revealed that various genera in human milk were correlated with the infant gut microbiota. The Linear discriminant analysis effect size (LEfSe) showed that 6 bacteria in the human milk microbiota in the control group were significantly more abundant than those in the probiotic group. Nine bacteria were significantly more abundant in the human milk microbiota in the probiotic group than the control group. According to the LEfSe results, 11 bacteria in the infant gut microbiota in the control group were significantly more abundant than those in the probiotic group. Fourteen bacteria were significantly more abundant in the infant gut microbiota in the probiotic group than in the control group.

Conclusion

The infant gut microbiota at 6 months has a complicated relationship with the maternal human milk microbiota. Oral probiotic supplementation can change the composition of the human milk microbiota and the infant gut microbiota.

Keywords

Gut microbiota
Infants
Milk microbiota
Pregnant women
Probiotic supplementation
==== Body
pmc1 Introduction

Human milk is a source of nutrients for neonates. It is an important factor affecting the early development of the gut microbiota in offspring and plays an important role in infant growth and immune development. Breastfeeding has been associated with a reduced incidence of neonatal respiratory and gastrointestinal infections and may also decrease the risk of obesity, asthma, Crohn's disease, ulcerative colitis, and diabetes, among other conditions [[1], [2], [3], [4]]. In addition, exclusive or partial breastfeeding can promote infant cognitive development [[5], [6], [7]]. The microbiota in human milk plays an important role in stimulating the immune system in newborns, which is conducive to the development of the intestinal microbiota with immunomodulatory activity in newborns [8].

The maternal gut microbiota may be the source of the milk microbiota. Some scholars have proposed the concept of “the bacterial entero-mammary pathway”. During pregnancy, maternal intestinal permeability increases, leading to gut bacterial translocation to the lamina propria and Peyer's plaques along with mucosal dendritic cells. Then, gut bacteria can migrate from Peyer's plaques to mesenteric lymph nodes and enter the breast via lymph vessels and peripheral blood [[9], [10], [11], [12]]. Human peripheral blood mononuclear cells (PBMCs) and milk have been found to contain bacteria and their genetic material. These bacteria are phagocytosed by PBMCs and transferred to the mammary gland intact. Currently, it is widely believed that human milk is nonsterile.

Probiotics are defined as live microorganisms that are beneficial to human health [[13], [14], [15]]. Previous work has shown that probiotic supplementation during pregnancy can reduce maternal depression or anxiety and reduce the risk of gestational obesity and gestational diabetes [[16], [17], [18]]. Probiotics can also reduce the incidence of allergic diseases such as eczema or special dermatitis during pregnancy [19,20]. The gut microbiota during pregnancy can affect the composition of the early lactation microbiota in the first two weeks of the postpartum period [21]. Human milk is the best source of nutrition for infants, and the World Health Organization advocates exclusive breastfeeding until six months of age [22]. The aim of this study was to observe the effect of oral probiotics administered during the last trimester on microbiota in human milk and infant gut. We also wanted to determine the relationship between human milk and the infant gut microbiota.

2 Materials and methods

2.1 Participants

All methods were performed in accordance with relevant guidelines and regulations for clinical trials (NCT06241222). Normal singleton pregnant women who had undergone a systematic birth examination at the First Affiliated Hospital of Jinan University and who were scheduled to give birth at our hospital were selected. The inclusion criteria were as follows: Chinese nationality, singleton pregnancy, pre-pregnancy BMI of 18.5–24 kg/m2, first normal pregnancy, infants exclusively breastfed for 6 months, and infants delivered vaginally. The exclusion criteria were women of an pregnant women over 35 years old, with a history of gastrointestinal tract issues, with a history of vaginal infection before pregnancy, and with hypertension, diabetes, immune diseases, gestational hypertension, gestational diabetes, or other pregnancy complications and those taking probiotics or antibiotics during lactation. Informed consent was obtained from every participant before 32 weeks of pregnancy.

A total of 27 pregnant women meeting the inclusion criteria were recruited and randomly assigned to either the probiotic group (14 participants) or the control group (13 participants). Pregnant women in the probiotic group were advised to consume live combined Bifidobacterium, Lactobacillus, and Streptococcus tablets twice a day after 32 weeks of pregnancy (Inner Mongolia Shuangqi Pharmaceutical, 0.5 g/tablet; each tablet contained live bacteria of Bifidobacterium longum of no less than 0.5 × 107 CFU and live Lactobacillus bulgaricus and Streptococcus thermophilus of no less than 0.5 × 106 CFU), hereinafter referred to as triple-bacteria tablets, until childbirth.

Postpartum human milk and infant fecal specimens were collected for a period of 6 months. In the control group, due to personal reasons, one infant fecal specimen was not obtained, resulting in a total collection of 13 human milk samples and 12 infant fecal samples. In the probiotic group, two infant fecal specimens were not collected, while 14 human milk samples and 12 infant fecal samples were successfully obtained. The human milk samples from the control group were designated CM, and the infant samples were designated CI; however, in the probiotic group, the human milk samples were labeled PM, and the infant samples were labeled PI. The detailed sample collection information is presented in Table 1, Table 2.Table 1 Collection of human milk samples and microbial detection status.

Table 1Milk Sample Name	Samples Collected or Not	Detection of Microbial Communities	
CM 1	Yes	Yes	
CM 2	Yes	Yes	
CM 3	Yes	Not	
CM 4	Yes	Yes	
CM 5	Yes	Not	
CM 6	Yes	Not	
CM 7	Yes	Not	
CM 8	Yes	Yes	
CM 9	Yes	Yes	
CM 10	Yes	Yes	
CM 11	Yes	Yes	
CM 12	Yes	Not	
CM 13	Yes	Yes	
PM 1	Yes	Yes	
PM 2	Yes	Yes	
PM 3	Yes	Not	
PM 4	Yes	Not	
PM 5	Yes	Not	
PM 6	Yes	Not	
PM 7	Yes	Yes	
PM 8	Yes	Not	
PM 9	Yes	Not	
PM 10	Yes	Yes	
PM 11	Yes	Not	
PM 12	Yes	Yes	
PM 13	Yes	Not	
PM 14	Yes	Yes	

Table 2 Collection of infant fecal samples and microbial detection status.

Table 2Fecal Sample Name	Samples Collected or Not	Detection of Microbial Communities	
CI 1	Yes	Yes	
CI 2	Yes	Yes	
CI 3	Yes	Yes	
CI 4	Yes	Yes	
CI 5	Yes	Yes	
CI 6	Yes	Yes	
CI 7	Yes	Yes	
CI 8	Yes	Yes	
CI 9	Yes	Yes	
CI 10	Not	–	
CI 11	Yes	Yes	
CI 12	Yes	Yes	
CI 13	Yes	Yes	
PI 1	Not	–	
PI 2	Yes	Yes	
PI 3	Yes	Yes	
PI 4	Yes	Yes	
PI 5	Yes	Yes	
PI 6	Yes	Yes	
PI 7	Yes	Yes	
PI 8	Yes	Yes	
PI 9	Yes	Yes	
PI 10	Not	–	
PI 11	Yes	Yes	
PI 12	Yes	Yes	
PI 13	Yes	Yes	
PI 14	Yes	Yes	

2.2 Sample collection

When infants were 6 months old, postpartum milk was collected by two full-time members of the research group. The human milk samples were collected according to a standardized protocol [23,24]. Before collecting the samples, the nursing mothers were required to clean the nipple and areola with soap and water. A total sample volume of 6 ml was collected in a sterile tube. After collection, the samples were stored with dry ice for transportation to the laboratory immediately. Upon arrival at the laboratory, the samples were briefly thawed to allow for accurate aliquoting. Each sample was then divided equally into smaller frozen storage tubes. The divided samples were promptly refrozen and stored in a freezer at −80 °C within 4 h after collection. This process ensured that the samples were handled with minimal temperature fluctuations to preserve their integrity.

Fecal samples from the infants were collected internally in sterile and dedicated collecting boxes, thereby avoiding any contamination with foreign material. These samples were stored in a domestic freezer (−20 °C) and transferred to a laboratory freezer at −80 °C within 24 h using a cryogenic transport container to ensure they remained frozen during the transfer.

2.3 Detection of milk and fecal microbiota

For the fecal samples, DNA was extracted by using a Magnetic Soil and Stool DNA Kit (TianGen, China, Catalog #: DP712). For the human milk samples, DNA was extracted by using the CTAB extraction method.

16S rRNA genes of distinct regions (16SV3-V4) were amplified using specific primers (515F- 806R) with barcodes. Each PCR batch included a positive control containing a known 16S rRNA gene sequence to verify the amplification efficiency and a negative control devoid of template DNA to monitor for potential contamination or false-positive results. All PCRs were conducted with 15 μL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs), 0.2 μM forward and reverse primers, and approximately 10 ng of template DNA. Thermal cycling consisted of initial denaturation at 98 °C for 1 min, followed by 30 cycles of denaturation at 98 °C for 10 s, annealing at 50 °C for 30 s, elongation at 72 °C for 30 s and extension at 72 °C for 5 min. The same volume of 1X loading buffer (containing SYB green) was mixed with the PCR products, and electrophoresis was performed on a 2 % agarose gel for detection. The PCR products were mixed in ratios with equal densities. The mixture of PCR products was subsequently purified with a Universal DNA Purification Kit (TianGen, China, Catalog #: DP214).

2.4 Data processing and analysis

Sequencing libraries were generated using the NEBNext® Ultra™ II FS DNA PCR-free Library Prep Kit (New England Biolabs, USA, Catalog#: E7430L) in accordance with the manufacturer's recommendations, and indices were added. The library was checked with a Qubit fluorometer and real-time PCR for quantification and a bioanalyzer for size distribution detection. The quantified libraries were pooled and sequenced on Illumina platforms according to the needed effective library concentration and data amount. Single-end reads assembly and quality control. Single-end reads were assigned to samples based on their unique barcode sequence. Each read was examined for the presence of its corresponding barcode, which was used to map the read to its respective sample. After assignment, the barcode and primer sequence were removed from each read using a sliding window approach with a minimum quality score of 20 (Phred quality score) and a minimum length of 50 bases to ensure precise truncation. This process resulted in high-quality reads that were ready for downstream analysis. Quality filtering of the raw reads was performed using Cutadapt (Martin M., 2011) (V1.9.1, http://cutadapt.readthedocs.io/en/stable/) quality control process to obtain high-quality clean reads. The following parameters were used: a quality cutoff of 30 (Phred quality score) to trim low-quality bases from the ends of the reads, a minimum length of 50 bases to retain trimmed reads, and adapter trimming using the specific adapter sequences provided by the sequencing platform. All other parameters were set to default values. The reads were compared with the reference database (Gold database, http://drive5.com/uchime/uchime_download.html) using the UCHIME algorithm (UCHIME Algorithm, http://www.drive5.com/usearch/manual/uchime_algo.html) to detect chimeric sequences, after which the chimeric sequences were removed. Then, the effective tags were obtained.

Sequence analysis was performed with Uparse software (Uparse v7.0.1001, http://drive5.com/uparse/). Sequences with ≥97 % similarity were assigned to the same OTUs. A representative sequence for each OTU was screened for further annotation. For each representative sequence, the Silva138 database (https://www.arb-silva.de/) was used based on the RDP classifier (version 2.2, http://sourceforge.net/projects/rdp-classifier/) algorithm to annotate taxonomic information.

To study the phylogenetic relationships of the different OTUs and the differences in the dominant species in the different samples (groups), multiple sequence alignment was conducted using MUSCLE software (version 3.8.31, http://www.drive5.com/muscle/). OTU abundance information was normalized using a standard sequence number corresponding to the sample with the least number of sequences. Subsequent analyses of alpha diversity and beta diversity were performed based on these output normalized data. Alpha diversity was applied to analyze the complexity of species diversity for a sample and was calculated with QIIME version 1.9.1 (http://qiime.org/scripts/split_libraries_fastq.html) [12] and displayed with R software (version 2.15.3).

2.5 Data normalization

The absolute abundance of ASVs was normalized using a standard sequence number corresponding to the sample with the fewest sequences. Subsequent analyses of alpha diversity and beta diversity were all performed based on the output normalized data.

2.6 Statistical analysis

Comparisons of alpha diversity indices between groups were performed using the Wilcoxon test with R software (version 2.15.3). The P-values were adjusted for false discovery rate (FDR) using the Benjamini–Hochberg method. Nonmetric multidimensional scaling (NMDS) was performed to obtain principal coordinates and visualize complex, multidimensional data. NMDS based on the unweighted UniFrac distance was calculated with QIIME software. NMDS analysis was performed with the WGCNA package, stat packages and ggplot2 package in R software. The significance of the differences in the sample communities was determined by linear discriminant analysis effect size (LEfSe) [25]. For the LEfSe analysis, we extracted the p-values and applied the Benjamini-Hochberg FDR correction. P < 0.05 after FDR correction was considered to indicate statistical significance. To estimate the correlation between genera, we calculated the Spearman coefcients for the top 100 genera by relative abundance. To control for multiple testing, we applied the Benjamini-Hochberg FDR correction to the p-values obtained from the Spearman correlation results. Genera with absolute Spearman correlation coefficients greater than 0.6 and FDR-adjusted P < 0.05 were considered for construction of the network. Network construction and visualization analysis were carried out using cytoscape software (Version 3.6). The Spearman correlation analysis was used to observe the correlation between human milk and the genera of infant gut microbiota, with a correlation coefficient threshold set at |r| > 0.6 and a significance level of FDR-adjusted P < 0.05.

3 Results

3.1 Status of 6-month postpartum milk and fecal microbiota library construction

In the CM, bacterial microbiota was detected in 8 of the 13 milk samples, and the detection rate for the milk microbiota was 61.5 %. In the PM, bacterial microbiota was detected in only 6 of the 14 milk samples, and the detection rate for the milk microbiota was 42.9 %. The remaining milk samples were not successfully established in the database because of insufficient bacterial content, they were not ideal for amplification and sequencing. All the fecal samples were successfully subjected to library construction for microbial analysis.

3.2 Constitution of the microbiota in human milk and the intestinal tract of offspring

3.2.1 Human milk microbiota at 6 months

Proteobacteria (63.2 %), Firmicutes (18.7 %), and Bacteroidetes (6.3 %) were the three most common phyla in the human milk microbiota at 6 months in the CM. Similarly, Proteobacteria (61.7 %), Firmicutes (17.1 %), and Bacteroidetes (7.5 %) were the three most common components in the human milk microbiota in the PM (Fig. 1A) (Table 3).Fig. 1 Relative abundance bar charts of the human milk microbiota species. (A) Relative abundance bar charts of the top 10 phyla in the human milk microbiota; (B) Relative abundance bar charts of the top 10 genera in the human milk microbiota. The x-axis represents sample names; the y-axis indicates the relative abundance.

Fig. 1

Table 3 Relative OTUs percentages of the top ten phyla in the human milk microbiome.

Table 3Microbiome abundance	CM (%)	PM (%)	
Proteobacteria	63.1746	61.7189	
Firmicutes	18.7010	17.0882	
Cyanobacteria	3.1978	5.4046	
Bacteroidetes	6.3067	7.4923	
Actinobacteria	4.3088	2.6032	
Acidobacteria	0.6727	1.0089	
unidentified_Bacteria	0.4942	0.3797	
Spirochaetes	0.1283	0.3827	
Tenericutes	0.5020	0.6267	
Chloroflexi	0.4280	0.6172	

Stenotrophomonas (9.2 %), Phyllobacterium (4.4 %), unidentified_Cyanobacteria (3.2 %), Streptococcus (3.0 %), and Sphingomonas (2.9 %) were the five most common genera in the human milk microbiota in the CM. In the PM, the five most common genera were Stenotrophomonas (8.9 %), Pseudomonas (8.2 %), unidentified_Cyanobacteria (5.3 %), Phyllobacterium (3.9 %) and Lactobacillus (3.3 %) (Fig. 1B) (Table 4).Table 4 Relative OTUs percentages of the top ten genera in the human milk microbiome.

Table 4Microbiome abundance	CM (%)	PM (%)	
Stenotrophomonas	9.2440	8.9274	
Pseudomonas	1.3712	8.2043	
unidentified_Cyanobacteria	3.1606	5.3248	
unidentified_Clostridiales	2.7689	0.8930	
Phyllobacterium	4.4248	3.8858	
Staphylococcus	1.7610	0.1758	
Lactobacillus	1.7672	3.3436	
Sphingomonas	2.8781	3.1777	
Streptococcus	3.0220	2.1532	
Bifidobacterium	1.4588	0.5590	

3.3 Microbiota of the infant gut at 6 months

Actinobacteria (43.1 %), Firmicutes (32.0 %), and Proteobacteria (19.0 %) were the three most common phyla in the infant gut microbiota at 6 months in the CI. Similarly, Actinobacteria (45.2 %), Proteobacteria (27.5 %), and Firmicutes (22.7 %) were the three most common phyla in the PI (Fig. 2A) (Table 5).Fig. 2 Relative abundance bar charts of the gut microbiota species. (A) Relative abundance bar charts of the top 10 phyla in the gut microbiota; (B) Relative abundance bar charts of the top 10 genera in the gut microbiota. The x-axis represents sample names; the y-axis indicates the relative abundance.

Fig. 2

Table 5 Relative OTUs percentages of the top ten phyla in the infant gut microbiota.

Table 5Microbiome abundance	CI (%)	PI (%)	
Proteobacteria	18.9798	27.5008	
Actinobacteria	43.0870	45.1788	
Firmicutes	32.0205	22.7048	
Bacteroidetes	4.9400	4.0103	
Cyanobacteria	0.6835	0.2603	
Acidobacteria	0.0388	0.0741	
Fusobacteria	0.0241	0.0381	
Chloroflexi	0.0295	0.0474	
Euryarchaeota	0.0170	0.0004	
Planctomycetes	0.0190	0.0194	

Bifidobacterium (41.6 %), Enterococcus (8.6 %), Bacteroides (4.4 %), Lactobacillus (4.2 %), and Erysipelatoclostridium (3.3 %) were the five most common genera in the infant gut microbiota in the CI. In the PI, the five most common genera were Bifidobacterium (43.9 %), Stenotrophomonas (5.0 %), Enterococcus (4.2 %), Blautia (3.2 %) and Streptococcus (3.0 %) (Fig. 2B) (Table 6).Table 6 Relative OTUs percentages of the top ten genera in the infant gut microbiota.

Table 6Microbiome abundance	CI (%)	PI (%)	
Bifidobacterium	41.6348	43.8852	
Stenotrophomonas	0.5868	5.0257	
Lactobacillus	4.1713	1.4593	
Blautia	1.6408	3.1628	
Bacteroides	4.4040	1.5141	
Enterococcus	8.5820	4.2153	
Erysipelatoclostridium	3.3076	2.0327	
Phyllobacterium	2.0424	1.3720	
Streptococcus	1.1507	3.0244	
unidentified_Clostridiales	2.4592	1.4554	

3.4 Α-Diversity of each microbiota

3.4.1 Venn graphs

The OTU results were obtained by clustering. There were 3364 OTUs shared by the human milk microbiota in both the CM and the PM (Fig. 3A). The human milk microbiota in the CM had 2425 unique OTUs, while that in the PM had 1781 unique OTUs. There were 1509 OTUs shared by the gut microbiota in both the CI and PI (Fig. 3B). In total, 941 unique OTUs were identified in the gut microbiota of the CI. The PI had 1353 unique OTUs.Fig. 3 Venn diagram. A Venn diagram shows the OTUs unique to each sample and shared by different samples. Each circle represents a group, and the numbers in the overlapping regions of the circles represent the number of shared OTUs between the groups. The numbers outside the overlapping regions represent the unique number of OTUs in each group.

Fig. 3

3.5 Rarefaction and rank abundance curves

Rarefaction curves and rank abundance curves were used to analyze the α diversity of all the groups. The rarefaction curve directly reflects the rationality of the sequencing data volume; the horizontal axis represents the randomly sampled sequencing data quantity, and the vertical axis represents the observed OTU count. With increasing sequencing quantity, the trend of increase in the rarefaction curves became more gradual (Fig. 4A). In this study, the sequencing data met the sequencing requirements, and the sequencing data volume per sample was reasonable. The rank abundance curve reflects the evenness of species in the samples. In the horizontal direction, the width of the curve reflects the abundance of the species, with a wider range on the x-axis indicating greater species abundance. The shape (smoothness) of the curve reflects the evenness of the species in the sample; the smoother the curve is, the more evenly distributed the species are. The rank abundance curve showed a gradual decrease, indicating a uniform distribution of species (Fig. 4B).Fig. 4 Rarefaction curve and rank abundance curve. (A) Rarefaction curve. Each curve in the diagram represents a group sample; different colors represent different groups. A flat curve indicates that a greater sequence number will lead to a lower OTU number. (B) Rank abundance curve. Each curve in the diagram represents a group sample; different colors represent different groups.

Fig. 4

3.5.1 Chao1 and simpson

The Chao1 index, which is used to estimate the number of species in the community, showed no significant differences between the CM and PM in the human milk microbiota (P = 0.850) (Fig. 5A). The Simpson index, which is used to evaluate the species diversity index of the microbiota, showed no significant differences between the CM and PM in terms of the human milk microbiota (P = 0.315) (Fig. 5B). The Chao1 index did not significantly differ between the CI and PI for the infant gut microbiota (P = 0.420) (Fig. 5C). The Simpson index did not significantly differ between the CI and PI for the infant gut microbiota (P = 0.120) (Fig. 5D).Fig. 5 (A, C) Chao1 index and (B, D) Simpson index. The Chao1 index is used to compare the species richness of gut microbiota between different groups. A higher Chao1 index indicates greater species richness. The Simpson index is used to assess the diversity of gut microbiota between different groups. A lower Simpson index value indicates higher diversity, meaning a more even distribution of species.

Fig. 5

3.6 NMDS analysis

The NMDS analysis (stress = 0.053) indicated that the human milk microbiota compositions of the CM and PM were distinct, while the intestinal microbiota compositions of the CI and PI were similar (Fig. 6).Fig. 6 NMDS analysis of the microbiota composition. NMDS plots showing the differences in gut microbiota composition between groups. Each point represents a sample, and the distance between points reflects the similarity in microbiota composition. Stress values indicate the goodness-of-fit of the NMDS representation. Different colors represent different groups.

Fig. 6

3.7 Differences between the groups at 6 months

LEfSe was used to determine the significant differences in abundance between the different groups. As shown in Fig. 7, Hyphomonadacea, Bifidobacterium_bifidum, Stenotrophobacter, Epulopiscium, Niameybacter_massiliensis, Erysipelatoclostridium and Erysipelatoclostridium_ ramosum were significantly more abundant in the human milk microbiota of the CM than in that of the PM (P < 0.05) (Fig. 7A). Moraxella, Moraxellaceae, Moraxella_catarrhalis, Mycoplasma, Mycoplasma_hyorhinis, Mycoplasma_pneumoniae, Bergeyella, Roseburia and Roseburia_intestinalis were significantly more abundant in the human milk microbiota in the PM than in the CM (P < 0.05) (Fig. 7A).Fig. 7 Analysis of differences in the microbiota between the two groups. (A). LDA analysis results for the bacterial milk microbiota; (B). LDA analysis results for the gut microbiota. The bar chart in the LDA score distribution displays species with LDA scores greater than the set threshold; these are biomarkers indicating statistically significant differences between groups. The length of the bars represents the magnitude of the impact of the differential species (i.e., LDA score). p__: Phylum, c__: Class, o__: Order, f__: Family, g__: Genus, s__: Species.

Fig. 7

Firmicutes, Anaerostipess_caccae, Pygmaiobacter_massiliensis, Pygmaiobacter, Fonticella, TM7_phylum_sp_oral_clone_FR058, Armatimonadetes, Intestinibacter, Lactobacillus_backii, Megasphaera_micronuciformis, and Megasphaera were significantly more abundant in the infant gut microbiota in the CI than in the PI (P < 0.05) (Fig. 7B). Desulfomicrobium, Desulfomicrobiaceae, Oscillibacter_sp_ER4, Peptostreptococcus, Chryseolinea, Unidentified_Verrucomicrobiae, Peptoniphilus_sp_EL1, Peptoniphilus, Pararhizobium_giardinii, Streptosporangiales, Microlunatus, Dehalococcoidia, Steptococcus, and Streptococcaceae were significantly more abundant in the infant gut microbiota in the PI than in the CI (P < 0.05) (Fig. 7B).

3.8 Human milk and infant gut microbiotas co-occurrence network analysis

Microbial co-occurrence networks demonstrate the relationships of mutual influence and restriction between microbial communities (Fig. 8, Fig. 9). To further prove whether probiotic supplementation contributes to changes in the bacterial interaction network, we calculated the network index for each group (Supplemental Table 1). The network index can provide insight into the overall structure and complexity of the microbial communities, which can be influenced by probiotic intervention.Fig. 8 Co-occurrence networks of milk microbiota at the genus level. (A) Co-occurrence network of the bacterial milk microbiota in CM. (B) Co-occurrence network of the bacterial milk microbiota in PM. The nodes in the above figure are identified by different colors, representing each dominant genus. The line between the nodes indicates a correlation between the two genera, the red line indicates a positive correlation, the blue line indicates a negative correlation, and the thickness of the line indicates the intensity of correlation. The node size represents the content of the OTUs, and the dots of the same color represent the genus classified into the same phylum. The more connections through a node, the more associations the genus has with other members of the genus.

Fig. 8

Fig. 9 Co-occurrence networks of gut microbiota at the genus level. (A) Co-occurrence network of the infant gut microbiota in CI. (B) Co-occurrence network of the infant gut microbiota in PI. The nodes in the above figure are identified by different colors, representing each dominant genus. The line between the nodes indicates a correlation between the two genera, the red line indicates a positive correlation, the blue line indicates a negative correlation, and the thickness of the line indicates the intensity of correlation. The node size represents the content of the OTUs, and the dots of the same color represent the genus classified into the same phylum. The more connections through a node, the more associations the genus has with other members of the genus.

Fig. 9

Generic relationships refer to the number of direct interactions or edges between nodes (microorganisms) in the network. Probiotic supplementation increases the number of interactions in infant microbiotas (PI vs. CI) but decreases them in maternal milk (PM vs. CM) (Fig. 10A).Fig. 10 Network Index of Human Milk and Infant Gut Microbiotas. (A) Generic relationships. Shows the overall structure and connections among different microbiota species or genera. (B) Network diameter. Represents the longest shortest path between any two nodes, indicating the maximum extent of the network. (C) Clustering coefficient. Indicates the tendency of nodes to cluster together. A high coefficient shows tightly knit groups. (D) Graph density. The ratio of actual edges to the maximum possible edges, reflecting overall connectivity. (E) Average degree. The average number of connections per node, showing the level of interaction among microbiota. (F) Average path length. The average number of steps along the shortest paths between all node pairs, indicating the efficiency of interaction flow within the network.

Fig. 10

The network diameter is the longest of all the shortest paths between any two nodes in the network, providing a measure of the network's overall size and connectedness. A smaller network diameter suggests a more compact network. Probiotic groups (PI and PM) had smaller diameters compared to their control counterparts, indicating a more closely connected network (Fig. 10B).

The clustering coefficient is a measure of the degree to which nodes in a network tend to cluster together. Higher clustering coefficient indicate a more complex network with stronger interactions among microorganisms. Based on the clustering coefficients (CM: 0.513, PM: 0.541, CI: 0.568, and PI: 0.572), the clustering coefficient increased from 0.5128 in the control group to 0.5406 in the probiotic group for maternal milk, indicating an enhancement in network complexity due to probiotic supplementation (Fig. 10C). However, the impact of probiotics on the network structure of the infant microbiota was less pronounced. The clustering coefficient showed a slight increase from 0.568 in the control group to 0.572 in the probiotic group, suggesting that while probiotics do influence the infant microbiota network, the effect was not as significant as in maternal milk.

Graph density measures the proportion of potential connections in a network that are actual connections. Higher density indicates a more interconnected network. Graph density was higher in the probiotic groups, particularly in PI, indicating more interconnected networks (Fig. 10D).

The average degree measures the average number of connections each node (representing a microorganism) has in the network. In the PI, the average degree was 18.845, higher than the 15.250 in the CI, indicating more interactions per microorganism (Fig. 10E). In the PM, the average degree was 8.020, slightly lower than the 9.734 in the CM, suggesting balanced but robust interactions.

The average path length represents the average number of steps along the shortest paths for all possible pairs of nodes. In the PI, the average path length was 2.976, slightly shorter than the 3.025 in the CI, indicating more direct connections (Fig. 10F). In the PM, the average path length was 3.533, shorter than the 4.763 in the CM, suggesting a more efficiently connected network.

3.9 Correlation analysis between human milk microbiota and infant gut microbiota

Using Spearman correlation analysis, a heatmap was created to identify the top 20 genera of the human milk microbiota and the top 20 genera of the intestinal microbiota in 6-month-old infants, in order to explore the regulatory effects of the human milk microbiota on the infant gut microbiota.

The Spearman ranking for the top 20 genera between PM and PI did not find any genera showing correlation between the milk microbiome and the infant gut microbiome (Fig. 11A). However, various genera in human milk were correlated with the infant gut microbiota between CM and CI (Fig. 11B). We found that the abundance of Bifidobacterium in human milk was negatively correlated with the abundance of Streptococcus and Veillonella in infant feces (Table 7). The abundance of Streptococcus in human milk was positively correlated with the abundance of Blautia in infant feces. The abundance of Sphingomonas in human milk was positively correlated with the abundance of Bifidobacterium in infant feces.Fig. 11 Spearman Correlation Heatmap Analysis of the Top 20 Genera Levels in the Milk Microbiome and Infant Gut Microbiome in Breastfed Infants. (A) Correlation Heatmap between CM and CI. (B) Correlation Heatmap between PM and PI. M: Human milk with bacteria at 6 months postpartum. I: Infant gut microbiotas.

Fig. 11

Table 7 Correlation analysis of microbiota OTUs between the CM and the CI according to the top 20 genera (n = 7).

Table 7CM	CI	R	p Value	
Lachnoclostridium	Romboutsia	−0.847	0.016	
Delftia	Stenotrophomonas	0.847	0.016	
Delftia	Flavonifractor	0.786	0.048	
Vibrio	Erysipelatoclostridium	0.821	0.034	
Vibrio	Phyllobacterium	0.793	0.033	
Vibrio	Flavonifractor	0.893	0.012	
unidentified_Lachnospiraceae	unidentified_Clostridiales	−0.964	0.003	
unidentified_Rhizobiaceae	Enterococcus	−0.857	0.024	
unidentified_Rhizobiaceae	Erysipelatoclostridium	0.821	0.034	
unidentified_Rhizobiaceae	Romboutsia	0.857	0.024	
Acinetobacter	Sellimonas	0.889	0.007	
Enterococcus	Phyllobacterium	−0.77	0.041	
Veillonella	Streptococcus	−0.829	0.021	
Streptococcus	Blautia	0.893	0.012	
Bifidobacterium	Streptococcus	−0.857	0.024	
Bifidobacterium	Veillonella	−0.829	0.021	
Phyllobacterium	Stenotrophomonas	0.811	0.027	
Sphingomonas	Bifidobacterium	−0.857	0.024	
Staphylococcus	Streptococcus	0.893	0.012	
Stenotrophomonas	Enterococcus	−0.929	0.007	
Stenotrophomonas	Erysipelatoclostridium	0.857	0.024	
Stenotrophomonas	Phyllobacterium	0.793	0.033	
Stenotrophomonas	unidentified_Ruminococcaceae	0.786	0.048	
Stenotrophomonas	Flavonifractor	0.857	0.024	
unidentified_Clostridiales	Enterococcus	0.829	0.021	
unidentified_Clostridiales	Erysipelatoclostridium	−0.775	0.041	
unidentified_Clostridiales	Romboutsia	−0.775	0.041	
unidentified_Clostridiales	Streptococcus	−0.811	0.027	

4 Discussion

Previous studies have shown that human milk is nonsterile. The presence of bacteria in sterile-collected milk samples has already been confirmed through bacterial culture methods [26]. In recent years, the development of research technologies such as Illumina sequencing has shown that highly diverse bacterial communities do exist in human milk samples through metagenomic studies [27]. The results of this experiment showed that among the 13 CM milk samples, microbiota was detected in 8 of them. In the PM samples, microbiota was detected in 6 out of 14 samples. The remaining milk samples could not be successfully sequenced, possibly due to an insufficient microbial content for optimal amplification and sequencing. In our study, Proteobacteria, Firmicutes, and Bacteroidetes were the three most common phyla in the human milk microbiota at 6 months in the CM. Stenotrophomonas, Phyllobacterium, unidentified_Cyanobacteria, Streptococcus, and Sphingomonas were the five most common genera in the human milk microbiota in the CM. The human milk microbiota in the CM group had 2425 unique OTUs. Research by Fitzstevens et al. suggested that facultative anaerobic or aerobic bacteria are the predominant colonizers of the human milk ecosystem. Streptococcus and Staphylococcus are the most common and abundant bacterial classes in human milk samples, followed by Propionibacterium and Bifidobacterium, which originate from the skin or the environment; additionally, well-known probiotic species such as Lactobacillus and Bifidobacterium can be detected in the gut [28]. Previous work showed similar results to our results, with small differences. These differences may be influenced by race, region, dietary habits, the environment and many other factors.

The intestinal microbiota of infants undergoes a dynamic and evolving process within the first year of life [[29], [30], [31]]. A study by Raspini et al. [32] revealed that the dominant phyla in the intestinal microbiota of exclusively breastfed 6-month-old infants were Firmicutes, Actinobacteria and Proteobacteria. In this study, at the phylum level, Actinobacteria, Firmicutes, and Proteobacteria were identified as the dominant bacterial groups in the CI, contrary to previous findings. At the genus level, Bifidobacterium, Enterococcus, Bacteroides, Lactobacillus, and Erysipelatoclostridium were the five most common genera in the infant gut microbiota in the CI group. The previous study included 6-month-old infants from Italy, and we speculate that the potential factors contributing to the differences in the infant gut microbiota could be variations in complementary feeding practices, regional disparities, and ethnic diversity.

The early establishment of the gut microbiota may directly impact future health. After birth, the intestinal microbiota of newborns is susceptible to various factors, including the mode of delivery, feeding method, gestational age, and use of antibiotics by pregnant women during pregnancy [33]. Human milk provides various important symbiotic and mutually beneficial dominant bacterial groups, such as Bifidobacterium, Phyllobacterium, Lactobacillus, and Clostridium, for the establishment of the intestinal microbiota in newborns. Bifidobacterium, the primary dominant colonizer, collaborates with other bacteria in the intestines to resist harmful microorganisms such as Enterobacter and Salmonella, stimulating the production of secretory immunoglobulin A (IgA) and activating T regulatory cells [34]. This collaboration promotes the development of the newborn's immune system, prevents infections, and ensures the stability of the intestinal microbiota in newborns. Human milk is the best source of nutrition for infants and is also one of the carriers of bacteria that can be transmitted from mother to baby [35,36]. In our study, we found that the infant gut microbiota at 6 months has a complicated relationship with the maternal human milk microbiota. For instance, the abundance of Bifidobacterium in human milk is negatively correlated with the abundance of Streptococcus and Veillonella in infant feces. The network graph showed that although Probiotic supplementation did not directly alter the abundance of probiotics, they could modify the microbial microenvironment. To further demonstrate the impact of probiotic supplementation on bacterial interaction networks, we calculated the network index for each group. This index provides insights into the structure and complexity of microbial communities [37,38]. Network index analysis reveals that probiotic supplementation significantly enhances microbial network complexity and interactions in maternal milk. This is evidenced by a smaller network diameter, indicating a more compact network, and a higher clustering coefficient, suggesting stronger node interactions [39,40]. Additionally, the network exhibits a smaller average degree and graph density, reflecting fewer interactions between nodes and less dense connectivity [41]. The shorter average path length further supports the increased complexity and robustness of microbial interactions in maternal milk due to probiotics [42]. In the infant gut microbiota, while the impact of probiotics is noticeable, it is less pronounced compared to maternal milk. The network in the infant gut microbiota shows higher graph density and average degree, indicating more interconnected and interactive nodes. It also has a shorter average path length, signifying more efficient communication between nodes. The clustering coefficient is slightly higher, suggesting somewhat stronger interactions among nodes. Overall, although probiotic supplementation does enhance microbial network complexity and interactions in the infant gut, the changes are less dramatic than in maternal milk, highlighting a more substantial impact of probiotics on the maternal milk microbiome. Not all studies have found changes in specific microbiota following probiotic supplementation. The conventional view holds that taking probiotics can increase the presence of beneficial bacteria. Research suggested that supplementing with oral probiotics during late pregnancy and the breastfeeding period can positively impact the beneficial microbiota in human milk. In a group taking VSL#3 (which includes Lactobacillus, Bifidobacterium, and Streptococcus thermophilus), significant increases in the amounts of Bifidobacterium and Lactobacillus were observed in both colostrum and mature milk compared to the control group [43]. However, in most cases, taking probiotics could change the corresponding microbiota, but there were some instances where it cannot be changed [44]. This relationship showed that the microbiota in human milk does not directly influence the infant gut microbiota by changing the abundance of the bacteria themselves. This influence may be reflected in variations in the infant gut environment, inducing gut microbiota variation. Our study aim was to evaluate the probiotic effect on mother milk and infant gut microbiota, as well as the correlation between mother milk and infant gut microbiota. So, we want to exclude other confounding factor to make our manuscript more rigorous. Further study is needed to explore the other factors which will influence infant gut microbiota. In addition, this study collected human milk samples from disinfected breasts, eliminating the influence of skin microbiota. However, before collecting milk, infants suckling human milk make contact with the nipple and surrounding skin. Is it possible that the combination of low CFU content in human milk and high levels of skin microbiome means that maternal skin is a more important source for infant gut colonization than human milk? The relevant mechanism of this relationship requires further investigation. The establishment of a stable gut microbiota occurs during two critical periods of transition [45]. The first transition takes place shortly after birth when newborns shift from maternal dependence to relying on human milk. The second transition happens during weaning, marked by the end of breastfeeding and the introduction of solid foods, and these changes persist until the age of three. After this point, children develop a gut microbiota structure similar to that of adults. Many studies showed that the influence of human milk microbiota, as well as the skin microbiome, was much more apparent in the early stages (2–8 weeks postpartum). However, between 3 and 6 months, particularly when early weaning begins, a significant shift in the gut microbiome occurs. This study collected human milk samples at 6 months postpartum, and infants were exclusively breastfed during this period. It might be too late to study the effect of milk microbiota on gut microbiota establishment at 6 months postpartum. The dominant genera in the infant gut can vary depending on the feeding method. In 6-month-old infants who were formula-fed, the dominant genera included Streptococcus, Klebsiella, and Lactobacillus. In mixed-fed infants, the dominant genera were Streptococcus, Klebsiella, and Bifidobacterium [46]. In this study, Bifidobacterium was the dominant genus in the infant gut microbiota in both the CI and PI groups, while in another similar study with exclusively breastfed infants, Bacteroides was the dominant genus, with Bifidobacterium ranking third [46]. The reason for the difference could be the variation in geographic location and environmental factors.

Previously, we investigated the impact of taking probiotics during pregnancy on the gut microbiota of pregnant women. Probiotics were taken from 32 weeks of gestation until delivery. Our study found that taking probiotics during pregnancy did not alter the alpha diversity of the gut microbiota but did affect beta diversity. Additionally, the control group showed significant depletion of Turicibacter and Clostridium_sensu_stricto, while the probiotic group experienced a significant reduction in Phascolarctobacterium [47]. Does supplementation with probiotics affect the establishment of milk at 6 months postpartum? In our study, at the phylum level, Proteobacteria, Firmicutes, and Bacteroidetes were the three most common phyla in the human milk microbiota in the PI, similar to those in the CI. The five most common genera were Stenotrophomonas, Pseudomonas, unidentified_Cyanobacteria, Phyllobacterium and Lactobacillus in the PI. However, the five most common genera were Bifidobacterium, Stenotrophomonas, Enterococcus, Blautia and Streptococcus in the CI. At the phylum level, the consumption of probiotics by pregnant women had a relatively minor impact on the composition of the infant gut microbiota at 6 months. However, at the genus level, it exerted a more significant influence on the composition of the infant gut microbiota. LEfSe was used to determine the significant differences in abundance between the different groups. The six bacteria in the human milk microbiota in the CM were significantly more abundant than those in the PM. Nine bacteria were significantly more abundant in the human milk microbiota in the PM than in the CM. According to the LEfSe results, 11 bacteria in the infant gut microbiota in the CI were significantly more abundant than those in the PI. Fourteen bacteria were significantly more abundant in the infant gut microbiota in the PI than in the CI. Probiotics cannot change the microbiota constitution directly but can change the microenvironment through their metabolites or receptor-mediated genes.

The supplementation of probiotics is a controversial topic. There have been varying opinions and research findings on whether probiotics have an impact on the structure of the microbiome. Research suggested that supplementing with oral probiotics during late pregnancy and the breastfeeding period can positively impact the beneficial microbiota in human milk. In a group taking VSL#3 (which includes Lactobacillus, Bifidobacterium, and Streptococcus thermophilus), significant increases in the amounts of Bifidobacterium and Lactobacillus were observed in both colostrum and mature milk compared to the control group [43]. A survey included 24 randomized controlled trials with a total of 2761 mothers and 1756 infants. The overall effect of probiotics on the detection rate of beneficial bacteria in human milk showed a risk difference of 24 % [48]. However, some research suggested that maternal probiotic supplementation did not significantly affect the overall composition of the human milk microbiota transferred during breastfeeding [[49], [50], [51]]. Our study did not find beneficial effects of probiotic supplementation on the beneficial bacteria in human milk and infant feces through analysis of LEfSe differential results. This may be related to differences in probiotic strains, dosage, administration timing, and the timing of sample collection. Further research is needed to determine the specific impact of probiotic supplementation on the abundance of beneficial bacteria in human milk and infant feces.

Our study had several limitations. We acknowledge that a major limitation of this study is the small sample size. Most ongoing studies in this area involve several hundred or even over 1000 participants to achieve sufficient statistical power for studying the high complexity of neonatal gut microbiome development. This limitation may impact the generalizability of our findings. Future studies should aim to recruit larger populations. We acknowledge the importance of a double-blind placebo-controlled design for the study's robustness. Because pregnant women are a special vulnerable group, using placebos may increase their anxiety. Therefore, we chose not to use placebos in our study. In future research, we will introduce double-blind animal experiments to better study the mechanisms involved. We recognize that another limitation of this study is the selection of a timepoint at 6 months postpartum. The influence of milk microbiota may be more direct in the earlier stages of lactation. Given that the probiotic intervention occurred during the third trimester of pregnancy, and previous studies have demonstrated that changes in the microbiota due to intervention can last several weeks before potentially reverting back to their original state, there is a question of whether any alterations in the maternal gut microbiome persist into the lactational stage, particularly up to 6 months postpartum. Spearman correlation analysis can help identify correlations between microbial communities. However, microbiome data is typically compositional, meaning each sample's total relative abundance sums to one. This constraint can lead to spurious correlations, as increases in the relative abundance of one bacterium must correspond to decreases in others, regardless of biological significance. Thus, Spearman correlations may capture spurious associations driven by these inherent constraints. For example, two unrelated bacterial species might show significant negative correlation due to the sum constraint. Additionally, microbiome data often contain many zero values, indicating that some bacteria are not detected in certain samples. Spearman correlation can be biased by these zeros, as it cannot distinguish between true absence and zeros due to detection limits. While Spearman correlation has limitations and may not always be reliable, many studies still use it to analyze relationships between microbial communities [52]. We will seek better methods in future studies to explore the relationship between human milk and the infant gut microbiota. This study lacks data on the maternal gut microbiota at 6 months postpartum, which is a key area for further investigation.

5 Conclusion

In this study, we demonstrated the constitution of both the human milk and infant gut microbiotas at six months postpartum. The infant gut microbiota at 6 months has a complicated relationship with the maternal human milk microbiota. At the phylum level, the consumption of probiotics by pregnant women has a relatively minor impact on the composition of the infant gut microbiota at 6 months. However, at the genus level, it exerts a more significant influence on the composition of the infant gut microbiota. However, further investigations are needed to determine the underlying mechanism involved.

Funding

This study was supported by the Department of Science and Technology of Guangzhou (Grant no. SL2022A04J00792 ) and the Research Collaboration Agreement from the Third Affiliated Hospital of Sun Yat-sen University (Grant no. 202403282179 ), and the 10.13039/100014717 National Natural Science Foundation of China (Grant no. 81771664 ).

Ethics approval statement

The study project was authorized by the Institutional Review Board (IRB) for Human Subject Research at the First Affiliated Hospital of Jinan University (2019–011) and performed following the tenets of the Declaration of Helsinki. Written informed consent was obtained from every pregnant woman.

Consent to participate

Not applicable.

Data availability statement

All data generated or analyzed during this study are included in this published article. The sequence data reported in this study was archived in the Sequence Read Archive with the accession number PRJNA1067545.

Code availability

Not applicable.

CRediT authorship contribution statement

Guangyu Ma: Writing – review & editing, Writing – original draft, Software, Methodology. Yimi Li: Supervision, Investigation. Kian Deng Tye: Resources, Data curation. Ting Huang: Software, Investigation. Xiaomei Tang: Software, Resources. Huijuan Luo: Visualization, Supervision. Dongju Wang: Formal analysis, Conceptualization. Juan Zhou: Methodology, Investigation. Zhe Li: Writing – review & editing, Writing – original draft, Resources, Methodology, Conceptualization. Xiaomin Xiao: Supervision, Project administration, Investigation, Funding acquisition, Conceptualization.

Declaration of competing interest

The authors declare the following financial interests/personal relationships which may be considered as potential competing interests:

Xiaomin Xiao reports financial support was provided by 10.13039/501100001809 National Natural Science Foundation of China . If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A Supplementary data

The following is the Supplementary data to this article:Multimedia component 1

Multimedia component 1

Acknowledgments

We would like to thank Biobank of the first affiliated hospital of Jinan University for sample storage.

Appendix A Supplementary data to this article can be found online at https://doi.org/10.1016/j.heliyon.2024.e37157.
==== Refs
References

1 Horta B.L. Loret de Mola C. Victora C.G. Long-term consequences of breastfeeding on cholesterol, obesity, systolic blood pressure and type 2 diabetes: a systematic review and meta-analysis Acta Paediatr. Oslo Nor 104 2015 30 37 10.1111/apa.13133 1992
2 Klopp A. Vehling L. Becker A.B. Subbarao P. Mandhane P.J. Turvey S.E. Lefebvre D.L. Sears M.R. Daley D. Silverman F. Hayglass K. Kobor M. Turvey S. Kollmann T. Brook J. Ramsey C. Macri J. Sandford A. Pare P. Tebbutt S. Brauer M. Denburg J. Cyr M. Kozyrskyj A. Becker A. Chen E. Miller G. Takaro T. Ratjen F. Grasemann H. Dell S. To T. Moraes T. Subbarao P. Holness L. Befus A.D. Mandhane P. Laprise C. Scott J. Hegele R. Lou W. Anand S. Sears M. Simons E. Azad M. Eiwegger T. Duan Q. Hystad P. Scott J. Azad M.B. Modes of infant feeding and the risk of childhood asthma: a prospective birth cohort study J. Pediatr. 190 2017 192 199.e2 10.1016/j.jpeds.2017.07.012 29144244
3 Wang L. Collins C. Ratliff M. Xie B. Wang Y. Breastfeeding reduces childhood obesity risks, child Obes. Print 13 2017 197 204 10.1089/chi.2016.0210
4 Xu L. Lochhead P. Ko Y. Claggett B. Leong R.W. Ananthakrishnan A.N. Systematic review with meta-analysis: breastfeeding and the risk of Crohn's disease and ulcerative colitis, Aliment Pharmacol. Ther. 46 2017 780 789 10.1111/apt.14291
5 Michaelsen K.F. Lauritzen L. Mortensen E.L. Effects of breast-feeding on cognitive function Adv. Exp. Med. Biol. 639 2009 199 215 10.1007/978-1-4020-8749-3_15 19227544
6 Ibrahim S.H. Ghany S.M.A.E.- Shafie T.M.E. Hady M.E. Cognitive functions in breastfed versus artificially fed in preschool children, Egypt J. Hosp. Med. 77 2019 5742 5751 https://www.semanticscholar.org/paper/Cognitive-Functions-in-Breastfed-versus-Fed-in-Ibrahim-Ghany/1239438c98e19b8cd4b9c93c6687643a9bfaeda3
7 Pang W.W. Tan P.T. Cai S. Fok D. Chua M.C. Lim S.B. Shek L.P. Chan S.-Y. Tan K.H. Yap F. Gluckman P.D. Godfrey K.M. Meaney M.J. Broekman B.F.P. Kramer M.S. Chong Y.-S. Rifkin-Graboi A. Nutrients or nursing? Understanding how breast milk feeding affects child cognition Eur. J. Nutr. 59 2020 609 619 10.1007/s00394-019-01929-2 30809702
8 Toscano M. De Grandi R. Grossi E. Drago L. Role of the human breast milk-associated microbiota on the newborns' immune system: a mini review Front. Microbiol. 8 2017 2100 10.3389/fmicb.2017.02100 29118752
9 Gueimonde M. Laitinen K. Salminen S. Isolauri E. Breast milk: a source of bifidobacteria for infant gut development and maturation? Neonatology 92 2007 64 66 10.1159/000100088 17596738
10 Perez P.F. Doré J. Leclerc M. Levenez F. Benyacoub J. Serrant P. Segura-Roggero I. Schiffrin E.J. Donnet-Hughes A. Bacterial imprinting of the neonatal immune system: lessons from maternal cells? Pediatrics 119 2007 e724 e732 10.1542/peds.2006-1649 17332189
11 Donnet-Hughes A. Perez P.F. Doré J. Leclerc M. Levenez F. Benyacoub J. Serrant P. Segura-Roggero I. Schiffrin E.J. Potential role of the intestinal microbiota of the mother in neonatal immune education Proc. Nutr. Soc. 69 2010 407 415 10.1017/S0029665110001898 20633308
12 Rodríguez J.M. The origin of human milk bacteria: is there a bacterial entero-mammary pathway during late pregnancy and lactation? Adv. Nutr. Bethesda Md 5 2014 779 784 10.3945/an.114.007229
13 Sanders M.E. Probiotics: definition, sources, selection, and uses Clin. Infect. Dis. Off. Publ. Infect. Dis. Soc. Am. 46 Suppl 2 2008 S58 10.1086/523341 61; discussion S144-151
14 Morelli L. Yogurt, living cultures, and gut health Am. J. Clin. Nutr. 99 2014 1248S 10.3945/ajcn.113.073072 50S 24695895
15 Zok C. Die intestinale Mikrobiota: Das Faktotum im menschlichen Darm, DMW - Dtsch Med. Wochenschr. 139 2014 1282 1283 10.1055/s-0034-1374690
16 Wickens K.L. Barthow C.A. Murphy R. Abels P.R. Maude R.M. Stone P.R. Mitchell E.A. Stanley T.V. Purdie G.L. Kang J.M. Hood F.E. Rowden J.L. Barnes P.K. Fitzharris P.F. Crane J. Early pregnancy probiotic supplementation with Lactobacillus rhamnosus HN001 may reduce the prevalence of gestational diabetes mellitus: a randomised controlled trial Br. J. Nutr. 117 2017 804 813 10.1017/S0007114517000289 28367765
17 Browne P.D. Bolte A. Claassen E. de Weerth C. Probiotics in pregnancy: protocol of a double-blind randomized controlled pilot trial for pregnant women with depression and anxiety (PIP pilot trial) Trials 20 2019 440 10.1186/s13063-019-3389-1 31315657
18 Okesene-Gafa K.A. Moore A.E. Jordan V. McCowan L. Crowther C.A. Probiotic treatment for women with gestational diabetes to improve maternal and infant health and well-being Cochrane Database Syst. Rev. 6 2020 CD012970 10.1002/14651858.CD012970.pub2 32575163
19 Huurre A. Laitinen K. Rautava S. Korkeamäki M. Isolauri E. Impact of maternal atopy and probiotic supplementation during pregnancy on infant sensitization: a double-blind placebo-controlled study Clin. Exp. Allergy J. Br. Soc. Allergy Clin. Immunol. 38 2008 1342 1348 10.1111/j.1365-2222.2008.03008.x
20 Boyle R.J. Ismail I.H. Kivivuori S. Licciardi P.V. Robins-Browne R.M. Mah L.-J. Axelrad C. Moore S. Donath S. Carlin J.B. Lahtinen S.J. Tang M.L.K. Lactobacillus GG treatment during pregnancy for the prevention of eczema: a randomized controlled trial Allergy 66 2011 509 516 10.1111/j.1398-9995.2010.02507.x 21121927
21 Huang T. Zeng Z. Liang X. Tang X. Luo H. Wang D. Zhou J. Xiao X. Effect of breast milk with or without bacteria on infant gut microbiota BMC Pregnancy Childbirth 22 2022 595 10.1186/s12884-022-04930-6 35883060
22 Infant and Young Child Feeding: Model Chapter for Textbooks for Medical Students and Allied Health Professionals 2009 World Health Organization Geneva http://www.ncbi.nlm.nih.gov/books/NBK148965/
23 Cortes-Macías E. Selma-Royo M. García-Mantrana I. Calatayud M. González S. Martínez-Costa C. Collado M.C. Maternal diet shapes the breast milk microbiota composition and diversity: impact of mode of delivery and antibiotic exposure J. Nutr. 151 2020 330 340 10.1093/jn/nxaa310
24 Yeruva L. Mulakala B.K. Rajasundaram D. Gonzalez S. Cabrera-Rubio R. Martínez-Costa C. Collado M.C. Human milk miRNAs associate to maternal dietary nutrients, milk microbiota, infant gut microbiota and growth Clin. Nutr. 42 2023 2528 2539 10.1016/j.clnu.2023.10.011 37931372
25 Chang F. He S. Dang C. Assisted selection of biomarkers by linear discriminant analysis effect size (LEfSe) in microbiome data J. Vis. Exp. JoVE 2022 10.3791/61715
26 Ward T.L. Hosid S. Ioshikhes I. Altosaar I. Human milk metagenome: a functional capacity analysis BMC Microbiol. 13 2013 116 10.1186/1471-2180-13-116 23705844
27 Zivkovic A.M. German J.B. Lebrilla C.B. Mills D.A. Human milk glycobiome and its impact on the infant gastrointestinal microbiota Proc. Natl. Acad. Sci. U. S. A. 108 Suppl 1 2011 4653 4658 10.1073/pnas.1000083107 20679197
28 Fitzstevens J.L. Smith K.C. Hagadorn J.I. Caimano M.J. Matson A.P. Brownell E.A. Systematic review of the human milk microbiota Nutr. Clin. Pract. Off. Publ. Am. Soc. Parenter. Enter. Nutr. 32 2017 354 364 10.1177/0884533616670150
29 Tapiainen T. Paalanne N. Tejesvi M.V. Koivusaari P. Korpela K. Pokka T. Salo J. Kaukola T. Pirttilä A.M. Uhari M. Renko M. Maternal influence on the fetal microbiome in a population-based study of the first-pass meconium Pediatr. Res. 84 2018 371 379 10.1038/pr.2018.29 29538354
30 Perez-Muñoz M.E. Arrieta M.-C. Ramer-Tait A.E. Walter J. A critical assessment of the “sterile womb” and “in utero colonization” hypotheses: implications for research on the pioneer infant microbiome Microbiome 5 2017 48 10.1186/s40168-017-0268-4 28454555
31 Ma G. Shi Y. Meng L. Fan H. Tang X. Luo H. Wang D. Zhou J. Xiao X. Factors affecting the early establishment of neonatal intestinal flora and its intervention measures Front. Cell. Infect. Microbiol. 13 2023 10.3389/fcimb.2023.1295111
32 Raspini B. Porri D. De Giuseppe R. Chieppa M. Liso M. Cerbo R.M. Civardi E. Garofoli F. Monti M.C. Vacca M. De Angelis M. Cena H. Prenatal and postnatal determinants in shaping offspring's microbiome in the first 1000 days: study protocol and preliminary results at one month of life, Ital J. Pediatr. 46 2020 45 10.1186/s13052-020-0794-8
33 Kolokotroni O. Middleton N. Gavatha M. Lamnisos D. Priftis K.N. Yiallouros P.K. Asthma and atopy in children born by caesarean section: effect modification by family history of allergies - a population based cross-sectional study BMC Pediatr. 12 2012 179 10.1186/1471-2431-12-179 23153011
34 Ahrné S. Lönnermark E. Wold A.E. Aberg N. Hesselmar B. Saalman R. Strannegård I.-L. Molin G. Adlerberth I. Lactobacilli in the intestinal microbiota of Swedish infants Microbes Infect 7 2005 1256 1262 10.1016/j.micinf.2005.04.011 16002310
35 Heikkilä M.P. Saris P.E.J. Inhibition of Staphylococcus aureus by the commensal bacteria of human milk J. Appl. Microbiol. 95 2003 471 478 10.1046/j.1365-2672.2003.02002.x 12911694
36 Saturio S. Nogacka A.M. Alvarado-Jasso G.M. Salazar N. de Los Reyes-Gavilán C.G. Gueimonde M. Arboleya S. Role of bifidobacteria on infant health Microorganisms 9 2021 2415 10.3390/microorganisms9122415 34946017
37 Park T. Kim H.-J. Myeong N.R. Lee H.G. Kwack I. Lee J. Kim B.J. Sul W.J. An S. Collapse of human scalp microbiome network in dandruff and seborrhoeic dermatitis Exp. Dermatol. 26 2017 835 838 10.1111/exd.13293 28094891
38 Li H. Li T. Berasategui A. Rui J. Zhang X. Li C. Xiao Z. Li X. Gut region influences the diversity and interactions of bacterial communities in pikas (Ochotona curzoniae and Ochotona daurica) FEMS Microbiol. Ecol. 93 2017 10.1093/femsec/fix149
39 Cornell C.R. Zhang Y. Ning D. Xiao N. Wagle P. Xiao X. Zhou J. Land use conversion increases network complexity and stability of soil microbial communities in a temperate grassland ISME J. 17 2023 2210 2220 10.1038/s41396-023-01521-x 37833523
40 Chen X. Li Q. Chen D. Zhao L. Xiao C. Restoration measures of fencing after tilling guided succession of grassland soil microbial community structure to natural grassland in the sanjiangyuan agro-pasture ecotone of the qinghai-Tibetan plateau Microb. Ecol. 86 2023 2870 2881 10.1007/s00248-023-02287-7 37620628
41 Li H. Wang Y. Yu Q. Feng T. Zhou R. Shao L. Qu J. Li N. Bo T. Zhou H. Elevation is associated with human skin microbiomes Microorganisms 7 2019 611 10.3390/microorganisms7120611 31771258
42 Guo B. Zhang L. Sun H. Gao M. Yu N. Zhang Q. Mou A. Liu Y. Microbial co-occurrence network topological properties link with reactor parameters and reveal importance of low-abundance genera NPJ Biofilms Microbiomes 8 2022 3 10.1038/s41522-021-00263-y 35039527
43 Mastromarino P. Capobianco D. Miccheli A. Praticò G. Campagna G. Laforgia N. Capursi T. Baldassarre M.E. Administration of a multistrain probiotic product (VSL#3) to women in the perinatal period differentially affects breast milk beneficial microbiota in relation to mode of delivery Pharmacol. Res. 95–96 2015 63 70 10.1016/j.phrs.2015.03.013
44 Ferro L.E. Crowley L.N. Bittinger K. Friedman E.S. Decker J.E. Russel K. Katz S. Kim J.K. Trabulsi J.C. Effects of prebiotics, probiotics, and synbiotics on the infant gut microbiota and other health outcomes: a systematic review Crit. Rev. Food Sci. Nutr. 63 2023 5620 5642 10.1080/10408398.2021.2022595 37667870
45 Tanaka M. Nakayama J. Development of the gut microbiota in infancy and its impact on health in later life Allergol. Int. Off. J. Jpn. Soc. Allergol. 66 2017 515 522 10.1016/j.alit.2017.07.010
46 Fan W. Huo G. Li X. Yang L. Duan C. Impact of Diet in Shaping Gut Microbiota Revealed by a Comparative Study in Infants during the Six Months of Life vol. 24 2014 133 143 10.4014/jmb.1309.09029
47 Chen Y. Li Z. Tye K.D. Luo H. Tang X. Liao Y. Wang D. Zhou J. Yang P. Li Y. Su Y. Xiao X. Probiotic supplementation during human pregnancy affects the gut microbiota and immune status Front. Cell. Infect. Microbiol. 9 2019 254 10.3389/fcimb.2019.00254 31380297
48 Alemu B.K. Azeze G.G. Wu L. Lau S.L. Wang C.C. Wang Y. Effects of maternal probiotic supplementation on breast milk microbiome and infant gut microbiome and health: a systematic review and meta-analysis of randomized controlled trials Am. J. Obstet. Gynecol. MFM 5 2023 101148 10.1016/j.ajogmf.2023.101148
49 Simpson M.R. Avershina E. Storrø O. Johnsen R. Rudi K. Øien T. Breastfeeding-associated microbiota in human milk following supplementation with Lactobacillus rhamnosus GG, Lactobacillus acidophilus La-5, and Bifidobacterium animalis ssp. lactis Bb-12 J. Dairy Sci. 101 2018 889 899 10.3168/jds.2017-13411 29248229
50 Bisanz J.E. Enos M.K. PrayGod G. Seney S. Macklaim J.M. Chilton S. Willner D. Knight R. Fusch C. Fusch G. Gloor G.B. Burton J.P. Reid G. Microbiota at multiple body sites during pregnancy in a rural Tanzanian population and effects of moringa-supplemented probiotic yogurt Appl. Environ. Microbiol. 81 2015 4965 4975 10.1128/AEM.00780-15 25979893
51 Azagra-Boronat I. Tres A. Massot-Cladera M. Franch À. Castell M. Guardiola F. Pérez-Cano F.J. Rodríguez-Lagunas M.J. Lactobacillus fermentum CECT5716 supplementation in rats during pregnancy and lactation affects mammary milk composition J. Dairy Sci. 103 2020 2982 2992 10.3168/jds.2019-17384 32008776
52 Jia Y. He T. Wu D. Tong J. Zhu J. Li Z. Dong J. The treatment of Qibai Pingfei Capsule on chronic obstructive pulmonary disease may be mediated by Th17/Treg balance and gut-lung axis microbiota J. Transl. Med. 20 2022 281 10.1186/s12967-022-03481-w 35729584
