Array-based assay detects genome-wide 5-mC and 5-hmC in the brains of humans, non-human primates, and mice
© Chopra et al.; licensee BioMed Central Ltd. 2014
Received: 17 September 2013
Accepted: 5 February 2014
Published: 13 February 2014
Methylation on the fifth position of cytosine (5-mC) is an essential epigenetic mark that is linked to both normal neurodevelopment and neurological diseases. The recent identification of another modified form of cytosine, 5-hydroxymethylcytosine (5-hmC), in both stem cells and post-mitotic neurons, raises new questions as to the role of this base in mediating epigenetic effects. Genomic studies of these marks using model systems are limited, particularly with array-based tools, because the standard method of detecting DNA methylation cannot distinguish between 5-mC and 5-hmC and most methods have been developed to only survey the human genome.
We show that non-human data generated using the optimization of a widely used human DNA methylation array, designed only to detect 5-mC, reproducibly distinguishes tissue types within and between chimpanzee, rhesus, and mouse, with correlations near the human DNA level (R2 > 0.99). Genome-wide methylation analysis, using this approach, reveals 6,102 differentially methylated loci between rhesus placental and fetal tissues with pathways analysis significantly overrepresented for developmental processes. Restricting the analysis to oncogenes and tumor suppressor genes finds 76 differentially methylated loci, suggesting that rhesus placental tissue carries a cancer epigenetic signature. Similarly, adapting the assay to detect 5-hmC finds highly reproducible 5-hmC levels within human, rhesus, and mouse brain tissue that is species-specific with a hierarchical abundance among the three species (human > rhesus >> mouse). Annotation of 5-hmC with respect to gene structure reveals a significant prevalence in the 3'UTR and an association with chromatin-related ontological terms, suggesting an epigenetic feedback loop mechanism for 5-hmC.
Together, these data show that this array-based methylation assay is generalizable to all mammals for the detection of both 5-mC and 5-hmC, greatly improving the utility of mammalian model systems to study the role of epigenetics in human health, disease, and evolution.
KeywordsEpigenetics DNA methylation 5-hydroxymethylcytosine (5-hmC) Evolution
DNA methylation and histone modifications are essential epigenetic components in the establishment of the transcriptional state of eukaryotic genes throughout the genome. The best understood of these epigenetic modifications is DNA methylation, which occurs primarily at cytosines located 5' to guanosine in the CpG dinucleotide of differentiated eukaryotic nuclei . This modification is less common in CpG-rich areas, known as CpG islands, which are often located in the promoter regions of many genes and involved in transcriptional regulation. While DNA methylation is a hallmark of human cancer , the identification of a growing number of human diseases resulting from epigenetic disruptions emphasizes the importance of this mark [3–7].
The recent identification of another modified form of cytosine, 5-hydroxymethylcytosine (5-hmC), in both stem cells and post-mitotic neurons, raises new questions as to the role of this base in mediating epigenetic effects. 5-hmC is mediated by members of the ten-eleven translocation (TET) family of proteins, which each contain a C-terminal oxidase domain. Recently, a highly efficient and selective chemical approach to label and capture 5-hmC revealed the first distribution map of 5-hmC in a mammalian brain genome (mouse) and showed that 5-hmC is enriched in the gene bodies of active or highly transcribed genes . Furthermore, acquisition of 5-hmC in specific gene bodies occurs during neuronal differentiation and maturation, indicating that this mark may be distributed in a spatial and temporal manner in brain during development.
Traditionally, non-human studies of DNA methylation have relied heavily on the ‘gold standard’ for DNA methylation detection: treatment of genomic DNA with sodium bisulfite, amplification of the gene(s)/region(s) of interest, and sequencing the resulting amplicons . While this method is robust, it is subject to amplification biases and considered very laborious, allowing only for a limited number of sites to be interrogated at a time. The introduction of genome-wide arrays and ‘Next-Generation’ sequencing technologies permits genome-wide interrogation of 5-methylcytosine (5-mC) but these methodologies either focus solely on human studies (arrays), do not provide base resolution (enrichment techniques), or are cost prohibitive (whole genome sequencing). In addition, these methods could not distinguish between 5-mC and 5-hmC; however, numerous sequencing-based methods have recently emerged that now allow for the quantitative discrimination of the two modified bases [10, 11]. In this study we combined these recent 5-hmC detection methods with standard array-based methods to develop a novel use for the HumanMethylation450 BeadChip from Illumina to measure both 5-mC and 5-hmC in mammalian model systems.
Non-human primates (e.g. rhesus monkeys) and humans have structural and functional similarities in the neural circuits mediating neuropsychiatric disease. Access to relevant brain regions associated with neuropsychiatric disorders from a monkey model improves our ability to translate research findings to humans and provides opportunities to test hypotheses that cannot be tested in humans. However, genome-wide molecular approaches, such as array based assays, in non-human genomes have not matched resources aimed at the human genome. Here, we optimized an established human array that was originally designed to detect 5-mC levels at nearly half a million loci throughout the human genome and show that it is also capable of measuring genome-wide DNA methylation levels in other species, including non-human primates and mice.
To overcome the paucity of molecular tools for DNA methylation detection in non-human primates, we optimized a highly reproducible and widely used human array [12–16] to measure genome-wide DNA methylation levels in rhesus macaque (Macaca mulatta). The human array is designed to provide a quantitative measure of 5-methylcytocine (DNA methylation; 5-mC) at 485,512 CpG dinucleotides located proximal to nearly all human RefSeq genes (N = 21,231). We interrogated each fifty nucleotide probe sequence on the array to determine the fraction that had either an exact match or up to four mismatches to the rhesus genome and found that approximately 7% and 30%, respectively, of the array probes met this criteria, (rhesus-competent probes; exact match: N = 35,901; mismatch: N = 154,030; Additional file 1; see Methods). We hypothesized that both sets of these rhesus-competent probes should reproducibly provide a quantitative measure of DNA methylation at the majority (exact match: > 55%, N = 11,804; mismatch > 90%, N = 19,362) of the human annotated RefSeq genes in the rhesus genome. Indeed filtering the array probes to only the rhesus-competent probes reduces the median number of CpG dinucleotides interrogated per gene from fifteen (human) to two or six (rhesus exact match or mismatch, respectively; Additional file 2A).
Genome-wide DNA methylation analysis in rhesus macaque
To further explore the potential of this assay, we used the rhesus-competent probes to conduct a genome-wide regression analysis of methylation levels between placental and fetal tissue and identified 6,102 differentially methylated loci (Bonferroni <0.05; Additional file 4). To statistically validate our methods and adjust for possible failure of asymptotic assumptions due to the small sample size, we performed one hundred permutations of the data and found that there was a high correlation between the asymptotic P-values and the permuted P-values (Additional file 2B and Additional file 4; Methods). Annotation of the differentially methylated loci using gene ontology and Kyoto encyclopedia of genes and genomes (KEGG) pathways revealed biologically relevant terms that were significantly overrepresented for developmental processes (Additional file 5). However, the KEGG analysis did not find an overrepresentation of cancer related terms, suggesting that normal tissue differentiation is the primary purpose of the different methylation marks. Nonetheless, we restricted the probes used in the regression analysis to those annotated to oncogenes and tumor suppressor genes and identified cancer-specific differentially methylated loci (N = 76; Additional file 6), suggesting placental tissue carries a cancer epigenetic signature.
Species-specific methylation marks
These findings led us to test the hypothesis that other species-competent probe sets could be identified and used to study genome-wide DNA methylation. Thus we conducted a similar probe analysis to identify competent probes for the mouse (N = 9,734; Additional file 1) and the chimpanzee (N = 360,491; Additional file 1). Using the mouse-competent probes, we tested the reproducibility of this array and found that each mouse’s replicates (N = 3) were highly correlated to each other, with R2-values similar to those found among human and rhesus replicates for this assay (R2 > 0.99; Figure 1). In addition, using each set of species-competent probes (human, chimpanzee, rhesus, and mouse), we compared the consistency of this array between individuals of the same species (N = 3) and showed that there is a similarly high consistency among individuals of the same species (mean R2 > 0.90; Additional file 7), and that a set of common probes (N = 8,189) can clearly distinguish the species from one another (Additional file 8A). Furthermore, we investigated the gene structure distribution of the species-specific probes (i.e. monkey and mouse; Additional file 9). To test the significance of this probe distribution in monkey and mouse we compared the proportion of CpG loci that reside in each structure to the total number of CpG loci interrogated on the DNA methylation panel via permutation testing. We found that the optimization in both monkey and mouse is significantly biased (P < 0.0001) to retain more conserved gene structures (i.e. 5’ UTR, 1st exon, and gene body regions). This finding is consistent with the probe selection method used to optimize the assay, which requires that the species-specific probes be a near match to the human array (allowing for only a four base mismatch), meaning that conserved genomic structures are going to be preferentially retained and that the analyses of non-human genomes are going to be limited to these regions. Together, these data suggest that this optimization method is generalizable to any sequenced mammalian species, making genome-wide DNA methylation detection possible for a much broader group of species.
Novel utility of the array
This finding prompted us to investigate whether the array could accurately detect 5-hmC in other mammalian genomes and we found that it provides a consistent quantitative measure of 5-hmC throughout both the monkey and mouse brain genomes (Additional file 8C-D); again suggesting that optimizing this assay is generalizable to any sequenced mammalian species. Comparison with the human 5-hmC profiles also revealed several interesting biological findings, including that 5-hmC profiles clearly distinguish between species (N = 8,189 probes; Additional file 8B) and that the primate median 5-hmC levels are highly correlated (Figure 4C). Thus, despite the fact that there are species-distinguishing 5-hmC marks, the majority of marks have a consistent level of 5-hmC, at least in primates (Figure 4D), suggesting that there is likely to be a significant overlap of primate genes with abundant 5-hmC levels. Consistent with this rationale, we found that primates indeed have a significant overlap of probes with > 20% 5-hmC levels (Figure 4E). In addition, primates have significantly more probes with > 20% 5-hmC than mice (53% (human), 42% (monkey), and 17% (mouse); Figure 4D), suggesting that this mark may play a larger role in primates than mice. Together, these data show that this array-based assay has a broad utility in DNA methylation studies of both 5-mC and 5-hmC.
Hydroxymethylation and genomic context
Locations of DNA methylation relative to a gene refer to regions that include the distance to the gene TSS (TSS1500 (i.e. within 1500 bp of the TSS) and TSS200), 5’UTR, 1st exon, body, and 3’UTR. Previous studies of 5-hmC levels in mouse hippocampal and cerebellum tissue revealed that 5-hmC is enriched in gene bodies . Here we found a similar enrichment in the gene body of all species (P < 0.0001) but our data also revealed that the 3’ UTR has more 5-hmC than would be expected by chance alone (P < 0.0001; Figure 5F-K and Additional file 12). Since we also found a significant underrepresentation in all other genic regions (P < 0.0001), these data suggest that a non-stochastic mechanism may prevent 5-hmC from the 5’ end of genes.
X-linked 5-hmC loci have functions in chromatin modifications and cell death
In humans the X chromosome was deficient of loci with high levels of 5-hmC. Annotation of X-linked CpG loci with > 40% 5-hmC in primates or > 20% in mice using gene ontology revealed several biologically relevant ontologies that were significantly overrepresented among all species. These findings predominately included terms associated with epigenetic functions including covalent chromatin modifications and histone acetylation. (Additional file 5; see Methods). Together, these data suggest that 5-hmC may play a role in an epigenetic feedback loop. Interestingly, annotation of the human X-linked CpG loci with > 60% 5-hmC using gene ontology revealed several biologically relevant ontologies that were significantly overrepresented for cell death and apoptotic processes, suggesting that 5-hmC may target a specific subset of genes on the X chromosome. Notably, gene ontological analysis of the entire array (N = 21,231 genes) also finds significant terms, suggesting that the full array may have a bias toward a subset of genes and that conclusions from ontological analysis should be taken with caution. Clearly further studies are needed to determine the diverse functional contributions of 5-hmC in the mammalian brain.
This study optimizes an established human array-based DNA methylation (5-mC) assay to allow for accurate detection of 5-mC in non-human primates and mice and, perhaps more importantly, for the detection of genome-wide 5-hmC in mammals. Optimization of this assay dramatically improves its utility because now it can be used to simultaneously investigate genome-wide 5-mC and 5-hmC in any sequenced mammalian genome. Moreover, since the analysis of data generated by this array is well established, this assay will reduce the cost and time of DNA methylation studies, allowing for a greater utilization of mammalian model systems to investigate the role of epigenetics in mammalian biological processes.
This optimization achieved high reproducibility among 5-mC data generated for chimpanzee, rhesus, and mouse by utilizing an approach that relied on filtering the data down to the probes with a near match (allowing for up to 4 mismatches) between the human genome and the species of interest. This filter for species-specific probes still provides data for > 90% of the human annotated genes in the monkey and nearly 25% for the mouse. Finding that the chimpanzee analysis is able to retain far more probes (N = 360,491) is consistent with known evolutionary similarity with humans but also may reflect the higher quality reference sequence of the human and chimpanzee genomes compared to the rhesus genome. Recently, a similar approach identified 13,715 mouse-competent probes on the human array , which contains nearly 75% of the mouse-competent probes reported here. Together these independent findings strongly support the use of this array with other species. The only other array-based method available for 5-mC detection in mice (Mouse CpG Island Microarray; AgilentTechnologies) relies heavily on the efficiency/reproducibility of restriction enzymes and the data is only reports 5-mC methylation at CpG islands, which our data here clearly shows is lacking the majority of the variable DNA methylation information (Figure 5).
Having the ability to conduct an epigenome study in wide variety of sequenced mammals represents a significant new methodology. Here we compared the output of the optimized genome-wide array to that of a previously reported locus-specific study by Chiu et al., both investigating the epigenetic differences between monkey placenta and fetal tissue. Clearly the sheer magnitude of data from the array is not comparable between studies and neither are the conclusions. Notably, in 2007, when the Chui et al. study was conducted, the HumanMethylation450 BeadChips had not been developed; meaning Chiu and colleagues were using the molecular tools available at that time. This fact notwithstanding, while Chui et al. only found methylation differences in one of the nine interrogated genes (RASSF1), the much more comprehensive array approach developed here found differences in four out of the nine genes. In addition, the array identified differences in a total of 2,825 genes throughout the genome, including hypermethylated in cancer 1 (HIC1) and TIMP metallopeptidase inhibitor 3 (TIMP3), both of which were previously shown, using a locus-specific approach, to not be differentially methylated in normal placental tissue . Together, these data highlight the dramatic advantage to using the array-based approach and suggest that studies of placental methylation need to be revisited.
It was previously shown that β-glucosyltransferase could be used to modify an existing array-based protocol for the detection of 5-hmC. However, in that case the array was not specifically designed for methylation detection, the method relied on the efficiency/reproducibility of restriction enzymes, and the data was biased toward CpG islands. In contrast, using a modified protocol and the HumanMethylation450 BeadChip largely overcomes these limitations and produces a more reliable unbiased profile of 5-hmC throughout the genome, providing a new resource for epigenetic studies.
The density of 5-hmC probes in all three species profiles suggests that the bulk of 5-hmC levels are from 20% - 60%, especially for primates (Figure 4D). While it was expected that the majority of probes would have a low level of 5-hmC (< 40%), we did not anticipate that, in the brain, nearly 50% of the primate loci would be > 20% hydroxymethylated and nearly 20% of human loci would be > 50% hydroxymethylated. These data suggest that 5-hmC is much more prevalent than previously reported and that it may play a more significant role in primates than mice. However, the gene ontological analysis suggests that this role is in a feedback loop for epigenetic processes such as histone lysine methylation and that these functions are common to all mammals.
In the brain tissues of all three species we found the distribution of 5-hmC with relation to gene structure to be consistent with previous reports that found 5-hmC is especially prevalent in the regions adjacent to CpG islands (in the shelves), indicating that the highest levels of 5-hmC are at regions of the genome with the lowest CpG density. Also consistent with previous work was the finding that 5-hmC is more prevalent in gene bodies compared to the 5’ end of a gene; however, we find that the 3’ UTR is particularly enriched for this mark, a location that may facilitate a proposed role in a epigenetic feedback loop. At the chromosomal level 5-hmC was reported elsewhere to be deficient from the X chromosome in mouse hippocampal and cerebellum tissue. In contrast, we didn’t find a deficiency on any mouse chromosomes but we did find it is significantly deficient on the human X and approaching significance on the monkey X. Although the synteny between human and mouse X chromosomes is high, the smaller number of loci from the mouse may have resulted in a higher genome-wide variance that precluded the mouse X data to reach significance.
As we learned with genome-wide association studies, successful epigenome-wide association studies of complex phenotypes are going to require the analysis of a very large number of individuals (N > 1000). The development of this assay improves our ability to quickly and reproducibly analyze the dynamic role of DNA methylation in a large cohort of individuals. This assay also augments the use of model systems for the study of epigenetic mechanisms and phenomena that are difficult to study in humans, including other recently identified modifications to methylcytosine, namely 5-formylcytosine and 5-carboxylcytosine. In addition, data from model systems will make it quite tempting to draw evolutionary conclusions between the different mammalian species and the gene ontology analysis shown here suggests that there are more similarities among species than differences. In short, this optimized assay holds great promise to improve our understanding of DNA methylation and its contributions to health, disease, and evolution.
Here we have developed an array-based assay to detect 5-mC and 5-hmC in mammalian genomes and show that it is generalizable to all mammals. Using this assay we find that placental tissue carries a cancer epigenetic signature, that 5-hmC has a hierarchical abundance among mammals, and that X-linked loci and the 3’UTR of genes have a heightened prevalence for 5-hmC. In short, this optimized assay greatly improves the utility of mammalian model systems for studying the role of epigenetics in human health, disease, and evolution.
Tissue acquisition and DNA extraction
Except for the mouse brain tissues (UW IACUC protocol #M02529), all other tissues were supplied by other research groups, including monkey skin and brain tissue (Dr. Kalin’s group; UW IACUC protocol #G00181), monkey placental and fetal tissue (Dr. Golos’ group; UW IACUC protocol #G00245), and human brain tissue (Emory Neuroscience NINDS Core Facilities; approved by the Emory Internal Review Board on August 18, 2013 and given the approval number CR2_IRB00045782). All groups follow animal housing and experimental procedures that are in accordance with institutional guidelines. Once acquired, approximately thirty milligrams of tissue was homogenized with glass beads (Sigma) and DNA extraction was performed using AllPrep DNA/RNA mini kit (Qiagen).
DNA methylation profiling
Genomic DNA was split into the following two aliquots: a sodium bisulfite-only (SB-only) aliquot; and a Tet1 treated followed by SB treated aliquot (TAB). The Tet1 treatment was conducted using previously described methods . Briefly, one microgram (μg) of genomic DNA was spiked with four nanograms (ng) of synthetically methylated control DNAs (N = 3). Together, these DNAs were sonicated to 400 basepairs (bp) using a Diagenode Bioruptor Plus (Diagenode) and concentrated to 19 μl (Amicon Ultacel 100 k, Millipore). To mark 5-hmC, the concentrated DNA was then glucosylated with T4-β-glucosyltransferase (1 μM; Wisegene) and UDP-glucose (200 μM, Wisegene) at 37°C for 1 hour, and purified using the QIAquick Nucleotide Removal Kit (Qiagen). To mark 5-mC, five hundred ng of the glucosylated DNA was oxidized in the presence of recombinant mouse Tet1 (0.48 μg/μl; Wisegene) at 37°C for 1 hour followed by enzyme inactivation with protease K at 50°C for 1 hour. The oxidized DNA was twice purified (Micro Bio-Spin 30, Bio-Rad; QIAquick PCR Purification Kit, Qiagen) and eluted in 20 μl of water. SB-treatment of both aliquots: Five hundred nanograms of genomic DNA or Tet-treated DNA was sodium bisulfite–treated for cytosine (C) to thymine (T) conversion using the EZ DNA Methylation-Gold kit (Zymo Research). The converted DNA was purified and prepped for analysis on the Illumina HumanMethylation450 BeadChips following the manufacturer’s guidelines. Briefly, converted DNA was amplified, fragmented, and hybridized to the HumanMethylation450 pool of allele-differentiating oligonucleotides. After a series of extension, ligation, and cleanup reactions, the DNA was labeled with a fluorescent dye. The labeled DNA was then scanned using an Illumina BeadArray Reader or iScan. Image analysis and signal determination were performed using the GenomeStudio software, Methylation Module (Illumina).
Interpretation and QC of DNA methylation data
CpG DNA methylation data were interpreted using GenomeStudio to quantify methylated (M) and unmethylated (U) signal intensities for genomic DNA. The signals were quantile normalized, and overall methylation levels (β) were calculated as the ratio of methylated to total signal [i.e. β = M/(M + U + 100)], where β ranges from 0 (unmethylated) to 1 (methylated). For human samples, quality control of data resulted in removal of samples with aberrantly low signal intensity (mean < 2000) or with fewer than 99% of CpG loci detected, where a given locus was deemed not detected if the detection P-value was > 0.01 (detection P-value provided by Genome Studio and calculated relative to background signal). For rhesus, chimp and mouse samples, the data for only the ‘competent’ probe-set (ref. section ‘Assay optimization’) for each species was taken. Furthermore, any probes having more than 20% of samples with detected P-values > 0.01 were discarded from the analysis. Assay controls were inspected to remove samples with poor bisulfite conversion, staining, extension (single nucleotide extension assay), hybridization, or specificity. Outliers identified by hierarchical clustering and/or dissimilarity matrices were removed. Additionally, one control DNA replicate was run on each BeadChip to assess overall assay reproducibility. Methylation profiles of the control DNA correlated well, with an average Pearson correlation coefficient (R) of 0.990 between human replicates.
Analysis of cancer-associated CpG loci
To analyze DNA methylation differences associated with cancer, we fit a separate regression for each CpG site. Although samples were randomly distributed across BeadChips and experiments with respect to disease, BeadChip was also included as a random effect covariate in all analyses to account for potential batch effects. The package “nlme” in R (Cran) was used for the mixed effect model. Fixed effects included intensity in the model. To correct for multiple hypothesis testing, we applied a Benjamini-Hochberg False Discovery Rate (FDR) correction using the R function “p.adjust,” but to avoid false positives due to the small sample size, we used conservative Bonferroni adjustment for our ultimate determination of significance.
Permutation analyses for rhesus differentially methylated loci
All permutation analyses were conducted in R using the same linear model as the actual analysis, where BeadChip was treated as a mixed-effects covariate, but in each permutation the tissue type of the sample was randomly reassigned. In total, 100 permutations were conducted. Permutation P-values for each CpG locus were calculated by assessing the number of times each locus was more significantly associated with fetal tissue in the 100 permuted data sets than the actual.
Analysis of genomic location
Loci were mapped to chromosome location using the Illumina annotation files for the HumanMethylation450 panel based on the NCBI Human Genome (build 37). Over- and underrepresentation of CpG loci with > 20% median 5-hmC were assessed using both Fisher’s exact test and permutation analysis, which corresponded well (data not shown). Permutation P-values were determined as described previously , by comparing the proportion of >20% 5-hmC loci on a chromosome (termed “actual proportion”) to the proportions obtained using the permuted loci (termed “permuted proportions”). To correct for multiple hypotheses (i.e. 24 chromosomes and both over and underrepresentation) the actual proportion of each chromosome was compared to the permuted proportions of any chromosome for each permutation. Thus the permutation P-value is the number of times a permuted proportion is more significant than the actual proportion, divided by the number of permutations (10,000).
We have submitted the data generated from the monkey, chimp, human and mouse samples for this study to the Gene Expression Omnibus (GEO), which can be found under the Gene Series: GSE49177.
The authors would like to thank Dr. Peng Jin for critical comments on the manuscript and Julie Mowrey and Brian Lynch for technical assistance. This work was supported in part by the University of Wisconsin-Madison department of Psychiatry (RSA), the Alzheimer’s Disease Research Center grant number P50 AG025688 (Emory University human brain tissue), the Conti Center, and National Institutes of Health (NIH) grants MH081884, MH084051, and MH046729, all to Ned H. Kalin and NIH grant MH089606 to S.T.W. Emory University’s research IT service center and its high performance computer cluster also supported this research. This research project was supported in part by the Neuropathology core of the Emory Neuroscience NINDS Core Facilities grant, P30NS055077.
- Jaenisch R, Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003, 33 (Suppl): 245-254.PubMedView ArticleGoogle Scholar
- Feinberg AP, Vogelstein B: Hypomethylation distinguishes genes of some human cancers from their normal counterparts. Nature. 1983, 301: 89-92. 10.1038/301089a0.PubMedView ArticleGoogle Scholar
- Moulton T, Crenshaw T, Hao Y, Moosikasuwan J, Lin N, Dembitzer F, Hensle T, Weiss L, McMorrow L, Loew T, et al: Epigenetic lesions at the H19 locus in Wilms’ tumour patients. Nat Genet. 1994, 7: 440-447. 10.1038/ng0794-440.PubMedView ArticleGoogle Scholar
- Steenman MJ, Rainier S, Dobry CJ, Grundy P, Horon IL, Feinberg AP: Loss of imprinting of IGF2 is linked to reduced expression and abnormal methylation of H19 in Wilms’ tumour. Nat Genet. 1994, 7: 433-439. 10.1038/ng0794-433.PubMedView ArticleGoogle Scholar
- Sutcliffe JS, Nakao M, Christian S, Orstavik KH, Tommerup N, Ledbetter DH, Beaudet AL: Deletions of a differentially methylated CpG island at the SNRPN gene define a putative imprinting control region. Nat Genet. 1994, 8: 52-58. 10.1038/ng0994-52.PubMedView ArticleGoogle Scholar
- Arima T, Drewell RA, Arney KL, Inoue J, Makita Y, Hata A, Oshimura M, Wake N, Surani MA: A conserved imprinting control region at the HYMAI/ZAC domain is implicated in transient neonatal diabetes mellitus. Hum Mol Genet. 2001, 10: 1475-1483. 10.1093/hmg/10.14.1475.PubMedView ArticleGoogle Scholar
- Orr HT, Zoghbi HY: Trinucleotide repeat disorders. Annu Rev Neurosci. 2007, 30: 575-621. 10.1146/annurev.neuro.29.051605.113042.PubMedView ArticleGoogle Scholar
- Szulwach KE, Li X, Li Y, Song CX, Wu H, Dai Q, Irier H, Upadhyay AK, Gearing M, Levey AI, et al: 5-hmC-mediated epigenetic dynamics during postnatal neurodevelopment and aging. Nat Neurosci. 2011, 14: 1607-1616. 10.1038/nn.2959.PubMed CentralPubMedView ArticleGoogle Scholar
- Frommer M, McDonald LE, Millar DS, Collis CM, Watt F, Grigg GW, Molloy PL, Paul CL: A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. Proc Natl Acad Sci USA. 1992, 89: 1827-1831. 10.1073/pnas.89.5.1827.PubMed CentralPubMedView ArticleGoogle Scholar
- Booth MJ, Branco MR, Ficz G, Oxley D, Krueger F, Reik W, Balasubramanian S: Quantitative sequencing of 5-methylcytosine and 5-hydroxymethylcytosine at single-base resolution. Science. 2012, 336: 934-937. 10.1126/science.1220671.PubMedView ArticleGoogle Scholar
- Yu M, Hon GC, Szulwach KE, Song CX, Zhang L, Kim A, Li X, Dai Q, Shen Y, Park B, et al: Base-resolution analysis of 5-hydroxymethylcytosine in the mammalian genome. Cell. 2012, 149: 1368-1380. 10.1016/j.cell.2012.04.027.PubMed CentralPubMedView ArticleGoogle Scholar
- Bork S, Pfister S, Witt H, Horn P, Korn B, Ho AD, Wagner W: DNA methylation pattern changes upon long-term culture and aging of human mesenchymal stromal cells. Aging Cell. 2010, 9: 54-63. 10.1111/j.1474-9726.2009.00535.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Weisenberger DJ, Shen H, Campan M, Noushmehr H, Bell CG, Maxwell AP, et al: Age-dependent DNA methylation of genes that are suppressed in stem cells is a hallmark of cancer. Genome Res. 2010, 20: 440-446. 10.1101/gr.103606.109.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen YA, Choufani S, Ferreira JC, Grafodatskaya D, Butcher DT, Weksberg R: Sequence overlap between autosomal and sex-linked probes on the Illumina HumanMethylation27 microarray. Genomics. 2011, 97: 214-222. 10.1016/j.ygeno.2010.12.004.PubMedView ArticleGoogle Scholar
- Koch CM, Suschek CV, Lin Q, Bork S, Goergens M, Joussen S, Pallua N, Ho AD, Zenke M, Wagner W: Specific age-associated DNA methylation changes in human dermal fibroblasts. PLoS One. 2011, 6: e16679-10.1371/journal.pone.0016679.PubMed CentralPubMedView ArticleGoogle Scholar
- Alisch RS, Barwick BG, Chopra P, Myrick LK, Satten GA, Conneely KN, Warren ST: Age-associated DNA methylation in pediatric populations. Genome Res. 2012, 22: 623-632. 10.1101/gr.125187.111.PubMed CentralPubMedView ArticleGoogle Scholar
- Nazor KL, Altun G, Lynch C, Tran H, Harness JV, Slavin I, Garitaonandia I, Muller FJ, Wang YC, Boscolo FS, et al: Recurrent variations in DNA methylation in human pluripotent stem cells and their differentiated derivatives. Cell Stem Cell. 2012, 10: 620-634. 10.1016/j.stem.2012.02.013.PubMed CentralPubMedView ArticleGoogle Scholar
- Bischof P, Meisser A, Campana A: Biochemistry and molecular biology of trophoblast invasion. Ann N Y Acad Sci. 2001, 943: 157-162.PubMedView ArticleGoogle Scholar
- Janneau JL, Maldonado-Estrada J, Tachdjian G, Miran I, Motte N, Saulnier P, Sabourin JC, Cote JF, Simon B, Frydman R, et al: Transcriptional expression of genes involved in cell invasion and migration by normal and tumoral trophoblast cells. J Clin Endocrinol Metab. 2002, 87: 5336-5339. 10.1210/jc.2002-021093.PubMedView ArticleGoogle Scholar
- Zhou Y, Genbacev O, Fisher SJ: The human placenta remodels the uterus by using a combination of molecules that govern vasculogenesis or leukocyte extravasation. Ann N Y Acad Sci. 2003, 995: 73-83. 10.1111/j.1749-6632.2003.tb03211.x.PubMedView ArticleGoogle Scholar
- Ohlsson R: Growth factors, protooncogenes and human placental development. Cell Differ Dev. 1989, 28: 1-15. 10.1016/0922-3371(89)90018-X.PubMedView ArticleGoogle Scholar
- Strickland S, Richards WG: Invasion of the trophoblasts. Cell. 1992, 71: 355-357. 10.1016/0092-8674(92)90503-5.PubMedView ArticleGoogle Scholar
- Chiu RW, Chim SS, Wong IH, Wong CS, Lee WS, To KF, Tong JH, Yuen RK, Shum AS, Chan JK, et al: Hypermethylation of RASSF1A in human and rhesus placentas. Am J Pathol. 2007, 170: 941-950. 10.2353/ajpath.2007.060641.PubMed CentralPubMedView ArticleGoogle Scholar
- Ito S, D’Alessio AC, Taranova OV, Hong K, Sowers LC, Zhang Y: Role of Tet proteins in 5mC to 5hmC conversion, ES-cell self-renewal and inner cell mass specification. Nature. 2010, 466: 1129-1133. 10.1038/nature09303.PubMed CentralPubMedView ArticleGoogle Scholar
- Ficz G, Branco MR, Seisenberger S, Santos F, Krueger F, Hore TA, Marques CJ, Andrews S, Reik W: Dynamic regulation of 5-hydroxymethylcytosine in mouse ES cells and during differentiation. Nature. 2011, 473: 398-402. 10.1038/nature10008.PubMedView ArticleGoogle Scholar
- Williams K, Christensen J, Pedersen MT, Johansen JV, Cloos PA, Rappsilber J, Helin K: TET1 and hydroxymethylcytosine in transcription and DNA methylation fidelity. Nature. 2011, 473: 343-348. 10.1038/nature10066.PubMed CentralPubMedView ArticleGoogle Scholar
- Wu H, Zhang Y: Mechanisms and functions of Tet protein-mediated 5-methylcytosine oxidation. Genes Dev. 2011, 25: 2436-2452. 10.1101/gad.179184.111.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang T, Wu H, Li Y, Szulwach KE, Lin L, Li X, Chen IP, Goldlust IS, Chamberlain SJ, Dodd A, et al: Subtelomeric hotspots of aberrant 5-hydroxymethylcytosine-mediated epigenetic modifications during reprogramming to pluripotency. Nat Cell Biol. 2013, 15: 700-711. 10.1038/ncb2748.PubMed CentralPubMedView ArticleGoogle Scholar
- Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, et al: The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009, 41: 178-186. 10.1038/ng.298.PubMed CentralPubMedView ArticleGoogle Scholar
- Szulwach KE, Li X, Li Y, Song CX, Han JW, Kim S, Namburi S, Hermetz K, Kim JJ, Rudd MK, et al: Integrating 5-hydroxymethylcytosine into the epigenomic landscape of human embryonic stem cells. PLoS Genet. 2011, 7: e1002154-10.1371/journal.pgen.1002154.PubMed CentralPubMedView ArticleGoogle Scholar
- Pastor WA, Pape UJ, Huang Y, Henderson HR, Lister R, Ko M, McLoughlin EM, Brudno Y, Mahapatra S, Kapranov P, et al: Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells. Nature. 2011, 473: 394-397. 10.1038/nature10102.PubMed CentralPubMedView ArticleGoogle Scholar
- Wong NC, Ng J, Hall NE, Lunke S, Salmanidis M, Brumatti G, Ekert PG, Craig JM, Saffery R: Exploring the utility of human DNA methylation arrays for profiling mouse genomic DNA. Genomics. 2013, 102: 38-46. 10.1016/j.ygeno.2013.04.014.PubMedView ArticleGoogle Scholar
- Xue WC, Chan KY, Feng HC, Chiu PM, Ngan HY, Tsao SW, Cheung AN: Promoter hypermethylation of multiple genes in hydatidiform mole and choriocarcinoma. J Mol Diagn. 2004, 6: 326-334. 10.1016/S1525-1578(10)60528-4.PubMed CentralPubMedView ArticleGoogle Scholar
- Khare T, Pai S, Koncevicius K, Pal M, Kriukiene E, Liutkeviciute Z, Irimia M, Jia P, Ptak C, Xia M, et al: 5-hmC in the brain is abundant in synaptic genes and shows differences at the exon-intron boundary. Nat Struct Mol Biol. 2012, 19: 1037-1043. 10.1038/nsmb.2372.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.