Search

DNA 5-hydroxymethylcytosine (5hmC) modification is known to be associated with gene transcription and frequently used as a mark to investigate dynamic DNA methylation conversion during mammalian development and in human diseases. However, the lack of genome-wide 5hmC profiles in different human tissue types impedes drawing generalized conclusions about how 5hmC is implicated in transcription activity and tissue specificity. To meet this need, we describe the development of a 5hmC tissue map by characterizing the genomic distributions of 5hmC in 19 human tissues derived from ten organ systems. Subsequent sequencing results enabled the identification of genome-wide 5hmC distributions that uniquely separates samples by tissue type. Further comparison of the 5hmC profiles with transcriptomes and histone modifications revealed that 5hmC is preferentially enriched on tissue-specific gene bodies and enhancers. Taken together, the results provide an extensive 5hmC map across diverse human tissue types that suggests a potential role of 5hmC in tissue-specific development; as well as a resource to facilitate future studies of DNA demethylation in pathogenesis and the development of 5hmC as biomarkers.

Applications associated with nucleobase protonation events are grounded in their fundamental impact on DNA thermodynamics, structure, and hybridization dynamics. Of the canonical nucleobases, N3 protonation of cytosine (C) is the most widely utilized in both biology and nanotechnology. Naturally occurring C derivatives that shift the N3 p K a introduce an additional level of tunability. The epigenetic nucleobase 5-carboxylcytosine (caC) presents a particularly interesting example since this derivative forms Watson-Crick base pairs of similar stability and displays pH-dependent behavior over the same range as the canonical nucleobase. However, the titratable group in caC corresponds to the exocyclic carboxyl group rather than N3, and the implications of these divergent protonation events toward DNA hybridization thermodynamics, kinetics, and base pairing dynamics remain poorly understood. Here, we study the pH dependence of these physical properties using model oligonucleotides containing C and caC with FTIR and temperature-jump IR spectroscopy. We demonstrate that N3 protonation of C completely disrupts duplex stability, leading to large shifts in the duplex/single-strand equilibrium, a reduction in the cooperativity of melting, and an acceleration in the rate of duplex dissociation. In contrast, while increasing 5-carboxyl protonation in caC-containing duplexes induces an increase in base pair fluctuations, the DNA duplex can tolerate substantial protonation without significant perturbation to the duplex/single-strand equilibrium. However, 5-carboxyl protonation has a large impact on hybridization kinetics by reducing the transition state free energy. Our thermodynamic and kinetic analysis provides new insight on the impact of two divergent protonation mechanisms in naturally occurring nucleobases on the biophysical properties of DNA.

Background: Directed DNA methylation on N6-adenine (6mA), N4-cytosine (4mC), and C5-cytosine (5mC) can potentially increase DNA coding capacity and regulate a variety of biological functions. These modifications are relatively abundant in bacteria, occurring in about a percent of all bases of most bacteria. Until recently, 5mC and its oxidized derivatives were thought to be the only directed DNA methylation events in metazoa. New and more sensitive detection techniques (ultra-high performance liquid chromatography coupled with mass spectrometry (UHPLC-ms/ms) and single molecule real-time sequencing (SMRTseq)) have suggested that 6mA and 4mC modifications could be present in a variety of metazoa., Results: Here, we find that both of these techniques are prone to inaccuracies, which overestimate DNA methylation concentrations in metazoan genomic DNA. Artifacts can arise from methylated bacterial DNA contamination of enzyme preparations used to digest DNA and contaminating bacterial DNA in eukaryotic DNA preparations. Moreover, DNA sonication introduces a novel modified base from 5mC that has a retention time near 4mC that can be confused with 4mC. Our analyses also suggest that SMRTseq systematically overestimates 4mC in prokaryotic and eukaryotic DNA and 6mA in DNA samples in which it is rare. Using UHPLC-ms/ms designed to minimize and subtract artifacts, we find low to undetectable levels of 4mC and 6mA in genomes of representative worms, insects, amphibians, birds, rodents and primates under normal growth conditions. We also find that mammalian cells incorporate exogenous methylated nucleosides into their genome, suggesting that a portion of 6mA modifications could derive from incorporation of nucleosides from bacteria in food or microbiota. However, gDNA samples from gnotobiotic mouse tissues found rare (0.9-3.7 ppm) 6mA modifications above background., Conclusions: Altogether these data demonstrate that 6mA and 4mC are rarer in metazoa than previously reported, and highlight the importance of careful sample preparation and measurement, and need for more accurate sequencing techniques.

Cytosine can undergo modifications, forming 5-methylcytosine (5-mC) and its oxidized products 5-hydroxymethylcytosine (5-hmC), 5-formylcytosine (5-fC) and 5-carboxylcytosine (5-caC). Despite their importance as epigenetic markers and as central players in cellular processes, it is not well understood how these modifications influence physical properties of DNA and chromatin. Here we report a comprehensive survey of the effect of cytosine modifications on DNA flexibility. We find that even a single copy of 5-fC increases DNA flexibility markedly. 5-mC reduces and 5-hmC enhances flexibility, and 5-caC does not have a measurable effect. Molecular dynamics simulations show that these modifications promote or dampen structural fluctuations, likely through competing effects of base polarity and steric hindrance, without changing the average structure. The increase in DNA flexibility increases the mechanical stability of the nucleosome and vice versa, suggesting a gene regulation mechanism where cytosine modifications change the accessibility of nucleosomal DNA through their effects on DNA flexibility.

In the active cytosine demethylation pathway, 5-methylcytosine (5mC) is oxidized sequentially to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). Thymine DNA glycosylase (TDG) selectively excises 5fC and 5caC but not cytosine (C), 5mC, and 5hmC. We propose that the electron-withdrawing properties of -CHO and -COOH in 5fC and 5caC increase N3 acidity, leading to weakened hydrogen bonding and reduced base pair stability relative to C, 5mC, and 5hmC, thereby facilitating the selective recognition of 5fC and 5caC by TDG. Through (13)C NMR, we measured the pKa at N3 of 5fC as 2.4 and the two pKa's of 5caC as 2.1 and 4.2. We used isotope-edited IR spectroscopy coupled with density functional theory (DFT) calculations to site-specifically assign the more acidic pKa of 5caC to protonation at N3, indicating that N3 acidity is increased in 5fC and 5caC relative to C. IR and UV melting studies of self-complementary DNA oligomers confirm reduced stability for 5fC-G and 5caC-G base pairs. Furthermore, while the 5fC-G base pair stability is insensitive to pH, the 5caC-G stability is reduced as pH decreases and the carboxyl group is increasingly protonated. Despite suggestions that 5fC and 5caC may exist in rare tautomeric structures which form wobble GC base pairs, our two-dimensional infrared (2D IR) spectroscopy of 5fC and 5caC free nucleosides confirms that both bases are predominantly in the canonical amino-keto form. Taken together, these findings support our model that weakened base pairing ability for 5fC and 5caC in dsDNA contributes to their selective recognition by TDG.

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.

Restriction-modification (R-M) systems pose a major barrier to DNA transformation and genetic engineering of bacterial species. Systematic identification of DNA methylation in R-M systems, including N(6)-methyladenine (6mA), 5-methylcytosine (5mC) and N(4)-methylcytosine (4mC), will enable strategies to make these species genetically tractable. Although single-molecule, real time (SMRT) sequencing technology is capable of detecting 4mC directly for any bacterial species regardless of whether an assembled genome exists or not, it is not as scalable to profiling hundreds to thousands of samples compared with the commonly used next-generation sequencing technologies. Here, we present 4mC-Tet-assisted bisulfite-sequencing (4mC-TAB-seq), a next-generation sequencing method that rapidly and cost efficiently reveals the genome-wide locations of 4mC for bacterial species with an available assembled reference genome. In 4mC-TAB-seq, both cytosines and 5mCs are read out as thymines, whereas only 4mCs are read out as cytosines, revealing their specific positions throughout the genome. We applied 4mC-TAB-seq to study the methylation of a member of the hyperthermophilc genus, Caldicellulosiruptor, in which 4mC-related restriction is a major barrier to DNA transformation from other species. In combination with MethylC-seq, both 4mC- and 5mC-containing motifs are identified which can assist in rapid and efficient genetic engineering of these bacteria in the future., (© Published by Oxford University Press on behalf of Nucleic Acids Research 2015. This work is written by (a) US Government employee(s) and is in the public domain in the US.)

Active DNA demethylation in mammals involves oxidation of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC). However, genome-wide detection of 5fC at single-base resolution remains challenging. Here we present fC-CET, a bisulfite-free method for whole-genome analysis of 5fC based on selective chemical labeling of 5fC and subsequent C-to-T transition during PCR. Base-resolution 5fC maps showed limited overlap with 5hmC, with 5fC-marked regions more active than 5hmC-marked ones.

DNA methylation at selective cytosine residues (5-methylcytosine (5mC)) and their removal by TET-mediated DNA demethylation are critical for setting up pluripotent states in early embryonic development. TET enzymes successively convert 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC), with 5fC and 5caC subject to removal by thymine DNA glycosylase (TDG) in conjunction with base excision repair. Early reports indicate that 5fC and 5caC could be stably detected on enhancers, promoters and gene bodies, with distinct effects on gene expression, but the mechanisms have remained elusive. Here we determined the X-ray crystal structure of yeast elongating RNA polymerase II (Pol II) in complex with a DNA template containing oxidized 5mCs, revealing specific hydrogen bonds between the 5-carboxyl group of 5caC and the conserved epi-DNA recognition loop in the polymerase. This causes a positional shift for incoming nucleoside 5'-triphosphate (NTP), thus compromising nucleotide addition. To test the implication of this structural insight in vivo, we determined the global effect of increased 5fC/5caC levels on transcription, finding that such DNA modifications indeed retarded Pol II elongation on gene bodies. These results demonstrate the functional impact of oxidized 5mCs on gene expression and suggest a novel role for Pol II as a specific and direct epigenetic sensor during transcription elongation.

5-hydroxymethylcytosine (5hmC), an oxidized derivative of 5-methylcytosine (5mC), has been implicated as an important epigenetic regulator of mammalian development. Current procedures use DNA sequencing methods to discriminate 5hmC from 5mC, limiting their accessibility to the scientific community. Here we report a method that combines TET-assisted bisulfite conversion with Illumina 450K DNA methylation arrays for a low-cost high-throughput approach that distinguishes 5hmC and 5mC signals at base resolution. Implementing this approach, termed "TAB-array", we assessed DNA methylation dynamics in the differentiation of human pluripotent stem cells into cardiovascular progenitors and neural precursor cells. With the ability to discriminate 5mC and 5hmC, we identified a large number of novel dynamically methylated genomic regions that are implicated in the development of these lineages. The increased resolution and accuracy afforded by this approach provides a powerful means to investigate the distinct contributions of 5mC and 5hmC in human development and disease., (Copyright © 2014 Elsevier Inc. All rights reserved.)

DNA methylation (5-methylcytosine, 5mC) plays critical biological functions in mammals and plants as a vital epigenetic marker. The Ten-Eleven translocation dioxygenases (TET1, 2, and 3) have been found to oxidize 5mC to 5-hydroxymethylcytosine (5hmC) and then to 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC) in mammalian cells. We report herein three mushroom TET homologues from Coprinopsis cinerea that can mediate 5mC oxidation. Specifically, one homologue (CC1G_05589, CcTET) shows similar activity to its mammalian TET homologues. Biochemically, CcTET actively converts 5mC to 5hmC, 5fC, and 5caC under natural conditions (pH 7.0). Interestingly, CcTET also converts the majority of 5mC to 5fC under slightly acidic (pH 5.8) and neutral conditions. Kinetics analyses of the oxidation by CcTET under neutral conditions indicate that conversion of 5mC to 5hmC and 5hmC to 5fC are faster than that of 5fC to 5caC, respectively. Our results provide an example of a TET homologue in a non-mammalian organism that exhibits full 5mC-to-5caC oxidation activity and a slight preference to producing 5fC. The preferential accumulation of 5fC in the in vitro oxidation reactions under both neutral and acidic conditions may have biological implications for 5mC oxidation in fungi species.

Background: 5-methylcytosine (mC) can be oxidized by the tet methylcytosine dioxygenase (Tet) family of enzymes to 5-hydroxymethylcytosine (hmC), which is an intermediate of mC demethylation and may also be a stable epigenetic modification that influences chromatin structure. hmC is particularly abundant in mammalian brains but its function is currently unknown. A high-resolution hydroxymethylome map is required to fully understand the function of hmC in the human brain., Results: We present genome-wide and single-base resolution maps of hmC and mC in the human brain by combined application of Tet-assisted bisulfite sequencing and bisulfite sequencing. We demonstrate that hmCs increase markedly from the fetal to the adult stage, and in the adult brain, 13% of all CpGs are highly hydroxymethylated with strong enrichment at genic regions and distal regulatory elements. Notably, hmC peaks are identified at the 5'splicing sites at the exon-intron boundary, suggesting a mechanistic link between hmC and splicing. We report a surprising transcription-correlated hmC bias toward the sense strand and an mC bias toward the antisense strand of gene bodies. Furthermore, hmC is negatively correlated with H3K27me3-marked and H3K9me3-marked repressive genomic regions, and is more enriched at poised enhancers than active enhancers., Conclusions: We provide single-base resolution hmC and mC maps in the human brain and our data imply novel roles of hmC in regulating splicing and gene expression. Hydroxymethylation is the main modification status for a large portion of CpGs situated at poised enhancers and actively transcribed regions, suggesting its roles in epigenetic tuning at these regions.

DNA methylation is implicated in mammalian brain development and plasticity underlying learning and memory. We report the genome-wide composition, patterning, cell specificity, and dynamics of DNA methylation at single-base resolution in human and mouse frontal cortex throughout their lifespan. Widespread methylome reconfiguration occurs during fetal to young adult development, coincident with synaptogenesis. During this period, highly conserved non-CG methylation (mCH) accumulates in neurons, but not glia, to become the dominant form of methylation in the human neuronal genome. Moreover, we found an mCH signature that identifies genes escaping X-chromosome inactivation. Last, whole-genome single-base resolution 5-hydroxymethylcytosine (hmC) maps revealed that hmC marks fetal brain cell genomes at putative regulatory regions that are CG-demethylated and activated in the adult brain and that CG demethylation at these hmC-poised loci depends on Tet2 activity.

Labeling chemically for a "single" result: A one-step, nonenzymatic 5hmC thiolation method mediated by bisulfite is highlighted. The 5hmC base can be chemically labeled in single-stranded DNA and detected by using nanopore sequencing., (Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

5-Methylcytosine (5mC) in DNA can be oxidized stepwise to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) by the TET family proteins. Thymine DNA glycosylase can further remove 5fC and 5caC, connecting 5mC oxidation with active DNA demethylation. Here, we present a chemical modification-assisted bisulfite sequencing (CAB-Seq) that can detect 5caC with single-base resolution in DNA. We optimized 1-ethyl-3-[3-dimethylaminopropyl]carbodiimide hydrochloride (EDC)-catalyzed amide bond formation between the carboxyl group of 5caC and a primary amine group. We found that the modified 5caC can survive the bisulfite treatment without deamination. Therefore, this chemical labeling coupled with bisulfite treatment provides a base-resolution detection and sequencing method for 5caC.

Mammalian somatic cells can be directly reprogrammed into induced pluripotent stem cells (iPSCs) by introducing defined sets of transcription factors. Somatic cell reprogramming involves epigenomic reconfiguration, conferring iPSCs with characteristics similar to embryonic stem cells (ESCs). Human ESCs (hESCs) contain 5-hydroxymethylcytosine (5hmC), which is generated through the oxidation of 5-methylcytosine by the TET enzyme family. Here we show that 5hmC levels increase significantly during reprogramming to human iPSCs mainly owing to TET1 activation, and this hydroxymethylation change is critical for optimal epigenetic reprogramming, but does not compromise primed pluripotency. Compared with hESCs, we find that iPSCs tend to form large-scale (100 kb-1.3 Mb) aberrant reprogramming hotspots in subtelomeric regions, most of which exhibit incomplete hydroxymethylation on CG sites. Strikingly, these 5hmC aberrant hotspots largely coincide (~80%) with aberrant iPSC-ESC non-CG methylation regions. Our results suggest that TET1-mediated 5hmC modification could contribute to the epigenetic variation of iPSCs and iPSC-hESC differences.

TET proteins oxidize 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). 5fC and 5caC are excised by mammalian DNA glycosylase TDG, implicating 5mC oxidation in DNA demethylation. Here, we show that the genomic locations of 5fC can be determined by coupling chemical reduction with biotin tagging. Genome-wide mapping of 5fC in mouse embryonic stem cells (mESCs) reveals that 5fC preferentially occurs at poised enhancers among other gene regulatory elements. Application to Tdg null mESCs further suggests that 5fC production coordinates with p300 in remodeling epigenetic states of enhancers. This process, which is not influenced by 5hmC, appears to be associated with further oxidation of 5hmC and commitment to demethylation through 5fC. Finally, we resolved 5fC at base resolution by hydroxylamine-based protection from bisulfite-mediated deamination, thereby confirming sites of 5fC accumulation. Our results reveal roles of active 5mC/5hmC oxidation and TDG-mediated demethylation in epigenetic tuning at regulatory elements., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Little is known about how patterns of DNA methylation change during mammalian spermatogenesis. 5 hmC has been recognized as a stable intermediate of DNA demethylation with potential regulatory functions in the mammalian genome. However, its global pattern in germ cells has yet to be addressed. Here, we first conducted absolute quantification of 5 hmC in eight consecutive types of mouse spermatogenic cells using liquid chromatography-tandem mass spectrometry, and then mapped its distributions in various genomic regions using our chemical labeling and enrichment method coupled with deep sequencing. We found that 5 hmC mapped differentially to and changed dynamically in genomic regions related to expression regulation of protein-coding genes, piRNA precursor genes and repetitive elements. Moreover, 5 hmC content correlated with the levels of various transcripts quantified by RNA-seq. These results suggest that the highly ordered alterations of 5 hmC in the mouse genome are potentially crucial for the differentiation of spermatogenic cells.

5-Hydroxymethylcytosine (5-hmC) is a newly discovered modified form of cytosine that has been suspected to be an important epigenetic modification in neurodevelopment. While DNA methylation dynamics have already been implicated during neurodevelopment, little is known about hydroxymethylation in this process. Here, we report DNA hydroxymethylation dynamics during cerebellum development in the human brain. Overall, we find a positive correlation between 5-hmC levels and cerebellum development. Genome-wide profiling reveals that 5-hmC is highly enriched on specific gene regions including exons and especially the untranslated regions (UTRs), but it is depleted on introns and intergenic regions. Furthermore, we have identified fetus-specific and adult-specific differentially hydroxymethylated regions (DhMRs), most of which overlap with genes and CpG island shores. Surprisingly, during development, DhMRs are highly enriched in genes encoding mRNAs that can be regulated by fragile X mental retardation protein (FMRP), some of which are disrupted in autism, as well as in many known autism genes. Our results suggest that 5-hmC-mediated epigenetic regulation may broadly impact the development of the human brain, and its dysregulation could contribute to the molecular pathogenesis of neurodevelopmental disorders. Accession number: Sequencing data have been deposited to GEO with accession number GSE40539.

A complete understanding of the potential function of 5-hydroxymethylcytosine (5-hmC), a DNA cytosine modification in mammalian cells, requires an accurate single-base resolution sequencing method. Here we describe a modified bisulfite-sequencing method, Tet-assisted bisulfite sequencing (TAB-seq), which can identify 5-hmC at single-base resolution, as well as determine its abundance at each modification site. This protocol involves β-glucosyltransferase (β-GT)-mediated protection of 5-hmC (glucosylation) and recombinant mouse Tet1(mTet1)-mediated oxidation of 5-methylcytosine (5-mC) to 5-carboxylcytosine (5-caC). After the subsequent bisulfite treatment and PCR amplification, both cytosine and 5-caC (derived from 5-mC) are converted to thymine (T), whereas 5-hmC reads as C. The treated genomic DNA is suitable for both whole-genome and locus-specific sequencing. The entire procedure (which does not include data analysis) can be completed in 14 d for whole-genome sequencing or 7 d for locus-specific sequencing.

Genome-wide 5-hydroxymethylome analysis of a rodent hepatocarcinogen model reveals that 5-hydroxymethylcytosine-dependent active DNA demethylation may be functionally important in the early stages of carcinogenesis.

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(1). 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(11). 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(12). 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.

Although the roles of 5-methylcytosine and 5-hydroxymethylcytosine in epigenetic regulation of gene expression are well established, the functional effects of 5-formylcytosine and 5-carboxylcytosine on the process of transcription are not clear. Here we report a systematic study of the effects of five different forms of cytosine in DNA on mammalian and yeast RNA polymerase II transcription, providing new insights into potential functional interplay between cytosine methylation status and transcription.

The study of 5-hydroxylmethylcytosines (5hmC) has been hampered by the lack of a method to map it at single-base resolution on a genome-wide scale. Affinity purification-based methods cannot precisely locate 5hmC nor accurately determine its relative abundance at each modified site. We here present a genome-wide approach, Tet-assisted bisulfite sequencing (TAB-Seq), that when combined with traditional bisulfite sequencing can be used for mapping 5hmC at base resolution and quantifying the relative abundance of 5hmC as well as 5mC. Application of this method to embryonic stem cells not only confirms widespread distribution of 5hmC in the mammalian genome but also reveals sequence bias and strand asymmetry at 5hmC sites. We observe high levels of 5hmC and reciprocally low levels of 5mC near but not on transcription factor-binding sites. Additionally, the relative abundance of 5hmC varies significantly among distinct functional sequence elements, suggesting different mechanisms for 5hmC deposition and maintenance., (Copyright © 2012 Elsevier Inc. All rights reserved.)

The discovery of 5-hydroxymethyl-cytosine (5hmC) in mammalian cells prompted us to look for this base in the DNA of Arabidopsis thaliana (thale cress), and to ask how well the Arabidopsis Variant in Methylation 1 (VIM1) protein, an essential factor in maintaining 5-cytosine methylation (5mC) homeostasis and epigenetic silencing in this plant, recognizes this novel base. We found that the DNA of Arabidopsis' leaves and flowers contain low levels of 5hmC. We also cloned and expressed in Escherichia coli full-length VIM1 protein, the archetypal member of the five Arabidopsis VIM gene family. Using in vitro binding assays, we observed that full-length VIM1 binds preferentially to hemi-methylated DNA with a single modified 5mCpG site; this result is consistent with its known role in preserving DNA methylation in vivo following DNA replication. However, when 5hmC replaces one or both cytosine residues at a palindromic CpG site, VIM1 binds with approximately ≥10-fold lower affinity. These results suggest that 5hmC may contribute to VIM-mediated passive loss of cytosine methylation in vivo during Arabidopsis DNA replication., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Human thymine DNA glycosylase (hTDG) efficiently excises 5-carboxylcytosine (5caC), a key oxidation product of 5-methylcytosine in genomic DNA, in a recently discovered cytosine demethylation pathway. We present here the crystal structures of the hTDG catalytic domain in complex with duplex DNA containing either 5caC or a fluorinated analog. These structures, together with biochemical and computational analyses, reveal that 5caC is specifically recognized in the active site of hTDG, supporting the role of TDG in mammalian 5-methylcytosine demethylation.

One of the recent advances in the epigenetic field is the demonstration that the Tet family of proteins are capable of catalyzing conversion of 5-methylcytosine (5mC) of DNA to 5-hydroxymethylcytosine (5hmC). Interestingly, recent studies have shown that 5hmC can be further oxidized by Tet proteins to generate 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC), which can be removed by thymine DNA glycosylase (TDG). To determine whether Tet-catalyzed conversion of 5mC to 5fC and 5caC occurs in vivo in zygotes, we generated antibodies specific for 5fC and 5caC. By immunostaining, we demonstrate that loss of 5mC in the paternal pronucleus is concurrent with the appearance of 5fC and 5caC, similar to that of 5hmC. Importantly, instead of being quickly removed through an enzyme-catalyzed process, both 5fC and 5caC exhibit replication-dependent dilution during mouse preimplantation development. These results not only demonstrate the conversion of 5mC to 5fC and 5caC in zygotes, but also indicate that both 5fC and 5caC are relatively stable and may be functional during preimplantation development. Together with previous studies, our study suggests that Tet-catalyzed conversion of 5mC to 5hmC/5fC/5caC followed by replication-dependent dilution accounts for paternal DNA demethylation during preimplantation development.

We describe strand-specific, base-resolution detection of 5-hydroxymethylcytosine (5-hmC) in genomic DNA with single-molecule sensitivity, combining a bioorthogonal, selective chemical labeling method of 5-hmC with single-molecule, real-time (SMRT) DNA sequencing. The chemical labeling not only allows affinity enrichment of 5-hmC-containing DNA fragments but also enhances the kinetic signal of 5-hmC during SMRT sequencing. We applied the approach to sequence 5-hmC in a genomic DNA sample with high confidence.

DNA methylation dynamics influence brain function and are altered in neurological disorders. 5-hydroxymethylcytosine (5-hmC), a DNA base that is derived from 5-methylcytosine, accounts for ∼40% of modified cytosine in the brain and has been implicated in DNA methylation-related plasticity. We mapped 5-hmC genome-wide in mouse hippocampus and cerebellum at three different ages, which allowed us to assess its stability and dynamic regulation during postnatal neurodevelopment through adulthood. We found developmentally programmed acquisition of 5-hmC in neuronal cells. Epigenomic localization of 5-hmC-regulated regions revealed stable and dynamically modified loci during neurodevelopment and aging. By profiling 5-hmC in human cerebellum, we found conserved genomic features of 5-hmC. Finally, we found that 5-hmC levels were inversely correlated with methyl-CpG-binding protein 2 dosage, a protein encoded by a gene in which mutations cause Rett syndrome. These data suggest that 5-hmC-mediated epigenetic modification is critical in neurodevelopment and diseases.

DNA is not merely a combination of four genetic codes, namely A, T, C, and G. It also contains minor modifications that play crucial roles throughout biology. For example, the fifth DNA base, 5-methylcytosine (5-mC), which accounts for ∼1% of all the nucleotides in mammalian genomic DNA, is a vital epigenetic mark. It impacts a broad range of biological functions, from development to cancer. Recently, an oxidized form of 5-methylcytosine, 5-hydroxymethylcytosine (5-hmC), was found to constitute the sixth base in the mammalian genome; it was believed to be another crucial epigenetic mark. Unfortunately, further study of this newly discovered DNA base modification has been hampered by inadequate detection and sequencing methods, because current techniques fail to differentiate 5-hmC from 5-mC. The immediate challenge, therefore, is to develop robust methods for ascertaining the positions of 5-hmC within the mammalian genome. In this Account, we describe our development of the first bioorthogonal, selective labeling of 5-hmC to specifically address this challenge. We utilize β-glucosyltransferase (βGT) to transfer an azide-modified glucose onto 5-hmC in genomic DNA. The azide moiety enables further bioorthogonal click chemistry to install a biotin group, which allows for detection, affinity enrichment, and, most importantly, deep sequencing of the 5-hmC-containing DNA. With this highly effective and selective method, we revealed the first genome-wide distribution of 5-hmC in the mouse genome and began to shed further light on the biology of 5-hmC. The strategy lays the foundation for developing high-throughput, single-base-resolution sequencing methods for 5-hmC in mammalian genomes in the future. DNA and RNA are not static inside cells. They interact with protein and other DNA and RNA in fundamental biological processes such as replication, transcription, translation, and DNA and RNA modification and repair. The ability to investigate these interactions will also be enhanced by developing and utilizing bioorthogonal probes. We have chosen the photoreactive diazirine photophore as a bioorthogonal moiety to develop nucleic acid probes. The small size and unique photo-cross-linking activity of diazirine enabled us to develop a series of novel cross-linking probes to streamline the study of protein-nucleic acid and nucleic acid-nucleic acid interactions. In the second half of this Account, we highlight a few examples of these probes., (© 2011 American Chemical Society)

The prevalent DNA modification in higher organisms is the methylation of cytosine to 5-methylcytosine (5mC), which is partially converted to 5-hydroxymethylcytosine (5hmC) by the Tet (ten eleven translocation) family of dioxygenases. Despite their importance in epigenetic regulation, it is unclear how these cytosine modifications are reversed. Here, we demonstrate that 5mC and 5hmC in DNA are oxidized to 5-carboxylcytosine (5caC) by Tet dioxygenases in vitro and in cultured cells. 5caC is specifically recognized and excised by thymine-DNA glycosylase (TDG). Depletion of TDG in mouse embyronic stem cells leads to accumulation of 5caC to a readily detectable level. These data suggest that oxidation of 5mC by Tet proteins followed by TDG-mediated base excision of 5caC constitutes a pathway for active DNA demethylation.

5-methylcytosine (5mC) in DNA plays an important role in gene expression, genomic imprinting, and suppression of transposable elements. 5mC can be converted to 5-hydroxymethylcytosine (5hmC) by the Tet (ten eleven translocation) proteins. Here, we show that, in addition to 5hmC, the Tet proteins can generate 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC) from 5mC in an enzymatic activity-dependent manner. Furthermore, we reveal the presence of 5fC and 5caC in genomic DNA of mouse embryonic stem cells and mouse organs. The genomic content of 5hmC, 5fC, and 5caC can be increased or reduced through overexpression or depletion of Tet proteins. Thus, we identify two previously unknown cytosine derivatives in genomic DNA as the products of Tet proteins. Our study raises the possibility that DNA demethylation may occur through Tet-catalyzed oxidation followed by decarboxylation.

5-Hydroxymethylcytosine (5-hmC) is a newly discovered DNA base in mammalian cells that is believed to be another important epigenetic modification. Here we report the use of a methylation-insensitive restriction enzyme Taq(α)I coupled with selective chemical labeling of 5-hmC in a combined glycosylation restriction analysis (CGRA) to detect 5-hmC in TCGA sequences. This method, differentiates fully versus hemi-hydroxymethylated cytosine in the CpG dinucleotide, adds a new tool to facilitate biological studies of 5-hmC., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

To investigate the potential oxidation products of 5-hydroxymethylcytosine (5-hmC)-containing DNA, we present here efficient syntheses of 5-formyl- and 5-methoxycarbonyl-2'-deoxycytidine phosphoramidites. The 5-formyl group in III was easy to introduce and was compatible with phosphoramidite and DNA syntheses. An additional treatment of ODN1 with NaBH(4) produced the corresponding ODN2 quantitatively. Phosphoramidite V was also incorporated into DNA, and the methyl ester could be hydrolyzed under mild basic conditions to afford ODN3.

Covalent modification of DNA distinguishes cellular identities and is crucial for regulating the pluripotency and differentiation of embryonic stem (ES) cells. The recent demonstration that 5-methylcytosine (5-mC) may be further modified to 5-hydroxymethylcytosine (5-hmC) in ES cells has revealed a novel regulatory paradigm to modulate the epigenetic landscape of pluripotency. To understand the role of 5-hmC in the epigenomic landscape of pluripotent cells, here we profile the genome-wide 5-hmC distribution and correlate it with the genomic profiles of 11 diverse histone modifications and six transcription factors in human ES cells. By integrating genomic 5-hmC signals with maps of histone enrichment, we link particular pluripotency-associated chromatin contexts with 5-hmC. Intriguingly, through additional correlations with defined chromatin signatures at promoter and enhancer subtypes, we show distinct enrichment of 5-hmC at enhancers marked with H3K4me1 and H3K27ac. These results suggest potential role(s) for 5-hmC in the regulation of specific promoters and enhancers. In addition, our results provide a detailed epigenomic map of 5-hmC from which to pursue future functional studies on the diverse regulatory roles associated with 5-hmC., Competing Interests: The authors have declared that no competing interests exist.

5-Hydroxymethylcytosine (5-hmC) is a newly discovered DNA base modification in mammalian genomic DNA that is proposed to be a major epigenetic mark. We report here the syntheses of two new versions of phosphoramidites III and IV from 5-iodo-2'-deoxyuridine in 18% and 32% overall yields, respectively, with TBDMS as the 5-hydroxyl protecting group. Phosphoramidites III and IV allow efficient incorporation of 5-hmC into DNA and a "one-step" deprotection procedure to cleanly remove all the protecting groups. A "two-step" deprotection strategy is compatible with ultramild DNA synthesis, which enables the synthesis of 5hmC-containing DNA with additional modifications.

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.

N(1)-meA and N(3)-meC are cytotoxic DNA base methylation lesions that can accumulate in the genomes of various organisms in the presence of S(N)2 type methylating agents. We report here the structural characterization of these base lesions in duplex DNA using a cross-linked protein-DNA crystallization system. The crystal structure of N(1)-meA:T pair shows an unambiguous Hoogsteen base pair with a syn conformation adopted by N(1)-meA, which exhibits significant changes in the opening, roll and twist angles as compared to the normal A:T base pair. Unlike N(1)-meA, N(3)-meC does not establish any interaction with the opposite G, but remains partially intrahelical. Also, structurally characterized is the N(6)-meA base modification that forms a normal base pair with the opposite T in duplex DNA. Structural characterization of these base methylation modifications provides molecular level information on how they affect the overall structure of duplex DNA. In addition, the base pairs containing N(1)-meA or N(3)-meC do not share any specific characteristic properties except that both lesions create thermodynamically unstable regions in a duplex DNA, a property that may be explored by the repair proteins to locate these lesions.

In mammals, cytosine methylation (5mC) is widely distributed throughout the genome but is notably depleted from active promoters and enhancers. While the role of DNA methylation in promoter silencing has been well documented, the function of this epigenetic mark at enhancers remains unclear. Recent experiments have demonstrated that enhancers are enriched for 5-hydroxymethylcytosine (5hmC), an oxidization product of the Tet family of 5mC dioxygenases and an intermediate of DNA demethylation. These results support the involvement of Tet proteins in the regulation of dynamic DNA methylation at enhancers. By mapping DNA methylation and hydroxymethylation at base resolution, we find that deletion of Tet2 causes extensive loss of 5hmC at enhancers, accompanied by enhancer hypermethylation, reduction of enhancer activity, and delayed gene induction in the early steps of differentiation. Our results reveal that DNA demethylation modulates enhancer activity, and its disruption influences the timing of transcriptome reprogramming during cellular differentiation.

Despite explosive growth in genomic datasets, the methods for studying epigenomic mechanisms of gene regulation remain primitive. Here we present a model-based approach to systematically analyze the epigenomic functions in modulating transcription factor-DNA binding. Based on the first principles of statistical mechanics, this model considers the interactions between epigenomic modifications and a cis-regulatory module, which contains multiple binding sites arranged in any configurations. We compiled a comprehensive epigenomic dataset in mouse embryonic stem (mES) cells, including DNA methylation (MeDIP-seq and MRE-seq), DNA hydroxymethylation (5-hmC-seq), and histone modifications (ChIP-seq). We discovered correlations of transcription factors (TFs) for specific combinations of epigenomic modifications, which we term epigenomic motifs. Epigenomic motifs explained why some TFs appeared to have different DNA binding motifs derived from in vivo (ChIP-seq) and in vitro experiments. Theoretical analyses suggested that the epigenome can modulate transcriptional noise and boost the cooperativity of weak TF binding sites. ChIP-seq data suggested that epigenomic boost of binding affinities in weak TF binding sites can function in mES cells. We showed in theory that the epigenome should suppress the TF binding differences on SNP-containing binding sites in two people. Using personal data, we identified strong associations between H3K4me2/H3K9ac and the degree of personal differences in NFκB binding in SNP-containing binding sites, which may explain why some SNPs introduce much smaller personal variations on TF binding than other SNPs. In summary, this model presents a powerful approach to analyze the functions of epigenomic modifications. This model was implemented into an open source program APEG (Affinity Prediction by Epigenome and Genome, http://systemsbio.ucsd.edu/apeg).

The reprogramming of cellular memory in specific cell types, and in visceral adipocytes in particular, appears to be a fundamental aspect of obesity and its related negative health outcomes. We explored the hypothesis that adipose tissue contains epigenetically distinct subpopulations of adipocytes that are differentially potentiated to record cellular memories of their environment. Adipocytes are large, fragile, and technically difficult to efficiently isolate and fractionate. We developed fluorescence nuclear cytometry (FNC) and fluorescence activated nuclear sorting (FANS) of cellular nuclei from visceral adipose tissue (VAT) using the levels of the pan-adipocyte protein, peroxisome proliferator-activated receptor gamma-2 (PPARg2), to distinguish classes of PPARg2-Positive (PPARg2-Pos) adipocyte nuclei from PPARg2-Negative (PPARg2-Neg) leukocyte and endothelial cell nuclei. PPARg2-Pos nuclei were 10-fold enriched for most adipocyte marker transcripts relative to PPARg2-Neg nuclei. PPARg2-Pos nuclei showed 2- to 50-fold higher levels of transcripts encoding most of the chromatin-remodeling factors assayed, which regulate the methylation of histones and DNA cytosine (e.g., DNMT1, TET1, TET2, KDM4A, KMT2C, SETDB1, PAXIP1, ARID1A, JMJD6, CARM1, and PRMT5). PPARg2-Pos nuclei were large with decondensed chromatin. TAB-seq demonstrated 5-hydroxymethylcytosine (5hmC) levels were remarkably dynamic in gene bodies of various classes of VAT nuclei, dropping 3.8-fold from the highest quintile of expressed genes to the lowest. In short, VAT-derived adipocytes appear to be more actively remodeling their chromatin than non-adipocytes. [ABSTRACT FROM AUTHOR]

Background Peripheral blood leukocytes are the most commonly used surrogates to study epigenome-induced risk and epigenomic response to disease-related stress. We considered the hypothesis that the various classes of peripheral leukocytes differentially regulate the synthesis of 5-methylcytosine (5mCG) and its removal via Ten-Eleven Translocation (TET) dioxygenase catalyzed hydroxymethylation to 5-hydroxymethylcytosine (5hmCG), reflecting their responsiveness to environment. Although it is known that reductions in TET1 and/or TET2 activity lead to the over-proliferation of various leukocyte precursors in bone marrow and in development of chronic myelomonocytic leukemia and myeloproliferative neoplasms, the role of 5mCG hydroxymethylation in peripheral blood is less well studied. Results We developed simplified protocols to rapidly and reiteratively isolate non-overlapping leukocyte populations from a single small sample of fresh or frozen whole blood. Among peripheral leukocyte types we found extreme variation in the levels of transcripts encoding proteins involved in cytosine methylation (DNMT1, 3A, 3B), the turnover of 5mC by demethylation (TET1, 2, 3), and DNA repair (GADD45A, B, G) and in the global and gene-region-specific levels of DNA 5hmCG (CD4 + T cells ≫ CD14 + monocytes > CD16 + neutrophils > CD19 + B cells > CD56 + NK cells > Siglec8 + eosinophils > CD8 + T cells). Conclusions Our data taken together suggest a potential hierarchy of responsiveness among classes of leukocytes with CD4 +, CD8 + T cells and CD14 + monocytes being the most distinctly poised for a rapid methylome response to physiological stress and disease. [ABSTRACT FROM AUTHOR]

Summary: The reprogramming of parental methylomes is essential for embryonic development. In mammals, paternal 5-methylcytosines (5mCs) have been proposed to be actively converted to oxidized bases. These paternal oxidized bases and maternal 5mCs are believed to be passively diluted by cell divisions. By generating single-base resolution, allele-specific DNA methylomes from mouse gametes, early embryos, and primordial germ cell (PGC), as well as single-base-resolution maps of oxidized cytosine bases for early embryos, we report the existence of 5hmC and 5fC in both maternal and paternal genomes and find that 5mC or its oxidized derivatives, at the majority of demethylated CpGs, are converted to unmodified cytosines independent of passive dilution from gametes to four-cell embryos. Therefore, we conclude that paternal methylome and at least a significant proportion of maternal methylome go through active demethylation during embryonic development. Additionally, all the known imprinting control regions (ICRs) were classified into germ-line or somatic ICRs. [ABSTRACT FROM AUTHOR]

Abstract: As an important step of the active demethylation of 5-methylcytosine (5mC), human thymine DNA glycosylase (hTDG) efficiently excises 5-carboxylcytosine (5caC) from double-stranded DNA (dsDNA). Here, we present synthesis of DNA oligos containing a 2′-deoxy-2′-fluoro-d-arabinofuranosyl-5-carboxylcytidine (F-5caC) modification that act as hTDG inhibitors. The glycosylase activity assay showed that F-5caC oligos were resistant to excision by the hTDG catalytic domain (hTDGcat, residues 111–308) and they could inhibit the excision of DNA oligos containing 5caC. The electrophoretic mobility shift assay confirmed that DNA oligos containing F-5caC could bind well with unmodified hTDGcat to form a stable complex, which makes it possible to obtain the crystal structure of the complex to reveal details on how hTDGcat recognizes the DNA substrate. [Copyright &y& Elsevier]