Site-specific quantification of 5-carboxylcytosine in DNA by chemical conversion coupled with ligation-based PCR

A disposable electrochemical magnetic immunosensor for the rapid and sensitive detection of 5-formylcytosine and 5-carboxylcytosine in DNA

5-formylcytosine (5 fC) and 5-carboxylcytosine (5caC) serve as key intermediates in DNA demethylation process with significant implications for gene regulation and disease progression. In this study, we introduce a novel electrochemical sensing platform specifically designed for the sensitive and selective detection of 5 fC and 5caC in DNA. Protein A-modified magnetic beads (ProtA-MBs) coupled with specific antibodies facilitate the immunorecognition and enrichment of these modified bases. Signal amplification is achieved through several chemical reactions involving the interaction between N3-kethonaxl and guanine, copper-free click chemistry for the attachment of dibenzocyclooctyne (DBCO)-Biotin, and the subsequent recognition by streptavidin-conjugated horseradish peroxidase (SA-HRP). The assay's readout is performed on a disposable laser-induced graphene (LIG) electrode, modified with the bead-antibody-DNA complex in a magnetic field, and analyzed using differential pulse voltammetry in a system employing hydroquinone (HQ) as the redox mediator and H2O2 as the substrate. This immunosensor displayed excellent sensitivity, with detection limits of 14.8 fM for 5 fC across a 0.1-1000 pM linear range and 87.4 fM for 5caC across a 0.5-5000 pM linear range, and maintained high selectivity even in the presence of interferences from other DNA modifications. Successful application in quantifying 5 fC and 5caC in genomic DNA from cell extracts, with recovery rates between 97.7% to 102.9%, underscores its potential for clinical diagnostics. N3-kethoxal was used for the first time in an electrochemical sensor. This work not only broadens the toolkit for detecting DNA modifications but also provides a fresh impetus for the development of point-of-care testing (POCT) technologies.

Elongation and Ligation-Mediated Differential Coding for Label-Free and Locus-Specific Analysis of 8-Oxo-7,8-dihydroguanine in DNA

Oxidative DNA damage is closely associated with the occurrence of numerous human diseases and cancers. 8-Oxo-7,8-dihydroguanine (8-oxoG) is the most prevalent form of DNA damage, and it has become not only an oxidative stress biomarker but also a new epigenetic-like biomarker. However, few approaches are available for the locus-specific detection of 8-oxoG because of the low abundance of 8-oxoG damage in DNA and the limited sensitivity of existing assays. Herein, we demonstrate the elongation and ligation-mediated differential coding for label-free and locus-specific analysis of 8-oxoG in DNA. This assay is very simple without the involvement of any specific labeled probes, complicated steps, and large sample consumption. The utilization of Bsu DNA polymerase can specifically initiate a single-base extension reaction to incorporate dATP into the opposite position of 8-oxoG, endowing this assay with excellent selectivity. The introduction of cascade amplification reaction significantly enhances the sensitivity. The proposed method can monitor 8-oxoG with a limit of detection of 8.21 × 10-19 M (0.82 aM), and it can identify as low as 0.001% 8-oxoG damage from a complex mixture with excessive undamaged DNAs. This method can be further applied to measure 8-oxoG levels in the genomic DNA of human cells under diverse oxidative stress, holding prospect potential in the dynamic monitoring of critical 8-oxoG sites, early clinical diagnosis, and gene damage-related biomedical research.

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.

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-Base Resolution Detection of N6-Methyladenosine in RNA by Adenosine Deamination Sequencing

N⁶-Methyladenosine (m⁶A) is one of the most abundant and prevalent natural modifications occurring in diverse RNA species. m⁶A plays a wide range of roles in physiological and pathological processes. Revealing the functions of m⁶A relies on the faithful detection of individual m⁶A sites in RNA. However, developing a simple method for the single-base resolution detection of m⁶A is still a challenging task. Herein, we report an adenosine deamination sequencing (AD-seq) technique for the facile detection of m⁶A in RNA at single-base resolution. The AD-seq approach capitalizes on the selective deamination of adenosine, but not m⁶A, by the evolved tRNA adenosine deaminase (TadA) variant of TadA8e or the dimer protein of TadA-TadA8e. In AD-seq, adenosine is deaminated by TadA8e or TadA-TadA8e to form inosine, which pairs with cytidine and is read as guanosine in sequencing. m⁶A resists deamination due to the interference of the methyl group at the N6 position of adenosine. Thus, the m⁶A base pairs with thymine and is still read as adenosine in sequencing. The differential readouts from A and m⁶A in sequencing can achieve the single-base resolution detection of m⁶A in RNA. Application of the proposed AD-seq successfully identified individual m⁶A sites in Escherichia coli 23S rRNA. Taken together, the proposed AD-seq allows simple and cost-effective detection of m⁶A at single-base resolution in RNA, which provides a valuable tool to decipher the functions of m⁶A in RNA.

Blocker-dUThiophene poly tailing-based method for assessing methyl transferase activity

In this report, we present a method for the selective and sensitive detection of methyl transferase activity. The method uses a dsDNA probe that contains C3 spacers and is coupled with dUThioTP-TdT polymerase-based poly-tailing. The short dsDNA probe is designed with C3 spacers at both 3' ends to prevent any type of tailing reaction. However, the probe contains a methyl transferase recognition sequence that can methylate adenosines in the palindromic part of both strands. When a specific DpnI endonuclease is introduced, it selectively cleaves the dsDNA probe such that both strands are methylated, unblocking the probe into two separate dsDNA forms with exposed 3' OH groups. This makes the probe susceptible to tailing in the presence of a TdT tailing polymerase. The unblocked probe is then subjected to fluorescent dUThioTP-based tailing, which produces a strong fluorescent signal that indicates the presence of methyl transferase activity. In the absence of methyl transferase, the probe remains in the blocked state and does not undergo fluorescence. This method has a limit of detection of 0.049 U/mL with good selectivity and the potential for accurate MTase analysis.

Enzymatic Cleavage-Mediated Extension Stalling Enables Accurate Recognition and Quantification of Locus-Specific Uracil Modification in DNA

Chemical modifications in DNA have profound influences on the structures and functions of DNA. Uracil, a naturally occurring DNA modification, can originate from the deamination of cytosine or arise from misincorporation of dUTP into DNA during DNA replication. Uracil in DNA will imperil genomic stability due to their potential in producing detrimental mutations. An in-depth understanding of the functions of uracil modification requires the accurate determination of its site as well as content in genomes. Herein, we characterized that a new member of the uracil-DNA glycosylase (UDG) family enzyme (UdgX-H109S) could selectively cleave both uracil-containing single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA). Based on this unique property of UdgX-H109S, we developed an enzymatic cleavage-mediated extension stalling (ECES) method for the locus-specific detection and quantification of uracil in genomic DNA. In the ECES method, UdgX-H109S specifically recognizes and cleaves the N-glycosidic bond of uracil from dsDNA and generates an apurinic/apyrimidinic (AP) site, which could be broken by APE1 to form a one-nucleotide gap. The specific cleavage by UdgX-H109S is then evaluated and quantified by qPCR. With the developed ECES approach, we demonstrated that the level of uracil at position Chr4:50566961 in genomic DNA of breast cancer tissues was significantly decreased. Collectively, the ECES method has been proved to be accurate and reproducible in the locus-specific quantification of uracil in genomic DNA from biological and clinical samples.

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.

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.

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.

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.

Highly Sensitive Fluorescence Detection of Global 5-Hydroxymethylcytosine from Nanogram Input with Strongly Emitting Copper Nanotags

Quantitative analysis of 5-hydroxymethylcytosine (5hmC) has remarkable clinical significance to early cancer diagnosis; however, it is limited by the requirement in current assays for large amounts of starting material and expensive instruments requring expertise. Herein, we present a highly sensitive fluorescence method, termed hmC-TACN, for global 5hmC quantification from several nanogram inputs based on terminal deoxynucleotide transferase (TdT)-assisted formation of fluorescent copper (Cu) nanotags. In this method, 5hmC is labeled with click tags by T4 phage β-glucosyltransferase (β-GT) and cross-linked with a random DNA primer via click chemistry. TdT initiates the template-free extension along the primer at the modified 5hmC site and then generates a long polythymine (T) tail, which can template the production of strongly emitting Cu nanoparticles (CuNPs). Consequently, an intensely fluorescent tag containing numerous CuNPs can be labeled onto the 5hmC site, providing the sensitive quantification of 5hmC with a limit of detection (LOD) as low as 0.021% of total nucleotides (S/N = 3). With only a 5 ng input (∼1000 cells) of genomic DNA, global 5hmC levels were accurately determined in mouse tissues, human cell lines (including normal and cancer cells of breast, lung, and liver), and urines of a bladder cancer patient and healthy control. Moreover, as few as 100 cells can also be distinguished between normal and cancer cells. The hmC-TACN method has great promise of being cost effective and easily mastered, with low-input clinical utility, and even for the microzone analysis of tumor models.

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

DNA cytosine methylation (5-methylcytosine, 5mC) is the most important epigenetic mark in higher eukaryotes. 5mC in genomes is dynamically controlled by writers and erasers. DNA (cytosine-5)-methyltransferases (DNMTs) are responsible for the generation and maintenance of 5mC in genomes. Active demethylation of 5-methylcytosine (5mC) is achieved by ten-eleven translocation (TET) dioxygenase-mediated oxidation of 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). 5fC and 5caC are further processed by thymine DNA glycosylase (TDG)-initiated base excision repair (BER) to restore unmodified cytosines. The TET-TDG-BER pathway could cause the production of DNA strand breaks and therefore jeopardize the integrity of genomes. Here, we investigated the direct decarboxylation of 5caC in mammalian genomes by using metabolic labeling with 2'-fluorinated 5caC (F-5caC) and mass spectrometry analysis. Our results clearly demonstrated the decarboxylation of 5caC occurring in mammalian genomes, which unveiled that, in addition to the TET-TDG-BER pathway, the direct decarboxylation of TET-produced 5caC constituted a new pathway for active demethylation of 5mC in mammalian genomes.