ORIGINAL RESEARCH article
Sec. Systems Microbiology
Volume 12 - 2021 | https://doi.org/10.3389/fmicb.2021.691312
The Influence of Cesarean Section on the Composition and Development of Gut Microbiota During the First 3 Months of Life
- Department of Gastroenterology, Children’s Hospital, Zhejiang University School of Medicine, National Clinical Research Center for Child Health, National Children’s Regional Medical Center, Zhejiang, China
The intestinal microbiota has emerged as a critical regulator of growth and development in the early postnatal period of life. Cesarean section (CS) delivery is one of the strongest disrupting factors of the normal colonization process and has been reported as a risk factor for disorders in later life. In this study, we dynamically and longitudinally evaluated the impact of CS on the initial colonization pattern and development of gut microbiota by 16 healthy Chinese infants with fecal samples collected at 9 time points (day 5, day 8, day 11, week 2, week 4, week 6, week 7, month 2, and month 3) during the first 3 months of life. The V3–V4 regions of 16S rRNA gene were analyzed by Illumina sequencing. In comparison with vaginally delivered (VD) infants, infants born by CS showed decreased relative abundance of Bacteroides and Parabacteroides and enrichment of Clostridium_sensu_stricto_1, Enterococcus, Klebsiella, Clostridioides, and Veillonella. Most interestingly, Firmicutes/Bacteroidetes ratio was found to be significantly higher in the CS group than in the VD group from day 5 until month 3. Besides, the results of microbial functions showed that the VD group harbored significantly higher levels of functional genes in vitamin B6 metabolism at day 5, day 8, week 2, week 4, week 6, week 7, month 2, and month 3 and taurine and hypotaurine metabolism at day 5, while the phosphotransferase system and starch and sucrose metabolism involved functional genes were plentiful in the CS group at day 11, week 2, week 4, week 6, week 7, and month 2 and at week 2, week 7, and month 2, respectively. Our results establish a new evidence that CS affected the composition and development of gut microbiota in the first 3 months and provide a novel insight into strategies for CS-related disorders in later life.
Accumulating evidence has shown that the gut microbiota plays a fundamental role in the health and disease by assisting in the synthesis and absorption of nutrients, strengthening gut integrity, protecting against enteropathogens, modulating the immune system, and exerting control over the gut–brain axis (Bäumler and Sperandio, 2016; Adak and Khan, 2019; Ahern and Maloy, 2020). The gastrointestinal tract contains 1014 microorganisms, which encompasses 10 times more bacterial cells than the total number of human cells (Thursby and Juge, 2017). The gut is of particular relevance to human health, as it contains the majority and most diverse set of human commensal bacteria (Savage, 1977). In the gut microbiota, bacteria are an essential part, and they are also the most commonly studied. The 16S rRNA high-throughput sequencing has become a popular technology for detecting intestinal flora in recent years due to its characteristics of fast sequencing, detailed results, and high accuracy (Fouhy et al., 2012). It is well known that the dynamic infant gut microbiota symbiosis is established from birth and developed rapidly in the first 3 years of life, particularly during infancy (Robertson et al., 2019; Shao et al., 2019). The early postnatal period of life is a vital stage for growth and development. The intestinal microbiota has emerged as a critical regulator of growth and development in the early postnatal period of life (Lee et al., 2021). The intestinal microbiota is deemed to coincide with the development of our body systems (Cowan et al., 2020; Wilkins and Reimer, 2021). Accordingly, perturbations influenced by any detrimental factors in early life may cause adversely long-lasting consequences for host health (Klancic and Reimer, 2020), ranging from gastrointestinal diseases, such as inflammatory bowel disease (Zuo et al., 2018) and irritable bowel syndrome (Rincel and Darnaudéry, 2020), immune reactivity, and immunopathology (Al Nabhani and Eberl, 2020) to allergies, including asthma, food allergy, and atopic dermatitis (Zhuang et al., 2019; Rachid et al., 2021; Renz and Skevaki, 2021).
The maturation of gut microbiota is strongly influenced by both internal host properties and external factors, such as ethnicity, gestational age, delivery mode, feeding pattern, and antibiotic usage, of which cesarean section (CS) is a crucial one (Heiss and Olofsson, 2019; Kim et al., 2019; Shao et al., 2019). It is reported that CS delivery is one of the strongest disrupting factors of the normal colonization process (Korpela et al., 2020). The microbial population of neonates born vaginally resembles the maternal vagina and perianal microbes, while CS-delivered neonates predominantly acquire bacteria derived from the maternal skin and the surrounding environment (Dominguez-Bello et al., 2010). Maternal microbiota exerts a significant effect on the neonatal microbiome and contributes to regulating the development of offspring immunity, metabolism, brain function, and behavior (Jašarević and Bale, 2019; García-Mantrana et al., 2020). These advantage effects of maternal microbiota can be vertically transmitted to offspring through vaginal delivery (Mueller et al., 2016). Despite CS can be a life-saving intervention when medically indicated, this procedure can lead to short-term and long-term health effects for children with the absence of transmission from maternal microbiota to offspring (Sandall et al., 2018). There is growing evidence that CS delivery is a risk factor for the onset of several clinical diseases later in life including asthma (Black et al., 2015), atopic disease (Seo et al., 2015), allergies (Sandall et al., 2018), and obesity (Kuhle et al., 2015). Furthermore, CS delivery has been reported to increase the risk of respiratory infections in later childhood (Unger and Bogaert, 2017) and impair cognitive capabilities in school-age children (Polidano et al., 2017). The correlation between CS delivery and these increased risks may be interpreted by the subsequent disturbance in immune response regulation after the absence of contact with maternal vaginal microbes at birth (Neu and Rushing, 2011).
Overall, the first 3 months of early life is the first peak of growth and development. Given the effects of CS delivery on the developing microbial colonization process in infants may pose risks for morbidity in later life, it is indispensable to characterize early gut microbiota profiles based on the influence of delivery mode for elucidating the mechanisms underlying their relationship with affected diseases. However, few studies have focused on the influence of CS on the dynamic of gut microbiota, particularly on this critical period. Therefore, we longitudinally collected fecal samples from a cohort of 16 healthy full-term Chinese infants at 9 time points throughout the first 3 months of life and aimed to identify the effects of CS delivery on the composition of gut microbiota by high-throughput sequencing of the V3–V4 regions of 16S rRNA gene in such period. To predict the functional potential of intestinal flora and point out the direction for the subsequent research, the PICRUSt software was used based on its high accuracy (>90% in most cases) (Douglas et al., 2018).
Materials and Methods
Study Design and Participants
The study was approved by the Ethics Committee of the Children’s Hospital of Zhejiang University School of Medicine. Written informed consent was obtained from mothers before collecting samples. The longitudinal cohort study was conducted between November 6, 2017 and December 5, 2018. Healthy infant volunteers who met the following criteria were recruited in our study: (1) whose mothers worked and were residents in the city of Hangzhou, China and (2) full-term healthy infants born between gestational age of 37 and 41 weeks; any subjects containing one of the following conditions were excluded: (1) infants born from mothers who had any infection during the third trimester of pregnancy and chronic diseases, such as diabetes, hypertensive disorders, or autoimmune disease and (2) infants who presented any abnormal conditions, such as neonatal asphyxia, hyperbilirubinemia, and any other conditions that needed extensive care. Parents were educated for fecal sampling and storage. Fecal samples were collected by parents at 9 time points during the first 3 months of life (day 5, day 8, day 11, week 2, week 4, week 6, week 7, month 2, and month 3 after birth). Clinical data regarding both mothers and infants including maternal age, gravidity and parity history, mode of delivery, sex, gestational age, birth weight, and other factors, such as family smoking and pet keeping, were recorded. A total of 16 infant volunteers were recruited in the longitudinal cohort study, including 6 of vaginally delivered (VD) and 10 of CS. 123 fecal samples were collected at 9 time points in the first 3 months of life.
Fecal Sample Collection
Fecal samples were collected into sterile containers by parents and immediately frozen in a −20°C refrigerator. Samples were then transported to the laboratory using dry ice by investigators and were frozen at −80°C for subsequent experimentation.
Genomic DNA Extraction
Genomic DNA was isolated from fecal samples using the QIAmp Fast DNA Stool Mini Kit (51604; Qiagen, Valencia, CA, United States) in accordance with the manufacturer’s instructions. DNA was eluted with 50 μl elution buffer and stored at −80°C before use. The DNA extract was checked on 1% agarose gel, and DNA concentration and purity were assessed by NanoDrop 2000 UV–vis spectrophotometer (Thermo Scientific, Wilmington, DE, United States).
Sequencing of 16S rRNA Gene
The sequencing of 16S rRNA gene was performed as previously reported method (Zheng et al., 2019, 2020). Concretely, the V3–V4 hypervariable regions of the bacteria 16S rRNA gene were amplified with barcode-indexed primers 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) by an ABI GeneAmp® 9700 PCR thermocycler (ABI, Foster City, CA, United States). The PCR amplification mixture contained: 4 μl 5 × TransStart FastPfu Buffer, 2 μl 2.5 mM dNTPs, 0.8 μl forward primer (5 μM), 0.8 μl reverse primer (5 μM), 0.4 μl TransStart FastPfu DNA Polymerase, 0.2 μl BSA, 10 ng template DNA, and finally ddH2O up to 20 μl. The PCR conditions for amplification were as follows: initial denaturation for 3 min at 95°C, followed by 27 cycles of denaturation at 95°C for 30 s, annealing at 55°C for 30 s and extension at 72°C for 45 s, and single extension at 72°C for 10 min, and end at 10°C. PCR reactions were performed in triplicate. PCR products were extracted from 2% agarose gel, purified by the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, United States), and quantified using QuantusTM Fluorometer (Promega, Madison, WI, United States) according to the manufacturer’s instructions. Purified amplicons were pooled in equimolar and paired-end sequenced on an Illumina MiSeq platform (Illumina, San Diego, CA, United States) according to the standard protocols by Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China).
Sequence Data Processing and Statistical Analysis
The protocol of sequence data processing was similar to the previously published study (Zheng et al., 2020). Briefly, the raw 16S rRNA gene sequencing reads were demultiplexed, quality-filtered by Trimmomatic (Bolger et al., 2014), and merged by FLASH (Magoc̀ and Salzberg, 2011) with the following criteria: (i) the 300 bp reads were truncated at any site receiving an average quality score of < 20 over a 50-bp sliding window, and the truncated reads shorter than 50 bp were discarded; reads containing ambiguous characters were also discarded; (ii) only overlapping sequences longer than 10 bp were assembled according to their overlapped sequence. The maximum mismatch ratio of the overlap region was 0.2. Reads that could not be assembled were discarded; (iii) samples were distinguished according to the barcode and primers, and the sequence direction was adjusted. The allowed maximum barcode mismatch number was 0, and the primer mismatch number was 2. Operational taxonomic units (OTUs) with 97% similarity cutoff were clustered using UPARSE (Edgar, 2013) (version 7.1,1), and chimeric sequences were identified and removed. To minimize the effects of sequencing depth on alpha and beta diversity measures, the number of reads from each sample was rarefied to 29,101, which still yielded an average Good’s coverage of 99.60%. The taxonomy of each OTU representative sequence was analyzed by RDP Classifier (Wang et al., 2007)2 against the 16S rRNA database (Silva 132) using the confidence threshold of 0.7. Taxonomic relative abundance profiles at taxa levels (domain, kingdom, phylum, class, order, family, and genus) were generated based on OTU annotation. Principal coordinates analysis (PCoA) was performed based on Weighted Unifrac distance to visualize differences in bacterial structure of infants under delivery modes and sampling time points. Before Unifrac analysis, the phylogenetic tree was constructed with representative OTU sequences by FastTree version 2.1.3 and MUSCLE (multiple sequence comparison by log-expectation). Differences in microbial composition and abundance between groups of PCoA were tested using the analysis of similarities (ANOSIM) on the beta diversity matrix. Generalized Linear Models with participant identity as a random effect were used to determine the variation tendency of gut microbial abundance with age. And the Wilcoxon rank-sum test was used to determine the significant differences between groups at both phylum and genus levels for abundance of gut microbiota, and in the microbial functional gene composition. The microbiota enterotype analysis was performed based on R (version 3.3.1) combined with Jensen–Shannon divergence (JSD) distance to evaluate the distinct growth pattern of gut microbiota according to delivery mode. To predict the microbial functions, functional prediction and annotations were conducted using the Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) (Douglas et al., 2018). The distribution of data, such as infant information (gestational age and birth weight), maternal information (age, gravidity, parity, pre-pregnancy weight, and pregnancy weight), and Firmicutes/Bacteroidetes (F/B) ratio, were analyzed by Kolmogorov–Smirnov test. Thereinto, normally distributed data (gestational age, birth weight, maternal age, gravidity, parity, pre-pregnancy weight, and pregnancy weight) were expressed as mean ± SD, and data between two groups were compared using the Student’s t-test. Non-normally distributed data (F/B ratio) were presented as median (minimum, maximum) unless specifically highlighted, and data between two groups were compared using the Mann–Whitney U test. Countable data, such as infant sex and other information (family smoking, pet keeping), were presented as number (%) and tested using chi-square test.
General Characteristics of the Subjects
The flow diagram of the study is generalized in Figure 1. The general characteristics of the study population according to delivery mode are summarized in Table 1. There were no significant differences in the baseline characteristics between VD and CS-delivered infants, such as infant information (gestational age and birth weight) and maternal information (age, gravidity, parity, pre-pregnancy weight, and pregnancy weight) (Student’s t-test, all p > 0.05), as well as infant information (sex) and other information (family smoking and keeping pets) (chi-square test, all p > 0.05).
Figure 1. Flow diagram of the study. The diagram showed the fecal samples collected at 9 time points. x, fecal samples were not collected by parents as the infants had no defecation.
Comparison of Gut Microbial Diversity
The rarefaction curves showed that curves have reached a flat level, indicating that the sequencing depth was sufficient (Supplementary Figure 1). The Chao1 (Figure 2A) and Shannon (Figure 2B) indexes, which represent microbial alpha diversity and microbiota richness, showed non-coordinated dynamics for all individuals in the infants both born with VD and born with CS. The results showed that there was no statistical difference both in the Chao1 (Mann–Whitney U test, all p > 0.05; Figure 2C) and in the Shannon index (Mann–Whitney U test, all p > 0.05; Figure 2D) at 9 time points between the two groups.
Figure 2. Alpha diversity dynamic of the gut microbiota based on 9 time points between infants born with vaginal delivery and born with cesarean section. (A) Chao1 index; (B) Shannon index; (C) the comparison of Chao1 index between VD and CS at 9 time points (Mann–Whitney U test, all p > 0.05); (D) the comparison of Shannon index between VD and CS at 9 time points (Mann–Whitney U test, all p > 0.05). The results were reported as median values (minimum–maximum). VD, vaginal delivery; CS, cesarean section.
Impact of CS on Temporal Changes of Gut Microbial Composition During the First 3 Months of Life
In two delivery mode groups, we determined the most predominant 4 phyla for the first 3 months of life: Firmicutes, Proteobacteria, Bacteroidetes, and Actinobacteria (Figure 3A). Despite the 4 dominant phyla were present in both groups, we observed a clear distinction of temporal changes in their proportions at the phylum level between the two groups. Within the CS group, Firmicutes was observed as the most dominant phylum from day 5 to week 4 after birth. The relative abundance of Firmicutes gradually decreased from day 5 to month 3 after birth (Supplementary Table 1, Generalized Linear Models, p = 0.000). In contrast, the prevalence of Proteobacteria continuously increased from day 5 to month 3 after birth (Generalized Linear Models, p = 0.006) and dominated from week 6 to month 3 in the CS group. Additionally, Bacteroidetes had a lower abundance from day 5 to week 4 and then increased gradually from week 6 to month 3 in the CS group. Whereas in the VD group, the predominant phyla from day 5 to day 11 were Proteobacteria and Firmicutes. The relative proportion of Firmicutes significantly decreased from day 5 to month 3 (Generalized Linear Models, p = 0.021). At the genus level, several genera demonstrated differences in variation tendency between two delivery mode groups (Figure 3B). For instance, the relative abundance of Escherichia–Shigella displayed an increasing trend from day 5 to month 3 in CS infants (Supplementary Table 2, Generalized Linear Models, p = 0.001). Clostridioides and Streptococcus showed the decreasing tendency in abundance from day 5 to month 3 in the CS group (Generalized Linear Models, p < 0.01), whereas genera in the VD group did not show obvious variation tendency.
Figure 3. Dynamic development of microbiota in the vaginal delivery and cesarean section groups. (A) Dynamic development of microbiota in phylum level; (B) dynamic development of microbiota in genus level; CS, cesarean section; VD, vaginal delivery.
Impact of CS on Temporal Change of F/B Ratio
Cesarean section delivery significantly affected F/B ratio from day 5 to month 3 compared with VD delivery (Figure 4, Mann–Whitney U test, all p < 0.05). F/B ratio was significantly higher in the CS group than in the VD group at every time point (Figure 4, Mann–Whitney U test, all p < 0.05).
Figure 4. Cesarean section significantly increased Firmicutes/Bacteroidetes ratio in the first 3 months of life. Data were representative with median. F/B ratio: Firmicutes/Bacteroidetes ratio; CS: cesarean section; VD: vaginal delivery. *p < 0.05; **p < 0.01 (Mann–Whitney U test).
Comparison of the Gut Microbiota Composition by the Mode of Delivery
The overall difference in microbiota community between two groups (beta diversity) was compared using PCoA based on Weighted Unifrac distance combined with ANOSIM, which distinctly revealed the effect of delivery mode on the gut microbiota at 9 time points (Figure 5, ANOSIM, all p < 0.05). The PCoA plot showed differences in bacterial phylogenetic structure between the CS and VD groups at every time point during the first 3 months in life. And the ANOSIM results revealed that such differences between two groups at every time point were statistically significant (Figure 5, ANOSIM, all p < 0.05). Thereafter, the significance test of the difference between groups was conducted by Wilcoxon rank-sum test at both phylum and genus levels. Concerning the phylum level, Firmicutes was found to be significantly more abundant in the CS group than in the VD group at day 5 and week 2, while Bacteroidetes exhibited significant superiority in the VD group at day 5, day 8, day 11, week 2, week 4, week 6, week 7, and month 2 (Figure 6, Wilcoxon rank-sum test, p < 0.05). Concerning the genus level, the microbiota of the VD group was characterized by significantly high levels of Escherichia–Shigella at day 5, day 8, day 11, and week 2 and persistent high levels of Bacteroides at 9 time points in comparison with the CS group. Besides, Phascolarctobacterium and Parabacteroides were also significantly more enriched in VD infants at day 5 and week 7 and at day 11, week 4, week 7, and month 3, respectively. We also observed that the VD group had a higher relative proportion of Collinsella at week 6 and week 7 than the CS group. As for infants delivered by CS, Clostridium_sensu_stricto_1 (at day 5, day 11, week 2, week 4, and week 7), Enterococcus (at day 5, week 2, week 4, and week 6), Klebsiella (week 2, week 6, and week 7), Clostridioides (week 6, month 3), Citrobacter (week 6, month 3), Streptococcus (month 2), Veillonella (week 6), and Enterobacter (week 7) displayed significantly higher relative abundance (Figure 7, Wilcoxon rank-sum test, p < 0.05).
Figure 5. Principal coordinates analysis of the infant’s gut microbiota among different time points based on delivery mode. R and p were calculated using ANOSIM.
Figure 6. Cesarean section significantly affected the development of gut microbiota at the phylum level in the first 3 months of life. *p < 0.05; **p < 0.01 (Wilcoxon rank-sum test).
Figure 7. Cesarean section significantly affected the development of gut microbiota at the genus level in the first 3 months of life. *p < 0.05; **p < 0.01 (Wilcoxon rank-sum test).
The Microbiota Enterotype Analysis Based on Delivery Mode
To evaluate the distinct growth pattern of gut microbiota according to delivery mode, the microbiota enterotype was analyzed at specific time points. As shown in Figure 8, enterotypes mainly represented by Clostridioides and Clostridium_sensu_stricto_1 were predominant in the CSgroup at day 5, while the VD group possessed enterotype dominated by Escherichia–Shigella. In addition, Clostridioides and Escherichia–Shigella-dominated enterotypes were abundant in the CS and VD groups at day 11, respectively. Furthermore, we found the enrichment of Clostridium_sensu_stricto_1 and Bacteroides-represented enterotypes at week 4 in the CS and VD groups, respectively. Besides, the CS group was predominated by Klebsiella and Clostridium_sensu_stricto_1-represented enterotypes, whereas Bifidobacterium-dominated enterotype was pronounced in VD infants at week 6.
Comparison of Microbial Functions According to the Delivery Mode
Changes in microbial functions by different delivery modes using PICRUSt were predicted. We found significant differences in the functional gene composition in the metabolic pathway between two delivery mode groups (Supplementary Figures 2–4, Wilcoxon rank-sum test, p < 0.05). The results showed that the VD group harbored significantly higher levels of functional genes in vitamin B6 metabolism at day 5, day 8, week 2, week 4, week 6, week 7, month 2, and month 3 and taurine and hypotaurine metabolism at day 5 (Supplementary Figures 2–4, Wilcoxon rank-sum test, p < 0.05), whereas the phosphotransferase system (PTS) and starch and sucrose metabolism involved functional genes were plentiful in the CS group at day 11, week 2, week 4, week 6, week 7, and month 2 and at week 2, week 7, and month 2, respectively (Supplementary Figures 2–4, Wilcoxon rank-sum test, p < 0.05).
In this study, we explored the effects of CS on the development of gut microbiota in healthy Chinese infants by longitudinally collected fecal samples analyzed by 16S rRNA gene sequencing during the first 3 months of life. Our results demonstrated the significant perturbation of the early colonization of gut microbiota caused by the CS compared with VD. Concordant with previous reports, we found that infants delivered by CS exhibited a decreased colonization of Bacteroidetes and an increased abundance of the Firmicute phylum (Jakobsson et al., 2014; Rutayisire et al., 2016; Yang R. et al., 2019), which might be interpreted by the lack of contact with maternal vaginal microbiota during childbirth (Rutayisire et al., 2016). Of particular importance, the colonization with Bacteroides was evidently delayed in infants born by CS, which was in line with previous findings (Hesla et al., 2014; Lee et al., 2016; Rutayisire et al., 2016; Kim et al., 2019; Shao et al., 2019). Bacteroides genus was thought to be capable of activating T-cell-dependent immune responses, thereby modulating the development and homeostasis of the host immune system (Troy and Kasper, 2010). Coincidentally, the delayed Bacteroidetes colonization and reduced type 1 of helper T lymphocyte (Th1) responses were illustrated in infants delivered by CS (Jakobsson et al., 2014). Furthermore, Bacteroides has been reported as a major component of gut microbiota in adults (Eckburg et al., 2005). In this regard, the presence of the maternal transmission in gut microbiota during delivery might contribute to the high abundance of Bacteroides in infants born vaginally (Vaishampayan et al., 2010; Lee et al., 2016). This indicated the essential role of maternal transmission during vaginal delivery in the establishment of early gut microbiota in infants.
Most interestingly, we found that the F/B ratio was higher in the CS group than in the VD group from day 5 until month 3. The F/B ratio was widely accepted to have an important influence in maintaining normal intestinal homeostasis. Increased F/B ratio was regarded as dysbiosis, and was usually observed with obesity (Stojanov et al., 2020). A previous systematic review and meta-analysis revealed that children born by CS were at higher risk of developing obesity in childhood, despite that findings were limited by a moderate heterogeneity among studies and the potential for residual confounding and publication bias (Kuhle et al., 2015). Moreover, a recent study from Mexico showed that infants born by CS had a strikingly high F/B ratio (Murata et al., 2020). High F/B ratio in CS delivery infants may help explain why Mexico that had one of the highest global average annual rate increases in CS led the world in childhood obesity, which may be associated with the high F/B ratio in CS delivery infants. Also, another recent study from Mexico demonstrated that Firmicutes and Bacteroidetes phyla were positively and negatively associated with obesity, respectively (Vazquez-Moreno et al., 2020). By intensively monitoring the F/B ratio at 9 time points, the CS group had a permanently higher F/B ratio than the VD group that might provide crucial evidence for the causality of high F/B ratio in the first 3 months of life and the obesity in the later life. As one study has demonstrated that mice in the CS group gained more body mass after weaning than those in the VD group and established a causal relationship between CS and increased body weight (Martinez et al., 2017). Despite CS delivery infants in our study seems to have a similar body weight with VD delivery infants at the short-term observation, it is necessary to implement long-term follow-up and determine the relationship between higher F/B ratio of early life and the risk for obesity of later life. In parallel, our results also provide new insight into the prevention and management of obesity through early regulation of the normal F/B ratio by probiotics usage or other treatments, such as fecal bacteria transplantation (Korpela et al., 2020). However, some studies revealed no significant differences in phyla abundances or F/B ratios observed between normal-weight and obese children (López-Contreras et al., 2018; Chen et al., 2020). Therefore, more studies are required to confirm the findings.
In addition, in agreement with several studies performed before (Liu et al., 2015; Lee et al., 2016; Hill et al., 2017), we observed a significantly higher level of Clostridium in the CS group, which has been reported to be ascribed to nosocomial infection (Penders et al., 2006). What is more, the correlation between the colonization of Clostridium in the gut and the increased risks of wheeze and eczema has been proven (van Nimwegen et al., 2011). CS delivery also elevated the abundances of potential pathogenic and proinflammatory bacteria associated with the hospital environment including Klebsiella, Enterococcus, and Enterobacter in infants, which were in accordance with previous studies (Liu et al., 2015; Reyman et al., 2019; Shao et al., 2019; Yang R. et al., 2019). And the increased proportion of Klebsiella in early life has been validated to be connected with later development of allergy in children (Low et al., 2017). Also, Klebsiella, as a common opportunistic pathogenic bacterium in hospitals, is usually associated with nosocomial infection and acts as a reservoir for a myriad of antimicrobial resistance genes (Rosenthal et al., 2012; Gibson et al., 2016), which might suggest the possible correlation between delivery mode and the increased risk of antibiotic-resistant infections (Reyman et al., 2019). In our study, Streptococcus, as a member of maternal skin bacteria, was more often colonized in CS-delivered infants, as other studies reported (Bäckhed et al., 2015; Moore and Townsend, 2019), while Escherichia–Shigella has also been reported to be more common in VD infants before (Bäckhed et al., 2015; Yang R. et al., 2019). Besides, Parabacteroides was significantly overrepresented in VD infants, which was also discovered in previous findings (Bäckhed et al., 2015; Hill et al., 2017). The lower relative abundance of Parabacteroides was found to be associated with a high risk for eczema in infants (Wopereis et al., 2018). Nevertheless, infants born by CS showed significantly higher relative levels of Veillonella in our results, which was consistent with other reports (Hesla et al., 2014; Bäckhed et al., 2015). Moreover, Veillonella was ascertained to be relevant to the penetrating complications of Crohn’s disease in children (Kugathasan et al., 2017). Based on the above mentioned, we conclude that gut microbial dysbiosis and delayed maturation are correlated with CS.
Our functional prediction revealed that delivery mode might have a certain impact on microbial functions. It has also been previously reported that the functional genes of taurine and hypotaurine metabolism were significantly higher in the VD group, while PTS was enriched in the CS group (Yang B. et al., 2019). Taurine was bound to bile acid to form taurocholic acid, which was involved in the fat emulsification process, for promoting the digestion of mother’s milk fat, while PTS and starch and sucrose metabolism were involved in the digestion and absorption of carbohydrates. This indicated that the gut microbiota in different delivery groups might associate with the metabolism of diverse nutrients.
The limits of this study were the small sample size and short follow-up period. Moreover, we did not collect all the fecal samples at different time points because of no defecation of infants. In spite of this, we took intensive sampling within the first 3 months of life, and the total 9 time points provided the dynamic evaluation of the early gut microflora during this early life period. Our results decipher the crucial role of delivery mode in shaping the early gut microbiota of infants.
In conclusion, CS delivery mode significantly perturbs the early colonization and development of gut microbiota. Furthermore, CS causes significantly higher F/B ratio in the first 3 months of life, which may be associated with obesity in later life. Importantly, our findings establish new evidence that CS affected the composition and development of gut microbiota in the first 3 months and provide novel insight into strategies for CS-related disorders in later life.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: SRA, SRP283175.
The studies involving human participants were reviewed and approved by Ethics Committee of the Children’s Hospital of Zhejiang University School of Medicine. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin. Written informed consent was obtained from the individual(s), and minor(s)’ legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
GL and MJ designed the present study and drafted the manuscript. GL, YH, XS, and WZ performed the literature search. GL, ET, and BC collected the samples and performed the analyses. GL, ET, and MJ edited the manuscript. All authors have read and approved the manuscript.
This research was funded by grants from the Biostime International Holdings Limited Maternal and Infant Nutrition and Health Research Projects of the National Center for Women and Children’s Health Chinese Center for Disease Control and Prevention (No. 2018FYH011).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.691312/full#supplementary-material
Adak, A., and Khan, M. R. (2019). An insight into gut microbiota and its functionalities. Cell. Mol. Life Sci. 76, 473–493. doi: 10.1007/s00018-018-2943-4
Ahern, P. P., and Maloy, K. J. (2020). Understanding immune-microbiota interactions in the intestine. Immunology 159, 4–14. doi: 10.1111/imm.13150
Al Nabhani, Z., and Eberl, G. (2020). Imprinting of the immune system by the microbiota early in life. Mucosal Immunol. 13, 183–189. doi: 10.1038/s41385-020-0257-y
Bäckhed, F., Roswall, J., Peng, Y., Feng, Q., Jia, H., Kovatcheva-Datchary, P., et al. (2015). Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe 17:852. doi: 10.1016/j.chom.2015.05.012
Bäumler, A. J., and Sperandio, V. (2016). Interactions between the microbiota and pathogenic bacteria in the gut. Nature 535, 85–93. doi: 10.1038/nature18849
Black, M., Bhattacharya, S., Philip, S., Norman, J. E., and McLernon, D. J. (2015). Planned cesarean delivery at term and adverse outcomes in childhood health. JAMA 314, 2271–2279. doi: 10.1001/jama.2015.16176
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics (Oxford, England) 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Chen, X., Sun, H., Jiang, F., Shen, Y., Li, X., Hu, X., et al. (2020). Alteration of the gut microbiota associated with childhood obesity by 16S rRNA gene sequencing. PeerJ 8:e8317. doi: 10.7717/peerj.8317
Cowan, C. S. M., Dinan, T. G., and Cryan, J. F. (2020). Annual research review: critical windows - the microbiota-gut-brain axis in neurocognitive development. J. Child Psychol. Psychiatry 61, 353–371. doi: 10.1111/jcpp.13156
Dominguez-Bello, M. G., Costello, E. K., Contreras, M., Magris, M., Hidalgo, G., Fierer, N., et al. (2010). Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc. Natl. Acad. Sci. U.S.A. 107, 11971–11975. doi: 10.1073/pnas.1002601107
Douglas, G. M., Beiko, R. G., and Langille, M. G. I. (2018). Predicting the functional potential of the microbiome from marker genes using PICRUSt. Methods Mol. Biol. 1849, 169–177. doi: 10.1007/978-1-4939-8728-3_11
Eckburg, P. B., Bik, E. M., Bernstein, C. N., Purdom, E., Dethlefsen, L., Sargent, M., et al. (2005). Diversity of the human intestinal microbial flora. Science (New York, N.Y.) 308, 1635–1638. doi: 10.1126/science.1110591
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Fouhy, F., Ross, R. P., Fitzgerald, G. F., Stanton, C., and Cotter, P. D. (2012). Composition of the early intestinal microbiota: knowledge, knowledge gaps and the use of high-throughput sequencing to address these gaps. Gut Microbes 3, 203–220. doi: 10.4161/gmic.20169
García-Mantrana, I., Selma-Royo, M., González, S., Parra-Llorca, A., Martínez-Costa, C., and Collado, M. C. (2020). Distinct maternal microbiota clusters are associated with diet during pregnancy: impact on neonatal microbiota and infant growth during the first 18 months of life. Gut microbes 11, 962–978. doi: 10.1080/19490976.2020.1730294
Gibson, M. K., Wang, B., Ahmadi, S., Burnham, C.-A. D., Tarr, P. I., Warner, B. B., et al. (2016). Developmental dynamics of the preterm infant gut microbiota and antibiotic resistome. Nat. Microbiol. 1:16024. doi: 10.1038/nmicrobiol.2016.24
Heiss, C. N., and Olofsson, L. E. (2019). The role of the gut microbiota in development, function and disorders of the central nervous system and the enteric nervous system. J. Neuroendocrinol. 31:e12684. doi: 10.1111/jne.12684
Hesla, H. M., Stenius, F., Jäderlund, L., Nelson, R., Engstrand, L., Alm, J., et al. (2014). Impact of lifestyle on the gut microbiota of healthy infants and their mothers—the ALADDIN birth cohort. FEMS Microbiol. Ecol. 90, 791–801. doi: 10.1111/1574-6941.12434
Hill, C. J., Lynch, D. B., Murphy, K., Ulaszewska, M., Jeffery, I. B., O’Shea, C. A., et al. (2017). Evolution of gut microbiota composition from birth to 24 weeks in the INFANTMET Cohort. Microbiome 5:4. doi: 10.1186/s40168-016-0213-y
Jakobsson, H. E., Abrahamsson, T. R., Jenmalm, M. C., Harris, K., Quince, C., Jernberg, C., et al. (2014). Decreased gut microbiota diversity, delayed Bacteroidetes colonisation and reduced Th1 responses in infants delivered by caesarean section. Gut 63, 559–566. doi: 10.1136/gutjnl-2012-303249
Jašarević, E., and Bale, T. L. (2019). Prenatal and postnatal contributions of the maternal microbiome on offspring programming. Front. Neuroendocrinol. 55:100797. doi: 10.1016/j.yfrne.2019.100797
Kim, H., Sitarik, A. R., Woodcroft, K., Johnson, C. C., and Zoratti, E. (2019). Birth mode, breastfeeding, pet exposure, and antibiotic use: associations with the gut microbiome and sensitization in children. Curr. Allergy Asthma Rep. 19:22. doi: 10.1007/s11882-019-0851-9
Klancic, T., and Reimer, R. A. (2020). Gut microbiota and obesity: impact of antibiotics and prebiotics and potential for musculoskeletal health. J. Sport Health Sci. 9, 110–118. doi: 10.1016/j.jshs.2019.04.004
Korpela, K., Helve, O., Kolho, K.-L., Saisto, T., Skogberg, K., Dikareva, E., et al. (2020). Maternal fecal microbiota transplantation in cesarean-born infants rapidly restores normal gut microbial development: a proof-of-concept study. Cell 183, 324–334.e5. doi: 10.1016/j.cell.2020.08.047
Kugathasan, S., Denson, L. A., Walters, T. D., Kim, M.-O., Marigorta, U. M., Schirmer, M., et al. (2017). Prediction of complicated disease course for children newly diagnosed with Crohn’s disease: a multicentre inception cohort study. Lancet (London, England) 389, 1710–1718. doi: 10.1016/S0140-6736(17)30317-3
Kuhle, S., Tong, O. S., and Woolcott, C. G. (2015). Association between caesarean section and childhood obesity: a systematic review and meta-analysis. Obes. Rev. 16, 295–303. doi: 10.1111/obr.12267
Lee, E., Kim, B. J., Kang, M. J., Choi, K. Y., Cho, H. J., Kim, Y., et al. (2016). Dynamics of gut microbiota according to the delivery mode in healthy Korean Infants. Allergy Asthma Immunol. Res. 8, 471–477. doi: 10.4168/aair.2016.8.5.471
Lee, Y. M., Mu, A., Wallace, M., Gengatharan, J. M., Furst, A. J., Bode, L., et al. (2021). Microbiota control of maternal behavior regulates early postnatal growth of offspring. Sci. Adv. 7:eabe6563. doi: 10.1126/sciadv.abe6563
Liu, D., Yu, J., Li, L., Ai, Q., Feng, J., Song, C., et al. (2015). Bacterial community structure associated with elective cesarean section versus vaginal delivery in Chinese newborns. J. Pediatr. Gastroenterol. Nutr. 60, 240–246. doi: 10.1097/MPG.0000000000000606
López-Contreras, B. E., Morán-Ramos, S., Villarruel-Vázquez, R., Macías-Kauffer, L., Villamil-Ramírez, H., León-Mimila, P., et al. (2018). Composition of gut microbiota in obese and normal-weight Mexican school-age children and its association with metabolic traits. Pediatr. Obes. 13, 381–388. doi: 10.1111/ijpo.12262
Low, J. S. Y., Soh, S. E., Lee, Y. K., Kwek, K. Y. C., Holbrook, J. D., Van der Beek, E. M., et al. (2017). Ratio of Klebsiella/Bifidobacterium in early life correlates with later development of paediatric allergy. Benef. Microbes 8, 681–695. doi: 10.3920/BM2017.0020
Magoc̀, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics (Oxford, England) 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Martinez, K. A. II, Devlin, J. C., Lacher, C. R., Yin, Y., Cai, Y., Wang, J., et al. (2017). Increased weight gain by C-section: functional significance of the primordial microbiome. Sci Adv 3:eaao1874. doi: 10.1126/sciadv.aao1874
Moore, R. E., and Townsend, S. D. (2019). Temporal development of the infant gut microbiome. Open Biol. 9:190128. doi: 10.1098/rsob.190128
Mueller, N. T., Shin, H., Pizoni, A., Werlang, I. C., Matte, U., Goldani, M. Z., et al. (2016). Birth mode-dependent association between pre-pregnancy maternal weight status and the neonatal intestinal microbiome. Sci. Rep. 6:23133. doi: 10.1038/srep23133
Murata, C., Gutiérrez-Castrellón, P., Pérez-Villatoro, F., García-Torres, I., Enríquez-Flores, S., de la Mora-de la Mora, I., et al. (2020). Delivery mode-associated gut microbiota in the first 3 months of life in a country with high obesity rates: a descriptive study. Medicine 99:e22442. doi: 10.1097/MD.0000000000022442
Neu, J., and Rushing, J. (2011). Cesarean versus vaginal delivery: long-term infant outcomes and the hygiene hypothesis. Clin. Perinatol. 38, 321–331. doi: 10.1016/j.clp.2011.03.008
Penders, J., Thijs, C., Vink, C., Stelma, F. F., Snijders, B., Kummeling, I., et al. (2006). Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics 118, 511–521. doi: 10.1542/peds.2005-2824
Polidano, C., Zhu, A., and Bornstein, J. C. (2017). The relation between cesarean birth and child cognitive development. Sci. Rep. 7:11483. doi: 10.1038/s41598-017-10831-y
Rachid, R., Stephen-Victor, E., and Chatila, T. A. (2021). The microbial origins of food allergy. J. Allergy Clin. Immunol. 147, 808–813. doi: 10.1016/j.jaci.2020.12.624
Renz, H., and Skevaki, C. (2021). Early life microbial exposures and allergy risks: opportunities for prevention. Nat. Rev. Immunol. 21, 177–191. doi: 10.1038/s41577-020-00420-y
Reyman, M., van Houten, M. A., van Baarle, D., Bosch, A. A. T. M., Man, W. H., Chu, M. L. J. N., et al. (2019). Impact of delivery mode-associated gut microbiota dynamics on health in the first year of life. Nat. Commun. 10:4997. doi: 10.1038/s41467-019-13014-7
Rincel, M., and Darnaudéry, M. (2020). Maternal separation in rodents: a journey from gut to brain and nutritional perspectives. Proc. Nutr. Soc. 79, 113–132. doi: 10.1017/S0029665119000958
Robertson, R. C., Manges, A. R., Finlay, B. B., and Prendergast, A. J. (2019). The human microbiome and child growth - first 1000 days and beyond. Trends Microbiol. 27, 131–147. doi: 10.1016/j.tim.2018.09.008
Rosenthal, V. D., Bijie, H., Maki, D. G., Mehta, Y., Apisarnthanarak, A., Medeiros, E. A., et al. (2012). International Nosocomial Infection Control Consortium (INICC) report, data summary of 36 countries, for 2004-2009. Am. J. Infect. Control 40, 396–407. doi: 10.1016/j.ajic.2011.05.020
Rutayisire, E., Huang, K., Liu, Y., and Tao, F. (2016). The mode of delivery affects the diversity and colonization pattern of the gut microbiota during the first year of infants’ life: a systematic review. BMC Gastroenterol. 16:86. doi: 10.1186/s12876-016-0498-0
Sandall, J., Tribe, R. M., Avery, L., Mola, G., Visser, G. H., Homer, C. S., et al. (2018). Short-term and long-term effects of caesarean section on the health of women and children. Lancet (London, England) 392, 1349–1357. doi: 10.1016/S0140-6736(18)31930-5
Savage, D. C. (1977). Microbial ecology of the gastrointestinal tract. Annu. Rev. Microbiol. 31, 107–133. doi: 10.1146/annurev.mi.31.100177.000543
Seo, J. H., Kim, H. Y., Jung, Y. H., Lee, E., Yang, S. I., Yu, H. S., et al. (2015). Interactions between innate immunity genes and early-life risk factors in allergic rhinitis. Allergy Asthma Immunol. Res. 7, 241–248. doi: 10.4168/aair.2015.7.3.241
Shao, Y., Forster, S. C., Tsaliki, E., Vervier, K., Strang, A., Simpson, N., et al. (2019). Stunted microbiota and opportunistic pathogen colonization in caesarean-section birth. Nature 574, 117–121. doi: 10.1038/s41586-019-1560-1
Stojanov, S., Berlec, A., and Štrukelj, B. (2020). The influence of probiotics on the Firmicutes/Bacteroidetes ratio in the treatment of obesity and inflammatory bowel disease. Microorganisms 8:1715. doi: 10.3390/microorganisms8111715
Thursby, E., and Juge, N. (2017). Introduction to the human gut microbiota. Biochem. J. 474, 1823–1836. doi: 10.1042/bcj20160510
Troy, E. B., and Kasper, D. L. (2010). Beneficial effects of Bacteroides fragilis polysaccharides on the immune system. Front. Biosci. (Landmark Ed.) 15, 25–34. doi: 10.2741/3603
Unger, S. A., and Bogaert, D. (2017). The respiratory microbiome and respiratory infections. J. Infect. 74(Suppl. 1), S84–S88. doi: 10.1016/S0163-4453(17)30196-2
Vaishampayan, P. A., Kuehl, J. V., Froula, J. L., Morgan, J. L., Ochman, H., and Francino, M. P. (2010). Comparative metagenomics and population dynamics of the gut microbiota in mother and infant. Genome Biol. Evol. 2, 53–66. doi: 10.1093/gbe/evp057
van Nimwegen, F. A., Penders, J., Stobberingh, E. E., Postma, D. S., Koppelman, G. H., Kerkhof, M., et al. (2011). Mode and place of delivery, gastrointestinal microbiota, and their influence on asthma and atopy. J. Allergy Clin. Immunol. 128, 948–55.e1–3. doi: 10.1016/j.jaci.2011.07.027
Vazquez-Moreno, M., Perez-Herrera, A., Locia-Morales, D., Dizzel, S., Meyre, D., Stearns, J. C., et al. (2020). Association of gut microbiome with fasting triglycerides, fasting insulin and obesity status in Mexican children. Pediatr. Obes. 16:e12748. doi: 10.1111/ijpo.12748
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/aem.00062-07
Wilkins, A. T., and Reimer, R. A. (2021). Obesity, early life gut microbiota, and antibiotics. Microorganisms 9:413. doi: 10.3390/microorganisms9020413
Wopereis, H., Sim, K., Shaw, A., Warner, J. O., Knol, J., and Kroll, J. S. (2018). Intestinal microbiota in infants at high risk for allergy: effects of prebiotics and role in eczema development. J. Allergy Clin. Immunol. 141, 1334–1342.e5. doi: 10.1016/j.jaci.2017.05.054
Yang, B., Chen, Y., Stanton, C., Ross, R. P., Lee, Y.-K., Zhao, J., et al. (2019). Bifidobacterium and Lactobacillus Composition at species level and gut microbiota diversity in infants before 6 weeks. Int. J. Mol. Sci. 20:3306. doi: 10.3390/ijms20133306
Yang, R., Gao, R., Cui, S., Zhong, H., Zhang, X., Chen, Y., et al. (2019). Dynamic signatures of gut microbiota and influences of delivery and feeding modes during the first 6 months of life. Physiol. Genomics 51, 368–378. doi: 10.1152/physiolgenomics.00026.2019
Zheng, P., Li, Y., Wu, J., Zhang, H., Huang, Y., Tan, X., et al. (2019). Perturbed microbial ecology in myasthenia gravis: evidence from the gut microbiome and fecal metabolome. Adv. Sci. (Weinh.) 6:1901441. doi: 10.1002/advs.201901441
Zheng, P., Yang, J., Li, Y., Wu, J., Liang, W., Yin, B., et al. (2020). Gut microbial signatures can discriminate unipolar from bipolar depression. Adv. Sci. (Weinh.) 7:1902862. doi: 10.1002/advs.201902862
Zhuang, L., Chen, H., Zhang, S., Zhuang, J., Li, Q., and Feng, Z. (2019). Intestinal microbiota in early life and its implications on childhood health. Genomics Proteomics Bioinformatics 17, 13–25. doi: 10.1016/j.gpb.2018.10.002
Zuo, T., Kamm, M. A., Colombel, J.-F., and Ng, S. C. (2018). Urbanization and the gut microbiota in health and inflammatory bowel disease. Nat. Rev. Gastroenterol. Hepatol. 15, 440–452. doi: 10.1038/s41575-018-0003-z
Keywords: gut microbiota, cesarean section, Firmicutes/Bacteroidetes ratio, development, obesity
Citation: Long G, Hu Y, Tao E, Chen B, Shu X, Zheng W and Jiang M (2021) The Influence of Cesarean Section on the Composition and Development of Gut Microbiota During the First 3 Months of Life. Front. Microbiol. 12:691312. doi: 10.3389/fmicb.2021.691312
Received: 07 April 2021; Accepted: 26 July 2021;
Published: 18 August 2021.
Edited by:Junjun Wang, China Agricultural University, China
Reviewed by:María Victoria Quiroga, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Argentina
Casper Kamutando, University of Zimbabwe, Zimbabwe
Copyright © 2021 Long, Hu, Tao, Chen, Shu, Zheng and Jiang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mizu Jiang, email@example.com