5-Hydroxymethylcytosine: the many faces of the sixth base of mammalian DNA.

Epigenetic phenomena play a central role in cell regulatory processes and are important factors for understanding complex human disease. One of the best understood epigenetic mechanisms is DNA methylation. In the mammalian genome, cytosines (C) in CpG dinucleotides were long known to undergo methylation at the 5-position of the pyrimidine ring (mC). Later it was found that mC can be oxidized to 5-hydroxymethylcytosine (hmC) or even further to 5-formylcytosine (fC) and to 5-carboxylcytosine (caC) by the action of 2-oxoglutarate-dependent dioxygenases of the TET family. These findings unveiled a long elusive mechanism of active DNA demethylation and bolstered a wave of studies in the area of epigenetic regulation in mammals. This review is dedicated to critical assessment of recent data on biochemical and chemical aspects of the formation and conversion of hmC in DNA, analytical techniques used for detection and mapping of this nucleobase in mammalian genomes as well as epigenetic roles of hmC in DNA replication, transcription, cell differentiation and human disease.

One-pot trimodal mapping of unmethylated, hydroxymethylated, and open chromatin sites unveils distinctive 5hmC roles at dynamic chromatin loci.

We present a method, named Mx-TOP, for profiling of three epigenetic regulatory layers-chromatin accessibility, general DNA modification, and DNA hydroxymethylation-from a single library. The approach is based on chemo-enzymatic covalent tagging of unmodified CG sites and hydroxymethylated cytosine (5hmC) along with GC sites in chromatin, which are then mapped using tag-selective base-resolution TOP-seq sequencing. Our in-depth validation of the approach revealed its sensitivity and informativity in evaluating chromatin accessibility and DNA modification interactions that drive transcriptional regulation. We employed the technology in a study of chromatin and DNA demethylation dynamics during in vitro neuronal differentiation. The study highlighted the involvement of gene body 5hmC in modulating an extensive decoupling between promoter accessibility and transcription. The importance of 5hmC in chromatin remodeling was further demonstrated by the observed resistance of the developmentally acquired open loci to the global 5hmC erasure in neuronal progenitors.

Chemical Expansion of the Methyltransferase Reaction: Tools for DNA Labeling and Epigenome Analysis.

ConspectusDNA is the genetic matter of life composed of four major nucleotides which can be further furnished with biologically important covalent modifications. Among the variety of enzymes involved in DNA metabolism, AdoMet-dependent methyltransferases (MTases) combine the recognition of specific sequences and covalent methylation of a target nucleotide. The naturally transferred methyl groups play important roles in biological signaling, but they are poor physical reporters and largely resistant to chemical derivatization. Therefore, an obvious strategy to unlock the practical utility of the methyltransferase reactions is to enable the transfer of "prederivatized" (extended) versions of the methyl group.However, previous enzymatic studies of extended AdoMet analogs indicated that the transalkylation reactions are drastically impaired as the size of the carbon chain increases. In collaborative efforts, we proposed that, akin to enhanced SN2 reactivity of allylic and propargylic systems, addition of a π orbital next to the transferable carbon atom might confer the needed activation of the reaction. Indeed, we found that MTase-catalyzed transalkylations of DNA with cofactors containing a double or a triple C-C bond in the ß position occurred in a robust and sequence-specific manner. Altogether, this breakthrough approach named mTAG (methyltransferase-directed transfer of activated groups) has proven instrumental for targeted labeling of DNA and other types of biomolecules (using appropriate MTases) including RNA and proteins.Our further work focused on the propargylic cofactors and their reactions with DNA cytosine-5 MTases, a class of MTases common for both prokaryotes and eukaryotes. Here, we learned that the 4-X-but-2-yn-1-yl (X = polar group) cofactors suffered from a rapid loss of activity in aqueous buffers due to susceptibility of the triple bond to hydration. This problem was remedied by synthetically increasing the separation between X and the triple bond from one to three carbon units (6-X-hex-2-ynyl cofactors). To further optimize the transfer of the bulkier groups, we performed structure-guided engineering of the MTase cofactor pocket. Alanine replacements of two conserved residues conferred substantial improvements of the transalkylation activity with M.HhaI and three other engineered bacterial C5-MTases. Of particular interest were CpG-specific DNA MTases (M.SssI), which proved valuable tools for studies of mammalian methylomes and chemical probing of DNA function.Inspired by the successful repurposing of bacterial enzymes, we turned to more complex mammalian C5-MTases (Dnmt1, Dnmt3A, and Dnmt3B) and asked if they could ultimately lead to mTAG labeling inside mammalian cells. Our efforts to engineer mouse Dnmt1 produced a variant (Dnmt1*) that enabled efficient Dnmt1-directed deposition of 6-azide-hexynyl groups on DNA in vitro. CRISPR-Cas9 editing of the corresponding codons in the genomic Dnmt1 alleles established endogenous expression of Dnmt1* in mouse embryonic stem cells. To circumvent the poor cellular uptake of AdoMet and its analogs, we elaborated their efficient internalization by electroporation, which has finally enabled selective catalysis-dependent azide tagging of natural Dnmt1 targets in live mammalian cells. The deposited chemical groups were then exploited as "click" handles for reading adjoining sequences and precise genomic mapping of the methylation sites. These findings offer unprecedented inroads into studies of DNA methylation in a wide range of eukaryotic model systems.

DNA Labeling Using DNA Methyltransferases.

DNA methyltransferases (MTases) uniquely combine the ability to recognize and covalently modify specific target sequences in DNA using the ubiquitous cofactor S-Adenosyl-L-methionine (AdoMet). Although DNA methylation plays important roles in biological signaling, the transferred methyl group is a poor reporter and is highly inert to further biocompatible derivatization. To unlock the biotechnological power of these enzymes, extended cofactor AdoMet analogs have been developed that enable targeted MTase-directed attachment of larger moieties containing functional or reporter groups onto DNA. As the enlarged cofactors are not always compatible with the active sites of native MTases, steric engineering of the active site has been employed to optimize their alkyltransferase activity. In addition to the described cofactor analogs, recently discovered atypical reactions of DNA cytosine-5 MTases involving non-cofactor-like compounds can also be exploited for targeted derivatization and labeling of DNA. Altogether, these approaches offer new powerful tools for sequence-specific covalent DNA labeling, leading to a variety of useful techniques in DNA research, diagnostics and nanotechnologies, and have already proven practical utility for optical DNA mapping and high-throughput epigenome studies.

Distribution and regulatory roles of oxidized 5-methylcytosines in DNA and RNA of the basidiomycete fungi Laccaria bicolor and Coprinopsis cinerea.

The formation of three oxidative DNA 5-methylcytosine (5mC) modifications (oxi-mCs)-5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC)-by the TET/JBP family of dioxygenases prompted intensive studies of their functional roles in mammalian cells. However, the functional interplay of these less abundant modified nucleotides in other eukaryotic lineages remains poorly understood. We carried out a systematic study of the content and distribution of oxi-mCs in the DNA and RNA of the basidiomycetes Laccaria bicolor and Coprinopsis cinerea, which are established models to study DNA methylation and developmental and symbiotic processes. Quantitative liquid chromatography-tandem mass spectrometry revealed persistent but uneven occurrences of 5hmC, 5fC and 5caC in the DNA and RNA of the two organisms, which could be upregulated by vitamin C. 5caC in RNA (5carC) was predominantly found in non-ribosomal RNA, which potentially includes non-coding, messenger and small RNA species. Genome-wide mapping of 5hmC and 5fC using the single CG analysis techniques hmTOP-seq and foTOP-seq pointed at involvement of oxi-mCs in the regulation of gene expression and silencing of transposable elements. The implicated diverse roles of 5mC and oxi-mCs in the two fungi highlight the epigenetic importance of the latter modifications, which are often neglected in standard whole-genome bisulfite analyses.

Multiomics Analysis of Neuroblastoma Cells Reveals a Diversity of Malignant Transformations.

Neuroblastoma (NB) is a pediatric cancer of the developing sympathetic nervous system that exhibits significant variation in the stage of differentiation and cell composition of tumors. Global loss of DNA methylation and genomic 5-hydroxymethylcytosine (5hmC) is a hallmark of human cancers. Here, we used our recently developed single-base resolution approaches, hmTOP-seq and uTOP-seq, for construction of 5hmC maps and identification of large partially methylated domains (PMDs) in different NB cell subpopulations. The 5hmC profiles revealed distinct signatures characteristic to different cell lineages and stages of malignant transformation of NB cells in a conventional and oxygen-depleted environment, which often occurs in tumors. The analysis of the cell-type-specific PMD distribution highlighted differences in global genome organization among NB cells that were ascribed to the same lineage identity by transcriptomic networks. Collectively, we demonstrated a high informativeness of the integrative epigenomic and transcriptomic research and large-scale genome structure in investigating the mechanisms that regulate cell identities and developmental stages of NB cells. Such multiomics analysis, as compared with mutational studies, open new ways for identification of novel disease-associated features which bring prognostic and therapeutic value in treating this aggressive pediatric disease.

Identification of fetal unmodified and 5-hydroxymethylated CG sites in maternal cell-free DNA for non-invasive prenatal testing.

BACKGROUND: Massively parallel sequencing of maternal cell-free DNA (cfDNA) is widely used to test fetal genetic abnormalities in non-invasive prenatal testing (NIPT). However, sequencing-based approaches are still of high cost. Building upon previous knowledge that placenta, the main source of fetal circulating DNA, is hypomethylated in comparison to maternal tissue counterparts of cfDNA, we propose that targeting either unmodified or 5-hydroxymethylated CG sites specifically enriches fetal genetic material and reduces numbers of required analytical sequencing reads thereby decreasing cost of a test. METHODS: We employed uTOPseq and hmTOP-seq approaches which combine covalent derivatization of unmodified or hydroxymethylated CG sites, respectively, with next generation sequencing, or quantitative real-time PCR. RESULTS: We detected increased 5-hydroxymethylcytosine (5hmC) levels in fetal chorionic villi (CV) tissue samples as compared with peripheral blood. Using our previously developed uTOP-seq and hmTOP-seq approaches we obtained whole-genome uCG and 5hmCG maps of 10 CV tissue and 38 cfDNA samples in total. Our results indicated that, in contrast to conventional whole genome sequencing, such epigenomic analysis highly specifically enriches fetal DNA fragments from maternal cfDNA. While both our approaches yielded 100% accuracy in detecting Down syndrome in fetuses, hmTOP-seq maintained such accuracy at ultra-low sequencing depths using only one million reads. We identified 2164 and 1589 placenta-specific differentially modified and 5-hydroxymethylated regions, respectively, in chromosome 21, as well as 3490 and 2002 Down syndrome-specific differentially modified and 5-hydroxymethylated regions, respectively, that can be used as biomarkers for identification of Down syndrome or other epigenetic diseases of a fetus. CONCLUSIONS: uTOP-seq and hmTOP-seq approaches provide a cost-efficient and sensitive epigenetic analysis of fetal abnormalities in maternal cfDNA. The results demonstrated that T21 fetuses contain a perturbed epigenome and also indicated that fetal cfDNA might originate from fetal tissues other than placental chorionic villi. Robust covalent derivatization followed by targeted analysis of fetal DNA by sequencing or qPCR presents an attractive strategy that could help achieve superior sensitivity and specificity in prenatal diagnostics.

A Bisulfite-free Approach for Base-Resolution Analysis of Genomic 5-Carboxylcytosine.

Due to an extreme rarity of 5-carboxylcytosine (5caC) in the mammalian genome, investigation of its role brings a considerable challenge. Methods based on bisulfite sequencing have been proposed for genome-wide 5caC analysis. However, bisulfite-based sequencing of scarcely abundant 5caC demands significant experimental and computational resources, increasing sequencing cost. Here, we present a bisulfite-free approach, caCLEAR, for high-resolution mapping of 5caCGs. The method uses an atypical activity of the methyltransferase eM.SssI to remove a carboxyl group from 5caC, generating unmodified CGs, which are localized by uTOP-seq sequencing. Validation of caCLEAR on model DNA systems and mouse ESCs supports the suitability of caCLEAR for analysis of 5caCGs. The 5caCG profiles of naive and primed pluripotent ESCs reflect their distinct demethylation dynamics and demonstrate an association of 5caC with gene expression. Generally, we demonstrate that caCLEAR is a robust economical approach that could help provide deeper insights into biological roles of 5caC.

Precise genomic mapping of 5-hydroxymethylcytosine via covalent tether-directed sequencing.

5-hydroxymethylcytosine (5hmC) is the most prevalent intermediate on the oxidative DNA demethylation pathway and is implicated in regulation of embryogenesis, neurological processes, and cancerogenesis. Profiling of this relatively scarce genomic modification in clinical samples requires cost-effective high-resolution techniques that avoid harsh chemical treatment. Here, we present a bisulfite-free approach for 5hmC profiling at single-nucleotide resolution, named hmTOP-seq (5hmC-specific tethered oligonucleotide-primed sequencing), which is based on direct sequence readout primed at covalently labeled 5hmC sites from an in situ tethered DNA oligonucleotide. Examination of distinct conjugation chemistries suggested a structural model for the tether-directed nonhomologous polymerase priming enabling theoretical evaluation of suitable tethers at the design stage. The hmTOP-seq procedure was optimized and validated on a small model genome and mouse embryonic stem cells, which allowed construction of single-nucleotide 5hmC maps reflecting subtle differences in strand-specific CG hydroxymethylation. Collectively, hmTOP-seq provides a new valuable tool for cost-effective and precise identification of 5hmC in characterizing its biological role and epigenetic changes associated with human disease.

Epigenomic technologies for deciphering circulating tumor DNA.

Tumor-derived DNA, found in body fluids (liquid biopsy) of cancer patients as part of cell-free DNA (cfDNA), lends itself for noninvasive cancer detection and monitoring. Advantages of cfDNA as analytical target have evoked a burst of sophisticated techniques, providing clinically relevant information. Each cell type carries a unique DNA modification profile consisting mainly of patterns of 5-methylcytosine in CpG dinucleotides, which are critical for establishing and maintaining cellular identity and which are frequently disturbed in cancer. Assessment of the tumor-derived cfDNA modifications combined with high-throughput analysis techniques holds promise for developing highly specific noninvasive diagnostic tests. This review highlights recent advances in locus-specific and whole-genome analysis of cfDNA, with a specific focus on epigenetic phenomena and their clinical value.

Tethered Oligonucleotide-Primed Sequencing, TOP-Seq: A High-Resolution Economical Approach for DNA Epigenome Profiling.

Modification of CG dinucleotides in DNA is part of epigenetic regulation of gene function in vertebrates and is associated with complex human disease. Bisulfite sequencing permits high-resolution analysis of cytosine modification in mammalian genomes; however, its utility is often limited due to substantial cost. Here, we describe an alternative epigenome profiling approach, named TOP-seq, which is based on covalent tagging of individual unmodified CG sites followed by non-homologous priming of the DNA polymerase action at these sites to directly produce adjoining regions for their sequencing and precise genomic mapping. Pilot TOP-seq analyses of bacterial and human genomes showed a better agreement of TOP-seq with published bisulfite sequencing maps as compared to widely used MBD-seq and MRE-seq and permitted identification of long-range and gene-level differential methylation among human tissues and neuroblastoma cell types. Altogether, we propose an affordable single CG-resolution technique well suited for large-scale epigenome studies.

DNA Labeling Using DNA Methyltransferases.

DNA methyltransferases (MTases) uniquely combine the ability to recognize and covalently modify specific target sequences in DNA using the ubiquitous cofactor S-adenosyl-L-methionine (AdoMet). Although DNA methylation plays important roles in biological signaling, the transferred methyl group is a poor reporter and is highly inert to further biocompatible derivatization. To unlock the biotechnological power of these enzymes, two major types of cofactor AdoMet analogs were developed that permit targeted MTase-directed attachment of larger moieties containing functional or reporter groups onto DNA. One such approach (named sequence-specific methyltransferase-induced labeling, SMILing) uses reactive aziridine or N-mustard mimics of the cofactor AdoMet, which render targeted coupling of a whole cofactor molecule to the target DNA. The second approach (methyltransferase-directed transfer of activated groups, mTAG) uses AdoMet analogs with a sulfonium-bound extended side chain replacing the methyl group, which permits MTase-directed covalent transfer of the activated side chain alone. As the enlarged cofactors are not always compatible with the active sites of native MTases, steric engineering of the active site has been employed to optimize their alkyltransferase activity. In addition to the described cofactor analogs, recently discovered atypical reactions of DNA cytosine-5 MTases involving non-cofactor-like compounds can also be exploited for targeted derivatization and labeling of DNA. Altogether, these approaches offer new powerful tools for sequence-specific covalent DNA labeling, which not only pave the way to developing a variety of useful techniques in DNA research, diagnostics, and nanotechnologies but have already proven practical utility for optical DNA mapping and epigenome studies.

Lactase nonpersistence is directed by DNA-variation-dependent epigenetic aging.

The inability to digest lactose, due to lactase nonpersistence, is a common trait in adult mammals, except in certain human populations that exhibit lactase persistence. It is not known how the lactase gene is dramatically downregulated with age in most individuals but remains active in some individuals. We performed a comprehensive epigenetic study of human and mouse small intestines, by using chromosome-wide DNA-modification profiling and targeted bisulfite sequencing. Epigenetically controlled regulatory elements accounted for the differences in lactase mRNA levels among individuals, intestinal cell types and species. We confirmed the importance of these regulatory elements in modulating lactase mRNA levels by using CRISPR-Cas9-induced deletions. Genetic factors contribute to epigenetic changes occurring with age at the regulatory elements, because lactase-persistence and lactase-nonpersistence DNA haplotypes demonstrated markedly different epigenetic aging. Thus, genetic factors enable a gradual accumulation of epigenetic changes with age, thereby influencing phenotypic outcome.

Direct decarboxylation of 5-carboxylcytosine by DNA C5-methyltransferases.

S-Adenosylmethionine-dependent DNA methyltransferases (MTases) perform direct methylation of cytosine to yield 5-methylcytosine (5mC), which serves as part of the epigenetic regulation mechanism in vertebrates. Active demethylation of 5mC by TET oxygenases produces 5-formylcytosine (fC) and 5-carboxylcytosine (caC), which were shown to be enzymatically excised and then replaced with an unmodified nucleotide. Here we find that both bacterial and mammalian C5-MTases can catalyze the direct decarboxylation of caC yielding unmodified cytosine in DNA in vitro but are inert toward fC. The observed atypical enzymatic C-C bond cleavage reaction provides a plausible precedent for a direct reversal of caC to the unmodified state in DNA and offers a unique approach for sequence-specific analysis of genomic caC.

DNA unmethylome profiling by covalent capture of CpG sites.

Dynamic patterns of cytosine-5 methylation and successive hydroxylation are part of epigenetic regulation in eukaryotes, including humans, which contributes to normal phenotypic variation and disease risk. Here we present an approach for the mapping of unmodified regions of the genome, which we call the unmethylome. Our technique is based on DNA methyltransferase-directed transfer of activated groups and covalent biotin tagging of unmodified CpG sites followed by affinity enrichment and interrogation on tiling microarrays or next generation sequencing. Control experiments and pilot studies of human genomic DNA from cultured cells and tissues demonstrate that, along with providing a unique cross-section through the chemical landscape of the epigenome, the methyltransferase-directed transfer of activated groups-based approach offers high precision and robustness as compared with existing affinity-based techniques.

5-hmC in the brain is abundant in synaptic genes and shows differences at the exon-intron boundary.

The 5-methylcytosine (5-mC) derivative 5-hydroxymethylcytosine (5-hmC) is abundant in the brain for unknown reasons. Here we characterize the genomic distribution of 5-hmC and 5-mC in human and mouse tissues. We assayed 5-hmC by using glucosylation coupled with restriction-enzyme digestion and microarray analysis. We detected 5-hmC enrichment in genes with synapse-related functions in both human and mouse brain. We also identified substantial tissue-specific differential distributions of these DNA modifications at the exon-intron boundary in human and mouse. This boundary change was mainly due to 5-hmC in the brain but due to 5-mC in non-neural contexts. This pattern was replicated in multiple independent data sets and with single-molecule sequencing. Moreover, in human frontal cortex, constitutive exons contained higher levels of 5-hmC relative to alternatively spliced exons. Our study suggests a new role for 5-hmC in RNA splicing and synaptic function in the brain.

5-Hydroxymethylcytosine--the elusive epigenetic mark in mammalian DNA.

Over the past decade, epigenetic phenomena claimed a central role in cell regulatory processes and proved to be important factors for understanding complex human diseases. One of the best understood epigenetic mechanisms is DNA methylation. In the mammalian genome, cytosines (C) were long known to exist in two functional states: unmethylated or methylated at the 5-position of the pyrimidine ring (5mC). Recent studies of genomic DNA from the human and mouse brain, neurons and from mouse embryonic stem cells found that a substantial fraction of 5mC in CpG dinucleotides is converted to 5-hydroxymethyl-cytosine (hmC) by the action of 2-oxoglutarate- and Fe(ii)-dependent oxygenases of the TET family. These findings provided important clues in a long elusive mechanism of active DNA demethylation and bolstered a fresh wave of studies in the area of epigenetic regulation in mammals. This review is dedicated to critical assessment of the most popular techniques with respect to their suitability for analysis of hmC in mammalian genomes. It also discusses the most recent data on biochemical and chemical aspects of the formation and further conversion of this nucleobase in DNA and its possible biological roles in cell differentiation, embryogenesis and brain function.

Bioinformatic and partial functional analysis of pEspA and pEspB, two plasmids from Exiguobacterium arabatum sp. nov. RFL1109.

The complete nucleotide sequences of two plasmids from Exiguobacterium arabatum sp. nov. RFL1109, pEspA (4563bp) and pEspB (38,945bp), have been determined. Five ORFs were identified in the pEspA plasmid, and putative functions were assigned to two of them. Using deletion mapping approach, the Rep-independent replication region of pEspA, which functions in Bacillus subtilis, was localized within a 0.6kb DNA region. Analysis of the pEspB sequence revealed 42 ORFs. From these, function of two genes encoding enzymes of the Lsp1109I restriction-modification system was confirmed experimentally, while putative functions of another 18 ORFs were suggested based on comparative analysis. Three functional regions have been proposed for the pEspB plasmid: the putative conjugative transfer region, the region involved in plasmid replication and maintenance, and the region responsible for transposition of the IS21 family-like transposable elements.

Generation of DNA cleavage specificities of type II restriction endonucleases by reassortment of target recognition domains.

Type II restriction endonucleases (REases) cleave double-stranded DNA at specific sites within or close to their recognition sequences. Shortly after their discovery in 1970, REases have become one of the primary tools in molecular biology. However, the list of available specificities of type II REases is relatively short despite the extensive search for them in natural sources and multiple attempts to artificially change their specificity. In this study, we examined the possibility of generating cleavage specificities of REases by swapping putative target recognition domains (TRDs) between the type IIB enzymes AloI, PpiI, and TstI. Our results demonstrate that individual TRDs recognize distinct parts of the bipartite DNA targets of these enzymes and are interchangeable. Based on these properties, we engineered a functional type IIB REase having previously undescribed DNA specificity. Our study suggests that the TRD-swapping approach may be used as a general technique for the generation of type II enzymes with predetermined specificities.