Abstract
The microbiome is a regulator of host immunity, metabolism, neurodevelopment, and behavior. During early life, bacterial communities within maternal gut and vaginal compartments can have an impact on directing these processes. Maternal stress experience during pregnancy may impact offspring development by altering the temporal and spatial dynamics of the maternal microbiome during pregnancy. To examine the hypothesis that maternal stress disrupts gut and vaginal microbial dynamics during critical prenatal and postnatal windows, we used high-resolution 16S rRNA marker gene sequencing to examine outcomes in our mouse model of early prenatal stress. Consistent with predictions, maternal fecal communities shift across pregnancy, a process that is disrupted by stress. Vaginal bacterial community structure and composition exhibit lasting disruption following stress exposure. Comparison of maternal and offspring microbiota revealed that similarities in bacterial community composition was predicted by a complex interaction between maternal body niche and offspring age and sex. Importantly, early prenatal stress influenced offspring bacterial community assembly in a temporal and sex-specific manner. Taken together, our results demonstrate that early prenatal stress may influence offspring development through converging modifications to gut microbial composition during pregnancy and transmission of dysbiotic vaginal microbiome at birth.
Similar content being viewed by others
Introduction
Early environmental perturbations, such as chronic stress or infection during pregnancy, are associated with profound disruption to the developing central nervous system, including circuits that regulate stress, affect, and cognition1,2,3,4,5. Rodent and primate models of early life adversity demonstrate alterations to the gut-brain axis that mimic endophenotypes of neurodevelopmental disorders, including increased intestinal permeability, visceral hypersensitivity, stress pathway dysregulation, altered development of the enteric nervous system, and changes to behavior6,7,8,9,10,11,12,13,14,15,16,17,18,19. The mechanisms through which fetal antecedents contribute to gut-brain dysfunction likely involve dynamic interactions between the maternal and fetal compartments1. One potential interaction that has yet to be explored is the prolonged effect of stress on bacterial communities during pregnancy.
During pregnancy a continuous maternal supply of nutrients and substrates to the fetus is essential to support normal growth and development20. This nutritional and metabolic demand of the developing fetus is likely modulated by the capacity of the maternal gut microbiome to harvest and produce these substrates21,22. Indeed, metabolites of maternal gut microbial origin, such as short chain fatty acids, saccharides, and amino acids cross the placental barrier, enter fetal circulation, and contribute to the formation of the blood-brain barrier and innate immune development23,24. Previous work in our laboratory has shown that maternal stress experience during the first week of pregnancy results in alterations to placental transfer of nutrients in a sex-specific manner25,26,27. As maternal nutritional status influences placental function, alterations to the structure and metabolic potential of the maternal bacterial microbiome may represent a currently uncharacterized association through which maternal stress experience impacts offspring outcomes28,29.
During parturition the bacterial communities within the maternal vagina also play a significant role in offspring postnatal development6,7,8,9,13,14,15,18,30,31. As maternal vaginal bacterial communities provide the primary inoculum that colonize the neonate gut at birth, the composition of these pioneer communities contributes disproportionally to long-term health outcomes32,33,34,35,36. Maternal factors, such as maternal immune activation, diet, antibiotics, and mode of delivery impact vertical transmission of vaginal communities33,34,36,37,38,39,40. As a result, transmission of dysbiotic vaginal bacterial communities has been proposed to disrupt neonate-microbe interactions necessary for immune education, metabolism, and neurodevelopment37,41,42,43,44. Consistent with this hypothesis, previous work in our laboratory has shown that early prenatal stress decreased the relative abundance of the commensal Lactobacillus in the maternal vagina and neonate gut at colonization. Further, these paralleled changes to offspring gut and plasma metabolite profiles related to impaired energy metabolism and availability and sex- and brain region-specific changes in amino acid transport18. Transplantation of maternal vaginal microbiota to cesarean delivered mice was sufficient in recapitulating bacterial load observed in vaginally delivered offspring, recapitulating recent evidence in infants18,34. While these results show that maternal stress experience impacts microbial community composition and elicits metabolic consequences at birth, the lasting consequences of transmission of dysbiotic vaginal communities on the assembly of offspring microbial communities remains unclear.
Therefore, to examine the hypothesis that maternal stress experience disrupts gut and vaginal microbial dynamics during critical prenatal and postnatal windows, we used a high-resolution sampling strategy coupled with 16S rRNA marker gene sequencing to examine outcomes in our mouse model of early prenatal stress, in which male, but not female, offspring demonstrate significant neurodevelopmental changes in hypothalamic and limbic circuits and in the regulation of stress responsivity, cognitive dysfunction, and post-pubertal growth18,25,26,45,46,47,48. This approach provides an opportunity to draw novel associations between the effects of stress during pregnancy on the structure and functional potential of maternal bacterial communities and sex-specific phenotypes that emerge in offspring as adults. Further, unique assemblages of microbiota that differentiate treatment groups as a function of time, sex, and their interaction were identified using a supervised machine learning technique. Together, these studies provide an in-depth timeline reconstruction of temporal and spatial kinetics of maternal microbial communities following stress exposure and subsequent influence on the assembly of offspring microbial communities.
Results
Stress alters maternal fecal microbial diversity, community, and composition during pregnancy
Mammalian pregnancy exhibits two distinct metabolic states: an anabolic stage characterized by a linear phase in maternal and fetal growth and a catabolic stage characterized by ketogenesis, gluconeogenesis, and accelerated placental transfer of nutrients, all of which support the exponential growth rate of the developing fetus20. In mice the linear growth phase is restricted to the first week of pregnancy, followed by an exponential growth phase the remaining two weeks of pregnancy (r2 = 0.94, p < 0.01, Fig. 1a). Previous work has shown that the maternal gut microbiome undergoes dynamic remodeling during pregnancy, favoring taxa with an increased capacity for extracting nutrients necessary to meet the nutritional and metabolic demand of pregnancy and in preparation for lactation21. To determine the impact of stress during the first week of pregnancy on murine gut microbial communities, fecal samples were obtained starting one day prior to mating and then daily during the course of pregnancy (nineteen days in mice) (Supplementary Fig. 1A). To determine temporal patterns in the bacterial community structure during pregnancy, alpha and beta diversity measures were calculated and compared as a function of day of pregnancy. Linear regression modeling revealed that community diversity across pregnancy was best predicted by a sigmoid function, characterized by an initial decrease during the first eight days and followed by an increase and stabilization during the remaining eleven days of pregnancy (Fig. 1b). Further, maternal body weight was significantly correlated Shannon diversity (r2 = 0.11, p = 0.011, Supplementary Fig. 1B). To determine the temporal patterning of bacterial taxa across pregnancy, unsupervised hierarchical clustering of 97%-identity OTUs was conducted. Consistent with Shannon diversity results, taxa clustered into two distinct stages (Fig. 1c, Supplemental Table 1).
(a) Percent increase in total body weight across pregnancy demonstrates a slow and exponential growth phase during pregnancy in mice. There was no effect of stress on total body weight across pregnancy. N = 4–8 females per time point. (b) Shannon diversity index plotted against day of pregnancy, revealing that dynamic effects of pregnancy on community diversity. N = 4–8 females sampled daily per treatment, yielding 228 unique fecal samples in total; Non-parametric t-Test *p < 0.05, **p < 0.01, ***p < 0.001. (c) Heatmap depicting unsupervised clustering of samples (columns) of 97%-identity OTUs that exhibit at least 1% abundance in at least 2 females, demonstrating that taxa form two distinct clusters. The taxonomic annotation of each OTU is indicated to the left of rows. (d) Distinct structure of fecal bacterial communities during murine pregnancy. Communities clustered using PCoA of the unweighted UniFrac distance matrix, and each point corresponds to a sample that was collected daily during pregnancy. The percentage of variation explained by the PC is indicated on the axes, and colors correspond to samples collected either during early (white) and late (grey) pregnancy. N = 4–8 females sampled daily per treatment; PERMANOVA p-values based on 999 Monte Carlo simulations; PC, principal coordinate; PCoA, principal coordinates analysis.(e) Stress alterations to microbial community structure are limited to period of stress exposure. Average (±sem) unweighted UniFrac distances comparing phylogenetic configurations of fecal microbial communities of nonpregnant females, and control and stress-exposed females during early and late pregnancy. Mann Whitney U, ***p < 0.001. (f) Correlation matrix analysis of Random Forests-identified taxa that discriminatory across stage of pregnancy, demonstrating that stress disrupts co-occurrence and co-exclusion relationships among taxa. Empty boxes indicate nonsignificant correlations, while filled boxes are indicative of a significant correlation with a P-value ≤ 0.05 by Kendall’s test. The shading intensity of the bubble, along with size, is indicative of the Kendall rank correlation coefficient between matrices. Blue designates a negative correlation while red designates a positive correlation. N = 4–8 females sampled daily per treatment, yielding 228 unique fecal samples in total.
As the Shannon Diversity Index and unsupervised clustering of 97%-identity OTUs indicated two distinct phases of pregnancy-associated community structure, the longitudinal data was binned into early (e.g., day 1–8) and late (e.g., day 9–19) stages of pregnancy and the bacterial community structures were assessed. Unweighted UniFrac distance analysis confirmed differences in community structure between these two periods (PERMANOVA, pseudo-F = 35.55, p = 0.001) (Fig. 1d). Inter-individual variation was significantly greater during early pregnancy than late pregnancy in control females, while stress decreased UniFrac distances among exposed females during both stages of pregnancy (Fig. 1e). Similarly, the community structure between control and stress-exposed females was significantly different during the early and late period of pregnancy (Fig. 1e).
To identify the bacterial signature that accounts for the observed shift in community structure following stress exposure, we applied Random Forests models to our longitudinal data as previously described49. This approach identified 13 discriminatory taxa whose relative abundance was differentially represented between early and late pregnancy in control females (Supplemental Fig. 2A). Members of the order Clostridiales of the phylum Firmicutes (Lachnospiraceae and Ruminococcaceae) increase during late pregnancy, and is paralleled by decreases in taxa within the Bacteroidetes phyla (Prevotella and S24–7) (Supplemental Fig. 2B). Reconfiguration of taxa within the Bacteroidetes phylum was observed in stress-exposed females, including enrichment of Rikenellaceae and Odoribacter with a concomitant loss of Bacteroides that emerged early in pregnancy and remained stable until parturition. Enrichment of the Desulfovibrionaceae family within the Proteobacteria phylum and Mucispirillum within the Deferribacteres phylum was observed during late pregnancy. Stress increased Mucispirillum abundance during the early stage of pregnancy and remained elevated throughout pregnancy (Supplemental Fig. 2B). As taxa exhibit significant co-occurrence and co-exclusion relationships, we determined whether stress exposure during pregnancy impacted these interactions by generating correlation matrices of Random Forests-identified taxa50. Significant correlations were detected between taxa across pregnancy that were disrupted by stress. For instance, the significant positive correlation between Mucispirillum and Lachnospiraceae in control dams disappeared in dams exposed to stress during pregnancy, while new correlations emerged in stress-exposed females, including a significant positive correlation between Mucispirillum and Flexispira (Fig. 1g). Taken together, these results suggest that stress during the first week of pregnancy altered community dynamics of maternal fecal bacteria.
Given the well-established functional redundancy of bacterial communities, shifts in community composition may not always parallel functional alterations51. To determine whether stress influenced the functional potential of bacterial communities, 16S rRNA gene sequencing data were used to predict the abundance of KEGG functional orthologs using phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt), which were then filtered, and discriminatory features were identified with Random Forests regressions52. Random Forests models revealed 30 functional categories that discriminate between early and late pregnancy in control females, including pathways involved in gluconeogenesis, lipid biosynthesis, amino acids and fatty acid metabolism, and N-glycan degradation. (Supplemental Fig. 2A,B). Together, these results show that maternal stress experience during the first week of pregnancy exerts lasting disruption to community structure, composition, and functional capacity.
Chronic stress experience alters maternal vaginal bacteria composition
As the maternal vaginal bacteria provide the primary inoculum of bacterial communities that colonize the neonate gut at birth, we next assessed the impact of stress on these communities. To determine whether temporal patterns of the maternal vaginal bacteria are impacted by stress during pregnancy and the postpartum period, vaginal fluid samples were obtained a day prior to mating, at the end of the first week of pregnancy (i.e., end of stress exposure), and 48 h postpartum. A significant increase in the vaginal Shannon Diversity Index was observed on gestational day 7.5 (E7.5) in control females relative to nonpregnant and postpartum control females (ps < 0.05, t test), while markedly blunted the effect of early pregnancy on vaginal community structure (Fig. 2a). Similarly, a difference in UniFrac distances emerged between nonpregnant and the E7.5 vaginal communities of control females. The vaginal community structure was also significantly different between control and stress-exposed females (ps < 0.05, permutation test) (Fig. 2b).
(a) Shannon diversity index plotted against sampling time point demonstrating that stress prevents the increase in microbial diversity at E7.5 as observed in control females. N = 5–8/day/treatment; Mann Whitney U, *p < 0.05. (b) Stress alters microbial community structure in the vagina. Average (±sem) unweighted UniFrac distances comparing phylogenetic configurations of vaginal microbial communities of nonpregnant females, and control and stress-exposed females during early pregnancy and postpartum. N = 5–8/day/treatment; Mann Whitney U, *p < 0.05, ***p < 0.001. (c,d) Distribution of dominant taxa in control and stress-exposed females reveals an impact of stress on vaginal microbial community composition. (c) Left panel, mean phyla-level relative abundance of nonpregnant, pregnant, and postpartum vaginal samples. (d) Right panel, mean relative abundance of the most abundant taxa in vaginal communities sampled from nonpregnant, pregnant (E7.5), and postpartum (PN2) females. E7.5, gestation day 7.5; PN2, postnatal day 2; N = 5–8/day/treatment. (e) Vaginal Lactobacillus relative abundance is significantly decreased at PN2 in females exposed to stress compared with control females PN2, postnatal day 2; N = 5–8/day/treatment. Mann Whitney U, **p < 0.01.
Given the stress alterations to community structure, we next examined whether these alterations were lasting or transient and therefore inconsequential to offspring reprogramming. The phylum Proteobacteria dominated the nonpregnant and postpartum murine vagina regardless of treatment (Fig. 2c) and stress-exposed females failed to exhibit a bloom of the Firmicutes and Bacteroidetes phyla at E7.5 compared with control females (Fig. 2c). At the genus level, members of Aggregatibacter constituted the dominant vaginal taxa in nonpregnant and postpartum murine vaginas (Fig. 2c). The Firmicutes and Bacteroidetes bloom in control females at E7.5 was driven by an expansion of members of unclassified Lachnospiraceae, Clostridiales, S24-7, and Prevotella (Fig. 2c). Consistent with alpha diversity measures, stress-exposed females showed a further expansion of Proteobacteria at E7.5, predominantly by members of the Helicobacter genus (Fig. 2c). In addition, relative abundance of vaginal Lactobacillus was decreased at PN2 in females exposed stress (p = 0.0044, Mann Whitney U, Fig. 2d). These results show that stress during pregnancy affects bacterial community structure and composition within the vagina, suggesting the possibility that these dysbiotic communities are transmitted to offspring during birth.
Shared overlap between maternal and offspring communities is determined by offspring age
As intestinal microbiota acquisition begins at birth through exposure of the maternal microbiota, we next determined whether these primary colonizers shared similarities with bacterial communities within the maternal gut, vagina, or a combination. Unweighted UniFrac matrices were calculated to determine distances between offspring gut microbiota and communities within the maternal gut and vagina, and then visualized using PCoA. Consistent with previous reports in humans, maternal vaginal bacterial profiles showed the greatest overlap with neonate gut profiles, while maternal gut bacterial profiles showed the greatest overlap with gut profiles of offspring at weaning (Fig. 3a)33,35,37. Additional unweighted UniFrac distance comparisons revealed that the community structure in the colon of control offspring at PN2 was more similar to PN2 control maternal vaginal samples than PN2 EPS maternal vaginal samples (p = 0.0027, permutation test; Fig. 3b). Conversely, colonic community structure in PN2 EPS offspring was more similar to PN2 EPS maternal vaginal samples compared with PN2 control maternal vaginal samples (p = 0.0046, permutation test; Fig. 3b).
(a) Newborn microbiota profiles overlap maternal vaginal microbial profiles while the weaning microbiota are similar to maternal gut microbial communities, providing an important proof-of-concept that colonization and subsequent assembly of the murine microbiota follows a conserved program whereby the colonizing microbiome reflects the maternal vaginal microbiome, and become more similar to the maternal gut microbiome as offspring mature. Communities clustered using PCoA of the unweighted UniFrac distance matrix. The percentage of variation explained by the PC is indicated on the axes. Each point corresponds to a single community sample from a maternal body site or offspring colon, and color-coded by the indicated metadata. PCoA, principle coordinates analysis; PC, principle coordinates. N = 4–8 dams/site/day/treatment; N = 5–12 offspring/sex/day/treatment. (b) Average (±sem) unweighted UniFrac distance for pairwise comparisons between maternal gut and vaginal microbial communities at PN2. N = 4–8 dams/site/day/treatment; N = 5–12 offspring/sex/day/treatment. **p < 0.01. (c) Heatmap of mean relative abundances of 24-development discriminatory taxa identified by Random Forests. Cluster I corresponds to a microbial signature associated with early colonizers while Cluster II corresponds to a weaning microbial signature. Hierarchical clustering performed using the Spearman rank correlation distance matrix. See Methods for more detail on the Random Forests approach. (d) Random Forests identified development-discriminatory taxa are altered of EPS exposure. Differentially abundant taxa in colon samples from male and female control and EPS offspring at PN2, PN6, and PN28. See Supplemental Fig. 5B for additional comparisons. Barplots indicate median ± IQR. FCON, female control; MCON, male control; FEPS, female EPS; MEPS, male EPS; EPS, early prenatal stress; PN, postnatal day; IQR, Interquartile range. N = 5–12 offspring/sex/day/treatment; Mann Whitney U, *p < 0.05, **p < 0.001, ***p < 0.001. PC, principal coordinate; PCoA, principal coordinates analysis. PN, postnatal day; M, male; F, female. FCON, female control; MCON, male control; FEPS, female EPS; MEPS, male EPS; EPS, early prenatal stress; PN, postnatal day. N = 5–12 offspring/sex/day/treatment.
Effects of EPS on bacterial communities are sex-specific in offspring
Assembly of bacterial communities in the postnatal gut has been shown to parallel chronological age, and disruption in this synchrony has been associated with metabolic dysfunction49,53. As key features of the mouse model of early prenatal stress is male-specific metabolic dysfunction, we next determined assembly of bacterial communities is sexually dimorphic and altered by early prenatal stress. Colonic samples were obtained from control and EPS male and female pups at birth, during active breastfeeding, and weaning. Permutational multivariate analysis of variance using UniFrac distance measurements revealed a significant main effect of age (PERMANOVA pseudo-F = 17.61, p = 0.001), EPS by age interaction (PERMANOVA pseudo-F = 8.40, p = 0.001), and EPS by sex by age interaction (PERMANOVA pseudo-F = 4.33, p = 0.001). Post-hoc analysis revealed significant differences in the community structure of EPS males at PN28 compared with control males and EPS females at the same age (ps < 0.01, Supplemental Fig. 4B). Community diversity increased as a function of age in all offspring regardless of age and treatment across multiple measures of alpha diversity (Supplemental Fig. 4C).
Given the significant interaction between sex and chronological age on community structure, we next determined whether the assembly of bacterial communities differed between male and female offspring. To account for potential sex differences, two separate Random Forests models were regressed against chronological age. Constructing these age-discriminatory models for male and female offspring revealed that while many of the age-discriminatory taxa were shared between male and female offspring, their ranked importance in driving sex-specific maturation of the colonic bacteria differed (Supplemental Fig. 5A, Supplemental Table 2). Using this approach, we identified 27 age-discriminatory taxa that formed two distinct clusters: cluster 1 resembled an early life microbial community consisting of nine taxa, while cluster 2 resembled a weaning microbial community consisting of the remaining eighteen taxa (Fig. 3c,d). Further, eight age-discriminatory taxa exhibited significant sex differences (Supplemental Fig. 5B).
Using the microbial signatures identified from Random Forests analysis, we determined the sex-specific EPS effects on postnatal colonic microbiota maturation. Early colonizers were comprised of taxa within the Proteobacteria and Firmicutes phyla. At the genus level, Aggregatibacter, Paenibacillus, and Lactobacillus abundances were similar among male and female offspring, while sex differences in the relative abundance of Streptococcus, Haemophilus, and Sphingobium were observed (ps < 0.05, Mann Whitney U) (Supplemental Fig. 5B). EPS exposure disrupted the composition of the early colonizers Lactobacillus and Streptococcus independent of sex (ps < 0.05, Mann Whitney U) (Fig. 3d). These early colonizers were rapidly replaced by exclusive dominance of Lactobacillus and Aggregatibacter, regardless of sex and treatment (Fig. 3d, Supplemental Fig. 5B). By weaning, the abovementioned expansion in microbial diversity and richness during this period was related to the replacement of facultative anaerobes such as Aggregatibacter, Lactobacillus, and Streptococcus by a vast community of strictly obligate anaerobes (Fig. 3d, Supplemental Fig. 5B). Interestingly, sex differences also reemerged at weaning (e.g., PN28). Mucispirillum, Odoribacter, and Desulfovibrio were enriched in females at weaning, while Dehalobacterium and Flexispira were enriched at weaning in males (Fig. 3d, Supplemental Fig. 5B). Early prenatal stress exposure disrupted sex differences in Odoribacter, Desulfovibrio, Flexispira, and Mucispirillum abundance, while Lachnospiraceae and Clostridiales abundance were increased only in EPS males relative to control males (Fig. 3d, Supplemental Fig. 5B). In addition to the EPS sex-specific effects on the gut microbiota composition at weaning, a main effect of EPS on the relative abundance of the Bacteroidetes Rikenella was also observed (p < 0.001) (Fig. 3d, Supplemental Fig. 5B). Collectively, our results demonstrate that maternal stress during the first week of pregnancy shapes sex-specific maturation of the colonic microbiota.
To relate composition shifts to functional capacity, PICRUSt analysis was used to predict KEGG functional categories, and similar to the approach above, separate Random Forests regressions were constructed to identify functional shifts associated with maturation (Supplemental Fig. 6A). Thirty-four age-discriminatory pathways were identified (Supplemental Fig. 6A). Early colonizers were enriched in functional pathways related to metabolism and degradation of amino acids, ketone bodies, and the short chain fatty acids butyrate and propionate (Supplemental Fig. 6B). Despite the limited representation of microbial communities during the lactation (i.e., PN 6), numerous unique microbial metabolic pathways were enriched only during this period, including primary bile acid biosynthesis, ascorbate and aldarate metabolism, alpha-Linolenic metabolism, tryptophan metabolism, steroid hormone biosynthesis, and lysine degradation (Supplemental Fig. 6B). Consistent with the expansion of community diversity as a function of chronological age, bacterial functional capacity was characterized at weaning by extensive enrichment of carbohydrate and amino acid metabolic pathways (Supplemental Fig. 6B).
Discussion
Early life host-microbe interactions disproportionally contribute to later health outcomes including susceptibility and resistance to infection, predisposition to inflammatory and metabolic disorders, and changes to neural circuits that control stress and affect, yet remarkably little is known about the maternal factors that may influence this process19,54,55,56,57,58,59,60. Although maternal stress experience during pregnancy is associated with profound dysregulation of pathways along the gut-brain axis, far less is known regarding the direct impact of stress on maternal bacterial communities that, in turn, shape these health outcomes in offspring15,16,17,18,30,61,62,63,64,65. Thus, alterations to temporal and spatial dynamics of maternal gut and vaginal bacterial communities during pregnancy and postpartum may represent a novel mode of stress reprogramming of offspring outcomes.
To determine whether early gestation stress exhibits a lasting impact on maternal fecal bacterial communities across pregnancy, fecal samples were collected at daily intervals from the initiation of pregnancy to term. Time-series profiling revealed that the murine maternal fecal bacteria undergo remodeling over the course of pregnancy. Surprisingly, two distinct windows demarcate community structure, composition, and functional capacity during pregnancy that parallel distinct metabolic adaptations, such as weight gain, that occur during pregnancy20. One question raised by our observations relates to host factors that may contribute to maternal gut bacterial community temporal patterns. Complex hormonal, endocrine, and immune interactions in the pregnant female emerge between gestation day 8 and 11, suggesting factors specific to pregnancy may mediate community dynamics. Feto-placental vascularization occurs during this window, a process during which the placenta begins to invade the uterus, gains contact with maternal circulation, and substrate transport from maternal circulation to the fetus is initiated66. A possible explanation is that maternal and fetal signals converge upon certain groups of bacteria that provide a pool of substrates important for fetal growth and development23. In this regard, selection and subsequent stabilization of these microbial communities may represent an adaptation to ensure sequestration of substrates important for offspring development. Thus, the capacity of such signals to impose selective pressures upon maternal microbial composition and function during pregnancy is intriguing and warrants future study.
Our observations of lasting stress alterations to maternal gut bacterial diversity and composition reveal a novel mode of prenatal stress reprogramming. Bacterial diversity and composition change abruptly in response to stress, which may be related to previous observations that chronic stress exaggerates the normal physiological immunosuppression and metabolic syndrome that occurs during pregnancy67. Thus, stress-related changes to the nutritional, endocrine, or immune environment of the maternal gut might explain the simultaneous bloom of facultative anaerobes Mucispirillum (phylum Deferribacteres) and Desulfovibrionaceae (phylum Proteobacteria) at the expense of some obligatory anaerobes31,68,69. The combined capacity to degrade mucin (Mucispirillum) and production of hydrogen sulfide (Desulfovibrionaceae) may provide a putative link between stress during pregnancy and intestinal inflammation70,71,72. Taken together, these results suggest that stress during the first week of pregnancy induced a coordinated shift in the maternal fecal bacteria community structure, particularly during periods of pregnancy when community structure is most variable.
Dynamic shifts in the functional capacity of fecal bacterial communities parallel well-established metabolic adaptations in pregnancy, implicating pregnancy as an important period of host-microbe synchrony20. Accumulation of nutrients and lipid stores occurs during the early phase of pregnancy when offspring growth is limited20. Similarly, we observe microbial pathways involved in the biosynthesis of amino acids and lipids enriched during the early phase of pregnancy. To support the demand underlying the exponential growth phase of fetal development during the latter portion of pregnancy, maternal metabolic processes shift to a catabolic state that is characterized by accelerated lipolysis, gluconeogenesis, and ketogenesis20. The enrichment of predicted functional pathways involved in gluconeogenesis and metabolism of fatty acids during the late phase of pregnancy may highlight an important role of bacterial communities in modulating the metabolic demand of fetal growth and development. Enrichment of pathways involved in amino acid and glycerophospholipid metabolism observed during the early pregnancy window was disrupted by maternal stress experience. Similarly, upregulation of lysine biosynthesis pathways during late gestation was disrupted in females exposed to stress. In contrast, stress increased the abundance of pathways involved in biosynthesis of fatty acids and degradation of branched chain amino acids, suggesting the possibility that stress during pregnancy may induce an alternate metabolic state of the bacterial communities, that, in turn, may impact availability of metabolites. Building on these community-level analyses, additional mechanistic studies are now needed to identify stress-altered bacterial strains, the metabolites produced by these strains, and their importance in normal fetal growth and development.
In addition to maternal gut communities, stress may also alter vaginal bacterial communities that provide the primary inoculum during birth. To determine the lasting effect of maternal stress experience on vaginal bacterial communities, vaginal fluid samples were collected prior to mating, at day 7 of pregnancy, and 24 h postpartum. Surprisingly, Aggregatibacter of the Pasteurellaceae family within the phylum Proteobacteria dominated the communities of the mouse vagina. Although communities in the human vagina are dominated by the acidophilic lactobacilli, this appears to be largely impacted by vaginal pH: as vaginal pH shifts towards neutral, the communities reflect higher diversity and loss of Lactobacillus dominance without consequences on vaginal health73. Further, lactobacilli dominance does not appear to reflect the vaginal ecosystem of non-human primates and livestock, thus, it is not surprising that mouse vaginal communities were not dominated by Lactobacillus74,75. Nevertheless, the vaginal microbial composition shifted considerably during the first week of pregnancy, from near complete dominance of Aggregatibacter in nonpregnant females to a distinct community state type composed of Aggregatibacter, Prevotella, Clostridiales, and Lachnospiraceae. Dominance of Aggregatibacter reemerged in the early postpartum period, collectively highlighting temporal dynamics of the mouse vaginal microbiota during pregnancy. Building upon our previous work, we now show that stress during pregnancy alters bacterial community dynamics of the maternal vagina during pregnancy and the postpartum period18. Moreover, consistent with previous reports in rodents and humans, significant overlap was observed between the maternal vaginal and neonate colonic microbiota32,33,34,35,76. Similarities between mother-offspring microbiota were dependent on maternal body site and age of offspring, whereby neonate colonic microbiota resemble the maternal vaginal microbiota immediately following birth and become more similar to the maternal gut microbiota as a function of age33,34,35.
To determine whether maternal stress experience shapes assembly of offspring colonic microbiota, community composition was evaluated in offspring colon samples at birth, during active breastfeeding, and at weaning. Microbial diversity increased as a function of age and independent of sex or treatment. Conversely, the structure of bacterial communities appears to be mediated through a complex interaction between age, sex, and treatment. EPS exposure significantly altered community structure in males at weaning compared with control males, and this alteration parallels the emergence of key EPS phenotypes such as changes in metabolism, stress responsivity, and cognition25,26,27,45,46,47,48.
To further examine the lasting impact of maternal stress exposure on offspring microbial communities, we used Random Forests models to identify sex- and age-discriminatory taxa that define a normal program of microbial assembly. This approach showed that postnatal maturation in control mice followed a patterned progression, characterized by colonization of
As with any descriptive dataset, there are limitations here in making causal inferences. The current data has a reliance on data generated from 16S rRNA maker gene sequencing and predictive algorithms to drive conclusions regarding the impact of the maternal and offspring bacterial communities as a novel contributor to the EPS phenotype. Additional work that combines time-series profiling with bacterial strain-level resolution sequencing may identify temporal and spatial dynamics of bacterial strains and functional categories by stress in dam and offspring. Additionally, the current study does not provide functional links between these alterations in bacterial community composition and metabolite availability in dam and offspring. While previous work in our laboratory has demonstrated significant shifts in metabolite availability in EPS offspring at postnatal day 2, additional work is required to identify whether stress-mediated alterations to maternal and fetal metabolites during the prenatal window represent a causal link to key aspects of the EPS phenotype18.
With these important caveats in mind, our data provide novel insight into the role of chronic stress on prenatal and postnatal development via the maternal gut and vaginal microbiome, respectively. Stress during the first week of pregnancy induces rapid and lasting alterations to the microbial landscape within the maternal gut and vagina. Further, maternal stress experience impacts assembly of the postnatal colonic microbiome in a sex-specific manner and these alterations precede the emergence of key EPS phenotypes, such as stress pathway dysregulation, cognitive deficits, and metabolic changes in these males as adults. Collectively, our results suggest stress reprogramming of offspring development may occur via the maternal microbiota through two novel mechanisms. First, stress alterations to maternal gut microbiota composition and function during pregnancy may impact the available pool of microbe derived metabolites and substrates that are necessary for normal prenatal growth and development86. Second, transmission of stress-altered vaginal microbiota may disrupt neonate-microbe interactions necessary for normal immune education, metabolism, and neurodevelopment37,41,42,43,44. Gut colonization by these stress-altered vaginal microbiota may, in turn, prime the intestinal niche for successors that are associated with greater disease risk later in life, particularly those disorders that exhibit a strong sex bias in onset, severity, and treatment outcomes. Our results open possibilities for future studies on defining mechanisms of these associations, and producing studies in which causality can be demonstrated between the neonate microbiome and the developing brain.
Methods
Experimental animals
Mice used in these studies were on a mixed C57BL/6:129 background that are maintained in-house as an outbred mixed colony. A total of 30 mixed colony breeding pairs were used for these studies, and females were checked daily at 700 EST for copulation plugs. Dams were housed under a 12 h light/day photoperiod with a standard chow diet (Purina Rodent Chow, St. Louis, MO; 28.1% protein, 59.8% carbohydrate, 12.1% fat) and ad libitum access to water. Noon on the day that the plug was observed was considered to be embryonic day 0.5 (E0.5). Samples from one male and one female per litter were used for all subsequent analyses. All experiments were approved by the University of Pennsylvania Institutional Animal Care and Use Committee and performed in accordance with National Institutes of Health Animal Care and Use Guidelines.
Early prenatal stress
Administration of chronic variable stress was performed as previously described25,26,27,45,46,47,48. Dams were randomly assigned to treatment groups to receive stress during days 1–7 of gestation (EPS; n = 8) or to a control (n = 5) unstressed group. Following confirmation of a copulation plug, pregnant mice assigned to the EPS group experienced each of the following stressors on different days: 60 min of fox odor exposure (1:5,000 2,4,5-trimethylthiazole; Acros Organics), 15 min of restraint (beginning at 1300 EST) in a modified 50-mL conical tube, 36 h of constant light, novel noise (White Noise/Nature Sound-Sleep Machine; Brookstone) overnight, three cage changes throughout the light cycle, saturated bedding (700 mL, 23 °C water) overnight, and novel object (eight marbles of similar shape and color) exposure overnight. Stressors were selected to be non-habituating and did not induce pain or directly influence maternal food intake, weight gain, behavior, litter size, or sex ratio25,26,27,45,46,47,48.
Sampling strategy
The fecal, vaginal, and colonic high-resolution sampling schedule is shown in Fig. 1a. Following confirmation of the mating plug, females were randomly assigned to either control or stress (EPS) groups and daily fecal samples were collected at the same time each day to control for diurnal confounds. A vaginal fluid sample was collected following cessation of the chronic variable stress protocol (i.e., gestation day 7.5). Twenty-four hours after birth, an additional vaginal fluid sample was collected to assess similarity between maternal microbiota and pup gut microbiota profiles. On postnatal day (PN) 2, 6 and 28 colon samples from 1 male and 1 female offspring per litter were collected for sequencing. For these experiments in offspring, bacterial enumeration by 16S rRNA sequencing was conducted from colonic scrapes that include the loose colonic contents and mucus. These time points were selected to assess neonate microbial structure and function during key periods of assembly: colonization (PN2), breastfeeding (PN2), and weaning (PN28). Maternal fecal (n = 4–8 females per treatment sampled every day, yielding a total of 228 unique fecal samples), maternal vaginal (n = 5–8 /treatment/time point; 39 unique vaginal fluid samples in total), and offspring distal colon samples (n = 5–8 offspring/sex/treatment/time point, for a total of 80 colon samples) comprise the total number of samples used in this study. To control for potential confounds of the housing and the cage environment in the assessment of EPS effects on the microbiota, several precautions were taken: first, in the bacterial community composition assessment in maternal fecal and vaginal samples, females were singly housed during the entire duration of the pregnancy; and, second, in the bacterial communities assessment within offspring colonic samples, one male and one female offspring were collected at each time point from a different dam, such that the sample size correspond to cage number.
DNA isolation, amplification, and sequencing of the 16S rRNA gene
Fecal, vaginal fluid, and colonic samples were rapidly frozen in liquid nitrogen and subsequently stored at −80 °C before DNA extraction. Genomic DNA from fecal, vaginal and colon samples was isolated using the Stratec PSP Spin Stool DNA Plus kit using the difficult to lyse bacteria protocol from the manufacturer (STRATEC Molecular GmbH, Berlin, Germany). For each set of extractions, one blank swab exposed to laboratory air was processed as a negative control. The V4 region of the 16S rRNA gene was amplified using barcoded primers for the Illumina platform as previously described87. Sequencing was performed on a MiSeq instrument (Illumina, San Diego, CA) using 250 base paired-end chemistry at the University of Pennsylvania Next Generation Sequencing Core. Mock communities, consisting of genomic DNA from 12 ATCC strains, were included as a positive control and for run-to-run quality control. Quality filtering and chimera checking yielded 15,471,769 quality-filtered sequences with a mean ± SD depth of 24,326 ± 4,1875 sequences per sample.
16S rRNA gene analysis
Paired-end reads were assembled and quality filtered to include sequences with a Q score ≥30. mothur v. 1.36.1 was employed to remove sequences <248 bp and >255 bp in length and sequences with homopolymers >10 bp in length88. QIIME v. 1.8 was used for further downstream processing and analyses89. OTUs were defined using 97% sequence similarity with CD-HIT, and a representative sequence from each OTU containing ≥10 sequences was chosen for downstream analyses (based on the most abundant sequence). Chimeric sequences were removed using ChimeraSlayer. Representative sequences were assigned to genera using the Ribosomal Database Project (RDP) classifier v 2.2, multiple sequence alignment was performed using PyNAST, and a phylogeny was built with FastTree. The samples were rarified to 1,000 sequences per sample for calculating alpha- and beta-diversity metrics.
Bioinformatics quality assurance
For quality control purposes, water and processed blank samples were sequenced and analyzed through the bioinformatics pipeline. Taxa identified as cyanobacteria or ‘unclassified’ to the phylum level were removed. For further quality control assurance and to ensure run-to-run reproducibility, genomic DNA from abovementioned mock communities were sequenced and the expected sequences were compared to the obtained sequences.
Identifying microbial communities involved with pregnancy and offspring postnatal maturation using Random Forests
Random Forests, a machine-based-learning algorithm, was used to regress microbiota relative abundance in the time-series profiling of mothers and offspring against day of pregnancy or chronologic age, respectively. The R package ‘randomForest’, ntree = 1000 was used with a default mtry parameter of p/3 where p is the number of input 97%-identity OTUs (features)) as previously described49. This approach allows for the identification of taxa that discriminate between different periods of interest, including stages of pregnancy and age of development in offspring. OTU importance was ranked by the percent decrease in prediction accuracy of the model that occurs when that OTU is removed. To estimate the minimal number of top ranking discriminatory taxa required for prediction, fivefold cross-validation using the “rfcv” function implemented in the ‘randomForest’ package was applied over 100 iterations49. A rarefied OTU table at 1,000 sequences per sample served as input data. The model derived from control maternal and offspring microbiota was then applied to microbiota from stress-exposed mothers and male and female offspring. The R package “Heatplus” was used to plot OTU relative abundances as heatmaps.
Predicted metagenomic analysis based on 16s rRNA sequencing data
The predicted functional capacity of the maternal gut and vaginal microbiota, as well offspring gut microbiota, was predicted using phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt)52. This algorithm extended ancestral-state reconstruction algorithm to predict which gene families are present and then combines gene families to estimate the composite metagenome. The analysis includes the following steps. (i) Obtaining a closed reference OTU table; (ii) Normalization of the OTU table; (ii) Predicting the functions for metagenomes. The last step produces the actual metagenome functional predictions (KEGG Orthologs) for a given OTU table, which multiplies each normalized OTU abundance by each predicted functional trait abundance to produce a table of functions by samples. As our a priori hypotheses are directly related to the metabolic contribution of the microbiota to the host, the predicted metagenome dataset was filtered to only include KEGG Orthologs that mapped on to the “Metabolism” category. Random Forests were applied to this filtered metagenome dataset to identify functional categories that shift as a function of stage of pregnancy or postnatal age, and then impact of stress on the proportional representation of these predicted microbial metabolic pathways was assessed.
Additional statistical methods
Alpha diversity and beta diversity are reported as mean ± SEM. Taxonomic relative abundances are reported as the median ± interquartile range (IQR). Permutational multivariate analysis of variance using distance matrices (PERMANOVA) was used to analyzed effects of age, sex, treatment, and interactions on unweighted UniFrac distances90. Linear regression models were used to predict the best fit for community diversity as a function of day of pregnancy. Repeated measures were used to assess significant changes in alpha diversity across gestation, followed by pairwise comparisons adjusted for multiple comparisons. For pairwise comparisons, P values were calculated using the Kruskal-Wallis test, Wilcoxon rank-sum test, or Mann Whitney U in GraphPad.
Additional Information
How to cite this article: Jašarević, E. et al. Stress during pregnancy alters temporal and spatial dynamics of the maternal and offspring microbiome in a sex-specific manner. Sci. Rep. 7, 44182; doi: 10.1038/srep44182 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
Bale, T. L. et al. Early life programming and neurodevelopmental disorders. Biological psychiatry 68, 314–319, doi: 10.1016/j.biopsych.2010.05.028 (2010).
de Kloet, E. R., Sibug, R. M., Helmerhorst, F. M. & Schmidt, M. V. Stress, genes and the mechanism of programming the brain for later life. Neuroscience and biobehavioral reviews 29, 271–281, doi: 10.1016/j.neubiorev.2004.10.008 (2005).
Malkova, N. V., Yu, C. Z., Hsiao, E. Y., Moore, M. J. & Patterson, P. H. Maternal immune activation yields offspring displaying mouse versions of the three core symptoms of autism. Brain Behavior and Immunity 26, 607–616, doi: 10.1016/J.Bbi.2012.01.011 (2012).
Nestler, E. J. et al. Neurobiology of depression. Neuron 34, 13–25 (2002).
McEwen, B. S. & Sapolsky, R. M. Stress and cognitive function. Curr Opin Neurobiol 5, 205–216 (1995).
Golubeva, A. V. et al. Prenatal stress-induced alterations in major physiological systems correlate with gut microbiota composition in adulthood. Psychoneuroendocrinology 60, 58–74, doi: 10.1016/j.psyneuen.2015.06.002 (2015).
Moloney, R. D., O’Mahony, S. M., Dinan, T. G. & Cryan, J. F. Stress-induced visceral pain: toward animal models of irritable-bowel syndrome and associated comorbidities. Front Psychiatry 6, 15, doi: 10.3389/fpsyt.2015.00015 (2015).
Moloney, R. D., Stilling, R. M., Dinan, T. G. & Cryan, J. F. Early-life stress-induced visceral hypersensitivity and anxiety behavior is reversed by histone deacetylase inhibition. Neurogastroenterology and motility; the official journal of the European Gastrointestinal Motility Society, doi: 10.1111/nmo.12675 (2015).
O’Mahony, S. M. et al. Early life stress alters behavior, immunity, and microbiota in rats: implications for irritable bowel syndrome and psychiatric illnesses. Biological psychiatry 65, 263–267, doi: 10.1016/j.biopsych.2008.06.026 (2009).
Pusceddu, M. M. et al. N-3 Polyunsaturated Fatty Acids (PUFAs) Reverse the Impact of Early-Life Stress on the Gut Microbiota. PloS One 10, e0139721, doi: 10.1371/journal.pone.0139721 (2015).
Pusceddu, M. M. et al. n-3 PUFAs have beneficial effects on anxiety and cognition in female rats: Effects of early life stress. Psychoneuroendocrinology 58, 79–90, doi: 10.1016/j.psyneuen.2015.04.015 (2015).
Bailey, M. T. The contributing role of the intestinal microbiota in stressor-induced increases in susceptibility to enteric infection and systemic immunomodulation. Hormones and behavior 62, 286–294, doi: 10.1016/J.Yhbeh.2012.02.006 (2012).
Bailey, M. T. et al. Exposure to a social stressor alters the structure of the intestinal microbiota: Implications for stressor-induced immunomodulation. Brain Behavior and Immunity 25, 397–407, doi: 10.1016/J.Bbi.2010.10.023 (2011).
Bailey, M. T., Engler, H. & Sheridan, J. F. Stress induces the translocation of cutaneous and gastrointestinal microflora to secondary lymphoid organs of C57BL/6 mice. Journal of neuroimmunology 171, 29–37, doi: 10.1016/j.jneuroim.2005.09.008 (2006).
Bailey, M. T., Lubach, G. R. & Coe, C. L. Prenatal stress alters bacterial colonization of the gut in infant monkeys. Journal of pediatric gastroenterology and nutrition 38, 414–421 (2004).
Bailey, M. T. & Coe, C. L. Maternal separation disrupts the integrity of the intestinal microflora in infant rhesus monkeys. Developmental psychobiology 35, 146–155, doi: 10.1002/(Sici)1098-2302(199909)35:2<146::Aid-Dev7>3.0.Co;2-G (1999).
De Palma, G. et al. Microbiota and host determinants of behavioural phenotype in maternally separated mice. Nature communications 6, 7735, doi: 10.1038/ncomms8735 (2015).
Jasarevic, E., Howerton, C. L., Howard, C. D. & Bale, T. L. Alterations in the Vaginal Microbiome by Maternal Stress Are Associated With Metabolic Reprogramming of the Offspring Gut and Brain. Endocrinology 156, 3265–3276, doi: 10.1210/en.2015-1177 (2015).
Hsiao, E. Y. et al. Microbiota modulate behavioral and physiological abnormalities associated with neurodevelopmental disorders. Cell 155, 1451–1463, doi: 10.1016/j.cell.2013.11.024 (2013).
Herrera, E. Metabolic adaptations in pregnancy and their implications for the availability of substrates to the fetus. European journal of clinical nutrition 54 Suppl 1, S47–51 (2000).
Koren, O. et al. Host remodeling of the gut microbiome and metabolic changes during pregnancy. Cell 150, 470–480, doi: 10.1016/j.cell.2012.07.008 (2012).
Jasarevic, E., Morrison, K. E. & Bale, T. L. Sex differences in the gut microbiome-brain axis across the lifespan. Philos Trans R Soc Lond B Biol Sci 371, doi: 10.1098/rstb.2015.0122 (2016).
Gomez de Aguero, M. et al. The maternal microbiota drives early postnatal innate immune development. Science 351, 1296–1302, doi: 10.1126/science.aad2571 (2016).
Braniste, V. et al. The gut microbiota influences blood-brain barrier permeability in mice. Science translational medicine 6, 263ra158, doi: 10.1126/scitranslmed.3009759 (2014).
Howerton, C. L., Morgan, C. P., Fischer, D. B. & Bale, T. L. O-GlcNAc transferase (OGT) as a placental biomarker of maternal stress and reprogramming of CNS gene transcription in development. Proceedings of the National Academy of Sciences of the United States of America 110, 5169–5174, doi: 10.1073/pnas.1300065110 (2013).
Mueller, B. R. & Bale, T. L. Sex-specific programming of offspring emotionality after stress early in pregnancy. The Journal of neuroscience; the official journal of the Society for Neuroscience 28, 9055–9065, doi: 10.1523/JNEUROSCI.1424-08.2008 (2008).
Howerton, C. L. & Bale, T. L. Targeted placental deletion of OGT recapitulates the prenatal stress phenotype including hypothalamic mitochondrial dysfunction. Proceedings of the National Academy of Sciences of the United States of America 111, 9639–9644, doi: 10.1073/pnas.1401203111 (2014).
Gluckman, P. D. & Hanson, M. A. Maternal constraint of fetal growth and its consequences. Seminars in fetal & neonatal medicine 9, 419–425, doi: 10.1016/j.siny.2004.03.001 (2004).
Belkacemi, L., Nelson, D. M., Desai, M. & Ross, M. G. Maternal undernutrition influences placental-fetal development. Biology of reproduction 83, 325–331, doi: 10.1095/biolreprod.110.084517 (2010).
De Palma, G., Collins, S. M., Bercik, P. & Verdu, E. F. The microbiota-gut-brain axis in gastrointestinal disorders: stressed bugs, stressed brain or both? The Journal of physiology 592, 2989–2997, doi: 10.1113/jphysiol.2014.273995 (2014).
Lyte, M., Vulchanova, L. & Brown, D. R. Stress at the intestinal surface: catecholamines and mucosa-bacteria interactions. Cell Tissue Res 343, 23–32, doi: 10.1007/s00441-010-1050-0 (2011).
Palmer, C., Bik, E. M., DiGiulio, D. B., Relman, D. A. & Brown, P. O. Development of the human infant intestinal microbiota. PLoS biology 5, e177, doi: 10.1371/journal.pbio.0050177 (2007).
Dominguez-Bello, M. G. et al. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proceedings of the National Academy of Sciences of the United States of America 107, 11971–11975, doi: 10.1073/pnas.1002601107 (2010).
Dominguez-Bello, M. G. et al. Partial restoration of the microbiota of cesarean-born infants via vaginal microbial transfer. Nat Med 22, 250–253, doi: 10.1038/nm.4039 (2016).
Yatsunenko, T. et al. Human gut microbiome viewed across age and geography. Nature 486, 222–227, doi: 10.1038/nature11053 (2012).
Backhed, F. et al. 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 (2015).
Dominguez-Bello, M. G., Blaser, M. J., Ley, R. E. & Knight, R. Development of the human gastrointestinal microbiota and insights from high-throughput sequencing. Gastroenterology 140, 1713–1719, doi: 10.1053/j.gastro.2011.02.011 (2011).
Cox, L. M. et al. Altering the intestinal microbiota during a critical developmental window has lasting metabolic consequences. Cell 158, 705–721, doi: 10.1016/j.cell.2014.05.052 (2014).
Deshmukh, H. S. et al. The microbiota regulates neutrophil homeostasis and host resistance to Escherichia coli K1 sepsis in neonatal mice. Nat Med 20, 524–530, doi: 10.1038/nm.3542 (2014).
Ma, J. et al. High-fat maternal diet during pregnancy persistently alters the offspring microbiome in a primate model. Nature communications 5, 3889, doi: 10.1038/ncomms4889 (2014).
Blaser, M. J. & Falkow, S. What are the consequences of the disappearing human microbiota? Nature reviews. Microbiology 7, 887–894, doi: 10.1038/nrmicro2245 (2009).
Mueller, N. T., Bakacs, E., Combellick, J., Grigoryan, Z. & Dominguez-Bello, M. G. The infant microbiome development: mom matters. Trends in molecular medicine 21, 109–117, doi: 10.1016/j.molmed.2014.12.002 (2015).
LaRosa, P. S. et al. Patterned progression of bacterial populations in the premature infant gut. Proceedings of the National Academy of Sciences of the United States of America(2014).
Subramanian, S. et al. Cultivating healthy growth and nutrition through the gut microbiota. Cell 161, 36–48, doi: 10.1016/j.cell.2015.03.013 (2015).
Bronson, S. L. & Bale, T. L. Prenatal stress-induced increases in placental inflammation and offspring hyperactivity are male-specific and ameliorated by maternal antiinflammatory treatment. Endocrinology 155, 2635–2646, doi: 10.1210/en.2014-1040 (2014).
Mueller, B. R. & Bale, T. L. Impact of prenatal stress on long term body weight is dependent on timing and maternal sensitivity. Physiology & behavior 88, 605–614, doi: 10.1016/j.physbeh.2006.05.019 (2006).
Mueller, B. R. & Bale, T. L. Early prenatal stress impact on coping strategies and learning performance is sex dependent. Physiology & behavior 91, 55–65, doi: 10.1016/j.physbeh.2007.01.017 (2007).
Pankevich, D. E., Mueller, B. R., Brockel, B. & Bale, T. L. Prenatal stress programming of offspring feeding behavior and energy balance begins early in pregnancy. Physiology & behavior 98, 94–102, doi: 10.1016/j.physbeh.2009.04.015 (2009).
Subramanian, S. et al. Persistent gut microbiota immaturity in malnourished Bangladeshi children. Nature 510, 417–421, doi: 10.1038/nature13421 (2014).
Faust, K. et al. Microbial co-occurrence relationships in the human microbiome. PLoS computational biology 8, e1002606, doi: 10.1371/journal.pcbi.1002606 (2012).
Turnbaugh, P. J. et al. The human microbiome project. Nature 449, 804–810, doi: 10.1038/nature06244 (2007).
Langille, M. G. et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nature biotechnology 31, 814–821, doi: 10.1038/nbt.2676 (2013).
Nobel, Y. R. et al. Metabolic and metagenomic outcomes from early-life pulsed antibiotic treatment. Nature communications 6, 7486, doi: 10.1038/ncomms8486 (2015).
Round, J. L. & Mazmanian, S. K. The gut microbiota shapes intestinal immune responses during health and disease. Nature reviews. Immunology 9, 313–323, doi: 10.1038/nri2515 (2009).
Jostins, L. et al. Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature 491, 119–124, doi: 10.1038/nature11582 (2012).
Cryan, J. F. & Dinan, T. G. Mind-altering microorganisms: the impact of the gut microbiota on brain and behaviour. Nature reviews. Neuroscience 13, 701–712, doi: 10.1038/nrn3346 (2012).
Nicholson, J. K. et al. Host-gut microbiota metabolic interactions. Science 336, 1262–1267, doi: 10.1126/science.1223813 (2012).
Tremaroli, V. & Backhed, F. Functional interactions between the gut microbiota and host metabolism. Nature 489, 242–249, doi: 10.1038/nature11552 (2012).
Hooper, L. V., Littman, D. R. & Macpherson, A. J. Interactions between the microbiota and the immune system. Science 336, 1268–1273, doi: 10.1126/science.1223490 (2012).
Costello, E. K., Stagaman, K., Dethlefsen, L., Bohannan, B. J. & Relman, D. A. The application of ecological theory toward an understanding of the human microbiome. Science 336, 1255–1262, doi: 10.1126/science.1224203 (2012).
Foster, J. A. & McVey Neufeld, K. A. Gut-brain axis: how the microbiome influences anxiety and depression. Trends in neurosciences 36, 305–312, doi: 10.1016/j.tins.2013.01.005 (2013).
Mayer, E. A., Tillisch, K. & Gupta, A. Gut/brain axis and the microbiota. The Journal of clinical investigation 125, 926–938, doi: 10.1172/JCI76304 (2015).
O’Mahony, S. M., Hyland, N. P., Dinan, T. G. & Cryan, J. F. Maternal separation as a model of brain-gut axis dysfunction. Psychopharmacology 214, 71–88, doi: 10.1007/s00213-010-2010-9 (2011).
Desbonnet, L. et al. Effects of the probiotic Bifidobacterium infantis in the maternal separation model of depression. Neuroscience 170, 1179–1188, doi: 10.1016/j.neuroscience.2010.08.005 (2010).
Rhee, S. H., Pothoulakis, C. & Mayer, E. A. Principles and clinical implications of the brain-gut-enteric microbiota axis. Nat Rev Gastro Hepat 6, 306–314, doi: 10.1038/nrgastro.2009.35 (2009).
Todros, T. et al. Review: Feto-placental vascularization: a multifaceted approach. Placenta 32 Suppl 2, S165–169, doi: 10.1016/j.placenta.2010.12.020 (2011).
Lupien, S. J., McEwen, B. S., Gunnar, M. R. & Heim, C. Effects of stress throughout the lifespan on the brain, behaviour and cognition. Nature reviews. Neuroscience 10, 434–445, doi: 10.1038/nrn2639 (2009).
Winter, S. E. & Baumler, A. J. Why related bacterial species bloom simultaneously in the gut: principles underlying the ‘Like will to like’ concept. Cell Microbiol 16, 179–184, doi: 10.1111/cmi.12245 (2014).
Stecher, B. et al. Like will to like: abundances of closely related species can predict susceptibility to intestinal colonization by pathogenic and commensal bacteria. PLoS pathogens 6, e1000711, doi: 10.1371/journal.ppat.1000711 (2010).
Robertson, B. R. et al. Mucispirillum schaedleri gen. nov., sp. nov., a spiral-shaped bacterium colonizing the mucus layer of the gastrointestinal tract of laboratory rodents. Int J Syst Evol Microbiol 55, 1199–1204, doi: 10.1099/ijs.0.63472-0 (2005).
Berry, D. et al. Phylotype-level 16S rRNA analysis reveals new bacterial indicators of health state in acute murine colitis. The ISME journal 6, 2091–2106, doi: 10.1038/ismej.2012.39 (2012).
Carbonero, F., Benefiel, A. C. & Gaskins, H. R. Contributions of the microbial hydrogen economy to colonic homeostasis. Nature reviews . Gastroenterology & hepatology 9, 504–518, doi: 10.1038/nrgastro.2012.85 (2012).
Ravel, J. et al. Vaginal microbiome of reproductive-age women. Proceedings of the National Academy of Sciences of the United States of America 108 Suppl 1, 4680–4687, doi: 10.1073/pnas.1002611107 (2011).
Swartz, J. D. et al. Characterization of the Vaginal Microbiota of Ewes and Cows Reveals a Unique Microbiota with Low Levels of Lactobacilli and Near-Neutral pH. Front Vet Sci 1, 19, doi: 10.3389/fvets.2014.00019 (2014).
Yildirim, S. et al. Primate vaginal microbiomes exhibit species specificity without universal Lactobacillus dominance. The ISME journal 8, 2431–2444, doi: 10.1038/ismej.2014.90 (2014).
Pantoja-Feliciano, I. G. et al. Biphasic assembly of the murine intestinal microbiota during early development. The ISME journal 7, 1112–1115, doi: 10.1038/ismej.2013.15 (2013).
Mirpuri, J. et al. Proteobacteria-specific IgA regulates maturation of the intestinal microbiota. Gut microbes 5, 28–39, doi: 10.4161/gmic.26489 (2014).
Rakoff-Nahoum, S. et al. Analysis of gene-environment interactions in postnatal development of the mammalian intestine. Proceedings of the National Academy of Sciences of the United States of America 112, 1929–1936, doi: 10.1073/pnas.1424886112 (2015).
Markle, J. G. et al. Sex differences in the gut microbiome drive hormone-dependent regulation of autoimmunity. Science 339, 1084–1088, doi: 10.1126/science.1233521 (2013).
Yurkovetskiy, L. et al. Gender bias in autoimmunity is influenced by microbiota. Immunity 39, 400–412, doi: 10.1016/j.immuni.2013.08.013 (2013).
Rooks, M. G. et al. Gut microbiome composition and function in experimental colitis during active disease and treatment-induced remission. The ISME journal 8, 1403–1417, doi: 10.1038/ismej.2014.3 (2014).
Bolnick, D. I. et al. Individual diet has sex-dependent effects on vertebrate gut microbiota. Nature communications 5, doi: ARTN 450010.1038/ncomms5500 (2014).
Klein, S. L., Marriott, I. & Fish, E. N. Sex-based differences in immune function and responses to vaccination. Trans R Soc Trop Med Hyg 109, 9–15, doi: 10.1093/trstmh/tru167 (2015).
Oertelt-Prigione, S. The influence of sex and gender on the immune response. Autoimmun Rev 11, A479–485, doi: 10.1016/j.autrev.2011.11.022 (2012).
Nance, D. M. Sex differences in the hypothalamic regulation of feeding behavior in the rat. Adv Psychobiol 3, 75–123 (1976).
Jasarevic, E., Morrison, K. E. & Bale, T. L. Sex differences in the gut microbiome-brain axis across the lifespan. Philos Trans R Soc Lond B Biol Sci 371, 20150122, doi: 10.1098/rstb.2015.0122 (2016).
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K. & Schloss, P. D. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Applied and environmental microbiology 79, 5112–5120, doi: 10.1128/AEM.01043-13 (2013).
Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and environmental microbiology 75, 7537–7541, doi: 10.1128/AEM.01541-09 (2009).
Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7, 335–336, doi: 10.1038/nmeth.f.303 (2010).
Anderson, M. J. A new method for non-parametric multivariate analysis of variance. Austral Ecology 26, 32–56 (2001).
Acknowledgements
The authors thank the University of Pennsylvania Diabetes Research Center for use of the NGSC (P30-DK19525). The research reported in this publication was supported by a pilot award from the PennVet Center for Host-Microbial Interactions at the University of Pennsylvania and by the National Institutes of Mental Health under Award Numbers P50-MH099910, MH104184, MH091258, MH087597, MH073030 and MH108286 (T.L.B.); and F32MH109298 (E.J.).
Author information
Authors and Affiliations
Contributions
E.J. and T.L.B. designed experiments and drafted the manuscript. E.J. and C.D.H. conducted experiments. A.M. and D.P.B. assisted with marker gene sequencing and data analysis and visualization. All authors contributed equally to revision of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Jašarević, E., Howard, C., Misic, A. et al. Stress during pregnancy alters temporal and spatial dynamics of the maternal and offspring microbiome in a sex-specific manner. Sci Rep 7, 44182 (2017). https://doi.org/10.1038/srep44182
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep44182
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.