Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2009, p. 1121-1128, Vol. 75, No. 4
0099-2240/09/$08.00+0 doi:10.1128/AEM.02359-07
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Division of Nutritional Sciences,1 Institute for Genomic Biology,2 Department of Animal Sciences,3 Department of Pathobiology, University of Illinois at Urbana-Champaign, Urbana, Illinois,4 Mead Johnson & Company, Evansville, Indiana,5 Bristol-Myers Squibb Co., Princeton, New Jersey6
Received 19 October 2007/ Accepted 14 November 2008
|
|
|---|
|
|
|---|
Clinical investigations of infant formulas supplemented with galactooligosaccharide (GOS) and fructooligosaccharide (FOS) at a range of concentrations and a range of blends have revealed increases in the numbers of fecal bifidobacteria in preterm infants (7, 18) and in the numbers of both bifidobacteria and lactobacilli in term infants (2, 3, 6, 7, 13, 14, 18, 19, 23, 24, 31) compared with infants receiving unsupplemented formulas. The Bifidobacterium levels (3, 6, 7, 14) and stool consistencies (7, 23, 24) have often been comparable for infants in both prebiotic-supplemented formula groups and human breast milk-fed reference groups. Nondigestible polydextrose (PDX) and lactulose (LOS) also stimulate bifidobacterial growth (8, 16, 29). In a recent safety and tolerance study, healthy term infants who received a formula supplemented with prebiotic blends of PDX, GOS, and LOS when they were 14 to 120 days old exhibited normal growth and stool characteristics more similar to those of breast-fed infants than to those of infants fed an unsupplemented control formula (35).
In the present study, the effects of the same prebiotic blends on the fecal bacterial populations of formula-fed infants were investigated during a 28-day feeding period using 16S rRNA gene-based real-time quantitative PCR (qPCR), 16S rRNA-based fluorescence in situ hybridization (FISH), and PCR-denaturing gradient gel electrophoresis (DGGE). These culture-independent molecular techniques are now commonly used to identify, quantify, and determine the profiles of bacteria in complex microbial communities that are difficult to assess using standard culture methods, such as the communities in the gastrointestinal tract (2, 6, 9, 13-15, 19, 22, 28, 31, 36). In this study, the Clostridium botulinum group (cluster I), the Clostridium lituseburense group (cluster XI), the Clostridium coccoides-Eubacterium rectale group (cluster XIV), and the genera Lactobacillus and Bifidobacterium in infant fecal samples were targeted with 16S rRNA gene group-specific primers, and changes in the subpopulations were evaluated by qPCR. Potential changes in the bifidobacterial subpopulation and total bacteria were also evaluated by using FISH, and Enterobacteriaceae were enumerated on selective media. Specific changes in the sum of the levels of Clostridium clusters I and XI, which include the pathogenic organisms Clostridium perfringens and Clostridium difficile, respectively, were chosen for the primary comparison in this study since consistently lower levels of Clostridium spp. have been found in breast-fed infants than in formula-fed infants (10). Differences in bacterial community composition between feeding groups and between individual participants were analyzed by using DGGE. Infant growth, stool characteristics, and tolerance of formulas were also monitored. The results were compared with the results for a parallel breast-fed reference group.
|
|
|---|
Fecal collection and sampling.
Study personnel collected infant stool specimens at baseline, on day 14, and on day 28 no more than 2 h following defecation by using chemical-free diapers (Tushies; TenderCare International, Eau Claire, WI). Samples were evenly divided and placed into two separate 15-ml transport containers, and the aliquots were treated as follows: (i) shipped under ambient conditions to a central clinical laboratory (North Coast Clinical Laboratory, Sandusky, OH) for enumeration of Enterobacteriaceae or (ii) immediately frozen at
–20°C and shipped on cold packs to a second central laboratory (Department of Animal Sciences, University of Illinois, Urbana) for DNA isolation and microbial community analysis. Following storage for approximately 3 years at –20°C, the second set of fecal samples was shipped on dry ice to a third laboratory (Department of Medical Microbiology, University Medical Center, Groningen, The Netherlands) for relative quantification of bifidobacteria by FISH.
Bacterial strains and culture conditions.
Bacterial strains used to obtain the qPCR standard curves are listed in Table 1. All Clostridium strains and Bacteroides fragilis were cultured in reinforced clostridial medium broth, Lactobacillus delbrueckii was cultured in de Man-Rogosa-Sharpe broth, and Bifidobacterium adolescentis was cultured in Belitsky minimal broth supplemented with hemin and menadione. All strains were grown at 37°C and harvested during exponential growth. Enterobacteriaceae were cultured on MacConkey agar at 37°C for at least 48 h, and the results were expressed in CFU.
|
View this table: [in a new window] |
TABLE 1. qPCR standards and detection limits
|
Real-time qPCR.
Assays were performed using the ABI PRISM 7900HT sequence detection system (Applied Biosystems, Foster City, CA). The PCR amplification protocols were the protocols described previously for total bacteria (9) and the genera Bifidobacterium (22), Lactobacillus (9), and Clostridium (32). All reaction mixtures contained 1 ng template DNA, a specific primer set (Table 1), and bovine serum albumin (0.01 mg/ml). SYBR green PCR master mixture (2x; Applied Biosystems) was used for detection of total bacteria, Bifidobacterium, and Lactobacillus. For detection of Clostridium, TaqMan PCR master mixture was used, and the internal probe was labeled with 6-carboxyfluorescein as the 5' reporter dye and 6-carboxytetramethylrhodamine as the 3' quencher dye. Standard curves were generated using serial dilutions of genomic DNA extracted from pure cultures of representative bacteria (Table 1). The initial data were expressed as target DNA concentrations, which were normalized by expression relative to the total bacterial DNA concentration.
16S rRNA gene PCR-DGGE.
Template DNA (10 ng) was amplified using universal bacterial 16S-V3 primers (Table 2) (a 40-bp GC clamp was added to the forward primer) as described previously (25). PCR amplicons were loaded in 35 to 60% linear DNA-denaturing gradient gels, and DGGE was performed using the Bio-Rad D-code system (Bio-Rad, Los Angeles, CA). The reference ladder contained PCR amplicons from pure cultures of B. fragilis, E. rectale, L. delbrueckii, and B. adolescentis from the top to the bottom (both the fourth and fifth bands were derived from B. adolescentis and consistently appeared on the gel). Following electrophoresis, gels were silver stained and scanned using a GS-710 calibrated imaging densitometer (Bio-Rad). Images were analyzed using Diversity Database software (Bio-Rad). Each band on the gel was identified regardless of the band intensity. Pairwise similarity was determined using the Dice coefficient (Cs), where patterns with no common bands received a score of 0 and identical band patterns received a score of 1.0 (21). Dendrograms were generated by using Ward's algorithm, in which samples were clustered by the degree of similarity to visualize community relatedness.
|
View this table: [in a new window] |
TABLE 2. PCR primers and probes used for PCR-DGGE, qPCR, and FISH
|
Statistical analysis.
The change in the sum for Clostridium clusters I and XI from baseline to day 28 was the primary outcome used to evaluate study formulas. A sample size of 17 infants per formula group was used to detect a difference of one standard deviation at 80% power (
= 0.05). qPCR data for bacterial subpopulations expressed as percentages of total bacteria were not normally distributed; therefore, data were expressed as median values. Changes from the baseline values among groups of formula-fed infants were analyzed using the Wilcoxon rank sum test; the sign test was used for within-group comparisons. Comparisons of values for all feeding groups at baseline, day 14, and day 28 were performed using the Wilcoxon rank sum test. Comparisons of changes in the levels of Enterobacteriaceae (log10 CFU/g stool) for formula groups compared with the baseline level were performed using analysis of covariance. A repeated-measures analysis that included all feeding groups was performed using data obtained at baseline, day 14, and day 28. Cs and Ward's algorithm calculations were performed using Diversity Database software (Bio-Rad). Cs were analyzed by using analysis of variance (ANOVA). Feeding group and age were factored into the analysis. Anthropometric measurements were converted to z-scores based on U.S. National Center for Health Statistics reference data (http:/www.cdc.gov/nccdphp/dnpa/growthcharts/resources/sas.htm) and were analyzed by using ANOVA. Race, ethnicity, gender, and study discontinuance were analyzed by using Fisher's exact test. Analyses were performed by using SAS, version 8 (Cary, NC).
|
|
|---|
The feeding groups were similar with respect to gender distribution and age of the participants. There was a statistically significant difference in race only between the BF group (83% white, 3% black, and 13% more than one race) and the control group (56% white and 44% black) (P = 0.002). The mean formula intake values were similar during the 7-day formula run-in period (195.0 to 214.8 ml/kg per day) and during the 28-day feeding trial (206.3 to 235.0 ml/kg per day). The means for weight-, length-, and head circumference-for-age z-scores at enrollment and day 28 are shown in Table 3. Statistically significant differences in mean weight- and length-for-age z-scores among feeding groups were detected at enrollment and day 28 (Table 3 shows individual P values for all pairwise comparisons of z-scores). Mean weight- and length-for-age z-scores for each group either remained the same or increased between the time of enrollment and day 28. There were no significant differences in mean head circumference-for-age z-scores at enrollment or at study day 28.
|
View this table: [in a new window] |
TABLE 3. Mean weight-, length-, and head circumference-for-age z-scores at enrollment and day 28
|
0.018) during the run-in period; there were no significant differences among feeding groups in average mean stool frequency from baseline to day 28 (1.9 to 2.8 ± 0.3 stools/day). The daily stool consistency scores (mean ± standard error) were significantly higher for the BF group (3.8 ± 0.1) than for all of the formula groups (control group, 3.3 ± 0.1; PG4 group, 3.3 ± 0.1; PGL4 group, 3.2 ± 0.1; PGL8 group, 3.1 ± 0.1) (P
0.004) during the run-in period; the scores for the BF group (3.8 ± 0.1) and each formula group were also significantly different when the averages from baseline to day 28 were compared (control group, 3.2 ± 0.1; PG4, PGL4, or PGL8 group, 3.4 ± 0.1) (P < 0.05). The only significant difference in the amount of gas (mean ± standard error) was the difference between the BF group (0.9 ± 0.1) and the PGL8 group (1.4 ± 0.1) (P = 0.006) from day 15 to day 21. There were no significant differences among the feeding groups in the amount of fussiness (mean ± standard error) during the run-in period or feeding trial (1.0 to 1.3 ± 0.1). There were no differences among the study groups in the number of participants who experienced at least one adverse event, and there were no differences in the number of participants who experienced adverse events when body systems or individual events were analyzed. There were no adverse events that had particular clinical significance in any feeding group (data not shown).
Quantification of Enterobacteriaceae.
Selective plating of fecal samples on MacConkey agar was used to enumerate fecal Enterobacteriaceae populations in all feeding groups (log10 CFU/g stool). No significant differences in the mean changes in the Enterobacteriaceae populations from baseline to day 28 were detected among the formula groups (control, PG4, or PGL4 group, 0.2 ± 0.1 log10 CFU/g stool; PGL8 group, 0.3 ± 0.1 log10 CFU/g stool). Within the PGL8 group there was a slight but statistically significant increase in the mean level of Enterobacteriaceae (P = 0.027) from baseline to day 28, but there were no significant changes within the other feeding groups. The means for Enterobacteriaceae for the different feeding groups at baseline, day 14, and day 28 are shown in Table 4. The mean for the BF group was significantly higher than the means for all other groups at baseline (P < 0.001). Significantly higher levels of Enterobacteriaceae were also observed on day 14 for the BF group than for the prebiotic-supplemented groups (P
0.002). The Enterobacteriaceae levels for the feeding groups at day 28 were similar.
|
View this table: [in a new window] |
TABLE 4. Quantification of Enterobacteriaceae in fecal samples from different feeding groups by selective enumeration on MacConkey agar
|
Median values for Clostridium clusters, Lactobacillus spp., and Bifidobacterium spp. for each feeding group at baseline, day 14, and day 28 are shown in Table 5. No significant differences among feeding groups in Clostridium cluster I values were observed at any time point. The median values for Clostridium cluster XI were significantly lower for the BF group than for the PG4 group (P = 0.006), the PGL4 group (P < 0.001), and the PGL8 group (P = 0.016) at baseline, for all other groups at day 14 (control group, P = 0.011; PG4 group, P < 0.001; PGL4 group, P = 0.003; PGL8 group, P = 0.019), and for the PG4 group (P = 0.005) and the PGL4 group (P = 0.015) at day 28. No significant differences among feeding groups for Clostridium cluster XIV, the sum for Clostridium clusters I and XI, Lactobacillus spp., or Bifidobacterium spp. were observed at any time point.
|
View this table: [in a new window] |
TABLE 5. Evaluation of Clostridium clusters I, XI, and XIVa, Lactobacillus spp., and Bifidobacterium spp. at baseline, day 14, and day 28 by qPCR and FISH
|
Profiling of bacterial communities using DGGE.
Bacterial community profiles were generated by 16S rRNA gene PCR-DGGE. Samples collected from the same subject at baseline, day 14, and day 28 were run on the same gel to maintain electrophoretic consistency (representative results are shown in Fig. 1). Cs were calculated to quantify the pairwise similarity of DNA banding patterns within (intragroup) and among (intergroup) feeding groups at baseline, day 14, and day 28. The levels of intergroup similarity at baseline (Cs range, 18.0 to 19.2), day 14 (Cs range, 19.2 to 21.7), and day 28 (Cs range, 16.4 to 19.9) were low. The intragroup Cs values were within similar ranges at all time points (data not shown). The intragroup and intergroup Cs values were not significantly different at any time point. A cluster analysis was performed using Ward's algorithm, and dendrograms were constructed to visualize relationships among banding patterns for samples. Samples obtained from the same subject at all time points examined tended to cluster together. No specific clustering patterns were observed for feeding groups, ages, or races (data not shown).
![]() View larger version (118K): [in a new window] |
FIG. 1. Representative DGGE profiles for individual infants at three study time points, baseline (lanes 1), day 14 (lanes 2), and day 28 (lanes 3). Lanes L contained a standard ladder consisting of B. fragilis, E. rectale, L. delbrueckii, and B. adolescentis (fourth and fifth bands) from top to bottom. Each collection of samples is grouped by study diet.
|
|
View this table: [in a new window] |
TABLE 6. Cs of individual DGGE profiles between baseline and day 28
|
|
|
|---|
Small but statistically significant differences in mean weight- and length-for-age z-scores were noted among feeding groups at enrollment and study day 28. For all groups, however, the mean weight- and length-for-age z-scores either remained the same or increased between enrollment and day 28, indicating that each group was growing normally. Consequently, differences among feeding groups were not considered clinically significant. The current study was not designed to evaluate growth as a primary clinical outcome. Ziegler et al. (35), however, fed infants from days 14 to 120 old the formulas given to the control, PG4, and PGL8 groups in a study designed specifically to evaluate infant growth according to the American Academy of Pediatrics Task Force recommendations for clinical testing of infant formulas (11). These workers found no differences in the rates of increases in weight or length over the duration of the 120-day study and concluded that infants fed these prebiotic-containing formulas exhibited normal growth (35). Stools from infants in the BF group in the current study were found to be softer and looser than stools from infants in groups that received formula. However, the stool consistency scores changed for the groups of infants who received prebiotic-supplemented formula over the 28-day feeding trial and approached the stool consistency scores for the breast-fed infants. Ziegler et al. (35) also reported that infants fed formula supplemented with prebiotics had stool characteristics more similar to those of breast-fed infants than to those of infants fed unsupplemented formula. Overall, infant formulas containing blends of PDX, GOS, and LOS were generally well tolerated and safe.
Few significant differences in the overall bacterial composition were found from baseline to day 28 within feeding groups or among feeding groups at any time point. Compared to the traditional culture methods used in previous prebiotic feeding studies, the more sensitive molecular methods for quantifying bacteria in complex microbial communities may make broad differences between breast-fed and formula-fed infants less apparent. For example, a recent examination of fecal samples from breast-fed and formula-fed infants using qPCR demonstrated that there were no significant differences in the Bifidobacterium levels (28). Although the levels of bifidobacteria tended to be higher in the BF group throughout the current study, the levels of neither bifidobacteria nor lactobacilli were significantly different in the formula-fed infants and the infants fed breast milk at any time point, with one exception. The exception occurred in the baseline phase, when fecal samples from infants in the PGL8 group contained significantly less bifidobacteria (20.7%) than fecal samples from infants in the human milk reference group (83.5%), as determined by FISH. The general lack of significance could have been due to the great variation in the data, which could be attributed to the wider range of infant ages (21 to 104 days at baseline) and the shorter feeding period (28 days) compared with other similar prebiotic studies.
In addition, subtle differences between the microbiota of breast-fed infants and the microbiota of formula-fed infants might occur only at the species level, where species in certain subpopulations of the gut microbiota shift toward the pattern observed in breast-fed infants (28, 30). Rinne et al. (30) demonstrated that the total number of bifidobacteria was significantly lower in infants who received control formula than in infants who received breast milk or a prebiotic-supplemented formula (GOS and FOS at a ratio of 9:1); additionally, species-specific PCR revealed that the Bifidobacterium species composition obtained for the breast-fed infants was found in infants who received prebiotic-supplemented formula. Using a species-specific qPCR assay, Penders et al. (28) found a significantly lower prevalence and significantly lower numbers of C. difficile, a member of Clostridium cluster XI, in breast-fed infants than in formula-fed infants. In the current study, a significantly greater number of subjects in the PGL8 group exhibited a decrease of the number of Clostridium cluster XI organisms between baseline and day 28. However, no species-specific prebiotic effects, such as a decrease in the level of a pathogen or similarity in species diversity, could be attributed to supplementation of infant formula with GOS, PDX, and LOS with the group-specific primers used.
The level of bifidobacteria observed expressed as a percentage of the total bacterial DNA in this study appeared to be low compared with the common assumption that Bifidobacterium spp. dominate the infant fecal microbiota. However, the results obtained by using a DNA concentration-based standard curve in the qPCR analysis may not have accurately reflected the number of bacterial cells in the complete microbial community. To estimate the number of bifidobacterial cells in the fecal microbiota, we performed post hoc FISH using Bifidobacterium-specific and universal bacterial probes. The percentages of bifidobacterial cells in the original samples indicated by the FISH results were much higher than the percentages estimated by DNA concentration-based qPCR. Consistent with the qPCR results, the percentages of bifidobacteria were higher in BF group than in the other groups, but there were few statistically significant differences among groups due to the great variation among subjects within groups. The difference between the percentages of bifidobacteria estimated by FISH and qPCR can probably be partially explained by inherent differences in methodology. Bifidobacterium spp. have relatively low numbers of 16S rRNA gene copies compared with other major intestinal bacterial species, which leads to potential underestimation by 16S rRNA gene-targeted qPCR. Palmer et al. (27) also observed a relatively low frequency and low abundance of bifidobacteria in infant fecal samples based on the results of molecular methods, including qPCR (as in this study) and microarray analysis, but in the latter analysis the primers used for broad-range amplification of rRNA genes were suboptimal and underestimated bifidobacteria eightfold. Another explanation for the low level of bifidobacteria detected by qPCR in the current study may be the difference in DNA extraction efficiency between gram-positive and gram-negative bacteria, which may result in a bacterial DNA composition bias. The differences in the percentages of bifidobacteria reported in different studies emphasize the importance of reevaluating PCR-related biases, including DNA extraction bias, in order to obtain more accurate estimates of the contribution of this bacterial group to the intestinal microbiota of infants.
An unexpected finding was that the BF group had significantly higher levels of Enterobacteriaceae than the formula-fed groups at the beginning of the study, since the levels of Enterobacteriaceae are reportedly lower in breast-fed infants (28). However, by the end of the study period, the levels of Enterobacteriaceae were comparable for all feeding groups, suggesting that neither prebiotic carbohydrates nor human milk oligosaccharides influenced the growth of these bacteria.
Analysis of the 16S rRNA gene PCR-DGGE results for infant fecal samples revealed little similarity in the bacterial community profiles among subjects regardless of the feeding practice. The calculated Cs for DGGE banding patterns within or among feeding groups were low, and significant differences were not observed. Cluster analysis of DGGE profiles demonstrated that over time the individual bacterial community profiles clustered by subject rather than by feeding group. Significant differences in bacterial community stability were not observed among feeding groups during the 28-day feeding period, whereas a significant age effect was observed. A post hoc analysis of younger and older infants within each group revealed that the shifting microbial populations in younger infants in the PG4, PGL4, and PGL8 groups, as well as the BF reference group, were less stable from baseline to day 28. In contrast, the Cs values for the control group appeared to be relatively high for both younger and older infants, indicating that the overall microbial profile was reasonably stable over the 28-day feeding period regardless of age. These age-related differences might have been due to the presence of prebiotics in the formulas and the presence of human milk oligosaccharides in breast milk, as corresponding age-related differences were not observed for the control group. Although further age-based statistical analyses were not performed, the results suggest that the prebiotic blends had a greater effect on the fecal bacterial communities of younger infants than on the fecal bacterial communities of older infants and imply that infant age is a variable that is more important than originally expected.
With the current study design, few statistically significant differences in the composition of fecal bacterial populations were observed in the infants who received formula supplemented with GOS, PDX, or LOS prebiotic blends compared with the infants who received a control formula or the reference BF group. However, the age-related differences in fecal bacterial community stability observed within the PG4, PGL4, and PGL8 feeding groups indicate that these prebiotic blends may have a greater impact on infant fecal bacterial populations in younger infants than in older infants. Further analysis is necessary before definitive conclusions concerning this possibility can be reached.
We thank study investigators and their staffs for their cooperation, including Michael Brown (Hill Top Research, St. Petersburg, FL), Ann Grandjean (The Center for Human Nutrition, Omaha, NE), William Geffen (Hillcrest Medical Group, Tulsa, OK), Robert Lewine (Hill Top Research, Scottsdale, AZ), Derek Lewis (Little Rock, AR), Susanne Matthews (Metro Pediatrics, Birmingham, AL), William Parker (Portico Pediatrics, Shreveport, LA), Ekhard Ziegler (Children's Hospital of Iowa, Iowa City, IA), and Stacy Jones (Arkansas Children's Hospital, Little Rock, AR). The participation of parents and infants in this study is gratefully acknowledged. Gerwin C. Raangs and Gjalt W. Welling (University Medical Center, Groningen, The Netherlands) are acknowledged for conducting the FISH analysis of bifidobacteria. Julia Boettcher is acknowledged for assistance with the manuscript.
Published ahead of print on 16 December 2008. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»