Selective Capture of 5-hydroxymethylcytosine from Genomic DNA

5-methylcytosine (5-mC) constitutes ~2-8% of the total cytosines in human genomic DNA and impacts a broad range of biological functions, including gene expression, maintenance of genome integrity, parental imprinting, X-chromosome inactivation, regulation of development, aging, and cancer¹. Recently, the presence of an oxidized 5-mC, 5-hydroxymethylcytosine (5-hmC), was discovered in mammalian cells, in particular in embryonic stem (ES) cells and neuronal cells^2-4. 5-hmC is generated by oxidation of 5-mC catalyzed by TET family iron (II)/α-ketoglutarate-dependent dioxygenases^{2, 3}. 5-hmC is proposed to be involved in the maintenance of embryonic stem (mES) cell, normal hematopoiesis and malignancies, and zygote development^{2, 5-10}. To better understand the function of 5-hmC, a reliable and straightforward sequencing system is essential. Traditional bisulfite sequencing cannot distinguish 5-hmC from 5-mC¹¹. To unravel the biology of 5-hmC, we have developed a highly efficient and selective chemical approach to label and capture 5-hmC, taking advantage of a bacteriophage enzyme that adds a glucose moiety to 5-hmC specifically¹².Here we describe a straightforward two-step procedure for selective chemical labeling of 5-hmC. In the first labeling step, 5-hmC in genomic DNA is labeled with a 6-azide-glucose catalyzed by β-GT, a glucosyltransferase from T4 bacteriophage, in a way that transfers the 6-azide-glucose to 5-hmC from the modified cofactor, UDP-6-N3-Glc (6-N3UDPG). In the second step, biotinylation, a disulfide biotin linker is attached to the azide group by click chemistry. Both steps are highly specific and efficient, leading to complete labeling regardless of the abundance of 5-hmC in genomic regions and giving extremely low background. Following biotinylation of 5-hmC, the 5-hmC-containing DNA fragments are then selectively captured using streptavidin beads in a density-independent manner. The resulting 5-hmC-enriched DNA fragments could be used for downstream analyses, including next-generation sequencing.Our selective labeling and capture protocol confers high sensitivity, applicable to any source of genomic DNA with variable/diverse 5-hmC abundances. Although the main purpose of this protocol is its downstream application (i.e., next-generation sequencing to map out the 5-hmC distribution in genome), it is compatible with single-molecule, real-time SMRT (DNA) sequencing, which is capable of delivering single-base resolution sequencing of 5-hmC.     

Biochemical characterization of recombinant β-glucosyltransferase and analysis of global 5-hydroxymethylcytosine in unique genomes

5-Hydroxymethylcytosine (5-hmC) is an enzymatic oxidative product of 5-methylcytosine (5-mC). The Ten Eleven Translocation (TET) family of enzymes catalyze the conversion of 5-mC to 5-hmC. Phage-encoded glucosyltransferases are known to glucosylate 5-hmC, which can be utilized to detect and analyze the 5-hmC as an epigenetic mark in the mammalian epigenome. Here we have performed a detailed biochemical characterization and steady-state kinetic parameter analysis of T4 phage β-glucosyltransferase (β-GT). Recombinant β-GT glucosylates 5-hmC DNA in a nonprocessive manner, and binding to either 5-hmC DNA or uridine diphosphoglucose (UDP-glucose) substrates is random, with both binary complexes being catalytically competent. Product inhibition studies with β-GT demonstrated that UDP is a competitive inhibitor with respect to UDP-glucose and a mixed inhibitor with respect to 5-hmC DNA. Similarly, the glucosylated-5-hmC (5-ghmC) DNA is a competitive inhibitor with respect to 5-hmC DNA and mixed inhibitor with respect to UDP-glucose. 5-hmC DNA binds ~10 fold stronger to the β-GT enzyme when compared to its glucosylated product. The numbers of 5-hmC on target sequences influenced the turnover numbers for recombinant β-GT. Furthermore, we have utilized recombinant β-GT to estimate global 5-hmC content in a variety of genomic DNAs. Most of the genomic DNAs derived from vertebrate tissue and cell lines contained 5-hmC. DNA from mouse, human, and bovine brains displayed 0.5-0.9% of the total nucleotides as 5-hmC, which was higher compared to the levels found in other tissues. A comparison between cancer and healthy tissue genomes suggested a lower percentage of 5-hmC in cancer, which may reflect the global hypomethylation of 5-mC observed during oncogenesis.     

The Behaviour of 5-Hydroxymethylcytosine in Bisulfite Sequencing

Background

We recently showed that enzymes of the TET family convert 5-mC to 5-hydroxymethylcytosine (5-hmC) in DNA. 5-hmC is present at high levels in embryonic stem cells and Purkinje neurons. The methylation status of cytosines is typically assessed by reaction with sodium bisulfite followed by PCR amplification. Reaction with sodium bisulfite promotes cytosine deamination, whereas 5-methylcytosine (5-mC) reacts poorly with bisulfite and is resistant to deamination. Since 5-hmC reacts with bisulfite to yield cytosine 5-methylenesulfonate (CMS), we asked how DNA containing 5-hmC behaves in bisulfite sequencing.

Methodology/Principal Findings

We used synthetic oligonucleotides with different distributions of cytosine as templates for generation of DNAs containing C, 5-mC and 5-hmC. The resulting DNAs were subjected in parallel to bisulfite treatment, followed by exposure to conditions promoting cytosine deamination. The extent of conversion of 5-hmC to CMS was estimated to be 99.7%. Sequencing of PCR products showed that neither 5-mC nor 5-hmC undergo C-to-T transitions after bisulfite treatment, confirming that these two modified cytosine species are indistinguishable by the bisulfite technique. DNA in which CMS constituted a large fraction of all bases (28/201) was much less efficiently amplified than DNA in which those bases were 5-mC or uracil (the latter produced by cytosine deamination). Using a series of primer extension experiments, we traced the inefficient amplification of CMS-containing DNA to stalling of Taq polymerase at sites of CMS modification, especially when two CMS bases were either adjacent to one another or separated by 1–2 nucleotides.

Conclusions

We have confirmed that the widely used bisulfite sequencing technique does not distinguish between 5-mC and 5-hmC. Moreover, we show that CMS, the product of bisulfite conversion of 5-hmC, tends to stall DNA polymerases during PCR, suggesting that densely hydroxymethylated regions of DNA may be underrepresented in quantitative methylation analyses.     

Association of 5-Methylcytosine and 5-Hydroxymethylcytosine with Mitochondrial DNA Content and Clinical and Biochemical Parameters in Hepatocellular Carcinoma

Increasing epidemiological evidence has indicated that inherited variations of mitochondrial DNA (mtDNA) copy number affect the genetic susceptibility of many malignancies in a tumour-specific manner and that DNA methylation also plays an important role in controlling gene expression during the differentiation and development of hepatocellular carcinoma (HCC). Our previous study demonstrated that HCC tissues showed a lower 5-hydroxymethylcytosine (5-hmC) content when compared to tumour-adjacent tissues, but the relationship among 5-hmC, 5-methylcytosine (5-mC) and mtDNA content in HCC patients is still unknown. This study aimed to clarify the correlation among mtDNA content, 5-mC and 5-hmC by quantitative real-time PCR and liquid chromatography tandem mass spectrometry analysis. We demonstrated that 5-hmC correlated with tumour size [odds ratio (OR) 0.847, 95% confidence interval (CI) 0.746–0.962, P = 0.011], and HCC patients with a tumour size ≥5.0 cm showed a lower 5-hmC content and higher levels of fasting plasma aspartate aminotransferase, the ratio of alanine amiotransferase to aspartate aminotransferase, γ-glutamyltransferase, alpha-fetoprotein than those with a tumour size <5 cm (all P<0.05). We further revealed that the mtDNA content of HCC tumour tissues was 225.97(105.42, 430.54) [median (25th Percentile, 75th Percentile)] and was negatively correlated with 5-mC content (P = 0.035), but not 5-hmC content, in genomic DNA from HCC tumour tissues.     

Crystal structures of B-DNA dodecamer containing the epigenetic modifications 5-hydroxymethylcytosine or 5-methylcytosine

5-Hydroxymethylcytosine (5-hmC) was recently identified as a relatively frequent base in eukaryotic genomes. Its physiological function is still unclear, but it is supposed to serve as an intermediate in DNA de novo demethylation. Using X-ray diffraction, we solved five structures of four variants of the d(CGCGAATTCGCG) dodecamer, containing either 5-hmC or 5-methylcytosine (5-mC) at position 3 or at position 9. The observed resolutions were between 1.42 and 1.99 Å. Cytosine modification in all cases influences neither the whole B-DNA double helix structure nor the modified base pair geometry. The additional hydroxyl group of 5-hmC with rotational freedom along the C5-C5A bond is preferentially oriented in the 3′ direction. A comparison of thermodynamic properties of the dodecamers shows no effect of 5-mC modification and a sequence-dependent only slight destabilizing effect of 5-hmC modification. Also taking into account the results of a previous functional study [Münzel et al. (2011) (Improved synthesis and mutagenicity of oligonucleotides containing 5-hydroxymethylcytosine, 5-formylcytosine and 5-carboxylcytosine. Chem. Eur. J., 17, 13782−13788)], we conclude that the 5 position of cytosine is an ideal place to encode epigenetic information. Like this, neither the helical structure nor the thermodynamics are changed, and polymerases cannot distinguish 5-hmC and 5-mC from unmodified cytosine, all these effects are making the former ones non-mutagenic.     

Lineage-specific distribution of high levels of genomic 5-hydroxymethylcytosine in mammalian development

Methylation of cytosine is a DNA modification associated with gene repression. Recently, a novel cytosine modification, 5-hydroxymethylcytosine (5-hmC) has been discovered. Here we examine 5-hmC distribution during mammalian development and in cellular systems, and show that the developmental dynamics of 5-hmC are different from those of 5-methylcytosine (5-mC); in particular 5-hmC is enriched in embryonic contexts compared to adult tissues. A detectable 5-hmC signal appears in pre-implantation development starting at the zygote stage, where the paternal genome is subjected to a genome-wide hydroxylation of 5-mC, which precisely coincides with the loss of the 5-mC signal in the paternal pronucleus. Levels of 5-hmC are high in cells of the inner cell mass in blastocysts, and the modification colocalises with nestin-expressing cell populations in mouse post-implantation embryos. Compared to other adult mammalian organs, 5-hmC is strongly enriched in bone marrow and brain, wherein high 5-hmC content is a feature of both neuronal progenitors and post-mitotic neurons. We show that high levels of 5-hmC are not only present in mouse and human embryonic stem cells (ESCs) and lost during differentiation, as has been reported previously, but also reappear during the generation of induced pluripotent stem cells; thus 5-hmC enrichment correlates with a pluripotent cell state. Our findings suggest that apart from the cells of neuronal lineages, high levels of genomic 5-hmC are an epigenetic feature of embryonic cell populations and cellular pluri- and multi-lineage potency. To our knowledge, 5-hmC represents the first epigenetic modification of DNA discovered whose enrichment is so cell-type specific.     

DNA methylation in mouse cells in culture as measured by restriction enzymes

Methylation of DNA in normal mouse cultured 3T3 cells and in their virally or chemically transformed derivatives was studied. DNA methylation was studied by restriction with HpaII, MspI, or HpaII plus MspI. DNA from the chemically transformed cells was cleaved about twice as often with HpaII than was the DNA of normal and virally transformed cells. Digests with MspI and HpaII plus MspI were identical in all cell lines studied. Densitometry of the restriction patterns allowed an estimate of total DNA methylation from the weight average lengths. The chemically transformed cell line showed 25% reduction in methylation compared to the other cell lines. Southern blot hybridization using satellite DNA showed that these sequences followed a pattern of modification similar to that of total DNA.     

Restriction and modification in Bacillus subtilis: sequence specificities of restriction/modification systems BsuM, BsuE, and BsuF.

The sequence specificities of three Bacillus subtilis restriction/modification systems were established: (i) BsuM (CTCGAG), an isoschizomer to XhoI; (ii) BsuE (CGCG), an isoschizomer to FnuDII; and (iii) BsuF (CCGG), an isoschizomer to MspI, HpaII. The BsuM modification enzyme methylates the 3' cytosine of the recognition sequence. The BsuF modification enzyme methylates the 5' cytosine of the sequence, rendering such sites resistant to MspI degradation and leaving the majority of sites sensitive to HpaII degradation.     

Methylation-sensitive restriction endonuclease digestion patterns revealed in Vicia faba L. chromosomes by in situ nick-translation

Restriction endonucleases sensitive to cytosine methylation (HpaII, MspI and HhaI) and 5-azacitidine were used to study the localization of target sequences in Vicia faba metaphase chromosomes by in situ digestion and radioactive or non-radioactive nick-translation. In control experiments, neither isolated DNA nor chromosomes in situ were digested by HpaII and MspI. Pretreatment with demethylating agent, 5-azacitidine resulted both in increased effectiveness of in situ digestion and nick-translation. In 5-azacitidine-treated material, negative bands in M chromosomes appeared. HhaI cleaved isolated DNA, digested it in situ and gave positive signals as a result of nick-translation procedure in metaphase chromosomes. In S chromosomes containing heterochromatin without target sequences for HpaII and MspI, negative bands were shown after nick-translation. Such heterochromatin contains FokI sequences and in situ nick-translation driven by that restriction enzyme resulted in positive bands.     

Selective chemical labeling reveals the genome-wide distribution of 5-hydroxymethylcytosine

In contrast to 5-methylcytosine (5-mC), which has been studied extensively, little is known about 5-hydroxymethylcytosine (5-hmC), a recently identified epigenetic modification present in substantial amounts in certain mammalian cell types. Here we present a method for determining the genome-wide distribution of 5-hmC. We use the T4 bacteriophage β-glucosyltransferase to transfer an engineered glucose moiety containing an azide group onto the hydroxyl group of 5-hmC. The azide group can be chemically modified with biotin for detection, affinity enrichment and sequencing of 5-hmC-containing DNA fragments in mammalian genomes. Using this method, we demonstrate that 5-hmC is present in human cell lines beyond those previously recognized. We also find a gene expression level-dependent enrichment of intragenic 5-hmC in mouse cerebellum and an age-dependent acquisition of this modification in specific gene bodies linked to neurodegenerative disorders.     

Changes in methylation of tissue cultured soybean cells detected by digestion with the restriction enzymes HpaII and MspI

Each of the tandemly arranged 5S RNA genes of soybean contain two CCGG sites which, if unmethylated, can be digested by both MspI and HpaII. Methylation of the internal cytosine (CmeCGG) prevents digestion by HpaII but allows digestions by MspI.Suspension cultures were prepared from soybean plants and the DNA from these cultures was examined for the susceptibility of 5S RNA genes to digestion by MspI and HpaII. 5S genes from DNA extracted from intact plants can be partially digested with MspI but not at all by HpaII. In contrast, shortly after cells were cultured the 5S RNA could be hydrolyzed by both HpaII and MspI. After prolonged cell culture, the 5S genes from some cell lines were found to have become partially or even completely resistant to HpaII digestion. The results suggest that lack of methylation can occur when cells are cultured and that such methylation may play a role in the heritable changes observed in cell culture.Research supported by Grant 01498 from the National Institutes of Environmental Health Sciences     

TET2-mediated 5-hydroxymethylcytosine induces genetic instability and mutagenesis

The family of Ten-Eleven Translocation (TET) proteins is implicated in the process of active DNA demethylation and thus in epigenetic regulation. TET 1, 2 and 3 proteins are oxygenases that can hydroxylate 5-methylcytosine (5-mC) into 5-hydroxymethylcytosine (5-hmC) and further oxidize 5-hmC into 5-formylcytosine (5-fC) and 5-carboxylcytosine (5-caC). The base excision repair (BER) pathway removes the resulting 5-fC and 5-caC bases paired with a guanine and replaces them with regular cytosine. The question arises whether active modification of 5-mC residues and their subsequent elimination could affect the genomic DNA stability. Here, we generated two inducible cell lines (Ba/F3-EPOR, and UT7) overexpressing wild-type or catalytically inactive human TET2 proteins. Wild-type TET2 induction resulted in an increased level of 5-hmC and a cell cycle defect in S phase associated with higher level of phosphorylated P53, chromosomal and centrosomal abnormalities. Furthermore, in a thymine-DNA glycosylase (Tdg) deficient context, the TET2-mediated increase of 5-hmC induces mutagenesis characterized by GC > AT transitions in CpG context suggesting a mutagenic potential of 5-hmC metabolites. Altogether, these data suggest that TET2 activity and the levels of 5-hmC and its derivatives should be tightly controlled to avoid genetic and chromosomal instabilities. Moreover, TET2-mediated active demethylation might be a very dangerous process if used to entirely demethylate the genome and might rather be used only at specific loci.     

DNA binding of the p21 repressor ZBTB2 is inhibited by cytosine hydroxymethylation

Recent studies have demonstrated that the modified base 5-hydroxymethylcytosine (5-hmC) is detectable at various rates in DNA extracted from human tissues. This oxidative product of 5-methylcytosine (5-mC) constitutes a new and important actor of epigenetic mechanisms. We designed a DNA pull down assay to trap and identify nuclear proteins bound to 5-hmC and/or 5-mC. We applied this strategy to three cancerous cell lines (HeLa, SH-SY5Y and UT7-MPL) in which we also measured 5-mC and 5-hmC levels by HPLC-MS/MS. We found that the putative oncoprotein Zinc finger and BTB domain-containing protein 2 (ZBTB2) is associated with methylated DNA sequences and that this interaction is inhibited by the presence of 5-hmC replacing 5-mC. As published data mention ZBTB2 recognition of p21 regulating sequences, we verified that this sequence specific binding was also alleviated by 5-hmC. ZBTB2 being considered as a multifunctional cell proliferation activator, notably through p21 repression, this work points out new epigenetic processes potentially involved in carcinogenesis.     

Modulation of Ten-Eleven Translocation 1 (TET1), Isocitrate Dehydrogenase (IDH) Expression, α-Ketoglutarate (α-KG), and DNA Hydroxymethylation Levels by Interleukin-1β in Primary Human Chondrocytes

5-Hydroxymethylcytosine (5-hmC) generated by ten-eleven translocation 1–3 (TET1–3) enzymes is an epigenetic mark present in many tissues with different degrees of abundance. IL-1β and TNF-α are the two major cytokines present in arthritic joints that modulate the expression of many genes associated with cartilage degradation in osteoarthritis. In the present study, we investigated the global 5-hmC content, the effects of IL-1β and TNF-α on 5-hmC content, and the expression and activity of TETs and isocitrate dehydrogenases in primary human chondrocytes. The global 5-hmC content was found to be ∼0.1% of the total genome. There was a significant decrease in the levels of 5-hmC and the TET enzyme activity upon treatment of chondrocytes with IL-1β alone or in combination with TNF-α. We observed a dramatic (10–20-fold) decrease in the levels of TET1 mRNA expression and a small increase (2–3-fold) in TET3 expression in chondrocytes stimulated with IL-1β and TNF-α. IL-1β and TNF-α significantly suppressed the activity and expression of IDHs, which correlated with the reduced α-ketoglutarate levels. Whole genome profiling showed an erasure effect of IL-1β and TNF-α, resulting in a significant decrease in hydroxymethylation in a myriad of genes including many genes that are important in chondrocyte physiology. Our data demonstrate that DNA hydroxymethylation is modulated by pro-inflammatory cytokines via suppression of the cytosine hydroxymethylation machinery. These data point to new mechanisms of epigenetic control of gene expression by pro-inflammatory cytokines in human chondrocytes.     

Methylation of nuclear DNA in Physarum polycephalum.

The restriction endonucleases HpaII and HhaI, whose action is inhibited by the presence of methylated base analogues at the recognition sequences in the DNA substrate, were used to investigate the distribution of 5-methylcytosine in nuclear DNA from Physarum polycephalum. Physarum DNA is digested into two fractions by these enzymes: a low-molecular-weight (M--) compartment comprising 80% of the DNA, and a high-molecular-weight (M+) compartment containing 20% of the DNA. The DNA fraction showing resistance to digestion by restriction endonuclease HpaII is cleaved by its isoschizomer MspI, indicating that methylated endonuclease-HpaII-specific sites are present in M + DNA. Additional properties of sequences in the M+ compartment were investigated.     

Cloning of the MspI modification enzyme. The site of modification and its effects on cleavage by MspI and HpaII

The gene for the MspI modification enzyme from Moraxella was cloned in Escherichia coli using the plasmid vector pBR322. Selection of transformants carrying the gene was based on the resistance of the modified plasmid encoding the enzyme to cleavage by MspI. Both chromosomal and plasmid DNA were modified in the selected clones. None of the clones obtained produced the cognate restriction enzyme which suggests that in this system the genes for the restriction enzyme and methylase are not closely linked. Crude cell extracts prepared from the recombinant strains, but not the host (E. coli HB101), contain an S-adenosylmethionine-dependent methyltransferase specific for the MspI recognition site, CCGG. Production of the enzyme is 3-4-fold greater in the transformants than in the original Moraxella strain. 5-Methylcytosine was identified as the product of the reaction chromatographically. The outer cytosine of the recognition sequence, *CCGG, was shown to be the site of methylation by DNA-sequencing methods. This modification blocks cleavage by both MspI and its isoschizomer HpaII. HpaII, but not MspI, is able to cleave the unmethylated strand of a hemimethylated substrate. The relevance of these results to the use of MspI and HpaII to analyze patterns of methylation in genomic DNA is discussed.     

A simple method for high-throughput quantification of genome-wide DNA methylation by fluorescence polarization

To rapidly determine DNA methylation levels from a large number of biological or clinical samples, we have developed an accurate and sensitive method for high-throughput quantification of global methylation of 5′-Cm⁵CGG-3′ sites in the genome, visualized by fluorescence polarization (FP) based measurement of DNA methylation (FPDM). In FPDM, the methyl-sensitive HpaII and methyl-insensitive MspI restriction enzymes were employed to achieve DNA cleavage, followed by incorporation of fluorescent dCMP into the enzyme-cleavage products through polymerase chain extension, yielding an FP-ratio between the HpaII- and MspI-restricted preparations as a measure of methylation. FPDM provided stable estimates of methylation level of submicrograms of lambda or human DNA, and of a 255-bp DNA segment containing a single HpaII/MspI restriction site in accord with, and more accurate than, determination by gel electrophoresis. FPDM was also applied to measure dose-dependent DNA hypomethylation in human embryonic kidney 293T cells treated with the DNA-methyltransferase inhibitor 5-aza-dC.     

Discrimination between 5-hydroxymethylcytosine and 5-methylcytosine in DNA by selective chemical labeling

Detection of DNA damage has been greatly improved following the development of equipment and techniques, however, discrimination between 5-hydroxymethylcytosine (5-hmC) and 5-methylcytosine (5-mC) is still a thorny problem. In the present study, an approach to oxidize and selective label (Ox-Labeling) 5-hmC in native DNA has been reported, which conveniently distinguishes 5-hmC from 5-mC using simple and effective processes.     

MspI, an isoschizomer of hpaII which cleaves both unmethylated and methylated hpaII sites.

The cleavage of DNA by restriction endonucleases HpaII and HapII is prevented by the presence of a 5-methyl group at the internal C residue of its recognition sequence CCGG. MspI, an isoschizomer of HpaII available from New England Biolabs, cleaves DNA irrespective of the presence of a methyl group at this position. This enzyme cleaves DNA from Haemophilus parainfluenzae and Haemophilus aphrophilus readily while HpaII and HapII cannot degrade these DNAs. Practically all HpaII sites in mammalian sperm DNA are also protected by methylation at the internal C position since HpaII and HapII barely cleave this DNA (average molecular weight 40 kb). MspI, however, cleaves the DNA to an average size of about 5 kb.