Comparison of the gut microbiota profile in breast-fed and formula-fed Korean infants using pyrosequencing.

BACKGROUND/
OBJECTIVES: Feeding in infancy is the most significant determinant of the intestinal microbiota in early life. The aim of this study was to determine the gut microbiota of Korean infants and compare the microbiota obtained between breast-fed and formula-fed Korean infants. SUBJECTS/
METHODS: We analyzed the microbial communities in fecal samples collected from twenty 4-week old Korean (ten samples in each breast-fed or formula-fed) infants using pyrosequencing.
RESULTS: The fecal microbiota of the 4-week-old Korean infants consisted of the three phyla Actinobacteria, Firmicutes, and Proteobacteria. In addition, five species, including Bifidocbacterium longum, Streptococcus salivarius, Strepotococcus lactarius, Streptococcus pseudopneumoniae, and Lactobacillus gasseri were common commensal intestinal microbiota in all infants. The predominant intestinal microbiota in the breast-fed infants (BFI) included the phylum Actinobacteria (average 70.55%), family Bifidobacteriacea (70.12%), genus Bifidobacterium (70.03%) and species Bifidobacterium longum (69.96%). In the microbiota from the formula-fed infants (FFI), the proportion of the phylum Actinobacteria (40.68%) was less, whereas the proportions of Firmicutes (45.38%) and Proteobacteria (13.85%) as well as the diversity of each taxonomic level were greater, compared to those of the BFI. The probiotic species found in the 4-week-old Korean infants were Bifidobacterium longum, Streptococcus salivarius, and Lactobacillus gasseri. These probiotic species accounted for 93.81% of the microbiota from the BFI, while only 63.80% of the microbiota from the FFI. In particular, B. longum was more abundant in BFI (69.96%) than in FFI (34.17%).
CONCLUSIONS: Breast milk supports the growth of B. longum and inhibits others. To the best of our knowledge, this study was the first attempt to analyze the gut microbiota of healthy Korean infants according to the feeding type using pyrosequencing. Our data can be used as a basis for further studies to investigate the development of intestinal microbiota with aging and disease status.

INTRODUCTION

The human intestinal microbiota is a complex and dynamic ecosystem with 300-500 different species of bacteria [1] and plays an important role in maintaining host health since it is involved in nutrition, pathogenesis, and immunology [2].

The gut microbiota of an individual is influenced by different dietary habits during life. The gut microbiota of children fed a modern western diet compared to a rural diet has been characterized [3], and the effect of dietary polyphenols on the human gut microbiota has been reported [4]. Arumugam et al. have proposed that human gut microbial communities represent three predominant variants, "enterotypes," dominated by Bacteroides, Prevotella, and Ruminococcus, respectively [5]. In addition, the enterotype is related to disease states, and long-term dietary interventions may allow modulation of an individual's microbiota, which subsequently improves health [6].

The neonatal period is a crucial stage for gastrointestinal colonization [7], which influences the adult intestinal microbiota and lifelong health [8]. In general, the fetus is known to be sterile, and gut colonization in newborn infants starts rapidly within a few hours by contact with the maternal vagina and the environment [9]. The critical stage of gut colonization occurs within the first few weeks after birth and is influenced by a variety of factors including the environment, diet, gestational age, hospitalization, antibiotics, delivery mode, and feeding method [10]. In early life, one of the most significant determinants of the intestinal microbiota is the type of feeding.

Breast milk is the gold standard for growing infants, because it is a complex species-specific biological fluid adapted to satisfy the nutritional requirements and develop the immune system [11]. In addition, breast milk has been reported to protect against different infectious diseases [12,13,14]. Colostrum and mature milk contain potential antimicrobial factors such as immunoglobulins, cytokines, oligosaccharides, lysozyme, lactoferrin, other glycoproteins, and antimicrobial peptides that can inactivate pathogens individually or synergistically [15].

Recently barcoded pyrosequencing based on the small subunit ribosomal RNA (16S rRNA) gene has been used to analyze the entire intestinal microbiota. This powerful technique provides information regarding the complex microbiota, which had not been previously detected by culture methods [7,16,17]. In this study, we compared the gut microbiota of Korean infants fed breast milk versus formula milk by using pyrosequencing. This pilot study will allow us to better understand the relationship between the gut microbiota and feeding type, and more importantly, to Korean infant health.

SUBJECTS AND METHODS

Subjects and sampling

Fecal samples were collected from twenty healthy newborn infants at 4 weeks of age who were delivered from healthy mothers by normal spontaneous vaginal delivery at Arante Women's Hospital in Seoul, Korea, from May 2012 to February 2013. Ten infants were exclusively breast-fed, and the other ten infants were formula-fed (some exposure to breast milk but fed 70-100% with formula milk). Meconium and fecal samples were collected from the infants, which were defecated spontaneously on disposable diapers at 4 weeks after birth. The specimens were immediately sealed in a plastic bag and were stored at -80℃. The fecal samples from infants who were delivered by Caesarean section, premature, or give antibiotic prophylaxis or therapy as well as those whose mothers had infections during pregnancy, clinical illness, antibiotic administration, or probiotic supplementation were excluded from this study. This study was approved by the Institutional Review Board of Ewha Womans University (IRB No. 2012-3-10; Seoul, Korea), and all of the participants provided informed consent.

Genomic DNA extraction

Metagenomic DNA was isolated from fecal samples of infants using a Fast DNA SPIN Kit (MP BIO, Santa Ana, CA, USA). Extraction was performed according to the manufacturer's instructions. DNA was eluted in 50 µl of elution buffer and stored at -20℃ prior to use. The genomic DNA concentration was determined by the optimal density at 260 nm, and the purity was assessed by the ratio of the absorbance at 260 nm and 280 nm.

Polymerase chain reaction (PCR)

The PCR mixture was prepared with Red-Taq Ready Mix (Sigma-Aldrich, St. Louis, MO), 40 µM primers (COSMO Genetech, Korea), and 100 ng of genomic DNA in a total volume of 25 µl. For PCR, the primer set 27F and 518R was used to amplify the 27-518 region of the bacterial 16S rRNA containing the V1 and V3 regions [18].

Pyrosequencing and sequence analysis

The PCR products were purified with the QIAquick PCR purification kit (Qiagen), and quantified using the PicoGreen dsDNA Assay kit (Invitrogen, Carlsbad, CA). Equimolar concentration of each amplicon were mixed and sequenced on a Roche/454 GS Junior system at ChunLab, Inc (Seoul, Korea).

The sequences obtained from the pyrosequencer were analyzed according to previously described methods [19]. Briefly, the sequences of each sample were sorted by a unique barcode, and low quality reads (average quality score < 25 or read length < 300bp) were removed for further analysis. Primer sequences were trimmed by pairwise alignment, and trimmed sequences were clustered to correct for sequencing errors. Representative sequences in each cluster were identified using the EzTaxon-e database [20]. Chimeric sequences were removed using the UCHIME program [21] and the diversity indices were calculated using the Mothur program [22]. The pyrosequences obtained in this study are available in the EMBL SRA database under the study PRJEB6615 (http://www.ebi.ac.uk/ena/data/view/PRJEB6615).

Statistical analysis

The CLcommunity™ program was used to analyze the pyrosequencing sequences of the intestinal microbiota. The operational taxonomic unit (OTU) refers to a group of sequences that are mathematically defined with a sequence similarity of 97%. This diversity measure was carried out using CD-HIT methods. Alpha-diversity indices such as ACE, Chao1, and Shannon diversity were used to estimate species richness. All data were expressed as mean ± standard error of the mean (SEM) for each group. Comparison between the two groups of infants was performed with the Wilcoxon-Mann-Whitney test. All statistical analyses were performed with PASW Statistics 18.0 (SPSS Inc., Chicago, IL, USA). P-values < 0.05 were considered statistically significant.

RESULTS

General characteristics of the participants

There were no significant differences in terms of the characteristics between the test groups of breast-fed infants (BFI) and formula-fed infants (FFI) (Table 1). The mean maternal age was 31.4 ± 3.3 years old in the BFI group and 32.3 ± 2.3 years old in the FFI group. There were more females than males in the BFI group, whereas there were equal numbers of females and males in the FFI group. The mean birth weights were 3.1 ± 0.4 kg and 3.3 ± 0.5 kg, respectively. All infants were 4 weeks of age. No infants received any additional foods and no mothers took antibiotics after the first week of delivery.

Table 1

General characteristic of mothers and infants

1)BFI: Breast-Fed Infants

2)FFI: Formula-Fed Infants

3)Values are mean ± SD

Pyrosequences obtained from the fecal samples of BFI and FFI

A total of 207,596 sequences were used for the final analysis (Table 2). The average numbers of sequences were 11,100 from the BFI and 9,600 from the FFI, indicating no significant difference (P > 0.05). The Good's coverages of all samples were greater than 0.98, and the total number of OTUs obtained was 737 (average, 73.7 per sample) in the BFI samples and 1,364 (average, 136.4 per sample) in the FFI samples (Table 3). The Shannon diversity indices in the BFI ranged from 0.82 to 1.91, while those in the FFI were from 1.56 to 2.21. These diversity indices indicated that the intestinal microbiota obtained from the FFI was more diverse that from the BFI.

Table 2

Summary of pyrosequences obtained from fecal samples of BFI and FFI

1)BFI: Breast-Fed Infants, FFI: Formula-Fed Infants

2)The OTUs are defined based on the sequence similarity of 97% results

Table 3

The operational taxonomic units (OTUs) and Shannon diversity index of the fecal microbiota of 4 wks old BFI and FFI

The number in the parentheses ( ) is the median value and the numbers in the brackets [ ] are the minimum and maximum value.

1)BFI: Breast-Fed Infants

2)FFI: Formula-Fed Infants

3)Significantly different by independent t-test between two groups (P < 0.05, P < 0.001)

Comparison of taxonomic compositions between BFI and FFI

The analyzed sequences of microbiota were classified from the species up to the phylum level. The microbiota of the infant fecal samples consisted of the three phyla Actinobacteria, Firmicutes, and Proteobacteria (Table 4 & Fig. 1). Actinobacteria (average 70.55% of total reads) was the predominant phylum in the BFI, followed by Firmicutes (26.15%) and Proteobacteria (3.29%). In contrast, in the FFI, the proportions of Actinobacteria (40.68%) and Firmicutes (45.38%) were similar, followed by Proteobacteria (13.85%). The relative abundance of Firmicutes and Proteobacteria was higher in the FFI than in the BFI.

Table 4

Composition of fecal microbiota in 4 wks old BFI and FFI1)

1)Cut-off: 0.05

2)BFI: Breast-Fed Infants

3)FFI: Formula-Fed infants

4)Significantly different by independent t-test between two groups (P < 0.05, P < 0.001)

Fig. 1

Aggregate microbiota composition at pylum/family level from fecal samples of 4 wks old infants.

(A) Breast-fed infants, (B) Formula-fed infants.

Analysis at the family level revealed that the microbiota from the FFI was more diverse than that from the BFI. Bifidobacteriaceae was the most predominant family within the phylum Actinobacteria (Table 4 & Fig. 1), and comprised 70.12% of the BFI samples, compared to 36.61% of the FFI samples. In Actinobacteria from the FFI, Actinomycetaceae (3.73%) was the second most predominant family after Bifidobacteriaceae. The phylum Firmicutes from the BFI consisted mainly of Streptococcaceae (17.31%) and Lactobacillaceae (8.02%), while the five families Streptococcaceae (28.86%), Veillonellaceae (6.12%), Enterococcaceae (4.96%), Lactobacillaceae (2.56%), and Coprobacillus (1.94%) were abundantly detected within Firmicutes from the FFI. Within the phylum Proteobacteria, Enterobacteriaceae was the major family in both the BFI and FFI samples. The proportion of Enterobacteriaceae from the FFI (12.64%) was higher than that from the BFI (3.17%).

The fecal microbiota from individual BFI and FFI was analyzed to the genus level (Tables 5 & 6). Bifidobacterium was remarkably abundant in the BFI, with an average percentage of 70.03%. In BFI 7, 8, 9, and 10, the percentage of Bifidobacterium was below the average amount; instead, the genera Streptococcus and Lactobacillus were more abundant, indicating that total amount of probiotics was stable (Table 5). The main genera in the FFI were Bifidobacterium (average, 35.99%) and Streptococcus (average, 23.10%); however, their percentages were quite varied in each individual (Table 6). The proportion of Bifidobacterium was lower than in the FFI than in the BFI. In addition, the microbiota of the FFI contained a significant amount of the genera Escherichia, Veillonella, Enterococcus, and Enterobacter, whereas the content of Lactobacillus was low.

Table 5

The proportion of the most ubiquitous genera identified from the fecal microbiota of breast-fed infants (BFI)

1)Number of fecal samples containing the bacterial group

2)Values are %

Table 6

The proportion of the most ubiquitous genera identified from the fecal microbiota of formula-fed infants (FFI)

1)Number of fecal samples containing the bacterial group

2)Values are %

The species of fecal microbiota in individual infants were compared by heatmap analysis (Fig. 2). The species compositions from the FFI were more diverse than those from the BFI. The five species Bifidocbacterium longum, Streptococcus salivarius, Strepotococcus lactarius, Streptococcus pseudopneumoniae, and Lactobacillus gasseri were found in all infants; therefore, these five species are common commensal intestinal microbiota in 4-week-old Korean infants, independently of feeding type. In addition to these five species, Escherichia coli, Enterococcus faecalis, and Veillonella atypical were abundant in most of the FFI microbiota.

Fig. 2

Heatmap analysis of the species detected from pyrosequencing of fecal samples from 4 wks old BFI and FFI.

The data contained in the matrix are represented by black and white colors. The minimum ratio of a taxon is displayed as an individual sample. The cut-off value was set to 1%. BFI, Breast-fed infants; FFI, Formula-fed infants

Comparison of probiotics in the fecal microbiota between BFI and FFI

Probiotics are defined as living microorganisms that confer a health benefit on the host when consumed [23]. S. salivarius, L. gasseri, and species of Bifidobacterium were determined to be the dominant probiotics from the infant microbiota in this study. Comparison of probiotics between the BFI and FFI showed that the probiotics in the BFI samples (93.8%) were more abundant than the FFI samples (65.6%) (Fig. 3). The proportion of B. longum was higher in the BFI (70.0%) than in the FFI (34.2%), and the proportion of L. gasseri was also higher in the BFI (7.9%) than in the FFI (2.5%). On the other hand, the proportion of S. salivarius from the FFI (27.2%) was higher than that from the BFI (16.0%). Diverse species of Bifidobacteria such as B. breve, B. bifidum, L. brevis, and L. reuteri were detected only in the FFI. The relative abundance of the genus Bifidobacterium showed a significant difference between the BFI (average 70.0%) and the FFI (average 36.4%) (Wilcoxon-Mann-Whitney test, P = 0.0015) (Fig. 4).

Fig. 3

The compositions of probiotics at the species level in fecal samples from 4-week-old BFI and FFI.

etc., bacteria other than probiotics; BFI, Breast-fed infants; FFI, Formula-fed infants

Fig. 4

Scatter plots showing the proportions of Bifidobacterium sequences with respect to the total sequences obtained from pyrosequencing 16S rRNA genes associated with feeding type.

Mean values (horizontal lines) and significance values (P, Mann-Whitney) are shown. BFI, Breast-fed infants; FFI, Formula-fed infants

Overall, the dominant fecal bacteria from the 4-week-old Korean infants studied were the probiotics Bifidobacterium longum, Streptococcus salivarius, and Lactobacillus gasseri. In the BFI, probiotics consisted of 93.81% of the total microbiota, of which B. longum was the most predominant species.

DISCUSSION

In the past few years, determination of the composition of infant intestinal microbiota has made progress through the use of pyrosequencing-based monitoring [24]. Accordingly, the role of the infant intestinal microbiota has been emphasized not only for intestinal health but also lifelong health in adulthood. The present study showed that the fecal microbiota of 4-week-old Korean infants consisted of the three phyla Actinobacteria, Firmicutes, and Proteobacteria in both the BFI and FFI groups. Staphylococcus epidermis was also detected in the BFI, which apparently originated from the mother's skin during breast-feeding [25]. In the microbiota from the BFI, the phylum Actinobacteria (70.55%), family Bifidobacteriacea (70.12%), genus Bifidobacterium (70.03%), and species Bifidobacterium longum (69.96%) were predominant. In the microbiota from the FFI, the proportion of Actinobacteria (40.68%) was less, whereas the proportion of Firmicutes (45.38%) and Proteobacteria (13.85%) as well as the diversity of each taxonomic level were greater than in the BFI. These results demonstrated a specific gut microbiota community for Korean 4-week-old infants that is different from 6-week-old European infants [7], in which the genera Bifidobacterium (40.0%), Bacteroides (11.4%) and enterobacteria (7.5%) were detected. Fallani et al. (2010) have described that there is a geographic difference in the proportion of Bacteroides, Bifidobacteria, and Enterobacteria: higher proportions of Bifidobacteria were found in 6-week-old infants from Northern European countries, whereas a more diverse microbiota with a higher proportion of Bacteroides was found in infants from Southern European countries [7]. They found that Bifidobacteria was dominant in BFI, whereas FFI harbored higher proportions of Bacteroides, Clostridium, and Lactobacillus. In addition, Chinese infants have Firmicutes and Proteobacteria as dominant microbiota in all BFI, FFI, and mixed-fed infants aged 1 to 6 months [26]. They found a significant depletion of Bacteroides and Actinobacteria in the FFI, whereas Bifidobacteriaceae comprised only 8.16% of the microbiota in the BFI. The differences between these results and our results may be explained by the different geographical origins, feeding types, delivery modes, and ages of the infants.

The probiotic species found in this study were Bifidobacterium longum, Streptococcus salivarius, and Lactobacillus gasseri. The sum of the probiotic species comprised 93.81% of the microbiota from the BFI, while only 63.80% of that from the FFI. In particular, B. longum was significantly more abundant in the BFI (69.96%) compared with the FFI (34.17%). The FFI exhibited a more diverse microbiota community composed of various bacterial species compared to the BFI. These results are consistent with previous studies showing that the dominant bacteria in BFI are Bifidobacteriaceae and Lactobacillaceae; however, the FFI had multiple intestinal microbiota composed of Bifidobacteria, Bacteroides, Clostridia, Enterococcus, and Staphylococcus [27,28,29]. The increase in the microbiota diversity of the FFI probably depends on the contents of the formula milk; therefore, further studies are required to investigate the composition of formula milk.

The predominance of Bifidobacteriaceae due to the high proportion of B. longum in the BFI is likely related to the higher diversity of microbiota in the FFI. Breast milk contains a low buffering capacity, casein, calcium, and phosphate, but a high content of lactose, which seem to favor the growth of Bifidobacterium compared to other bacteria [30]. Breast milk differs from cow milks as it contains sialylated and fucosylated oligosaccharides, which are utilized for growth of Bifidobacteria [31,32]. The colonization of Lactobacilli and Bifidobacteria in the intestine inhibits the growth of pathogenic microorganisms, such as Staphylococcus aureus, Salmonella typhimurium, Yersinia enterocolitica, and Clostridium perfringens, by nutrient competition and preventing the adhesion of pathogens [33]. Therefore, the relationship between various contents of formula with the infant intestinal microbiota should be studied in detail.

The bacterial microbiota formed in infants depends on the diet of the infants [34], but the environment during birth, prematurity, and hygiene measures also influence the intestinal microbiota community in early infancy [10]. To minimize the effects of factors other than feeding, we only included subjects who did not take antibiotics after the first week of birth and were delivered via normal spontaneous vaginal delivery. In addition, all subjects were full-term infants.

There were several limitations in this study. Only a relatively small number of samples were analyzed, since pyrosequencing requires a high cost. In addition, samples were analyzed only at a single time point. Meconium was also analyzed, but there were not enough bacteria for DNA extraction. Lastly, some FFI were exposed to breast milk during the 4 weeks of life, and the formula compositions were not investigated for the FFI group.

Maternal dietary supplementation, e.g., with probiotics, may provide a health benefit to infants via improvement of the maternal gastrointestinal tract [35]. However, little is known about the microbiota of infants breast fed by mothers with different diets. Furthermore, infant formula resembling breast milk or the direct addition of Bifidobacterium to infant formula leads to a gut microbiota of FFI similar to that of BFI [36]. Therefore, studying the influence of the mother's dietary behaviors and types of formula on the infant's intestinal microbiota for a longer duration is required.

In conclusion, the present study was the first attempt to analyze the gut microbiota in healthy Korean 4-week-old infants by using pyrosequencing. Furthermore, this study provides a basis for a more comprehensive understanding of the diversity of the intestinal microbiota and valuable information for future studies on the role of the intestinal microbiota in lifelong health.