Abstract
Recently, the cytosine modifications 5-hydroxymethylcytosine (5hmC) and 5-formylcytosine (5fC) were found to exist in the genomic deoxyribonucleic acid (DNA) of a wide range of mammalian cell types. It is now important to understand their role in normal biological function and disease. Here we introduce reduced bisulfite sequencing (redBS-Seq), a quantitative method to decode 5fC in DNA at single-base resolution, based on a selective chemical reduction of 5fC to 5hmC followed by bisulfite treatment. After extensive validation on synthetic and genomic DNA, we combined redBS-Seq and oxidative bisulfite sequencing (oxBS-Seq) to generate the first combined genomic map of 5-methylcytosine, 5hmC and 5fC in mouse embryonic stem cells. Our experiments revealed that in certain genomic locations 5fC is present at comparable levels to 5hmC and 5mC. The combination of these chemical methods can quantify and precisely map these three cytosine derivatives in the genome and will help provide insights into their function.