Loss of maternal EED results in postnatal overgrowth
Clinical Epigenetics volume 10, Article number: 95 (2018)
Investigating how epigenetic information is transmitted through the mammalian germline is the key to understanding how this information impacts on health and disease susceptibility in offspring. EED is essential for regulating the repressive histone modification, histone 3 lysine 27 tri-methylation (H3K27me3) at many developmental genes.
In this study, we used oocyte-specific Zp3-Cre recombinase (Zp3Cre) to delete Eed specifically in mouse growing oocytes, permitting the study of EED function in oocytes and the impact of depleting EED in oocytes on outcomes in offspring. As EED deletion occurred only in growing oocytes and females were mated to normal wild type males, this model allowed the study of oocyte programming without confounding factors such as altered in utero environment. Loss of EED from growing oocytes resulted in a significant overgrowth phenotype that persisted into adult life. Significantly, this involved increased adiposity (total fat) and bone mineral density in offspring. Similar overgrowth occurs in humans with Cohen-Gibson (OMIM 617561) and Weaver (OMIM 277590) syndromes, that result from de novo germline mutations in EED or its co-factor EZH2, respectively. Consistent with a role for EZH2 in human oocytes, we demonstrate that de novo germline mutations in EZH2 occurred in the maternal germline in some cases of Weaver syndrome. However, deletion of Ezh2 in mouse oocytes resulted in a distinct phenotype compared to that resulting from oocyte-specific deletion of Eed.
This study provides novel evidence that altering EED-dependent oocyte programming leads to compromised offspring growth and development in the next generation.
Factors regulating oocyte (egg) and sperm programming and early embryonic development have been associated with the fetal origins of disease, including reduced cognitive ability and increased chronic diseases, such as type-2 diabetes, obesity, heart disease and behavioural anomalies [1,2,3,4,5]. The causes of these defects are poorly understood, but are likely to be in part due to altered epigenetic programming of oocytes or sperm that significantly impact on embryonic development and underlie the fetal origin of some of these disorders [1, 2, 4]. Defining how inherited epigenetic information regulates fetal development and postnatal phenotypic outcomes is therefore important for understanding how inherited epigenetic information impacts on human health and disease. While the establishment of aberrant epigenetic states has been associated with disease, functional analyses of epigenetic inheritance that isolate in vivo effects generated in the germline are not possible in humans. Here, we have used genetic mouse models to mediate oocyte-specific deletion of Eed and Ezh2 to understand how changes in epigenetic information established in the oocyte can impact on phenotypic outcomes in offspring.
H3K27me3 is a critical epigenetic modification that is catalysed by polycomb repressive complex 2 (PRC2), a highly conserved epigenetic modifying complex. PRC2 is comprised of three core protein subunits: Embryonic Ectoderm Development (EED), Enhancer of Zeste 1/2 (EZH1/2) and Suppressor of Zeste 12 (SUZ12). Global loss of function to any one of these components results in drastically compromised enzymatic activity of PRC2, substantial loss of H3K27me3 and embryonic lethality in mice [6,7,8]. Recent studies of germ cells have demonstrated enrichment of H3K27me3 at key developmental genes in male and female germ cells [9,10,11] and the maintenance of some H3K27me3 marked histones in mature sperm [12,13,14]. Furthermore, oocyte-specific deletion of Ezh2 results in loss of H3K27me3 in the zygote and ~ 40% growth restriction in maternal offspring . In humans, de novo germline mutations in EED or EZH2 lead to Cohen-Gibson (OMIM 617561) or Weaver (OMIM 28229590) syndromes, which are characterised by overgrowth, skeletal defects and advanced bone age [16,17,18,19,20,21,22,23], indicating that PRC2 activity may be required in the human germline for regulating outcomes in offspring. Together, these studies raise the possibility that PRC2 and H3K27me3 may underpin epigenetic inheritance effects on offspring development and postnatal outcomes.
During female fetal development, germ cells commit to oogenesis and enter meiotic prophase . Folliculogenesis occurs in early postnatal life and results in establishment of a finite pool of quiescent primordial follicles that underpin the female reproductive lifespan. Subsequently, primordial follicles are continuously released in reproductively mature females, initiating a prolonged period of oocyte growth and maturation that takes around 21 days in mice and up to 12 months in humans . While most of this period involves oocyte growth while in a diploid state, it is completed by rapid maturation of each oocyte through meiosis I and meiosis II, ultimately producing a haploid oocyte at fertilisation . Before completion of the first meiotic division, the diploid oocyte enters the germinal vesicle (GV) stage, characterised by decondensed chromatin and a period of high transcription that includes the production of maternal factors (proteins and RNAs) that are required for directing preimplantation development in the offspring . In addition to maternal factors, the mature oocyte carries specific epigenetic information required for offspring development, but the nature, extent and effects of this information on offspring development are not yet fully understood [28, 29].
Germ cell-specific knockout mouse models have been valuable for analysing the function of maternal factors in the absence of confounding effects such as in utero environment and nutritional influences derived from the mother [15, 28, 30]. Zp3Cre, is a well-established model for generating oocyte-specific gene deletion mediated by transcription of Cre recombinase under control of the Zona Pellucida 3 (Zp3) promoter . This model allows the production of offspring from oocytes that lack specific genes only during their maturation and can be effectively used for functional analyses of maternal inheritance [15, 31,32,33].
Although deletion of the PRC2 component, Ezh2, in the growing oocyte leads to restricted growth in offspring , little is known about the role of Eed in the female germline. However, the PRC1 components, Ring finger protein 2 (Rnf2) and Ring finger protein 1 (Ring1), are regulated throughout oogenesis and absence of maternal RNF2 and RING1 proteins leads to developmental arrest at the two-cell stage of embryogenesis . Moreover, injection of the H3K27me3-specific demethylase, Kdm6b, into zygotes resulted in ectopic maternal expression of specific genes that are normally expressed only from the paternal allele, demonstrating that maternal H3K27me3 regulates DNA methylation-independent imprinting . These data suggest oocyte-specific roles for polycomb group proteins in the establishment of maternal factors and/or epigenetic modifications that are required for correct development of offspring.
In this study, we used Zp3Cre to delete Eed specifically from growing oocytes. The resulting oocytes lacked H3K27me3 and produced offspring with significant overgrowth that involved increased adiposity (percentage fat) and increased bone mineral density. These data demonstrate that EED is required in the oocyte for programming developmental outcomes that affect life-long outcomes in offspring. This contrasted with growth restriction that resulted from deletion of Ezh2 in growing oocytes, indicating that while EED mediates epigenetic inheritance through the maternal germline, the function of EZH2 is more complex.
Immunofluorescent studies of oocyte development demonstrated that H3K27me3 is increasingly enriched in the nucleus of oocytes during follicle growth . In this study, we used a Zp3Cre transgene to mediate oocyte-specific deletion of Eed, in order to determine the effect loss of maternal PRC2 function had on offspring (Fig. 1a). The Eedfl/fl model was originally developed by Orkin and colleagues, and results in deletion of exons 3–6 of Eed and loss of the ability for PRC2 to catalyse H3K27 methylation . Genotyping of progeny from Eed floxed females carrying Zp3Cre mated to wild type males demonstrated that Eed deletion was 100% efficient in oocytes (Additional file 1: Figure S1; Chi-square test, nsd). Mouse oocyte growth encompasses 22–24 days, after which they complete meiosis during maturation within a single day, ultimately producing haploid ova [37, 38]. Deletion using Zp3Cre occurs early during oocyte growth, resulting in loss of the target gene during the majority of the diploid (4C) oocyte growth phase . Critically, immunofluorescence demonstrated that Eeddel/del oocytes contained markedly reduced nuclear H3K27me3, consistent with substantial loss of PRC2 function in these cells (Fig. 1b). Initially, we determined whether deleting Eed in oocytes affected average litter size produced from Eedwt/wt Eedwt/del and Eeddel/del growing oocytes. Despite severe depletion of H3K27me3, Eeddel/del oocytes produced live born pups in normal sex ratios. However, the average litter size produced by females producing Eeddel/del oocytes was significantly reduced compared to females producing Eedwt/wt or Eedwt/del oocytes, although the causes for the litter size reduction remain unknown (Fig. 1c, n = 13–20 litters per genotype). This model therefore provides an opportunity to study offspring derived from oocytes that lack EED and H3K27me3, specifically during oocyte growth, maturation and preimplantation development, prior to activation of the paternal Eed allele.
Mating females with oocyte-specific deletion of Eed to wild type C57BL/6 males provided the opportunity to compare isogenic offspring in the absence of confounding maternal in utero effects (Fig. 2). Importantly, oocytes developing in Eedwt/fl;Zp3-Cre transgenic mothers have one intact copy of Eed during the oocyte growth period, when epigenetic modifications including DNA methylation and H3K27me3 are established. In contrast, oocytes from Eedfl/fl;Zp3-Cre transgenic females lack both copies of Eed (Eeddel/del) and have no EED function in their growing oocytes. Therefore, in this model, we expect the epigenome of Eeddel haploid oocytes produced by Eedfl/wt;Zp3-Cre females to be relatively normal as they developed in the presence of reduced, yet sufficient, EED function during their growth and maturation. However, Eeddel haploid oocytes produced by Eedfl/fl;Zp3-Cre females lacked EED protein during oocyte growth and maturation (Fig. 2).
This model allowed the production of heterozygous (HET) offspring from females producing Eeddel/del (hom) or Eedwt/del (het) growing oocytes mated to wild type males. Although HET offspring from Eeddel/del or Eedwt/del growing oocytes were isogenic, they have two critical differences. Firstly, we expect these offspring to be epigenetically different, reflecting the presence or absence of EED function in the growing oocyte. Secondly, these offspring lacked maternal EED during the zygotic to paternal genome activation phase of preimplantation development. Therefore, we proposed that differences between the HET offspring from females producing Eeddel/del or Eedwt/del growing oocytes could be attributed to altered programming in Eeddel/del growing oocytes and during the earliest stages of preimplantation development. Similar comparisons were made between wild type (WT) offspring from Eedwt/fl;Zp3-Cre and Eedwt/wt females, which respectively produced Eedwt/del and Eedwt/wt growing oocytes (Fig. 2).
To identify postnatal differences in offspring due to compromised EED-dependent programming in the oocyte and early development, we initially weighed postnatal day (PND) 2 pups generated from Eeddel/del, Eedwt/del and Eedwt/wt growing oocytes and isogenic wild type sperm (Figs. 2 and 3a). Remarkably, complete loss of EED specifically in the growing oocyte resulted in a 29% increase in weight of PND2 HET offspring from Eeddel/del growing oocytes (average weight 2.19 g, n = 51) compared to PND2 HET offspring from Eedwt/del growing oocytes (average weight 1.69 g, n = 52) (Fig. 3a P < 0.0001, one-way ANOVA). In addition, HET offspring from Eeddel/del growing oocytes had increased nose to rump length, compared to isogenic age-matched HET offspring from Eedwt/del growing oocytes (Fig. 3b–d). Similarly, HET offspring from Eeddel/del growing oocytes were significantly longer and heavier than WT offspring from Eedwt/del growing oocytes, but there was no difference in length or weight between HET and WT offspring generated from Eedwt/del or Eedwt/wt growing oocytes (Fig. 3d). There was no bias in this PND2 overgrowth phenotype between the sexes (Additional file 1: Figure S2).
One possible explanation for the increased pup weight at PND2 was the reduced average litter size observed for females producing Eeddel/del growing oocytes (Fig. 1c). To account for this, we analysed the weight data using linear mixed models to statistically evaluate the litter size dependent, and litter size independent contributions to offspring weight. To account for non-independence of pups from individual litters, we included litter identification as a random effect in the analysis. Initially, we ensured that there were no differences in offspring weight due to litter size differences that were attributable individual oocyte genotypes. This revealed that the relationship between litter size and offspring weight among WT and HET pups produced from Eedwt/wt, Eedwt/del and Eeddel/del oocytes was consistent (litter size by offspring group interaction: F3,48.4 = 0.048, P = 0.99; variance components: litter ID = 0.0423, residual = 0.0346). Next, we determined the relationship between litter size and offspring weight for all genotypes. This revealed that as litter size decreased, there was a minor increase in offspring weight (t33.5 = − 1.32, P = 0.20; variance components: litter ID = 0.0398, residual = 0.0345; Fig. 3e). However, this effect was not significant and was insufficient to account for the substantial increase in the weight of HET offspring generated from Eeddel/del growing oocytes compared to HET or WT offspring generated from Eedwt/del or Eedwt/wt oocytes. In contrast, after removing the weight change attributable to litter size, there remained a highly significant increase in weight of HET pups produced from Eeddel/del growing oocytes compared to HET pups produced from Eedwt/del oocytes (t34.8 = 3.44, P = 0.002) and WT pups produced from either Eedwt/del (t36.3 = 2.81, P = 0.008) or Eedwt/wt (t33.1 = 2.29, P = 0.03) growing oocytes (Fig. 3e). Given that this increase in offspring weight only occurred in the HET pups generated from Eeddel/del growing oocytes and was not accounted for by litter size differences, we concluded that loss of EED in the oocyte and early embryo resulted in substantial, postnatal overgrowth.
As the Eed HET pups produced were isogenic, the simplest explanation for this difference is that loss of maternal EED in the mouse oocyte and zygote led to a significant early developmental programming effect that impacted on postnatal weight in offspring from Eeddel/del compared to Eedwt/del growing oocytes. As EED is known only to regulate epigenetic outcomes, this postnatal effect is likely to result from altered epigenetic programming in the oocyte and early embryo.
Recent studies in humans have demonstrated that de novo germline mutations in EED lead to Cohen-Gibson (OMIM 617561) syndrome, which is characterised by overgrowth, skeletal defects and advanced bone age [16,17,18,19,20]. DEXA scanning of the PND2 mouse offspring produced in this study revealed increased bone mineral density (BMD) in HET offspring from Eeddel/del growing oocytes compared to HET controls generated from Eedwt/del growing oocytes (Eedwt/wt n = 28, Eedwt/del n = 14, Eeddel/del n = 19; P < 0.05, one-way ANOVA; Fig. 4a). Moreover, fat content was also significantly increased in PND2 HET offspring from Eeddel/del oocytes compared to HET offspring from Eedwt/del growing oocytes, and WT offspring from Eedwt/del and Eedwt/wt oocytes (Fig. 4b, P < 0.01, one-way ANOVA). Lean muscle content was significantly reduced in HET offspring from Eeddel/del growing oocytes, compared to WT offspring from Eedwt/wt oocytes (Fig. 4c, P < 0.05, one-way ANOVA), but there was no significant difference between HET offspring from Eedwt/del and Eeddel/del oocytes. The ponderal index (weight/crown-length3) of WT and HET offspring from Eedwt/wt and Eedwt/del growing oocytes, respectively, was not significantly different indicating that there was no substantial effect on corpulence in these offspring. However, while there was an unexpected marginal decrease in ponderal index in WT pups from Eedwt/del growing oocytes compared to all other treatment groups, the significance of this change remains unknown (Fig. 4d, *P < 0.05, **P < 0.01 one-way ANOVA).
To determine the weight phenotype observed in PND2 offspring persisted through adult life, female and male offspring from Eedwt/wt, Eedwt/del and Eed del/del growing oocytes were weighed at PND30, PND49 and PND130 (Fig. 5a–b, Additional file 1: Figure S2). HET offspring from Eeddel/del growing oocytes remained significantly heavier than age-matched WT offspring from Eedwt/wt oocytes at PND30 (females), PND49 and PND130 (Additional file 1: Figure S2). However, there was no statistical difference in male weights at PND30, a growth period encompassing the transition of adolescence to adulthood in mice. Although the reason for this discrepancy remains unknown, individual variation in growth and maturation at this time point may explain the greater individual variability evident in the data. Notwithstanding this limitation, our findings are consistent with Cohen-Gibson patients , as the weight differences between animals were partially ameliorated through time and between sexes such that there was no longer a consistent statistically significant difference between HET offspring from Eedwt/del and Eed del/del growing oocytes at PND30 and at older ages (Fig. 5a–b, Additional file 1: Figure S2). Calculation of the change in average sex-specific weights in each group at PND30 and PND49 relative to PND2, detected no difference in the growth rates of WT and HET offspring produced from Eed wt/wt and Eed del/del growing oocytes (Fig. 5c–d).
In humans, de novo germline mutations in EZH2 lead to Weaver Syndrome (OMIM 28229590) and de novo germline mutations in EED lead to Cohen-Gibson syndrome, both of which are characterised by overgrowth [16,17,18, 21,22,23]. To determine the germline origin of the EZH2/EED de novo mutation, comparative SNPs were analysed in the EZH2 and EED genomic regions in DNA samples from Weaver syndrome patients and their parents (Fig. 6a). Long range, targeted sequencing revealed that the causative EZH2 mutations in human patients occurred either in the maternal or paternal germlines (Fig. 6b). Although patient numbers are quite restricted in this study, parent-of-origin appeared to have no effect on birth weight in human patients Unfortunately, analysis of two patients carrying EED mutations were uninformative, one patient lacked informative SNPs in the sequenced region and data from the other patient was inconclusive as amplification of the template failed for this individual.
An earlier study in mice found that deletion of Ezh2 in growing oocytes resulted in reduced birth weight in maternal offspring . As the offspring produced after deletion of Eed in growing oocytes resulted in overgrown offspring, we independently generated offspring from Ezhfl/fl (Zp3-Cre negative), Ezh2wt/fl;Zp3-Cre and Ezh2fl/fl;Zp3-Cre mothers. Consistent with deletion of Eed in growing oocytes, deletion of Ezh2 using the same Zp3Cre strategy was highly efficient (Additional file 1: Figure S3) and resulted in substantially reduced H3K27me3 in growing oocytes. The degree to which H3K27me3 was lost in Ezh2del/del and Eeddel/del oocytes was comparable, indicating that loss of EZH2 is unlikely to be compensated for by the related protein EZH1. In contrast to Eed, oocyte-specific deletion of Ezh2 resulted in significantly reduced PND2 weight in HET offspring from Ezh2del/del growing oocytes compared to WT pups from Ezh2wt/del and Ezh2wt/wt growing oocytes. However, there was no difference in PND2 weight between HET offspring from Ezh2del/del and Ezh2wt/del growing oocytes (Fig. 7). Significantly, this indicated that deletion of Ezh2 in growing oocytes resulted in a genetically defined, heterozygous vs wild type reduction in offspring weight, rather than the epigenetically defined, heterozygous vs heterozygous difference in offspring weight observed for deletion of Eed in growing oocytes. A caveat to this is that Ezh2 HET offspring from Ezh2wt/del growing oocytes were not significantly different in weight to WT offspring from Ezh2wt/del or Ezh2wt/wt growing oocytes. Therefore, it is likely that this phenotype results from a combination of genetic and epigenetic effects mediated through the maternal germline.
We have established a model in which oocyte-specific deletion of Eed results in overgrowth in postnatal offspring. Heterozygous offspring from Eeddel/del growing oocytes had increased body weight, length, fat content and bone mineral density compared to isogenic heterozygous control offspring from Eedwt/del growing oocytes. Our findings are reminiscent of weight gains observed in offspring in other epigenetic models, including mice with disrupted maternal imprinting of H19 . This is an interesting parallel considering recent evidence has demonstrated a role for H3K27me3 in DNA methylation-independent maternal imprinting .
Growth trajectories were similar between offspring from all Eed oocyte genotypes and mice from Eeddel/del growing oocytes remained moderately larger throughout adulthood (18 weeks of age), although this overgrowth phenotype was ameliorated through time. These phenotypic effects demonstrate the importance of maternal PRC2 for establishing fetal growth patterns and provide evidence that maternal EED regulates offspring growth and postnatal outcomes. Linear regression analysis indicated that litter size did not strongly correlate with weight at PND2. Therefore, while litter size may make a minor confounding contribution to PND2 offspring size, the majority of the overgrowth phenotype was attributable to differences in EED-dependent regulation in the oocyte.
In contrast to the offspring overgrowth observed for loss of EED in the oocyte, we found that oocytes lacking EZH2 resulted in reduced offspring birth weight. The observed EZH2-dependent growth restriction was consistent with a previous study . However, the pattern of inheritance for EED-mediated effects to offspring was not consistent with Ezh2 regulating a clear epigenetic inheritance effect from the oocyte as HET offspring generated from Ezh2 wt/del and Ezh2 del/del oocytes were equivalent, but WT and HET offspring differed in birth weight. In contrast, loss of Eed in oocytes resulted in postnatal overgrowth in HET offspring generated from Eed del/del oocytes compared to HET offspring from Eed wt/del oocytes. As these offspring were isogenic and derived from oocytes that lacked EED during oocyte growth and preimplantation development, it appears that these effects are mediated by altered epigenetic patterning in the oocyte and early life. Furthermore, this Eed-mediated, inherited phenotype persisted into the later stages of life (> 18 weeks old). Moreover, while in the previous study, the EZH2-mediated growth restriction was resolved by 4 weeks of age , the increased offspring size initiated by lack of EED in the oocyte was more persistent and affected animals through adulthood.
Despite the dependence of PRC2 on EZH2 and EED for catalysing H3K27me3, offspring generated from Eeddel/del and Ezh2del/del growing oocytes produced remarkably different overgrowth and growth restriction phenotypes in offspring. Growth restricted offspring were produced from Ezh2del/del oocytes, but there was no difference in weight between HET and HET offspring from Ezh2del/wt and Ezh2del/del oocytes, or between WT and WT offspring from Ezh2del/wt and Ezh2wt/wt oocytes. However, there was a significant difference in weight between genetically different HET and WT offspring. These observations indicate that the growth restriction observed in Ezh2 HET vs WT offspring is dependent on Ezh2 heterozygosity rather than on epigenetic or maternal factor differences identified for EED. One possible explanation is that haploinsufficiency of PRC2 occurs due to compensation for the loss of EZH2 by EZH1 in Ezh2del/del growing oocytes, whereas PRC2 function is entirely lost in Eeddel/del oocytes. However, this scenario seems unlikely as H3K27me3 levels were similar in Eeddel/del and Ezh2del/del oocytes, but offspring phenotypes differed substantially for these two models. Moreover, although growth restriction was related to Ezh2 heterozygosity in offspring, Eed heterozygosity did not result in a growth phenotype.
An alternative explanation is that the functions of EZH2 and EED differ in these models. EED is known to function only in PRC2 and is essential for PRC2 function . However, EZH2 can methylate non-histone protein targets [41,42,43] and result in phenotypic change that is PRC2 independent. For example, EZH2 methylates PLZF in T cells, leading to the ubiquitination of PLZF and its subsequent degradation . Furthermore, T cell-specific deletion of Ezh2 did not affect H3K27me3 levels and induced rapid expansion of natural killer T (NKT) cells. In contrast, conditional deletion of Suz12 or Eed halted NKT cell development due to PRC2 destabilisation . Furthermore, PLZF exerts growth-suppressive activities  and regulates limb-axial skeletal patterning in mice . Interestingly, skeletal development is compromised when Ezh2 is deleted from mesenchymal stem cells in mice, resulting in reduced skeletal size and reduced body weight . Therefore, although further analyses are required to determine the role of Ezh2 in this system, it remains possible that the growth restriction observed in Ezh2 heterozygotes derived from Ezh2del/del oocytes may result from PRC2 independent EZH2 effects in the embryo rather than reduced H3K27me3 in the oocyte.
Significantly, de novo germline mutations in EED or EZH2 in humans lead to Cohen-Gibson and Weaver syndromes, characterised by overgrowth, skeletal defects and learning/cognitive disabilities [16,17,18, 21, 22]. As deletion of Eed in the maternal germline also resulted in offspring overgrowth, it appears likely that this mouse model, at least in part, reflects the developmental defects typified in Cohen-Gibson syndrome patients. In mice, this appears to be either an epigenetically inherited phenomenon, or to result from loss of a maternal factor activity of EED, for which the only known activity is to regulate epigenetic state through H3K27 methylation. De novo germline mutations in human EED lead to Cohen-Gibson, also characterised by fetal overgrowth. Moreover, EED is very highly conserved between mice and humans and loss of EED in the oocyte results in similar offspring phenotypes in both species. Therefore, the EedZP3Cre mouse model appears highly relevant for both the study of maternal epigenetic inheritance and Cohen-Gibson syndrome in humans.
An apparent anomaly in the mouse and human models is that murine offspring produced from oocytes lacking Ezh2 were born with growth restriction, but patients with de novo missense mutations in either EZH2 or EED were born with overgrowth [16,17,18]. The sequencing data presented here demonstrated that in two affected families, de novo p.Pro132Ser (c.394 C>T) and p.Arg684Cys (c.2050 C>T) changes in EZH2 were maternally inherited. Both patients presented with increased length at birth (94th and 95th percentiles) and the patient with the p.Arg684Cys (c.2050 C>T) mutation was in the 95th percentile for weight and was delivered preterm. Consistent with the mouse model, which lacks the EZH2 SET domain and is catalytically inactive , the maternally derived EZH2 mutations in both patients reduced histone methyltransferase activity, although the impact of Pro132Ser was less severe than that mediated by Arg684Cys . The reasons for the differing phenotypes in mouse and human (i.e., growth restriction in mouse and overgrowth in Weaver patients) are not clear. However, mesenchymal stem cell-specific Ezh2 deletion in mice resulted in large offspring with increased size and weight in heterozygous animals, but reduced size and weight in homozygous animals, demonstrating that partial loss of Ezh2 activity (via heterozygosity of a loss-of-function mutation) and complete loss of Ezh2 (via deletion) have differing impacts on skeletal size and birth weight . Whether similar effects of EZH2 dosage in human patients result in variable skeletal development remains to be determined. However, human data in DECIPHER associates copy number variations in either EZH2 or EED suggesting that both over- and under-growth can be related to alterations of PRC2 dosage .
Epigenetic modifications are established and removed through the activity of specific chromatin modifying complexes, which can be altered by mutations, dysregulated transcription, drugs or other environmental influences. The most prominent examples of epigenetic dysregulation have been found in cancer, with mutations in human tumours commonly detected in genes that regulate chromatin organisation . For example, gain-of-function mutations in EZH2, EED or SUZ12 occur in multiple cancer types, and both EZH2 and EED are current targets for cancer therapy . A prominent example is provided by Tazemetostat (Epizyme, Inc.), an EZH1/2 inhibitor in stage I/II clinical trials for treatment of refractory malignant mesothelioma, lymphomas and a number of other tumours . Significantly, as observed in Eed and Ezh2 deleted oocytes, H3K27me3 is substantially and rapidly reduced in growing oocytes of mice treated with pre-clinical doses of Tazemetostat for just 10 days . While the impacts of maternal Tazemetostat treatment on offspring outcomes remain unknown, our study demonstrates that loss of H3K27me3 in growing oocytes substantially affects life-long offspring outcomes.
Importantly, a range of cancers affect men and women of reproductive age, including tumours carrying gain-of-function mutations in EZH2. Our mouse model provides evidence that loss of maternal PRC2 function in the germline of adult females disrupts H3K27me3 and results in defects in offspring growth. Based on these findings, it is likely that drugs targeting PRC2 will alter the epigenome in the oocytes of patients undergoing cancer therapy, raising the possibility that these drugs will affect health outcomes in offspring should the patient conceive during or soon after treatment . Indeed, since the oocyte pool in all women is finite, and oocyte growth and maturation spans 12 months in humans, these drugs may have effects in patients’ oocytes long after the termination of treatment. Further evaluation of the impacts of drugs that target PRC2 and other epigenetic modifier enzymes on oocytes and sperm is essential to assess the risks of these drugs in reproductive biology. Such studies would reveal the potential risks of drugs that target epigenetic modifier complexes and would facilitate development of more informed clinical approaches for patients of reproductive age , potentially including germline preservation or avoidance of pregnancy after treatment has been completed.
Using a genetic approach that deleted Eed specifically in the growing oocyte, we have demonstrated that this conserved and developmentally important epigenetic modifier mediates programming effects in the oocyte and in the earliest stages of development that are important for life-long outcomes in mice. Further understanding of this and other similar models is essential for determining how epigenetic modifiers regulate early life and long-term developmental outcomes in human health and the developmental origins of disease.
Mouse strains, animal housing, breeding and ethics
Mice were housed at Monash Medical Centre Animal Facility using a 12-h light-dark cycle. Food and water were available ad libitum and room temperature was 21–23 °C with controlled humidity. With the exception of breeding pairs and newborn mice (up to 3 weeks old), males and females were weaned at 21 days and kept in cages of up to five individuals. All animal work was undertaken in accordance with Monash University Animal Ethics Committee (AEC) approvals. Zp3Cre mice (C57BL/6-Tg 93knw/J; Jackson Labs line 003651) were constructed by Professor Barbara Knowles and obtained from The Jackson Laboratory. Eed floxed mice (Eedfl/fl) (B6; 129S1-Eedtm1Sho/J; Jackson Labs line 0022727) were constructed by Stuart Orkin and colleagues  and obtained from the Jackson Laboratory. Ezh2 floxed mice (Ezh2fl/fl) were constructed by Alexander Tarakhovsy and colleagues . The Eed and Ezh2 lines were backcrossed to a pure C57BL6/J and shared with us by Rhys Allen and Marnie Blewitt, Walter and Eliza Hall Institute for Medical Research, Melbourne.
Colony maintenance animals were genotyped via tail collection at PND2 or ear punch at weaning by Transnetyx (Cordova, TN) using real-time PCR assays designed for each gene (details available on request). Assays were designed based on the genomic structure of Eed and Ezh2 in relation to the conditional genetic modifications established in  and , respectively.
Tissue fixation and embedding
Ovaries were fixed in 4% paraformaldehyde (PFA) in PBS overnight at 4 °C. Samples were then washed in PBS and left in a 30% sucrose solution overnight at 4 °C. Samples were then placed in disposable cryostat moulds (Sakura Finetek, #4565) filled with OCT (Sakura Finetek, #4583) and frozen in dry ice. Blocks were stored at − 80 °C.
Eight micron sections were cut from OCT embedded ovaries fixed in 4% PFA, mounted on Superfrost Plus slides and dried for 5 min before immersing in 1 × PBS. Sections were then permeabilised by incubation in 1% Triton × 100 (Sigma, #T8787) in PBS for 10 min at room temperature (RT). Slides were washed in PBS. Sections were blocked in PBS containing 5% BSA (Sigma, #A9647) and 10% donkey serum (Sigma, #D9663) and incubated for 45 min at RT. Blocking solution was replaced by PBS containing 1% BSA and appropriately diluted H3K27me3 antibody (1:400, rabbit anti-H3K27me3 Cell Signaling Technologies #C36B11) and incubated for 1 h at RT. Slides were washed three times for 5 min in PBS and secondary antibodies diluted in 1% BSA in 1 × PBS according to antibody dilutions (1:300, donkey anti-rabbit 594, Alexa Fluor Life Technologies #A21207. Secondary antibody incubation was carried out in a dark box for 1 h at RT. Slides were washed three times in PBS (5 min each wash) and mounted in ProLong Gold® containing DAPI (Life Technologies, #P36931) and left in a dark box overnight to dry. For control slides, only a secondary antibody was applied. Confocal images were taken as single optical sections using a Nikon® C1 inverted Confocal microscope. All pictures were taken at × 80, using a × 40 oil immersion lens.
Phenotypic analysis of offspring
Offspring were weighed at PND2 and/or 30, 49 and 130. At PND30, 49 and 130, female and male mice were analysed separately to account for age-related sex-specific differences. Crown to rump (rump defined by base of tail) and nose to rump measurements were taken by one individual to minimise variances in method. Maternal genotypes were concealed to ensure no bias during collection of data.
Lunar PIXImus DEXA (bone mineral density, lean muscle, fat content)
PND2 offspring were measured, euthanised and stored at − 20 °C. Lunar PIXImus for small animals was used to measure bone and tissue composition using dual-energy X-ray absorptiometry (DEXA). Quality control is ensured by calibrating Lunar PIXImus to a QC phantom upon each run. Animal is place inside the region of interest (ROI), and scan is completed, n = 61. All data was collected using Lunar Piximus 2.0 software.
Long-range PCR and sequencing of patient and parental DNA samples
Regions spanning the defined mutations in the EZH2 gene of patients were amplified using long-range PCR in patient and parental DNA samples and processed for single molecule real-time (SMRT) sequencing (Pacific Bioscience) by the Leiden Genome Technology Centre, Department of Human Genetics, Leiden University Medical Center, the Netherlands. All samples were obtained by Professor William Gibson under University of British Columbia and British Columbia Children’s Hospital Human Ethics approval numbers H08-00784, H09-01228 and H10-03215, University of British Columbia, Vancouver, Canada. Sample preparation and workflow was carried out as previously described . Primer sequences are included in Additional file 1.
One-way ANOVA plus Tukey’s post hoc test was used to statistically analyse all quantitative data. Pearson r correlation test used to analyse weight vs. litter size. PRISM software v6.0e; GraphPad Prism 7 was used to analyse and graph data sets. Where data are presented graphically, statistically significant (P < 0.05) post hoc outcomes are represented by an asterisk.
A linear mixed model with litter ID as a random effect was used to test for differences among offspring groups while accounting for the relationship between offspring weight and litter size. Tests of significance for fixed effects were undertaken using Satterthwaite approximation for degrees of freedom and type III sums of squares for ANOVA. Linear mixed models were implemented in R v3.3.2 (R Core ) using lme4 v1.1-13 and lmerTest v2.0-33 [55, 56].
Embryonic Ectoderm Development
Enhancer of Zeste 1
Enhancer of Zeste 2
Trimethylated lysine 27 on histone 3
Mouse vasa homologue (also known as DDX4)
- Nsd :
No significant difference
Polycomb repressive complex 2
Supressor of Zeste 12
Zona pellucida 3 Cre recombinase
Stegemann R, Buchner DA. Transgenerational inheritance of metabolic disease. Semin Cell DevBiol. 2015;43:131–40.
Fernandez-Twinn DS, Constancia M, Ozanne SE. Intergenerational epigenetic inheritance in models of developmental programming of adult disease. Semin Cell DevBiol. 2015;43:85–95.
Gapp K, Woldemichael BT, Bohacek J, Mansuy IM. Epigenetic regulation in neurodevelopment and neurodegenerative diseases. Neuroscience. 2014;264:99–111.
Keverne EB. Significance of epigenetics for understanding brain development, brain evolution and behaviour. Neuroscience. 2014;264:207–17.
Barker DJ. The origins of the developmental origins theory. J Intern Med. 2007;261(5):412–7.
Pasini D, Bracken AP, Jensen MR, Denchi EL, Helin K. Suz12 is essential for mouse development and for EZH2 histone methyltransferase activity. EMBO J. 2004;23(20):4061–71.
Faust C, Schumacher A, Holdener B, Magnuson T. The Eed mutation disrupts anterior mesoderm production in mice. Development. 1995;121(2):273–85.
O'Carroll D, Erhardt S, Pagani M, Barton SC, Surani MA, Jenuwein T. The polycomb-group gene Ezh2 is required for early mouse development. Mol Cell Biol. 2001;21(13):4330–6.
Lesch BJ, Dokshin GA, Young RA, McCarrey JR, Page DC. A set of genes critical to development is epigenetically poised in mouse germ cells from fetal stages through completion of meiosis. Proc Natl Acad Sci U S A. 2013;110(40):16061–6.
Sachs M, Onodera C, Blaschke K, Ebata KT, Song JS, Ramalho-Santos M. Bivalent chromatin marks developmental regulatory genes in the mouse embryonic germline in vivo. Cell Rep. 2013;3(6):1777–84.
Ng JH, Kumar V, Muratani M, Kraus P, Yeo JC, Yaw LP, Xue K, Lufkin T, Prabhakar S, Ng HH. In vivo epigenomic profiling of germ cells reveals germ cell molecular signatures. Dev Cell. 2013;24(3):324–33.
Erkek S, Hisano M, Liang CY, Gill M, Murr R, Dieker J, Schubeler D, van der Vlag J, Stadler MB, Peters AH. Molecular determinants of nucleosome retention at CpG-rich sequences in mouse spermatozoa. Nat Struct Mol Biol. 2013;20(7):868–75.
Brykczynska U, Hisano M, Erkek S, Ramos L, Oakeley EJ, Roloff TC, Beisel C, Schubeler D, Stadler MB, Peters AH. Repressive and active histone methylation mark distinct promoters in human and mouse spermatozoa. Nat Struct Mol Biol. 2010;17(6):679–87.
Hammoud SS, Nix DA, Zhang H, Purwar J, Carrell DT, Cairns BR. Distinctive chromatin in human sperm packages genes for embryo development. Nature. 2009;460(7254):473–8.
Erhardt S, Su IH, Schneider R, Barton S, Bannister AJ, Perez-Burgos L, Jenuwein T, Kouzarides T, Tarakhovsky A, Surani MA. Consequences of the depletion of zygotic and embryonic enhancer of zeste 2 during preimplantation mouse development. Development. 2003;130(18):4235–48.
Cohen ASA, Tuysuz B, Shen YQ, Bhalla SK, Jones SJM, Gibson WT. A novel mutation in EED associated with overgrowth. J Hum Genet. 2015;60(6):339–42.
Cohen ASA, Gibson WT. EED-associated overgrowth in a second male patient. J Hum Genet. 2016;61(9):831–4.
Cooney E, Bi WM, Schlesinger AE, Vinson S, Potocki L. Novel EED mutation in patient with Weaver syndrome. Am J Med Genet A. 2017;173(2):541–5.
Imagawa E, Higashimoto K, Sakai Y, Numakura C, Okamoto N, Matsunaga S, Ryo A, Sato Y, Sanefuji M, Ihara K, et al. Mutations in genes encoding polycomb repressive complex 2 subunits cause Weaver syndrome. Hum Mutat. 2017;38(6):637–48.
Tatton-Brown K, Loveday C, Yost S, Clarke M, Ramsay E, Zachariou A, Elliott A, Wylie H, Ardissone A, Rittinger O, et al. Mutations in epigenetic regulation genes are a major cause of overgrowth with intellectual disability. Am J Hum Genet. 2017;100(5):725–36.
Gibson WT, Hood RL, Zhan SH, Bulman DE, Fejes AP, Moore R, Mungall AJ, Eydoux P, Babul-Hirji R, An JH, et al. Mutations in EZH2 cause Weaver syndrome. Am J Hum Genet. 2012;90(1):110–8.
Tatton-Brown K, Hanks S, Ruark E, Zachariou A, Duarte SD, Ramsay E, Snape K, Murray A, Perdeaux ER, Seal S, et al. Germline mutations in the oncogene EZH2 cause Weaver syndrome and increased human height. Oncotarget. 2011;2(12):1127–33.
Cohen AS, Yap DB, Lewis ME, Chijiwa C, Ramos-Arroyo MA, Tkachenko N, Milano V, Fradin M, McKinnon ML, Townsend KN, et al. Weaver syndrome-associated EZH2 protein variants show impaired histone methyltransferase function in vitro. Hum Mutat. 2016;37(3):301–7.
Smith P, Wilhelm D, Rodgers RJ. Development of mammalian ovary. J Endocrinol. 2014;221(3):R145–61.
Eppig JJ. Oocyte control of ovarian follicular development and function in mammals. Reproduction. 2001;122(6):829–38.
Petronczki M, Siomos MF, Nasmyth K. Un menage a quatre: the molecular biology of chromosome segregation in meiosis. Cell. 2003;112(4):423–40.
Kim K-H, Lee K-A. Maternal effect genes: findings and effects on mouse embryo development. Clin ExpReprodMed. 2014;41(2):47–61.
Prokopuk L, Western PS, Stringer JM. Transgenerational epigenetic inheritance: adaptation through the germline epigenome? Epigenomics. 2015;7(5):829–46.
Hogg K, Western PS. Refurbishing the germline epigenome: out with the old, in with the new. Semin Cell DevBiol. 2015;45:104–13.
Posfai E, Kunzmann R, Brochard V, Salvaing J, Cabuy E, Roloff TC, Liu ZC, Tardat M, van Lohuizen M, Vidal M, et al. Polycomb function during oogenesis is required for mouse embryonic development. Genes Dev. 2012;26(9):920–32.
Lewandoski M, Wassarman KM, Martin GR. Zp3-cre, a transgenic mouse line for the activation or inactivation of loxP-flanked target genes specifically in the female germ line. Curr Biol. 1997;7(2):148–51.
Bao S, Tang F, Li X, Hayashi K, Gillich A, Lao K, Surani MA. Epigenetic reversion of post-implantation epiblast to pluripotent embryonic stem cells. Nature. 2009;461(7268):1292–5.
Kemler R, Hierholzer A, Kanzler B, Kuppig S, Hansen K, Taketo MM, de Vries WN, Knowles BB, Solter D. Stabilization of beta-catenin in the mouse zygote leads to premature epithelial-mesenchymal transition in the epiblast. Development. 2004;131(23):5817–24.
Inoue A, Jiang L, Lu F, Suzuki T, Zhang Y. Maternal H3K27me3 controls DNA methylation-independent imprinting. Nature. 2017;547(7664):419–24.
Prokopuk L, Stringer JM, Hogg K, Elgass KD, Western PS. PRC2 is required for extensive reorganization of H3K27me3 during epigenetic reprogramming in mouse fetal germ cells. Epigenetics Chromatin. 2017;10:1–20.
Xie H, Xu J, Hsu JH, Nguyen M, Fujiwara Y, Peng C, Orkin SH. Polycomb repressive complex 2 regulates normal hematopoietic stem cell function in a developmental-stage-specific manner. Cell Stem Cell. 2014;14(1):68–80.
Racki WJ, Richter JD. CPEB controls oocyte growth and follicle development in the mouse. Development. 2006;133(22):4527–37.
Li R, Albertini DF. The road to maturation: somatic cell interaction and self-organization of the mammalian oocyte. Nat Rev Mol Cell Biol. 2013;14(3):141–52.
Leighton PA, Ingram RS, Eggenschwiler J, Efstratiadis A, Tilghman SM. Disruption of imprinting caused by deletion of the H19 gene region in mice. Nature. 1995;375(6526):34–9.
Cao Q, Wang X, Zhao M, Yang R, Malik R, Qiao Y, Poliakov A, Yocum AK, Li Y, Chen W, et al. The central role of EED in the orchestration of polycomb group complexes. Nat Commun. 2014;5:3127.
Yan JL, Ng SB, Tay JLS, Lin BH, Koh TL, Tan J, Selvarajan V, Liu SC, Bi CL, Wang S, et al. EZH2 overexpression in natural killer/T-cell lymphoma confers growth advantage independently of histone methyltransferase activity. Blood. 2013;121(22):4512–20.
Bracken AP, Pasini D, Capra M, Prosperini E, Colli E, Helin K. EZH2 is downstream of the pRB-E2F pathway, essential for proliferation and amplified in cancer. EMBO J. 2003;22(20):5323–35.
Vasanthakumar A, Xu DK, Lun ATL, Kueh AJ, van Gisbergen KPJM, Iannarella N, Li XF, Yu L, Wang D, Williams BRG, et al. A non-canonical function of Ezh2 preserves immune homeostasis. EMBO Rep. 2017;18(4):619–31.
Shaknovich R, Yeyati PL, Ivins S, Melnick A, Lempert C, Waxman S, Zelent A, Licht JD. The promyelocytic leukemia zinc finger protein affects myeloid cell growth, differentiation, and apoptosis. Mol Cell Biol. 1998;18(9):5533–45.
Barna M, Hawe N, Niswander L, Pandolfi PP. Plzf regulates limb and axial skeletal patterning. Nat Genet. 2000;25(2):166–72.
Hemming S, Cakouros D, Codrington J, Vandyke K, Arthur A, Zannettino A, Gronthos S. EZH2 deletion in early mesenchyme compromises postnatal bone micro architecture and structural integrity and accelerates remodeling. FASEB J. 2017;31(3):1011.
Su IH, Basavaraj A, Krutchinsky AN, Hobert O, Ullrich A, Chait BT, Tarakhovsky A. Ezh2 controls B cell development through histone H3 methylation and Igh rearrangement. Nat Immunol. 2003;4(2):124–31.
Firth HV, Richards SM, Bevan AP, Clayton S, Corpas M, Rajan D, Van Vooren S, Moreau Y, Pettett RM, Carter NP. DECIPHER: database of chromosomal imbalance and phenotype in humans using Ensembl resources. Am J Hum Genet. 2009;84(4):524–33.
Jones PA, Issa JPJ, Baylin S. Targeting the cancer epigenome for therapy. Nat Rev Genet. 2016;17(10):630–41.
Knutson SK, Kawano S, Minoshima Y, Warholic NM, Huang KC, Xiao YH, Kadowaki T, Uesugi M, Kuznetsov G, Kumar N, et al. Selective inhibition of EZH2 by EPZ-6438 leads to potent antitumor activity in EZH2-mutant non-Hodgkin lymphoma. Mol Cancer Ther. 2014;13(4):842–54.
Prokopuk L, Hogg K, Western PS. Pharmacological inhibition of EZH2 disrupts the female germline epigenome. Clin Epigenetics. 2018; In press
Western PS. Epigenomic drugs and the germline: collateral damage in the home of heritability? Mol Cell Endocrinol. 2018; In press
Buermans HPJ, Vossen RHAM, Anvar SY, Allard WG, Guchelaar HJ, White SJ, den Dunnen JT, Swen JJ, van der Straaten T. Flexible and scalable full-length CYP2D6 long amplicon PacBio sequencing. Hum Mutat. 2017;38(3):310–6.
Team RC. R: a language and environment for statistical computing [computer software], v3.3.2 edn. Vienna: R Foundation for Statistical Computing; 2016.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48.
Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest package: tests in linear mixed effects models. J Stat Softw. 2017;82(13):1–26.
We thank Monash Animal Research Platform staff for assistance with mouse care; MMC Histology Platform and MMI Micro Imaging Facility for technical advice. We thank Stuart Orkin, Barbara Knowles and Alexander Tarakhovsy for respectively sharing the Eedfl/fl, Ezh2fl/fl, and ZP3Cre lines and Marnie Blewitt, Rhys Allen, John Carroll and Eileen McLaughlin for sharing mice. We also thank the patients and families involved in the study.
This work was supported by the National Health and Medical Research Grants 1144966 and 1144887 awarded to PW, funding from the Monash University Faculty of Medicine, Nursing and Health Sciences funding granted to PW and the Victorian Government’s Operational Infrastructure Support Program. LP was supported by an Australian Postgraduate Award.
Availability of data and materials
All data is available from the corresponding author.
LP is now located at the Stem Cells and Cancer Division of the Walter and Eliza Hall Institute of Medical Research, Parkville, Victoria, Australia.
Ethics approval and consent to participate
Consent was obtained for use of all patient samples by Professor William Gibson under University of British Columbia and British Columbia Children’s Hospital Human Ethics approval numbers H08-00784, H09-01228 and H10-03215, University of British Columbia, Vancouver, Canada. Animal work was undertaken in accordance with Monash University Animal Ethics Committee (AEC) approval MMCA-2016-18.
Consent for publication
All authors have approved this manuscript and consented to publication of the data it contains.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Prokopuk, L., Stringer, J.M., White, C.R. et al. Loss of maternal EED results in postnatal overgrowth. Clin Epigenet 10, 95 (2018). https://doi.org/10.1186/s13148-018-0526-8
- Epigenetic inheritance