You are here

  • Home
  • Archive
  • Volume 63, Issue 4
  • Decreased gut microbiota diversity, delayed Bacteroidetes colonisation and reduced Th1 responses in infants delivered by Caesarean section

Article Text

Article menu
Gut microbiota
Original article
Decreased gut microbiota diversity, delayed Bacteroidetes colonisation and reduced Th1 responses in infants delivered by Caesarean section
  1. Hedvig E Jakobsson1,2,
  2. Thomas R Abrahamsson3,
  3. Maria C Jenmalm3,4,
  4. Keith Harris5,
  5. Christopher Quince5,
  6. Cecilia Jernberg1,
  7. Bengt Björkstén6,7,
  8. Lars Engstrand2,
  9. Anders F Andersson8
  1. 1Department of Preparedness, Swedish Institute for Communicable Disease Control, Solna, Sweden
  2. 2Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
  3. 3Department of Clinical and Experimental Medicine, Division of Pediatrics, Linköping University, Linköping, Sweden
  4. 4Department of Clinical and Experimental Medicine, Division of Inflammation Medicine, Linköping University, Linköping, Sweden
  5. 5School of Engineering, University of Glasgow, Glasgow, UK
  6. 6Institute of Environmental Medicine, Karolinska Institutet, Stockholm, Sweden
  7. 7School of Health and Medical Sciences, Örebro University, Örebro, Sweden
  8. 8KTH Royal Institute of Technology, Science for Life Laboratory, School of Biotechnology, Division of Gene Technology, Solna, Sweden
  1. Correspondence to Dr Anders F Andersson: KTH Royal Institute of Technology, Science for Life Laboratory, School of Biotechnology, Division of Gene Technology, PO Box 1031, Solna SE-171 21, Sweden; anders.andersson{at}scilifelab.se and Professor Lars Engstrand, Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm SE-171 82, Sweden; lars.engstrand{at}ki.se

Abstract

Objective The early intestinal microbiota exerts important stimuli for immune development, and a reduced microbial exposure as well as caesarean section (CS) has been associated with the development of allergic disease. Here we address how microbiota development in infants is affected by mode of delivery, and relate differences in colonisation patterns to the maturation of a balanced Th1/Th2 immune response.

Design The postnatal intestinal colonisation pattern was investigated in 24 infants, born vaginally (15) or by CS (nine). The intestinal microbiota were characterised using pyrosequencing of 16S rRNA genes at 1 week and 1, 3, 6, 12 and 24 months after birth. Venous blood levels of Th1- and Th2-associated chemokines were measured at 6, 12 and 24 months.

Results Infants born through CS had lower total microbiota diversity during the first 2 years of life. CS delivered infants also had a lower abundance and diversity of the Bacteroidetes phylum and were less often colonised with the Bacteroidetes phylum. Infants born through CS had significantly lower levels of the Th1-associated chemokines CXCL10 and CXCL11 in blood.

Conclusions CS was associated with a lower total microbial diversity, delayed colonisation of the Bacteroidetes phylum and reduced Th1 responses during the first 2 years of life.

  • Chemokines
  • Infant Gut
  • Intestinal Bacteria
  • Intestinal Microbiology
  • Molecular Biology

https://doi.org/10.1136/gutjnl-2012-303249

Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Significance of this study

    What is already known on this subject?

    • The infant gut microbiota diversity increases during the first years of life.

    • The microbiota composition differs between infants born by caesarean section (CS) or vaginal delivery with a delayed colonisation of the genus Bacteroides.

    • Bacterial colonisation is necessary for the development of the immune system and immune regulation.

    • An association between CS delivery and the development of allergic disease has been observed in several studies.

    What are the new findings?

    • The total microbiota diversity is lower in CS than vaginally delivered infants through the first 2 years of life.

    • The diversity of the Bacteroidetes phylum is lower in CS born infants during the first 2 years of life.

    • Vaginal delivery is associated with increased circulating levels of Th1-associated chemokines during infancy.

    How might it impact on clinical practice in the foreseeable future?

    • Deeper knowledge of the impact of delivery mode on microbiota composition and immune regulation may lead to improved allergy preventive strategies.

    Introduction

    The gastrointestinal tract of the newborn infant is considered to be sterile. Bacteria from the environment, mainly from the mother, colonise the infant gut immediately following birth. Dominant members of anaerobic Firmicutes and Bacteroidetes do not appear to grow outside the gut and hence need to be transmitted between human hosts.1 To what extent the transmission occurs from mother to offspring is not clear, but differences in microbiota composition depending on delivery mode indicate a mother–child transmission during vaginal delivery. A recent study based on pyrosequencing of 16S rRNA genes demonstrated that the microbiota of vaginally delivered (VD) neonates (<24 h postdelivery) resembled the vaginal microbiota of their own mother and was similar across multiple body habitats (skin, oral, nasopharynx and faeces), while in neonates born by caesarean section (CS), it resembled the mother's skin microbiota.2 While this study provided evidence that microbiota from the birth canal is transferred from mother to child, providing an inoculum for the initial microbiota, it remains to be shown that specific gut microbes are successfully transmitted during vaginal delivery.

    The incidence of caesarean delivery has increased from 5% in the 1970 s to more than 60% in some hospitals in China according to recent reports.3 The early colonisation pattern differs between VD infants and those delivered by CS, including a delayed colonisation of, for example, Bacteroides and Bifidobacterium spp. in CS infants.4 ,5 The influence of delivery mode on gut microbiota development has not previously been longitudinally characterised using powerful cultivation-independent microbiological methods, however. The intestinal microbiota is important for the development of the immune system and immunological tolerance.6 Differences in the postnatal microbial colonisation may explain the higher incidence of immune mediated diseases such as allergy in children born by CS as compared with those born vaginally.7 ,8 Indeed, allergic disease has been associated with a low prevalence of Bacteroides and Bifidobacterium,9 ,10 and a low intestinal microbiota biodiversity in early infancy appears to have an impact on the development of allergic disease later in life.11–13 A failure of Th2-silencing during maturation of the immune system may underlie development of Th2-mediated allergic disease.14 Appropriate microbial stimulation may be required to avoid this pathophysiological process, as early differences in the gut microbiota may shape later immune development.6 ,15 The influence of CS on immune development is largely unknown, however.16 The aim of the present study was to monitor the development of the infant intestinal microbiota in babies born vaginally and through CS, and to relate the findings to the maternal microbiota and to Th1- and Th2-associated chemokine levels during infancy.

    Materials and methods

    Ethics

    The human ethic committee at Linköping University, Linköping, Sweden, approved the study. Informed consent was obtained from both parents before inclusion.

    Subjects and sample collection

    The study group comprised 24 healthy women and their infants. Nine of the infants were born by CS and 15 by VD (see online supplementary table S1 regarding, sex, delivery, birth weight, the use of antimicrobials and length of breast feeding for the different infants). Seven out of nine mothers, who gave birth through CS, were given antibiotics prophylactically during the surgery (see online supplementary table S2). This was done after the delivery, however, and thus the infants were not exposed to antibiotics via the placenta. No children were treated with antibiotics during the neonatal period. Twenty of the infants (83%) were partly breast fed until at least 6 months of age. The women and children included in this study were part of a larger study assessing the prevention of allergic disease by probiotics17 (ClinicalTrials.gov ID NCT01285830) and they all received placebo. Stool samples were collected from the mothers 1 week after delivery and from the children at 1 week, 1, 3, 6, 12 and 24 months. The faecal samples were immediately frozen at −20°C following collections and later stored at −70°C. Samples were collected in 2002–2005 and stored in −70°C until DNA extractions were conducted in 2009. No systematic differences in storage times existed between the VD and CS samples (mean months: VD 80±7, CS 78±4, Student t test, p=0.42).

    DNA extraction and 16S rRNA gene amplification

    The DNA extraction and 16S rRNA gene amplification were performed as described previously18 with the following modifications; the primer pair used, targeting the variable regions 3 and 4 of the 16S rRNA gene, was 341f 5′CCTACGGGNGGCWGCAG with adaptor B and 805r 5′GACTACHVGGGTATCTAATCC with adaptor A and a sample-specific barcode sequence consisting of five nucleotides.19 The barcodes contained no homopolymers and a pair of barcodes differed in at least two positions. The 341f-805r primer pair was shown to be the least biased among 512 primer pairs evaluated in silico for bacterial amplification and was experimentally shown to give a taxonomic composition similar to shotgun metagenomics.20 The primer pair has good coverage of bacterial groups typically found in the human lower intestine. For phylum Bacteroidetes, 121 862 of 132 120, for phylum Firmicutes 376 912 of 406 649, for phylum Proteobacteria 336 471 of 368 375 and for family Bifidobacteriaceae 1112 of the 1239 sequences are matched, when considering Ribosomal Database Project (RDP) sequences that span the region. A negative PCR reaction without template was also included for all primer pairs in each run. The PCR products with approximate lengths of 450 bp were purified with AMPure beads (Becton Dickinson, Franklin Lakes, New Jersey, USA) using a Magnet Particle Separator (Invitrogen, Carlsbad, California, USA). The concentrations of the purified products were measured by Qubit fluorometer (Invitrogen), the quality was assessed on a Bioanalyzer 2100 (Agilent, Santa Clara, California, USA), and the samples were amplified in PCR-mixture-in-oil emulsions and sequenced from the 805r primer on different lanes of a 2-lane PicoTiterPlate on a Genome Sequencer FLX system (Roche, Basel, Switzerland) at the Swedish Institute for Communicable Disease Control (Solna, Sweden). Sequence processing was carried out with the AmpliconNoise software package21 correcting for errors introduced in the PCR and pyrosequencing as well as removing chimeric sequences. Also, reads lacking a correct primer and/or having less than 360 successful pyrosequencing flows were excluded.21 Denoised sequences were trimmed to 198 bp after primer and barcode removal and clustered by complete linkage clustering into operational taxonomic units (OTUs) at the 97% similarity level using AmpliconNoise.

    Taxonomic classification

    Each denoised sequence, as well as the most abundant sequence for each OTU, was BLAST searched with default parameters against a local BLAST database comprising 836 814 near full-length bacterial 16S rRNA gene sequences from the RDP v.10.10.22 The sequences inherited the taxonomic annotation (down to genus level) of the best scoring RDP hit fulfilling the criteria of ≥95% identity over an alignment of length ≥180 bp. If no such hit was found the sequence was classified as ‘no match’. If multiple best hits (of same score) were found, the taxonomy was set to the most detailed level of taxonomy shared by the best hits. The majority of reads had an RDP relative within 95% sequence similarity and were hence of bacterial origin. After removal of pyrosequencing noise and chimeric sequences using the AmpliconNoise package,21 357 685 high quality, typically 198 bp long, sequence reads remained, with 828–4395 reads per sample (mean=2129). These corresponded to 3048 unique sequences and 1818 OTUs, clustered at 97% similarity level using complete linkage clustering.

    Sample clustering

    The online version of Fast Unifrac (http://bmf2.colorado.edu/fastunifrac/)23 was used to calculate weighted sample distances by mapping our OTU sequences with BLAST onto the Greengenes reference sequences (downloaded from the Fast Unifrac web page, May 2009) and using the corresponding Greengenes tree. A principal coordinates analysis plot based on all pair-wise sample distance was created on the Fast Unifrac web page. Our OTU sequences were mapped onto 154 Greengenes sequences.

    Statistical testing for over-representation or under-representation of bacterial lineages

    Statistical tests of over-representation or under-representation of bacterial lineages among sample groups were made at the phylum and genus levels using Wilcoxon rank-sum test. To correct for multiple testing, the p values were converted to false discovery rate values (Q values).

    Diversity estimations

    Shannon diversity index was calculated as H=–Σ log(pi)pi, where pi denotes the relative frequency of OTU i,24 Pielou's evenness index as –Σ log(pi)pi/log(Sobs), where Sobs denotes the number of observed OTUs in the sample,25 and Chao1 richness estimate as Sobs + n1(n1 − 1)/n2(n2 − 1), where n1 and n2 are the number of observed singleton and doubleton OTUs, respectively.26 Since these parameters are influenced by sequencing depth and the sequencing depth differed between samples, we subsampled (with replacement) 1400 reads from each sample, counted the occurrences of the corresponding OTUs and performed the diversity calculations on these counts. This was repeated 10 times and averages of the diversity parameters calculated and used for further analysis. Four (out of 168) samples had fewer than 1400 reads and were excluded from this part of the analysis (three VD infants: one at 1 week, one at 3 months and one at 12 months, and one CS infant at 1 month). Diversity calculations and statistics were done with the R software (http://www.r-project.org/) and the R package vegan (http://cran.r-project.org/web/packages/vegan/). Repeated measures ANOVA were employed, using the Shannon diversity index at the different time points (1 week, 1, 3, 6, 12 and 24 months) as the repeated measures.

    Statistical testing of mother–child overlap in sequence types

    For each time point and for each infant, we calculated the number of specific sequences (sequence types) shared with its mother/number of sequence types observed in the infant (Rown). Likewise, we calculated the average number of sequence types shared with other mothers/number of sequence types observed in the infant (Rother). We then compared the Rown and Rother values pair-wise for all infants within each group (VD or CS) with the Wilcoxon signed rank test.

    Chemokine analyses in venous blood and association with mode of delivery

    Venous blood was collected at 6 (n=24), 12 (n=24) and 24 months (n=24) and stored as plasma or serum in −20°C pending analysis. The Th1-associated chemokines CXCL10 and CXCL11 and the Th2-associated chemokines CCL17 and CCL22 were analysed with an inhouse multiplexed Luminex assay.27 ,28 The limit of detection was 6 pg/mL for CXCL10 and CXCL11 and 2 pg/mL for CCL17 and for CCL22. All samples were analysed in duplicates and the sample was reanalysed if the coefficient of variance was >15%. The chemokine levels were compared between infants being born vaginally or by CS by repeated measures ANOVA using log transformed chemokine levels at the different time points (6, 12 and 24 months) as the repeated measures.

    Results

    Microbiota succession in infants

    At the phylum level, the microbiota developed in a similar fashion in infants delivered vaginally and by CS, with a gradual decline in Proteobacteria from 1 week to 24 months, a peak of Actinobacteria at 3 months, an expansion of Firmicutes from 3 months and onward, and the emergence of Verrucomicrobia at around 6 months of age (figure 1, see online supplementary table S3). However, a notable difference between the VD and CS infants was the higher proportion of Bacteroidetes in VD infants during the first 12 months (significantly higher at 1 week, 3 months and 12 months; see online supplementary table S3).

    Figure 1
    Figure 1

    Phylum level microbiota composition in mothers and their infants at 1 week, 1, 3, 6, 12 and 24 months. The mean relative abundances (%) of the most abundant bacterial phyla in the 15 vaginally delivered infants (A) and nine caesarean section infants (B) as well as in their mothers are shown.

    The maternal microbiota resembled the typical adult flora as demonstrated in several previous studies18 ,29 ,30 and was independent of delivery mode. The Firmicutes was the dominating phylum, representing 74% and 71% in mean relative abundance for the 15 VD and nine CS mothers, respectively, followed by Bacteroidetes (16% and 13%, respectively), Actinobacteria (7% and 12%, respectively), Proteobacteria (3% and 2%, respectively) and Verrucomicrobia (1% and 2%, respectively) (figure 1, see online supplementary table S3).

    The relative abundance of the major genera found in the infants and mothers are illustrated in figure 2A,B (see also online supplementary table S4). Infants in the VD group were colonised by Bacteroides to a greater extent than in the CS group (significantly more at 1 week (11/15 in VD vs 1/9 in CS; Fischer's exact test p=0.005), 3 months (11/15 vs 1/9; p=0.005) and 12 months (14/15 vs 4/9; p=0.015)). At 1 month of age, Bifidobacterium dominated the microbiota in both groups. The Enterococcus genus was found in significantly higher relative abundance in the CS compared with the VD infants at 1 month (p<0.0001; see online supplementary table S4). Following 6 months of age, there was a gradual increase in previous low abundant genera in both VD and CS infants. At 24 months of age, Bacteroides and several genera belonging to the Clostridia class, for example Ruminococcus, a dominant member of the adult microbiota, also dominated the infant microbiota.

    Figure 2
    Figure 2

    Genus level microbiota composition in mothers and their infants at 1 week, 1, 3, 6, 12 and 24 months. The mean relative abundances (%) of the most abundant bacterial genera in the 15 vaginally delivered infants (A) and nine caesarean section infants (B) as well as in their mothers are shown. Only genera comprising ≥1% of the total community were included. Ac, Actinobacteria; Ba, Bacteroides; Fi, Firmicutes; Pr, Proteobacteria.

    The gradual shift in community composition was accompanied by an increase in α-diversity over time, with a significant increase in Shannon diversity index between each pair of succeeding time points from 3 months of age onwards (figure 3). Similar results were obtained for evenness (Pielou's index; see online supplementary figure S1A) and estimated richness (Chao1; see online supplementary figure S1B). The low increase in diversity during the first 3 months may be related to that most (83%) infants were exclusively breast fed up to this age (see online supplementary table S1). CS delivery was associated with significantly lower total microbiota diversity when considering all time points in the infants (p=0.047 with repeated measures ANOVA; see online supplementary table S5). At individual time points, the total microbiota diversity was significantly lower in the CS delivered infants at 12 months (see online supplementary table S5; figure 3). The total microbiota diversity did not differ significantly between the VD and CS mothers (see online supplementary table S6). Narrowing the analysis to the phylum level, CS delivery was associated with a lower diversity of the Bacteroidetes phylum when considering all time points (p=0.002; see online supplementary table S5). For individual time points, the Bacteroidetes diversity was significantly lower in the CS infants than VD infants at 1, 3, 12 and 24 months (figure 4; see online supplementary table S6). The other phyla, Firmicutes, Proteobacteria and Actinobacteria, did not display any consistent differences in Shannon diversity between the groups, although diversity of Firmicutes was significantly lower at 12 months and Proteobacteria at 24 months in the CS infants (see online supplementary table S6).

    Figure 3
    Figure 3

    Increase in faecal microbiota α-diversity over time. Distributions of Shannon diversity indices displayed for the 15 vaginally delivered (VD) infants and nine caesarean section (CS) infants at 1 week, 1, 3, 6, 12 and 24 months, and for their mothers. In all, 50% of the data points reside within boxes, 75% within whiskers and median values are indicated by horizontal lines within boxes (circles indicate individual values). Wilcoxon signed rank tests were conducted to compare Shannon diversity between adjacent time points, and Wilcoxon rank-sum tests to compare diversity between delivery modes within time points; ***p<0.001, **p<0.01 and *p<0.05.

    Figure 4
    Figure 4

    Increase in Bacteroidetes α-diversity over time. Distributions of Shannon diversity indices displayed for the 15 vaginally delivered (VD) infants and nine caesarean section (CS) infants at 1 week, 1, 3, 6, 12 and 24 months, and for their mothers. **p<0.01 and *p<0.05. In all, 50% of the data points reside within boxes, 75% within whiskers and median values are indicated by horizontal lines within boxes (circles indicate individual values).

    A principal coordinates analysis plot based on pair-wise sample community differences calculated with UniFrac23 is illustrated in figure 5. As shown in a previous study31 the microbial communities became more uniform across infants over time, which is also evident when comparing the distributions of pair-wise community differences at each time point (see online supplementary figure S2). By 24 months of age, the communities closely resembled those of the mothers (figure 5). Although community composition converged to an adult-type microbiota, diversity estimates were still significantly lower at 24 months than in the mothers (figure 3; see online supplementary figure S1). Birth weight, antibiotic intake during the time course, sex and breast feeding had no apparent impact on the microbiota composition at any time point (see online supplementary figures S3–S6).

    Figure 5
    Figure 5

    Individuality and convergence of infant microbiota. Principal coordinates analysis was performed on all pair-wise community differences (calculated with UniFrac23), and samples from infants at 1 week, 1, 3, 6, 12 and 24 months, and from the mothers are highlighted in the different boxes. Vaginally delivered and caesarean section infants/mothers are displayed in red and blue, respectively.

    In order to investigate possible mother–child transmission of bacteria, the presence of specific unique sequences (sequence types) was compared between infants and their mothers, as well as between infants and other mothers. For each mother–child pair, we calculated the fraction of the sequence types found in the child that were also found in the mother (number of shared sequence types/number of child sequence types). The VD infants shared a significantly higher proportion of sequence types with their own mother than with the other mothers at 1 week and 24 months when considering all bacterial taxa (see online supplementary table S7). Considering one phylum at a time, Bacteroidetes displayed a significantly higher overlap at 1 and 6 months, and Firmicutes at 6 and 24 months. No significant overlap was observed for the CS delivered infants at any time point.

    CS was associated with moderately lower levels of the Th1-associated chemokines CXCL10 and CXCL11 (table 1). The levels of the Th2-chemokines CCL17 and CCL22 did not differ significantly between the birth modes.

    View this table:
    Table 1

    Repeated measures ANOVAs to test whether there were any significant differences in CXCL10 and CXCL11 levels during the first 2 years of life between caesarean section (CS) and vaginally delivered (VD) infants

    Discussion

    Microbial colonisation of the infant gut gastrointestinal tract is important for the postnatal development of the immune system. In this study, CS delivered infants who are not entering the birth canal of the mother either lacked or displayed a delayed colonisation of one of the major gut phylum, the Bacteroidetes. The colonisation of this phylum was delayed by up to 1 year for some infants. The total microbiota diversity was also lower in the CS infants, probably largely as a consequence of the lack of this phylum. This was not a consequence of antimicrobial treatment, as none of the CS mothers were given antibiotics before surgery and the microbial diversity did not differ between mothers who were given antibiotics prophylactically and those who did not receive antibiotics. Comparisons of intestinal microbiota have not conclusively confirmed bacterial transmission,31 ,32 although the genus Bacteroides has been proposed to be transmitted from the maternal gut.33 ,34 Our study corroborates earlier studies reporting a delayed colonisation of Bacteroides in babies delivered by CS.4 ,5 In addition, our study provides evidence that specific lineages of the intestinal microbiota, as defined by 16S rRNA gene sequences, are transmitted from mother to child during vaginal delivery.

    It is important to note that bacterial composition changes as a consequence of freezing the faecal samples35 ,36 and that PCR amplification can induce taxon-specific biases. However, there were no significant differences in storage times between the VD and CS samples and since the same primer pair and PCR conditions were used for all samples, these effects should not contribute to the observed differences in microbiota composition and α-diversity between the sample groups.

    The genus Enterococcus, which is a typical faecal bacterium, is usually acquired during the first week of life.37 Colonisation has previously not been shown to depend on delivery mode, suggesting other sources in addition to the maternal intestinal microbiota,5 such as the environment38 and breast milk.39 We found that CS infants had a higher relative abundance of Enterococcus at 1 month of age, suggesting that the lack of bacteria transmitted through vaginal delivery favours the growth and colonisation of enterococci.

    Appropriate microbial stimulation during infancy is required for the development of a more balanced immune phenotype, including maturation of Th1-like responses and appropriate development of regulatory T cell responses.6 ,40 ,41 It is well known that early life events occurring during critical windows of immune development can have long-term impact on immune-mediated diseases such as allergy,15 ,42 ,43 diabetes and inflammatory bowel disease. We hypothesised that early differences in the gut microbiota could shape later immune responsiveness, influencing the Th1 maturation trajectory. Our findings of a lower microbial diversity in the CS infants and lower circulating levels of the Th1-related chemokines CXCL10 and CXCL11 support this view. Previous studies have shown that Bacteroides fragilis exert strong effects on the immune system. This is mediated by the capsular polysaccharide, which enhances T cell mediated immune responses and affects the Th1/Th2 balance.44 ,45 Furthermore, Bacteroides thetaiotamicron is also known to affect the immune system.46 Thus, the lower abundance of Bacteroides among the CS infants may be a contributing factor to the observed differences in Th1-associated chemokines. Future studies with larger sample sizes will be able to address the effects of individual microbes on chemokine levels.

    With few exceptions,31 previous studies have reported Bifidobacterium to be one of the dominant genera of the early infant intestinal microbiota,4 ,33 ,47––49 and changes in relative abundance of this genus have been related to delivery mode.4 ,16 ,34 Also in the present study, Bifidobacterium was the dominant genus from 1 to 12 months of age with a gradual decline following weaning. The abundance was not affected by delivery mode, however, and we could not detect any significant overlap in the mothers’ and babies’ rRNA sequences. Hence, Bifidobacterium could primarily be transmitted from the breast milk, and to a lesser extent from the intestinal micobiota, as suggested but not confirmed previously.50

    In accordance with recent studies,31 ,47 ,48 our results demonstrate considerable individual differences in the microbial succession during the first year of life. This is probably a result of differences in time of weaning, and incidental exposure of bacteria from the environment. Community composition converges to an adult-like state within 2 years. However, even at 2 years the microbiota appears not to be fully developed, since the diversity was significantly lower than in mothers. This was also evident in a recent study reporting a lower microbial diversity in a 2.5-year-old child than in its mother.51

    An association between CS delivery and the development of allergic disease has been observed in several studies52 ,53 and a lower microbial diversity has been observed in allergic infants before onset of disease.13 We conclude that CS is associated with a lower bacterial diversity during the first 2 years of life, a lower abundance and diversity of the phylum Bacteroidetes, and lower circulating levels of Th1-associated chemokines during infancy.

    Acknowledgments

    We thank Mrs Lena Lindell, Linköping, Mrs Elisabeth Andersson, Norrköping, Mrs Linnea Andersson, Jönköping, and Mrs Eivor Folkesson, Motala, Dr Göran Oldaeus, Jönköping, and Dr Ted Jacobsson, Linköping, for their brilliant and enthusiastic work guiding the families through the study and all the sampling procedures. We also thank Ms Martina Abelius, Mrs Anne-Marie Fornander and Ms Anna Forsberg in Linköping for excellent technical assistance.

    References

    1. ↵
      1. Ley RE,
      2. Peterson DA,
      3. Gordon JI
      . Ecological and evolutionary forces shaping microbial diversity in the human intestine. Cell 2006;124:837–48.
      OpenUrlCrossRefPubMedWeb of Science
    2. ↵
      1. Dominguez-Bello MG,
      2. Costello EK,
      3. Contreras M,
      4. et al
      . Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci USA 2010;107:11971–5.
      OpenUrlAbstract/FREE Full Text
    3. ↵
      1. Sufang G,
      2. Padmadas SS,
      3. Fengmin Z,
      4. et al
      . Delivery settings and caesarean section rates in China. Bull World Health Organ 2007;85:755–62.
      OpenUrlCrossRefPubMedWeb of Science
    4. ↵
      1. Penders J,
      2. Thijs C,
      3. Vink C,
      4. et al
      . Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics 2006;118:511–21.
      OpenUrlAbstract/FREE Full Text
    5. ↵
      1. Adlerberth I,
      2. Lindberg E,
      3. Åberg N,
      4. et al
      . Reduced enterobacterial and increased staphylococcal colonization of the infantile bowel: an effect of hygienic lifestyle? Pediatr Res 2006;59:96–101.
      OpenUrlCrossRefPubMedWeb of Science
    6. ↵
      1. McLoughlin RM,
      2. Mills KH
      . Influence of gastrointestinal commensal bacteria on the immune responses that mediate allergy and asthma. J Allergy Clin Immunol 2011;127:1097–107.
      OpenUrlCrossRefPubMedWeb of Science
    7. ↵
      1. Eggesbø M,
      2. Botten G,
      3. Stigum H,
      4. et al
      . Is delivery by cesarean section a risk factor for food allergy? J Allergy Clin Immunol 2003;112:420–6.
      OpenUrlCrossRefPubMedWeb of Science
    8. ↵
      1. Bager P,
      2. Wohlfahrt J,
      3. Westergaard T
      . Caesarean delivery and risk of atopy and allergic disease: meta-analyses. Clin Exp Allergy 2008;38:634–42.
      OpenUrlCrossRefPubMedWeb of Science
    9. ↵
      1. Watanabe S,
      2. Narisawa Y,
      3. Arase S,
      4. et al
      . Differences in fecal microflora between patients with atopic dermatitis and healthy control subjects. J Allergy Clin Immunol 2003;111:587–91.
      OpenUrlCrossRefPubMedWeb of Science
    10. ↵
      1. Björkstén B,
      2. Naaber P,
      3. Sepp E,
      4. et al
      . The intestinal microflora in allergic Estonian and Swedish 2-year-old children. Clin Exp Allergy 1999;29:1–5.
      OpenUrlPubMedWeb of Science
    11. ↵
      1. Wang M,
      2. Karlsson C,
      3. Olsson C,
      4. et al
      . Reduced diversity in the early fecal microbiota of infants with atopic eczema. J Allergy Clin Immunol 2008;121:129–34.
      OpenUrlCrossRefPubMedWeb of Science
      1. Forno E,
      2. Onderdonk AB,
      3. McCracken J,
      4. et al
      . Diversity of the gut microbiota and eczema in early life. Clin Mol Allergy 2008;6:11.
      OpenUrlCrossRefPubMed
    13. ↵
      1. Abrahamsson TR,
      2. Jakobsson HE,
      3. Andersson AF,
      4. et al
      . Low diversity of the gut microbiota in infants with atopic eczema. J Allergy Clin Immunol 2012;129:434–40, 40 e1–2.
      OpenUrlCrossRefPubMedWeb of Science
    14. ↵
      1. Böttcher MF,
      2. Jenmalm MC,
      3. Björkstén B
      . Immune responses to birch in young children during their first 7 years of life. Clin Exp Allergy 2002;32:1690–8.
      OpenUrlCrossRefPubMedWeb of Science
    15. ↵
      1. Jenmalm MC
      . Childhood immune maturation and allergy development: regulation by maternal immunity and microbial exposure. Am J Reprod Immunol 2011;66(Suppl 1):75–80.
      OpenUrlCrossRefPubMed
    16. ↵
      1. Huurre A,
      2. Kalliomäki M,
      3. Rautava S,
      4. et al
      . Mode of delivery—effects on gut microbiota and humoral immunity. Neonatology 2008;93:236–40.
      OpenUrlCrossRefPubMedWeb of Science
    17. ↵
      1. Abrahamsson TR,
      2. Jakobsson T,
      3. Böttcher MF,
      4. et al
      . Probiotics in prevention of IgE-associated eczema: a double-blind, randomized, placebo-controlled trial. J Allergy Clin Immunol 2007;119:1174–80.
      OpenUrlCrossRefPubMedWeb of Science
    18. ↵
      1. Andersson AF,
      2. Lindberg M,
      3. Jakobsson H,
      4. et al
      . Comparative analysis of human gut microbiota by barcoded pyrosequencing. PLoS ONE 2008;3:e2836.
      OpenUrlCrossRefPubMed
    19. ↵
      1. Herlemann DP,
      2. Labrenz M,
      3. Jurgens K,
      4. et al
      . Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. Isme J 2011;5:1571–9.
      OpenUrlCrossRefPubMedWeb of Science
    20. ↵
      1. Klindworth A,
      2. Pruesse E,
      3. Schweer T,
      4. et al
      . Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res 2013;41:e1.
      OpenUrlAbstract/FREE Full Text
    21. ↵
      1. Quince C,
      2. Lanzen A,
      3. Davenport RJ,
      4. et al
      . Removing noise from pyrosequenced amplicons. BMC Bioinformatics 2011;12:38.
      OpenUrlCrossRefPubMed
    22. ↵
      1. Maidak BL,
      2. Cole JR,
      3. Lilburn TG,
      4. et al
      . The RDP (Ribosomal Database Project) continues. Nucleic Acids Res 2000;28:173–4.
      OpenUrlAbstract/FREE Full Text
    23. ↵
      Hamady M, Lozupone C, Knight R. Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J 2010;4:17–27.
    24. ↵
      1. Hayek LCaMAB
      . Surveying natural populations. New York: Columbia University Press, 1996.
    25. ↵
      1. Pielou EC
      . The measurement of diversity in different types of biological collections. J Theor Biol 1966;13:131–44.
      OpenUrlCrossRefWeb of Science
    26. ↵
      1. Chao A
      . Non-parametric estimation of the number of classes in a population. Scand J Stat 1984;11:265–70.
      OpenUrlWeb of Science
    27. ↵
      1. de Jager W,
      2. Prakken BJ,
      3. Bijlsma JW,
      4. et al
      . Improved multiplex immunoassay performance in human plasma and synovial fluid following removal of interfering heterophilic antibodies. J Immunol Methods 2005;300:124–35.
      OpenUrlCrossRefPubMedWeb of Science
    28. ↵
      1. Abrahamsson TR,
      2. Abelius M,
      3. Forsberg A,
      4. et al
      . A Th1/Th2-associated chemokine imbalance during infancy in children developing eczema, wheeze and sensitization. Clin Exp Allergy 2011;41:1729–39.
      OpenUrlCrossRefPubMed
    29. ↵
      1. Dethlefsen L,
      2. Huse S,
      3. Sogin ML,
      4. et al
      . The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing. PLoS Biol 2008;6:e280.
      OpenUrlCrossRefPubMed
    30. ↵
      1. Eckburg PB,
      2. Bik EM,
      3. Bernstein CN,
      4. et al
      . Diversity of the human intestinal microbial flora. Science 2005;308:1635–8.
      OpenUrlAbstract/FREE Full Text
    31. ↵
      1. Palmer C,
      2. Bik EM,
      3. Digiulio DB,
      4. et al
      . Development of the Human Infant Intestinal Microbiota. PLoS Biol 2007;5:e177.
      OpenUrlCrossRefPubMed
    32. ↵
      1. Turnbaugh PJ,
      2. Hamady M,
      3. Yatsunenko T,
      4. et al
      . A core gut microbiome in obese and lean twins. Nature 2009;457:480–4.
      OpenUrlCrossRefPubMedWeb of Science
    33. ↵
      1. Vaishampayan PA,
      2. Kuehl JV,
      3. Froula JL,
      4. et al
      . Comparative metagenomics and population dynamics of the gut microbiota in mother and infant. Genome Biol Evol 2010;2010:53–66.
      OpenUrl
    34. ↵
      1. Adlerberth I,
      2. Strachan DP,
      3. Matricardi PM,
      4. et al
      . Gut microbiota and development of atopic eczema in 3 European birth cohorts. J Allergy Clin Immunol 2007;120:343–50.
      OpenUrlCrossRefPubMedWeb of Science
    35. ↵
      1. Maukonen J,
      2. Simoes C,
      3. Saarela M
      . The currently used commercial DNA-extraction methods give different results of clostridial and actinobacterial populations derived from human fecal samples. FEMS Microbiol Ecol 2012;79:697–708.
      OpenUrlAbstract/FREE Full Text
    36. ↵
      1. Bahl MI,
      2. Bergstrom A,
      3. Licht TR
      . Freezing fecal samples prior to DNA extraction affects the Firmicutes to Bacteroidetes ratio determined by downstream quantitative PCR analysis. FEMS Microbiol Lett 2012;329:193–7.
      OpenUrlAbstract/FREE Full Text
    37. ↵
      1. Rotimi VO,
      2. Duerden BI
      . The development of the bacterial flora in normal neonates. J Med Microbiol 1981;14:51–62.
      OpenUrlAbstract/FREE Full Text
    38. ↵
      1. Singh N,
      2. Leger MM,
      3. Campbell J,
      4. et al
      . Control of vancomycin-resistant enterococci in the neonatal intensive care unit. Infect Control Hosp Epidemiol 2005;26:646–9.
      OpenUrlCrossRefPubMed
    39. ↵
      1. Martin R,
      2. Langa S,
      3. Reviriego C,
      4. et al
      . Human milk is a source of lactic acid bacteria for the infant gut. J Pediatr 2003;143:754–8.
      OpenUrlCrossRefPubMedWeb of Science
    40. ↵
      1. Lloyd CM,
      2. Hawrylowicz CM
      . Regulatory T cells in asthma. Immunity 2009;31:438–49.
      OpenUrlCrossRefPubMed
    41. ↵
      1. Vuillermin PJ,
      2. Ponsonby AL,
      3. Saffery R,
      4. et al
      . Microbial exposure, interferon gamma gene demethylation in naive T-cells, and the risk of allergic disease. Allergy 2009;64:348–53.
      OpenUrlCrossRefPubMedWeb of Science
    42. ↵
      1. Björkstén B
      . The intrauterine and postnatal environments. J Allergy Clin Immunol 1999;104:1119–27.
      OpenUrlCrossRefPubMedWeb of Science
    43. ↵
      1. Prescott SL
      . Early origins of allergic disease: a review of processes and influences during early immune development. Curr Opin Allergy Clin Immunol 2003;3:125–32.
      OpenUrlCrossRefPubMed
    44. ↵
      1. Mazmanian SK,
      2. Liu CH,
      3. Tzianabos AO,
      4. et al
      . An immunomodulatory molecule of symbiotic bacteria directs maturation of the host immune system. Cell 2005;122:107–18.
      OpenUrlCrossRefPubMedWeb of Science
    45. ↵
      1. Mazmanian SK,
      2. Round JL,
      3. Kasper DL
      . A microbial symbiosis factor prevents intestinal inflammatory disease. Nature 2008;453:620–5.
      OpenUrlCrossRefPubMedWeb of Science
    46. ↵
      1. Kelly D,
      2. Campbell JI,
      3. King TP,
      4. et al
      . Commensal anaerobic gut bacteria attenuate inflammation by regulating nuclear-cytoplasmic shuttling of PPAR-gamma and RelA. Nat Immunol 2004;5:104–12.
      OpenUrlCrossRefPubMedWeb of Science
    47. ↵
      1. Eggesbö M,
      2. Moen B,
      3. Peddada S,
      4. et al
      . Development of gut microbiota in infants not exposed to medical interventions. APMIS 2011;119:17–35.
      OpenUrlCrossRefPubMed
    48. ↵
      1. Yatsunenko T,
      2. Rey FE,
      3. Manary MJ,
      4. et al
      . Human gut microbiome viewed across age and geography. Nature 2012;486:222–7.
      OpenUrlCrossRefPubMedWeb of Science
    49. ↵
      1. Turroni F,
      2. Peano C,
      3. Pass DA,
      4. et al
      . Diversity of Bifidobacteria within the Infant Gut Microbiota. PLoS ONE 2012;7:e36957.
      OpenUrlCrossRefPubMed
    50. ↵
      1. Tannock GW,
      2. Fuller R,
      3. Smith SL,
      4. et al
      . Plasmid profiling of members of the family Enterobacteriaceae, lactobacilli, and bifidobacteria to study the transmission of bacteria from mother to infant. J Clin Microbiol 1990;28:1225–8.
      OpenUrlAbstract/FREE Full Text
    51. ↵
      1. Koenig JE,
      2. Spor A,
      3. Scalfone N,
      4. et al
      . Succession of microbial consortia in the developing infant gut microbiome. Proc Natl Acad Sci USA 2011;108(Suppl 1):4578–85.
      OpenUrlAbstract/FREE Full Text
    52. ↵
      1. Salam MT,
      2. Margolis HG,
      3. McConnell R,
      4. et al
      . Mode of delivery is associated with asthma and allergy occurrences in children. Ann Epidemiol 2006;16:341–6.
      OpenUrlCrossRefPubMedWeb of Science
    53. ↵
      1. Thavagnanam S,
      2. Fleming J,
      3. Bromley A,
      4. et al
      . A meta-analysis of the association between Caesarean section and childhood asthma. Clin Exp Allergy 2008;38:629–33.
      OpenUrlCrossRefPubMedWeb of Science
    View Abstract

    Supplementary materials

    Footnotes

    • Contributors Conception and design: HEJ, TRA, MCJ, BB, CJ, LE, AFA. Analysis and interpretation of data: HEJ, TRA, MCJ, KH, CQ, AFA. Drafting the article: HEJ, TRA, MCJ, BB, AFA. Final approval of submitted version: HEJ, TRA, MCJ, KH, CQ, CJ, BB, LE, AFA.

    • Funding This work was supported by the Ekhaga Foundation and the Söderbergs Foundation to LE and by the Swedish Research Council (K2011-56X-21854-01-06), the Research Council for the South-East Sweden, the Swedish Asthma and Allergy Association, the Olle Engkvist Foundation, the Vårdal Foundation–for Health Care Sciences and Allergy Research to MJ and by the Swedish Research Councils VR (621-2011-5689) and FORMAS (215-2009-1174) to AA. KH was funded through a direct grant from Unilever.

    • Competing interests None.

    • Ethics approval The human ethic committee at Linköping University, Linköping, Sweden.

    • Provenance and peer review Not commissioned; externally peer reviewed.