Examination of the specificity of DNA methylation profiling techniques towards 5-methylcytosine and 5-hydroxymethylcytosine

DNA cytosine-5 methylation is a well-studied epigenetic pathway implicated in gene expression control and disease pathogenesis. Different technologies have been developed to examine the distribution of 5-methylcytosine (5mC) in specific sequences of the genome. Recently, substantial amounts of 5-hydroxymethylcytosine (5hmC), most likely derived from enzymatic oxidation of 5mC by TET1, have been detected in certain mammalian tissues. Here, we have examined the ability of several commonly used DNA methylation profiling methods to distinguish between 5mC and 5hmC. We show that techniques based on sodium bisulfite treatment of DNA are incapable of distinguishing between the two modified bases. In contrast, techniques based on immunoprecipitation with anti-5mC antibody (methylated DNA immunoprecipitation, MeDIP) or those based on proteins that bind to methylated CpG sequences (e.g. methylated-CpG island recovery assay, MIRA) do not detect 5hmC and are specific for 5mC unless both modified bases occur in the same DNA fragment. We also report that several methyl-CpG binding proteins including MBD1, MBD2 and MBD4 do not bind to sequences containing 5hmC. Selective mapping of 5hmC will require the development of unique tools for the detection of this modified base.     

MIRA-SNuPE,a quantitative,multiplex method for measuring allele-specific DNA methylation

5-methyl-C (5mC) and 5-hydroxymethyl-C (5hmC) are epigenetic marks with well-known and putative roles in gene regulation, respectively. These two DNA covalent modifications cannot be distinguished by bisulfite sequencing or restriction digestion, the standard methods of 5mC detection. The methylated CpG island recovery assay (MIRA), however, specifically detects 5mC but not 5hmC. We further developed MIRA for the analysis of allele-specific CpG methylation at differentially methylated regions (DMRs) of imprinted genes. MIRA specifically distinguished between the parental alleles by capturing the paternally methylated H19/Igf2 DMR and maternally methylated KvDMR1 in mouse embryo fibroblasts (MEFs) carrying paternal and maternal duplication of mouse distal Chr7, respectively. MIRA in combination with multiplex single nucleotide primer extension (SNuPE) assays specifically captured the methylated parental allele from normal cells at a set of maternally and paternally methylated DMRs. The assay correctly recognized aberrant biallelic methylation in a case of loss of imprinting. The MIRA-SNuPE assays revealed that placenta exhibited less DNA methylation bias at DMRs compared to yolk sac, amnion, brain, heart, kidney, liver and muscle. This method should be useful for the analysis of allele-specific methylation events related to genomic imprinting, X chromosome inactivation and for verifying and screening haplotype-associated methylation differences in the human population.Key words: epigenetics, imprinting, DMR, MIRA, MBD, DNA methylation, SNuPE     

表观遗传学新标记——5-羟甲基胞嘧啶检测方法的研究进展

被称为"第六种碱基"的5-羟甲基胞嘧啶(5-hydroxymethylcytosine, 5hmC),广泛分布于多种哺乳动物的组织和细胞中,与胚胎发育,神经系统功能以及肿瘤研究高度相关.与5-甲基胞嘧啶(5-methylcytosine, 5mC)相比,5hmC在组织中含量更低,难以精确的检测.随着研究的深入,5hmC参与的重要生物学作用逐渐被人们发现,同时也促使着5hmC的检测和定量方法不断发展.为了区分5hmC与其他胞嘧啶衍生物,很多利用化学或者酶学修饰实现靶向检测或非靶向富集5hmC的方法应运而生.因此,选择并发展灵敏,准确,可靠的5hmC检测技术对于表观遗传研究至关重要.本文重点综述了近年来发展起来的5hmC检测和测序技术,通过比较分析各种方法的优缺点,为研究人员选择特定合适的方法开展相关研究提供重要的参考.     

Effects of particulate matter exposure on blood 5-hydroxymethylation: results from the Beijing truck driver air pollution study

Previous studies have reported epigenetic changes induced by environmental exposures. However, previous investigations did not distinguish 5-methylcytosine (5mC) from a similar oxidative form with opposite functions, 5-hydroxymethylcytosine (5hmC). Here, we measured blood DNA global 5mC and 5hmC by ELISA and used adjusted mixed-effects regression models to evaluate the effects of ambient PM₁₀ and personal PM_2.5 and its elemental components—black carbon (BC), aluminum (Al), calcium (Ca), potassium (K), iron (Fe), sulfur (S), silicon (Si), titanium (Ti), and zinc (Zn)—on blood global 5mC and 5hmC levels. The study was conducted in 60 truck drivers and 60 office workers in Beijing, China from The Beijing Truck Driver Air Pollution Study at 2 exams separated by one to 2 weeks. Blood 5hmC level (0.08%) was ∼83-fold lower than 5mC (6.61%). An inter-quartile range (IQR) increase in same-day PM₁₀ was associated with increases in 5hmC of 26.1% in office workers (P = 0.004), 20.2% in truck drivers (P = 0.014), and 21.9% in all participants combined (P < 0.001). PM₁₀ effects on 5hmC were increasingly stronger when averaged over 4, 7, and 14 d preceding assessment (up to 132.6% for the 14-d average in all participants, P < 0.001). PM₁₀ effects were also significant after controlling for multiple testing (family-wise error rate; FWER < 0.05). 5hmC was not correlated with personal measures of PM_2.5 and elemental components (FWER > 0.05). 5mC showed no correlations with PM₁₀, PM_2.5, and elemental components measures (FWER > 0.05). Our study suggests that exposure to ambient PM₁₀ affects 5hmC over time, but not 5mC. This finding demonstrates the need to differentiate 5hmC and 5mC in environmental studies of DNA methylation.     

Altering TET dioxygenase levels within physiological range affects DNA methylation dynamics of HEK293 cells

The TET family of dioxygenases (TET1/2/3) can convert 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) and has been shown to be involved in active and passive DNA demethylation. Here, we demonstrate that altering TET dioxygenase levels within physiological range can affect DNA methylation dynamics of HEK293 cells. Overexpression of TET1 increased global 5hmC levels and was accompanied by mild DNA demethylation of promoters, gene bodies and CpG islands. Conversely, the simultaneous knockdown of TET1, TET2, and TET3 led to decreased global 5hmC levels and mild DNA hypermethylation of above-mentioned regions. The methylation changes observed in the overexpression and knockdown studies were mostly non-reciprocal and occurred with different preference depending on endogenous methylation and gene expression levels. Single-nucleotide 5hmC profiling performed on a genome-wide scale revealed that TET1 overexpression induced 5mC oxidation without a distribution bias among genetic elements and structures. Detailed analysis showed that this oxidation was related to endogenous 5hmC levels. In addition, our results support the notion that the effects of TET1 overexpression on gene expression are generally unrelated to its catalytic activity.     

Loss of 5-hydroxymethylcytosine is linked to gene body hypermethylation in kidney cancer

Both 5-methylcytosine (5mC) and its oxidized form 5-hydroxymethylcytosine (5hmC) have been proposed to be involved in tumorigenesis. Because the readout of the broadly used 5mC mapping method, bisulfite sequencing (BS-seq), is the sum of 5mC and 5hmC levels, the 5mC/5hmC patterns and relationship of these two modifications remain poorly understood. By profiling real 5mC (BS-seq corrected by Tet-assisted BS-seq, TAB-seq) and 5hmC (TAB-seq) levels simultaneously at single-nucleotide resolution, we here demonstrate that there is no global loss of 5mC in kidney tumors compared with matched normal tissues. Conversely, 5hmC was globally lost in virtually all kidney tumor tissues. The 5hmC level in tumor tissues is an independent prognostic marker for kidney cancer, with lower levels of 5hmC associated with shorter overall survival. Furthermore, we demonstrated that loss of 5hmC is linked to hypermethylation in tumors compared with matched normal tissues, particularly in gene body regions. Strikingly, gene body hypermethylation was significantly associated with silencing of the tumor-related genes. Downregulation of IDH1 was identified as a mechanism underlying 5hmC loss in kidney cancer. Restoring 5hmC levels attenuated the invasion capacity of tumor cells and suppressed tumor growth in a xenograft model. Collectively, our results demonstrate that loss of 5hmC is both a prognostic marker and an oncogenic event in kidney cancer by remodeling the DNA methylation pattern.     

Genome-wide antagonism between 5-hydroxymethylcytosine and DNA methylation in the adult mouse brain

Mounting evidence points to critical roles for DNA modifications, including 5-methylcytosine （5mC） and its oxidized forms, in the development, plasticity and disorders of the mammalian nervous system. The novel DNA base 5- hydroxymethylcytosine （5hmC） is known to be capable of initiating passive or active DNA demethylation, but whether and how extensively 5hmC functions in shaping the post-mitotic neuronal DNA methylome is unclear. Here we report the genome-wide distribution of 5hmC in dentate granule neurons from adult mouse hippocampus in vivo. 5hmC in the neuronal genome is highly enriched in gene bodies, especially in exons, and correlates with gene expression. Direct genome-wide comparison of 5hmC distribution between embryonic stem cells and neurons reveals extensive differences, reflecting the functional disparity between these two cell types. Importantly, integrative analysis of 5hmC, overall DNA methylation and gene expression profiles of dentate granule neurons in vivo reveals the genome-wide antagonism between these two states of cytosine modifications, supporting a role for 5hmC in shaping the neuronal DNA methylome by promoting active DNA demethylation.     

Dynamic reprogramming of 5-hydroxymethylcytosine during early porcine embryogenesis

DNA active demethylation is an important epigenetic phenomenon observed in porcine zygotes, yet its molecular origins are unknown. Our results show that 5-methylcytosine (5mC) converts into 5-hydroxymethylcytosine (5hmC) during the first cell cycle in porcine in vivo fertilization (IVV), IVF, and SCNT embryos, but not in parthenogenetically activated embryos. Expression of Ten-Eleven Translocation 1 (TET1) correlates with this conversion. Expression of 5mC gradually decreases until the morula stage; it is only expressed in the inner cell mass, but not trophectoderm regions of IVV and IVF blastocysts. Expression of 5mC in SCNT embryos is ectopically distinct from that observed in IVV and IVF embryos. In addition, 5hmC expression was similar to that of 5mC in IVV cleavage-stage embryos. Expression of 5hmC remained constant in IVF and SCNT embryos, and was evenly distributed among the inner cell mass and trophectoderm regions derived from IVV, IVF, and SCNT blastocysts. Ten-Eleven Translocation 3 was highly expressed in two-cell embryos, whereas TET1 and TET2 were highly expressed in blastocysts. These data suggest that TET1-catalyzed 5hmC may be involved in active DNA demethylation in porcine early embryos. In addition, 5mC, but not 5hmC, participates in the initial cell lineage specification in porcine IVV and IVF blastocysts. Last, SCNT embryos show aberrant 5mC and 5hmC expression during early porcine embryonic development.     

Age-related epigenome-wide DNA methylation and hydroxymethylation in longitudinal mouse blood

DNA methylation at cytosine-phosphate-guanine (CpG) dinucleotides changes as a function of age in humans and animal models, a process that may contribute to chronic disease development. Recent studies have investigated the role of an oxidized form of DNA methylation – 5-hydroxymethylcytosine (5hmC) – in the epigenome, but its contribution to age-related DNA methylation remains unclear. We tested the hypothesis that 5hmC changes with age, but in a direction opposite to 5-methylcytosine (5mC), potentially playing a distinct role in aging. To characterize epigenetic aging, genome-wide 5mC and 5hmC were measured in longitudinal blood samples (2, 4, and 10 months of age) from isogenic mice using two sequencing methods – enhanced reduced representation bisulfite sequencing and hydroxymethylated DNA immunoprecipitation sequencing. Examining the epigenome by age, we identified 28,196 unique differentially methylated CpGs (DMCs) and 8,613 differentially hydroxymethylated regions (DHMRs). Mouse blood showed a general pattern of epigenome-wide hypermethylation and hypo-hydroxymethylation with age. Comparing age-related DMCs and DHMRs, 1,854 annotated genes showed both differential 5mC and 5hmC, including one gene – Nfic – at five CpGs in the same 250 bp chromosomal region. At this region, 5mC and 5hmC levels both decreased with age. Reflecting these age-related epigenetic changes, Nfic RNA expression in blood decreased with age, suggesting that age-related regulation of this gene may be driven by 5hmC, not canonical DNA methylation. Combined, our genome-wide results show age-related differential 5mC and 5hmC, as well as some evidence that changes in 5hmC may drive age-related DNA methylation and gene expression.     

Lead exposure induces changes in 5-hydroxymethylcytosine clusters in CpG islands in human embryonic stem cells and umbilical cord blood

Prenatal exposure to neurotoxicants such as lead (Pb) may cause stable changes in the DNA methylation (5mC) profile of the fetal genome. However, few studies have examined its effect on the DNA de-methylation pathway, specifically the dynamic changes of the 5-hydroxymethylcytosine (5hmC) profile. Therefore, in this study, we investigate the relationship between Pb exposure and 5mC and 5hmC modifications during early development. To study the changes in the 5hmC profile, we use a novel modification of the Infinium™ HumanMethylation450 assay (Illumina, Inc.), which we named HMeDIP-450K assay, in an in vitro human embryonic stem cell model of Pb exposure. We model Pb exposure-associated 5hmC changes as clusters of correlated, adjacent CpG sites, which are co-responding to Pb. We further extend our study to look at Pb-dependent changes in high density 5hmC regions in umbilical cord blood DNA from 48 mother-infant pairs from the Early Life Exposure in Mexico to Environmental Toxicants (ELEMENT) cohort. For our study, we randomly selected umbilical cord blood from 24 male and 24 female children from the 1st and 4th quartiles of Pb levels. Our data show that Pb-associated changes in the 5hmC and 5mC profiles can be divided into sex-dependent and sex-independent categories. Interestingly, differential 5mC sites are better markers of Pb-associated sex-dependent changes compared to differential 5hmC sites. In this study we identified several 5hmC and 5mC genomic loci, which we believe might have some potential as early biomarkers of prenatal Pb exposure.     

Hydroxymethylcytosine and demethylation of the γ-globin gene promoter during erythroid differentiation

The mechanism responsible for developmental stage-specific regulation of γ-globin gene expression involves DNA methylation. Previous results have shown that the γ-globin promoter is nearly fully demethylated during fetal liver erythroid differentiation and partially demethylated during adult bone marrow erythroid differentiation. The hypothesis that 5-hydroxymethylcytosine (5hmC), a known intermediate in DNA demethylation pathways, is involved in demethylation of the γ-globin gene promoter during erythroid differentiation was investigated by analyzing levels of 5-methylcytosine (5mC) and 5hmC at a CCGG site within the 5′ γ-globin gene promoter region in FACS-purified cells from baboon bone marrow and fetal liver enriched for different stages of erythroid differentiation. Our results show that 5mC and 5hmC levels at the γ-globin promoter are dynamically modulated during erythroid differentiation with peak levels of 5hmC preceding and/or coinciding with demethylation. The Tet2 and Tet3 dioxygenases that catalyze formation of 5hmC are expressed during early stages of erythroid differentiation and Tet3 expression increases as differentiation proceeds. In baboon CD34+ bone marrow-derived erythroid progenitor cell cultures, γ-globin expression was positively correlated with 5hmC and negatively correlated with 5mC at the γ-globin promoter. Supplementation of culture media with Vitamin C, a cofactor of the Tet dioxygenases, reduced γ-globin promoter DNA methylation and increased γ-globin expression when added alone and in an additive manner in combination with either DNA methyltransferase or LSD1 inhibitors. These results strongly support the hypothesis that the Tet-mediated 5hmC pathway is involved in developmental stage-specific regulation of γ-globin expression by mediating demethylation of the γ-globin promoter.     

Genome-wide regulation of 5hmC, 5mC, and gene expression by Tet1 hydroxylase in mouse embryonic stem cells

DNA methylation at the 5 position of cytosine (5mC) in the mammalian genome is a key epigenetic event critical for various cellular processes. The ten-eleven translocation (Tet) family of 5mC-hydroxylases, which convert 5mC to 5-hydroxymethylcytosine (5hmC), offers a way for dynamic regulation of DNA methylation. Here we report that Tet1 binds to unmodified C or 5mC- or 5hmC-modified CpG-rich DNA through its CXXC domain. Genome-wide mapping of Tet1 and 5hmC reveals mechanisms by which Tet1 controls 5hmC and 5mC levels in mouse embryonic stem cells (mESCs). We also uncover a comprehensive gene network influenced by Tet1. Collectively, our data suggest that Tet1 controls DNA methylation both by binding to CpG-rich regions to prevent unwanted DNA methyltransferase activity, and by converting 5mC to 5hmC through hydroxylase activity. This Tet1-mediated antagonism of CpG methylation imparts differential maintenance of DNA methylation status at Tet1 targets, ultimately contributing to mESC differentiation and the onset of embryonic development.     

5‐Hydroxymethylcytosine: a new kid on the epigenetic block?

The discovery of the Ten‐Eleven‐Translocation (TET) oxygenases that catalyze the hydroxylation of 5‐methylcytosine (5mC) to 5‐hydroxymethylcytosine (5hmC) has triggered an avalanche of studies aiming to resolve the role of 5hmC in gene regulation if any. Hitherto, TET1 is reported to bind to CpG‐island (CGI) and bivalent promoters in mouse embryonic stem cells, whereas binding at DNAseI hypersensitive sites (HS) had escaped previous analysis. Significant enrichment/accumulation of 5hmC but not 5mC can indeed be detected at bivalent promoters and at DNaseI‐HS. Surprisingly, however, 5hmC is not detected or present at very low levels at CGI promoters notwithstanding the presence of TET1. Our meta‐analysis of DNA methylation profiling points to potential issues with regard to the various methodologies that are part of the toolbox used to detect 5mC and 5hmC. Discrepancies between published studies and technical limitations prevent an unambiguous assignment of 5hmC as a ‘true’ epigenetic mark, that is, read and interpreted by other factors and/or as a transiently accumulating intermediary product of the conversion of 5mC to unmodified cytosines.