Direct decarboxylation of ten-eleven translocation-produced 5-carboxylcytosine in mammalian genomes forms a new mechanism for active DNA demethylation

Direct single-nucleotide resolution sequencing of DNA 5-methylcytosine using engineered DNA methyltransferase-mediated CMD-seq

5-Methylcytosine (5mC) is a crucial epigenetic modification in the mammalian genome, primarily occurring at CG dinucleotides. Accurate localization of 5mC is essential for understanding its functional significance. In this study, we discovered a novel DNA methyltransferase, designated M.MedI, from the bacterium Mycoplasmopsis edwardii. M.MedI exhibits carboxymethylation activity towards cytosines in CG sites in DNA. We further engineered a variant of M.MedI by mutating its critical active site residue 377 asparagine (N) to lysine (K), resulting in M.MedI-N377K. This engineered M.MedI-N377K enzyme demonstrated superior carboxymethylation activity towards cytosines in CG sites in both unmethylated and hemi-methylated DNA. Utilizing the newly identified M.MedI-N377K methyltransferase, we developed a novel method, engineered DNA methyltransferase-mediated carboxymethylation deamination sequencing (CMD-seq), for the stoichiometric detection of 5mC in DNA at single-nucleotide resolution. In CMD-seq, M.MedI-N377K efficiently transfers a carboxymethyl group to cytosines in CG sites in the presence of carboxy-S-adenosyl-l-methionine (caSAM), generating 5-carboxymethylcytosine (5camC). Subsequent treatment with the deaminase A3A deaminates 5mC to form thymine (T), which pairs with adenine (A) and is read as T, while 5camC remains unchanged, pairing with guanine (G) and being read as cytosine (C) during sequencing. We successfully applied CMD-seq to quantify 5mC sites in the promoters of tumor suppressor genes RASSF1A and SHOX2 in human lung cancer tissue and adjacent normal tissue. The quantification results were highly comparable to those obtained using traditional bisulfite sequencing. Overall, CMD-seq provides a valuable tool for bisulfite-free, single-nucleotide resolution, and quantitative detection of 5mC in limited DNA samples.

Temporal epigenome modulation enables efficient bacteriophage engineering and functional analysis of phage DNA modifications

Lytic bacteriophages hold substantial promise in medical and biotechnological applications. Therefore a comprehensive understanding of phage infection mechanisms is crucial. CRISPR-Cas systems offer a way to explore these mechanisms via site-specific phage mutagenesis. However, phages can resist Cas-mediated cleavage through extensive DNA modifications like cytosine glycosylation, hindering mutagenesis efficiency. Our study utilizes the eukaryotic enzyme NgTET to temporarily reduce phage DNA modifications, facilitating Cas nuclease cleavage and enhancing mutagenesis efficiency. This approach enables precise DNA targeting and seamless point mutation integration, exemplified by deactivating specific ADP-ribosyltransferases crucial for phage infection. Furthermore, by temporally removing DNA modifications, we elucidated the effects of these modifications on T4 phage infections without necessitating gene deletions. Our results present a strategy enabling the investigation of phage epigenome functions and streamlining the engineering of phages with cytosine DNA modifications. The described temporal modulation of the phage epigenome is valuable for synthetic biology and fundamental research to comprehend phage infection mechanisms through the generation of mutants.

Simultaneous detection of 5-methylcytosine and 5-hydroxymethylcytosine at specific genomic loci by engineered deaminase-assisted sequencing

Cytosine modifications, particularly 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), play crucial roles in numerous biological processes. Current analytical methods are often constrained to the separate detection of either 5mC or 5hmC, or the combination of both modifications. The ability to simultaneously detect C, 5mC, and 5hmC at the same genomic locations with precise stoichiometry is highly desirable. Herein, we introduce a method termed engineered deaminase-assisted sequencing (EDA-seq) for the simultaneous quantification of C, 5mC, and 5hmC at the same genomic sites. EDA-seq utilizes a specially engineered protein, derived from human APOBEC3A (A3A), known as eA3A-M5. eA3A-M5 exhibits distinct deamination capabilities for C, 5mC, and 5hmC. In EDA-seq, C undergoes complete deamination and is sequenced as T. 5mC is partially deaminated resulting in a mixed readout of T and C, and 5hmC remains undeaminated and is read as C. Consequently, the proportion of T readouts (P T) reflects the collective occurrences of C and 5mC, regulated by the deamination rate of 5mC (R 5mC). By determining R 5mC and P T values, we can deduce the precise levels of C, 5mC, and 5hmC at particular genomic locations. We successfully used EDA-seq to simultaneously measure C, 5mC, and 5hmC at specific loci within human lung cancer tissue and their normal counterpart. The results from EDA-seq demonstrated a strong concordance with those obtained from the combined application of BS-seq and ACE-seq methods. EDA-seq eliminates the need for bisulfite treatment, DNA oxidation or glycosylation and uniquely enables simultaneous quantification of C, 5mC and 5hmC at the same genomic locations.

Whole-Genome Mapping of Epigenetic Modification of 5-Formylcytosine at Single-Base Resolution by Chemical Labeling Enrichment and Deamination Sequencing

DNA cytosine methylation (5-methylcytosine, 5mC) is a predominant epigenetic modification that plays a critical role in a variety of biological and pathological processes in mammals. In active DNA demethylation, the 10-11 translocation (TET) dioxygenases can sequentially oxidize 5mC to generate three modified forms of cytosine, 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). Beyond being a demethylation intermediate, recent studies have shown that 5fC has regulatory functions in gene expression and chromatin organization. While some methods have been developed to detect 5fC, genome-wide mapping of 5fC at base resolution is still highly desirable. Herein, we propose a chemical labeling enrichment and deamination sequencing (CLED-seq) method for detecting 5fC in genomic DNA at single-base resolution. The CLED-seq method utilizes selective labeling and enrichment of 5fC-containing DNA fragments, followed by deamination mediated by apolipoprotein B mRNA-editing catalytic polypeptide-like 3A (APOBEC3A or A3A) and sequencing. In the CLED-seq process, while all C, 5mC, and 5hmC are interpreted as T during sequencing, 5fC is still read as C, enabling the precise detection of 5fC in DNA. Using the proposed CLED-seq method, we accomplished genome-wide mapping of 5fC in mouse embryonic stem cells. The mapping study revealed that promoter regions enriched with 5fC overlapped with H3K4me1, H3K4me3, and H3K27ac marks. These findings suggest a correlation between 5fC marks and active gene expression in mESCs. In conclusion, CLED-seq is a straightforward, bisulfite-free method that offers a valuable tool for detecting 5fC in genomes at a single-base resolution.

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.

Whole-Genome Sequencing of 5-Hydroxymethylcytosine at Base Resolution by Bisulfite-Free Single-Step Deamination with Engineered Cytosine Deaminase

The epigenetic modification 5-hydroxymethylcytosine (5hmC) plays a crucial role in the regulation of gene expression. Although some methods have been developed to detect 5hmC, direct genome-wide mapping of 5hmC at base resolution is still highly desirable. Herein, we proposed a single-step deamination sequencing (SSD-seq) method, designed to precisely map 5hmC across the genome at single-base resolution. SSD-seq takes advantage of a screened engineered human apolipoprotein B mRNA-editing catalytic polypeptide-like 3A (A3A) protein, known as eA3A-v10, to selectively deaminate cytosine (C) and 5-methylcytosine (5mC) but not 5hmC. During sequencing, the deaminated C and 5mC are converted to uracil (U) and thymine (T), read as T in the sequencing data. However, 5hmC remains unaffected by eA3A-v10 and is read as C during sequencing. Consequently, the presence of C in the sequence reads indicates the original 5hmC. We applied SSD-seq to generate a base-resolution map of 5hmC in human lung tissue. Our findings revealed that 5hmC was predominantly localized to CpG dinucleotides. Furthermore, the base-resolution map of 5hmC generated by SSD-seq demonstrated a strong correlation with prior ACE-seq results. The advantages of SSD-seq are its single-step process, absence of bisulfite treatment or DNA glycosylation, cost effectiveness, and ability to detect and quantify 5hmC directly at single-base resolution.

Single-Nucleotide Resolution Mapping of N6-Methyladenine in Genomic DNA

N6-Methyladenine (6mA) is a naturally occurring DNA modification in both prokaryotes and eukaryotes. Herein, we developed a deaminase-mediated sequencing (DM-seq) method for genome-wide mapping of 6mA at single-nucleotide resolution. The method capitalizes on the selective deamination of adenine, but not 6mA, in DNA mediated by an evolved adenine deaminase, ABE8e. By employing this method, we achieved genome-wide mapping of 6mA in Escherichia coli and in mammalian mitochondrial DNA (mtDNA) at single-nucleotide resolution. We found that the 6mA sites are mainly located in the GATC motif in the E. coli genome. We also identified 17 6mA sites in mtDNA of HepG2 cells, where all of the 6mA sites are distributed in the heavy strand of mtDNA. We envision that DM-seq will be a valuable tool for uncovering new functions of 6mA in DNA and for exploring its potential roles in mitochondria-related human diseases.

Substrate DNA length regulates the activity of TET 5-methylcytosine dioxygenases

The ten-eleven translocation (TET) isoforms (TET1-3) play critical roles in epigenetic transcription regulation. In addition, mutations in the TET2 gene are frequently detected in patients with glioma and myeloid malignancies. TET isoforms can oxidize 5-methylcytosine to 5-hydroxymethylcytosine, 5-formylcytosine, and 5-carboxylcytosine, by iterative oxidation. The in vivo DNA demethylation activity of TET isoforms may depend on many factors including enzyme's structural features, its interaction with DNA-binding proteins, chromatin context, DNA sequence, DNA length, and configuration. The rationale for this study is to identify the preferred DNA length and configuration in the substrates of TET isoforms. We have used a highly sensitive LC-MS/MS-based method to compare the substrate preference of TET isoforms. To this end, four DNA substrate sets (S1, S2, S3, S4) of different sequences were chosen. In addition, in each set, four different lengths of DNA substrates comprising 7-, 13-, 19-, and 25-mer nucleotides were synthesized. Each DNA substrate was further used in three different configurations, that is, double stranded symmetrically-methylated, double stranded hemi-methylated, and single stranded single-methylated to evaluate their effect on TET-mediated 5mC oxidation. We demonstrate that mouse TET1 (mTET1) and human TET2 (hTET2) have highest preference for 13-mer dsDNA substrates. Increasing or decreasing the length of dsDNA substrate reduces product formation. In contrast to their dsDNA counterparts, the length of ssDNA substrates did not have a predictable effect on 5mC oxidation. Finally, we show that substrate specificity of TET isoforms correlates with their DNA binding efficiency. Our results demonstrate that mTET1 and hTET2 prefer 13-mer dsDNA as a substrate over ssDNA. These results may help elucidate novel properties of TET-mediated 5mC oxidation and help develop novel diagnostic tools to detect TET2 function in patients.

Dynamic changes of DNA methylation induced by benzo(a)pyrene in cancer

Benzo(a)pyrene (BaP), the earliest and most significant carcinogen among polycyclic aromatic hydrocarbons (PAHs), has been found in foods, tobacco smoke, and automobiles exhaust, etc. Exposure to BaP induced DNA damage directly, or oxidative stress-related damage, resulting in cell apoptosis and carcinogenesis in human respiratory system, digestive system, reproductive system, etc. Moreover, BaP triggered genome-wide epigenetic alterations by methylation, which might cause disturbances in regulation of gene expression, and thereby induced cancer. It has been proved that BaP reduced genome-wide DNA methylation, and activated proto-oncogene by hypomethylation in the promoter region, but silenced tumor suppressor genes by promoter hypermethylation, resulting in cancer initiation and progression. Here we summarized the changes in DNA methylation in BaP exposure, and revealed the methylation of DNA plays a role in cancer development.

DNA Glycosylases Define the Outcome of Endogenous Base Modifications

Chemically modified nucleic acid bases are sources of genomic instability and mutations but may also regulate gene expression as epigenetic or epitranscriptomic modifications. Depending on the cellular context, they can have vastly diverse impacts on cells, from mutagenesis or cytotoxicity to changing cell fate by regulating chromatin organisation and gene expression. Identical chemical modifications exerting different functions pose a challenge for the cell's DNA repair machinery, as it needs to accurately distinguish between epigenetic marks and DNA damage to ensure proper repair and maintenance of (epi)genomic integrity. The specificity and selectivity of the recognition of these modified bases relies on DNA glycosylases, which acts as DNA damage, or more correctly, as modified bases sensors for the base excision repair (BER) pathway. Here, we will illustrate this duality by summarizing the role of uracil-DNA glycosylases, with particular attention to SMUG1, in the regulation of the epigenetic landscape as active regulators of gene expression and chromatin remodelling. We will also describe how epigenetic marks, with a special focus on 5-hydroxymethyluracil, can affect the damage susceptibility of nucleic acids and conversely how DNA damage can induce changes in the epigenetic landscape by altering the pattern of DNA methylation and chromatin structure.

Engineered APOBEC3C Sequencing Enables Bisulfite-Free and Direct Detection of DNA Methylation at a Single-Base Resolution

DNA methylation (5-methylcytosine, 5mC) is the most important epigenetic modification in mammals. Deciphering the roles of 5mC relies on the quantitative detection of 5mC at the single-base resolution. Bisulfite sequencing (BS-seq) is the most often employed technique for mapping 5mC in DNA. However, bisulfite treatment may cause serious degradation of input DNA due to the harsh reaction conditions. Here, we engineered the human apolipoprotein B mRNA-editing catalytic polypeptide-like 3C (A3C) protein to endow the engineered A3C (eA3C) protein with differential deamination activity toward cytosine and 5mC. By the virtue of the unique property of eA3C, we proposed an engineered A3C sequencing (EAC-seq) method for the bisulfite-free and quantitative mapping of 5mC in DNA at the single-base resolution. In EAC-seq, the eA3C protein can deaminate C but not 5mC, which is employed to differentiate C and 5mC in sequencing. Using the EAC-seq method, we quantitatively detected 5mC in genomic DNA of lung cancer tissue. In contrast to the harsh reaction conditions of BS-seq, which could lead to significant degradation of DNA, the whole procedure of EAC-seq is carried out under mild conditions, thereby preventing DNA damage. Taken together, the EAC-seq approach is bisulfite-free and straightforward, making it an invaluable tool for the quantitative detection of 5mC in limited DNA at the single-base resolution.

Active DNA demethylation promotes cell fate specification and the DNA damage response

Neurons harbor high levels of single-strand DNA breaks (SSBs) that are targeted to neuronal enhancers, but the source of this endogenous damage remains unclear. Using two systems of postmitotic lineage specification-induced pluripotent stem cell-derived neurons and transdifferentiated macrophages-we show that thymidine DNA glycosylase (TDG)-driven excision of methylcytosines oxidized with ten-eleven translocation enzymes (TET) is a source of SSBs. Although macrophage differentiation favors short-patch base excision repair to fill in single-nucleotide gaps, neurons also frequently use the long-patch subpathway. Disrupting this gap-filling process using anti-neoplastic cytosine analogs triggers a DNA damage response and neuronal cell death, which is dependent on TDG. Thus, TET-mediated active DNA demethylation promotes endogenous DNA damage, a process that normally safeguards cell identity but can also provoke neurotoxicity after anticancer treatments.

Formation and removal of 1,N6-dimethyladenosine in mammalian transfer RNA

RNA molecules harbor diverse modifications that play important regulatory roles in a variety of biological processes. Over 150 modifications have been identified in RNA molecules. N⁶-methyladenosine (m⁶A) and 1-methyladenosine (m¹A) are prevalent modifications occurring in various RNA species of mammals. Apart from the single methylation of adenosine (m⁶A and m¹A), dual methylation modification occurring in the nucleobase of adenosine, such as N⁶,N⁶-dimethyladenosine (m^6,6A), also has been reported to be present in RNA of mammals. Whether there are other forms of dual methylation modification occurring in the nucleobase of adenosine other than m^6,6A remains elusive. Here, we reported the existence of a novel adenosine dual methylation modification, i.e. 1,N⁶-dimethyladenosine (m^1,6A), in tRNAs of living organisms. We confirmed that m^1,6A is located at position 58 of tRNAs and is prevalent in mammalian cells and tissues. The measured level of m^1,6A ranged from 0.0049% to 0.047% in tRNAs. Furthermore, we demonstrated that TRMT6/61A could catalyze the formation of m^1,6A in tRNAs and m^1,6A could be demethylated by ALKBH3. Collectively, the discovery of m^1,6A expands the diversity of RNA modifications and may elicit a new tRNA modification-mediated gene regulation pathway.

Single-Base Resolution Detection of Adenosine-to-Inosine RNA Editing by Endonuclease-Mediated Sequencing

RNA molecules contain diverse modifications that play crucial roles in a wide variety of biological processes. Adenosine-to-inosine (A-to-Ino) RNA editing is one of the most prevalent modifications among all types of RNA. Abnormal A-to-InoRNA editing has been demonstrated to be associated with many human diseases. Identification of A-to-Ino editing sites is indispensable to deciphering their biological roles. Herein, by employing the unique property of human endonuclease V (hEndoV), we proposed a hEndoV-mediated sequencing (hEndoV-seq) method for the single-base resolution detection of A-to-InoRNA editing sites. In this approach, the terminal 3'OH of RNA is first blocked by 3'-deoxyadenosine (3'-deoxy-A). Specific cleavage of Ino sites by hEndoV protein produces new terminal 3'OH, which can be identified by sequencing analysis, and therefore offers the site-specific detection of Ino in RNA. The principle of hEndoV-seq is straightforward and the analytical procedure is simple. No chemical reaction is involved in the sequencing library preparation. The whole procedure in hEndoV-seq is carried out under mild conditions and RNA is not prone to degradation. Taken together, the proposed hEndoV-seq method is capable of site-specific identification of A-to-Ino editing in RNA, which provides a valuable tool for elucidating the functions of A-to-Ino editing in RNA.

The Methylation Game: Epigenetic and Epitranscriptomic Dynamics of 5-Methylcytosine

DNA and RNA methylation dynamics have been linked to a variety of cellular processes such as development, differentiation, and the maintenance of genome integrity. The correct deposition and removal of methylated cytosine and its oxidized analogues is pivotal for cellular homeostasis, rapid responses to exogenous stimuli, and regulated gene expression. Uncoordinated expression of DNA/RNA methyltransferases and demethylase enzymes has been linked to genome instability and consequently to cancer progression. Furthermore, accumulating evidence indicates that post-transcriptional DNA/RNA modifications are important features in DNA/RNA function, regulating the timely recruitment of modification-specific reader proteins. Understanding the biological processes that lead to tumorigenesis or somatic reprogramming has attracted a lot of attention from the scientific community. This work has revealed extensive crosstalk between epigenetic and epitranscriptomic pathways, adding a new layer of complexity to our understanding of cellular programming and responses to environmental cues. One of the key modifications, m5C, has been identified as a contributor to regulation of the DNA damage response (DDR). However, the various mechanisms of dynamic m5C deposition and removal, and the role m5C plays within the cell, remains to be fully understood.

Identification of Inosine and 2'-O-Methylinosine Modifications in Yeast Messenger RNA by Liquid Chromatography-Tandem Mass Spectrometry Analysis

The discovery of reversible modifications in messenger RNA (mRNA) opens new research directions in RNA modification-mediated epigenetic regulation. Yeast is an extensively used model organism in molecular biology. Systematic investigation and profiling of modifications in yeast mRNA would promote our understanding of the physiological regulation mechanisms in yeast. However, due to the high abundance of ribosomal RNA (rRNA) and transfer RNA (tRNA) in total RNA, isolation of low abundance of mRNA frequently suffers from the contamination of rRNA and tRNA, which will lead to the false-positive determination and inaccurate quantification of modifications in mRNA. Therefore, obtaining high-purity mRNA is critical for precise determination and accurate quantification of modifications in mRNA, especially for studies that focus on discovering new ones. Herein, we proposed a successive orthogonal isolation method by combining polyT-based purification and agarose gel electrophoresis purification for extracting high-purity mRNA. With the extracted high-purity yeast mRNA, we systemically explored the modifications in yeast mRNA by liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) analysis. The results showed that in addition to the previously reported eight kinds of modifications, two novel modifications of inosine (Ino) and 2'-O-methylinosine (Im) were identified to be prevalent in yeast mRNA. It is worth noting that Im was reported for the first time, to the best of our knowledge, to exist in living organisms in the three domains of life. Moreover, we observed that the levels of 10 kinds of modifications including Ino and Im in yeast mRNA exhibited dynamic change at different growth stages of yeast cells. Furthermore, Im in mRNA showed a significant decrease while in response to H₂O₂ treatment. These results indicated that the two newly identified modifications in yeast mRNA were involved in yeast cell growth and response to environmental stress. Taken together, we reported two new modifications of Ino and Im in yeast mRNA, which expends the diversity of RNA modifications in yeast and also suggests new regulators for modulating yeast physiological functions.

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.

Bisulfite-free and single-nucleotide resolution sequencing of DNA epigenetic modification of 5-hydroxymethylcytosine by engineered deaminase

The discovery of 5-hydroxymethylcytosine (5hmC) in mammalian genomes is a landmark in epigenomics study. Similar to 5-methylcytosine (5mC), 5hmC is viewed a critical epigenetic modification. Deciphering the functions of 5hmC...

Genome-wide mapping of N4-methylcytosine at single-base resolution by APOBEC3A-mediated deamination sequencing

N⁴-methylcytosine (4mC) is a natural DNA modification occurring in thermophiles and plays important roles in restriction-modification (R-M) systems in bacterial genomes. However, the precise location and sequence context of 4mC in the whole genome are limited. In this study, we developed an APOBEC3A-mediated deamination sequencing (4mC-AMD-seq) method for genome-wide mapping of 4mC at single-base resolution. In the 4mC-AMD-seq method, cytosine and 5-methylcytosine (5mC) are deaminated by APOBEC3A (A3A) protein to generate uracil and thymine, both of which are read as thymine in sequencing, while 4mC is resistant to deamination and therefore read as cytosine. Thus, the readouts of cytosines from sequencing could manifest the original 4mC sites in genomes. With the 4mC-AMD-seq method, we achieved the genome-wide mapping of 4mC in Deinococcus radiodurans (D. radiodurans). In addition, we confirmed that 4mC, but not 5mC, was the major modification in the D. radiodurans genome. We identified 1586 4mC sites in the genome of D. radiodurans, among which 564 sites were located in the CCGCGG motif. The average methylation levels in the CCGCGG motif and non-CCGCGG sequence were 70.0% and 22.8%, respectively. We envision that the 4mC-AMD-seq method will facilitate the investigation of 4mC functions, including the 4mC-involved R-M systems, in uncharacterized but potentially useful strains.

Genome-wide mapping of N⁴-methylcytosine (4mC) at single-base resolution with APOBEC3A-mediated deamination sequencing (4mC-AMD-seq).

An enzyme-mediated bioorthogonal labeling method for genome-wide mapping of 5-hydroxymethyluracil

DNA 5-hydroxymethyluracil (5hmU) is a thymine modification existing in the genomes of various organisms. The post-replicative formation of 5hmU occurs via hydroxylation of thymine by ten-eleven translocation (TET) dioxygenases in mammals and J-binding proteins (JBPs) in protozoans, respectively. In addition, 5hmU can also be generated through oxidation of thymine by reactive oxygen species or deamination of 5hmC by cytidine deaminase. While the biological roles of 5hmU have not yet been fully explored, determining its genomic location will highly assist in elucidating its functions. Herein, we report a novel enzyme-mediated bioorthogonal labeling method for selective enrichment of 5hmU in genomes. 5hmU DNA kinase (5hmUDK) was utilized to selectively install an azide (N₃) group or alkynyl group into the hydroxyl moiety of 5hmU followed by incorporation of the biotin linker through click chemistry, which enabled the capture of 5hmU-containing DNA fragments via streptavidin pull-down. The enriched fragments were applied to deep sequencing to determine the genomic distribution of 5hmU. With this established enzyme-mediated bioorthogonal labeling strategy, we achieved the genome-wide mapping of 5hmU in Trypanosoma brucei. The method described here will allow for a better understanding of the functional roles and dynamics of 5hmU in genomes.

Intragenomic Decarboxylation of 5-Carboxy-2'-deoxycytidine

Cellular DNA is composed of four canonical nucleosides (dA, dC, dG and T), which form two Watson-Crick base pairs. In addition, 5-methylcytosine (mdC) may be present. The methylation of dC to mdC is known to regulate transcriptional activity. Next to these five nucleosides, the genome, particularly of stem cells, contains three additional dC derivatives, which are formed by stepwise oxidation of the methyl group of mdC with the help of Tet enzymes. These are 5-hydroxymethyl-dC (hmdC), 5-formyl-dC (fdC), and 5-carboxy-dC (cadC). It is believed that fdC and cadC are converted back into dC, which establishes an epigenetic control cycle that starts with methylation of dC to mdC, followed by oxidation and removal of fdC and cadC. While fdC was shown to undergo intragenomic deformylation to give dC directly, a similar decarboxylation of cadC was postulated but not yet observed on the genomic level. By using metabolic labelling, we show here that cadC decarboxylates in several cell types, which confirms that both fdC and cadC are nucleosides that are directly converted back to dC within the genome by C-C bond cleavage.