Skip to main content
  • Research article
  • Open access
  • Published:

Evaluation of stool microbiota signatures in two cohorts of Asian (Singapore and Indonesia) newborns at risk of atopy



Studies have suggested that demographic and lifestyle factors could shape the composition of fecal microbiota in early life. This study evaluated infant stool microbiota signatures in two Asian populations, Singapore (n = 42) and Indonesia (n = 32) with contrasting socioeconomic development, and examined the putative influences of demographic factors on these human fecal associated bacterial signatures.


Longitudinal analysis showed associations of geographical origin with Clostridium leptum, Atopobium and Bifidobacterium groups. Mode of delivery had the largest effect on stool microbiota signatures influencing the abundance of four bacterial groups. Significantly higher abundance of bacterial members belonging to the Bacteroides-Prevotella, Bifidobacterium and Atopobium groups, but lower abundance of Lactobacilli-Enterococci group members, were observed in vaginal delivered compared to caesarean delivered infants. Demographic factors influencing the structure of infants stool microbiota during the first year of life included breastfeeding, age of weaning, sibship size and exposure to antibiotics.


Differences in stool microbiota signatures were observed in relation to various demographic factors. These features may confound studies relating to the association of the structure of fecal microbiota and the predisposition to human modern disease.


The human gut microbiome is a complex ecosystem harbouring a rich diversity of commensal microorganisms. It is widely thought that the early life development of the neonatal intestinal microbiota plays an important role in the maturation of the host immune system and could in turn influence allergy development [13]. For example, germfree mice which lack the endemic intestinal microbiota showed impairment of intestinal mucosal and systemic immune system development. The impairment in the systemic immune system is reflected by poorly formed spleen and lymph nodes, hypoplastic Peyer's patches, reduced levels of secreted IgA and IgG, and lack of expansion of CD4+ T cell populations [2, 3]. Furthermore, these mice exhibited cytokine profiles that skewed towards Th2 [2], which is involved in the pathophysiology of allergic diseases.

Past studies have further reported that intestinal microbiota in subjects with allergy, particularly those with atopic eczema, differed from those of healthy controls [47]. Wang and colleagues showed that there is a reduced bacterial diversity in the early stool microbiota of infants with atopic eczema [7]. Recently, we further showed that the abundances of Bifidobacterium and Enterobacteriaceae were different among caesarean-delivered infants with and without eczema [5]. The association with specific gut microbiota signatures and the occurrence of allergy in infants has also been reported in studies comparing healthy and allergic infants from different countries [8, 9].

However, varying demographic and lifestyle characteristics at different geographical locations can pose as potential confounders in correlating multidimensional data generated from studies involving the diverse bacterial populations of the gut microbiota as "quantitative traits". For example, factors that have been shown to influence gut microbiota colonization in early life include the mode of delivery of the newborn, infant feeding pattern, and household factors such as sibship size [8, 1012]. Additionally, medication such as the use of antibiotics may also influence the pattern of intestinal microbiota colonization [10, 11]. Across geographical locations, socioeconomic and cultural differences would result in a significant variance in the mothers' choice of dietary regimen for their infants, the number of children born within a household (i.e., sibship size) and so on. Therefore, prior to examining the correlation between host health status and gut microbiota, it is essential to better elucidate how the gut microbiota would be affected by the various demographic and lifestyle factors arising from living in different geographic locations.

Our study aimed to investigate the influence of demographic factors on determining the microbial colonization of the infant colon in two Asian populations, Singapore (SG) and Yogyakarta, Indonesia (IN). SG represents an affluent and urbanized community, and IN being an urbanized but developing community. We employed molecular techniques: terminal restriction fragment length polymorphism (T-RFLP) and fluorescent in situ hybridization combined with flow cytometry (FISH-FC) targeting seven major bacterial groups to evaluate and monitor the structure of the colonic microbiota at four time points (i.e, 3 days, one month, three months and one year of age). This study would provide insight on the infant gut microbial succession pattern, as well as the demographic factors that influence stool microbiota signatures in these two Asian populations over the first year of life.


Demographic and Clinical Characteristics

The demographic and clinical characteristics are shown in both Singaporean (SG) and Indonesian (IN) populations (Table 1). Vaginal delivery was more common in SG compared to IN (p = 0.019). In early infancy till 6 months, 85.7% and 80.7% of the SG and IN cohorts, respectively, opted for partial breast and formula feeding. There were a higher percentage of Indonesian infants who were exclusively breastfed in the first 6 months (18.72%, 6/32). In contrast, none of the Singaporean infants were exclusively breastfed for that period of time (p = 0.004). Instead, more SG infants (14.3%) were exclusively formula fed in the first 6 months compared to none in the IN cohort (p = 0.035). Weaning to semisolids for IN cohort occurred later than SG cohort (6.72 months versus 3 months, respectively; p = 0.022). Prenatal antibiotics were administered only in IN cohort (p = 0.013), but the number of infants from both cohorts who received postnatal antibiotics were not significantly different. No significant differences arising from the geographic locations were observed for factors such as gender proportion, postnatal antibiotics consumption and sibling number.

Table 1 Demographic characteristics of Singapore (n = 42) and Indonesia (n = 32) children

Temporal change of relative abundance of seven bacterial groups

The relative abundance of seven bacterial groups was quantified (Figure 1). Although the proportions differed, the trends of bacterial colonization studied over the first year of life were similar for SG and IN cohorts (Figure 1). For example, in both SG and IN cohorts, members of the Enterobacteriaceae family, were one of the earliest colonizers and gradually decreased to an average 0.67% of total bacteria counts at 1 year of age. Colonization of Eubacterium rectale-Clostridium coccoides group increased gradually from 0.18% to 24.07% of total bacteria at 1 year old. The colonization pattern of Bifidobacterium showed an initial increase from a mean of 19.92% at 3 days to 49.50% at 3 months but later decreased to 27.34% at one year of age. A reversal of pattern was seen with Clostridium leptum group where a decrease in colonization from a mean of 5.88% to 1.59% occurred between 3 days and 3 months of age but increased subsequently at the age of one year. The other three bacterial groups such as Bacteroides-Prevotella, Atopobium and Lactobacilli-Enterococci group remained in relatively lower abundance throughout the first year of life, and each constituted less than 10% of the total bacteria detected in stool sample throughout all time points. The phylogenetic gap included the remaining bacterial members that were not targeted by our panel of probes, and the relative abundance of the phylogenetic gap ranged from 22.89% to 37.40% of total bacteria.

Figure 1
figure 1

Comparison of relative abundance of seven predominant bacterial groups between Singapore and Indonesia infants. Singapore cohort is represented by SG while Indonesia cohort is represented by IN. Enterobacteriaceae (Enter), Bifidobacterium (Bif), Clostridium leptum (Clep), Bacteroides-Prevotella (Bac), Lactobacilli-Enterococci (Lab), Atopobium (Ato), Eubacterium rectale-Clostridium coccoides (Erec), and phylogenetic gap (Gap) represent the remaining microbiota that was not targeted by our panel of probes.

Influence of Demographic and Lifestyle/Clinical Characteristics

We comparatively examined the abundances of bacterial groups in relation to demographic factors: geographical origin, mode of delivery, dietary regimen and weaning age, and sibship size. These demographic characteristics differed between the SG and IN cohorts (Table 1), and we aimed to determine if these demographic factors would result in corresponding differences in the abundance of specific fecal associated bacterial groups. The additional file 1 details the univariate analysis of relative abundance (%) of seven bacterial group members of the infant fecal microbiota in relation to geographical and clinical factors.

(A) Geographical Origin

Three bacterial groups, namely Clostridium leptum, Atopobium and Bifidobacterium, differed in abundance between the SG and IN cohorts. Linear mixed model revealed that the relative abundances of Clostridium leptum [coefficient (B): 7.758, 95% confidence interval (CI): 5.063-10.453, adj p < 0.001] and Atopobium group [B: 3.526, 95%CI: 1.102 - 5.949, adj p = 0.005] were higher in IN cohort (Figure 2). Conversely, lower relative abundance of Bifidobacterium was observed in the IN cohort [B: -13.950, 95%CI: -26.423 - -1.476), adj p = 0.029] (Figure 2).

Figure 2
figure 2

Longitudinal comparison of fecal microbiota by geographical origin of subjects. Linear mixed model analysis involved adjustment against other confounding factors (Mode of delivery, weaning age, sibling number, total breastfeeding up to 6 month, eczema and prenatal antibiotics). Indonesia cohort (IN) represented by dotted line, and Singapore cohort (SG) represented by solid line.

(B) Mode of Delivery

Longitudinal analyses over four time points showed that mode of delivery had the largest effect on the abundance of four bacterial groups (Figure 3). In vaginal delivered infants, significantly higher abundance of Bacteroides-Prevotella [B: 3.016, 95%CI: 0.639 - 5.394, adj p = 0.014], Bifidobacterium [B: 16.040, 95%CI: 5.667 - 26.414, adj p = 0.003] and Atopobium group [B: 2.531, 95%CI: 0.472 -4.589, adj p = 0.017] were observed. However, lower abundance of Lactobacilli - Enterococci group [B: -3.665, 95%CI: -5.949 - -1.381, adj p = 0.002] was observed in vaginal delivered infants.

Figure 3
figure 3

Longitudinal comparison of fecal microbiota by mode of delivery. Linear mixed model analysis adjusted with confounding factors (Location, weaning age, sibling number, total breastfeeding up to 6 month, eczema and prenatal antibiotics). Caesarean delivery (LSCS) represented by dotted line, and vaginal delivery (V) represented by solid line.

(C) Breastfeeding and Age of Weaning

The relative abundance of Lactobacilli-Enterococci group was significantly higher in infants who were exclusively breastfed throughout all time points [B: 5.236, 95%CI: 1.044 - 9.428, adj p = 0.015] (Table 2). Another factor that affected the abundance of stool microbiota was the age of weaning to semisolids. Linear mixed model showed a decrease in abundance of Clostridium leptum group for every month of increase in weaning age [B: -0.827, 95%CI: -1.5934 - -0.0602, adj p = 0.035].

Table 2 Feeding habits and demographic factors affecting the relative abundance of microbial groups

(D) Sibship Size

Relative abundance of Bifidobacterium increased by 3.873% with every increase in sibling number [B: 3.873, 95%CI: 0.112 -7.634, adj p = 0.044]. On the other hand, the abundance of Enterobacteriaceae decreased with each increase in sibship size [B: -0.526, 95%CI: -0.8725 - -0.1801, adj p = 0.004] (Table 2).

(E) Exposure to Antibiotics

The relative abundance of Clostridium leptum group at 1 year of age was significantly higher for infants that reported their postnatal antibiotic intake at period of 6 months to 1 year of age [B: 5.706; adj p = 0.025], as compared to the infants who did not consume antibiotics.

Stool Microbial Richness/Diversity

T-RFs of stool microbiota in SG and IN cohorts were obtained from three individual enzymatic digestions (i.e., AluI, MspI and RsaI), and compared for their microbial richness based on Shannon and Simpson Index. Microbial richness between the cohorts was considered different when both Shannon and Simpson Index from all three enzymatic digestions were significantly different. Table 3 shows that there were no observable differences in the microbial richness of SG and IN cohorts at both 3 months and 12 months of age, both before and after adjusting for demographic confounders. In contrast, when the infants from both geographical locations were grouped according to their mode of delivery, microbial richness of stool microbiota in vaginal-delivered infants had a significantly higher microbial richness compared to caesarean-delivered infants at 12 months of age (Table 3). The microbial richness of stool microbiota did not correlate with other lifestyle factors.

Table 3 Shannon and Simpson diversity index determined from T-RFLP profiles


In this study, we first examined the temporal succession of specific bacterial groups in the infant stool microbiota sampled from two geographical locations (i.e., Singapore and Yogyakarta, Indonesia). These two populations in Asia are considered diverse, in terms of geographical distance, cultural ethnicity and overall social economic development. Despite these differences between the two populations, our findings showed that the dynamic of colonization was similar in both cohorts. For example, Enterobacteriaceae and Bifidobacterium constitute the predominant bacterial groups in stool microbiota before three months of age, and were present at a relative abundance of up to 98% of total bacteria. This observation is in agreement with past reports which found that healthy infants from Netherlands, breastfed Indian infants from Guatemala, preterm infants from Nigeria and 6-week old infants across Europe also had a similar predominance of Enterobacteriaceae and Bifidobacterium [10, 1315]. As the infants age, our study also showed that Firmicutes represented by members of the Eubacterium rectale-Clostridium coccoides group increased in its abundance, and gradually resembled that of an adult stool microbiota i.e. mainly populated with members of the Firmicutes and Bacteroidetes phylum [16]. The similarities in the pattern of colonization from early till late infancy despite geographical differences may be related to multiple factors. For example, the prevalence of facultative anaerobes Enterobacteriaceae during early life may be due to a relatively aerobic gastrointestinal tract, and the need for the facultative anaerobes to deplete the oxygen content so as to provide an anoxic environment suitable for other commensal microbes to establish [17].

There remains no clear explanation for the predominance of Bifidobacterium in most infants, including those who were exclusively formula-fed, but not in adults. A possible reason may be related to the diet consumed by the human host at different stages of life. To illustrate, dietary carbohydrates that are consumed by infants comprise mainly disaccharides (lactose) and oligosaccharides [18, 19], which are in turn rapidly hydrolyzed to form galactose and glucose monosaccharides [20, 21]. A portion of these monosaccharides becomes available for the commensal microbiota, and because Bifidobacterium spp. produce more ATP per mole of glucose through the bifidus pathway [22], there remains a selective advantage for Bifidobacterium to out compete the other commensal bacterial groups fermenting carbohydrates through the conventional glycolysis and 6-phosphogluconate pathways. Subsequently, as the host matures and undergoes weaning, the dietary carbohydrates become more complex and eventually favour the establishment of other bacterial members belonging to the Bacteroides and Clostridium for instance, which are known to contain a wide repertoire of polysaccharides-utilizing gene clusters that can effectively degrade complex dietary carbohydrates [2325]. A recent study conducted in five European centres reported a decrease in the abundance of Bifidobacterium, Enterobacteriaceae and certain species of Clostridium, whereas the abundance of members belonging to the Clostridium coccoides and Clostridium leptum groups increased after weaning [26]. Variations in dietary habits between Singapore and Indonesia may explain the differences in rates of colonization of these bacterial groups between Singapore and Indonesian subjects and therefore the slopes of the curves with age for Bifidobacterium, Clostridium leptum and Bacteroides (Figure 2).

A low relative abundance of the Bacteroides-Prevotella group was observed throughout all time points up till the age of 12 month (mean 7.31%). Our previous publication based on 16S rRNA pyrosequencing reported similar proportion of Bacteroides (8.90%) in healthy infants at 12 months [5] and substantiates the findings in this current study. On the contrary in adult the Bacteroidetes co-inhabits with the Firmicutes and both phyla dominate the bacterial community of the human gut microbiome [16, 27, 28]. The structure of the infant gut microbiome is dynamic and evolves over the first years of life toward an adult-like microbiota [2931].

Besides monitoring for the temporal succession of stool microbiota, we further evaluate if demographic and lifestyle differences in the two studied geographical locations (Singapore, SG and Indonesia, IN) would influence the abundance of specific bacterial groups. A study conducted across Europe showed that the geographic origin had an impact on the composition of the gut microbiota [10], and it remains unknown if the structure of the microbiota is influenced to the same extent in Asia. In this study, both SG and IN differ in its extent of development and urbanization, and we observed a higher relative abundance of Bifidobacterium in the SG cohort compared to IN. This might be a common feature of urban populations, as it has also been reported previously for Northern European countries such as Stockholm to have a higher abundance of Bifidobacterium in infants stool microbiota as compared to those sampled in the Spanish province of Granada [10].

In addition, the two geographical locations in this study differ significantly in various aspects, for instance in mode of delivery, feeding history, occurrence of antibiotics consumption and sibling number. Interestingly, these factors studied have also been associated with the development of allergic diseases [3235]. It has been postulated that the influence of these factors have on atopic disease may at least be in part through the effects on profile of gut microbiota. When we examined the effects of demographic and lifestyle factors, we found that the mode of delivery had the largest effect on stool microbiota of infants. These observations are supported by previous studies, where higher numbers of bacterial members belonging to the genus Bifidobacterium [36, 37], Bacteroides and Atopobium group were observed for vaginal delivered infants compared to caesarean delivered infants [8, 10]. Compared to caesarean delivered infants, the general notion is that infants who are vaginally delivered are exposed to their mothers' vaginal microbiota, in particular the Lactobacillus spp. [38], and therefore would have a higher incidence of bacterial transmission to their gut. However, in contrast to previous report which detected a higher abundance of Lactobacillus spp. in vaginally delivered infants [39], we detected a lower abundance of Lactobacilli-Enterococci group in our studied cohort. This discrepancy may be due to the specificities of different oligonucleotide primers/probes used to target the Lactobacillus-Enterococci group. Alternatively, the close adherence of Lactobacillus spp. to mucosal layers might hinder its transmission to the infants while the other vaginal microbiota gets transmitted to the infant [40]. Future validations on a larger cohort of vaginally delivered infants residing in SG and IN will be needed to verify the associated low abundance of Lactobacillus.

Our study also showed that vaginal delivered infants had a significantly higher number of terminal restriction fragments (T-RFs) and microbial richness at 12 months of age. Previous studies had reported that the diversity of stool microbiota increased over time [41]. We postulate that the higher abundance of beneficial bacteria such as Bifidobacterium associated with vaginal delivery may promote the diversity of overall gut microbiota as the infant ages.

Our findings also suggest that antibiotics consumption and sibling number are potential factors that influence the bacterial composition of the human fecal microbiota. For example, the consumption of postnatal antibiotic exposure resulted in a higher relative abundance of members of the Clostridium leptum group at one year of age. Previous studies have also found that postnatal antibiotic intake were associated with decreased numbers of Bifidobacterium and Bacteroides [11, 42], further suggesting that antibiotics consumption can perturb the structure of the commensal microbiota. A higher abundance of Bifidobacterium was observed to be associated with the presence of older siblings [11]. Furthermore, we noted a corresponding decrease in the abundance of Enterobacteriaceae with the number of siblings. Interestingly, Lewis and colleagues have previously reported a decrease in the incidence of allergy with the number of siblings [34], while our past studies have found higher abundance of Bifidobacterium spp. and decreased abundance of Enterobacteriaceae in healthy infants compared to infants with eczema [5, 6]. It remains to be further established if these multitude of factors: the sibship size and abundance of Bifidobacterium spp. and Enterobacteriaceae are intricately linked with the development of allergy and its related disorders. Besides demographic and lifestyle characteristics, the genetic make-up of the host has been proposed to be an important contributing factor in shaping the composition of the gut microbiota. Ambiguous findings were reported for twin studies on the influence of host genetics on gut microbiota [27, 43], whereas absence, mutation or variation of certain single host genes were shown to affect the composition of gut microbiota [44].

FISH-FC approach showed a phylogenetic gap ranging from 22.89% to 37.40% of total bacteria for the four time points. A similar bacterial coverage was reported by Fallani et al using the same method, where the sum of bacterial cells detected were 72.7% ± 24.5% [10] and 74.3% ± 18.9% [45] with a panel of 10 non-overlapping probes. We acknowledge that the molecular techniques applied in this study do not permit a thorough description of the bacterial population inhabiting the human colon. Future studies would aim to utilize deep sequencing of the 16S rRNA genes so as to delve in depth the bacterial communities populating the human microbiome [46, 47]. Their greater depths of sampling offer the opportunity to explore within the phylogenetic gap and beyond, therefore allowing high-resolution association studies involving the bacterial populations of the human microbiome as "quantitative traits".


In conclusion, we have shown that variations in term of relative abundance in infant fecal microbiota are discernable for bacterial groups between two Asian populations of different geographical locations. The differences in the stool microbiota were partly explained by certain lifestyle and clinical factors. These features may confound studies relating to the association of stool microbiota and the predisposition to disease, and should be an important confounder to take note for comparative studies that enrol large population cohort across different geographical origins.


Subject recruitment and study design

The SG at risk of atopy cohort (n = 42) is a subgroup selected from the placebo arm (n = 112) of a randomized double-blind placebo controlled clinical trial on the administration of probiotics supplemented cow's milk-based infant formula for 6 months on the prevention of eczema and allergic diseases. The placebo group of the study received the same cow's milk-based infant formula without probiotics. This study was conducted at National University of Hospital, Singapore ( Identifier: NCT00318695) [48]. The Indonesia at risk of atopy cohort (n = 32) was selected from a birth cohort study (n = 66) recruited from expectant mothers who visited Gadjah Mada University Hospital, Yogyakarta. The inclusion criteria for both cohorts were 1) first-degree relative with a history of allergic disorder as confirmed by a doctor's diagnosis of asthma, allergic rhinitis, or eczema and a positive skin prick test to any of a panel of common dust mite allergens, which are the most important inhalant allergens in our atopic population [49]; 2) gestational age above 35 wk and birth weight above 2 kg; 3) absence of major congenital malformations or major illness at birth; 4) deemed to be in good health based on medical history and physical examination; and 5) the family assessed to be able to complete the trial. Consecutive cases with availability of stool sample with 3 or 4 time points for FISH-FC were selected from both SG and IN cohorts. The sample size was calculated based on study by Sepp et al [50], which reported a higher prevalence of Lactobacillus at 12 months of age in Estonian infants (63%) compared with Swedish infants (38%). We therefore anticipated the difference to be approximately 25% with a power of 90% and a two-sided test size of 5%, 49 subjects were required in each group. No probiotics and prebiotics consumption was reported for SG cohort during the early infancy. Four infants within the IN cohort were partially fed with milk formula that contained prebiotics. Written informed consent for participation in the study was obtained from the parents/guardians of all infants. The study was approved by the National University Hospital's ethics review committee (Ref Code: B/00/322).

Stool sampling

Stool samples were collected on day 3, and at 1, 3, and 12 month after birth based on collection and processing produce as described previously [51]. Stool samples were collected into sterile plastic vials by parents, stored in the freezer at -20°C and delivered to the laboratory within 20 hours. The samples were kept cool on a dry-ice pack during transport and, immediately upon arrival at the laboratory, diluted with 0.85% sodium chloride solution (saline) to give a 0.1 g/ml homogenate. After preparation of the homogenate, samples were fixed in 4% paraformaldehyde (PFA), and stored in TN (10 mM Tris-HCl [pH 8], 150 mM NaCl) buffer at -80°C for later FISH-FC and DNA extraction respectively. Stool samples stored in PFA and TN buffer from Indonesia were shipped on dry-ice pack to single location (laboratory in National University of Singapore) for analysis.

DNA extraction and T-RFLP Analysis

Bacterial DNA extractions from stool samples were carried out as described previously [51] using 0.3 g of 0.1 mm zirconia/silica beads (Biospec, Inc) and a mini bead-beater (Biospec, Inc). In brief, the aqueous supernatant containing DNA will be subsequently subjected to two phenol-chloroform (1:2) extractions and precipitated with 1 ml of ethanol and 50 μl of 3 M sodium acetate. Finally, DNA will be dissolved in 25 μl of sterile TE Buffer (pH 8.0) and stored at -20°C until analysis.

For T-RFLP analysis, 16S rRNA genes were amplified using primers 47F (5'- Cy5 -GCY TAA YAC ATG CAA GT -3') and 1492R (5'-GGY TAC CTT GTT ACG ACT T-3'). PCR reaction mix comprises of 50 ng of DNA template, 0.2 μM (each) of forward and reverse primers, 0.2 mM dNTP, 1.5 U of Ex Taq DNA polymerase (Takara Bio Inc., Japan) and the volume added up to 50 μl with molecular biological grade water. The PCR conditions were as follow: 3 mins at 95°C, 20 cycles (30 seconds at 95°C, 45 seconds at 50°C, 1 min at 72°C), and 5 mins at 72°C. After PCR amplification, mung bean digestion (New England Biolabs [NEB], MA) was carried out to digest single-stranded overhangs. PCR products were purified with QIAquick PCR purification kit (Qiagen, Germany) prior to individual enzymatic digestion by AluI, MspI and RsaI (NEB, MA). Terminal restricted fragments (T-RFs) were analyzed after capillary electrophoresis on CEQ8000 genetic analyzer (Beckman Coulter, CA) [52].

Fluorescence in situhybridization combined with flow cytometry (FISH-FC)

FISH-FC was performed as described previously [53]. A panel of seven bacterial phylogenetic probes was used as described previously by [51]. These probes were selected to target the Eubacterium rectale-Clostridium coccoides group (Erec 482), Clostridium leptum subgroup (Clep 866 and the corresponding competitor probes), Bacteroides-Prevotella group (Bac 303), Bifidobacterium genus (Bif 164), Atopobium group (Ato 291), Lactobacilli-Enterococci group (Lab 158) and Enterobacteriaceae family (Enter 1432). Eubacterial probe EUB338 was used as positive control, while NON 338 probe was used as negative control. Samples were analyzed with FACS Calibur flow cytometer (Becton-Dickinson, USA). A total of 100,000 cells were acquired for analysis per sample and bacterial concentrations adjusted to lower than 3,000 events/s. Subsequent analyses were conducted using the Cell Quest Software (Becton Dickinson, USA). True-positive counts were determined by subtracting double-labelled bacteria with background level evaluated with NON 338 probe. The relative abundance of each bacterial group was expressed as percentage of total EUB338 labelled bacterial cells.

Statistical analysis

All statistical analyses were carried out using SPSS v.18 (SPSS: IBM, Chicago III, USA). Linear mixed model was used to evaluate the demographic effect and time effect (i.e. 4 time points) while adjusting for other confounders. The relative abundance of bacterial groups was used as dependent variable in the model. Three variance-covariance structures (compound symmetry, 1st order autoregressive and unstructured) were used for linear mixed model, and the selection of covariance structure was based on Akaike's Information Criterion (AIC) and Schwarz's Bayesian Criterion (BIC). Linear regression analyse was used to analyze the effects of the postnatal antibiotic consumption on the relative abundance of bacterial groups, with adjustment for confounding factors at single time point. Shannon and Simpson Index were calculated from relative intensity of T-RF as described previously [7], and were used as dependent variable in the model of linear regression to investigate the effects of confounding factors. Confounding factors used for adjustment were based on the all demographics factors being studied in this study and standardized for all statistical analysis (i.e, Linear mixed model and linear regression), those factors were location, mode of delivery, weaning age, sibling number, total breastfeeding up to 6 month, eczema and prenatal antibiotics. Statistical significance was set at p < 0.05. All statistical significance tests and confidence intervals were two-sided.


  1. Kelly D, King T, Aminov R: Importance of microbial colonization of the gut in early life to the development of immunity. Mutat Res. 2007, 622 (1-2): 58-69.

    Article  PubMed  CAS  Google Scholar 

  2. Sekirov I, Russell SL, Antunes LC, Finlay BB: Gut microbiota in health and disease. Physiol Rev. 2010, 90 (3): 859-904. 10.1152/physrev.00045.2009.

    Article  PubMed  CAS  Google Scholar 

  3. Macpherson AJ, Harris NL: Interactions between commensal intestinal bacteria and the immune system. Nat Rev Immunol. 2004, 4 (6): 478-485. 10.1038/nri1373.

    Article  PubMed  CAS  Google Scholar 

  4. Bottcher MF, Nordin EK, Sandin A, Midtvedt T, Bjorksten B: Microflora-associated characteristics in faeces from allergic and nonallergic infants. Clin Exp Allergy. 2000, 30 (11): 1590-1596.

    Article  PubMed  CAS  Google Scholar 

  5. Hong PY, Lee BW, Aw M, Shek LP, Yap GC, Chua KY, Liu WT: Comparative analysis of fecal microbiota in infants with and without eczema. PLoS One. 2010, 5 (4): e9964-10.1371/journal.pone.0009964.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Mah KW, Bjorksten B, Lee BW, van Bever HP, Shek LP, Tan TN, Lee YK, Chua KY: Distinct pattern of commensal gut microbiota in toddlers with eczema. Int Arch Allergy Immunol. 2006, 140 (2): 157-163. 10.1159/000092555.

    Article  PubMed  CAS  Google Scholar 

  7. Wang M, Karlsson C, Olsson C, Adlerberth I, Wold AE, Strachan DP, Martricardi PM, Aberg N, Perkin MR, Tripodi S, et al: Reduced diversity in the early fecal microbiota of infants with atopic eczema. J Allergy Clin Immunol. 2008, 121 (1): 129-134. 10.1016/j.jaci.2007.09.011.

    Article  PubMed  Google Scholar 

  8. Adlerberth I, Strachan DP, Matricardi PM, Ahrne S, Orfei L, Aberg N, Perkin MR, Tripodi S, Hesselmar B, Saalman R, et al: Gut microbiota and development of atopic eczema in 3 European birth cohorts. J Allergy Clin Immunol. 2007, 120 (2): 343-350. 10.1016/j.jaci.2007.05.018.

    Article  PubMed  CAS  Google Scholar 

  9. Bjorksten B, Naaber P, Sepp E, Mikelsaar M: The intestinal microflora in allergic Estonian and Swedish 2-year-old children. Clin Exp Allergy. 1999, 29 (3): 342-346. 10.1046/j.1365-2222.1999.00560.x.

    Article  PubMed  CAS  Google Scholar 

  10. Fallani M, Young D, Scott J, Norin E, Amarri S, Adam R, Aguilera M, Khanna S, Gil A, Edwards CA, et al: Intestinal microbiota of 6-week-old infants across Europe: geographic influence beyond delivery mode, breast-feeding, and antibiotics. J Pediatr Gastroenterol Nutr. 2010, 51 (1): 77-84. 10.1097/MPG.0b013e3181d1b11e.

    Article  PubMed  Google Scholar 

  11. Penders J, Thijs C, Vink C, Stelma FF, Snijders B, Kummeling I, van den Brandt PA, Stobberingh EE: Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics. 2006, 118 (2): 511-521. 10.1542/peds.2005-2824.

    Article  PubMed  Google Scholar 

  12. Gore C, Munro K, Lay C, Bibiloni R, Morris J, Woodcock A, Custovic A, Tannock GW: Bifidobacterium pseudocatenulatum is associated with atopic eczema: a nested case-control study investigating the fecal microbiota of infants. J Allergy Clin Immunol. 2008, 121 (1): 135-140. 10.1016/j.jaci.2007.07.061.

    Article  PubMed  Google Scholar 

  13. Mata LJ, Urrutia JJ: Intestinal Colonization of Breast-Fed Children in a Rural Area of Low Socioeconomic Level. Ann Ny Acad Sci. 1971, 176 (Jan7): 93-

    Article  Google Scholar 

  14. Favier CF, de Vos WM, Akkermans AD: Development of bacterial and bifidobacterial communities in feces of newborn babies. Anaerobe. 2003, 9 (5): 219-229. 10.1016/j.anaerobe.2003.07.001.

    Article  PubMed  Google Scholar 

  15. Rotimi VO, Olowe SA, Ahmed I: The development of bacterial flora of premature neonates. J Hyg (Lond). 1985, 94 (3): 309-318. 10.1017/S0022172400061532.

    Article  CAS  Google Scholar 

  16. Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, Gill SR, Nelson KE, Relman DA: Diversity of the human intestinal microbial flora. Science. 2005, 308 (5728): 1635-1638. 10.1126/science.1110591.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Gaskins HR, Croix JA, Nakamura N, Nava GM: Impact of the intestinal microbiota on the development of mucosal defense. Clin Infect Dis. 2008, 46 (Suppl 2): S80-86. discussion S144-151

    Article  PubMed  Google Scholar 

  18. Ferreira IMPLVO, Gomes AMP, Ferreira MA: Determination of sugars, and some other compounds in infant formulae, follow-up milks and human milk by HPLC-UV/RI. Carbohydrate Polymers. 1998, 37: 225-229. 10.1016/S0144-8617(98)00064-2.

    Article  Google Scholar 

  19. Newburg DS: Oligosaccharides and glycoconjugates in human milk: their role in host defense. Journal of mammary gland biology and neoplasia. 1996, 1 (3): 271-283. 10.1007/BF02018080.

    Article  PubMed  CAS  Google Scholar 

  20. Mobassaleh M, Montgomery RK, Biller JA, Grand RJ: Development of carbohydrate absorption in the fetus and neonate. Pediatrics. 1985, 75 (1 Pt 2): 160-166.

    PubMed  CAS  Google Scholar 

  21. MacLean WC, Fink BB, Schoeller DA, Wong W, Klein PD: Lactose assimilation by full-term infants: relation of [13C] and H2 breath tests with fecal [13C] excretion. Pediatric research. 1983, 17 (8): 629-633.

    Article  PubMed  CAS  Google Scholar 

  22. Palframan RJ, Gibson GR, Rastall RA: Carbohydrate preferences of Bifidobacterium species isolated from the human gut. Current issues in intestinal microbiology. 2003, 4 (2): 71-75.

    PubMed  CAS  Google Scholar 

  23. Xu J, Bjursell MK, Himrod J, Deng S, Carmichael LK, Chiang HC, Hooper LV, Gordon JI: A genomic view of the human-Bacteroides thetaiotaomicron symbiosis. Science (New York, NY. 2003, 299 (5615): 2074-2076. 10.1126/science.1080029.

    Article  CAS  Google Scholar 

  24. Shah HN, Gharbia SE: Ecophysiology and taxonomy of Bacteroides and related taxa. Clin Infect Dis. 1993, 16 (Suppl 4): S160-167.

    Article  PubMed  CAS  Google Scholar 

  25. Pope PB, Denman SE, Jones M, Tringe SG, Barry K, Malfatti SA, McHardy AC, Cheng JF, Hugenholtz P, McSweeney CS, et al: Adaptation to herbivory by the Tammar wallaby includes bacterial and glycoside hydrolase profiles different from other herbivores. Proceedings of the National Academy of Sciences of the United States of America. 107 (33): 14793-14798.

  26. Fallani M, Amarri S, Uusijarvi A, Adam R, Khanna S, Aguilera M, Gil A, Vieites JM, Norin E, Young D, et al: Determinants of the human infant intestinal microbiota after the introduction of first complementary foods in infant samples from five European centres. Microbiology. 2011, 157 (Pt 5): 1385-1392.

    Article  PubMed  CAS  Google Scholar 

  27. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, et al: A core gut microbiome in obese and lean twins. Nature. 2009, 457 (7228): 480-484. 10.1038/nature07540.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Suau A, Bonnet R, Sutren M, Godon JJ, Gibson GR, Collins MD, Dore J: Direct analysis of genes encoding 16S rRNA from complex communities reveals many novel molecular species within the human gut. Appl Environ Microbiol. 1999, 65 (11): 4799-4807.

    PubMed  CAS  PubMed Central  Google Scholar 

  29. Koenig JE, Spor A, Scalfone N, Fricker AD, Stombaugh J, Knight R, Angenent LT, Ley RE: Succession of microbial consortia in the developing infant gut microbiome. Proc Natl Acad Sci USA. 2011, 108 (Suppl 1): 4578-4585.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  30. Favier CF, Vaughan EE, De Vos WM, Akkermans AD: Molecular monitoring of succession of bacterial communities in human neonates. Appl Environ Microbiol. 2002, 68 (1): 219-226. 10.1128/AEM.68.1.219-226.2002.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  31. Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO: Development of the human infant intestinal microbiota. PLoS Biol. 2007, 5 (7): e177-10.1371/journal.pbio.0050177.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Bager P, Wohlfahrt J, Westergaard T: Caesarean delivery and risk of atopy and allergic disease: meta-analyses. Clin Exp Allergy. 2008, 38 (4): 634-642. 10.1111/j.1365-2222.2008.02939.x.

    Article  PubMed  CAS  Google Scholar 

  33. Kummeling I, Stelma FF, Dagnelie PC, Snijders BE, Penders J, Huber M, van Ree R, van den Brandt PA, Thijs C: Early life exposure to antibiotics and the subsequent development of eczema, wheeze, and allergic sensitization in the first 2 years of life: the KOALA Birth Cohort Study. Pediatrics. 2007, 119 (1): e225-231. 10.1542/peds.2006-0896.

    Article  PubMed  Google Scholar 

  34. Lewis SA, Britton JR: Consistent effects of high socioeconomic status and low birth order, and the modifying effect of maternal smoking on the risk of allergic disease during childhood. Respir Med. 1998, 92 (10): 1237-1244. 10.1016/S0954-6111(98)90427-9.

    Article  PubMed  CAS  Google Scholar 

  35. Kull I, Bohme M, Wahlgren CF, Nordvall L, Pershagen G, Wickman M: Breast-feeding reduces the risk for childhood eczema. J Allergy Clin Immunol. 2005, 116 (3): 657-661. 10.1016/j.jaci.2005.04.028.

    Article  PubMed  Google Scholar 

  36. Biasucci G, Rubini M, Riboni S, Morelli L, Bessi E, Retetangos C: Mode of delivery affects the bacterial community in the newborn gut. Early Hum Dev. 2010, 86 (Suppl 1): 13-15.

    Article  PubMed  Google Scholar 

  37. Huurre A, Kalliomaki M, Rautava S, Rinne M, Salminen S, Isolauri E: Mode of delivery - effects on gut microbiota and humoral immunity. Neonatology. 2008, 93 (4): 236-240. 10.1159/000111102.

    Article  PubMed  Google Scholar 

  38. Zhou X, Bent SJ, Schneider MG, Davis CC, Islam MR, Forney LJ: Characterization of vaginal microbial communities in adult healthy women using cultivation-independent methods. Microbiology (Reading, England). 2004, 150 (Pt 8): 2565-2573.

    Article  CAS  Google Scholar 

  39. Matsumiya Y, Kato N, Watanabe K, Kato H: Molecular epidemiological study of vertical transmission of vaginal Lactobacillus species from mothers to newborn infants in Japanese, by arbitrarily primed polymerase chain reaction. J Infect Chemother. 2002, 8 (1): 43-49.

    Article  PubMed  CAS  Google Scholar 

  40. Velez MP, De Keersmaecker SC, Vanderleyden J: Adherence factors of Lactobacillus in the human gastrointestinal tract. FEMS microbiology letters. 2007, 276 (2): 140-148. 10.1111/j.1574-6968.2007.00908.x.

    Article  PubMed  CAS  Google Scholar 

  41. Roger LC, McCartney AL: Longitudinal investigation of the faecal microbiota of healthy full-term infants using fluorescence in situ hybridization and denaturing gradient gel electrophoresis. Microbiology. 2010, 156 (Pt 11): 3317-3328.

    Article  PubMed  CAS  Google Scholar 

  42. Dethlefsen L, Huse S, Sogin ML, Relman DA: The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing. PLoS biology. 2008, 6 (11): e280-10.1371/journal.pbio.0060280.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Stewart JA, Chadwick VS, Murray A: Investigations into the influence of host genetics on the predominant eubacteria in the faecal microflora of children. J Med Microbiol. 2005, 54 (Pt 12): 1239-1242.

    Article  PubMed  CAS  Google Scholar 

  44. Spor A, Koren O, Ley R: Unravelling the effects of the environment and host genotype on the gut microbiome. Nat Rev Microbiol. 2011, 9 (4): 279-290. 10.1038/nrmicro2540.

    Article  PubMed  CAS  Google Scholar 

  45. Fallani M, Amarri S, Uusijarvi A, Adam R, Khanna S, Aguilera M, Gil A, Vieites JM, Norin E, Young D, et al: Determinants of the human infant intestinal microbiota after the introduction of first complementary foods in infant samples from five European centres. Microbiology. 2011, 157 (Pt 5): 1385-1392.

    Article  PubMed  CAS  Google Scholar 

  46. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R: Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci USA. 2011, 108 (Suppl 1): 4516-4522.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  47. Huse SM, Dethlefsen L, Huber JA, Mark Welch D, Relman DA, Sogin ML: Exploring microbial diversity and taxonomy using SSU rRNA hypervariable tag sequencing. PLoS Genet. 2008, 4 (11): e1000255-10.1371/journal.pgen.1000255.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Soh SE, Aw M, Gerez I, Chong YS, Rauff M, Ng YP, Wong HB, Pai N, Lee BW, Shek LP: Probiotic supplementation in the first 6 months of life in at risk Asian infants--effects on eczema and atopic sensitization at the age of 1 year. Clin Exp Allergy. 2009, 39 (4): 571-578. 10.1111/j.1365-2222.2008.03133.x.

    Article  PubMed  CAS  Google Scholar 

  49. Chew FT, Lim SH, Goh DY, Lee BW: Sensitization to local dust-mite fauna in Singapore. Allergy. 1999, 54 (11): 1150-1159. 10.1034/j.1398-9995.1999.00050.x.

    Article  PubMed  CAS  Google Scholar 

  50. Sepp E, Julge K, Vasar M, Naaber P, Bjorksten B, Mikelsaar M: Intestinal microflora of Estonian and Swedish infants. Acta Paediatr. 1997, 86 (9): 956-961. 10.1111/j.1651-2227.1997.tb15178.x.

    Article  PubMed  CAS  Google Scholar 

  51. Mah KW, Chin VI, Wong WS, Lay C, Tannock GW, Shek LP, Aw MM, Chua KY, Wong HB, Panchalingham A, et al: Effect of a milk formula containing probiotics on the fecal microbiota of asian infants at risk of atopic diseases. Pediatr Res. 2007, 62 (6): 674-679. 10.1203/PDR.0b013e31815991d5.

    Article  PubMed  Google Scholar 

  52. Liu WT, Marsh TL, Cheng H, Forney LJ: Characterization of microbial diversity by determining terminal restriction fragment length polymorphisms of genes encoding 16S rRNA. Appl Environ Microbiol. 1997, 63 (11): 4516-4522.

    PubMed  CAS  PubMed Central  Google Scholar 

  53. Lay C, Sutren M, Rochet V, Saunier K, Dore J, Rigottier-Gois L: Design and validation of 16S rRNA probes to enumerate members of the Clostridium leptum subgroup in human faecal microbiota. Environ Microbiol. 2005, 7 (7): 933-946. 10.1111/j.1462-2920.2005.00763.x.

    Article  PubMed  CAS  Google Scholar 

Download references


The authors would like to thank Dr Shen Liang for statistical advice, Drs Dawn Lim and Irvin Gerez, Ms Hor Chuen Yee, Corinne Kwek Poh Lian and Judy Anthony who assisted with follow-up of the subjects. The voluntary participation of all subjects in this study is sincerely appreciated. This study was supported by A*STAR's Biomedical Research Council (BMRC) and the MOH's National Medical Research Council (BMRC/08/1/21/19/566).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Bee Wah Lee.

Additional information

Authors' contributions

GCY and KKC performed the experiments, data analysis and statistical analysis. GCY drafted the manuscript. PYH and BWL helped to revise the manuscript. CL helped in experimental techniques for FISH-FC. YDZ and DL participated in collation of clinical data and helped in statistical analysis. MA, LPCS and BWL conceived the study. DC, S, YS, MA, LPCS, KYC and BWL participated in study design and helped in coordination of sample and data collection. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Univariate analysis of relative abundance of seven predominant bacterial groups. Univariate analysis of relative abundance of seven predominant bacterial groups were performed for location, mode of delivery, total breastfeeding up to 6 month, eczema, prenatal antibiotics and postnatal antibiotics. Statistical significance were bold formatted (p value < 0.05). (XLS 50 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Yap, G.C., Chee, K.K., Hong, PY. et al. Evaluation of stool microbiota signatures in two cohorts of Asian (Singapore and Indonesia) newborns at risk of atopy. BMC Microbiol 11, 193 (2011).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: