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.     

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.     

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.     

Structure of Naegleria Tet-like dioxygenase (NgTet1) in complexes with a reaction intermediate 5-hydroxymethylcytosine DNA

The family of ten-eleven translocation (Tet) dioxygenases is widely distributed across the eukaryotic tree of life, from mammals to the amoeboflagellate Naegleria gruberi. Like mammalian Tet proteins, the Naegleria Tet-like protein, NgTet1, acts on 5-methylcytosine (5mC) and generates 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC) in three consecutive, Fe(II)- and α-ketoglutarate-dependent oxidation reactions. The two intermediates, 5hmC and 5fC, could be considered either as the reaction product of the previous enzymatic cycle or the substrate for the next cycle. Here we present a new crystal structure of NgTet1 in complex with DNA containing a 5hmC. Along with the previously solved NgTet1–5mC structure, the two complexes offer a detailed picture of the active site at individual stages of the reaction cycle. In the crystal, the hydroxymethyl (OH-CH₂-) moiety of 5hmC points to the metal center, representing the reaction product of 5mC hydroxylation. The hydroxyl oxygen atom could be rotated away from the metal center, to a hydrophobic pocket formed by Ala212, Val293 and Phe295. Such rotation turns the hydroxyl oxygen atom away from the product conformation, and exposes the target CH₂ towards the metal-ligand water molecule, where a dioxygen O₂ molecule would occupy to initiate the next round of reaction by abstracting a hydrogen atom from the substrate. The Ala212-to-Val (A212V) mutant profoundly limits the product to 5hmC, probably because the reduced hydrophobic pocket size restricts the binding of 5hmC as a substrate.     

Transient reduction of 5-methylcytosine and 5-hydroxymethylcytosine is associated with active DNA demethylation during regeneration of zebrafish fin

Although dedifferentiation, transformation of differentiated cells into progenitor cells, is a critical step in the regeneration of amphibians and fish, the molecular mechanisms underlying this process, including epigenetic changes, remain unclear. Dot blot assays and immunohistochemical analyses revealed that, during regeneration of zebrafish fin, the levels of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) are transiently reduced in blastema cells and cells adjacent to the amputation plane at 30 h post-amputation (hpa), and the level of 5mC, but not 5hmC, is almost restored by 72 hpa. We observed that the dedifferentiated cells showed reduced levels of 5mC and 5hmC independent of cell proliferation by 24 hpa. Furthermore, expressions of the proposed demethylation- and DNA repair-related genes were detected during fin regeneration. Taken together, our findings illustrate that the transient reduction of 5mC and 5hmC in dedifferentiated cells is associated with active demethylation during regeneration of zebrafish fin.     

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.     

Placental 5-methylcytosine and 5-hydroxymethylcytosine patterns associate with size at birth

Altered placental function as a consequence of aberrant imprinted gene expression may be one mechanism mediating the association between low birth weight and increased cardiometabolic disease risk. Imprinted gene expression is regulated by epigenetic mechanisms, particularly DNA methylation (5mC) at differentially methylated regions (DMRs). While 5-hydroxymethylcytosine (5hmC) is also present at DMRs, many techniques do not distinguish between 5mC and 5hmC. Using human placental samples, we show that the expression of the imprinted gene CDKN1C associates with birth weight. Using specific techniques to map 5mC and 5hmC at DMRs controlling the expression of CDKN1C and the imprinted gene IGF2, we show that 5mC enrichment at KvDMR and DMR0, and 5hmC enrichment within the H19 gene body, associate positively with birth weight. Importantly, the presence of 5hmC at imprinted DMRs may complicate the interpretation of DNA methylation studies in placenta; future studies should consider using techniques that distinguish between, and permit quantification of, both modifications.     

Chemical Display of Pyrimidine Bases Flipped Out by Modification-Dependent Restriction Endonucleases of MspJI and PvuRts1I Families

The epigenetic DNA modifications 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) in eukaryotes are recognized either in the context of double-stranded DNA (e.g., by the methyl-CpG binding domain of MeCP2), or in the flipped-out state (e.g., by the SRA domain of UHRF1). The SRA-like domains and the base-flipping mechanism for 5(h)mC recognition are also shared by the recently discovered prokaryotic modification-dependent endonucleases of the MspJI and PvuRts1I families. Since the mechanism of modified cytosine recognition by many potential eukaryotic and prokaryotic 5(h)mC “readers” is still unknown, a fast solution based method for the detection of extrahelical 5(h)mC would be very useful. In the present study we tested base-flipping by MspJI- and PvuRts1I-like restriction enzymes using several solution-based methods, including fluorescence measurements of the cytosine analog pyrrolocytosine and chemical modification of extrahelical pyrimidines with chloroacetaldehyde and KMnO₄. We find that only KMnO₄ proved an efficient probe for the positive display of flipped out pyrimidines, albeit the method required either non-physiological pH (4.3) or a substitution of the target cytosine with thymine. Our results imply that DNA recognition mechanism of 5(h)mC binding proteins should be tested using a combination of all available methods, as the lack of a positive signal in some assays does not exclude the base flipping mechanism.     

The role of 5-hydroxymethylcytosine in human cancer

The patterns of DNA methylation in human cancer cells are highly abnormal and often involve the acquisition of DNA hypermethylation at hundreds or thousands of CpG islands that are usually unmethylated in normal tissues. The recent discovery of 5-hydroxymethylcytosine (5hmC) as an enzymatic oxidation product of 5-methylcytosine (5mC) has led to models and experimental data in which the hypermethylation and 5mC oxidation pathways seem to be connected. Key discoveries in this setting include the findings that several genes coding for proteins involved in the 5mC oxidation reaction are mutated in human tumors, and that a broad loss of 5hmC occurs across many types of cancer. In this review, we will summarize current knowledge and discuss models of the potential roles of 5hmC in human cancer biology.     

Hydroxyl-radical-induced oxidation of 5-methylcytosine in isolated and cellular DNA

The methylation and oxidative demethylation of cytosine in CpG dinucleotides plays a critical role in the regulation of genes during cell differentiation, embryogenesis and carcinogenesis. Despite its low abundance, 5-methylcytosine (5mC) is a hotspot for mutations in mammalian cells. Here, we measured five oxidation products of 5mC together with the analogous products of cytosine and thymine in DNA exposed to ionizing radiation in oxygenated aqueous solution. The products can be divided into those that arise from hydroxyl radical (•OH) addition at the 5,6-double bond of 5mC (glycol, hydantoin and imidazolidine products) and those that arise from H-atom abstraction from the methyl group of 5mC including 5-hydroxymethylcytosine (5hmC) and 5-formylcytosine (5fC). Based on the analysis of these products, we show that the total damage at 5mC is about 2-fold greater than that at C in identical sequences. The formation of hydantoin products of 5mC is favored, compared to analogous reactions of thymine and cytosine, which favor the formation of glycol products. The distribution of oxidation products is sequence dependent in specific ODN duplexes. In the case of 5mC, the formation of 5hmC and 5fC represents about half of the total of •OH-induced oxidation products of 5mC. Several products of thymine, cytosine, 5mC, as well as 8-oxo-7,8-dihydroguanine (8oxoG), were also estimated in irradiated cells.     

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.     

A sensitive mass spectrometry method for simultaneous quantification of DNA methylation and hydroxymethylation levels in biological samples

The recent discovery of 5-hydroxymethyl-cytosine (5hmC) in embryonic stem cells and postmitotic neurons has triggered the need for quantitative measurements of both 5-methyl-cytosine (5mC) and 5hmC in the same sample. We have developed a method using liquid chromatography electrospray ionization tandem mass spectrometry with multiple reaction monitoring (LC–ESI–MS/MS–MRM) to simultaneously measure levels of 5mC and 5hmC in digested genomic DNA. This method is fast, robust, and accurate, and it is more sensitive than the current 5hmC quantitation methods such as end labeling with thin layer chromatography and radiolabeling by glycosylation. Only 50 ng of digested genomic DNA is required to measure the presence of 0.1% 5hmC in DNA from mouse embryonic stem cells. Using this procedure, we show that human induced pluripotent stem cells exhibit a dramatic increase in 5mC and 5hmC levels compared with parental fibroblast cells, suggesting a dynamic regulation of DNA methylation and hydroxymethylation during cellular reprogramming.     

Mouse olfactory bulb methylome and hydroxymethylome maps reveal noncanonical active turnover of DNA methylation

Hydroxylation of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) by TET enzymes presents a particular regulatory mechanism in the mammalian brain. However, although methylation and hydroxymethylation of cytosines in non-CpG contexts have been reported, these mechanisms remain poorly understood. Here, we applied TAB-seq and oxBS-seq selectively to detect 5hmC and 5mC at base resolution in olfactory bulb derived from female mice. We found that active turnover of 5mC to 5hmC occurred in both CpG and non-CpG contexts. Strikingly, we identified a different sequence preference for 5mC and 5hmC in a CH context, in which H = A, C, or T, TNCA/TC for 5mC and NNCA/T/CN for 5hmC. More importantly, we found that genes showing 5mC to 5hmC turnover showed only limited overlap in CpG and CH contexts, and that olfactory receptor genes were marked with higher turnover of 5mC to 5hmC in non-CpG context. Collectively, we identified an unexpected sequence preference for non-CpG hydroxymethylation and distinct target genes regulated by the turnover of 5mC to 5hmC in CpG and CH contexts.