Natural, modified DNA bases

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.

Harnessing non-Watson-Crick's base pairing to enhance CRISPR effectors cleavage activities and enable gene editing in mammalian cells

Genomic DNA of the cyanophage S-2L virus is composed of 2-aminoadenine (Z), thymine (T), guanine (G), and cytosine (C), forming the genetic alphabet ZTGC, which violates Watson-Crick base pairing rules. The Z-base has an extra amino group on the two position that allows the formation of a third hydrogen bond with thymine in DNA strands. Here, we explored and expanded applications of this non-Watson-Crick base pairing in protein expression and gene editing. Both ZTGC-DNA (Z-DNA) and ZUGC-RNA (Z-RNA) produced in vitro show detectable compatibility and can be decoded in mammalian cells, including Homo sapiens cells. Z-crRNA can guide CRISPR-effectors SpCas9 and LbCas12a to cleave specific DNA through non-Watson-Crick base pairing and boost cleavage activities compared to A-crRNA. Z-crRNA can also allow for efficient gene and base editing in human cells. Together, our results help pave the way for potential strategies for optimizing DNA or RNA payloads for gene editing therapeutics and give insights to understanding the natural Z-DNA genome.

Universal Design of Luciferase Fusion Proteins for Epigenetic Modifications Detection Based on Bioluminescence Resonance Energy Transfer

Global hypomethylation and promoter hypermethylation of tumor-suppressor genes are the hallmarks of cancer. We previously reported a global DNA methylation level sensing system based on dual-color bioluminescence resonance energy transfer (BRET) using methyl-CpG binding domain (MBD)-fused firefly luciferase (Fluc) and unmethyl-CpG binding domain (CXXC)-fused Oplophorus luciferase (Oluc). Moreover, BRET-based hydroxymethylation and hemi-methylation level sensing systems have been developed using hydroxymethyl-CpG and hemi-methyl-CpG binding domain-fused Fluc. These studies suggest that target epigenetic modifications can be simultaneously quantified using target-modification-binding protein-fused luciferases. In this study, we focused on the SnoopTag (SnT)/SnoopCatcher (SnC) protein ligation system to establish a universal design for fusion protein construction for any combination. SnT spontaneously forms an isopeptide bond with SnC; therefore, any kind of fusion protein would be constructed by the SnT/SnC system. To establish the proof of concept, MBD-SnT, CXXC-SnT, and SnC-Oluc were prepared and ligated MBD-SnT or CXXC-SnT to SnC-Oluc. The ligation products of MBD-SnT-SnC-Oluc and CXXC-SnT-SnC-Oluc showed luciferase activity and specific binding activity to methyl-CpG and unmethyl-CpG, respectively. The BRET signal using MBD-SnT-SnC-Oluc and CXXC-SnT-SnC-Oluc increased the amount of methyl-CpG and unmethyl-CpG in genomic DNA, respectively. There was a significant negative correlation between the BRET signals; therefore, the global DNA methylation level was quantified using the BRET signals (R² = 0.99, and R.S.D. <3.5%). These results indicate that the SnT/SnC protein ligation system can be utilized to construct target modification-binding protein-fused luciferases in any combination that detects target modifications in genomic DNA based on BRET.

MuLan-Methyl-multiple transformer-based language models for accurate DNA methylation prediction

Transformer-based language models are successfully used to address massive text-related tasks. DNA methylation is an important epigenetic mechanism, and its analysis provides valuable insights into gene regulation and biomarker identification. Several deep learning-based methods have been proposed to identify DNA methylation, and each seeks to strike a balance between computational effort and accuracy. Here, we introduce MuLan-Methyl, a deep learning framework for predicting DNA methylation sites, which is based on 5 popular transformer-based language models. The framework identifies methylation sites for 3 different types of DNA methylation: N6-adenine, N4-cytosine, and 5-hydroxymethylcytosine. Each of the employed language models is adapted to the task using the "pretrain and fine-tune" paradigm. Pretraining is performed on a custom corpus of DNA fragments and taxonomy lineages using self-supervised learning. Fine-tuning aims at predicting the DNA methylation status of each type. The 5 models are used to collectively predict the DNA methylation status. We report excellent performance of MuLan-Methyl on a benchmark dataset. Moreover, we argue that the model captures characteristic differences between different species that are relevant for methylation. This work demonstrates that language models can be successfully adapted to applications in biological sequence analysis and that joint utilization of different language models improves model performance. Mulan-Methyl is open source, and we provide a web server that implements the approach.

Epigenetic Pyrimidine Nucleotides in Competition with Natural dNTPs as Substrates for Diverse DNA Polymerases

Five 2'-deoxyribonucleoside triphosphates (dNTPs) derived from epigenetic pyrimidines (5-methylcytosine, 5-hydroxymethylcytosine, 5-formylcytosine, 5-hydroxymethyluracil, and 5-formyluracil) were prepared and systematically studied as substrates for nine DNA polymerases in competition with natural dNTPs by primer extension experiments. The incorporation of these substrates was evaluated by a restriction endonucleases cleavage-based assay and by a kinetic study of single nucleotide extension. All of the modified pyrimidine dNTPs were good substrates for the studied DNA polymerases that incorporated a significant percentage of the modified nucleotides into DNA even in the presence of natural nucleotides. 5-Methylcytosine dNTP was an even better substrate for most polymerases than natural dCTP. On the other hand, 5-hydroxymethyl-2'-deoxyuridine triphosphate was not the best substrate for SPO1 DNA polymerase, which naturally synthesizes 5hmU-rich genomes of the SPO1 bacteriophage. The results shed light onto the possibility of gene silencing through recycling and random incorporation of epigenetic nucleotides and into the replication of modified bacteriophage genomes.

Epigenetic Marks, DNA Damage Markers, or Both? The Impact of Desiccation and Accelerated Aging on Nucleobase Modifications in Plant Genomic DNA

Modifications of DNA nucleobases are present in all forms of life. The purpose of these modifications in eukaryotic cells, however, is not always clear. Although the role of 5-methylcytosine (m⁵C) in epigenetic regulation and the maintenance of stability in plant genomes is becoming better understood, knowledge pertaining to the origin and function of oxidized nucleobases is still scarce. The formation of 5-hydroxymetylcytosine (hm⁵C) in plant genomes is especially debatable. DNA modifications, functioning as regulatory factors or serving as DNA injury markers, may have an effect on DNA structure and the interaction of genomic DNA with proteins. Thus, these modifications can influence plant development and adaptation to environmental stress. Here, for the first time, the changes in DNA global levels of m⁵C, hm⁵C, and 8-oxo-7,8-dihydroguanine (8-oxoG) measured by ELISA have been documented in recalcitrant embryonic axes subjected to desiccation and accelerated aging. We demonstrated that tissue desiccation induces a similar trend in changes in the global level of hm⁵C and 8-oxoG, which may suggest that they both originate from the activity of reactive oxygen species (ROS). Our study supports the premise that m⁵C can serve as a marker of plant tissue viability whereas oxidized nucleobases, although indicating a cellular redox state, cannot.

1H NMR Chemical Exchange Techniques Reveal Local and Global Effects of Oxidized Cytosine Derivatives

5-Carboxycytosine (5caC) is a rare epigenetic modification found in nucleic acids of all domains of life. Despite its sparse genomic abundance, 5caC is presumed to play essential regulatory roles in transcription, maintenance and base-excision processes in DNA. In this work, we utilize nuclear magnetic resonance (NMR) spectroscopy to address the effects of 5caC incorporation into canonical DNA strands at multiple pH and temperature conditions. Our results demonstrate that 5caC has a pH-dependent global destabilizing and a base-pair mobility enhancing local impact on dsDNA, albeit without any detectable influence on the ground-state B-DNA structure. Measurement of hybridization thermodynamics and kinetics of 5caC-bearing DNA duplexes highlighted how acidic environment (pH 5.8 and 4.7) destabilizes the double-stranded structure by ∼10–20 kJ mol^–1 at 37 °C when compared to the same sample at neutral pH. Protonation of 5caC results in a lower activation energy for the dissociation process and a higher barrier for annealing. Studies on conformational exchange on the microsecond time scale regime revealed a sharply localized base-pair motion involving exclusively the modified site and its immediate surroundings. By direct comparison with canonical and 5-formylcytosine (5fC)-edited strands, we were able to address the impact of the two most oxidized naturally occurring cytosine derivatives in the genome. These insights on 5caC’s subtle sensitivity to acidic pH contribute to the long-standing questions of its capacity as a substrate in base excision repair processes and its purpose as an independent, stable epigenetic mark.

Evolved DNA Duplex Readers for Strand-Asymmetrically Modified 5-Hydroxymethylcytosine/5-Methylcytosine CpG Dyads

5-Methylcytosine (mC) and 5-hydroxymethylcytosine (hmC), the two main epigenetic modifications of mammalian DNA, exist in symmetric and asymmetric combinations in the two strands of CpG dyads. However, revealing such combinations in single DNA duplexes is a significant challenge. Here, we evolve methyl-CpG-binding domains (MBDs) derived from MeCP2 by bacterial cell surface display, resulting in the first affinity probes for hmC/mC CpGs. One mutant has low nanomolar affinity for a single hmC/mC CpG, discriminates against all 14 other modified CpG dyads, and rivals the selectivity of wild-type MeCP2. Structural studies indicate that this protein has a conserved scaffold and recognizes hmC and mC with two dedicated sets of residues. The mutant allows us to selectively address and enrich hmC/mC-containing DNA fragments from genomic DNA backgrounds. We anticipate that this novel probe will be a versatile tool to unravel the function of hmC/mC marks in diverse aspects of chromatin biology.

Homologues of epigenetic pyrimidines: 5-alkyl-, 5-hydroxyalkyl and 5-acyluracil and -cytosine nucleotides: synthesis, enzymatic incorporation into DNA and effect on transcription with bacterial RNA polymerase

Homologues of natural epigenetic pyrimidine nucleosides and nucleotides were designed and synthesized. They included 5-ethyl-, 5-propyl-, 5-(1-hydroxyethyl)-, 5-(1-hydroxypropyl)- and 5-acetyl- and 5-propionylcytosine and -uracil 2′-deoxyribonucleosides and their corresponding 5′-O-triphosphates (dN^XTPs). The epimers of 5-(1-hydroxyethyl)- and 5-(1-hydroxypropyl)pyrimidine nucleosides were separated and their absolute configuration was determined by a combination of X-ray and NMR analysis. The modified dN^XTPs were used as substrates for PCR synthesis of modified DNA templates used for the study of transcription with bacterial RNA polymerase. Fundamental differences in transcription efficiency were observed, depending on the various modifications. The most notable effects included pronounced stimulation of transcription from 5-ethyluracil-bearing templates (200% transcription yield compared to natural thymine) and an enhancing effect of 5-acetylcytosine versus inhibiting effect of 5-acetyluracil. In summary, these results reveal that RNA polymerase copes with dramatically altered DNA structure and suggest that these nucleobases could potentially play roles as artificial epigenetic DNA nucleobases.

Nucleotides derived from homologues of epigenetic pyrimidine bases were prepared and used for polymerase synthesis of modified DNA templates. Interesting effects of the substituents on PCR and transcription have been observed.

Epigenetic Quantification of DNA 5-Hydroxymethylcytosine Using DNA Hybridization-Based Single-Molecule Immunofluorescent Imaging

5-Hydroxymethylcytosine (5hmC) is a deoxyribonucleic acid (DNA) epigenetic modification that has an important function in embryonic development and human diseases. However, the numerous methods that have been developed to detect and quantify 5hmC, require large amounts of DNA sample to be modified via chemical reactions, which considerably limits their application with cell-free DNA (cfDNA). Meanwhile, other antibody-based methods of detecting 5hmC do not offer information about the DNA sequence. Here, in this article DNA hybridization-based single-molecule immunofluorescent imaging is presented, an ultrasensitive method of detecting 5hmC modification in DNA. Via using the probe DNA to capture the DNA fragment of interest and the 5hmC antibody to detect the 5hmC modification in DNA, the fluorescent response signal of the 5hmC modification from the secondary antibody at the single-molecule level is successfully detected. Using the method, one could determine the quantity of 5hmC in the gene of interest within 6 h. In addition, it requires only 3 pg of the DNA sample and minimal experience and training for operation and analysis.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Quantification and mapping of DNA modifications

Apart from the four canonical nucleobases, DNA molecules carry a number of natural modifications. Substantial evidence shows that DNA modifications can regulate diverse biological processes. Dynamic and reversible modifications of DNA are critical for cell differentiation and development. Dysregulation of DNA modifications is closely related to many human diseases. The research of DNA modifications is a rapidly expanding area and has been significantly stimulated by the innovations of analytical methods. With the recent advances in methods and techniques, a series of new DNA modifications have been discovered in the genomes of prokaryotes and eukaryotes. Deciphering the biological roles of DNA modifications depends on the sensitive detection, accurate quantification, and genome-wide mapping of modifications in genomic DNA. This review provides an overview of the recent advances in analytical methods and techniques for both the quantification and genome-wide mapping of natural DNA modifications. We discuss the principles, advantages, and limitations of these developed methods. It is anticipated that new methods and techniques will resolve the current challenges in this burgeoning research field and expedite the elucidation of the functions of DNA modifications.

Programmable tools for targeted analysis of epigenetic DNA modifications

Modifications of the cytosine 5-position are dynamic epigenetic marks of mammalian DNA with important regulatory roles in development and disease. Unraveling biological functions of such modified nucleobases is tightly connected with the potential of available methods for their analysis. Whereas genome-wide nucleobase quantification and mapping are first-line analyses, targeted analyses move into focus the more genomic sites with high biological significance are identified. We here review recent developments in an emerging field that addresses such targeted analyses via probes that combine a programmable, sequence-specific DNA-binding domain with the ability to directly recognize or cross-link an epigenetically modified nucleobase of interest. We highlight how such probes offer simple, high-resolution nucleobase analyses in vitro and enable in situ correlations between a nucleobase and other chromatin regulatory elements at user-defined loci on the single-cell level by imaging.

Active turnover of genomic methylcytosine in pluripotent cells

Epigenetic plasticity underpins cell potency, but the extent to which active turnover of DNA methylation contributes to such plasticity is not known, and the underlying pathways are poorly understood. Here we use metabolic labeling with stable isotopes and mass spectrometry to quantitatively address the global turnover of genomic 5-methyl-2'-deoxycytidine (mdC), 5-hydroxymethyl-2'-deoxycytidine (hmdC) and 5-formyl-2'-deoxycytidine (fdC) across mouse pluripotent cell states. High rates of mdC/hmdC oxidation and fdC turnover characterize a formative-like pluripotent state. In primed pluripotent cells, the global mdC turnover rate is about 3-6% faster than can be explained by passive dilution through DNA synthesis. While this active component is largely dependent on ten-eleven translocation (Tet)-mediated mdC oxidation, we unveil additional oxidation-independent mdC turnover, possibly through DNA repair. This process accelerates upon acquisition of primed pluripotency and returns to low levels in lineage-committed cells. Thus, in pluripotent cells, active mdC turnover involves both mdC oxidation-dependent and oxidation-independent processes.

Photocaged 5-(Hydroxymethyl)pyrimidine Nucleoside Phosphoramidites for Specific Photoactivatable Epigenetic Labeling of DNA

5-Hydroxymethylcytosine and uracil are epigenetic nucleobases, but their biological roles are still unclear. We present the synthesis of 2-nitrobenzyl photocaged 5-hydroxymethyl-2'-deoxycytidine and uridine 3'-O-phosphoramidites and their use in automated solid-phase synthesis of oligonucleotides (ONs) modified at specific positions. The ONs were used as primers for PCR to construct DNA templates modified in the promoter region that allowed switching of transcription through photochemical uncaging.

Discovery of an Unnatural DNA Modification Derived from a Natural Secondary Metabolite

Despite widespread interest for understanding how modified bases have evolved their contemporary functions, limited experimental evidence exists for measuring how close an organism is to accidentally creating a new, modified base within the framework of its existing genome. Here, we describe the biochemical and structural basis for how a single-point mutation in E. coli's naturally occurring cytosine methyltransferase can surprisingly endow a neomorphic ability to create the unnatural DNA base, 5-carboxymethylcytosine (5cxmC), in vivo. Mass spectrometry, bacterial genetics, and structure-guided biochemistry reveal this base to be exclusively derived from the natural but sparse secondary metabolite carboxy-S-adenosyl-L-methionine (CxSAM). Our discovery of a new, unnatural DNA modification reveals insights into the substrate selectivity of DNA methyltransferase enzymes, offers a promising new biotechnological tool for the characterization of the mammalian epigenome, and provides an unexpected model for how neomorphic bases could arise in nature from repurposed host metabolites.

Impact of 5-formylcytosine on the melting kinetics of DNA by 1H NMR chemical exchange

5-Formylcytosine (5fC) is a chemically edited, naturally occurring nucleobase which appears in the context of modified DNA strands. The understanding of the impact of 5fC on dsDNA physical properties is to date limited. In this work, we applied temperature-dependent 1H Chemical Exchange Saturation Transfer (CEST) NMR experiments to non-invasively and site-specifically measure the thermodynamic and kinetic influence of formylated cytosine nucleobase on the melting process involving dsDNA. Incorporation of 5fC within symmetrically positioned CpG sites destabilizes the whole dsDNA structure-as witnessed from the ∼2°C decrease in the melting temperature and 5-10 kJ mol-1 decrease in ΔG°-and affects the kinetic rates of association and dissociation. We observed an up to ∼5-fold enhancement of the dsDNA dissociation and an up to ∼3-fold reduction in ssDNA association rate constants, over multiple temperatures and for several proton reporters. Eyring and van't Hoff analysis proved that the destabilization is not localized, instead all base-pairs are affected and the transition states resembles the single-stranded conformation. These results advance our knowledge about the role of 5fC as a semi-permanent epigenetic modification and assist in the understanding of its interactions with reader proteins.

Synthetic Life with Alternative Nucleic Acids as Genetic Materials

DNA, the fundamental genetic polymer of all living organisms on Earth, can be chemically modified to embrace novel functions that do not exist in nature. The key chemical and structural parameters for genetic information storage, heredity, and evolution have been elucidated, and many xenobiotic nucleic acids (XNAs) with non-canonical structures are developed as alternative genetic materials in vitro. However, it is still particularly challenging to replace DNAs with XNAs in living cells. This review outlines some recent studies in which the storage and propagation of genetic information are achieved in vivo by expanding genetic systems with XNAs.