Early Developmental and Evolutionary Origins of Gene Body DNA Methylation Patterns in Mammalian Placentas

The placenta is vital for the proper development of the fetus, not only facilitating the exchange of nutrients, oxygen, and waste between the mother and the fetus but also acting as an interface to the maternal immune system and regulating fetal growth by excreting hormones and growth factors. DNA methylation is important for both placental and embryonic development as loss of proteins involved in DNA methylation can result in placental dysmorphology and early embryonic death. The human placenta has a unique DNA methylation landscape characterized by alternating regions of low methylation, covering silent genes with tissue-specific developmental functions, and high methylation, covering active genes. In order to better understand the significance of this DNA methylation landscape in the human placenta, we performed a cross-species comparison of DNA methylation in mammalian placentas, oocytes, and early embryos from this and other studies. Although the levels and extent of hypomethylation differed between mammalian placentas, what we found to be highly conserved was relatively higher methylation levels over active genes. These same genes also had high methylation in the opossum extraembryonic membrane, a primitive placenta, as well as oocytes and early embryos, suggesting that high methylation over these genes predated placental mammals and is established very early in development.

Published in the journal: . PLoS Genet 11(8): e32767. doi:10.1371/journal.pgen.1005442
Category: Research Article
doi: 10.1371/journal.pgen.1005442


The placenta is vital for the proper development of the fetus, not only facilitating the exchange of nutrients, oxygen, and waste between the mother and the fetus but also acting as an interface to the maternal immune system and regulating fetal growth by excreting hormones and growth factors. DNA methylation is important for both placental and embryonic development as loss of proteins involved in DNA methylation can result in placental dysmorphology and early embryonic death. The human placenta has a unique DNA methylation landscape characterized by alternating regions of low methylation, covering silent genes with tissue-specific developmental functions, and high methylation, covering active genes. In order to better understand the significance of this DNA methylation landscape in the human placenta, we performed a cross-species comparison of DNA methylation in mammalian placentas, oocytes, and early embryos from this and other studies. Although the levels and extent of hypomethylation differed between mammalian placentas, what we found to be highly conserved was relatively higher methylation levels over active genes. These same genes also had high methylation in the opossum extraembryonic membrane, a primitive placenta, as well as oocytes and early embryos, suggesting that high methylation over these genes predated placental mammals and is established very early in development.


In eutherian mammals the placenta plays a vital role in not only the transfer of nutrients and waste between mother and offspring but also as a protective layer between the maternal and fetal immune systems during fetal development. Despite this, the gross morphologies and cellular characteristics of the maternal/fetal interface are quite diverse [1,2] and have undergone multiple instances of both divergent and convergent evolution [3]. Marsupials also have an extraembryonic membrane (EEM) that, although short-lived, is also important for nutrient exchange, originates from a trophectoderm layer in the early embryo, and has been argued to be a true placenta [4,5,6,7].

DNA methylation is essential for proper embryo and placenta development. The offspring of Dnmt3a conditional knockout mothers die in utero by E11.5 [8]. The offspring of Dnmt3L null mothers die at E9.5 from placental abnormalities and/or imprinting defects [9,10,11]. Another DNMT important for normal placenta development is DNMT1o, an oocyte-specific isoform of DNMT1 that is present in eutherians as well as metatherians [12]. Loss of Dnmt1o in mice results in widespread placental dysmorphology [13,14].

It has long been known that human and mouse placentas are hypomethylated compared to other tissues [15,16,17,18,19]. However, recent analysis of human placenta has shown a large-scale pattern of PMDs and HMDs that are often over 200 kb in length and can cover entire genes and gene clusters [20,21]. Human placenta PMDs cover tissue-specific genes that are transcriptionally repressed. These findings provide a unique opportunity to use DNA methylation to study not only the evolution of transcriptional regulation in this developmentally important organ but also the molecular similarities between species with morphologically distinct placentas. One study has shown that gene-specific methylation in mammalian placentas tracks with phylogeny more than placental morphology [22], but otherwise little is known about the evolution of DNA methylation and transcriptional regulation in placenta.

The placenta derives from the trophectoderm layer of the blastocyst in the early embryo, before the implantation stage. The DNA methylation patterns of the early embryo have recently been elucidated. Human and mouse oocytes have low levels of methylation which decrease even further in the fertilized embryo until after the blastocyst stage [23,24,25,26,27]. Therefore, it is thought that the placenta may never undergo the wave of remethylation that occurs in cells with other somatic tissue fates [28].

We performed MethylC-seq on a representative set of mammalian placentas, both with respect to evolutionary relationship and placental morphology. We found that although PMD/HMDs were not found in the placentas of many species, what was conserved was high methylation over gene bodies, particularly over those of active genes. Thus, in many mammalian placentas high methylation is found over single genes instead of over entire clusters of genes, as is the case in placentas with PMD/HMDs. Low genome-wide methylation and high gene body methylation in transcriptionally active genes was also found in the opossum EEM, suggesting that this is a conserved feature in the evolution of the mammalian placenta. Finally, analysis of DNA methylation in oocytes and preimplantation embryos shows that low global methylation but higher methylation over gene bodies is present prior to fertilization and persists through the blastocyst stage


Mammalian placentas differ in their levels and patterns of global DNA methylation

Low coverage MethylC-seq was performed on rhesus, squirrel monkey, mouse, dog, cow, and horse placentas as well as opossum EEM (S1 Table). Fig 1A shows the phylogenetic relationship between the species in this study and the type of placental morphology in each.

Genome-wide methylation patterns in mammalian placentas show both large-scale divergence and gene-specific similarities.
Fig. 1. Genome-wide methylation patterns in mammalian placentas show both large-scale divergence and gene-specific similarities.
(A) Phylogenetic tree of the species studied and the classification of their placenta types. Branch lengths are not to scale. (B) Density curves of average percent methylation in non-overlapping 20 kb windows in mammalian placentas. For comparison, the methylation distribution for brain tissue is shown in beige for some species and human cord blood is shown in aquamarine. The opossum brain sample was fetal whereas the other brain samples were postnatal. The interquartile range and medians are shown as black bars and white dots, respectively, for the placenta samples. Mouse brain data is from Hon et al. [19] (GSE42836), human placenta data is from Schroeder et al [21] (GSE25930). (C) Comparison of global methylation patterns in select species after liftOver to the human genome and smoothing (for full figure, see S4 Fig). Human placenta PMDs, as determined by HMM, are shown in black bars at top. (D) Comparison of methylation patterns at the CNTNAP2 locus in select species (for full figure, see S5 Fig). Raw species CpG site methylation wig data were graphed on the UCSC Genome Browser without preprocessing. Black line represents 50% methylation. Genes of interest are orange. EEM = extra-embryonic membrane.

In order to compare the pattern of DNA hypomethylation across mammalian placentas, average methylation was calculated for non-overlapping 20 kb windows tiled across the autosomes for each species. To compare placenta methylation levels to those of a somatic tissue, we also performed MethylC-seq on adult cow and dog cerebrum samples as well as fetal opossum whole brain. Fig 1B shows that mammalian placentas vary greatly in their overall amount of DNA methylation, although all were globally hypomethylated compared to brain, including the opossum EEM. All placentas had lower than 66% average methylation (calculated as the average methylation of all CpG sites in all reads that aligned to the genome), but cow in particular had a remarkably low level of methylation in placenta (30%), even though cow brain global methylation is comparable to that of other mammals. In spite of this, methylation levels in regulatory elements such as CpG islands and promoters are similar between species (S1 Fig) and repetitive elements in general follow similar patterns, with LINEs and LTRs having methylation levels similar to those of non-repetitive regions and SINEs having relatively higher methylation (S2 Fig).

While human, squirrel monkey, and dog placentas show a bimodal distribution of methylation levels indicative of PMD and HMDs, other species showed a single peak or, like rhesus, had evidence of a small number of HMDs seen as a small secondary peak at higher methylation values. These results were similar for window lengths of 5–50 kb (S3 Fig). Since the species showing PMD/HMD divisions were not necessarily the most closely related to humans, we next examined the conservation of methylation patterns at the chromosome and gene level. Homologous chromosomal regions of each species were mapped to the human genome using the liftOver program (see Methods). Fig 1C shows that, despite the weak evidence for PMD/HMDs in the rhesus methylation density curves, methylation patterns across the rhesus chromosomes follow the patterns in human placenta remarkably well, with a Pearson correlation of 0.69 in 20 kb windows, the highest of all the species in this study (S2 Table). Interestingly, despite the vast differences in the medians and interquartile ranges of window methylation levels between the different species, the correlations are remarkably high, with mouse placenta having the least similarity with human placenta of the species studied (S2 Table and Figs 1C and S4).

Mapping a species genome to the human genome can introduce errors and artifacts, particularly in distantly-related species such as opossum (S4 Fig). To eliminate these we next examined two loci with high methylation conservation between the species without mapping to the human genome (Figs 1D, S5 and S6). Figs 1D and S5 show the CNTNAP2 locus which includes two genes involved in neuronal development, CNTNAP2 and DPP6, as well as the gene encoding Polycomb group protein EZH2. Since neuron-specific genes tend to be in PMDs in human placenta [21], we asked if hypomethylation was a conserved feature of this neuronal gene locus in the placenta. Figs 1D and S5 show that in most species both CNTNAP2 and DPP6 are in regions of relatively lower methylation compared to the surrounding locus. Most striking, however, are the short regions of relatively high methylation over individual genes or small clusters of genes that are highly conserved across all species, including the opossum EEM. EZH2, KMT2C (a histone methyltransferase expressed in placenta), PAXIP1 (a gene that maintains genome stability during mitosis), and RBM33 (a hypothetical RNA-binding protein) have higher than average methylation in all the species studied. Likewise in the DAB1 locus (S6 Fig), genes such as USP24 (a ubiquitin peptidase), FGGY (phosphorylates carbohydrates), INADL (scaffolding protein in the cell membrane), DOCK7 (a guanine nucleotide exchange factor), and ATG4A (cysteine protease required for autophagy) are highly methylated in most species.

Our previous MethylC-seq analysis of three full-term human placentas showed that PMDs and HMDs were highly reproducible across individuals and gestational ages [21]. Here we performed five additional experiments to test the sensitivity of placenta methylome data to sampling differences and cellular heterogeneity. First, MethylC-seq analysis of whole rhesus placental tissue and isolated rhesus trophoblast cells gave nearly identical results with a correlation of 0.89 (S7A, S7B, and S7C Fig). MethylC-seq results for E15.5 C57Bl6/J mouse placenta are very similar to those previously reported by Hon et al. [19] and a E11.5 C57BI6/J placenta sample (S7D, S7E and S7F Fig) with pairwise correlations over 0.87. Since cow placenta has such low methylation, to validate the results we sequenced three additional biological replicates using tissue material from a different source and found nearly identical results with correlations above 0.75 (S8A, S8B, and S8C Fig). We also sequenced two additional dog placentas from different breeds and developmental timepoints (S1 Table). Methylation patterns across the chromosomes were very similar and pairwise correlations were over 0.87 (S8D and S8F Fig). Although one of the dog placentas did not show a bimodal distribution whereas the other two did (S8E Fig), the relative pattern of methylation along chromosomes was highly similar and the pairwise correlations were high (0.89 and 0.87). Thus, the length of the range of methylation values in the curve may be more important than its bimodality. Finally, to confirm the reproducibility of human placenta MethylC-seq across samples, labs, and sequencing coverage, we compared our three previously sequenced placentas [21] to that of a higher-coverage MethylC-seq sample [29]. S9 Fig shows that all four human placenta samples have nearly identical global methylation patterns with pairwise correlations of >0.95, demonstrating that sequencing depth does not affect global methylation levels or patterns.

High methylation is associated with the gene bodies of actively transcribed genes

To determine the functions of genes with conserved high methylation in mammalian placentas, orthologous Ensembl genes in the seven species were clustered based on average gene body methylation (Fig 2A). We found a group of 3380 genes (branches A and B) that had relatively high methylation within each species. This group is enriched for functions such as cell cycle, protein localization, and protein ubiquitination (S3 Table). In contrast, genes with consistently low methylation across species tended to be transcription factors and have developmental functions (branch G), of which 30% are polycomb-regulated. Our previous study found that human placenta also has low methylation over some tissue-specific genes such as those involved in neuronal functions of synaptic transmission and ion transport [21]. Across species, however, this trend of lower methylation of neuronal genes is not as conserved since dog, horse, and opossum show mixed or higher levels of methylation in these genes (branch F).

Regions of high methylation in placentas cover gene bodies.
Fig. 2. Regions of high methylation in placentas cover gene bodies.
(A) Heatmap of average methylation in the gene bodies (introns and exons, excluding CpG islands and promoters) of orthologous genes. Only the top few GO biological processes with Benjamini p-values below 1.0E-3 are shown. For a complete list see S3 Table. Branches A and B were combined because they contain similar GO terms. (B) Comparison of percent methylation between human placenta (red curve) and rhesus placenta (purple curve). Rhesus methylation data was lifted over to the human genome. Vertical purple lines show large chromosomal breaks in synteny between the two species. The fourth ring in shows regions of higher (blue) and lower (red) methylation in human placenta compared to rhesus. The fourth circle in shows the locations of human genes in black. (C) Spinograms showing the probability that a 5 kb window is in a gene given that window's average percent methylation. Bars are color-coded by percent methylation and bar widths show the percentage of windows with that methylation level. Bars furthest from the 0.5 blue line marker show the most information about gene location.

An interesting question is whether conservation of long-range methylation patterns (eg. human and rhesus in Fig 1C) is due to conservation of large chromatin domains or whether such patterns arise at the gene level. To address this question, we compared regions with large chromosomal breaks in synteny with regions of differential relative methylation between human and other species (S10A Fig). Fig 2B shows a portion of human chromosome 7 with a relatively large number of chromosomal syntenic breaks between human and rhesus. Syntenic breaks are not enriched for differential methylation, suggesting that methylation patterns are established at a more local level.

Since we had shown that DNA methylation was most conserved in gene bodies, we examined more closely the relationship between DNA methylation and genes. S10B Fig shows that the boundaries of human placenta PMDs are closer to gene ends and CpG islands than expected by chance. We next asked if DNA methylation is predictive of gene location in the placentas of all species studied. Non-overlapping 5 kb windows were tiled across the autosomes of each species and classified as genic or intergenic. S11 Fig shows that intergenic regions have lower DNA methylation levels than genic regions. Figs 2C and S12 show the probability that a window is in a gene based on its average % methylation. The further the probability is from 0.5, the more informative the methylation level is in predicting the presence (above) or absence (below) of a gene. In all species studied, including opossum, high placental methylation was associated with genes, as was very low methylation since CpG islands were not removed. This is in contrast to brain where high methylation gives little to no information about the presence of a gene. Together these data indicate that DNA methylation patterns in placenta are being set at the gene level.

We next asked whether, as in human placenta, genes with high methylation are more likely to be expressed across a broad range of mammals with diverse placental anatomies [21]. Published polyA-selected RNA-seq data for human, mouse, and horse placenta and opossum EEM [30,31] were utilized to compare gene expression to gene body methylation in orthologous genes. In human placenta, expressed genes have clearly higher than average gene body methylation levels (Fig 3A, left column). To a lesser extent this can also be observed in mouse, horse, and opossum, although the range of gene body methylation values is smaller than human. To determine if these distributions are different than we would expect by chance given the marginal gene body methylation and gene expression distributions, we divided the x and y axes into 20 equally-spaced bins, counted the number of observations in each resulting quadrant, and compared that to the expected number of observations if gene body methylation and gene expression were independent (Fig 3A, right column). A co-independence test showed a statistically significant deviation from independence for all four species, but more importantly the patterns of deviation are remarkably similar between the species. In all species examined, genes with high gene body methylation are more likely to have intermediate expression than expected by chance and genes with low methylation are less likely to be expressed.

Regions of high methylation are enriched for actively transcribed genes.
Fig. 3. Regions of high methylation are enriched for actively transcribed genes.
(A) Relationship between average gene body percent methylation and gene expression in placentas. Only gene with orthologs were used. The right column shows a density scatterplot with gray and purple show areas of low and high density, respectively. Black lines show the marginal distribution of percent methylation in gene bodies. The left column shows the number of genes in each quadrant above/below what would be expected if gene expression and gene body methylation were independent. Human, mouse, and opossum expression data are from Necsulea et al. [30] (GSE43520) and horse expression data is from Wang et al. [31](GSE30243). (B) Comparison of methylation and gene expression patterns at the HECTD1 locus. Genes of interest are orange.

An example of the overlap between methylation and expression is the HECTD1 locus in Fig 3B. The HECTD1 gene encodes an E2 ubiquitin protein ligase that is important for proper placenta development in mouse [32]. HECTD1 has higher than average methylation levels and is expressed in human, mouse, and opossum. This is also true for the neighboring gene SCFD1, which encodes a protein involved in SNARE-pin assembly and vesicular transport.

Cow and mouse oocytes show gene body methylation patterns similar to those found in mammalian placentas

An important question is whether mammalian placentas inherit their distinct methylation patterns from the early embryo or the patterns instead emerge later during placental development. To answer this question we utilized a combination of mouse oocyte and early embryo methylation data [24] and human oocyte data [27] and also performed MethylC-seq on MII cow oocytes. Similar to what was observed in mouse and human oocytes [23,27], cow oocytes exhibit a bimodal methylation distribution (Figs 4A, left column, and S13), a surprising result since mouse and cow placentas do not show bimodal distributions. In addition, although both human oocytes and placenta have bimodal methylation distributions, the global patterns are quite dissimilar (S14A Fig).

Gene methylation and expression in oocytes compared to placenta in mouse and cow.
Fig. 4. Gene methylation and expression in oocytes compared to placenta in mouse and cow.
(A) Distribution of average methylation in 50 kb windows (first column) and gene bodies (second column) during early human, mouse, and cow development. The third column shows the relationship between average gene body methylation and gene expression in oocytes. Black lines show the marginal distribution of percent methylation in gene bodies. All protein-coding genes from each species were used. Human oocyte and sperm methylation data are from Okae et al. [27] (JGAS00000000006), human oocyte expression data are from Reich et al. [34] (GSE32689), mouse oocyte and sperm methylation and expression data are from Wang et al [24] (GSE56697), and cow oocyte expression data are from Graf et al [33] (GSE52415). (B) Heatmap of average gene body methylation in human, mouse, and cow oocytes. For a complete list of all GO and KEGG terms for each quadrant, see S3 Table.

We next looked at gene body methylation. A bimodal distribution of average gene body methylation was also observed in all three species' oocytes, although this is more prominent in the human and mouse oocytes (Fig 4A, center column). In mouse, where MethylC-seq data is available for multiple early developmental timepoints, gene body methylation shows bimodal distributions that persist through the inner cell mass (ICM) stage (S13 Fig) and gene body methylation patterns remain correlated between the oocyte and ICM stages (S15 Fig) suggesting that although global methylation levels drop, relative levels of methylation in gene bodies remains consistent.

Using oocyte expression data [24,33,34] we found that, similar to what was seen in placenta, methylation is enriched in gene bodies in oocytes and the mouse early embryo (S16 Fig) and genes with high gene body methylation are more likely to be expressed (Fig 4A, right column). Genes with high methylation in the oocytes of all three species are enriched for functions related to vesicle-mediated transport, protein phosphorylation, and chromosome organization (Fig 4B and S3 Table), similar to that seen in mammalian placentas (Fig 2A and S3 Table). In fact, genes with high expression in both placenta and oocytes are enriched for these same functions (S17A and S17B Fig) and there is a large amount of overlap between the genes that are highly methylated in placenta and oocytes and those that are highly expressed in placenta and oocytes. Thus, these genes appear to maintain a similar regulatory pattern in the placenta that was set up very early in development.


In this first comparative study of genome-wide methylation patterns in mammalian placentas, we have made several novel findings of relevance to understanding the association of transcription and methylation in early life. First, we confirm that hypomethylation of extraembryonic tissues compared to somatic tissues is observed across eutherian mammals and a metatherian mammal as well. Second, we demonstrate that while large differences exist between mammalian placentas in global methylation levels, higher relative methylation of active genes is an evolutionary conserved feature. Third, we show that relative methylation is predictive of genic location in placenta. Lastly, we establish that the pattern of higher gene body methylation of active genes is also observed in mammalian oocytes and persists in the preimplantation embryo.

While hypomethylation of placenta and EEM compared to somatic tissue was observed across mammalian species, global levels of placental methylation and the presence/absence of a bimodal PMD/HMD organization was diverse. The diversity of global methylation levels in mammalian placentas is not entirely surprising given the diversity of mammalian placenta morphologies. Based on current models of the mammalian radiation, it would appear that PMD/HMDs arose or were lost multiple times during placental evolution. However, our analysis of mouse and cow oocytes shows that even though some species do not have PMD/HMDs in their placentas, PMD/HMDs may still exist in the oocytes and early embryos, suggesting that these methylation patterns are later lost during placenta development in those species.

One potential limitation of this study is the inherent diversity of placental morphologies across species and issues of cellular heterogeneity. Differences in absolute methylation levels and the existence of PMD/HMDs could have to do with these confounding factors rather than real differences between species. However, analysis of placental biological replicates in multiple species shows that while absolute global levels of methylation may vary up to 20% due to tissue sampling and inter-individual differences, the methylation patterns across a chromosome remain remarkably similar. Thus, our conclusion of relatively higher methylation over active genes as a conserved feature of mammalian placentas is unaffected. Another potential limitation is the relatively low MethylC-seq coverage used in this study, as low as 1X coverage for most of the placentas used in this study (S1 Table). However, Ziller et al. experimentally determined coverage recommendations for whole-genome bisulfite sequencing, demonstrating that 1X coverage was sufficient for analysis of differentially methylated regions over 5 kb in length with methylation differences over 20% given two sequencing replicates [35].

One interesting question is why PMDs are hypomethylated. It was previously found in IMR90 fetal lung fibroblast cells that PMDs have a unique methylation distribution at individual CpG sites compared to HMDs [36]. High-coverage MethylC-seq data showed that within PMDs, each CpG site has a seemingly random level of methylation, but with some correlation in methylation between neighboring CpG sites less than 100 bp away. Moreover, the best predictor of the methylation level at a CpG site in PMDs was the distance to the nearest neighboring CpG site. In PMDs, unlike in CpG islands, CpG density is positively correlated with methylation. The authors hypothesize that PMDs may represent regions with reduced access to DNMTs and that in such an environment subtle DNMT sequence binding preferences may be easier to detect.

Why do extraembryonic and placental tissues exhibit low global methylation and high gene body methylation, but somatic tissues such as brain are globally highly methylated? One hypothesis stems from the observation that hypomethylation and well-defined PMD/HMD structures are also found in cancerous tumors [37,38,39]. Remarkably, the presence of hypomethylated PMDs is the most epigenetically defining feature of a wide range of human solid tumors [40], EBV transformation of B cells is characterized by the gain of hypomethylated PMDs [41], and brain tumors gain hypomethylated PMDs as they progress to higher metastases [42]. Since placenta shares many features in common with metastatic tumors, including rapid proliferation, invasiveness, and angiogenesis [43,44,45], perhaps a hypomethylated state would be evolutionary advantageous in a temporary tissue in which growth is highly regulated by pregnancy hormones. The temporary nature of the placenta may also make is less susceptible to the effects of retrotransposable elements. In contrast, somatic tissues in more long-lived mammalian species may have evolved globally higher methylation levels as a mechanism to prevent cancer in other tissues.


Tissue sources

The sources of the placenta and brain samples sequenced for this study are listed in S1 Table. Human subjects were approved by the UC Davis IRB (225645–17) and maternal written consent was obtained. All procedures involving animals were performed in accordance with the NIH Guide for the Care and Use of Laboratory Animals and under the approval of the University of California Davis, Animal Care and Use Committee (Animal Protocol #15639). Where possible, villus tissue was taken from the fetal side or interior (equidistant from the fetal and maternal sides) of the placenta, excluding large blood vessels, membranes, and any obvious connective tissue. Opossum EEM tissue was collected as described [46]. Samples were stored frozen at -80 degrees Celsius. Based on similarly processed and stored samples of human fetal side placenta, maternal cell contamination is expected to contribute <10% to methylation and sample location and inter-individual differences are expected to contribute to up to 20% differences in percent methylation over PMDs [21].

Rhesus trophoblast cells

Trophoblast cells were isolated from Rhesus monkey (Macaca mulatta) placental tissue (Gestation Day: 40–65 days) using procedures we have previously described [47,48,49]. All procedures involving animals were performed in accordance with the NIH Guide for the Care and Use of Laboratory Animals and under the approval of the University of California Davis, Animal Care and Use Committee (Animal Protocol #15639). These cells are 95% cytokeratin 7-positive and 5% vimentin-positive, consistent with a predominantly trophoblast population.

MII stage cow oocytes

MII oocytes were produced by in vitro maturation of GV oocytes collected from slaughter house derived ovaries according to standard protocols. Oocyte maturation was confirmed by presence of a polar body. Then, zona pellucida and first polar body were removed by incubation in 0.5% pronase solution for 2 minutes and vigorous pipetting. Zona-free oocytes were snap frozen in liquid nitrogen and stored at -80°C until DNA extraction.


DNA from rhesus, cow, horse, dog, squirrel monkey, and mouse placenta tissue was purified using Qiagen's Puregene kit. MethylC-seq libraries were made as described previously [21]. Briefly, the genomic DNA was sonicated to ~300 bp and methylated Illumina adapters were ligated to the ends. The library was bisulfite converted, amplified for 14 cycles, and sequenced on either an Illumina HiSeq or GAII. For the opossum EEM, rhesus trophoblast, human cord blood, and cow, dog, and opossum brain samples, MethylC-seq libraries were made using the Epicentre EpiGnome Methyl-Seq kit according to the manufacturer's recommendations except that 14 cycles of amplification were performed. For cow oocytes, DNA from 100 cells was purified using Zymo's Quick-gDNA MicroPrep kit and the MethylC-seq libraries were prepared using Zymo's Pico Methyl-seq Library Prep Kit according to the manufacturer's instructions. After sequencing on HiSeq2500/2000 machines, reads were mapped to the respective genomes using BS Seeker [50] and only one read per genomic position was kept to prevent clonal PCR amplification biases. CpG site methylation data were combined from both DNA strands. Because methylation data was analyzed over large genomic distances and/or smoothed, no minimum coverage was required for CpG sites used in this analysis [35].

Mapping MethylC-seq data to the human genome and determining chromosomal synteny

Individual CpGs were mapped to coordinates on the human genome using the liftOver program and species-specific liftOver chain files available on the UCSC Genome Browser [51]. Since our analysis focused on global methylation patterns, not individual CpG sites, mappings were done regardless of whether the CpG site was conserved in human. Since our analysis was focused on methylation in large syntenic regions, we removed CpG data in small inter- and intrachromosomal translocations, duplications, etc. using the species-specific synteny net files available on the UCSC Genome Browser. CpGs were removed if they were in the second level fill of the synteny net files and were less than 1 Mb in length. For graphing the cross-species global methylation patterns, these “cleaned” species liftOver data and raw human data were compressed into averages of non-overlapping 20 kb windows (windows with less than 20 CpG sites with methylation information were discarded) and smoothed in R using a kernel smoother.

For graphing the cross-species syntenic breaks, only the breaks between the large syntenic regions in the first level of the UCSC synteny net files and those between large (>1 Mb) syntenic regions in the second level were used since smaller non-syntenic regions had been removed from the methylation datasets. Due to the fact that the squirrel monkey genome (saiBol1) was not yet assembled into chromosomes, many of the “syntenic breaks” between human and squirrel monkey are actually between contigs that may or may not be on the same squirrel monkey chromosome.

Defining regions of differential methylation across species

After mapping species methylation data to the human genome, the data was smoothed and normalized. Average methylation was taken for non-overlapping 20 kb windows and a running median was computed using a width of 15 windows. The data was then scaled so that each species had the same mean and standard deviation. The smoothed and scaled methylation values were subtracted from the smoothed and scaled human methylation values and regions of differential methylation were defined as those with methylation differences over 1.5 standard deviations.

Methylation over orthologous gene bodies

Gene annotations for each species as well orthologous gene information were obtained from Ensembl's biomart. Promoter and CpG island sequences were removed from each gene before calculating average percent methylation over gene bodies. Genes with fewer than 20 remaining CpG sites with methylation information were removed from each species' dataset.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9

Attachment 10

Attachment 11

Attachment 12

Attachment 13

Attachment 14

Attachment 15

Attachment 16

Attachment 17

Attachment 18

Attachment 19

Attachment 20


1. Carter AM, Enders AC (2004) Comparative aspects of trophoblast development and placentation. Reproductive biology and endocrinology: RB&E 2: 46. 15236656

2. Turner (1876) Some General Observations on the Placenta, with especial reference to the Theory of Evolution. Journal of anatomy and physiology 11: 33–53. 17231137

3. Carter AM (2012) Evolution of placental function in mammals: the molecular basis of gas and nutrient transfer, hormone secretion, and immune responses. Physiol Rev 92: 1543–1576. doi: 10.1152/physrev.00040.2011 23073626

4. Padykula HA, Taylor JM (1982) Marsupial placentation and its evolutionary significance. Journal of reproduction and fertility Supplement 31: 95–104. 6762435

5. Selwood L, Johnson MH (2006) Trophoblast and hypoblast in the monotreme, marsupial and eutherian mammal: evolution and origins. BioEssays: news and reviews in molecular, cellular and developmental biology 28: 128–145.

6. Renfree MB (2010) Review: Marsupials: placental mammals with a difference. Placenta 31 Suppl: S21–26. doi: 10.1016/j.placenta.2009.12.023 20079531

7. Graves JA, Renfree MB (2013) Marsupials in the age of genomics. Annual review of genomics and human genetics 14: 393–420. doi: 10.1146/annurev-genom-091212-153452 23834319

8. Kaneda M, Okano M, Hata K, Sado T, Tsujimoto N, et al. (2004) Essential role for de novo DNA methyltransferase Dnmt3a in paternal and maternal imprinting. Nature 429: 900–903. 15215868

9. Arima T, Hata K, Tanaka S, Kusumi M, Li E, et al. (2006) Loss of the maternal imprint in Dnmt3Lmat-/- mice leads to a differentiation defect in the extraembryonic tissue. Dev Biol 297: 361–373. 16920095

10. Bourc'his D, Xu GL, Lin CS, Bollman B, Bestor TH (2001) Dnmt3L and the establishment of maternal genomic imprints. Science 294: 2536–2539. 11719692

11. Hata K, Okano M, Lei H, Li E (2002) Dnmt3L cooperates with the Dnmt3 family of de novo DNA methyltransferases to establish maternal imprints in mice. Development 129: 1983–1993. 11934864

12. Ding F, Patel C, Ratnam S, McCarrey JR, Chaillet JR (2003) Conservation of Dnmt1o cytosine methyltransferase in the marsupial Monodelphis domestica. Genesis 36: 209–213. 12929092

13. McGraw S, Oakes CC, Martel J, Cirio MC, de Zeeuw P, et al. (2013) Loss of DNMT1o disrupts imprinted X chromosome inactivation and accentuates placental defects in females. PLoS genetics 9: e1003873. doi: 10.1371/journal.pgen.1003873 24278026

14. Himes KP, Koppes E, Chaillet JR (2013) Generalized disruption of inherited genomic imprints leads to wide-ranging placental defects and dysregulated fetal growth. Developmental Biology 373: 72–82. doi: 10.1016/j.ydbio.2012.10.010 23085235

15. Ehrlich M, Gama-Sosa MA, Huang LH, Midgett RM, Kuo KC, et al. (1982) Amount and distribution of 5-methylcytosine in human DNA from different types of tissues of cells. Nucleic Acids Res 10: 2709–2721. 7079182

16. Fuke C, Shimabukuro M, Petronis A, Sugimoto J, Oda T, et al. (2004) Age related changes in 5-methylcytosine content in human peripheral leukocytes and placentas: an HPLC-based study. Annals of human genetics 68: 196–204. 15180700

17. Cotton AM, Avila L, Penaherrera MS, Affleck JG, Robinson WP, et al. (2009) Inactive X chromosome-specific reduction in placental DNA methylation. Hum Mol Genet 18: 3544–3552. doi: 10.1093/hmg/ddp299 19586922

18. Popp C, Dean W, Feng S, Cokus SJ, Andrews S, et al. (2010) Genome-wide erasure of DNA methylation in mouse primordial germ cells is affected by AID deficiency. Nature 463: 1101–1105. doi: 10.1038/nature08829 20098412

19. Hon GC, Rajagopal N, Shen Y, McCleary DF, Yue F, et al. (2013) Epigenetic memory at embryonic enhancers identified in DNA methylation maps from adult mouse tissues. Nature Genetics 45: 1198–1206. doi: 10.1038/ng.2746 23995138

20. Chu T, Handley D, Bunce K, Surti U, Hogge WA, et al. (2011) Structural and regulatory characterization of the placental epigenome at its maternal interface. PLoS One 6: e14723. doi: 10.1371/journal.pone.0014723 21373191

21. Schroeder DI, Blair JD, Lott P, Yu HO, Hong D, et al. (2013) The human placenta methylome. Proceedings of the National Academy of Sciences of the United States of America 110: 6037–6042. doi: 10.1073/pnas.1215145110 23530188

22. Ng HK, Novakovic B, Hiendleder S, Craig JM, Roberts CT, et al. (2010) Distinct patterns of gene-specific methylation in mammalian placentas: implications for placental evolution and function. Placenta 31: 259–268. doi: 10.1016/j.placenta.2010.01.009 20167366

23. Kobayashi H, Sakurai T, Imai M, Takahashi N, Fukuda A, et al. (2012) Contribution of intragenic DNA methylation in mouse gametic DNA methylomes to establish oocyte-specific heritable marks. PLoS genetics 8: e1002440. doi: 10.1371/journal.pgen.1002440 22242016

24. Wang L, Zhang J, Duan J, Gao X, Zhu W, et al. (2014) Programming and inheritance of parental DNA methylomes in mammals. Cell 157: 979–991. doi: 10.1016/j.cell.2014.04.017 24813617

25. Guo H, Zhu P, Yan L, Li R, Hu B, et al. (2014) The DNA methylation landscape of human early embryos. Nature 511: 606–610. doi: 10.1038/nature13544 25079557

26. Lee HJ, Hore TA, Reik W (2014) Reprogramming the methylome: erasing memory and creating diversity. Cell Stem Cell 14: 710–719. doi: 10.1016/j.stem.2014.05.008 24905162

27. Okae H, Chiba H, Hiura H, Hamada H, Sato A, et al. (2014) Genome-Wide Analysis of DNA Methylation Dynamics during Early Human Development. PLoS genetics 10: e1004868. doi: 10.1371/journal.pgen.1004868 25501653

28. Rugg-Gunn PJ (2012) Epigenetic features of the mouse trophoblast. Reproductive biomedicine online 25: 21–30. doi: 10.1016/j.rbmo.2012.01.012 22578826

29. Court F, Tayama C, Romanelli V, Martin-Trujillo A, Iglesias-Platas I, et al. (2014) Genome-wide parent-of-origin DNA methylation analysis reveals the intricacies of human imprinting and suggests a germline methylation-independent mechanism of establishment. Genome Research 24: 554–569. doi: 10.1101/gr.164913.113 24402520

30. Necsulea A, Soumillon M, Warnefors M, Liechti A, Daish T, et al. (2014) The evolution of lncRNA repertoires and expression patterns in tetrapods. Nature 505: 635–640. doi: 10.1038/nature12943 24463510

31. Wang X, Miller DC, Harman R, Antczak DF, Clark AG (2013) Paternally expressed genes predominate in the placenta. Proceedings of the National Academy of Sciences of the United States of America 110: 10705–10710. doi: 10.1073/pnas.1308998110 23754418

32. Sarkar AA, Nuwayhid SJ, Maynard T, Ghandchi F, Hill JT, et al. (2014) Hectd1 is required for development of the junctional zone of the placenta. Developmental Biology 392: 368–380. doi: 10.1016/j.ydbio.2014.05.007 24855001

33. Graf A, Krebs S, Zakhartchenko V, Schwalb B, Blum H, et al. (2014) Fine mapping of genome activation in bovine embryos by RNA sequencing. Proceedings of the National Academy of Sciences of the United States of America 111: 4139–4144. doi: 10.1073/pnas.1321569111 24591639

34. Reich A, Klatsky P, Carson S, Wessel G (2011) The transcriptome of a human polar body accurately reflects its sibling oocyte. The Journal of biological chemistry 286: 40743–40749. doi: 10.1074/jbc.M111.289868 21953461

35. Ziller MJ, Hansen KD, Meissner A, Aryee MJ (2015) Coverage recommendations for methylation analysis by whole-genome bisulfite sequencing. Nat Methods 12: 230–232, 231 p following 232. doi: 10.1038/nmeth.3152 25362363

36. Gaidatzis D, Burger L, Murr R, Lerch A, Dessus-Babus S, et al. (2014) DNA sequence explains seemingly disordered methylation levels in partially methylated domains of Mammalian genomes. PLoS genetics 10: e1004143. doi: 10.1371/journal.pgen.1004143 24550741

37. Hansen KD, Timp W, Bravo HC, Sabunciyan S, Langmead B, et al. (2011) Increased methylation variation in epigenetic domains across cancer types. Nature Genetics 43: 768–775. doi: 10.1038/ng.865 21706001

38. Hon GC, Hawkins RD, Caballero OL, Lo C, Lister R, et al. (2012) Global DNA hypomethylation coupled to repressive chromatin domain formation and gene silencing in breast cancer. Genome Research 22: 246–258. doi: 10.1101/gr.125872.111 22156296

39. Berman BP, Weisenberger DJ, Aman JF, Hinoue T, Ramjan Z, et al. (2012) Regions of focal DNA hypermethylation and long-range hypomethylation in colorectal cancer coincide with nuclear lamina-associated domains. Nature Genetics 44: 40–46.

40. Timp W, Bravo HC, McDonald OG, Goggins M, Umbricht C, et al. (2014) Large hypomethylated blocks as a universal defining epigenetic alteration in human solid tumors. Genome medicine 6: 61. doi: 10.1186/s13073-014-0061-y 25191524

41. Hansen KD, Sabunciyan S, Langmead B, Nagy N, Curley R, et al. (2014) Large-scale hypomethylated blocks associated with Epstein-Barr virus-induced B-cell immortalization. Genome Research 24: 177–184. doi: 10.1101/gr.157743.113 24068705

42. Marzese DM, Scolyer RA, Huynh JL, Huang SK, Hirose H, et al. (2014) Epigenome-wide DNA methylation landscape of melanoma progression to brain metastasis reveals aberrations on homeobox D cluster associated with prognosis. Human Molecular Genetics 23: 226–238. doi: 10.1093/hmg/ddt420 24014427

43. van Dijk M, Visser A, Posthuma J, Poutsma A, Oudejans CB (2012) Naturally occurring variation in trophoblast invasion as a source of novel (epigenetic) biomarkers. Frontiers in genetics 3: 22. doi: 10.3389/fgene.2012.00022 22363344

44. Rousseaux S, Debernardi A, Jacquiau B, Vitte AL, Vesin A, et al. (2013) Ectopic activation of germline and placental genes identifies aggressive metastasis-prone lung cancers. Science translational medicine 5: 186ra166.

45. Novakovic B, Saffery R (2013) Placental pseudo-malignancy from a DNA methylation perspective: unanswered questions and future directions. Frontiers in genetics 4: 285. doi: 10.3389/fgene.2013.00285 24368911

46. Wang X, Douglas KC, Vandeberg JL, Clark AG, Samollow PB (2014) Chromosome-wide profiling of X-chromosome inactivation and epigenetic states in fetal brain and placenta of the opossum, Monodelphis domestica. Genome Research 24: 70–83. doi: 10.1101/gr.161919.113 24065774

47. Douglas GC, King BF (1990) Isolation and morphologic differentiation in vitro of villous cytotrophoblast cells from rhesus monkey placenta. In Vitro Cell Dev Biol 26: 754–758. 2394671

48. Soghomonians A, Barakat AI, Thirkill TL, Blankenship TN, Douglas GC (2002) Effect of shear stress on migration and integrin expression in macaque trophoblast cells. Biochimica et biophysica acta 1589: 233–246. 12031791

49. Kumar P, Lindberg L, Thirkill TL, Ji JW, Martsching L, et al. (2012) The MUC1 extracellular domain subunit is found in nuclear speckles and associates with spliceosomes. PloS one 7: e42712. doi: 10.1371/journal.pone.0042712 22905162

50. Chen PY, Cokus SJ, Pellegrini M (2010) BS Seeker: precise mapping for bisulfite sequencing. BMC Bioinformatics 11: 203. doi: 10.1186/1471-2105-11-203 20416082

51. Karolchik D, Barber GP, Casper J, Clawson H, Cline MS, et al. (2014) The UCSC Genome Browser database: 2014 update. Nucleic Acids Research 42: D764–770. doi: 10.1093/nar/gkt1168 24270787

Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics

2015 Číslo 8

Nejčtenější v tomto čísle

Zvyšte si kvalifikaci online z pohodlí domova

Betablokátory a Ca antagonisté z jiného úhlu
nový kurz
Autoři: prof. MUDr. Michal Vrablík, Ph.D., MUDr. Petr Janský

Chronické žilní onemocnění a možnosti konzervativní léčby

Průvodce pomocnými prostředky při léčbě nemocí parodontu
Autoři: MUDr. Ladislav Korábek, CSc., MBA

Jak proměnil léčbu srdečního selhání nástup gliflozinů
Autoři: MUDr. Kristýna Kyšperská, MUDr. Jan Beneš

Autoři: doc. MUDr. Alena Šmahelová, Ph.D.

Všechny kurzy
Zapomenuté heslo

Nemáte účet?  Registrujte se

Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.


Nemáte účet?  Registrujte se