BACKGROUND
major microbial colonization of the gastrointestinal tract starts at delivery when an infant comes into contact with the environment.
the composition of developing microbiota is affected by factors such as mode of delivery  <cit> , dietary pattern  <cit>  and administration of probiotics or antibiotics  <cit> .
the early colonization events and the commensal intestinal microbiota shape the immune system and potentially affect the development of variety of diseases  <cit> .
previous studies have shown associations between the composition of intestinal microbiota and atopic diseases.
most of these have addressed the microbiota composition preceding the development of atopic disease, while microbiota aberrancies in infants already suffering from eczema have obtained less attention.
reduced diversity at early life  has been associated with an increased risk of developing atopic disease  <cit> .
the results on specific bacterial species or groups that would either increase or decrease the risk of developing allergy are still conflicting  <cit> .
few studies have observed microbiota alterations in allergic children  with also conflicting results  <cit> .
for example, faecal bifidobacterial counts have been reported to be both decreased  <cit>  or similar  <cit>  as compared to healthy children.
similarly, microbiota diversity in allergic children was found to be decreased in one study  <cit>  but not in another  <cit> .
notably, the studies were conducted by using traditional cultivation-based techniques or molecular techniques that target a sub-set of the intestinal microbiota.
finally, the administration of specific probiotic bacteria during pregnancy and/or during the first months of life has been shown to reduce the risk of atopy, especially atopic eczema  <cit> .
however, some studies have failed to find any connection between the microbiota composition and development of atopic eczema  <cit>  or to confirm the role of probiotics in prevention of atopic diseases  <cit> .
a variety of high-throughput methods based on 16s ribosomal rna  gene sequence analysis have been established to analyse the intestinal microbiota in a culture-independent way, including next -generation sequencing analysis and phylogenetic microarrays  <cit> .
the high-density phylogenetic microarray hitchip consists of  <dig> unique 16s rrna gene targeting oligonucleotide probes that selectively recognise microbes at different taxonomic levels  <cit> .
this and other microarrays have shown to be instrumental for the comprehensive and high-resolution analysis of the microbiota composition from microbial species  to phylum-like level  <cit> .
the objective of this study was to characterize the diversity and temporal changes of intestinal microbiota in early childhood and to identify specific bacterial groups associated with eczema.
by using the hitchip microarray and strategic qpcr analysis of early life fecal samples, we detected specific differences in microbiota composition between healthy children and those with eczema.
methods
study design, subjects and faecal samples
subjects of this study represent a sub-population from a prospective follow-up trial at turku university central hospital, finland, which has been described in detail previously  <cit> .
briefly, the inclusion criterion for the children was that they had a high risk of atopic diseases, i.e.
they had at least one close relative  with atopic eczema, allergic rhinitis or asthma.
further inclusion criteria for present study were vaginal delivery after full-term pregnancy , normal birth weight  and the availability of faecal samples taken at the ages of  <dig> and/or  <dig> months.
finally, all infants were exclusively or partially breast-fed for at least four months.
based on these criteria,  <dig> children from the original study population   <cit>  were included in this study.
the basic characteristics of the study subjects are shown in additional file  <dig>
mothers were randomized to receive capsules containing either placebo or  <dig> ×  <dig> colony-forming units of lactobacillus rhamnosus gg  daily for 2– <dig> weeks before expected delivery.
the intervention continued  <dig> months postnatally.
the capsule contents were consumed by mothers during the exclusive breastfeeding, otherwise infants received the agents.
the occurrence of eczema was diagnosed by the age of  <dig> years by typical skin lesions found in children and chronic relapsing course.
this last criterion was fulfilled if the child had had eczema for  <dig> month or longer at the 24-month study visit and on at least one previous visit .
eczema was considered atopic if it was associated with positive skin prick test at  <dig> and/or  <dig> -month study visit.
none of the study subjects included in present study suffered from asthma or allergic rhinitis.
also, all the infants were normal weight at the age of  <dig> and  <dig> months of age.
the study protocol was approved by the ethics committee of the hospital district of southwest finland and subjects were enrolled in the study after written informed consent was obtained.
faecal samples and dna extraction
the faecal samples were taken from children at age of  <dig> and  <dig> months.
the samples were aliquoted and frozen immediately after collection, and stored in −80°c.
dna was extracted from faecal samples using the repeated bead-beating method as described previously  <cit> .
16s rrna gene microarray analysis
the composition of total microbiota was assessed by using the phylogenetic human intestinal tract chip  as described previously  <cit> , except for the amplification step, where  <dig> cycles of end-point pcr were used.
microarray analysis of all samples were performed in at least two independent hybridizations until satisfactory reproducibility was achieved .
this study reports results of more than  <dig> independent microarray hybridizations.
the hitchip is a custom-made agilent microarray  designed to comprehensively cover the diversity of the human intestinal microbiota.
the array contains  <dig> unique oligonucleotide probes targeting the v <dig> and v <dig> hypervariable regions of the 16s rrna gene and covering over  <dig> intestinal bacterial phylotypes.
the hitchip allows the analysis at three phylogenetic levels: phylum-like level , genus-like level  and phylotype level .
the details of the hitchip have previously been described, including its validation for phylogenetic fingerprinting and quantification  <cit> .
microarray data extraction and microbiota diversity assessment
data were extracted from microarray images using the agilent feature extraction software, version  <dig> .
<dig> .
normalization of microarray data was performed as described earlier  <cit> .
further data processing was performed by using a custom designed relational database running under the mysql database management system  using r-based scripts  <cit> .
quantitative pcr
quantitative pcr  analysis of bifidobacterium genus and species was carried out in an applied biosystems  <dig> fast real-time pcr system in a 96-well format and by using sybr green chemistry .
the primers and their specificities are presented in additional file  <dig>  the pcr reactions and thermocycling conditions were as reported earlier  <cit> .
standards for qpcr were prepared as described in nermes et al.
<cit> .
samples were analysed in duplicate in at least two independent runs.
statistical and data analyses
statistical analysis of both qpcr and hitchip data was carried out with log-transformed data.
in qpcr data, non-detected values were imputed with the half of the theoretical detection limit.
for hitchip data, linear models with factors for treatment, health status, time point and breast-feeding with subsequent anova and contrast tests were used to determine the statistical differences between groups.
in microarray data, cut-off values for positive responding probes were calculated as described before  <cit> .
in hitchip data the analysed values were summary values on phylum-like and genus-like level, obtained by summing the intensities from all the probes assigned to the respective phylum-like or genus-like phylogetic groups.
totally  <dig> phylum-like and  <dig> genus-like level groups reached the detection threshold and were thus used in statistical analysis.
the data is presented as mean with standard deviation values.
redundancy analysis  was performed by using the multivariate statistical analysis package canoco  <cit> .
rda plot shows bacterial groups principally contributing to the difference between the groups of subjects.
the significance of separation in rda was assessed by monte carlo permutation procedure .
the diversity of the microbial community assessed by hitchip was expressed as simpson’s reciprocal index of diversity  as described before  <cit> .
RESULTS
temporal development of microbiota
the faecal microbiota of  <dig> children at age of  <dig> and  <dig> months was analysed using the hitchip phylogenetic microarray.
the diversity of total microbiota increased significantly with age, as the simpson’s the reciprocal diversity index has changed from  <dig> ±  <dig> to  <dig> ±  <dig> at age of  <dig> and  <dig> months, respectively .
at the phylum-like level, significant changes in the relative abundances of major bacterial groups were detected .
the most prominent decline in abundance was observed for actinobacteria that contributed  <dig> % and  <dig> % to the total signal at  <dig> and  <dig> months of age, respectively .
signal intensities for actinobacteria were almost entirely obtained from bifidobacteria .
this finding was consistent with quantitative pcr analysis, where total bifidobacteria counts decreased significantly with age .
at the species level, the amounts of b. longum/infantis group, b. breve, b. bifidum, b. catenulatum group and b. adolescentis decreased over time as assessed by qpcr.
in addition to actinobacteria, the relative abundance of bacilli decreased with age .
all genus-like groups belonging to bacilli decreased, most of which not significantly as individual groups, but the sum effect at the phylum-like level was significant .
an opposite trend was observed for the bacterial groups belonging to the major clostridium clusters, especially members of the clostridium cluster xiva which relative abundance increased with age from  <dig> % to  <dig> % .similar was observed for ruminococcus bromii et rel.
group from clostridium cluster iv that increased from  <dig> % to  <dig> % .
in total,  <dig> genus-like phylogenetic groups changed significantly with age, , which further highlights the extensive compositional changes that the microbiota is undergoing during this period of life.
genus-like groups with a p-value less than  <dig>  are presented in the table.
analysis of the intestinal microbiota composition in relation to the health status
when comparing the microbiota of the two groups of children at the age of  <dig> months, pronounced differences were observed both in the microbial composition and the diversity.
infants with eczema had a significantly more diverse total microbiota .
analysis at the species-like level showed that a large number of bacterial species have different abundance between healthy and eczematous infants, although the individual p-values are not particularly small .
the numerous, but mostly not significant, differences at the species-like level prompted us to look at the trends in microbiota differences at higher levels i.e.
at the phylum-like and genus-like levels.
analysis at the phylum-like level showed that the most remarkable difference between the groups was observed in the bacteroidetes, which were 3-fold more abundant in healthy children at the age of  <dig> months .
the effect of bacteroidetes is also clearly shown in the rda plot , which reveals the bacterial groups principally contributing to the difference among the groups of subjects.
the microbiota differences between the health groups shown in the rda are significant as assessed by mcpp  and a total of  <dig> % of the variation within the dataset could be related to the health status of the infants.
in contrast to the bacteroidetes, specific bacterial groups from the most abundant groups of the firmicutes phylum - clostridium clusters iv and xiva - were significantly more abundant in children with eczema .
in summary, the multiple differences in specific bacterial groups result in microbiota profiles that are significantly distinct between healthy and eczematous infants as assessed by mcpp .
* % of total hitchip signal
nd, below the detection level.
similar microbiota changes were observed between children with eczema and healthy children when only subjects in the placebo group were analysed, although the differences between the groups did not reach statistical difference due to the small number of subjects .
no significant differences were observed for the bifidobacterial sub-community between the two groups of children using both hitchip and qpcr analyses .
the comprehensive list of phylum-like and genus-like level data and p-values obtained by statistical analyses are presented in additional file  <dig> and additional file  <dig>  respectively.
notably, an indication towards altered microbiota composition in children with eczema was already identified at  <dig> months, although the difference did not reach the level of statistical significance .
a higher abundance of the clostridium cluster xiva bacteria was observed in infants with eczema than healthy controls .
l. rhamnosus gg supplementation in early infancy has minor long-term effects on the microbiota composition
when comparing the levels of hitchip signals between children from the placebo group and those who had received l. rhamnosus gg for their first  <dig> months of life, no statistically significant differences were observed at the age of  <dig> months.
however, the supplementation with l. rhamnosus gg showed effects on three genus-like bacterial groups at the age of  <dig> months i.e.
a year after the cessation of the probiotic supplementation.
the children that had received l. rhamnosus gg had higher levels of the butyrate-producing groups anaerostipes caccae et rel  and eubacterium ventriosum et rel  than those of placebo group .
moreover, the placebo group children had higher levels of clostridium difficile et rel at  <dig> months of age as compared to the lgg group children .
the comprehensive list of phylum-like and genus-like level data and p-values obtained by statistical analyses are presented in additional file  <dig> and additional file  <dig>  respectively.
the effect of the probiotic supplementation on the microbiota composition within the group of healthy children or the group of children with eczema was not addressed due to the small number of subjects.
discussion
we used a high-throughput phylogenetic microarray to reveal alterations in the gut microbiota composition throughout early childhood.
the used microarray has been developed and validated for determining the microbiota diversity and evaluating the relative proportions of genus-like or higher  phylogenetic groups  <cit> .
moreover, it has been demonstrated that this microarray has a power equal to or higher than deep new generation sequencing  <cit>  and it has previously been used to compare the microbiota diversity and composition in a variety of studies  <cit> .
despite the highly significant increase of microbiota diversity with age, the diversity indeces at  <dig> months of age are still relatively low  when compared to the approximately two-fold higher indexes  commonly observed in healthy adults  <cit> .
it has been suggested that by the age of  <dig> to  <dig> years the microbiota resembles that of an adult  <cit> .
our results show that microbiota succession continues at least until the age of  <dig> months and most likely even further, because the bacterial diversity has still not reached the diversity of an adult person.
thus, significant changes can be expected to occur in even after  <dig> months of age.
concerning the microbiota composition at  <dig> months of age, our results are in agreement with earlier studies  <cit> , except that we observed significant colonization by bifidobacteria in most of the children  while in the study of palmer et al.
<cit>  bifidobacteria were not detected, possibly due to differences in dna extraction, pcr primers, demographic and geographic origin, dietary patterns of the infants or other confounding factors.
primers used for pcr are often not so optimal for bifidobacteria than for other species and thus, high gc bacteria may perform less well in such pcrs.
further, in our previous studies we have shown that mechanical lysis of faecal bacteria is essential and improves the detection of especially gram-positive bacteria including bifidobacteria  <cit> .
in the palmer et al.
study  <cit> , mechanical lysis by bead-beating was not applied, which may have hampered the detection of bifidobacteria.
thus, we consider that the most likely explanation for the different results concerning bifidobacteria in our and palmer et al.
<cit>  study is the different dna extraction methods used.
when comparing healthy and eczematous children we found statistically significant differences in microbiota composition only at  <dig> months of age.
the total microbiota of children with eczema was found to become significantly more diverse than the microbiota of children who remained healthy by  <dig> months of age.
interestingly, the total microbiota and particularly firmicutes diversity was higher in the eczema group children, although the difference with the healthy subjects was not statistically significant.
abrahamsson et al.
described the infants as having atopic eczema during the first two years of life , but the age at the onset of symptoms was not clarified  <cit> .
however, it can be concluded from the abrahamsson et al.
<cit>  and our study, both taking advantage of the high resolution microbiota assessment techniques, that after  <dig> year of age the total microbiota diversity in children either developing or having eczema is comparable or even higher than that of healthy children.
secretory iga has been suggested to play a role in shaping the microbiota composition and diversity.
some early studies showed an association between the low levels of secretory iga and the risk of developing atopy  <cit>  and could suggest that the low iga levels permit establishment of a wider variety of bacteria and explain the higher bacterial diversity in children with eczema observed in this study.
however, more recent studies have shown a higher concentration of secretory iga in children with allergic sensitization during the first  <dig> years of life  <cit> .
another possible explanation for the increased bacterial diversity in children with eczema is the decreased levels or altered repertoire of antimicrobial peptides secreted into the gut lumen.
these peptides, such as alpha- and beta-defensins, have at least two key roles at the mucosal interface: contributing to the host defense against enteric bacterial attachment and homeostatic control of the intestinal bacterial ecosystem  <cit> .
recently, decreased alpha-defensin levels and increased beta-defensin levels were associated with increased risk of developing atopy  <cit> .
to our knowledge, the levels of faecal antimicrobial peptides in children already having eczema have not been studied.
however, a few studies have highlighted the role of alpha-defensins in microbiota composition and intestinal health.
for example, genetic mutations resulting in decreased alpha-defensin expression have been associated with the susceptibility and severity of inflammatory bowel disease in humans and decreased alpha-defensins may have an effect on the differences observed in microbiota composition between healthy and diseased subjects  <cit> .
interestingly, mice deficient in production of active alpha -defensins were shown to have a decrease in bacteroidetes  <cit> .
the reason for decreased bacteroidetes levels in children with eczema in this study remains unaccountable, but alpha-defensins provide one possible explanation for our observation.
also other host-dependent factors, such as the amount of mucus secretion and differences in mucus glycosylation  may have an influence on the microbiota diversity and composition, as recently reviewed by maynard et al.
<cit> .
clearly, the role of intestinal iga levels, antimicrobial peptides and mucus secretion in shaping the gut microbiota in healthy and eczematous children warrants for further investigation.
our results emphasize that the microbiota diversity in children with eczema should be further studied by using high-resolution techniques in order to define the favourable course of bacterial succession in early childhood and toddler age and to evaluate possible means to influence it.
it was observed that children with eczema harbour more bacteria belonging to the clostridium cluster iv and clostridium cluster xiva.
these bacteria are among the most abundant microbial groups detected in the healthy adult intestine  <cit> .
thus, prematurely occurring changes towards an adult-type microbiota were observed to take place in children with eczema.
it has been suggested that resident bacteria may shape the hosts’ physiology, among others, by modulating the expression of genes involved in intestinal functions, such as postnatal intestinal maturation and the maintenance of mucosal barrier  <cit> .
it may be speculated that an infant-type microbiota supports adequate gut barrier function and tolerance against food allergens in an immature gut.
infant-type microbiota may fortify the normal mucosal barrier function e.g.
by affecting the maturation of the gut epithelium and immune functions in an optimal way and decrease the low-grade intestinal inflammation observable in subjects with eczema  <cit> .
maintenance of adequate mucosal barrier function may also play a role in the level of sensitisation to food-derived compounds  <cit> .
the complex host-microbe interactions in the intestinal epithelium are only recently beginning to be understood  <cit> .
furthermore, we observed decreased relative abundances of bacteria belonging to bacteroidetes in children with eczema.
previous studies have reported an association between decreased amounts bacteroides spp.
and the development of atopy and increased risk for atopic sensitization  <cit> .
bacteria belonging to the bacteroidetes are among the first groups colonizing the gut  <cit>  and they are typical intestinal habitants in healthy adults  <cit> .
bacteroides spp.
are specialized in the breakdown of complex plant polysaccharides  <cit>  and their abundance has been associated with increased short-chain fatty acid concentrations in the infant gut after introduction of first solid foods  <cit> .
furthermore, b. fragilis polysaccharide has been shown in mice model to direct the cellular and physical maturation of the developing immune system via its ability to direct the development of cd4+ t cells, thus inducing the differentiation of th <dig> lineage and correction of the th1/th <dig> imbalance  <cit> .
together with our findings, these results suggest the significance of bacteroides spp.
in the development and maintenance of healthy infant gut and balanced mucosal immunity and necessitate the role of these bacteria to be considered in future studies.
when comparing healthy children with children with eczema we found statistically significant differences in microbiota composition only at  <dig> months, but not at  <dig> months of age.
breast-feeding is known as a major factor influencing the microbiota composition in infancy  <cit> .
at  <dig> months of age, the majority of children included in this study were still nursed and breast-feeding is likely to have had a strong influence on their microbiota composition at that age.
thus, it seems that breast-feeding could have evened up the microbiota differences between the healthy and eczematous children and masked the eczema-associated changes, which came apparent and measurable at  <dig> months of age after the withdrawal of breast-milk.
while many studies addressed the impact of l. rhamnosus gg on health parameters, the short and long-term effect on the intestinal microbiota has only received limited attention.
in the present intervention, the supplementation of l. rhamnosus gg continued until the age of  <dig> months.
interestingly, no significant effect on the microbiota composition was observed at the age of  <dig> months, but instead the supplementation of l. rhamnosus gg in early life was observed to a induce long-term effect and small but significant changes between the intervention groups were observed one year later at the age of  <dig> months.
the observation that the c. difficile et rel.
group bacteria were lower in the lgg groups as compared to placebo is of particular interest.
previously, clostridium difficile colonization at the age of  <dig> month has been associated with a higher risk of a diagnosis of atopic dermatitis at the age of  <dig> years  <cit> .
the higher anaerostipes caccae et rel levels in the children that had received the l. rhamnosus gg supplementation is also a potentially beneficial effect, because a. caccae produces butyrate, which is an energy source for epithelial cells of colonic mucosa  <cit> .
bacteria belonging to the eubacterium ventriosum et rel group that were higher in the children that received the probiotic supplementation, also have shown to produce butyrate but have been less investigated.
in mice, however, it has been shown that e. ventriosum was reduced in colitic mice as compared to non-colitic animals  <cit> .
to our knowledge this is the first high -throughput microbiota analysis study reporting the long-term effects of a probiotic strain on the microbiota composition in early life.
CONCLUSIONS
in conclusion, using a comprehensive microbial analysis approach we observed children with eczema to harbour a more diverse total microbiota and detected specific shifts in bacterial groups in different phylogenetic levels.
the results indicate that aberrancies in microbiota composition are associated with eczema.
our results also suggest that in children at high-risk for atopic disease, a diverse adult-type microbiota in too early childhood may be a potential risk factor and further strengthen the importance of early microbiota characterization and potential dietary modification.
abbreviations
et rel.
: and relatives; hitchip: human intestinal tract chip; mcpp: monte carlo permutation procedure; qpcr: quantitative real-time polymerase chain reaction; rrna: ribosomal rna.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
mk and es designed the original intervention study and organized the sample collection.
infants were clinically examined by mk.
ln, wmdv, rs and ss designed the current study.
ln performed faecal microbial dna extraction, qpcr analyses and hitchip experiments.
mr-s was involved in hitchip experiments.
jn performed bioinformatic analyses.
ln, rs and wmdv interpreted the results and wrote the paper.
all authors read and approved the final manuscript.
supplementary material
additional file 1
basic characteristics of the study subjects.
click here for file
additional file 2
primers targeting bifidobacterium genus and species used in this study.
click here for file
additional file 3
differences in bifidobacterial composition of all children at  <dig> and  <dig> months of age as assessed by using quantitative pcr.
click here for file
additional file 4
the histograms showing the distribution of p-values obtained from the statistical analyses of species-like level of hitchip data at  <dig> months.
each bar represents how many species-like groups gave a p-value in the given range when the effect of different factors on microbiota composition were analysed.
click here for file
additional file 5
the microbiota differences of healthy and eczematous children from placebo group as assessed by hitchip analysis.
click here for file
additional file 6
bifidobacterial sub-communities in infants with eczema and healthy controls as assessed by quantitative pcr and hitchip analyses.
click here for file
additional file 7
phylum-like  and genus-like  hitchip data used in this study.
data is presented as log-transformed values.
a letter a refers to  <dig> months samples and a letter d to  <dig> months samples, respectively.
click here for file
additional file 8
p-values obtained from the statistical analysis of phylum-like and genus-like groups of hitchip data at  <dig> months.
p-values are not corrected and therefore indicate trend-like differences in the abundance of individual bacterial groups between the groups of infants.
microbial groups that were over the detection level were included in the analysis.
click here for file
additional file 9
the microbiota differences between the intervention groups  at the age of  <dig> months as assessed by hitchip analysis.
click here for file
acknowledgements
this work was funded by finnish funding agency for technology and innovation .
in addition, the academy of finland is acknowledged for financial support .
hans heilig, outi immonen and alla kaljukivi are thanked for their excellent technical assistance.
we thank professor airi palva for valuable discussions and her support to carry out this study.

 major microbial colonization of the gastrointestinal tract starts at delivery when an infant comes into contact with the environment. the composition of developing microbiota is affected by factors such as mode of delivery, dietary pattern and administration of probiotics or antibiotics. the early colonization events and the commensal intestinal microbiota shape the immune system and potentially affect the development of variety of diseases. previous studies have shown associations between the composition of intestinal microbiota and atopic diseases. most of these have addressed the microbiota composition preceding the development of atopic disease, while microbiota aberrancies in infants already suffering from eczema have obtained less attention. reduced diversity at early life has been associated with an increased risk of developing atopic disease. the results on specific bacterial species or groups that would either increase or decrease the risk of developing allergy are still conflicting. few studies have observed microbiota alterations in allergic children with also conflicting results. for example, faecal bifidobacterial counts have been reported to be both decreased or similar as compared to healthy children. similarly, microbiota diversity in allergic children was found to be decreased in one study but not in another. notably, the studies were conducted by using traditional cultivation-based techniques or molecular techniques that target a sub-set of the intestinal microbiota. finally, the administration of specific probiotic bacteria during pregnancy and/or during the first months of life has been shown to reduce the risk of atopy, especially atopic eczema. however, some studies have failed to find any connection between the microbiota composition and development of atopic eczema or to confirm the role of probiotics in prevention of atopic diseases. a variety of high-throughput methods based on 16s ribosomal rna gene sequence analysis have been established to analyse the intestinal microbiota in a culture-independent way, including next -generation sequencing analysis and phylogenetic microarrays. the high-density phylogenetic microarray hitchip consists of unique 16s rrna gene targeting oligonucleotide probes that selectively recognise microbes at different taxonomic levels. this and other microarrays have shown to be instrumental for the comprehensive and high-resolution analysis of the microbiota composition from microbial species to phylum-like level. the objective of this study was to characterize the diversity and temporal changes of intestinal microbiota in early childhood and to identify specific bacterial groups associated with eczema. by using the hitchip microarray and strategic qpcr analysis of early life fecal samples, we detected specific differences in microbiota composition between healthy children and those with eczema. methods study design, subjects and faecal samples subjects of this study represent a sub-population from a prospective follow-up trial at turku university central hospital, finland, which has been described in detail previously. briefly, the inclusion criterion for the children was that they had a high risk of atopic diseases, i.e. they had at least one close relative with atopic eczema, allergic rhinitis or asthma. further inclusion criteria for present study were vaginal delivery after full-term pregnancy, normal birth weight and the availability of faecal samples taken at the ages of and/or months. finally, all infants were exclusively or partially breast-fed for at least four months. based on these criteria, children from the original study population were included in this study. the basic characteristics of the study subjects are shown in additional file mothers were randomized to receive capsules containing either placebo or × colony-forming units of lactobacillus rhamnosus gg daily for 2– weeks before expected delivery. the intervention continued months postnatally. the capsule contents were consumed by mothers during the exclusive breastfeeding, otherwise infants received the agents. the occurrence of eczema was diagnosed by the age of years by typical skin lesions found in children and chronic relapsing course. this last criterion was fulfilled if the child had had eczema for month or longer at the 24-month study visit and on at least one previous visit. eczema was considered atopic if it was associated with positive skin prick test at and/or -month study visit. none of the study subjects included in present study suffered from asthma or allergic rhinitis. also, all the infants were normal weight at the age of and months of age. the study protocol was approved by the ethics committee of the hospital district of southwest finland and subjects were enrolled in the study after written informed consent was obtained. faecal samples and dna extraction the faecal samples were taken from children at age of and months. the samples were aliquoted and frozen immediately after collection, and stored in −80°c. dna was extracted from faecal samples using the repeated bead-beating method as described previously. 16s rrna gene microarray analysis the composition of total microbiota was assessed by using the phylogenetic human intestinal tract chip as described previously, except for the amplification step, where cycles of end-point pcr were used. microarray analysis of all samples were performed in at least two independent hybridizations until satisfactory reproducibility was achieved. this study reports results of more than independent microarray hybridizations. the hitchip is a custom-made agilent microarray designed to comprehensively cover the diversity of the human intestinal microbiota. the array contains unique oligonucleotide probes targeting the v and v hypervariable regions of the 16s rrna gene and covering over intestinal bacterial phylotypes. the hitchip allows the analysis at three phylogenetic levels: phylum-like level, genus-like level and phylotype level. the details of the hitchip have previously been described, including its validation for phylogenetic fingerprinting and quantification. microarray data extraction and microbiota diversity assessment data were extracted from microarray images using the agilent feature extraction software, version.. normalization of microarray data was performed as described earlier. further data processing was performed by using a custom designed relational database running under the mysql database management system using r-based scripts. quantitative pcr quantitative pcr analysis of bifidobacterium genus and species was carried out in an applied biosystems fast real-time pcr system in a 96-well format and by using sybr green chemistry. the primers and their specificities are presented in additional file the pcr reactions and thermocycling conditions were as reported earlier. standards for qpcr were prepared as described in nermes et al.. samples were analysed in duplicate in at least two independent runs. statistical and data analyses statistical analysis of both qpcr and hitchip data was carried out with log-transformed data. in qpcr data, non-detected values were imputed with the half of the theoretical detection limit. for hitchip data, linear models with factors for treatment, health status, time point and breast-feeding with subsequent anova and contrast tests were used to determine the statistical differences between groups. in microarray data, cut-off values for positive responding probes were calculated as described before. in hitchip data the analysed values were summary values on phylum-like and genus-like level, obtained by summing the intensities from all the probes assigned to the respective phylum-like or genus-like phylogetic groups. totally phylum-like and genus-like level groups reached the detection threshold and were thus used in statistical analysis. the data is presented as mean with standard deviation values. redundancy analysis was performed by using the multivariate statistical analysis package canoco. rda plot shows bacterial groups principally contributing to the difference between the groups of subjects. the significance of separation in rda was assessed by monte carlo permutation procedure. the diversity of the microbial community assessed by hitchip was expressed as simpson’s reciprocal index of diversity as described before.  temporal development of microbiota the faecal microbiota of children at age of and months was analysed using the hitchip phylogenetic microarray. the diversity of total microbiota increased significantly with age, as the simpson’s the reciprocal diversity index has changed from ± to ± at age of and months, respectively. at the phylum-like level, significant changes in the relative abundances of major bacterial groups were detected. the most prominent decline in abundance was observed for actinobacteria that contributed % and % to the total signal at and months of age, respectively. signal intensities for actinobacteria were almost entirely obtained from bifidobacteria. this finding was consistent with quantitative pcr analysis, where total bifidobacteria counts decreased significantly with age. at the species level, the amounts of b. longum/infantis group, b. breve, b. bifidum, b. catenulatum group and b. adolescentis decreased over time as assessed by qpcr. in addition to actinobacteria, the relative abundance of bacilli decreased with age. all genus-like groups belonging to bacilli decreased, most of which not significantly as individual groups, but the sum effect at the phylum-like level was significant. an opposite trend was observed for the bacterial groups belonging to the major clostridium clusters, especially members of the clostridium cluster xiva which relative abundance increased with age from % to %.similar was observed for ruminococcus bromii et rel. group from clostridium cluster iv that increased from % to %. in total, genus-like phylogenetic groups changed significantly with age,, which further highlights the extensive compositional changes that the microbiota is undergoing during this period of life. genus-like groups with a p-value less than are presented in the table. analysis of the intestinal microbiota composition in relation to the health status when comparing the microbiota of the two groups of children at the age of months, pronounced differences were observed both in the microbial composition and the diversity. infants with eczema had a significantly more diverse total microbiota. analysis at the species-like level showed that a large number of bacterial species have different abundance between healthy and eczematous infants, although the individual p-values are not particularly small. the numerous, but mostly not significant, differences at the species-like level prompted us to look at the trends in microbiota differences at higher levels i.e. at the phylum-like and genus-like levels. analysis at the phylum-like level showed that the most remarkable difference between the groups was observed in the bacteroidetes, which were 3-fold more abundant in healthy children at the age of months. the effect of bacteroidetes is also clearly shown in the rda plot, which reveals the bacterial groups principally contributing to the difference among the groups of subjects. the microbiota differences between the health groups shown in the rda are significant as assessed by mcpp and a total of % of the variation within the dataset could be related to the health status of the infants. in contrast to the bacteroidetes, specific bacterial groups from the most abundant groups of the firmicutes phylum - clostridium clusters iv and xiva - were significantly more abundant in children with eczema. in summary, the multiple differences in specific bacterial groups result in microbiota profiles that are significantly distinct between healthy and eczematous infants as assessed by mcpp. * % of total hitchip signal nd, below the detection level. similar microbiota changes were observed between children with eczema and healthy children when only subjects in the placebo group were analysed, although the differences between the groups did not reach statistical difference due to the small number of subjects. no significant differences were observed for the bifidobacterial sub-community between the two groups of children using both hitchip and qpcr analyses. the comprehensive list of phylum-like and genus-like level data and p-values obtained by statistical analyses are presented in additional file and additional file respectively. notably, an indication towards altered microbiota composition in children with eczema was already identified at months, although the difference did not reach the level of statistical significance. a higher abundance of the clostridium cluster xiva bacteria was observed in infants with eczema than healthy controls. l. rhamnosus gg supplementation in early infancy has minor long-term effects on the microbiota composition when comparing the levels of hitchip signals between children from the placebo group and those who had received l. rhamnosus gg for their first months of life, no statistically significant differences were observed at the age of months. however, the supplementation with l. rhamnosus gg showed effects on three genus-like bacterial groups at the age of months i.e. a year after the cessation of the probiotic supplementation. the children that had received l. rhamnosus gg had higher levels of the butyrate-producing groups anaerostipes caccae et rel and eubacterium ventriosum et rel than those of placebo group. moreover, the placebo group children had higher levels of clostridium difficile et rel at months of age as compared to the lgg group children. the comprehensive list of phylum-like and genus-like level data and p-values obtained by statistical analyses are presented in additional file and additional file respectively. the effect of the probiotic supplementation on the microbiota composition within the group of healthy children or the group of children with eczema was not addressed due to the small number of subjects. discussion we used a high-throughput phylogenetic microarray to reveal alterations in the gut microbiota composition throughout early childhood. the used microarray has been developed and validated for determining the microbiota diversity and evaluating the relative proportions of genus-like or higher phylogenetic groups. moreover, it has been demonstrated that this microarray has a power equal to or higher than deep new generation sequencing and it has previously been used to compare the microbiota diversity and composition in a variety of studies. despite the highly significant increase of microbiota diversity with age, the diversity indeces at months of age are still relatively low when compared to the approximately two-fold higher indexes commonly observed in healthy adults. it has been suggested that by the age of to years the microbiota resembles that of an adult. our results show that microbiota succession continues at least until the age of months and most likely even further, because the bacterial diversity has still not reached the diversity of an adult person. thus, significant changes can be expected to occur in even after months of age. concerning the microbiota composition at months of age, our results are in agreement with earlier studies, except that we observed significant colonization by bifidobacteria in most of the children while in the study of palmer et al. bifidobacteria were not detected, possibly due to differences in dna extraction, pcr primers, demographic and geographic origin, dietary patterns of the infants or other confounding factors. primers used for pcr are often not so optimal for bifidobacteria than for other species and thus, high gc bacteria may perform less well in such pcrs. further, in our previous studies we have shown that mechanical lysis of faecal bacteria is essential and improves the detection of especially gram-positive bacteria including bifidobacteria. in the palmer et al. study, mechanical lysis by bead-beating was not applied, which may have hampered the detection of bifidobacteria. thus, we consider that the most likely explanation for the different results concerning bifidobacteria in our and palmer et al. study is the different dna extraction methods used. when comparing healthy and eczematous children we found statistically significant differences in microbiota composition only at months of age. the total microbiota of children with eczema was found to become significantly more diverse than the microbiota of children who remained healthy by months of age. interestingly, the total microbiota and particularly firmicutes diversity was higher in the eczema group children, although the difference with the healthy subjects was not statistically significant. abrahamsson et al. described the infants as having atopic eczema during the first two years of life, but the age at the onset of symptoms was not clarified. however, it can be concluded from the abrahamsson et al. and our study, both taking advantage of the high resolution microbiota assessment techniques, that after year of age the total microbiota diversity in children either developing or having eczema is comparable or even higher than that of healthy children. secretory iga has been suggested to play a role in shaping the microbiota composition and diversity. some early studies showed an association between the low levels of secretory iga and the risk of developing atopy and could suggest that the low iga levels permit establishment of a wider variety of bacteria and explain the higher bacterial diversity in children with eczema observed in this study. however, more recent studies have shown a higher concentration of secretory iga in children with allergic sensitization during the first years of life. another possible explanation for the increased bacterial diversity in children with eczema is the decreased levels or altered repertoire of antimicrobial peptides secreted into the gut lumen. these peptides, such as alpha- and beta-defensins, have at least two key roles at the mucosal interface: contributing to the host defense against enteric bacterial attachment and homeostatic control of the intestinal bacterial ecosystem. recently, decreased alpha-defensin levels and increased beta-defensin levels were associated with increased risk of developing atopy. to our knowledge, the levels of faecal antimicrobial peptides in children already having eczema have not been studied. however, a few studies have highlighted the role of alpha-defensins in microbiota composition and intestinal health. for example, genetic mutations resulting in decreased alpha-defensin expression have been associated with the susceptibility and severity of inflammatory bowel disease in humans and decreased alpha-defensins may have an effect on the differences observed in microbiota composition between healthy and diseased subjects. interestingly, mice deficient in production of active alpha -defensins were shown to have a decrease in bacteroidetes. the reason for decreased bacteroidetes levels in children with eczema in this study remains unaccountable, but alpha-defensins provide one possible explanation for our observation. also other host-dependent factors, such as the amount of mucus secretion and differences in mucus glycosylation may have an influence on the microbiota diversity and composition, as recently reviewed by maynard et al.. clearly, the role of intestinal iga levels, antimicrobial peptides and mucus secretion in shaping the gut microbiota in healthy and eczematous children warrants for further investigation. our results emphasize that the microbiota diversity in children with eczema should be further studied by using high-resolution techniques in order to define the favourable course of bacterial succession in early childhood and toddler age and to evaluate possible means to influence it. it was observed that children with eczema harbour more bacteria belonging to the clostridium cluster iv and clostridium cluster xiva. these bacteria are among the most abundant microbial groups detected in the healthy adult intestine. thus, prematurely occurring changes towards an adult-type microbiota were observed to take place in children with eczema. it has been suggested that resident bacteria may shape the hosts’ physiology, among others, by modulating the expression of genes involved in intestinal functions, such as postnatal intestinal maturation and the maintenance of mucosal barrier. it may be speculated that an infant-type microbiota supports adequate gut barrier function and tolerance against food allergens in an immature gut. infant-type microbiota may fortify the normal mucosal barrier function e.g. by affecting the maturation of the gut epithelium and immune functions in an optimal way and decrease the low-grade intestinal inflammation observable in subjects with eczema. maintenance of adequate mucosal barrier function may also play a role in the level of sensitisation to food-derived compounds. the complex host-microbe interactions in the intestinal epithelium are only recently beginning to be understood. furthermore, we observed decreased relative abundances of bacteria belonging to bacteroidetes in children with eczema. previous studies have reported an association between decreased amounts bacteroides spp. and the development of atopy and increased risk for atopic sensitization. bacteria belonging to the bacteroidetes are among the first groups colonizing the gut and they are typical intestinal habitants in healthy adults. bacteroides spp. are specialized in the breakdown of complex plant polysaccharides and their abundance has been associated with increased short-chain fatty acid concentrations in the infant gut after introduction of first solid foods. furthermore, b. fragilis polysaccharide has been shown in mice model to direct the cellular and physical maturation of the developing immune system via its ability to direct the development of cd4+ t cells, thus inducing the differentiation of th lineage and correction of the th1/th imbalance. together with our findings, these results suggest the significance of bacteroides spp. in the development and maintenance of healthy infant gut and balanced mucosal immunity and necessitate the role of these bacteria to be considered in future studies. when comparing healthy children with children with eczema we found statistically significant differences in microbiota composition only at months, but not at months of age. breast-feeding is known as a major factor influencing the microbiota composition in infancy. at months of age, the majority of children included in this study were still nursed and breast-feeding is likely to have had a strong influence on their microbiota composition at that age. thus, it seems that breast-feeding could have evened up the microbiota differences between the healthy and eczematous children and masked the eczema-associated changes, which came apparent and measurable at months of age after the withdrawal of breast-milk. while many studies addressed the impact of l. rhamnosus gg on health parameters, the short and long-term effect on the intestinal microbiota has only received limited attention. in the present intervention, the supplementation of l. rhamnosus gg continued until the age of months. interestingly, no significant effect on the microbiota composition was observed at the age of months, but instead the supplementation of l. rhamnosus gg in early life was observed to a induce long-term effect and small but significant changes between the intervention groups were observed one year later at the age of months. the observation that the c. difficile et rel. group bacteria were lower in the lgg groups as compared to placebo is of particular interest. previously, clostridium difficile colonization at the age of month has been associated with a higher risk of a diagnosis of atopic dermatitis at the age of years. the higher anaerostipes caccae et rel levels in the children that had received the l. rhamnosus gg supplementation is also a potentially beneficial effect, because a. caccae produces butyrate, which is an energy source for epithelial cells of colonic mucosa. bacteria belonging to the eubacterium ventriosum et rel group that were higher in the children that received the probiotic supplementation, also have shown to produce butyrate but have been less investigated. in mice, however, it has been shown that e. ventriosum was reduced in colitic mice as compared to non-colitic animals. to our knowledge this is the first high -throughput microbiota analysis study reporting the long-term effects of a probiotic strain on the microbiota composition in early life.  in conclusion, using a comprehensive microbial analysis approach we observed children with eczema to harbour a more diverse total microbiota and detected specific shifts in bacterial groups in different phylogenetic levels. the results indicate that aberrancies in microbiota composition are associated with eczema. our results also suggest that in children at high-risk for atopic disease, a diverse adult-type microbiota in too early childhood may be a potential risk factor and further strengthen the importance of early microbiota characterization and potential dietary modification. abbreviations et rel.: and relatives; hitchip: human intestinal tract chip; mcpp: monte carlo permutation procedure; qpcr: quantitative real-time polymerase chain reaction; rrna: ribosomal rna. competing interests the authors declare that they have no competing interests. authors’ contributions mk and es designed the original intervention study and organized the sample collection. infants were clinically examined by mk. ln, wmdv, rs and ss designed the current study. ln performed faecal microbial dna extraction, qpcr analyses and hitchip experiments. mr-s was involved in hitchip experiments. jn performed bioinformatic analyses. ln, rs and wmdv interpreted the results and wrote the paper. all authors read and approved the final manuscript. supplementary material additional file 1 basic characteristics of the study subjects. click here for file additional file 2 primers targeting bifidobacterium genus and species used in this study. click here for file additional file 3 differences in bifidobacterial composition of all children at and months of age as assessed by using quantitative pcr. click here for file additional file 4 the histograms showing the distribution of p-values obtained from the statistical analyses of species-like level of hitchip data at months. each bar represents how many species-like groups gave a p-value in the given range when the effect of different factors on microbiota composition were analysed. click here for file additional file 5 the microbiota differences of healthy and eczematous children from placebo group as assessed by hitchip analysis. click here for file additional file 6 bifidobacterial sub-communities in infants with eczema and healthy controls as assessed by quantitative pcr and hitchip analyses. click here for file additional file 7 phylum-like and genus-like hitchip data used in this study. data is presented as log-transformed values. a letter a refers to months samples and a letter d to months samples, respectively. click here for file additional file 8 p-values obtained from the statistical analysis of phylum-like and genus-like groups of hitchip data at months. p-values are not corrected and therefore indicate trend-like differences in the abundance of individual bacterial groups between the groups of infants. microbial groups that were over the detection level were included in the analysis. click here for file additional file 9 the microbiota differences between the intervention groups at the age of months as assessed by hitchip analysis. click here for file acknowledgements this work was funded by finnish funding agency for technology and innovation. in addition, the academy of finland is acknowledged for financial support. hans heilig, outi immonen and alla kaljukivi are thanked for their excellent technical assistance. we thank professor airi palva for valuable discussions and her support to carry out this study.
a diverse and adult-type microbiota in early childhood is associated with eczema and it may contribute to the perpetuation of eczema. probiotic lactobacillus rhamnosus gg supplementation in early infancy was observed to have minor long-term effects on the microbiota composition. deviations in composition and diversity of intestinal microbiota in infancy have been associated with both the development and recurrence of atopic eczema. thus, we decided to use a deep and global microarray-based method to characterize the diversity and temporal changes of the intestinal microbiota in infancy and to define specific bacterial signatures associated with eczema. faecal microbiota at and months of age were analysed from infants selected from a prospective follow-up study based on the availability of faecal samples. the infants were originally randomized to receive either lactobacillus rhamnosus gg or placebo. children with eczema harboured a more diverse total microbiota than control subjects as assessed by the simpson’s reciprocal diversity index of the microarray profiles. microbiota of children suffering from eczema had increased abundance of the clostridium clusters iv and xiva, which are typically abundant in adults.
 deviations in composition and diversity of intestinal microbiota in infancy have been associated with both the development and recurrence of atopic eczema. thus, we decided to use a deep and global microarray-based method to characterize the diversity and temporal changes of the intestinal microbiota in infancy and to define specific bacterial signatures associated with eczema. faecal microbiota at and months of age were analysed from infants selected from a prospective follow-up study based on the availability of faecal samples. the infants were originally randomized to receive either lactobacillus rhamnosus gg or placebo.  children with eczema harboured a more diverse total microbiota than control subjects as assessed by the simpson’s reciprocal diversity index of the microarray profiles. composition of the microbiota did not differ between study groups at age of months, but was significantly different at age of months as assessed by mcpp. at this age healthy children harboured -fold greater amount of members of the bacteroidetes. microbiota of children suffering from eczema had increased abundance of the clostridium clusters iv and xiva, which are typically abundant in adults. probiotic lactobacillus rhamnosus gg supplementation in early infancy was observed to have minor long-term effects on the microbiota composition.  a diverse and adult-type microbiota in early childhood is associated with eczema and it may contribute to the perpetuation of eczema. infantintestinal microbiotamicrobiota diversityphylogenetic microarrayeczema 