Ancestral exposure to stress epigenetically programs preterm birth risk and adverse maternal and newborn outcomes
© Yao et al.; licensee BioMed Central Ltd. 2014
Received: 10 January 2014
Accepted: 1 July 2014
Published: 7 August 2014
Chronic stress is considered to be one of many causes of human preterm birth (PTB), but no direct evidence has yet been provided. Here we show in rats that stress across generations has downstream effects on endocrine, metabolic and behavioural manifestations of PTB possibly via microRNA (miRNA) regulation.
Pregnant dams of the parental generation were exposed to stress from gestational days 12 to 18. Their pregnant daughters (F1) and grand-daughters (F2) either were stressed or remained as non-stressed controls. Gestational length, maternal gestational weight gain, blood glucose and plasma corticosterone levels, litter size and offspring weight gain from postnatal days 1 to 30 were recorded in each generation, including F3. Maternal behaviours were analysed for the first hour after completed parturition, and offspring sensorimotor development was recorded on postnatal day (P) 7. F0 through F2 maternal brain frontal cortex, uterus and placenta miRNA and gene expression patterns were used to identify stress-induced epigenetic regulatory pathways of maternal behaviour and pregnancy maintenance.
Progressively up to the F2 generation, stress gradually reduced gestational length, maternal weight gain and behavioural activity, and increased blood glucose levels. Reduced offspring growth and delayed behavioural development in the stress cohort was recognizable as early as P7, with the greatest effect in the F3 offspring of transgenerationally stressed mothers. Furthermore, stress altered miRNA expression patterns in the brain and uterus of F2 mothers, including the miR-200 family, which regulates pathways related to brain plasticity and parturition, respectively. Main miR-200 family target genes in the uterus, Stat5b, Zeb1 and Zeb2, were downregulated by multigenerational stress in the F1 generation. Zeb2 was also reduced in the stressed F2 generation, suggesting a causal mechanism for disturbed pregnancy maintenance. Additionally, stress increased placental miR-181a, a marker of human PTB.
The findings indicate that a family history of stress may program central and peripheral pathways regulating gestational length and maternal and newborn health outcomes in the maternal lineage. This new paradigm may model the origin of many human PTB causes.
KeywordsPreterm birth maternal stress prenatal stress transgenerational inheritance microRNA epigenetic regulation gestation maternal health behavioural development perinatal programming pregnancy
Preterm birth (PTB), which is associated with an intrauterine pro-inflammatory state, represents the leading cause of neonatal morbidity and mortality and one of the most critical factors for disease in later life. For example, infants born preterm, that is, born before 37 completed weeks of pregnancy, are at greater risk for mortality, developmental delay and health conditions than infants born at term . In spite of focused research efforts considering the drastic impact of PTB on health outcomes, in more than 50% of cases the causes of PTB remain unknown.
It has been difficult to demonstrate a clear causal relationship in humans , although PTB has been recognized as a consequence of severe maternal distress during pregnancy , or due to preconceptional factors . Earlier reports suggested that adverse perinatal programming by stress may increase the risk of PTB and low birth weight ,,,. Cumulative effects of stress seem to be of particular importance to PTB risk , which may include repeated exposure to stress across generations. Notably, elevated PTB risk was noted to propagate through generations , suggesting that factors determining PTB risk factors may be passed on to the offspring through the maternal lineage. Recent studies focusing on transmission through the male germ line in rodents have suggested that altered stress responses and associated emotional traits are linked to ancestral exposure to environmental toxins  and stressful experiences –. Furthermore, prenatal exposure to endocrine disruptors in female rats , or to maternal undernutrition in humans  have been associated with increased metabolic and endocrine disease risk in the offspring.
The molecular mechanisms leading to stress-induced pathologies in the maternal lineage occur through two different mechanisms. One mechanism occurs through direct exposure of fetal somatic cells in the female F1 and F2 generations ,. Alternatively, if phenotypic changes persist in the non-exposed F3 generation, truly transgenerational mechanisms involve changes in the germline that involve epigenetic mechanisms ,. Possible mechanisms of transgenerational transmission may be linked to a stress-associated epigenotype involving microRNAs (miRNAs) that are replicated in subsequent generations. MicroRNAs (miRNAs) are reasonable candidates for such a role since they are differentially regulated by progesterone during myometrial quiescence and initiation of parturition ,.
Here, we proposed that maternal stress or the cumulative effects of recurrent stress influence PTB risk and poor health outcomes across three generations. Using rats, we show that PTB risk, metabolic, endocrine and behavioural outcomes are affected by a single exposure to prenatal stress in one generation. In addition, the findings indicate that recurrent prenatal stress across multiple generations amplifies hypothalamic-pituitary-adrenal (HPA) axis responses to exacerbate variations in gestational lengths and adverse outcomes. We also show that stress-modulated gestational length is accompanied by miRNA expression changes and altered target gene pathways in somatic cells in F1 and F2 generations. Our data suggest that epigenomic programming of PTB risk factors may be an important mechanism involved in adverse pregnancy outcomes and altered maternal and offspring behaviours.
Four-hundred-and-eight Long-Evans hooded rats (Rattus norvegicus), bred and raised at the local vivarium, were used. Nulliparous, pair-housed female rats between the age of 100 and 160 days underwent timed pregnancy by being individually paired with a male for one hour per day until mating occurred. Rats were weighed prior to pregnancy and gestational weight was recorded daily. Pregnancy of the rats was confirmed by steady weight gain. Pregnant rats were housed individually from gestational day (GD) 19 until delivery. The hours of gestational length were monitored by an infrared security system (CCTV Cameras, Panasonic, Newark, NJ, USA). Pups remained with the dams until weaning on postnatal day (P) 21.
The rats were housed under a 12-hour light/day cycle with lights on at 7:30 AM. All procedures were performed in accordance with the guidelines of the Canadian Council for Animal Care and approved by the local Animal Welfare Committee.
Tissues for miRNA and mRNA expression analyses were collected from representative dams (F0-N, F0-S, F1-NN, F1-SN, F1-SS, F2-SSS, n = 3 per group) on GD21 (placenta) and after weaning of their offspring on LD21 (brain and uterus). From all other dams, uterine tissues were collected on LD21 to count the embryonic implantation sites to account for potential embryonic loss. Left and right uterine horns were pooled for further analyses.
Offspring were sexed at P1 and weighed on P1, P7, P15 and P30. Matched groups of two to three offspring from each sex in each litter were tested in sensorimotor behaviour on P7. Groups of offspring included F1-NN non-stress controls (n = 17), F1-SN stress animals (n = 48), F2-NNN (n = 10), F2-SNN (n = 33), F2-SSN (n = 36), F3-NNNN (n = 10), F3-SNNN (n = 31), F3-SSNN (n = 79) and F3-SSSN (n = 88) animals. Thus, ‘N’ was added to describe the tested offspring generations. Vivarium housing concerns prevented the analysis of pregnancy outcomes for the F3 generation. For the data shown in the figures, male and female animals were combined. All behavioural tests were performed by experimenters blind to the experimental groups.
Timed-pregnant rats were stressed daily from GD 12 to GD 18 by restraint and forced swimming. Restraint of the body for 20 minutes occurred between 8:00 and 9:00 AM. Animals were placed in a customized transparent Plexiglas container for a period of 20 minutes each day . The container had perforated ends to allow ventilation. The inner diameter of the container was adjusted to the size of the animals (6 cm inner diameter or larger) to prevent turning and maintain the animals in a standing position without compression of the body. Forced swimming occurred in a round water tank (45 cm diameter, 77 cm high, filled up to 50 cm with 21°C water) for five minutes between 4:00 and 5:00 PM .
Analysis of maternal postpartum behaviour
Postpartum tail chasing behavior was scored as the time spent engaged with the tail and the number of rotations was recorded from 24-hour infrared video recordings ,. Behaviour was videorecorded using an infrared video surveillance system (Panasonic WV-BP330, Panasonic, Minato-ku, Tokyo, Japan). Maternal behaviour was analysed for the first 60 minutes following the delivery of the last pup based on the video recorded data. The amount of time spent engaged in tail chasing and the total number of rotations performed were measured. Initiation of tail chasing behaviour was scored when the dam took interest in her tail followed by chasing or holding the tail with the mouth. Completion of a tail chasing event was scored once the rat disengaged with her tail and initiated a different activity . Tail grooming was not included in this analysis.
To test proprioceptive, musculoskeletal, and vestibular development on P7, offspring were placed head facing down on a custom-made Plexiglas 40° incline wedge covered with foam tack . Animals were video recorded for one minute and then returned to their mother. The time spent in the downward position until the initiation of a turn was recorded and averaged for three trials.
Blood samples (0.6 ml) were collected from the tail vein on GD18 and LD 1 in mothers between 8:00 and 9:00 AM under 4% isoflurane anesthesia . Blood glucose was measured using an Ascensia Breeze Blood Glucose Meter (Bayer, Toronto, ON, Canada) with test strips. The remaining blood was transferred to centrifuge tubes and plasma was obtained by centrifugation at 10,000 rpm for eight minutes. The samples were stored at –20°C. Plasma corticosterone (CORT) levels were determined by enzyme-linked immunosorbent assay (ELISA) using commercial kits (Cayman Chemical, Ann Arbor, MI, USA).
Brain, uterus and placenta
Dams received an intraperitoneal overdose of pentobarbital (Euthansol 100 mg/kg; CDMV Inc., Saint-Hyacinthe, QC, Canada). After rapid decapitation, tissues were dissected and flash-frozen for miRNA and transcriptomic analysis. Maternal brain and uterine tissues (n = 3/group) were collected at the time of weaning (three weeks postpartum). Placenta from female offspring was collected from dams (n = 3/group) on GD21.
mRNA and microRNA expression analysis
RNA extraction and microarrays
Total RNA was extracted using TRI Reagent Solution (Applied Biosystems, Foster City, CA, USA). Microarray assay was performed for F0-N, F0-S and F2-SSS frontal cortices using a service provider (LC Sciences, Houston, TX, USA). The assay started from 4 to 8 μg total RNA sample, which was size-fractionated using a YM-100 Microcon centrifugal filter (Millipore, Bedford, MA, USA) and the small RNAs (<300 nt) isolated were 3′-extended with a poly(A) tail using poly(A) polymerase. An oligonucleotide tag was then ligated to the poly(A) tail for later fluorescent dye staining; two different tags were used for the two RNA samples in dual-sample experiments. Hybridization was performed overnight on a μParaflo microfluidic chip using a micro-circulation pump (Atactic Technologies, Houston, TX, USA) ,. On the microfluidic chip, each detection probe consisted of a chemically modified nucleotide coding segment complementary to target miRNA or other RNA (control sequences) and a spacer segment of polyethylene glycol to extend the coding segment away from the substrate. The detection probes were made by in situ synthesis using photogenerated reagent (PGR) chemistry. The hybridization melting temperatures were balanced by chemical modifications of the detection probes. Hybridization used 100 μL 6xSSPE buffer (0.90 M NaCl, 60 mM Na2HPO4, 6 mM ethylenediaminetetraacetic acid (EDTA), pH 6.8) containing 25% formamide at 34°C. After RNA hybridization, tag-conjugating cyanine 3 (Cy3) and cyanine 5 (Cy5) dyes were circulated through the microfluidic chip for dye staining. Fluorescence images were collected using a laser scanner (GenePix 4000B, Molecular Device, Sunnyvale, CA, USA) and digitized using Array-Pro image analysis software (Media Cybernetics, Rockville, MD, USA). Data were analysed by first subtracting the background and then normalising the signals using a LOWESS filter (locally-weighted Regression) . For two-color experiments, the ratio of the two sets of detected signals (log2 transformed, balanced) and P-values of the t-test were calculated. Differentially detected signals were those with P-values of less than 0.10.
The putative gene targets for miRNAs were searched by computational analysis (TargetScan, Whitehead Institute for Biomedical Research MIT, Cambridge, MA, USA), which generated a list of predicted gene targets and related biological processes.
Quantitative real time PCR
In order to validate miRNAs, we performed quantitative real time PCR (qRT-PCR) analysis of these differentially regulated miRNAs (n = 3 per group for F0, F1 and F2 generations, three replicates per sample): miR-23b, miR-96, miR-141, miR-181a, miR-182, miR-183, miR-200a, miR-200b, miR-200c, miR429 and miR-451. Sno202, U6 and 5 s rRNA were used as references for calculation of the expression ratio. Reverse transcription oligos and amplification primers were designed according to an established protocol . The same samples of total RNA used for microarray analysis were used for qRT-PCR analysis. The generation of cDNAs from the total RNA samples was performed using M-MuLV Reverse Transcriptase, NEB#M0253S (New England Biolab, Ipswich, MA, USA; see Additional file 1: Table S1 for reverse transcription primers). For mRNA quantification, the cDNA were synthesized by using iScript cDNA synthesis kit (Bio-Rad, Mississauga, ON, Canada) following the supplier’s instructions. qRT-PCR reactions were conducted with Bio-Rad CFX96™ Real-Time PCR Systems, using SsoFas™ EvaGreen® Supermix (Bio-Rad) reaction premix added to the cDNAs templates and specific primers [see Additional file 1: Table S1 for primer sequences]. A total volume of 12 μl reaction mix was used, with 2.5 μl of cDNA template, 400 nM forward primer, 400 nM reverse primer and 6 μl of SsoFast™ EvaGreen® Supermix (Bio-Rad).
Phenotypic data, including gestational length, body weight, litter size, glucose, CORT and behavioural data, were analysed using a repeated measures analysis of variance (ANOVA) (group x sex) followed by two-way ANOVA (group x sex) at the different time points. Significant ANOVA results were explored further using either the post-hoc Fisher’s least significant difference (LSD) test or the Scheffe test for multiple comparisons. Plasma CORT values were transformed to normality. In addition, correlation analyses used Fisher’s R to Z transformations and Z-tests to calculate correlation coefficients. For all phenotypic data, a P-value of less than 0.05 was chosen as the significance level. Analyses were performed using Statview software version 5.0 (SAS Institute, Cary, NC, USA).
For miRNA microarray data t-values were calculated with P-values below a critical P-value (<0.10) selected for cluster analysis, which used a hierarchical method, average linkage and Euclidean distance metric . ANOVA was performed using Bio Rad CFX Manager for validation of miRNA expression by qRT-PCR. All data are presented as mean ± standard error of the mean (SEM).
Prenatal and multigenerational stress shortened gestational length across subsequent generations
Compared to the gestational length in non-stress controls (540.37 ± 3.8 hours), stress in the F0-S group did not significantly alter gestational length (537.38 ± 3.3 hours). Both a single generational or multigenerational exposure to prenatal stress modulated gestational length (F(6,34) = 3.48, P <0.05). The experience of prenatal stress significantly shortened gestational length in the F1 generation compared to non-stress controls (F(1,15) = 4.48, P <0.05; Figure 1B). Of interest, prenatal stress reduced gestational length in F1-SN animals (533.90 ± 2.95 hours) compared to non-stress controls (P <0.05). In subsequent generations, gestational length was further shortened by a stressful pregnancy in the F1-SS group (529.30 ± 0.87 hours, P <0.01) and among all groups that were exposed to prenatal and/or gestational stress (F(2,23) = 6.75, P <0.05; Figure 1B) in the F2 generation. In particular, non-stressed F2 dams whose grandmothers were stressed during gestation (F2-SNN; 530.10 ± 1.1) had a significantly shorter gestational length compared to that of non-stress controls (P <0.001). F2 dams whose grandmothers and mothers were stressed (F2-SSN; 528.2 ± 1.42 hours) or were stressed in each generation (F2-SSS; 527.78 ± 0.9 hours) also had shorter gestational lengths than those of non-stress controls (P <0.001).
Prenatal stress elevated blood glucose levels
The reductions in gestational lengths were accompanied by altered blood glucose concentrations in late pregnancy and postpartum. There was an overall effect of stress on glucose levels (F(3,33) = 3.71, P <0.05). In the F2 generation, SSN animals had higher basal blood glucose levels than non-stress F2-NNN (P <0.05) and F2-SSS (P <0.05) rats on GD 18 (Figure 1C). Furthermore, non-stressed dams whose grandmothers and mothers were stressed (F2-SSN) displayed higher blood glucose levels on GD18 (5.06 ± 0.25) compared to F2-NNN controls (4.0 ± 0.39, P <0.05; Figure 1C). Multigenerational stress in F2-SSS animals caused elevated gestational blood glucose levels on GD18 when compared to F2-NNN rats (5.19 ± 0.36, P <0.05) and after birth on LD 1 compared to F2-SNN and F2-SSN rats (7.0 ± 0.22, all Ps <0.05). Plasma CORT levels revealed that F2-SSS dams on GD18 had higher CORT levels than any F0 parental group (F2-SSS versus F0-S P <0.001, F2-SSS versus controls P <0.01).
Stress reduced gestational weight gain
Prenatal and transgenerational stress impeded offspring growth trajectories
In the developing offspring, prenatal stress delayed growth trajectories. There was an effect of group (F(1,8) = 6.176, P <0.001) and sex (F(1,8) = 64.756, P <0.001), but no interaction between the two factors. The effects of prenatal stress on offspring weight in the F1 generation became evident by P7 (F(1,64) = 57.97, P <0.0001), with F1-SS offspring (n = 48) being 3 g lighter than F1-NN controls (n = 17, P <0.001). This difference in weight remained significant into adulthood. Within each group and at all ages, female offspring were always lighter than male offspring (P <0.01).
Notably, the effects of prenatal stress on body weight were not noticeable at P1 until the F3 generation (see Figure 2D). All stress-treated F3 offspring groups were different from controls (F(3,241) = 5.12, P <0.001). Hence, the F3-SNNN (n = 31), F3-SSNN (n = 79) and F3-SSSN (n = 88) groups showed significant weight reduction of about 0.5 g compared to F3-NNNN animals (n = 44, all Ps <0.001). In F3-SNNN and F3-SSNN groups these effects remained significant throughout P7 (P <0.001), while F3-SSSN animals were not different from their F3-NNNN peers. At P15, this pattern was quite similar, revealing that the growth rate of F3-SSNN (P <0.001) and F3-SNNN (P <0.05) groups stayed behind that of F3-SSSN or F3-NNNN animals. Further, transgenerationally stressed F3-SNNN animals were lighter than their F3-SSNN counterparts (P <0.001) while multigenerationally stressed F3-SSSN animals were not different from F3-NNNN rats. At P30, F3-SSNN offspring were still lighter than F3-SSSN or F3-NNNN groups (P <0.001). Within all groups and at all ages, females were always lighter than males (P <0.01). Correlation analysis revealed that there was no influence of litter size on developmental trajectories.
Prenatal stress modified maternal behaviour across generations
Prenatal and transgenerational stress delayed offspring sensorimotor development
Proprioceptive, musculoskeletal, and vestibular development on P7 revealed an overall difference between groups (F(8,341) = 20.39, P <0.001). There was no effect of sex and no interaction between group and sex. Prenatally stressed F1-SN pups (n = 48), compared to non-stressed F1-NN pups (n = 17), showed a significantly delayed turning response (Figure 3C) when placed on an inclined plane (P <0.001; Figure 3D). The F2-SNN (n = 33) and F2-SSN pups (n = 36) showed further prolonged latencies compared to their F2-NNN counterparts (n = 10; P <0.001). The F3 generation showed a longer latency after grandmaternal stress in F3-SNNN animals (n = 31; P <0.001) and in F3-SSNN animals whose grandmothers and mothers were stressed (n = 79; P <0.001) or in F3-SSSN animals where all three generations experienced stress (n = 88; P <0.001; Figure 3D). These observations indicate that sensorimotor impairments resulted in slower response times.
MicroRNA (miRNA) profiles were altered in the F2 progenies of stressed animals
Target genes of altered miRNAs in brain include transcription regulators and mediators of neuropsychiatric disorders and endocrine pathways
Based on the validated list of altered miRNAs in the brain (frontal cortex), we compiled putative target genes using the mirSVR predicted target site scoring method . To explore the biological processes involving the putative gene targets, functional classification tools from DAVID  were used. Results are summarized in Supplementary Material [see Additional file 1: Table S1]. Between 17.1% to 23.7% of putative target genes were transcription regulators and an additional 2.3% to 5.5% were related to chromatin organization. Furthermore, a significant number of target genes assume roles in the genesis, reception or processing of endocrine functions including hormones, insulin, vitamins, carbohydrates, nutrients and drugs, or in embryonic development.
Possible involvement of miRNAs in disease pathways was suggested by the classification of putative gene target lists using PANTHER . Supplementary material [see Additional file 2: Table S2] summarizes the potential of these miRNA target genes that affect known disease pathways of metabolic, physiological, inflammatory, immunological, oncological, developmental and neuropsychiatric disorders.
Multigenerational programming by stress modulated uterine miRNA and gene expression involved in preterm birth
Stress across generations modified placental microRNA predictors of preterm birth
MiR-181a, which is altered in placentas of human preterm birth , remained unaltered in stressed F0 and F1 generations (Figure 5C). In the F2-SNN and F2-SSS groups, however, miR-181a was significantly upregulated compared to F2-NNN animals (n = 3, P <0.001 and P <0.01, respectively; Figure 5C), indicating programming by the cumulative effects of stress.
In about half of human cases, the causes of PTB remain unknown. Here, we provide evidence that gestational stress across generations of timed-pregnant rats has downstream effects on endocrine, metabolic and behavioural manifestations of PTB, and leads to shortened gestational length. Developmental trajectories across all generations of offspring were affected as early as P7. In terms of molecular mechanisms, stress in the parental F0 generation had minor effects on regulatory miRNA pathways in brain, uterus and placenta. By contrast, a history of stress in the F2 generation was associated with drastic changes in somatic tissue miRNA profiles and altered expression of genes that have been linked to PTB in humans. Notably, genuine transgenerational programming of developmental trajectories was observed in the F3 generation, in which gestational stress was imposed on the great-grandmaternal generation and was inherited to affect the developing embryo. These findings suggest that the mechanisms involved in the timing of parturition and associated behavioural and physiological signatures can be programmed through the maternal lineage.
A main finding of the present study is that gestational length is influenced by prenatal stress rather than by gestational stress. The impact of prenatal stress on phenotype in the present study was illustrated by delayed developmental milestones in the F1 generation, with an even stronger impact in subsequent F2 and F3 generations. Prenatal stress has been shown to program fetal brain development, HPA axis function and mental health –. Thus, early experiences can prime physiological and immunological processes that may lead to variations in gestational length , and susceptibility to altered glucose metabolism, such as type 2 diabetes  in adulthood. The timing and severity of the stressor is critical in that stress in early pregnancy may have greater effects on health outcomes than stress experienced in the last trimester . The present study induced stress from gestational days 12 to 18, representing a period thought to cover a large extent of the human second trimester , and in rats has been shown to be particularly susceptible to environmental influences, inflammatory processes and stress . While the present study controlled for adverse effects of stress, other factors, such as anesthetic administration, may still have affected the present outcomes . Altogether, it is possible that gestational and inter-generational programming of HPA axis responses may sensitize the response to environmental adversity thus resulting in gradually shortened gestation across generations and further reduction in the multigenerationally stressed cohorts (SSN and SSS).
Beyond fetal endocrine programming by an altered gestational endocrine milieu, maternal distress during pregnancy may also critically affect offspring brain development and physiology through variation in maternal behaviours ,,. The present data show that prenatal stress alters patterns of early postpartum maternal behaviours, which may be predictive of altered maternal care and stress coping at later times. The first hour after completed parturition may represent a critical transition phase in which the characteristic patterns of late antepartum behaviours, including tail chasing activity and nest building, convert into maternal care of the offspring . It has been demonstrated that the early postnatal environment, such as variations in maternal care, determines developmental and epigenetic outcomes ,. Patterns of altered maternal behaviour may transmit to subsequent generations ,. The contribution of endocrine and behavioural influences to generational programming is complex and likely reciprocally regulated by the epigenome.
In line with previous findings of stress-induced alteration in brain miRNA profiles  and according to the present observations of altered maternal behaviour, multigenerational stress in the F2-SSS group had prominent effects on miRNA expression patterns in the frontal cortex. Interestingly, F2-SSS dams showed upregulated miR-23b, which regulates oligodendrocyte development and myelination . miR-200 family members, including the downregulated miR-200a, are predicted to target genes that regulate synaptic function, neurodevelopment and neuronal survival . Stress also downregulated miRNAs that possess potential roles in the pathogenesis of psychiatric diseases, such as miR-96 , miR-182 and miR-183 . Furthermore, stress-induced downregulation concerned miR-429, which potentially influences development by altering cell proliferation and apoptosis . It is important to note that many neurodegenerative and psychiatric disorders share a pathology involving miRNA regulation ,, and that these miRNAs may in turn regulate central stress responses . Though not determined in this study, it is likely that these miRNA changes are not limited to the prefrontal cortex and uterus, thus indicating a potential intersection linking psychological stress to altered gestational length.
Across all generations, mechanisms of prenatal stress to modulate gestational length may include modulation of the complex pro-inflammatory state leading to PTB . Furthermore, stress may affect levels of hormones and neuropeptides, including prolactin, progesterone and oxytocin, which are involved in maintenance of pregnancy and timing of delivery . Increased fetal HPA axis activity may induce prostaglandin production by fetal membranes and decidua leading to uterine activation . In addition, stress may stimulate cytokines, which regulate the activity of placental 11-beta-hydroxysteroid dehydrogenase  to elevate PTB risk. These endocrine regulations have led to the notion that PTB risk may have roots in childhood . The present data confirm this notion and provide possible mechanistic links to epigenetic regulation of gene expression related to PTB risk.
Including the downregulated miR-200b, the miR-200 family may exert peripheral effects to control uterine quiescence and contractility during pregnancy and labour . Interestingly, miR-200b/200c/429 are induced at term labour in mice and humans and miR-200b/200c/429 are upregulated in mouse models of preterm labour . This group of miRNAs may largely interact with the endocrine cascade involved in pregnancy maintenance and termination, including progesterone and oxytocin . Moreover, miR-451 is expressed in the uterus , and regulated by estrogen and progesterone .
Target genes of the miR-200 family include three particular genes, Stat5b, Zeb1 and Zeb2, all involved in pregnancy maintenance . In the uterus, all three were downregulated by multigenerational stress in the F1 generation. Effects on Zeb2 expression were transmitted to the F2 generation. These findings concur with the reduction in gestational length. Accordingly, a decrease in Stat5b expression was linked to reduced progesterone activity and the initiation of labor, in particular in preterm birth . Furthermore, ZEB1 serves as transcription factor to inhibit the miR-200 family, thus enhancing Stat5b expression . As the myometrium transitions to term or preterm labor, reduced progesterone activity decreases ZEB1 and ZEB2 levels via a feed-forward mechanism ,, thus regulating the timing of parturition. The upregulation of uterine miR-200b may be causative for the suppression of Stat5b and ZEB1 and ZEB2; however, they may also reflect low postpartum progesterone levels due to timing of tissue sampling in the present study. Although the direction of these and the placental miR-181a changes are opposed to the downregulation found in human preterm birth , their differential expression across generations coincides with shortened gestational length and indicates a causal or, at least, predictive signature of preterm birth.
A role for genuine epigenetic inheritance of stress response is suggested by the present findings concerning the F3 generation. While context-dependent programming may have mainly determined the F1 and F2 phenotype, programming of the germ-line became evident by altered development in the F3 generation . Indeed, the most dramatic impact of prenatal stress on developmental trajectories was found in the F3 generation. Only in the F3 generation did the offspring display low body weight already on P1, which was associated with reduced growth trajectories and a drastic sensorimotor behaviour deficit. Since these phenotypic changes persisted to the F3-SNNN generation in the absence of direct somatic exposure, they are arguably mediated by genuine transgenerational programming of the female germline ,,. This suggests transgenerational epigenetic inheritance whereby the epigenetic modifications may have been passed on via the gametes that have escaped reprogramming ,,. Thus, the study of transgenerational programming of epigenetic signatures may provide a unique opportunity to identify predictive biomarkers and future therapeutic targets to promote maternal and child health.
The present findings show that prenatal stress is associated with an increased risk of shortened gestational length, poor pregnancy outcomes and delayed offspring development. Results from this study suggest that: 1) the mechanisms involved in the timing of parturition are vulnerable during early development; 2) there is a compounding effect of gestational stress on physiological and behavioural outcomes that propagate across subsequent generations; and 3) that these changes are accompanied by altered miRNA regulation in somatic cells. The identification of stress-induced epigenetic signatures in clinically accessible tissues, such as the placenta, offers an exciting potential for the prediction and prevention of PTB and poor pregnancy outcomes. The present findings concur with descriptions of inter-generational stress impacts by human migration, natural disasters and poverty, which may program maternal health preconceptionally via the maternal lineage. Although spontaneous PTB in humans is likely a multifactorial condition, the present data offer a potentially clinically relevant platform to study predictive factors and interventions for PTB and adverse developmental outcomes.
analysis of variance
quantitative real time-polymerase chain reaction
The authors are grateful to Erin Falkenberg and Rebecca Supina for assistance with the experiments, and to Amanda McRoberts for assistance with graphic illustrations. The authors would like to thank Dr. Bryan Kolb for helpful comments on the manuscript. The authors acknowledge support by the Alberta Innovates-Health Solutions Interdisciplinary Team Grant #200700595 ‘Preterm Birth and Healthy Outcomes’ (DO, GM), grants from Alberta Innovates-Health Solutions (AR, GM), the Natural Sciences and Engineering Research Council of Canada (GM), and the Canadian Institutes of Health Research #102652 (GM). GM is an AHFMR Senior Scholar. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
- Mikkola K, Ritari N, Tommiska V, Salokorpi T, Lehtonen L, Tammela O, Pääkkönen L, Olsen P, Korkman M, Fellman V: Neurodevelopmental outcome at 5 years of age of a national cohort of extremely low birth weight infants who were born in 1996-1997. Pediatrics. 2005, 116: 1391-1400. 10.1542/peds.2005-0171.View ArticlePubMed
- Kramer MS, Lydon J, Goulet L, Kahn S, Dahhou M, Platt RW, Sharma S, Meaney MJ, Séguin L: Maternal stress/distress, hormonal pathways and spontaneous preterm birth. Paediatr Perinat Epidemiol. 2013, 27: 237-246. 10.1111/ppe.12042.View ArticlePubMed
- Hobel CJ: Stress and preterm birth. Clin Obstet Gyn. 2004, 47: 856-880. 10.1097/01.grf.0000142512.38733.8c.View Article
- Zhu P, Tao F, Hao J, Sun Y, Jiang X: Prenatal life events stress: implications for preterm birth and infant birthweight. Am J Obstet Gynecol. 2010, 203: 34.e1-8. 10.1016/j.ajog.2010.02.023.View Article
- Emanuel I, Leisenring W, Wiliams MA, Kimpo C, Estee S, O’Brien W, Hale CB: The Washington state intergenerational study of birth outcomes: methodology and some comparisons of maternal birth weight and infant birthweight and gestation in four ethnic groups. Paediatr Perinat Epidemiol. 1999, 13: 352-371. 10.1046/j.1365-3016.1999.00184.x.View ArticlePubMed
- Rich-Edwards JW, Grizzard TA: Psychosocial stress and neuroendocrine mechanisms in preterm delivery. Am J Obstet Gynecol. 2005, 192: S30-S35. 10.1016/j.ajog.2005.01.072.View ArticlePubMed
- Wainstock T, Anteby E, Glasser S, Shoham-Vardi I, Lerner-Geva L: The association between prenatal maternal objective stress, perceived stress, preterm birth and low birthweight. J Matern Fetal Neonatal Med. 2013, 26: 973-977. 10.3109/14767058.2013.766696.View ArticlePubMed
- Porter F, Fraser AM, Hunter CY, Ward H, Varner MW: The risk of preterm birth across generations. Obstet Gynecol. 1997, 90: 63-67. 10.1016/S0029-7844(97)00215-9.View ArticlePubMed
- Crews D, Gillette R, Scarpino SV, Manikkam M, Savenkova MI, Skinner MK: Epigenetic transgenerational inheritance of altered stress responses. Proc Natl Acad Sci U S A. 2012, 109: 9143-9148. 10.1073/pnas.1118514109.PubMed CentralView ArticlePubMed
- Morgan CP, Bale TL: Early prenatal stress epigenetically programs dysmasculinization in second-generation offspring via the paternal lineage. J Neurosci. 2011, 31: 11748-11755. 10.1523/JNEUROSCI.1887-11.2011.PubMed CentralView ArticlePubMed
- Dias BG, Ressler KJ: Parental olfactory experience influences behavior and neural structure in subsequent generations. Nat Neurosci. 2014, 17: 89-96. 10.1038/nn.3594.PubMed CentralView ArticlePubMed
- Gapp K, Jawaid A, Sarkies P, Bohacek J, Pelczar P, Prados J, Farinelli L, Miska E, Mansuy IM: Implication of sperm RNAs in transgenerational inheritance of the effects of early trauma in mice. Nat Neurosci. 2014, 17: 667-669. 10.1038/nn.3695.PubMed CentralView ArticlePubMed
- Nilsson EE, Anway MD, Stanfield J, Skinner MK: Transgenerational epigenetic effects of the endocrine disruptor vinclozolin on pregnancies and female adult onset disease. Reproduction. 2008, 135: 713-721. 10.1530/REP-07-0542.View ArticlePubMed
- Skinner MK, Manikkam M, Tracey R, Guerrero-Bosagna C, Haque M, Nilsson EE: Ancestral dichlorodiphenyltrichloroethane (DDT) exposure promotes epigenetic transgenerational inheritance of obesity. BMC Med. 2013, 11: 228-10.1186/1741-7015-11-228.PubMed CentralView ArticlePubMed
- Veenendaal MV, Painter RC, de Rooij SR, Bossuyt PM, van der Post JA, Gluckman PD, Hanson MA, Roseboom TJ: Transgenerational effects of prenatal exposure to the 1944-45 Dutch famine. BJOG. 2013, 120: 548-553. 10.1111/1471-0528.12136.View ArticlePubMed
- Skinner MK: What is an epigenetic transgenerational phenotype? F3 or F2. Reprod Toxicol. 2008, 25: 2-6. 10.1016/j.reprotox.2007.09.001.PubMed CentralView ArticlePubMed
- Skinner MK, Manikkam M, Guerrero-Bosagna C: Epigenetic transgenerational actions of environmental factors in disease etiology. Trends Endocrinol Metab. 2010, 21: 214-222. 10.1016/j.tem.2009.12.007.PubMed CentralView ArticlePubMed
- Renthal NE, Chen CC, Williams KC, Gerard RD, Prange-Kiel J, Mendelson CR: miR-200 family and targets, ZEB1 and ZEB2, modulate uterine quiescence and contractility during pregnancy and labor. Proc Natl Acad Sci U S A. 2010, 107: 20828-20833. 10.1073/pnas.1008301107.PubMed CentralView ArticlePubMed
- Williams KC, Renthal NE, Condon JC, Gerard RD, Mendelson CR: MicroRNA-200a serves a key role in the decline of progesterone receptor function leading to term and preterm labor. Proc Natl Acad Sci U S A. 2012, 109: 7529-7534. 10.1073/pnas.1200650109.PubMed CentralView ArticlePubMed
- Metz GA, Jadavji NM, Smith LK: Modulation of motor function by stress: a novel concept of the effects of stress and corticosterone on behavior. Eur J Neurosci. 2005, 22: 1190-1200. 10.1111/j.1460-9568.2005.04285.x.View ArticlePubMed
- Ward ID, Zucchi FCR, Robbins J, Falkenberg EA, Olson DM, Benzies K, Metz GA: Transgenerational programming of maternal behaviour by prenatal stress. BMC Pregnancy Childbirth. 2013, 13: S9-10.1186/1471-2393-13-S1-S9.PubMed CentralView ArticlePubMed
- Zucchi FC, Yao Y, Ward ID, Ilnytskyy Y, Olson DM, Benzies K, Kovalchuk I, Kovalchuk O, Metz GA: Maternal stress induces epigenetic signatures of psychiatric and neurological diseases in the offspring. PLoS One. 2013, 8: e56967-10.1371/journal.pone.0056967.PubMed CentralView ArticlePubMed
- Ellenbroek BA, Derks N, Park HJ: Early maternal deprivation retards neurodevelopment in Wistar rats. Stress. 2005, 8: 247-257. 10.1080/10253890500404634.View ArticlePubMed
- Gao X, Gulari E, Zhou X: In situ synthesis of oligonucleotide microarrays. Biopolymers. 2004, 73: 579-596. 10.1002/bip.20005.View ArticlePubMed
- Zhu Q, Hong A, Sheng N, Zhang X, Jun KY, Srivannavit O, Gulari E, Gao X, Zhou X: Microfluidic biochip for nucleic acid and protein analysis. Methods Mol Biol. 2007, 382: 287-312. 10.1007/978-1-59745-304-2_19.View ArticlePubMed
- Bolstad BM, Irizarry RA, Astrandand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.View ArticlePubMed
- Wan G, Lim QE, Too HP: High-performance quantification of mature microRNAs by real-time RT-PCR using deoxyuridine-incorporated oligonucleotides and hemi-nested primers. RNA. 2010, 16: 1436-1445. 10.1261/rna.2001610.PubMed CentralView ArticlePubMed
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.PubMed CentralView ArticlePubMed
- Betel D, Koppal A, Agius P, Sander C, Leslie C: Comprehensive modeling of microRNA targets predicts functional non-conserved and non-canonical sites. Genome Biol. 2010, 11: R90-10.1186/gb-2010-11-8-r90. http://www.microrna.org/, [http://www.microrna.org/]PubMed CentralView ArticlePubMed
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009, 4: 44-57. 10.1038/nprot.2008.211. http://david.abcc.ncifcrf.gov/, [http://david.abcc.ncifcrf.gov/]View Article
- Mi H, Lazareva-Ulitsky B, Loo R, Kejariwal A, Vandergriff J, Rabkin S, Guo N, Muruganujan A, Doremieux O, Campbell MJ, Kitano H, Thomas PD: The PANTHER database of protein families, subfamilies, functions and pathways. Nucl Acids Res. 2005, 33: D284-D288. 10.1093/nar/gki078. http://www.pantherdb.org/, [http://www.pantherdb.org/]PubMed CentralView ArticlePubMed
- Mayor-Lynn K, Toloubeydokhti T, Cruz AC, Chegini N: Expression profile of microRNAs and mRNAs in human placentas from pregnancies complicated by preeclampsia and preterm labor. Reprod Sci. 2011, 18: 46-56. 10.1177/1933719110374115.PubMed CentralView ArticlePubMed
- Cottrell EC, Seckl JR: Prenatal stress, glucocorticoids and the programming of adult disease. Front Behav Neurosci. 2009, 3: 1-9. 10.3389/neuro.08.019.2009.View Article
- Schwab M: Intrauterine programming of disorders of brain function in later life. Gynakol Geburtshilfliche Rundsch. 2009, 49: 13-28. 10.1159/000184442. [in German]View ArticlePubMed
- Räikkönen K, Pesonen AK, Roseboom TJ, Eriksson JG: Early determinants of mental health. Best Pract Res Clin Endocrinol Metab. 2012, 26: 599-611. 10.1016/j.beem.2012.03.001.View ArticlePubMed
- Miller AE, Riegle GD: Progesterone and luteinizing hormone secretion following stress-induced interruption of constant estrus in aged rats. J Gerontol. 1985, 40: 129-132. 10.1093/geronj/40.2.129.View ArticlePubMed
- Arck PC: Stress and pregnancy loss: role of immune mediators, hormones and neurotransmitters. Am J Reprod Immunol. 2001, 46: 117-123. 10.1111/j.8755-8920.2001.460201.x.View ArticlePubMed
- Berends LM, Ozanne SE: Early determinants of type-2 diabetes. Best Pract Res Clin Endocrinol Metab. 2012, 26: 569-580. 10.1016/j.beem.2012.03.002.View ArticlePubMed
- Glynn LM, Wadhwa PD, Dunkel-Schetter C, Chicz-DeMet A, Sandman CA: When stress happens matters: effects of earthquake timing on stress responsivity in pregnancy. Obstetrics. 2000, 184: 637-642. 10.1067/mob.2001.111066.View Article
- Van den Hove DL, Blanco CE, Scheepens A, Desbonnet L, Myint AM, Leonard BE, Prickaerts J, Steinbusch HW: Prenatal maternal paroxetine treatment and neonatal mortality in the rat: a preliminary study. Neonatology. 2008, 93: 52-55. 10.1159/000106433.View ArticlePubMed
- Olivier JD, Akerud H, Kaihola H, Pawluski JL, Skalkidou A, Högberg U, Sundström-Poromaa I: The effects of maternal depression and maternal selective serotonin reuptake inhibitor exposure on offspring. Front Cell Neurosci. 2013, 7: 73-10.3389/fncel.2013.00073.PubMed CentralView ArticlePubMed
- Xu Z, Zhao J, Zhang H, Ke T, Xu P, Cai W, Katirai F, Ye D, Huang Y, Huang B: Spontaneous miscarriages are explained by the stress/glucocorticoid/lipoxin A4 axis. J Immuno. 2013, 190: 6051-6058. 10.4049/jimmunol.1202807.View Article
- Cheek TG, Baird E: Anesthesia for nonobstetric surgery: maternal and fetal considerations. Clin Obstet Gynecol. 2009, 52: 535-545. 10.1097/GRF.0b013e3181c11f60.View ArticlePubMed
- Uchida S, Hara K, Kobayashi A, Otsuki K, Hobara T, Yamagata H, Watanabe Y: Maternal and genetic factors in stress-resilient and -vulnerable rats: a cross-fostering study. Brain Res. 2010, 1316: 43-50. 10.1016/j.brainres.2009.11.070.View ArticlePubMed
- Weaver IC, Champagne FA, Brown SE, Dymov S, Sharma S, Meaney MJ, Szyf M: Reversal of maternal programming of stress responses in adult offspring through methyl supplementation: altering epigenetic marking later in life. J Neurosci. 2005, 25: 11045-11054. 10.1523/JNEUROSCI.3652-05.2005.View ArticlePubMed
- Zhang TY, Hellstrom IC, Bagot RC, Wen X, Diorio J, Meaney MJ: Maternal care and DNA methylation of a glutamic acid decarboxylase 1 promoter in rat hippocampus. J Neurosci. 2010, 30: 13130-13137. 10.1523/JNEUROSCI.1039-10.2010.View ArticlePubMed
- Champagne FA, Meaney MJ: Transgenerational effects of social environment on variations in maternal care and behavioral response to novelty. Behav Neurosci. 2007, 121: 1353-1363. 10.1037/0735-7044.121.6.1353.View ArticlePubMed
- Babenko O, Golubov A, Ilnytskyy Y, Kovalchuk I, Metz GA: Genomic and epigenomic responses to chronic stress involve miRNA-mediated programming. PLoS One. 2012, 7: e29441-10.1371/journal.pone.0029441.PubMed CentralView ArticlePubMed
- Lin ST, Fu YH: miR-23 regulation of lamin B1 is crucial for oligodendrocyte development and myelination. Dis Model Mech. 2009, 2: 178-188. 10.1242/dmm.001065.PubMed CentralView ArticlePubMed
- Jin J, Cheng Y, Zhang Y, Wood W, Peng Q, Hutchison E, Mattson MP, Becker KG, Duan W: Interrogation of brain miRNA and mRNA expression profiles reveals a molecular regulatory network that is perturbed by mutant huntingtin. J Neurochem. 2012, 123: 477-490. 10.1111/j.1471-4159.2012.07925.x.PubMed CentralView ArticlePubMed
- Wagner KV, Hartmann J, Mangold K, Wang XD, Labermaier C, Liebl C, Wolf M, Gassen NC, Holsboer F, Rein T, Müller MB, Schmidt MV: Homer1 mediates acute stress-induced cognitive deficits in the dorsal hippocampus. J Neurosci. 2013, 33: 3857-3864. 10.1523/JNEUROSCI.4333-12.2013.View ArticlePubMed
- Griggs EM, Young EJ, Rumbaugh G, Miller CA: MicroRNA-182 regulates amygdala-dependent memory formation. J Neurosci. 2013, 33: 1734-1740. 10.1523/JNEUROSCI.2873-12.2013.PubMed CentralView ArticlePubMed
- Wu Y, Xiao Y, Ding X, Zhuo Y, Ren P, Zhou C, Zhou J: A miR-200b/200c/429-binding site polymorphism in the 3′ untranslated region of the AP-2α gene is associated with cisplatin resistance. PLoS One. 2011, 6: e29043-10.1371/journal.pone.0029043.PubMed CentralView ArticlePubMed
- Haramati S, Navon I, Issler O, Ezra-Nevo G, Gil S, Zwang R, Hornstein E, Chen A: microRNA as repressors of stress-induced anxiety: the case of amygdalar miR-34. J Neurosci. 2011, 31: 14191-14203. 10.1523/JNEUROSCI.1673-11.2011.View ArticlePubMed
- Babenko O, Kovalchuk I, Metz GA: Environmental influences on epigenetic programming of neurodegenerative diseases. Brain Res. 2012, 1444: 96-111. 10.1016/j.brainres.2012.01.038.View ArticlePubMed
- Christiaens I, Zaragoza DB, Guilbert L, Robertson SA, Mitchell BF, Olson DM: Inflammatory processes in preterm and term parturition. J Reprod Immunol. 2008, 79: 50-57. 10.1016/j.jri.2008.04.002.View ArticlePubMed
- Kossintseva I, Wong S, Johnstone E, Guilbert L, Olson DM, Mitchell BF: Proinflammatory cytokines inhibit human placental 11beta-hydroxysteroid dehydrogenase type 2 activity through Ca2+ and cAMP pathways. Am J Physiol Endocrinol Metab. 2006, 290: E282-E288. 10.1152/ajpendo.00328.2005.View ArticlePubMed
- Nothnick WB, Healy C: Estrogen induces distinct patterns of microRNA expression within the mouse uterus. Reprod Sci. 2010, 17: 987-994. 10.1177/1933719110377472.PubMed CentralView ArticlePubMed
- Hou L, Lu Y, Li Y, Li L: MicroRNA-451 (miRNA-451) is a potential biomarker for estrogenicity in mouse uterus. Front Envir Sci Engineer. 2014, 1: 99-105. 10.1007/s11783-013-0490-7.View Article
- Heard E, Martienssen RA: Transgenerational epigeneticinheritance: myths and mechanisms. Cell. 2014, 157: 95-109. 10.1016/j.cell.2014.02.045.PubMed CentralView ArticlePubMed
- Migikovsky Z, Kovalchuk I: Epigenetic memory in mammals. Front Genet. 2011, 2: 28.
- Zucchi FC, Yao Y, Metz GA: The secret language of destiny: stress imprinting and transgenerational origins of disease. Front Genet. 2012, 3: 96-10.3389/fgene.2012.00096.PubMed CentralView ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public DomainDedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article,unless otherwise stated.