Assessing the impact of pregnancy and birth factors on the maternal and infant microbiota
1Canadian Centre for Human Microbiome and Probiotic Research, Lawson Health Research Institute, London, Ontario N6A4V2, Canada.
2Department of Microbiology and Immunology, Western University, London, Ontario N6A3K7, Canada.
3London Health Sciences Centre, London, Ontario N6A5W9, Canada.
4Department of Obstetrics and Gynaecology, Division of Maternal Fetal Medicine, Western University, London, Ontario N6H5W9, Canada.
5Division of Urology, Department of Surgery, Western University, London, Ontario N6A4V2, Canada.
6Children’s Health Research Institute, London, Ontario N6C 4V3, Canada.
Correspondence to: Dr. Barbra de Vrijer, Department of Obstetrics and Gynaecology, Division of Maternal Fetal Medicine, Western University, Victoria Hospital, Commissioners Rd E 800, London, Ontario N6A5W9, Canada. Email:
Background: The microbiota acquired at birth is known to play an intimate role in later life health and disease and has been shown to be affected by the mode of birth. There has been recent interest in microbiota correction by maternal vaginal seeding in Cesarean section-born infants; however, the safety of this practice has been debated. The aim of this study was to assess how other factors, such as timing of sampling, maternal obesity, vaginal Group B Streptococcus colonization (GBS), and antibiotic exposure, affect the maternal and infant microbiota.
Methods: Maternal vaginal and saliva samples were collected at three time periods: 35-37 weeks gestation (prenatal), within 24-36 hours after birth (birth), and at ~6 weeks postpartum. Infant saliva and stool samples were collected at ~6 weeks postpartum. 16S rRNA amplicon sequencing was utilized to assess the taxonomic and inferred functional compositions of the bacterial communities from both mothers and infants.
Results: Samples from 36 mothers and 32 infants were obtained. Gestational age, breastfeeding, mode of birth, and gravidity were associated with taxonomic alterations in the infant samples, while obesity, antibiotic use, and GBS status were not. Maternal samples were predominantly affected by time, whereby significant alterations including increased microbial diversity were seen at birth and persisted to 6 weeks postpartum.
Conclusion: This study provides information on the relationship between health and delivery factors and changes in vaginal and infant microbiota. These results may better direct clinicians and mothers in optimizing the infant microbiota towards health during infancy and later life.
The human body harbors abundant microbiota, which varies by anatomical site and influences health and disease, including metabolic, immune, and infectious conditions. The composition of a newborn’s microbiota is altered by a variety of host genetic and environmental factors, and primary colonization occurs at birth by vertical transmission from the mother. During pregnancy, the vaginal microbiota is important for the health of both mother and fetus, and changes in vaginal microbiota are associated with pregnancy complications such as preterm birth and late miscarriage.
The mode of birth has been shown to be a major determinant of infant microbiota composition with potentially extensive later-life health implications. A vaginal birth exposes the baby to maternal vaginal and fecal bacteria, while Cesarean section-born infants are enriched in skin-associated bacteria[3,4]. After birth, Cesarean section microbial acquisition in the infant has been associated with an increased prevalence of chronic immune conditions. There has been recent interest in the practice of exposing infants to maternal vaginal fluids immediately following a Cesarean birth, commonly known as vaginal seeding. However, this practice does not consider changes in maternal microbiota, potentially putting infants at an increased risk of infection, and the long-term health consequences of vaginal seeding have yet to be established.
In addition to the mode of birth, factors such as maternal body mass index (BMI), antibiotic use during pregnancy, and early-life antibiotic use in the infant are also believed to affect normal microbiota acquisition[7-9]. Obesity rates have been steadily increasing and affect women of reproductive age, with approximately 45% beginning their pregnancy either overweight or obese. Obesity, regardless of and in the case of pregnancy, affects the composition of the microbiota[7,10] and is associated with higher rates of antibiotic use. Obesity during pregnancy is also associated with longer labors, higher rates of Cesarean section, chorioamnionitis, and Group B Streptococcus (GBS) colonization[15,16], which in North America is managed with routine antibiotic administration during labor, all factors that may independently affect the acquisition of a normal microbiome. Furthermore, breastfeeding success rates are lower in women who have delivered by Cesarean section and in women with obesity regardless of mode of delivery[18,19]. Importantly, with few exceptions, most pregnancy-related microbiota studies have focused on normal weight mothers or do not follow the infant microbiota, and thus the impact of other pregnancy factors remains unclear[20-22].
Through an assessment of taxonomic and inferred functional alterations in the vaginal and salivary microbiota from mothers, as well as salivary and gut microbiota from infants, this study sought to highlight the correlations between pregnancy and delivery factors, such as obesity, GBS status, antibiotic use, and breastfeeding, as they relate to changes in maternal and infant microbiota. This information may better direct clinicians and mothers towards optimizing the infant microbiota for later life health.
Study design and sample collection
This study was approved by the Health Sciences Research Ethics Board (HSREB# 105630) and written informed consent was received from all participants. Individuals with singleton pregnancies were recruited from London Health Sciences Centre, Victoria Hospital in London, Ontario, Canada. Individuals who received antibiotics in labor for presumed chorioamnionitis or prior to labor for conditions affecting the immune system (e.g., auto-immune disease, severe anemia, cardiac conditions, and severe preeclampsia) were under the age of 18 years or over the age of 40 years, or had a fetus with suspected congenital anomalies, severe growth restriction or chromosomal disorders were excluded from the study. Data such as maternal age, height, pre-pregnancy weight, and gestational weight gain, as well as infant birth weight and type of feeding were collected.
Demographic information, GBS screen results, and pregnancy and birth data were obtained from patient charts. An Upset plot was generated in R (version 4.2.1) with the UpSetR package (version 1.4.0) to display the intersection of study participants based on eight birth factors. Additionally, maternal vaginal and saliva swab samples were collected during three perinatal periods: around 35-37 weeks gestation (prenatal), within 24-36 hours after birth (birth), and at ~6 weeks postpartum (postpartum). Infant saliva and stool swabs were collected at the postpartum appointment approximately 6-8 weeks after birth. All swabs were stored at -80 °C until downstream processing.
Bacterial DNA isolation, amplification, and sequencing
Swabs were cut directly into the wells of the bead plate from a 96-well DNeasy PowerSoil HTP 96 kit (QIAGEN, Toronto, ON) and the protocol was followed as per the manufacturer’s directions, as previously described. A Biomek® 3000 Laboratory Automation Workstation was utilized for automated PCR reagent set-up. Amplifications of the V4 region of the 16S ribosomal RNA gene were carried out with the primers (5'-3') ACACTCTTTCCCTACACGACGCTCTTCCGATCTNNNNxxxxxxxxxxxxGTGCCAGCMGCCGCGGTAA and (5'-3') CGGTCTCGGCATTCCTGCTGAACCGCTCTTCCGATCTNNNNxxxxxxxxxxxxGGACTACHVGGGTWTCT wherein "xxxxxxxxxxxx" represents a sample-specific nucleotide barcode and the preceding sequence is a portion of the Illumina adapter sequence for library construction. The Biomek® robot transferred 2 μL of the DNA template into a plate containing 10 μL of forward and reverse primers (3.2 pmol/µL). Then 20 μL of Promega GoTaq® Colourless Master Mix (Promega, Maddison, WI) was added to the DNA template and primers. The final plate was firmly sealed with a foil PCR plate cover. This plate was placed in the thermal cycler (Techne, Mississauga, ON) where the lid was kept at 105 °C. An initial warm-up temperature of 95 °C was used for 2 min to activate the GoTaq®. Afterwards, the volumes underwent 25 cycles of 95 °C for 1 min, 52 °C for 1 min, and 72 °C for 1 min. After completion, the temperature of the thermal cycler was held at 4 °C and amplicons were then stored at -20 °C.
DNA sequencing was conducted at the London Regional Genomics Centre at Robarts Research Institute (London, ON). Amplicons were quantified using PicoGreen (Quant-It; Life Technologies, Burlington, ON) and pooled at equimolar concentrations before cleanup (QIAquick PCR clean up; QIAGEN inc). The final samples were sequenced using the MiSeq by Illumina® platform, with 2 × 260 bp paired-end chemistry with a cluster density of 1100 and 5% PhiX.
Sequencing dataset processing
After sequencing, the paired reads were demultiplexed using custom scripts, then quality filtered, trimmed, denoised, and merged following the DADA2 pipeline (version 1.26.0) in R. The SILVA Database (version 138) was utilized in assigning taxonomy to the amplicon sequence variants (SVs). Further filtering was performed such that SVs were removed that did not comprise > 1% of the relative abundance in any sample.
Inferencing of microbial metagenomes was performed with PICRUSt2 (version 2.4.1). Downstream analysis was performed on both gene families (Enzyme classification numbers, EC) and functional pathways (BioCyc ID). Gene family and pathway counts were not filtered but were rounded to integers for downstream analyses.
Alpha diversity (Shannon’s Index) was calculated using custom R scripts. A Kruskall-Wallis test followed by a Dunn’s multiple comparison test (from the rstatix R package, version 0.7.2) was used to test for alpha diversity differences between sample subgroups. Sparsity in read count data was zero-adjusted with the cmultRepl function from the zCompositions R package (version 1.4.0-1). The zero-adjusted counts underwent centered log ratio (CLR) transformation and a principal components analysis (PCA) was performed. PCA biplots were generated with ggplot2 (version 3.4.2). Sample type comparisons of Aitchison distance (β-diversity) were calculated by PERMANOVA with the adonis2 function from the vegan R package (version 2.6-4). Potential covariates of microbiota composition were determined using the envfit function from the vegan R package with the Aitchison distance PCA as the input object. Differential abundance tests of taxonomy, inferred gene families, and functional pathways were determined using a generalized linear mixed model approach with the MaAsLin2 R package (version 1.12.0). Briefly, CLR-transformed counts were used as the input, and the fixed effects were potential covariate clinical metadata features with an envfit P-value < 0.1. For maternal samples, the longitudinal repeated measures were accounted for by treating the ‘participant’ as a random effect. For all differential abundance analyses, Benjamini-Hochberg adjusted P-values (Q-values) < 0.1 were considered significant. Co-occurrence network analysis was performed using CoNet software (version 1.1.1), as previously described[32,33], where only significant (Benjamini-Hochberg corrected Browns P-value < 0.05) edges were plotted. CoNet networks were visualized in Cytoscape (version 3.8.2), then the final copresence and mutual exclusion network figures were generated with Circos (version 0.69-9). Relative network positivity was calculated as the number of copresence edges divided by the number of mutual exclusion edges. Network property and baby stool microbiota origin figures were generated with GraphPad Prism (version 9.5.1).
A total of 36 mothers and 32 infants were recruited for this study. Patient demographic characteristics are summarized in Table 1. Briefly, the cohort was heterogeneous with a range of maternal age, BMI, antibiotic exposure, delivery method, infant feeding, and other aspects. Maternal obesity was associated with higher rates of Cesarean birth, GBS, and antibiotic use during the perinatal or postpartum period, as well as lower breastfeeding rates. [Figure 1].
Figure 1. Adverse birth factors are concomitant. Upset plot depicting the intersecting study participants between eight birth factors. Each column corresponds to a factor or set of factors. The number of participants within each set (Intersection size) is displayed above the column. Set size represents the number of participants fitting the specified binary factor (row). Mothers with BMI ≥ 30 were considered obese. LGA: large for gestational age; GA: gestational age.
Clinical characteristics of study participants
|Patient Characteristics||Mean ( ± SD) or n (%)|
|Age at delivery (years)||31.2 ± 5.5 (range 21-46)|
|Pre-pregnancy BMI||33.4 ± 10.6 (range 21-49)|
|BMI ≥ 30||20 (56)|
|Gravidity||2.4 ± 1.4 (range 1-7)|
|Gestational age at birth (weeks)||38.9 ± 1.2 (range 35.6-40.7)|
|Mode of delivery|
|Cesarean section||9 (25)|
|GBS positive||9 (25)|
|During delivery||18 (50)|
|Infant birth weight (g)||3498 ± 476 (range 2300-4370)|
|Infant birth length (cm)||51.8 ± 2.9 (range 46-58)|
16S rRNA gene sequencing was performed to elucidate the potential impacts of birth and delivery factors on the maternal and infant microbiota. In total, 256 samples were sequenced from the 68 individuals (36 mothers and 32 infants). Following quality filtering with DADA2, the 16S dataset yielded 11,820,704 reads; the average sample read depth was 46,175 (range 2499 -202,456). 1341 unique SVs were initially identified, with 271 SVs remaining after filtering out those that did not comprise > 1% of the relative abundance. For the principal component analysis (PCA) of all samples, 30.5% of the total variance was explained by the first two principal components [Figure 2A]. The samples clustered based on type, with maternal saliva, maternal vaginal, infant saliva, and infant stool samples being significantly distinct (PERMANOVA r2 = 0.159, P < 0.001; Figure 2A and B). Shannon’s index of alpha diversity also significantly differed between sample groups (Kruskal-Wallis test, P < 0.0001). The diversity of baby saliva was significantly lower than the mothers’ (Dunn’s pairwise tests, adjusted P < 0.001). Vaginal diversity was lowest in the prenatal samples and was significantly increased following birth and in the postpartum samples (Dunn’s pairwise tests, adjusted P < 0.0001; Figure 2C).
Figure 2. Microbiota composition diverges by sample type and perinatal period. (A) PCA biplot based on sample-wise CLR-transformed Aitchison distance. Each point represents a sample. Strength and association for SVs are depicted by the length and direction of the gray arrows, respectively. Samples are colored by type and ellipses represent 95% confidence intervals. (B) Average relative genus abundance by sample type and perinatal period. Each vertical bar represents the average relative abundance within the sample type subgroup, colored by genus. Genera representing < 1% relative abundance were grouped into the “Remainder” portion. (C) Alpha diversity plot displaying the Shannon index for maternal and infant samples at different perinatal periods. Box plot whiskers represent minimum and maximum. PN: Prenatal; B: Birth; PP: Postpartum; CLR: centered log ratio; PCA: principal components analysis; SVs: sequence variants.
Potential covariates (envfit P < 0.1, Supplementary Table 1] were accounted for and sample type subgroups were explored for trends based on clinical features in a multivariate manner [Supplementary Tables 1 and 2]. Perinatal period was the most significant covariate associated with taxonomic and functional changes in the maternal vaginal samples (P < 0.001), though BMI and gravidity were also significant factors. Vaginal samples were dominated by lactobacilli prenatally [Figure 2B, Figure 3A, Supplementary Table 2], whereas samples from the birth and postpartum periods were significantly altered in 74 SVs and 36 genera when SV counts were aggregated by genus [Supplementary Table 2]. These postpartum alterations included relative enrichment of the genera Dialister, Prevotella, Escherichia, and Streptococcus. BMI was negatively correlated with Gemella and Megasphaera, while primigravida was associated with a relative depletion of Escherichia.
Figure 3. Taxonomic composition and functional co-occurrence network structure are altered in vaginal samples over time. (A) PCA biplot based on sample-wise CLR-transformed Aitchison distance. Each point represents a vaginal sample. Strength and association for genera (SVs) are depicted by the length and direction of the gray arrows, respectively. Samples are colored by perinatal period and ellipses represent 95% confidence intervals. Plots of networks indicate (B) copresence and (C) mutual exclusion patterns of enzymes, and (D) copresence and (E) mutual exclusion patterns of functional pathways in prenatal vaginal samples. Only significant edges are plotted. Edges are colored based on the enzyme or pathway group with more total edges (of the two involved in the interaction). Relative network positivity at both the enzymatic (F) and pathway (G) levels was calculated as the number of copresence edges divided by the number of mutual exclusion edges for the network from vaginal samples at each perinatal period. (H) Total number of nodes and significant edges are plotted for vaginal samples at each perinatal period.
Over 1000 gene families/enzymes and 275 functional pathways were significantly altered in vaginal samples based on perinatal period, gravidity, and BMI [Supplementary Tables 3 and 4]. Co-occurrence network analysis was performed on the vaginal samples’ functional potential at both the gene family and pathway levels [Figure 3B-E]. At both levels, the network structure was significantly altered and became less robust over time, with the relative edge positivity increasing (copresence / mutual exclusion relationships) and the number of significant edges decreasing [Figure 3F-H].
Delivery method and infant sex were determined to be potential covariates of baby stool microbiota. After accounting for these factors, the genus Citrobacter was relatively increased in infants delivered by Cesarean section (Q-value = 0.04). In baby saliva, the covariates of breastfeeding, gestational age and primigravida were adjusted; infants born prior to 37 weeks harbored relatively more SV41: Limosilactobacillus fermentum than those born after 37 weeks (Q-value = 1 × 10-7), while infants born out of primigravid pregnancies harbored relatively less of the genus Escherichia (Q-value = 0.048), SV3: Escherichia (Q-value = 0.05), and SV14: Bacteroides (Q-value = 0.05) compared to infants of multigravidae.
When investigating microbiota function, baby stool was significantly shaped by the delivery method, whereby infants delivered vaginally were relatively depleted in 41 different gene families compared to those delivered by Cesarean section [Supplementary Table 3]. Several of these enzymes are involved in lipopolysaccharide (EC.220.127.116.11) or E. coli O-polysaccharide biosynthesis (EC.18.104.22.1684, EC.22.214.171.124). Infants delivered vaginally were also relatively depleted in 13 functional pathways compared to those delivered by Cesarean section [Supplementary Table 4]. Several of these functional pathways were Pseudomonadota-specific, including 3-phenylpropanoate degradation (PWY0.1277) and sulfoquinovose degradation (PWY.7446).
Baby salivary microbiota function was significantly shaped by gravidity, gestational age, and breastfeeding [Supplementary Tables 3 and 4]. For example, breastfed infants had a relatively lower abundance of EC.126.96.36.199, an NADH oxidase enzyme likely originating from Streptococcus mutans (Q = 0.009), and infants from primigravidae had a relatively lower abundance of the functional pathway for biosynthesis of the virulence factor enterobactin (Q = 0.005).
To track the potential source of baby stool taxa, dyad analysis was performed by determining shared SVs between i) baby stool and baby oral microbiota, and ii) baby stool and vaginal microbiota. This analysis was performed on the filtered SV counts table, so all non-zero SV counts in dyad samples were considered shared. On average, baby stool samples contained 54 SVs, of which 54%, 58%, 53%, and 60% were shared with the baby’s saliva, prenatal vaginal, birth vaginal, and postpartum vaginal samples, respectively (Tukey’s multiple comparison test, P = 0.007; Figure 4A). The shared taxa were generally of high abundance (for example, Lactobacillus iners, Lactobacillus crispatus, Lactobacillus jensenii, and Streptococcus spp.), whereas more rare taxa were not consistently shared.
Figure 4. Towards prediction of the origin of baby stool microbiota by dyad analysis. (A) The percentage of baby stool microbiota SVs shared between the same baby’s saliva samples, or their maternal vaginal samples. Box plot whiskers represent the minimum and maximum. (B) PCA biplot based on sample-wise CLR-transformed Aitchison distance. Each point represents a baby stool sample. Strength and association for genera (SVs) are depicted by the length and direction of the gray arrows, respectively. Samples are colored by delivery method and shaped by sex; ellipses represent the 95% confidence intervals of the delivery method. (C) Beta diversity (Aitchison distance) between baby stool samples and prenatal vaginal samples compared based on clinical factors. CLR: centered log ratio; PCA: principal components analysis; SVs: sequence variants.
Beta diversity (Aitchison distance) between the baby fecal and prenatal vaginal sample pairs was investigated as a marker of vertical transmission and microbiota engraftment. In these comparisons, the shorter the Aitchison distance, the more similar the baby fecal and prenatal vaginal sample microbiota, suggesting increased “engraftment” from the vaginal microbiota into the baby’s gut. This “engraftment” was compared between different birth factors to explore which, if any, influenced the diversity between the sample dyads. Vaginal delivery led to significantly more engraftment (Bonferroni corrected T-test, P = 0.008), as the Aitchison distance between the baby gut microbiota and prenatal vaginal samples was lower than the distance between dyads with Cesarean delivery [Figure 4B and C]. Trends were also observed with more engraftment in male babies and those born to non-obese mothers, although these did not reach statistical significance.
This study investigated how the maternal and infant microbiota are shaped by pregnancy and delivery factors including obesity, timing of sampling, GBS status, antibiotic use, and breastfeeding. It was demonstrated that the vaginal microbiota was significantly altered at the time of birth, regardless of the mode of birth. Other birth factors were less influential but still associated with significant alterations to the various microbial communities.
Several studies have demonstrated how these pregnancy and delivery factors individually affect the infant microbiome[3,4,7-9], but few have considered a more comprehensive approach. Some studies have accounted for multiple factors; for example, one study found that, after adjusting for maternal race, prenatal antifungal use and intrapartum antibiotic use, maternal GBS status was significantly associated with infant gut bacterial composition at 6 months. In another study, maternal obesity was only associated with altered infant gut microbiome composition in vaginally-delivered infants, whereas this association was not seen in Cesarean-delivered infants. In this cohort, many of the pregnant individuals with an elevated BMI were also on antibiotics, were GBS-positive, and were not breastfeeding.
Additionally, many of these factors are known to interact with one another, and to make a greater impact on infant health, it may be best to first address overall maternal health before and during pregnancy. For example, although some studies have found an association between antibiotic use in pregnancy and an increased risk of childhood obesity, one should also consider the reasons for prenatal antibiotics in the first place. For one, obesity is associated with an increased risk of bacterial infections during pregnancy and, thus, higher rates of antibiotic use. Furthermore, the increased risk of childhood obesity may be the result of maternal infection, rather than antibiotic use during pregnancy, as noted by one study.
In agreement with previous literature, this study determined that the mode of birth was associated with significant alterations to the baby stool microbiota, including enrichment of aerobic microbes in infants born by Cesarean section. This is unsurprising given the association of vaginal birth with anaerobic vaginal and gut microbes, rather than the Cesarean section seeding of microbes from the skin and external environment[3,4,42,43]. The current results contrast with previous findings where antibiotic exposure[44-46] and breastfeeding[47,48] have demonstrated significant association with baby stool microbiota; the lack of detectable signals in this cohort is likely due to the small and phenotypically heterogeneous population.
Numerous alterations to the infant salivary microbiota were associated with gestational age, breastfeeding, and primigravity. While first pregnancy is known to influence labor duration and a recent study suggests its association with altered gut microbiota and increased BMI in children between 1-3 years of age, we believe this to be the first instance detailing its influence on the infant’s oral microbiota. At the early timepoint measured in the current study, the oral cavity may be the first bacterial community to exhibit these distinct and novel microbial indicators, which, at a later age, could be consequential to health. Specifically, the changes determined here in oral lactobacilli may have future implications for childhood oral candidiasis and caries development[50,51]. The oral cavity also seeds the gut microbiota, and since many systemic diseases are associated with changes in the microbial community of both sites, gaining a better understanding of the factors shaping the oral microbiota could expand our ability to diagnose and predict both oral and systemic diseases.
Maternal obesity was significantly associated with variation in only the vaginal microbiota, indicating that it may not be as relevant to the infant microbiota as other birth factors including mode of delivery, primigravity, breastfeeding, and gestational age. Again, it is possible that our limited sample size and heterogeneous maternal population decreased power to detect significant microbiota-related factors, or that infant manifestations of the relevant factors were detectable beyond the postpartum timepoint measured in this study (6-8 weeks after birth). In this regard, the factors determined to be significant microbiota covariates (most of which were investigated in maternal and infant samples from multiple anatomical sites here for the first time) warrant further exploration in larger longitudinal studies to determine the role (if any) they have on microbiota development and overall health later in life.
The success of fecal microbiome transplant for microbially-involved and systemic diseases beyond Clostridioides difficile infection has led to the proposal of transplanting alternative materials and anatomical regions, including the vaginal microbiome for bacterial vaginosis and urinary tract infection[54-60]. Similarly, based on the known negative associations with Cesarean section microbial acquisition, interest in microbiota correction by maternal vaginal seeding has increased in recent years[43,61-63]. While this field is still evolving, it is important to characterize the optimally efficacious donor sample; studies diverge in this regard, with some utilizing maternal vaginal material, others suggesting maternal gut material, or even arguing that species-specific restoration is sufficient. While maternal stool was not investigated here, beyond any other birth factor, perinatal period was the most significant driver of vaginal microbiota community structure and function. The baby stool microbiota also shared the fewest microbial taxa with the vaginal sample collected at birth. Thus, if indeed vaginal seeding is to be undertaken, the sample should be collected prior to the onset of labor, after which the microbial community structure and overall function significantly change without recovery by 6 weeks postpartum. The extended displacement of vaginal lactobacilli during and after birth indicates that additional microbial therapeutic support may be advantageous to all perinatal individuals, and not just those with Cesarean births. Maternal probiotic supplementation during pregnancy and postpartum has shown promise in long-lasting health outcomes in both mother and child related to obesity[65,66], allergy, depression and anxiety, and necrotizing enterocolitis, among other conditions, warranting consideration for incorporation into standard perinatal care.
This study has several limitations. As previously mentioned, the sample size was modest, and the maternal cohort was heterogeneous in BMI, antibiotic exposure, GBS status, age, delivery, and infant feeding method. Although a multivariate approach was undertaken to account for such differences, a future study with a larger sample size could have better precision in detecting and adjusting for potential confounders. Maternal stool was also not investigated here. This study employed 16S rRNA amplicon sequencing of all samples; this technique provides relative compositional, but not absolute abundance information. It also did not provide taxonomic annotation to the species level in all cases, so caution should be taken when comparing these data and taxonomic annotations in future studies. The 16S rRNA amplicon data were used to infer functional metagenomes with PICRUSt2; although this technique demonstrates accuracy in samples of human origin[27,70], it is not a replacement for whole shotgun metagenomic sequencing which could be employed in the future to validate the current findings.
Overall, this study provides the first look at how a multitude of birth factors may influence maternal and infant microbiota. The maternal salivary and vaginal microbiota were associated most significantly with temporal factors related to birth, while the infant samples were influenced to a lesser degree by a range of factors including mode of delivery, breastfeeding, and gestational age. As obesity, GBS status, and antibiotic use continue to increase in our population, it is important to concurrently study these features and their interaction and effect on both the maternal and infant microbiota. Whether these differential features manifest in the infant’s later life disease risk requires additional study.
We thank Mrs. Laura McMurphy, who assisted with recruitment and sample collection. We thank David Carter at London Regional Genomics Centre for sequencing support. We acknowledge the Anishinaabek, Haudenosaunee, Lūnaapéewak and Chonnonton Nations, whose traditional territories are where this study was undertaken and the publication was produced.
Performed sample processing, data analysis and interpretation, figure generation, writing of original draft, and editing of the final draft: Al KF
Performed participant recruitment, sample and data acquisition, technical support, writing of original draft and editing final draft: Allen L
Contributed to data interpretation and editing final draft: Bedell S
Contributed to conception and study design, editing final draft: Burton JP
Acquired funding, contributed to conception and study design, sample and data acquisition, performed data interpretation, editing final draft: de Vrijer B
Availability of data and materials
All 16S rRNA gene sequencing data have been deposited into the NCBI Sequence Read Archive, BioProject ID # PRJNA801237.
Financial support and sponsorship
Funding was provided by the Western University Department of Obstetrics and Gynaecology Academic Enrichment Fund.
Conflicts of interest
All authors declare that there are no conflicts of interest.
Ethical approval and consent to participate
This study was approved by the Health Sciences Research Ethics Board (HSREB# 105630) of Western University, and written informed consent was received from all participants.
Consent for publication
© The Author(s) 2023.
1. Nuriel-Ohayon M, Neuman H, Koren O. Microbial changes during pregnancy, birth, and infancy. Front Microbiol 2016;7:1031.
2. Kindinger LM, Bennett PR, Lee YS, et al. The interaction between vaginal microbiota, cervical length, and vaginal progesterone treatment for preterm birth risk. Microbiome 2017;5:6.
3. Dominguez-Bello MG, Costello EK, Contreras M, et al. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci U S A 2010;107:11971-5.
4. Shao Y, Forster SC, Tsaliki E, et al. Stunted microbiota and opportunistic pathogen colonization in caesarean-section birth. Nature 2019;574:117-21.
5. Keag OE, Norman JE, Stock SJ. Long-term risks and benefits associated with cesarean delivery for mother, baby, and subsequent pregnancies: systematic review and meta-analysis. PLoS Med 2018;15:e1002494.
6. Wilson BC, Butler ÉM, Grigg CP, et al. Oral administration of maternal vaginal microbes at birth to restore gut microbiome development in infants born by caesarean section: A pilot randomised placebo-controlled trial. EBioMedicine 2021;69:103443.
7. Soderborg TK, Clark SE, Mulligan CE, et al. The gut microbiota in infants of obese mothers increases inflammation and susceptibility to NAFLD. Nat Commun 2018;9:4462.
8. Zou ZH, Liu D, Li HD, et al. Prenatal and postnatal antibiotic exposure influences the gut microbiota of preterm infants in neonatal intensive care units. Ann Clin Microbiol Antimicrob 2018;17:9.
9. Gunderson EP. Childbearing and obesity in women: weight before, during, and after pregnancy. Obstet Gynecol Clin North Am 2009;36:317-32, ix.
10. Castaner O, Goday A, Park YM, et al. The gut microbiome profile in obesity: a systematic review. Int J Endocrinol 2018;2018:4095789.
11. Merrill RM, Fowers R. To what extent does sex, age and BMI impact medical and pharmacy costs? A retrospective cohort study involving employees in a large school district in the USA. BMJ Open 2019;9:e024078.
12. Valente RP, Santos P, Ferraz T, Montenegro N, Rodrigues T. Valente RP, Santos P, Ferraz T, Montenegro N, Rodrigues T. Effect of obesity on labor duration among nulliparous women with epidural analgesia. J Matern Fetal Neonatal Med 2020;33:2195-201.
13. Marchi J, Berg M, Dencker A, Olander EK, Begley C. Risks associated with obesity in pregnancy, for the mother and baby: a systematic review of reviews. Obes Rev 2015;16:621-38.
14. Hadley EE, Discacciati A, Costantine MM, et al. Maternal obesity is associated with chorioamnionitis and earlier indicated preterm delivery among expectantly managed women with preterm premature rupture of membranes. J Matern Fetal Neonatal Med 2019;32:271-8.
15. Manzanares S, Zamorano M, Naveiro-Fuentes M, Pineda A, Rodríguez-Granger J, Puertas A. Maternal obesity and the risk of group B streptococcal colonisation in pregnant women. J Obstet Gynaecol 2019;39:628-32.
16. Kleweis SM, Cahill AG, Odibo AO, Tuuli MG. Maternal obesity and rectovaginal group b streptococcus colonization at term. Infect Dis Obstet Gynecol 2015;2015:586767.
17. Money D, Allen VM. The prevention of early-onset neonatal group b streptococcal disease. J Obstet Gynaecol Can 2016;38:S326-35.
18. Hobbs AJ, Mannion CA, McDonald SW, Brockway M, Tough SC. The impact of caesarean section on breastfeeding initiation, duration and difficulties in the first four months postpartum. BMC Pregnancy Childbirth 2016;16:90.
19. Babendure J, Reifsnider E, Mendias E, Moramarco MW, Davila YR. Reduced breastfeeding rates among obese mothers: a review of contributing factors, clinical considerations and future directions. Int Breastfeed J 2015;10:21.
20. Garcia-Mantrana I, Collado MC. Obesity and overweight: impact on maternal and milk microbiome and their role for infant health and nutrition. Mol Nutr Food Res 2016;60:1865-75.
21. Mueller NT, Shin H, Pizoni A, et al. Birth mode-dependent association between pre-pregnancy maternal weight status and the neonatal intestinal microbiome. Sci Rep 2016;6:23133.
22. Santacruz A, Collado MC, García-Valdés L, et al. Gut microbiota composition is associated with body weight, weight gain and biochemical parameters in pregnant women. Br J Nutr 2010;104:83-92.
23. Al KF, Bisanz JE, Gloor GB, Reid G, Burton JP. Evaluation of sampling and storage procedures on preserving the community structure of stool microbiota: a simple at-home toilet-paper collection method. J Microbiol Methods 2018;144:117-21.
24. Al KF, Burton JP. Processing human urine and ureteral stents for 16S rRNA amplicon sequencing. STAR Protoc 2021;2:100435.
25. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods 2016;13:581-3.
26. Quast C, Pruesse E, Yilmaz P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 2013;41:D590-6.
27. Douglas GM, Maffei VJ, Zaneveld JR, et al. PICRUSt2 for prediction of metagenome functions. Nat Biotechnol 2020;38:685-8.
28. Palarea-albaladejo J, Martín-fernández JA. zCompositions-R package for multivariate imputation of left-censored data under a compositional approach. Chemom Intell Lab Syst 2015;143:85-96.
30. Oksanen J, Simpson G, Blanchet, FG. vegan: an R package for community ecologists. Availble from: https://github.com/vegandevs/vegan. [Last accessed on 13 Jul].
31. Mallick H, Rahnavard A, McIver LJ, et al. Multivariable association discovery in population-scale meta-omics studies. PLoS Comput Biol 2021;17:e1009442.
32. Friedman J, Alm EJ. Inferring correlation networks from genomic survey data. PLoS Comput Biol 2012;8:e1002687.
33. Daisley BA, Chanyi RM, Abdur-Rashid K, et al. Abiraterone acetate preferentially enriches for the gut commensal Akkermansia muciniphila in castrate-resistant prostate cancer patients. Nat Commun 2020;11:4822.
34. Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003;13:2498-504.
35. Krzywinski M, Schein J, Birol I, et al. Circos: an information aesthetic for comparative genomics. Genome Res 2009;19:1639-45.
36. Saha P, Xiao X, Yeoh BS, et al. The bacterial siderophore enterobactin confers survival advantage to salmonella in macrophages. Gut Microbes 2019;10:412-23.
37. Cassidy-Bushrow AE, Sitarik A, Levin AM, et al. Maternal group B streptococcus and the infant gut microbiota. J Dev Orig Health Dis 2016;7:45-53.
38. Singh SB, Madan J, Coker M, et al. Does birth mode modify associations of maternal pre-pregnancy BMI and gestational weight gain with the infant gut microbiome? Int J Obes (Lond) 2020;44:23-32.
39. Mueller NT, Whyatt R, Hoepner L, et al. Prenatal exposure to antibiotics, cesarean section and risk of childhood obesity. Int J Obes 2015;39:665-70.
40. Ovalle A, Martínez MA, Fuentes A, et al. Obesity, a risk factor for ascending bacterial infection during pregnancy. Rev Med Chil 2016;144:476-82.
41. Li DK, Chen H, Ferber J, Odouli R. Maternal infection and antibiotic use in pregnancy and the risk of childhood obesity in offspring: a birth cohort study. Int J Obes 2020;44:771-80.
42. Kim G, Bae J, Kim MJ, et al. Delayed establishment of gut microbiota in infants delivered by cesarean section. Front Microbiol 2020;11:2099.
43. Dominguez-Bello MG, De Jesus-Laboy KM, Shen N, et al. Partial restoration of the microbiota of cesarean-born infants via vaginal microbial transfer. Nat Med 2016;22:250-3.
44. Azad MB, Konya T, Persaud RR, et al. CHILD study investigators. Impact of maternal intrapartum antibiotics, method of birth and breastfeeding on gut microbiota during the first year of life: a prospective cohort study. BJOG 2016;123:983-93.
45. Stearns JC, Simioni J, Gunn E, et al. Intrapartum antibiotics for GBS prophylaxis alter colonization patterns in the early infant gut microbiome of low risk infants. Sci Rep 2017;7:16527.
46. Zhang M, Differding MK, Benjamin-Neelon SE, Østbye T, Hoyo C, Mueller NT. Association of prenatal antibiotics with measures of infant adiposity and the gut microbiome. Ann Clin Microbiol Antimicrob 2019;18:18.
47. Azad MB, Konya T, Maughan H, et al. CHILD study investigators. Gut microbiota of healthy Canadian infants: profiles by mode of delivery and infant diet at 4 months. CMAJ 2013;185:385-94.
48. Fallani M, Young D, Scott J, et al. Other members of the INFABIO team. Intestinal microbiota of 6-week-old infants across Europe: geographic influence beyond delivery mode, breast-feeding, and antibiotics. J Pediatr Gastroenterol Nutr 2010;51:77-84.
49. Vu K, Lou W, Tun HM, et al. From birth to overweight and atopic disease: multiple and common pathways of the infant gut microbiome. Gastroenterology 2021;160:128-144.e10.
50. Xiao J, Fiscella KA, Gill SR. Oral microbiome: possible harbinger for children’s health. Int J Oral Sci 2020;12:12.
51. Twetman S, Stecksén-Blicks C. Probiotics and oral health effects in children. Int J Paediatr Dent 2008;18:3-10.
52. Kelly BJ, Gross R, Bittinger K, et al. Power and sample-size estimation for microbiome studies using pairwise distances and PERMANOVA. Bioinformatics 2015;31:2461-8.
53. Mitchell AM, Palettas M, Christian LM. Fetal sex is associated with maternal stimulated cytokine production, but not serum cytokine levels, in human pregnancy. Brain Behav Immun 2017;60:32-7.
54. Bakken JS, Borody T, Brandt LJ, et al. Fecal microbiota transplantation workgroup. Treating clostridium difficile infection with fecal microbiota transplantation. Clin Gastroenterol Hepatol 2011;9:1044-9.
55. Craven L, Rahman A, Nair Parvathy S, et al. Allogenic fecal microbiota transplantation in patients with nonalcoholic fatty liver disease improves abnormal small intestinal permeability: a randomized control trial. Am J Gastroenterol 2020;115:1055-65.
56. Baruch EN, Youngster I, Ben-Betzalel G, et al. Fecal microbiota transplant promotes response in immunotherapy-refractory melanoma patients. Science 2021;371:602-9.
57. de Groot P, Nikolic T, Pellegrini S, et al. Faecal microbiota transplantation halts progression of human new-onset type 1 diabetes in a randomised controlled trial. Gut 2021;70:92-105.
58. Lev-Sagie A, Goldman-Wohl D, Cohen Y, et al. Vaginal microbiome transplantation in women with intractable bacterial vaginosis. Nat Med 2019;25:1500-4.
59. DeLong K, Bensouda S, Zulfiqar F, et al. Conceptual design of a universal donor screening approach for vaginal microbiota transplant. Front Cell Infect Microbiol 2019;9:306.
60. Al KF, Craven LJ, Gibbons S, et al. Fecal microbiota transplantation is safe and tolerable in patients with multiple sclerosis: a pilot randomized controlled trial. Mult Scler J Exp Transl Clin 2022;8:20552173221086662.
61. Korpela K, Costea P, Coelho LP, et al. Selective maternal seeding and environment shape the human gut microbiome. Genome Res 2018;28:561-8.
62. Korpela K, Helve O, Kolho KL, et al. Maternal fecal microbiota transplantation in cesarean-born infants rapidly restores normal gut microbial development: a proof-of-concept study. Cell 2020;183:324-334.e5.
63. Duar RM, Kyle D, Tribe RM. Reintroducing B. infantis to the cesarean-born neonate: an ecologically sound alternative to “vaginal seeding”. FEMS Microbiol Lett 2020;367:fnaa032.
64. MacIntyre DA, Chandiramani M, Lee YS, et al. The vaginal microbiome during pregnancy and the postpartum period in a European population. Sci Rep 2015;5:8988.
65. Vähämiko S, Laiho A, Lund R, Isolauri E, Salminen S, Laitinen K. The impact of probiotic supplementation during pregnancy on DNA methylation of obesity-related genes in mothers and their children. Eur J Nutr 2019;58:367-77.
66. Ilmonen J, Isolauri E, Poussa T, Laitinen K. Impact of dietary counselling and probiotic intervention on maternal anthropometric measurements during and after pregnancy: a randomized placebo-controlled trial. Clin Nutr 2011;30:156-64.
67. Sulaiman T, Medi S, Erdem H, et al. The diagnostic utility of the “thwaites’ system” and “lancet consensus scoring system” in tuberculous vs. non-tuberculous subacute and chronic meningitis: multicenter analysis of 395 adult patients. BMC Infect Dis 2020;20:788.
68. Slykerman RF, Hood F, Wickens K, et al. Probiotic in pregnancy study group. effect of lactobacillus rhamnosus hn001 in pregnancy on postpartum symptoms of depression and anxiety: a randomised double-blind placebo-controlled trial. EBioMedicine 2017;24:159-65.
69. Benor S, Marom R, Ben Tov A, Armoni Domany K, Zaidenberg-Israeli G, Dollberg S. Probiotic supplementation in mothers of very low birth weight infants. Am J Perinatol 2014;31:497-504.
70. Sun S, Jones RB, Fodor AA. Inference-based accuracy of metagenome prediction tools varies across sample types and functional categories. Microbiome 2020;8:46.
Cite This Article
Al KF, Allen L, Bedell S, Burton JP, de Vrijer B. Assessing the impact of pregnancy and birth factors on the maternal and infant microbiota. Microbiome Res Rep 2023;2:29. http://dx.doi.org/10.20517/mrr.2023.32
Al KF, Allen L, Bedell S, Burton JP, de Vrijer B. Assessing the impact of pregnancy and birth factors on the maternal and infant microbiota. Microbiome Research Reports. 2023; 2(4): 29. http://dx.doi.org/10.20517/mrr.2023.32
Al, Kait F, Laura Allen, Samantha Bedell, Jeremy P Burton, Barbra de Vrijer. 2023. "Assessing the impact of pregnancy and birth factors on the maternal and infant microbiota" Microbiome Research Reports. 2, no.4: 29. http://dx.doi.org/10.20517/mrr.2023.32
Al, KF.; Allen L.; Bedell S.; Burton JP.; de Vrijer B. Assessing the impact of pregnancy and birth factors on the maternal and infant microbiota. Microbiome. Res. Rep. 2023, 2, 29. http://dx.doi.org/10.20517/mrr.2023.32
Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at email@example.com.