- Open Access
A screen for hydroxymethylcytosine and formylcytosine binding proteins suggests functions in transcription and chromatin regulation
© Iurlaro et al.; licensee BioMed Central Ltd. 2013
Received: 17 September 2013
Accepted: 24 October 2013
Published: 10 December 2013
DNA methylation (5mC) plays important roles in epigenetic regulation of genome function. Recently, TET hydroxylases have been found to oxidise 5mC to hydroxymethylcytosine (5hmC), formylcytosine (5fC) and carboxylcytosine (5caC) in DNA. These derivatives have a role in demethylation of DNA but in addition may have epigenetic signaling functions in their own right. A recent study identified proteins which showed preferential binding to 5-methylcytosine (5mC) and its oxidised forms, where readers for 5mC and 5hmC showed little overlap, and proteins bound to further oxidation forms were enriched for repair proteins and transcription regulators. We extend this study by using promoter sequences as baits and compare protein binding patterns to unmodified or modified cytosine using DNA from mouse embryonic stem cell extracts.
We compared protein enrichments from two DNA probes with different CpG composition and show that, whereas some of the enriched proteins show specificity to cytosine modifications, others are selective for both modification and target sequences. Only a few proteins were identified with a preference for 5hmC (such as RPL26, PRP8 and the DNA mismatch repair protein MHS6), but proteins with a strong preference for 5fC were more numerous, including transcriptional regulators (FOXK1, FOXK2, FOXP1, FOXP4 and FOXI3), DNA repair factors (TDG and MPG) and chromatin regulators (EHMT1, L3MBTL2 and all components of the NuRD complex).
0ur screen has identified novel proteins that bind to 5fC in genomic sequences with different CpG composition and suggests they regulate transcription and chromatin, hence opening up functional investigations of 5fC readers.
Levels of 5hmC in DNA (and where known 5fC and 5caC) vary between different mammalian tissues and are highest in ES cells and neural tissues [1–5]. In situations where oxidative derivatives of 5mC are implicated in demethylation of DNA, such as in pluripotent stem cells, early embryos and germ cells, there may be rapid turnover of these modifications through a combination of further oxidation, DNA replication, excision repair by TDG, and potentially deamination or decarboxylation [6–8]. In other tissues, especially those with non-dividing cells such as neural tissues, the modifications could potentially be more stable and might thus be used as epigenetic signals for genome function [9–11]. A variety of proteins that bind to histone modifications or to methylated DNA (methyl binding domain proteins (MBDs)) have been described and have a role in interpreting these epigenetic signals for the regulation of transcription, replication, DNA repair or other functions of the genome [12–14]. Recently, MBD3 and MECP2 have been shown to be able to bind 5hmC (MBD3 weakly so) in addition to 5mC, opening up the possibility that these proteins may also be able to interpret the 5hmC signal, for example, in the regulation of transcription or chromatin [15, 16]. A recently published unbiased screen  has identified and validated a number of proteins with specific binding to 5mC and its oxidised forms but the use of a single DNA probe overlooks the possibility that proteins in a cellular context might have a combined preference for both DNA modification and sequence context. Indeed some of the proteins identified as specific for a DNA modification are cell-type specific, suggesting a complex protein interaction network operating in modulating the intrinsic ability to bind to DNA modifications.
Results and discussion
Gene ontology term enrichment comparing modification specific binders to the full set of identified proteins showed highly significant groups enriching with relevance to gene transcription and chromatin regulation among 5fC binders on the Fgf15 probe (Figure 2b). Association of 5fC with repressive transcription complexes was a surprising finding where, notably, all members of the core NuRD complex were enriched in the group of 5fC specific binding proteins (Figure 2a), although it is likely that some of the members of the complex are not direct 5fC binders but are enriched by secondary protein-protein interactions. This indicates that 5fC is more likely to be associated with gene repression. Interestingly, many of the proteins enriched for 5fC at the Fgf15 probe were enriched for 5mC too, as seen by the hierarchical clustering, strengthening the potentially repressive properties of 5fC especially in the context of a CpG island sequence. This was not the case with the Pax6 probe, which is not a CpG island (Figure 3). It remains to be seen if the presence of 5fC in CGIs has inhibitory functions, especially in the process of cell differentiation. Clustering of proteins enriching on the Pax6 probe did not result in a similar grouping of 5fC and 5mC enriching repressive proteins and the GO analysis showed no significant enrichment for repressive complexes indicating that the DNA sequence of Pax6 might lack the DNA signatures of a typical CpG island therefore may not result in an inhibitory transcriptional signal in the presence of 5fC. While our experimental system made use of a promoter CpG island (in Fgf15) these insights may also be applicable to intragenic CpG islands, which can have higher levels of DNA modifications . The association between 5-formylcytosine and transcription has been investigated recently, resulting in its linkage variously with active or poised genes [25–27]. Our results potentially reinforce the idea that depending on context 5fC could have positive or negative effects on transcription. Nevertheless, some of the 5fC specific proteins were enriched with both DNA probes and are shown in Figure 4a. This comparison strongly suggests that Fork head box domain containing proteins have 5fC binding properties. Gene ontology results for the other cytosine modifications for the two probes are included as Additional files 3 and 4.
In order to validate some of these candidate proteins for 5fC binding specificity, we performed ELISA with purified recombinant proteins and differentially modified Fgf15 probes. His-tagged isoforms of MPG, L3MBTL2 and ZSCAN21 were expressed in Sf9 insect cells using a Baculovirus system, and purified by immobilised metal ion affinity chromatography (IMAC). We found that all three proteins bound with higher affinity to 5fC compared to the other modifications on the DNA (Figure 4b). MPG is one of the proteins common for both DNA targets and showed a strong binding preference for 5fC. In a recent study MPG was identified as a 5hmC specific binder but the data actually show some binding to 5fC as well , and considering different culturing conditions of ES cells (2i/LIF), post-translational modifications might modulate the binding of some proteins to their target . Finally, we considered the possibility that the 5fC binding proteins might have a role in the excision of 5fC similar to TDG. We therefore tested this hypothesis by RNAi in ES cells (Figure 4c, Additional files 5 and 6). While knockdown of TDG (which is known to excise 5fC and 5caC [29, 30]) resulted in increase of 5fC and 5caC (as measured by mass spectrometry), knockdown of the other candidates had no effect. We therefore conclude that the majority of 5fC binding proteins identified in this screen are less likely to metabolize 5fC, instead they are more likely to recognize 5fC as an epigenetic signal.
The preferential binding of TET1 to both 5mC (more strongly) and 5hmC, compared to C (Figure 3) was interesting since the CXXC domain of TET1 has been shown to differ from that of other CXXC domain-containing proteins, lacking a typical 'KFGG’ motif found in most of the family, with some studies showing its inability to bind DNA , and others suggesting that this peculiarity allows it to bind not only to unmodified and methylated DNA, but also to hydroxymethylated DNA [32, 33]. This opens the possibility that the binding could be influenced by sequence context or protein modifications.
It was of particular interest that our screen identified a higher number of proteins that appear to preferentially bind to 5fC (Figure 1c,d) rather than to other modifications, an observation also reported in Spruijit et al. . It is not immediately intuitive why there should be more proteins binding to 5fC than to 5hmC. Of course this could depend on the tissue analysed and there might be more 5hmC binding proteins in neural cell types, for example, where the modification is relatively prevalent. Intriguingly, FOXK2 in addition to being a member of the forkhead box transcription factor family has been shown to bind to T:G mismatches in DNA but no enzymatic activity has been identified . Another member of this family, FOXP1, a key transcriptional regulator in B cells and lung development was also identified as strong and specific 5fC binder in our screen. Recent reports have shown that an ES cell-specific isoform of FOXP1 is implicated in pluripotency regulation in ESCs by stimulating expression of pluripotency-related genes like Oct4, Nanog and Nr5a2. FOXP4, also enriched on both 5fC probes, is involved in development of the lung and is known to form homodimers and heterodimers with FOXP1, and to interact with NuRD components . FOXK1 is a transcriptional regulator involved in myogenic regulation , while relatively little is known about the function of mouse FOXI3. Another transcription factor that appears to bind specifically to 5fC in our screen is ZSCAN21, a strong transcriptional activator that plays a role in both male and female meiosis [38, 39]. The final protein in this category of transcriptional regulation linked with DNA repair is MPG, which is a base excision repair glycosylase known to excise modified bases resulting from alkylation damage. MPG was a highly specific binder for 5fC in our screen, while the human isoform bound strongly to 5fC in a HeLa sample extract providing an additional layer of confidence (data not shown); MPG has been identified as a interacting partner of MBD1  and, intriguingly, its methyl-purine glycosylase domain structurally resembles the formyl transferase, C-terminal-like domain (IPR011034).
The last category of 5fC binders makes interesting connections with chromatin regulation through the polycomb and histone methylation pathways. In addition to the previously mentioned correlation between 5fC and the NuRD complex, components of another chromatin regulator complex, E2F6.com-1, were also identified as 5fC binders. In addition to MGA and CBX3, we isolated and verified L3MBTL2 as a 5fC binder, which is a putative polycomb protein which may bind to modified histones, while EHMT1 is a euchromatin histone methyltransferase that methylates H3K9 to H3K9me1 and me2, potentially providing a link between modifications in euchromatin that are intermediates between transcriptional repression and activation [41, 42].
We have established a relatively simple and robust screen for proteins that bind 5hmC and 5fC in DNA. 5fC has so far been found in early embryos, embryonic stem cells and brain cortex, as well as in other major mouse organs like spleen, pancreas and liver . The distribution of 5fC in ESCs depends on TDG and recent studies have linked it with the regulation of transcription, variously associated with active or poised genes [25–27]. Our screen has identified 5fC-binding proteins with functions in transcription and in chromatin regulation, particularly involving forkhead box domain transcriptional regulators and the NuRD complex. This suggests that 5fC may be both an intermediate in demethylation and an epigenetic signal in its own right. The dual potential of some of the proteins we have identified (FOXK2 in transcription and DNA repair, EHMT1 mediating between 5fC and H3K9 methylation) is particularly interesting and warrants future functional investigations.
Cell lines and cell culture
E14 ES cells (derived from the E14 cell line strain 129P2/OlaHsd) were grown on a γ-irradiated pMEF feeder layer at 37°C and 5% CO2 in complete ES medium (DMEM 4,500 mg l-1 glucose, 4 mM l-glutamine and 110 mg l-1 sodium pyruvate, 15% fetal bovine serum, 100 U of penicillin/100 μg of streptomycin in 100 mL medium, 0.1 mM non-essential amino acids, 50 μM β-mercaptoethanol, 103 U LIF ESGRO).
Cells were washed with 1× PBS solution, detached adding trypsin at 37°C to the culture plate and centrifuged at 300 × g for 4 min. The pellet was then washed in ice-cold 1× PBS twice and resuspended gently in 5 volumes of ice-cold 1 Cytoplasmic Lysis Buffer (Chemicon International®) containing 0.5 mM DTT and 1/1,000 dilution of supplied protease inhibitor Cocktail. The solution was incubated on ice for 15 min, centrifuged at 300 × g for 5 min at 4°C, and the pellet was resuspended in two volumes of ice-cold 1× Cytoplasmic Lysis Buffer. Cells were lysed using a 27-gauge needle and the nuclear fraction was isolated from the cytosolic portion by centrifugation at 8,000 × g for 20 min at 4°C. Finally, the pellet was resuspended in two-thirds of the original cell pellet volume of ice-cold Nuclear Extraction Buffer (Chemicon International®) containing 0.5 mM DTT and 1/1,000 dilution of supplied protease inhibitor cocktail, incubated on orbital shaker for 60 min at 4°C, and centrifuged at 16,000 × g for 5 min at 4°C. The nuclear extract was then aliquoted and stored at -80°C.
The probes were obtained by PCR amplification of genomic region corresponding to the promoters of Pax6 (280 bp) and Fgf15 (248 bp) genes using DreamTaq™ DNA Polymerase (Fermentas). The primers used in the reaction were:
Pax6-F (Biotinylated): ATTCCCAAAGCAAGCAGAAG
Fgf15-F (Biotinylated): TTTCTTTCAGGCAGGGGAAT
The pull-down assay was carried out using Dynabeads® M-280 Streptavidin (Invitrogen™). For each sample, 2 μL of beads were washed in buffer PBT (1× PBS, 0.1% Triton X-100), and incubated with 50 ng of biotinylated DNA in 200 uL of PBS, overnight at 4°C. The beads were then washed three times in PBT and twice in buffer D-T (0.2 mM EDTA, 20% Glycerol, 20 mM Hepes-KOH pH 7.9, 0.1 M KCl, 1 mM DTT, 1 mM protease inhibitor PMSF, 0.1% Triton X-100), and incubated with 50 μg of nuclear extract for 15 min at 4°C in incubation buffer (0.05 mM EDTA, 5% Glycerol, 5 mM Hepes-KOH pH 7.9. 150 mM KCl, 1 mM DTT, 1 mM protease inhibitor PMSF, 0.025% Triton X-100 in PBS). The beads were washed six times in Buffer D-T, once in PBS and eluted in 1X LDS Loading buffer boiling at 95°C for 5 min. The eluted fraction was separated from the beads and finally analysed by mass spectrometry.
RNAi knockdown of Mpg, Tdg, L3mbtl2, Zscan21, Ehmt1, FoxK2 and FoxP1 in ES cells
Transfections of Dharmacon siGENOME SMARTpool against mouse Tdg (catalogue number M-040666-01; gaagugcaguauacauuug, gaguaaagguuaagaacuu, caaagaagauggcuguuaa, gcaaggaucugucuaguaa) and siGENOME ON-TARGETplus siRNA against Mpg (catalogue no. J-060513-11; ccggcuaggaccagaguuu), L3mbtl2 (catalogue no. J-065321-12; uuacugacuggaagagcua), FoxP1 (catalogue no. J-065400-09; gagcaugcgcuggacgaua), Ehmt1 (catalogue no. J-059041-12; gagcacagguggauccgaa), Zscan21/Zipro1 (catalogue no. J-048225-09; cuagagauaucccguaaga), FoxK2 (catalogue no. J-064514-12; ccagagcucaagcgaguua) were done with Lipofectamine 2000 according to the manufacturer’s instructions. Cells were harvested after three rounds of transfection for DNA/RNA isolation.
Eluted proteins were run a short distance (approximately 5 mm) into an SDS-PAGE gel, which was then stained with colloidal Coomassie stain (Imperial Blue, Invitrogen). The entire stained gel pieces were excised, then destained, reduced, carbamidomethylated and digested overnight with trypsin (Promega sequencing grade, 10 ng/μL in 25 mM ammonium bicarbonate) as previously described . Aliquots of each of the resulting tryptic digests were analysed by LC-MS/MS on a system comprising a nanoLC (Proxeon) coupled to a LTQ Orbitrap Velos mass spectrometer (Thermo). LC separation was achieved on a reversed-phase column (Reprosil C18AQ, 0.075 × 150 mm, 3 μm particle size), with an acetonitrile gradient (0-35% over 60 min, containing 0.1% formic acid, at a flow rate of 300 nL/min). The mass spectrometer was operated in data-dependent acquisition mode, with an acquisition cycle consisted of a high resolution precursor ion spectrum over the m/z range 350–1,500, followed by up to 20 CID spectra (with a 60 s dynamic exclusion of former target ions). Mass spectrometric data were searched against a database generated from the mammalian entries in Uniprot 2011.09 by concatenation of the forward and reversed sequences, using Mascot (Matrix Science) and the search results were processed using Scaffold software (Proteome Software Inc.). Criteria for protein identification were: minimum of two peptides, each with a probability of >50% and an overall protein probability of >99%, which gave a protein false discovery rate of 0.4%. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium  via the PRIDE partner repository  with the dataset identifier PXD000524.
Pulled-down proteins were eluted from beads in LDS Loading buffer, boiled and run on NuPAGE® Novex 4-12% Bis-Tris Gel 1.0 mm (Novex®). Proteins were transferred on a nitrocellulose membrane using iBlot® Blotting System (Life Technologies), membrane was blocked overnight in PBS-0.1%Tween (PBST) containing 5% BSA (blocking buffer). Primary antibody incubation was done at room temperature for 2 h with a rabbit polyclonal anti-UHRF1 Antibody (Santa Cruz M-132: sc-98817). Membrane was washed in PBST and incubated with HRP conjugated anti-rabbit secondary antibody in blocking buffer. HRP conjugates were detected with enhanced chemiluminescence (ECL, Amersham Biosciences).
Enzyme-linked immunosorbent assay (ELISA)
All binding reactions were carried out in buffer Z containing 20 mM TRIS HCL (pH 7.5), 150 mM NaCl, 20 mM KCl, 0.02% IGEPAL and 1 mM dithiothreitol. A Highbind Streptaplate (Roche) was blocked with 1 × PBS containing 3% BSA prior reaction. Subsequently, 50 μL of a 50 nM solution of biotinylated DNA were added per well and allowed to attach for 30 min at 37°C with gentle shaking. Wells were then washed three times with buffer Z. The proteins were diluted in buffer Z and 50 μL were added to each well. After incubation for 1 h at room temperature, plates were washed three times with buffer Z. For detection, 50 μL of mouse polyclonal anti-His tag antibody (Thermo Scientific) at 1:500 dilution in buffer Z were added per well and incubated for 1 h at room temperature. After washing three times with buffer Z, a polyclonal HRP-conjugated sheep anti-Mouse IgG antibody (GE Healthcare) diluted 1:2,000 in buffer Z was added and incubated for 30 min at room temperature. Wells were washed three times with buffer Z and peroxidase activity detected by adding 50 μL of TACS-Sapphire (Trevigen). Reactions were stopped by the addition of 50 uL of a 1 M HCL solution. Absorbance at 450 nm was measured using a SPECTROstar Nano (BMG Labtech). The equilibrium dissociation constants (Kd) for the protein-DNA interaction were determined by non-linear regression by fitting to a hyperbolic binding curve.
Purification of recombinant MPG, L3MBTL2 and ZSCAN21 from Baculovirus infected Sf9 cells
Coding sequences for the proteins MPG, L3MBTL2 and ZSCAN21 (Source BioScience) were cloned into Gateway® entry vector pENTR223.1 using SfiI restriction sites. CDS were then cloned into destination vector pDEST10 using Gateway® LR Clonase II mix (Invitrogen) and following manufacturer’s instructions. Resulting vectors were used to transform MAX Efficiency® DH10Bac™ cells (Invitrogen). Positive clones were selected by blue-complementation and correct insertion of sequence of interest was confirmed by PCR. Resulting bacmids were then transfected into Sf9 cells using Cellfectin® II Reagent (Invitrogen). Baculoviruses were then amplified and Sf9 cells expressing the proteins of interest were then harvested at 48, 72 and 96 h post infection for protein expression analysis. Cells pellets were resuspended in Lysis Buffer (50 mM NaH2PO4, 300 mM NaCl, 10 mM imidazole, 1% Triton and protease inhibitors), incubated on ice for 10 min and centrifuged at 10,000 × g for 10 min at 4°C. Cell lysates were filtered through a 0.2 μm filter and loaded on 1 mL HisTrap HP column (GE Healthcare) equilibrating with buffer A (50 mM NaH2PO4, 300 mM NaCl, 20 mM imidazole), washed with 10 column volumes of buffer A added with 40 mM imidazole. Proteins were eluted with a gradient of 40–500 mM imidazole over 20 column volumes. Protein samples were dialysed against storage buffer (25 mM Tris–HCl pH 7.5 10% glycerol, 150 mM NaCl, 1 mM DTT).
Spectral count values from LC-MS/MS were analysed by testing the sample groups using a non-parametric Kruskal Wallis t-test with a P value cutoff of 0.1, which was determined to be sufficient to identify any group where the most extreme values all fell within that group, regardless of how the values were distributed across the other groups.
Mass spectrometry of nucleosides
Quantitation of nucleosides in genomic DNA was done essentially as described previously  except that a Q-Exactive mass spectrometer (Thermo) fitted with an UltiMate 3000 RSLCnano HPLC (Dionex) was used and one additional transition 272.1 >156.0404 (caC) was monitored. Results are expressed as % or ppm of total unmodified and modified cytosines.
MI is supported by the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme FP7/2007-2013/under REA grant agreement no. 290123 and was supported by Unipharma-Graduates 7 Da Vinci Programme. MJB is supported by a BBRSC studentship. The WR lab is supported by BBSRC, MRC, the Wellcome Trust, EU EpiGeneSys and BLUEPRINT. The SB lab is supported by core funding from Cancer Research UK and a Wellcome Trust Senior Investigator Award. We would like to thank Judith Webster for the preparation of samples for mass spectrometry, Patrick Varga-Weisz and Sarah Elderkin for help with chromatography, Maureen Hamon for Baculovirus work, Phil Ewels for bioinformatic analysis.
- Ito S, D’Alessio AC, Taranova OV, Hong K, Sowers LC, Zhang Y: Role of Tet proteins in 5mC to 5hmC conversion, ES-cell self-renewal and inner cell mass specification. Nature. 2010, 466: 1129-1133. 10.1038/nature09303.View ArticlePubMedPubMed CentralGoogle Scholar
- Kriaucionis S, Heintz N: The nuclear DNA base 5-hydroxymethylcytosine is present in Purkinje neurons and the brain. Science. 2009, 324: 929-930. 10.1126/science.1169786.View ArticlePubMedPubMed CentralGoogle Scholar
- Szwagierczak A, Bultmann S, Schmidt CS, Spada F, Leonhardt H: Sensitive enzymatic quantification of 5-hydroxymethylcytosine in genomic DNA. Nucleic Acids Res. 2010, 38: e181-10.1093/nar/gkq684.View ArticlePubMedPubMed CentralGoogle Scholar
- Tahiliani M, Koh KP, Shen Y, Pastor WA, Bandukwala H, Brudno Y, Agarwal S, Iyer LM, Liu DR, Aravind L, Rao A: Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science. 2009, 324: 930-935. 10.1126/science.1170116.View ArticlePubMedPubMed CentralGoogle Scholar
- Khare T, Pai S, Koncevicius K, Pal M, Kriukiene E, Liutkeviciute Z, Irimia M, Jia PX, Ptak C, Xia MH, Tice R, Tochigi M, Morera S, Nazarians A, Belsham D, Wong AHC, Blencowe BJ, Wang SC, Kapranov P, Kustra R, Labrie V, Klimasauskas S, Petronis A: 5-hmC in the brain is abundant in synaptic genes and shows differences at the exon-intron boundary. Nat Struct Mole Biol. 2012, 19: 1037-U1094. 10.1038/nsmb.2372.View ArticleGoogle Scholar
- Branco MR, Ficz G, Reik W: Uncovering the role of 5-hydroxymethylcytosine in the epigenome. Nat Rev Genet. 2012, 13: 7-13.Google Scholar
- Zhu JK: Active DNA demethylation mediated by DNA glycosylases. Annu Rev Genet. 2009, 43: 143-166. 10.1146/annurev-genet-102108-134205.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu SC, Zhang Y: Active DNA demethylation: many roads lead to Rome. Nat Rev Mol Cell Biol. 2010, 11: 607-620. 10.1038/nrm2950.View ArticlePubMedPubMed CentralGoogle Scholar
- Jin SG, Wu X, Li AX, Pfeifer GP: Genomic mapping of 5-hydroxymethylcytosine in the human brain. Nucleic Acids Res. 2011, 39: 5015-5024. 10.1093/nar/gkr120.View ArticlePubMedPubMed CentralGoogle Scholar
- Szulwach KE, Li X, Li Y, Song CX, Wu H, Dai Q, Irier H, Upadhyay AK, Gearing M, Levey AI, Vasanthakumar A, Godley LA, Chang Q, Cheng X, He C, Jin P: 5-hmC-mediated epigenetic dynamics during postnatal neurodevelopment and aging. Nat Neurosci. 2011, 14: 1607-1616. 10.1038/nn.2959.View ArticlePubMedPubMed CentralGoogle Scholar
- Spruijt CG, Gnerlich F, Smits AH, Pfaffeneder T, Jansen PW, Bauer C, Munzel M, Wagner M, Muller M, Khan F, Eberl HC, Mensinga A, Brinkman AB, Lephikov K, Muller U, Walter J, Boelens R, van Ingen H, Leonhardt H, Carell T, Vermeulen M: Dynamic readers for 5-(hydroxy)methylcytosine and its oxidized derivatives. Cell. 2013, 152: 1146-1159. 10.1016/j.cell.2013.02.004.View ArticlePubMedGoogle Scholar
- Rando OJ: Combinatorial complexity in chromatin structure and function: revisiting the histone code. Curr Opin Genet Dev. 2012, 22: 148-155. 10.1016/j.gde.2012.02.013.View ArticlePubMedPubMed CentralGoogle Scholar
- Law JA, Jacobsen SE: Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010, 11: 204-220. 10.1038/nrg2719.View ArticlePubMedPubMed CentralGoogle Scholar
- Deaton AM, Bird A: CpG islands and the regulation of transcription. Genes Dev. 2011, 25: 1010-1022. 10.1101/gad.2037511.View ArticlePubMedPubMed CentralGoogle Scholar
- Yildirim O, Li R, Hung JH, Chen PB, Dong X, Ee LS, Weng Z, Rando OJ, Fazzio TG: Mbd3/NURD complex regulates expression of 5-Hydroxymethylcytosine marked genes in embryonic stem cells. Cell. 2011, 147: 1498-1510. 10.1016/j.cell.2011.11.054.View ArticlePubMedPubMed CentralGoogle Scholar
- Mellen M, Ayata P, Dewell S, Kriaucionis S, Heintz N: MeCP2 binds to 5hmC enriched within active genes and accessible chromatin in the nervous system. Cell. 2012, 151: 1417-1430. 10.1016/j.cell.2012.11.022.View ArticlePubMedPubMed CentralGoogle Scholar
- Meehan RR, Lewis JD, McKay S, Kleiner EL, Bird AP: Identification of a mammalian protein that binds specifically to DNA containing methylated CpGs. Cell. 1989, 58: 499-507. 10.1016/0092-8674(89)90430-3.View ArticlePubMedGoogle Scholar
- Ficz G, Branco MR, Seisenberger S, Santos F, Krueger F, Hore TA, Marques CJ, Andrews S, Reik W: Dynamic regulation of 5-hydroxymethylcytosine in mouse ES cells and during differentiation. Nature. 2011, 473: 398-402. 10.1038/nature10008.View ArticlePubMedGoogle Scholar
- Frauer C, Hoffmann T, Bultmann S, Casa V, Cardoso MC, Antes I, Leonhardt H: Recognition of 5-hydroxymethylcytosine by the Uhrf1 SRA domain. PLoS One. 2011, 6: e21306-10.1371/journal.pone.0021306.View ArticlePubMedPubMed CentralGoogle Scholar
- Bartels SJ, Spruijt CG, Brinkman AB, Jansen PW, Vermeulen M, Stunnenberg HG: A SILAC-based screen for Methyl-CpG binding proteins identifies RBP-J as a DNA methylation and sequence-specific binding protein. PLoS One. 2011, 6: e25884-10.1371/journal.pone.0025884.View ArticlePubMedPubMed CentralGoogle Scholar
- Hendrich B, Bird A: Identification and characterization of a family of mammalian methyl-CpG binding proteins. Mol Cell Biol. 1998, 18: 6538-6547.View ArticlePubMedPubMed CentralGoogle Scholar
- Martello G, Sugimoto T, Diamanti E, Joshi A, Hannah R, Ohtsuka S, Gottgens B, Niwa H, Smith A: Esrrb is a pivotal target of the Gsk3/Tcf3 axis regulating embryonic stem cell self-renewal. Cell Stem Cell. 2012, 11: 491-504. 10.1016/j.stem.2012.06.008.View ArticlePubMedPubMed CentralGoogle Scholar
- Hashimoto H, Zhang X, Cheng X: Excision of thymine and 5-hydroxymethyluracil by the MBD4 DNA glycosylase domain: structural basis and implications for active DNA demethylation. Nucleic Acids Res. 2012, 40: 8276-8284. 10.1093/nar/gks628.View ArticlePubMedPubMed CentralGoogle Scholar
- Illingworth RS, Gruenewald-Schneider U, Webb S, Kerr AR, James KD, Turner DJ, Smith C, Harrison DJ, Andrews R, Bird AP: Orphan CpG islands identify numerous conserved promoters in the mammalian genome. PLoS Genet. 2010, 6: e1001134-10.1371/journal.pgen.1001134.View ArticlePubMedPubMed CentralGoogle Scholar
- Raiber EA, Beraldi D, Ficz G, Burgess HE, Branco MR, Murat P, Oxley D, Booth MJ, Reik W, Balasubramanian S: Genome-wide distribution of 5-formylcytosine in embryonic stem cells is associated with transcription and depends on thymine DNA glycosylase. Genome Biol. 2012, 13: R69-10.1186/gb-2012-13-8-r69.View ArticlePubMedPubMed CentralGoogle Scholar
- Song CX, Szulwach KE, Dai Q, Fu Y, Mao SQ, Lin L, Street C, Li Y, Poidevin M, Wu H, Gao J, Liu P, Li L, Xu GL, Jin P, He C: Genome-wide profiling of 5-formylcytosine reveals its roles in epigenetic priming. Cell. 2013, 153: 678-691. 10.1016/j.cell.2013.04.001.View ArticlePubMedPubMed CentralGoogle Scholar
- Shen L, Wu H, Diep D, Yamaguchi S, D’Alessio AC, Fung HL, Zhang K, Zhang Y: Genome-wide analysis reveals TET- and TDG-dependent 5-methylcytosine oxidation dynamics. Cell. 2013, 153: 692-706. 10.1016/j.cell.2013.04.002.View ArticlePubMedPubMed CentralGoogle Scholar
- Ficz G, Hore TA, Santos F, Lee HJ, Dean W, Arand J, Krueger F, Oxley D, Paul YL, Walter J, Cook SJ, Andrews S, Branco MR, Reik W: FGF Signaling Inhibition in ESCs Drives Rapid Genome-wide Demethylation to the Epigenetic Ground State of Pluripotency. Cell Stem Cell. 2013, 13: 351-359. 10.1016/j.stem.2013.06.004.View ArticlePubMedPubMed CentralGoogle Scholar
- Maiti A, Drohat AC: Thymine DNA glycosylase can rapidly excise 5-formylcytosine and 5-carboxylcytosine: potential implications for active demethylation of CpG sites. J Biol Chem. 2011, 286: 35334-35338. 10.1074/jbc.C111.284620.View ArticlePubMedPubMed CentralGoogle Scholar
- Hashimoto H, Hong S, Bhagwat AS, Zhang X, Cheng X: Excision of 5-hydroxymethyluracil and 5-carboxylcytosine by the thymine DNA glycosylase domain: its structural basis and implications for active DNA demethylation. Nucleic Acids Res. 2012, 40: 10203-10214. 10.1093/nar/gks845.View ArticlePubMedPubMed CentralGoogle Scholar
- Frauer C, Rottach A, Meilinger D, Bultmann S, Fellinger K, Hasenoder S, Wang M, Qin W, Soding J, Spada F, Leonhardt H: Different binding properties and function of CXXC zinc finger domains in Dnmt1 and Tet1. PLoS One. 2011, 6: e16627-10.1371/journal.pone.0016627.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang H, Zhang X, Clark E, Mulcahey M, Huang S, Shi YG: TET1 is a DNA-binding protein that modulates DNA methylation and gene transcription via hydroxylation of 5-methylcytosine. Cell Res. 2010, 20: 1390-1393. 10.1038/cr.2010.156.View ArticlePubMedGoogle Scholar
- Xu Y, Wu F, Tan L, Kong L, Xiong L, Deng J, Barbera AJ, Zheng L, Zhang H, Huang S, Min J, Nicholson T, Chen T, Xu G, Shi Y, Zhang K, Shi YG: Genome-wide regulation of 5hmC, 5mC, and gene expression by Tet1 hydroxylase in mouse embryonic stem cells. Mol Cell. 2011, 42: 451-464. 10.1016/j.molcel.2011.04.005.View ArticlePubMedPubMed CentralGoogle Scholar
- Fujii Y, Nakamura M: FOXK2 transcription factor is a novel G/T-mismatch DNA binding protein. J Biochem. 2010, 147: 705-709. 10.1093/jb/mvq004.View ArticlePubMedGoogle Scholar
- Gabut M, Samavarchi-Tehrani P, Wang X, Slobodeniuc V, O’Hanlon D, Sung HK, Alvarez M, Talukder S, Pan Q, Mazzoni EO, Nedelec S, Wichterle H, Woltjen K, Hughes TR, Zandstra PW, Nagy A, Wrana JL, Blencowe BJ: An alternative splicing switch regulates embryonic stem cell pluripotency and reprogramming. Cell. 2011, 147: 132-146. 10.1016/j.cell.2011.08.023.View ArticlePubMedGoogle Scholar
- Chokas AL, Trivedi CM, Lu MM, Tucker PW, Li S, Epstein JA, Morrisey EE: Foxp1/2/4-NuRD interactions regulate gene expression and epithelial injury response in the lung via regulation of interleukin-6. J Biol Chem. 2010, 285: 13304-13313. 10.1074/jbc.M109.088468.View ArticlePubMedPubMed CentralGoogle Scholar
- Shi X, Wallis AM, Gerard RD, Voelker KA, Grange RW, DePinho RA, Garry MG, Garry DJ: Foxk1 promotes cell proliferation and represses myogenic differentiation by regulating Foxo4 and Mef2. J Cell Sci. 2012, 125: 5329-5337. 10.1242/jcs.105239.View ArticlePubMedPubMed CentralGoogle Scholar
- Noce T, Fujiwara Y, Sezaki M, Fujimoto H, Higashinakagawa T: Expression of a mouse zinc finger protein gene in both spermatocytes and oocytes during meiosis. Dev Biol. 1992, 153: 356-367. 10.1016/0012-1606(92)90120-6.View ArticlePubMedGoogle Scholar
- Chowdhury K, Goulding M, Walther C, Imai K, Fickenscher H: The ubiquitous transactivator Zfp-38 is upregulated during spermatogenesis with differential transcription. Mech Dev. 1992, 39: 129-142. 10.1016/0925-4773(92)90040-Q.View ArticlePubMedGoogle Scholar
- Watanabe S, Ichimura T, Fujita N, Tsuruzoe S, Ohki I, Shirakawa M, Kawasuji M, Nakao M: Methylated DNA-binding domain 1 and methylpurine-DNA glycosylase link transcriptional repression and DNA repair in chromatin. Proc Natl Acad Sci U S A. 2003, 100: 12859-12864. 10.1073/pnas.2131819100.View ArticlePubMedPubMed CentralGoogle Scholar
- Ogawa H, Ishiguro K, Gaubatz S, Livingston DM, Nakatani Y: A complex with chromatin modifiers that occupies E2F- and Myc-responsive genes in G0 cells. Science. 2002, 296: 1132-1136. 10.1126/science.1069861.View ArticlePubMedGoogle Scholar
- Trojer P, Cao AR, Gao Z, Li Y, Zhang J, Xu X, Li G, Losson R, Erdjument-Bromage H, Tempst P, Farnham PJ, Reinberg D: L3MBTL2 protein acts in concert with PcG protein-mediated monoubiquitination of H2A to establish a repressive chromatin structure. Mol Cell. 2011, 42: 438-450. 10.1016/j.molcel.2011.04.004.View ArticlePubMedPubMed CentralGoogle Scholar
- Ito S, Shen L, Dai Q, Wu SC, Collins LB, Swenberg JA, He C, Zhang Y: Tet proteins can convert 5-methylcytosine to 5-formylcytosine and 5-carboxylcytosine. Science. 2011, 333: 1300-1303. 10.1126/science.1210597.View ArticlePubMedPubMed CentralGoogle Scholar
- Webster J, Oxley D: Peptide mass fingerprinting: protein identification using MALDI-TOF mass spectrometry. Methods Mol Biol. 2005, 310: 227-240. 10.1007/978-1-59259-948-6_16.View ArticlePubMedGoogle Scholar
- The ProteomeXchange consortium. http://proteomecentral.proteomexchange.org,
- Vizcaino JA, Cote RG, Csordas A, Dianes JA, Fabregat A, Foster JM, Griss J, Alpi E, Birim M, Contell J, O’Kelly G, Schoenegger A, Ovelleiro D, Perez-Riverol Y, Reisinger F, Rios D, Wang R, Hermjakob H: The PRoteomics IDEntifications (PRIDE) database and associated tools: status in 2013. Nucleic Acids Res. 2013, 41: D1063-D1069. 10.1093/nar/gks1262.View ArticlePubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.View ArticlePubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37: 1-13. 10.1093/nar/gkn923.View ArticlePubMedGoogle Scholar
- The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7. http://david.abcc.ncifcrf.gov/home.jsp,
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.