Although enzymes that catalyze DNA methylation process are well-studied (1
), how DNA demethylation is achieved is less known, especially in animals (2
). A repair-based mechanism is used in DNA demethylation in plants, but whether a similar mechanism is also used in mammalian cells is unclear (3
). Identification of 5hmC as the 6th
base of the mammalian genome (5
) and the capacity of Tet (ten_eleven_translocation) proteins to convert 5mC to 5hmC in an Fe(II) and alpha-ketoglutarate (α-KG)-dependent oxidation reaction (6
) raised the possibility that a Tet-catalyzed reaction might be part of the DNA demethylation process.
A potential 5mC demethylation mechanism can be envisioned from similar chemistry for thymine to uracil conversion (3
) (Fig. S1A
) with the Tet proteins oxidizing 5mC not only to 5hmC, but also to the aldehyde (5fC) and potentially the carboxylic acid (5caC) forms (Fig. S1B
). The failure to detect such reaction products may simply be due to the limitations of the previous assay employed (6
). To determine if this might be the case, we synthesized 20mer DNA oligos with 5fC or 5caC in the internal C of an MspI site (10
) and found that although MspI is efficient in digesting the oligo DNAs with C/5mC/5hmC in the internal C, it failed to digest the DNA containing 5fC or 5caC (Fig. S2A-B
). Thus, if Tet proteins have the capacity to convert 5mC to 5fC or 5caC, these products would have evaded detection due to the inability of MspI to digest 5fC or 5caC-containing DNA. To overcome this problem, we identified and demonstrated that TaqI is capable of digesting DNA modified with 5mC, 5hmC (11
), 5fC or 5caC (Fig. S2C-D
In addition to restriction enzyme, TLC conditions can also affect the detection of 5fC and 5caC. Under previous TLC conditions (7
), 5hmC and 5fC have almost identical migration patterns (, lanes 4 and 5) and 5caC failed to migrate (, lane 6). Using a more acidic TLC buffer, all cytosine derivatives migrated (, lanes 4-7). However, 5mC and C cannot be separated under this condition (, lanes 1 and 7). Given that the TLC buffer used in can separate C from 5mC, two-dimensional TLC (2D-TLC) using the two buffer conditions should allow for separation of cytosine and its derivatives.
Optimization of conditions for detection of cytosine and its 5-position modified forms by TLC
Using TaqI digestion and 2D-TLC (Fig. S3
), we analyzed the enzymatic activity of the Tet proteins. Compared with the mutant control, incubation of the Tet1 protein with 5mC-containing substrate resulted in a decrease in the 5mC level concomitant with the appearance of a radioactive spot that correlates with 5hmC (, left two panels). Two additional radioactive spots, labeled “X” and “Y” whose appearance depend on Tet1 enzymatic activity, were observed. Similarly, Tet2 and Tet3 also generated three enzymatic activity-dependent radioactive spots that were detected in Tet1-catalyzed reaction although the signal that corresponds to the “Y” spot from the Tet3 reaction is extremely weak (, middle and right panels).
If our hypothetical model for DNA demethylation is correct (Fig. S1B
), the “X” and “Y” spots are likely to be 5fC and 5caC. We compared the migration patterns of 5fC and 5caC with that of Tet2-treated 5mC-containing DNA substrates and found that the “X” and “Y” spots match 5fC and 5caC with respect to their migration (, compare the first two panels). We further confirmed this by mixing radioactive 5fC (third panel) or 5caC (last panel) with the samples used in the first panel before performing 2D-TLC. To confirm the identities of the “X” and “Y” spots, we treated the Tet2-catalyzed reaction mixture with sodium borohydride (NaBH4
), which resulted in the disappearance of both “X” and “Y” spots concomitant with increase in 5hmC (, compare the first two panels) indicating that both are oxidation products of 5hmC, consistent with the notion that they are 5fC and 5caC.
Tet proteins are capable of converting 5mC to 5hmC, 5fC, and 5caC
O-ethylhydroxylamine hydrochloride (EHL) and 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide hydrochloride (EDC) react with formyl and carboxyl groups to generate oximes and amides, respectively (12
) (Fig. S4
). To determine the migration patterns of the reaction products, we performed reactions using standard 5fC and 5caC, and separated the products by 2D-TLC, establishing the migration pattern for oxime (Fig. S4A
) and amide (Fig. S4B
). Similar EHL treatment of the Tet2 reaction mixture specifically converted the “X” spot to a new spot that co-migrated with oxime (, compare panels 1, 3 and Fig. S4A
). In contrast, EDC treatment specifically converted the “Y” spot to a new signal that co-migrated with amide (, compare panels 1, 4, and Fig. S4B
). To unequivocally define the identities of “X” and Y”, we employed mass spectrometry. Having established the mass spectrometry fingerprints of standard 5fC and 5caC (, top panels), we extracted the “X” and “Y” spots and subjected them to mass spectrometric analysis. The “X” spot shows the same major fragment ions as that of 5fC, while the “Y” spot shows the same major fragment ions as that of 5caC. Collectively, 2D-TLC co-migration, chemical treatment, and mass spectrometry fingerprints demonstrate that Tet proteins not only can convert 5mC to 5hmC, but also can further oxidize 5hmC to 5fC and 5caC.
To determine if Tet proteins can use 5hmC or 5fC-containing DNA as substrates, 20mer DNA oligos with either 5hmC or 5fC in the TaqI site were incubated with Tet proteins. 2D-TLC analysis demonstrated that incubation with wild-type Tet proteins, but not the catalytic mutants, resulted in a decrease in the level of 5hmC/5fC concomitant with the appearance of 5fC and 5caC, or 5caC (, S5
) suggesting that Tet proteins can act upon 5hmC and 5fC-containing substrates. However, the 5caC signal generated by Tet3 is extremely weak.
Kinetic analysis of Tet2 using 5mC, 5hmC, and 5fC-containing oligo DNAs
We used a quantitative mass spectrometric assay to rule out the possibility that 5fC and 5caC are generated as a side reaction by Tet proteins. We generated a standard curve for each of the cytosine derivatives by mixing different amounts of each 5mC, 5hmC, 5fC, and 5caC followed by LC-MS (Fig. S6
). We then quantified the cytosine derivatives at different time points after incubating Tet2 with 5mC, 5hmC, or 5fC-containing DNA substrates. Quantification of the relative amount of the substrate and the various products during the reaction process demonstrated that the reaction plateaued after 10 min of incubation regardless whether 5mC, 5hmC, or 5fC-containing TaqI 20mer DNA is used as a substrate (). The reaction plateaued in 10 min due to the inactivation of the Tet2 enzyme during the incubation (Fig. S7
During this period, Tet2 is able to convert more than 95% of the 5mC to 5hmC (~60%), 5fC (~30%), and 5caC (5%), but it can only convert about 40% or 25% when 5hmC or 5fC-contianing DNA was used as a substrate (). From this data we calculated the initial reaction rate of Tet2 for 5mC, 5hmC, and 5fC-containing substrates to be 429 nM/min, 87.4 nM/min, and 56.6 nM/min, respectively (Fig. S8
). Although Tet2 has a clear preference for the 5mC-containing DNA substrate, its initial reaction rate for 5hmC and 5fC-containing substrate is only 4.9-7.6 fold lower. The fact that there is clear accumulation of 5fC and 5caC when 5mC is used as a substrate (, top panel) strongly suggests that Tet-catalyzed iterative oxidation is likely a kinetically relevant pathway.
To determine whether Tet-catalyzed iterative oxidation of 5mC can take place in vivo
, we transfected a mammalian expression construct containing the Tet2 catalytic domain fused to GFP into HEK293 cells. After FACS sorting, genomic DNA of GFP positive cells was analyzed for the presence of 5hmC, 5fC, and 5caC by 2D-TLC (Fig. S3
). Compared with the untransfected control, cells expressing Tet2 not only have increased 5hmC levels, but also contain two additional spots (), which correspond to 5fC and 5caC, respectively. In addition, we quantified the genomic content of 5hmC, 5fC, and 5caC following the procedure depicted in Fig. S9A
). After establishing the retention times for each of the cytosine derivatives on HPLC (Fig. S9B
, top panel), nucleosides derived from genomic DNA were subjected to the same HPLC conditions for fractionation. Fractions A and B (Fig. S9B
) that have the same retention times as that of 5caC and 5hmC or 5fC were collected. Mass spectrometry analysis demonstrates that both 5fC and 5caC are detected in the genomic DNA of cells overexpressing Tet2 (Fig. S10A
). By comparison to the standard curves (Fig. S11A
), overexpression of wild-type Tet2, but not a catalytic mutant, increased the genomic content of 5hmC, 5fC and 5caC ().
5fC and 5caC are present in genomic DNA and their abundance is regulated by Tet proteins
Next, we asked whether 5fC and 5caC are present in genomic DNA under physiological conditions. Using a similar approach as that used for the genomic DNA of Tet2-overexpressing HEK293, we show that not only 5hmC, but also 5fC and 5caC are present in the genomic DNA of mouse ES cells (Fig. S10B
). To quantify the genomic content of 5hmC, 5fC and 5caC in mouse ES cells, we generated standard curves for each of the 5mC derivatives at low concentrations and determined the limit of detection for 5fC and 5caC to be 5 fmol and 10 fmol, respectively (Fig. S11
). We then quantified the genomic content of these cytosine derivatives in mouse ES cells to be about 1.3×103
5hmC, 20 5fC, and 3 5caC in every 106
C (, Table S1
). Knockdown of Tet1 reduced the genomic content of 5hmC, as well as 5fC and 5caC () indicating that Tet1 is at least partially responsible for the generation of these cytosine derivatives. The presence of 5fC is not limited to ES cells as similar analysis also revealed their presence in genomic DNA of major mouse organs (). However, 5caC can be detected with confidence only in ES cells (, S10B
Here we demonstrate that the Tet family of proteins have the capacity to convert 5mC not only to 5hmC, but also to 5fC and 5caC in vitro
. In addition, we provide evidence for the presence of 5fC in the genomic DNA of mouse ES cells and organs and for the presence of 5caC in moue ES cells. We note that a similar study failed to detect their existence in genomic DNA of mouse organs (15
) likely due to the differences in the detection limits between the two studies (pmol vs fmol). The Tet-catalyzed oxidation reaction is reminiscent of the thymine hydroxylase catalyzed conversion of thymine to iso-orotate (8
) (Fig. S1
) raising the possibility that 5mC demethylation could be potentially achieved through a process similar to the conversion of thymine to uracil, which is achieved by conversion of thymine to iso-orotate followed by decarboxylation by the iso-orotate decarboxylase (8
). Although this hypothetic pathway for DNA demethylation is simple and appealing, the enzyme that is capable of decarboxylating 5caC-containing DNA has yet to be identified. Until such an enzyme is identified, we cannot rule out the possibility that the Tet family enzymes act together with other putative DNA demethylation pathways, such as the base excision DNA repair (BER) pathway. Indeed, recent studies have provided some supporting evidence for such a possibility (16