Ancestral exposure to stress epigenetically programs preterm birth risk and adverse maternal and newborn outcomes
BMC Medicine volume 12, Article number: 121 (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.
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.
Four successive generations of timed-pregnant female rats (n = 56) were bred under standardized conditions and split by treatment in each generation (see Figure 1A). Parental female rats (F0) were stressed during late gestation (S; n = 10). Their pregnant F1 daughters were split into either stressed (SS; n = 7) or non-stressed groups (SN; n = 5). Their pregnant F2 granddaughters were again either stressed (SSS to represent cumulative effects of stress; n = 9) or not stressed (SNN; n = 10, SSN; n = 7). Yoked controls were bred with each generation (N; n = 8). Data referring to offspring is provided in reference to postnatal days, data referring to mothers is provided in reference to GD and postpartum lactational days (LD).
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
In the absence of changes in litter size, stress during gestation reduced gestational weight gain in pregnant dams. Weight gain was analysed as a percentage of change compared to pre-gestational body weight. On GD11, F2-SSN and F2-SSS dams weighed on average 20 g less than F2-NNN dams (Ps <0.05). On GD21, the F0-S dams were lighter than F0-N dams (P <0.01) and F1-SS lighter than F1-SN and F1-NN (P <0.05), In the F2 generation, F2-SSS dams were lighter than F2-NNN (P <0.001) and F2-SNN dams (P <0.05; Figure 2A). There was no difference in litter size between groups, however (Figure 2B). There were no differences in uterine implantation sites between the groups in each generation; however, F2-SSN dams showed significantly more sites than F0-S and F1-SS dams (all Ps <0.05) in the absence of litter size differences. Notably, dams that showed higher blood glucose levels on LD1 also showed lower gestational weight gain (r = 0.36, P <0.05; Figure 2C) and elevated CORT levels on GD21 (r = 0.44, P <0.05; Figure 2C).
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
At one hour following delivery, the profile of motor activities in dams was used as an indicator of gestational and prenatal stress. There were no differences in the parental generation between non-stressed controls and stressed dams. Overall, the experience of prenatal stress reduced tail chasing across groups (P <0.01), which was further reduced in the F2 generation (P <0.001). In the F1 generation, tail chasing and rotational behaviours (Figure 3A) of F1-SN dams were reduced compared to F0-N controls (P <0.05, Figure 3B). In the F2 generation, SNN dams spent significantly less time in tail chasing compared to controls (P <0.01).
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
In consideration of the possible role of miRNA-mediated stress adaptation, we profiled miRNA from the frontal cortices of F0-N, F0-S and F2-SSS animals using a microarray-based approach. Compared to F0-N controls, rno-miR-138-1-3p* was significantly induced in the frontal cortex of F0-S dams whereas rno-miR-323-5p was significantly suppressed (P <0.01, n = 3). In addition, compared to F0-N rats, stress in F0-S dams induced one miRNA (rno-miR-466b-1-3p) and suppressed the expression of three miRNAs (rno-miR-145-3p, rno-miR-24-1-5p and rno-miR-375) (all Ps <0.10). Interestingly, ten miRNAs showed significant changes (P <0.10) between non-stressed F0-N and F2-SSS brains. However, since the signal level of these miRNAs was relatively low, we chose an additional set of miRNAs, which demonstrated high signal level and significantly altered expression levels based on t-test comparisons for qRT-PCR validation (Figure 4A). The qRT-PCR confirmed changes of the selected miRNAs (Figure 4B), decreased expression of miR-96, miR-141, miR-182, miR-183, miR-200a, miR-200b, miR-429 and miR-451 in F2-SSS compared to F0-S animals, whereas miR-23b and miR-200c showed increased expression levels. Thus, multigenerational stress in the F2-SSS cortex modulated miRNA profiles.
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
Cumulative multigenerational stress upregulated miR-200b and downregulated miR-429 expression levels in the uterus of F1-SS and F2-SSS generations (Figures 5A, B). Both miR-200b and miR-429 are known to modulate gestational length through interaction with their target genes Stat5b, Zeb1 and Zeb2. When upregulated, miR-200b may act to suppress Stat5b, Zeb1 and Zeb2 mRNA levels in the F1-SS and F2-SSS generations (Figures 5D-F), while reduced Zeb2 expression in particular was transmitted to the F2-SSS generation (Figure 5F). The findings suggest that miR-429 may not have a suppressive role on Stat5b, Zeb1 and Zeb2 in postpartum dams.
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
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.
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.
Hobel CJ: Stress and preterm birth. Clin Obstet Gyn. 2004, 47: 856-880. 10.1097/01.grf.0000142512.38733.8c.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Skinner MK: What is an epigenetic transgenerational phenotype? F3 or F2. Reprod Toxicol. 2008, 25: 2-6. 10.1016/j.reprotox.2007.09.001.
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.
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.
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.
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.
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.
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.
Ellenbroek BA, Derks N, Park HJ: Early maternal deprivation retards neurodevelopment in Wistar rats. Stress. 2005, 8: 247-257. 10.1080/10253890500404634.
Gao X, Gulari E, Zhou X: In situ synthesis of oligonucleotide microarrays. Biopolymers. 2004, 73: 579-596. 10.1002/bip.20005.
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.
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.
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.
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.
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/]
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/]
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/]
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.
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.
Schwab M: Intrauterine programming of disorders of brain function in later life. Gynakol Geburtshilfliche Rundsch. 2009, 49: 13-28. 10.1159/000184442. [in German]
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.
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.
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.
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.
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.
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.
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.
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.
Cheek TG, Baird E: Anesthesia for nonobstetric surgery: maternal and fetal considerations. Clin Obstet Gynecol. 2009, 52: 535-545. 10.1097/GRF.0b013e3181c11f60.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Nothnick WB, Healy C: Estrogen induces distinct patterns of microRNA expression within the mouse uterus. Reprod Sci. 2010, 17: 987-994. 10.1177/1933719110377472.
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.
Heard E, Martienssen RA: Transgenerational epigeneticinheritance: myths and mechanisms. Cell. 2014, 157: 95-109. 10.1016/j.cell.2014.02.045.
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.
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.
The authors declare that they have no competing interests.
YY, AMR, IK, DMO and GASM conceived and designed the study. YY, AMR, FCRZ, JCR and OB performed the experiments and acquired the data. All authors participated in the data analysis. YY, IK, DMO and GASM prepared the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1.: Percentage of genes participating in biological pathways that affect gestation. (DOCX 17 KB)
Additional file 2: Table S2.: Number of target genes involved in pathways screened by PANTHER (red font denotes pathways in neuropathology). (DOCX 30 KB)
About this article
Cite this article
Yao, Y., Robinson, A.M., Zucchi, F.C. et al. Ancestral exposure to stress epigenetically programs preterm birth risk and adverse maternal and newborn outcomes. BMC Med 12, 121 (2014). https://doi.org/10.1186/s12916-014-0121-6