Epigenetic regulation of histone modifications and Gata6 gene expression induced by maternal diet in mouse embryoid bodies in a model of developmental programming

Dietary interventions during pregnancy alter offspring fitness. We have shown mouse maternal low protein diet fed exclusively for the preimplantation period (Emb-LPD) before return to normal protein diet (NPD) for the rest of gestation, is sufficient to cause adult offspring cardiovascular and metabolic disease. Moreover, Emb-LPD blastocysts sense altered nutrition within the uterus and activate compensatory cellular responses including stimulated endocytosis within extra-embryonic trophectoderm and primitive endoderm (PE) lineages to protect fetal growth rate. However, these responses associate with later disease. Here, we investigate epigenetic mechanisms underlying nutritional programming of PE that may contribute to its altered phenotype, stabilised during subsequent development. We use embryonic stem (ES) cell lines established previously from Emb-LPD and NPD blastocysts that were differentiated into embryoid bodies (EBs) with outer PE-like layer. Emb-LPD EBs grow to a larger size than NPD EBs and express reduced Gata6 transcription factor (regulator of PE differentiation) at mRNA and protein levels, similar to Emb-LPD PE derivative visceral yolk sac tissue in vivo in later gestation. We analysed histone modifications at the Gata6 promoter in Emb-LPD EBs using chromatin immunoprecipitation assay. We found significant reduction in histone H3 and H4 acetylation and RNA polymerase II binding compared with NPD EBs, all markers of reduced transcription. Other histone modifications, H3K4Me2, H3K9Me3 and H3K27Me3, were unaltered. A similar but generally non-significant histone modification pattern was found on the Gata4 promoter. Consistent with these changes, histone deacetylase Hdac-1, but not Hdac-3, gene expression was upregulated in Emb-LPD EBs. First, these data demonstrate ES cells and EBs retain and propagate nutritional programming adaptations in vitro, suitable for molecular analysis of mechanisms, reducing animal use. Second, they reveal maternal diet induces persistent changes in histone modifications to regulate Gata6 expression and PE growth and differentiation that may affect lifetime health.


Background
Periconceptional environment, especially during oocyte maturation and preimplantation development, can influence the pattern of later gestation leading to permanent changes in offspring growth, physiology, health and disease risk through to adulthood [1][2][3]. Factors such as the quality and composition of maternal or paternal diet, parental metabolism and health, or specific conditions as used in assisted conception such as embryo culture, can all influence the developmental programme. This sensitive window in the lifecycle around conception can be viewed within the broader context of the Developmental Origins of Health and Disease (DOHaD) concept. This proposes that risk of adult onset diseases may derive from in utero conditions where nutrient availability may control fetal growth and metabolic homeostasis but which may predispose to adult disease, particularly cardiovascular dysfunction and metabolic syndrome, if homeostatic changes do not match postnatal environment. Epidemiological studies on human populations and various animal models show support for the DOHaD concept [4][5][6][7].
We have used a rodent maternal low protein diet model to study mechanisms of periconceptional programming whereby protein restriction is applied exclusively during the period from mating to blastocyst formation (Emb-LPD, 9% casein, E0-3.5 in mouse) with normal nutrition (NPD, 18% casein) provided for the remainder of gestation, and standard chow diet postnatally. This brief nutritional challenge is sufficient to induce cardiometabolic dysfunction, hypertension and abnormal behaviour in adulthood [8,9]. Emb-LPD changes the pattern of development by altering the composition of the uterine fluid which is detected by blastocysts via mTOR signalling [10]. The embryo responds to the nutrient challenge by activating several compensatory processes within extraembryonic lineages which collectively act to increase nutrient provision from the mother for the remainder of gestation to protect fetal growth. These responses include increased endocytosis and proliferation within the trophectoderm lineage (TE; progenitor of chorioallantoic placenta) and increased motility and invasiveness of outgrowing trophoblast at the time of implantation [10,11]. LPD treatment maintained throughout gestation leads to increased nutrient transport efficiency in the midand late-gestation placenta [12]. Stimulated endocytosis is also seen in response to Emb-LPD in the primitive endoderm (PE) extra-embryonic lineage formed from the blastocyst inner cell mass (ICM); this response is maintained until late gestation within the derivative visceral endoderm of the yolk sac placenta to promote nutrient uptake from the uterine milieu [9,11]. Nutrient provision and growth promotion resulting from these extraembryonic adaptations to poor maternal diet, whilst likely favouring competitive fitness of offspring during periods of limited food supply, also lead to later chronic disease when the diet improves, evidenced by the resulting perinatal weight correlating with adult CV and behavioural dysfunction [9].
Since extra-embryonic responses to Emb-LPD persist from early development throughout gestation and have important consequences for protecting conceptus growth and affecting adult disease risk, we anticipate conserved epigenetic mechanisms may be driving these physiological processes. Moreover, the compensatory changes persist well beyond the period of dietary challenge and reflect a 'memory' of an earlier environment. Periconceptional induction of epigenetic change has been demonstrated in other models of programming, such as following in vitro culture treatment of pre-implantation embryos [13][14][15][16][17]. However, clear evidence of epigenetic modifications driving physiological responses within an in vivo periconceptional programming model has not been forthcoming previously.
Here, we investigate the epigenetic status of histone modifications occurring within the PE lineage in response to Emb-LPD for evidence of the programming of altered phenotype. We have used embryoid bodies (EBs) derived from embryonic stem (ES) cell lines generated from Emb-LPD and NPD blastocysts since the PE-like layer on the surface of Emb-LPD EBs exhibit the enhanced endocytosis compensatory phenotype after at least six passages in standard culture [11].

Results
Effect of maternal diet on EB size ES cells derived from blastocysts collected from Emb-LPD and NPD females were maintained from passage 6 for 5.5 days in culture for EB formation in 96-well low adhesion plate culture. At this time point, EBs have formed primitive endoderm-like (PE) layer on their surface as demonstrated by the presence of Gata6 and Dab2 marker proteins [11]. EBs were imaged and diameters measured. EBs differentiated from Emb-LPD ES cells were significantly larger (~15%; p < 0.05) than NPD EBs ( Figure 1).

Effect of maternal diet on embryoid body Gata factor expression
The PE lineage is regulated through Gata6 and Gata4 transcription factors that activate PE specification and differentiation through several downstream target genes [18][19][20][21]. We investigated the expression of Gata4 and Gata6 and the downstream target gene, Dab2 in Emb-LPD and NPD EBs. Gata6 gene expression was significantly reduced in Emb-LPD EBs while Gata4 was reduced but only to trend level and Dab2 expression was unaffected by maternal dietary origin (Figure 2A). Reduced protein expression of Gata6 but not Dab2 was also evident in Emb-LPD EBs ( Figure 2B). The authenticity of the changed Gata6 expression in Emb-LPD EBs was supported by ex vivo analysis of mouse E17 visceral yolk sac tissue also showing reduced Gata6 protein expression compared with NPD control ( Figure 2C).

Effect of maternal diet on Gata6 promoter histone acetylation in embryoid bodies
Histone modifications have been shown to regulate Gata4 and Gata6 gene expression in other models [22,23]. We used ChIP assay to examine three upstream regions of the Gata4 and Gata6 genes ( Figure 3A) in EBs to compare levels of histone modifications within promoter domains with respect to maternal diet. A panel of antibodies was used to probe histone acetylation and methylation with targets for acetylated H3 and H4, H3K4Me3, H3K4Me2, H3K9Me3 and H3K27Me3. Our analysis revealed a distinct pattern of modifications dependent upon dietary origin of EBs. The Gata6 promoter at G6P1 site in Emb-LPD EBs exhibited significant hypoacetylation of both H3 and H4 and a trend of decrease in the density of the histone marker H3K4Me3 compared with NPD EBs ( Figure 3B). In addition to these three histone modification changes known to be associated with reduced gene expression [24], the G6P1 Gata6 promoter site in Emb-LPD EBs had reduced enrichment of RNA polymerase II ( Figure 3B), further supporting a suppressed state of expression. No significant changes in other histone modifications, including H3K4Me2, H3K9Me3 and H3K27Me3, were detected. The Gata4 promoter site at G4P1 showed a pattern of histone modifications in Emb-LPD EBs similar to those in the Gata6 promoter but not to the level of significance except that H3K9Me3 was significantly reduced compared with NPD EBs ( Figure 3C). At other sites upstream of the Gata6 and Gata4 promoter domains (G6P3, G6P5, G4P3, G4P5), maternal diet had little or no effects on histone modifications ( Figure 4A-D); however, H3K4Me3 on G4P3 showed a significant increase in Emb-LPD EBs ( Figure 4C). We used Gapdh as a control house-keeping gene in our ChIP analysis and found no difference in histone modifications on its promoter in EBs with respect to maternal diet ( Figure 4E).

Effect of maternal diet on expression of histone deacetylases in embryoid bodies
Given the pattern of histone hypoacetylation detected at the Gata6 promoter coinciding with reduced expression of this gene in Emb-LPD EBs, we evaluated whether the expression of histone deacetylases (HDACs) was altered in response to maternal diet. HDACs are expressed in early embryos and have been shown to modify their expression in response to in vitro culture [25]. We found Hdac-1 gene expression but not Hdac-3 was upregulated in Emb-LPD EBs ( Figure 5).

Discussion
In this study we have investigated the molecular mechanisms of developmental programming of the PE extraembryonic lineage within our Emb-LPD mouse model associated with adult-onset disease. The Emb-LPD PE shows enhanced endocytosis at ligand, lysosomal and receptor levels, a cellular modification that is sustained through to late gestation in the yolk sac visceral endoderm to support increased nutrient uptake despite poor maternal nutrition, thereby likely to protect fetal development [9,11]. Significantly, these changes are maintained after induction even if maternal diet is returned to control levels, a characteristic we have demonstrated also occurs in the trophectoderm (TE) extra-embryonic lineage [10,11] and is suggestive of epigenetic mechanisms, the focus of the current study. We chose to investigate potential PE epigenetic mechanisms using EBs since they maintain the enhanced endocytosis cellular phenotype induced by nutritional programming previously revealed from ex vivo tissues [9,11] and, in the current study, show similar downregulation of Gata6 as seen in ex vivo VYS, thereby confirming their authenticity. Also, the PE-like layer formed on EBs from undifferentiated ES cells is representative of the PE layer formed in the late blastocyst from the ICM, and PE-like cells in EBs as well as the core primitive ectoderm have been used extensively as a model for early cell lineage derivation in post-implantation development over many years [26,27].
Our results show significant change in gene expression and histone modifications in Emb-LPD EBs in comparison with NPD EBs. Emb-LPD EBs express reduced levels of Gata6 transcription factor, a key regulator of PE differentiation [18,19], consistent with Gata6 downregulation found in vivo within the Emb-LPD VYS. Reduced Gata6 expression coincided with histone H3 and H4 hypoacetylation and reduced recruitment of RNA polymerase II at the Gata6 promoter in Emb-LPD EBs, all factors consistent with reduced gene expression [24]. Moreover, the expression of Hdac-1 but not Hdac-3 was increased in Emb-LPD EBs. These HDACs are known to be expressed in early embryos and are sensitive to culture conditions so represent good candidates for coordinating histone epigenetic programming [25]. Collectively, these data reveal for the first time a persistent change in EB epigenetic status mediated through maternal diet and provide new clues to the origin of developmental programming mechanisms in this model. They also demonstrate that ES cells and derivative EBs maintain altered gene expression in vitro, long after the inductive dietary challenge, and provide a useful tool for analysis of underlying mechanisms, reducing the requirement for experimental animals.
The finding that Emb-LPD EBs, and ex vivo VYS, exhibit reduced Gata6 expression, apparently regulated through histone hypoacetylation in the EBs, may explain the increased size of Emb-LPD EBs. Gata6 and Gata4 are zinc-finger transcription factors that perform multiple roles both during development in the determination of cell lineages and in adult tissues in maintaining cell differentiation states [28,29]. Loss of Gata factor expression has been implicated in several forms of cancer and in ovarian cancer models where loss of Gata6 and Gata4 expression coincides with Gata6 and Gata4 promoter histone hypoacetylation in response to HDAC activity coupled with growth promotion and malignancy [23,30,31] indicating a similarity in epigenetic and cellular mechanisms to the current study.
Whilst Emb-LPD EBs had downregulated Gata6 expression, Gata4 expression was not significantly affected. In ES cells, absence of Gata6 gene leads to loss of Gata4 expression whilst absence of Gata4 gene does not inhibit Gata6 expression, indicating a hierarchical relationship [21,32,33]. However, functional redundancy exists between Gata6 and Gata4 expression in other models including pancreatic and ovarian germ cell differentiation and in EBs during myocyte differentiation with Gata4 expression not dependent upon Gata6 expression [34][35][36]. Thus, the distinction between Gata6 and Gata4 expression in the current study may be explained either by functional redundancy or by Gata6 expression, although reduced, being above the threshold required for Gata4 expression.
The Emb-LPD EBs also showed expression of the Gata6 downstream target gene, Dab2, required for epithelial function especially in receptor-mediated endocytosis [37]. Dab2 acts as a cargo-selective adaptor protein facilitating apical localisation of the megalin receptor (Lrp2 gene) and clathrin-mediated endocytosis [37]. The endocytic function in Emb-LPD EBs is stimulated together with increased expression of megalin as a compensatory response to maternal diet [11] and Dab2 expression and function is likely protected in the EB model to achieve this. Dab2 expression may be maintained either by a Gata4-dependent pathway [33,38], by Gata6 being above the threshold required for Dab2 promoter activation or by an alternative mechanism. We anticipate that the growth stimulation coinciding with reduced expression of Gata6 but not Dab2 in Emb-LPD EBs will be partly driven by the increased nutrient delivery provided by enhanced endocytosis. The resulting increase in growth in Emb-LPD EBs may extend to both surface PE-like layer and core epiblast cells and, like other compensatory response mechanisms, may reflect in vivo processes that safeguard fetal development against dietary deficiency. Indeed, the increased endocytic activity observed within the mature VYS in late gestation may be dependent upon this initial growth stimulation.

Conclusions
We have shown that maternal diet regulates the epigenetic status of the early embryo, reducing expression of the PE lineage regulator, Gata6 transcription factor, in EBs. Reduced Gata6 gene expression coincides with histone hypoacetylation and loss of RNA polymerase binding of the Gata6 promoter, and stimulation of Hdac-1 expression. These changes are associated with increased growth of the Emb-LPD EB which may contribute, alongside stimulation in endocytosis, as compensatory responses to support maintenance of nutrient delivery.

Ethics statement
All animal research was conducted under UK Home Office project license and local ethics approval (University of Southampton).

Animals, diet treatment and embryo collection
MF1 mice were bred in-house (University of Southampton Biomedical Research Facility) on a 0700-1900 light cycle with standard chow, under UK Home Office license and local ethics approval. Virgin females (7-8.5 weeks) were mated naturally overnight with MF1 males and plug positive females were housed individually the following morning and assigned randomly to either normal protein diet (18% casein, NPD) or isocaloric low protein diet (9% casein, Emb-LPD) until embryonic day 3.5 (E3.5); diet composition has been described elsewhere [8,9]. Embryos were collected at the blastocyst stage after cervical dislocation and uterine flushing with H6 medium with 4 mg/ml BSA (H6 + BSA) [39].

ES cell culture and embryoid body (EB) formation
Mouse embryonic stem (ES) cell lines were prepared using standard procedures from blastocysts derived from mothers fed NPD or Emb-LPD upon culture in knockout-Dulbecco's modified Eagle medium [high glucose] (DMEM [high glucose], Gibco). including 20% knock out serum replacement (Gibco) with other supplements and were subsequently maintained and expanded on mouse embryonic fibroblast feeder layers up to passage 5-7 as described in detail elsewhere [11]. A total of 18 ES lines were generated from Emb-LPD blastocysts and 38 lines from NPD blastocysts. Six clones, each from a different mother and of male gender and normal karyotype, from each diet treatment were selected for use in the current study and were used for all relevant experiments. For embryoid body (EB) formation, ES cells were dissociated with 0.05% trypsin-EDTA and suspended in ES cell culture medium without leukaemia inhibitory factor (LIF) supplementation for 1 h on gelatin-treated dishes. A cell suspension (4,000 in 200 μl) was subsequently pipetted into low-adherence 96-well plates and statically incubated at 37°C in humidified air with 5% CO 2 for 5.5 days to form individual EBs within each well using a method previously described that is optimised for uniform EB size production [11,40]. EB diameter at specific time intervals was measured using an Olympus microscope software and Cell sense®.

RNA isolation and real-time PCR
RNA isolation and quantitative real-time PCR (qRT-PCR) of EBs and tissues was performed as described previously [41]. Briefly, total RNA was extracted from EBs

Hdac1
Hdac3 Relative mRNA abundance NPD Emb-LPD * Figure 5 Gene expression of histone modification enzymic regulators Hdac-1 and Hdac-3 in embryoid bodies (EBs) at day 5.5 culture in relation to maternal diet. Expression presented as ratio to the geometric mean of Gapdh and Ppib transcripts. (n = 5 per treatment). *p < 0.05. and tissues using the RNeasy Mini kit (Qiagen, UK), with on-column DNase digestion. RNA was quantified using the Nanodrop ND-1000 spectrophotometer, and cDNA generated using a random priming strategy and the ImProm-II™ Reverse Transcription System (Promega, UK). cDNA was diluted to a concentration equivalent to 5 ng RNA per μl and used at 1 μl in a reaction volume of 20 μl with forward and reverse primers at final concentration of 300 nM each in qRT-PCR using the Chromo4 Real-Time Detector (BioRad, UK) with Opticon Monitor v3.1 software. Thermal cycling conditions were 95°C 10 min enzyme activation, then 40 cycles of 95°C for 15 s followed by 60°C for 1 min, with a final extension step of 10 min at 72°C. Primers used for qRT-PCR were designed by Primer3 software (Table 1). For EBs, Gapdh and Ppib were selected from 6 house-keeping candidates with GeNorm and NormFinder software for stability, showing no change in expression between Emb-LPD and NPD treatments ( Table 2). For quantification, efficiency of primers was determined by series 1/10 dilution and Ct value. Efficiency was calculated as E = 10 1/slope and qualified primer efficiency is between 1.9-2.1. Calculation of relative expression of target gene was calculated as E^-dCt then divided by the geometric mean of relative expression of the reference gene pair.

Protein isolation and western blotting
Approximately 100 EBs on day 5.5 were washed with PBS and lysed with 120 μl RIPA buffer (50 mM Tris-HCl [pH 7.5], 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 1% NP-40, 1% sodium deoxycholate, 1 mM Na 3 VO 4 , 50 mM NaF, cOmplete EDTA-free protease inhibitor cocktail (Roche), 0.5 mM PMSF). Lysate was sonicated on ice and protein content detected using the BCA assay (Pierce). 20 μl protein was mixed with 4x LDS sample buffer and DTT and boiled for 5 min before electrophoresis. After electrophoresis, protein was transferred to polyvinylidene fluoride membrane (0.45 μm) using wet transfer protocol. This was followed by blocking with 5% milk in TBS and incubation with primary antibody overnight at 4°C. The next day the membrane was washed with TBST and incubated with IRDye secondary antibody (Odyssey; 1:10,000) for 1 hour, washed and imaged with the Licor western detection system. Band intensity was quantified with Licor software. Antibodies used are shown in Table 3. α-tubulin was used as loading control and did not change in expression with respect to dietary treatment (data not shown).

Chromatin immunoprecipitation (ChIP)
A modified protocol was used to perform the Matrix ChIP assay [42]. EBs were dissociated with trypsin EDTA for 5 mins at 37°C, then serum was added to stop trypsinisation and the sample passed through a 21 g needle to   RT. Wells were then incubated with 0.25 μg antibody in 50 μl blocking buffer for 2 hours, RT. Chromatin (equivalent of 5x10 4 cells per assay) was diluted to 50 μl with blocking buffer, pre-incubated for 15 mins on ice and then incubated in the wells for 2.5 hours at 4°C. Wells were washed with 100 μl ice-cold IP buffer 7 times and twice with 100 μl ice cold TE buffer (10 mM Tris, 1 mM EDTA; pH 7.6). Finally, 50 μl elution buffer (0.1 mg/ml proteinase K, 25 mM Tris base, 1 mM EDTA, pH 10) was added per well. 1/10th of the sample chromatin in 50 μl of elution buffer was used as input.
Plates were incubated for 30 min at 55°C, followed by 10 min at 95°C in an ABI 7900HT thermocycler prior to quantitative PCR (qPCR). Antibodies used in ChIP are shown in Table 3. Each qPCR reaction used 2 μl immunoprecipitated DNA. The modified histone binding was expressed as fold enrichment ratio to IgG negative control (Cell Signaling). Standards and samples were simultaneously amplified in 10 μl reaction volume and primers were designed to amplify genomic sequences at the 5′UTR of Gata4, Gata6 and Gapdh genes (see Table 1).

Statistical analysis
Results were expressed as mean ± S.E.M. Comparison of mRNA expression, histone modification or Western blotting between control and treated groups was performed by one way ANOVA or t-test. EB diameters were compared across treatments using multilevel random effects regression model (SPSS). Significance testing was set at P < 0.05 (two-tailed).