The Behaviour of 5-Hydroxymethylcytosine in Bisulfite Sequencing

Background

We recently showed that enzymes of the TET family convert 5-mC to 5-hydroxymethylcytosine (5-hmC) in DNA. 5-hmC is present at high levels in embryonic stem cells and Purkinje neurons. The methylation status of cytosines is typically assessed by reaction with sodium bisulfite followed by PCR amplification. Reaction with sodium bisulfite promotes cytosine deamination, whereas 5-methylcytosine (5-mC) reacts poorly with bisulfite and is resistant to deamination. Since 5-hmC reacts with bisulfite to yield cytosine 5-methylenesulfonate (CMS), we asked how DNA containing 5-hmC behaves in bisulfite sequencing.

Methodology/Principal Findings

We used synthetic oligonucleotides with different distributions of cytosine as templates for generation of DNAs containing C, 5-mC and 5-hmC. The resulting DNAs were subjected in parallel to bisulfite treatment, followed by exposure to conditions promoting cytosine deamination. The extent of conversion of 5-hmC to CMS was estimated to be 99.7%. Sequencing of PCR products showed that neither 5-mC nor 5-hmC undergo C-to-T transitions after bisulfite treatment, confirming that these two modified cytosine species are indistinguishable by the bisulfite technique. DNA in which CMS constituted a large fraction of all bases (28/201) was much less efficiently amplified than DNA in which those bases were 5-mC or uracil (the latter produced by cytosine deamination). Using a series of primer extension experiments, we traced the inefficient amplification of CMS-containing DNA to stalling of Taq polymerase at sites of CMS modification, especially when two CMS bases were either adjacent to one another or separated by 1–2 nucleotides.

Conclusions

We have confirmed that the widely used bisulfite sequencing technique does not distinguish between 5-mC and 5-hmC. Moreover, we show that CMS, the product of bisulfite conversion of 5-hmC, tends to stall DNA polymerases during PCR, suggesting that densely hydroxymethylated regions of DNA may be underrepresented in quantitative methylation analyses.     

Reaction of cytidine with semicarbazide in the presence of bisulfite. A rapid modification specific for single-stranded polynucleotide.

Semicarbazide reacted rapidly with 5,6-dihydrocytidine-6-sulfonate, which was formed from cytidine by addition of bisulfite across the 5,6-double bond. The transaminated product, 5,6-dihydro-4-semicarbazido-2-ketotopyrimidine-6-sulfonate ribofuranoside, was identified by comparison with that formed by treatment of 4-semicarbazido-2-ketopyrimidine ribofuranoside with bisulfite. The progress of the transamination was monitored spectrophotometrically by use of a strong absorbance of the product in alkali. The reaction between cytidine and the semicarbazide-bisulfite mixture was optimal at pH 4.5. Complete transformation of cytidine into the product required only 5 min with the use of 3M semicarbazide-1M sodium bisulfite, pH 5.0, at the reaction temperature 37 degrees C. The product was stable in unbuffered solution but in phosphate buffers it underwent elimination of bisulfite to give 4-semicarbazido-2-ketopyrimidine ribofuranoside. The rate of the elimination at pH 7.0 and 37 degrees C increased proportionally with the increase of the phosphate concentration. Complete elimination was obtained by treatment with 1 M sodium phosphate for 2 h. When heat-denatured calf-thymus DNA was treated with 3 M semicarbazide-1 M bisulfite at 37 degrees C and pH 5.0 the transamination of reactive cytosine residues was completed by 10 min of incubation. At 20 degrees C, it required 85 min of incubation. Cytosine residues in native DNA did not react at all even by prolonged incubations. The modified DNA samples were further treated with a phosphate buffer at pH 7, producing 4-semicarbazido-2-ketopyrimidine residues in the DNA. Analysis of the base compositions of these samples by perchloric acid hydrolysis showed that the modification was selective to cytosine, which had been expected from studies with monomers. It also showed that the reactive cytosine residues in the denatured DNA, constitute about 80% of the total cytosine, which was consistent with the view that heat-denatured DNA still contains a considerable amount of secondary structure. The semicarbazide-bisulfite modification is expected to be a sensitive method to locate cytosine residues in single-stranded regions of polynucleotides.     

Cross-linking between 16S ribosomal RNA and protein S4 in Escherichia coli ribosomal 30S subunits effected by treatment with bisulfite/hydrazine and bromopyruvate

Cytosine in nucleic acids can be modified by treatment with a mixture of bisulfite and hydrazine. The reaction is specific for single-stranded regions of nucleic acids and the product is N4-aminocytosine. Bromopyruvate has been used for alkylation of protein SH groups and through its 2-oxo group it can form a hydrazone with N4-aminocytosine. Escherichia coli ribosomal 30S subunits were treated with 1 M sodium bisulfite + 2 M hydrazine in the presence of 10 mM MgCl2 at pH 7.0 and 37 degrees C for 30 min. By this treatment, 2.4 cytosine residues/molecule 16S rRNA were derivatized into N4-aminocytosines. 35S-labeled 30S subunits were modified in this way and then treated with 10 mM bromopyruvate at pH 8.0 and 37 degrees C for 5 min. Analysis in sodium dodecyl sulfate/sucrose density gradient centrifugation showed co-sedimentation of a part of the 35S radioactivity with the RNA. The co-sedimentation was dependent on both the bisulfite/hydrazine and the bromopyruvate treatments. The RNA-protein complex was prepared from unlabeled 30S subunits. The protein portion was labeled with 125I, the RNA portion was digested with nucleases, and then the hydrazone linkage between the protein and oligonucleotides was cleaved by treatment with 0.2 M HCl. The oligonucleotides formed were removed by dialysis and the protein was identified as S4 by two-dimensional electrophoresis and by sodium dodecyl sulfate/polyacrylamide gel electrophoresis. The results indicate that the cysteinyl residue of protein S4 at position 31 from the N-terminus is located close to a cytosine residue which is non-base-paired and easily accessible by the externally present bisulfite/hydrazine reagent.     

Rapid quantitation of methylation differences at specific sites using methylation-sensitive single nucleotide primer extension (Ms-SNuPE).

We have developed a rapid quantitative method (Ms-SNuPE) for assessing methylation differences at specific CpG sites based on bisulfite treatment of DNA followed by single nucleotide primer extension. Genomic DNA was first reacted with sodium bisulfite to convert unmethylated cytosine to uracil while leaving 5-methylcytosine unchanged. Amplification of the desired target sequence was then performed using PCR primers specific for bisulfite-converted DNA and the resulting product isolated and used as a template for methylation analysis at the CpG site(s) of interest. This methylation-sensitive technique has several advantages over existing methods used for detection of methylation changes because small amounts of DNA can be analyzed including microdissected pathology sections and it avoids utilization of restriction enzymes for determining the methylation status at CpG sites.     

Identification and resolution of artifacts in bisulfite sequencing

Bisulfite sequencing has become the most widely used application to detect 5-methylcytosine (5-MeC) in DNA, and provides a reliable way of detecting any methylated cytosine at single-molecule resolution in any sequence context. The process of bisulfite treatment exploits the different sensitivity of cytosine and 5-MeC to deamination by bisulfite under acidic conditions, in which cytosine undergoes conversion to uracil while 5-MeC remains unreactive. In this article, we address the more commonly encountered experimental artifacts associated with bisulfite sequencing, and provide methods for the detection and elimination of these artifacts. In particular, we focus on conditions that inhibit complete bisulfite-mediated conversion of cytosines in a target sequence, and demonstrate the necessity of complete protein removal from DNA samples prior to bisulfite treatment. We also include a brief summary of the experimental protocol for bisulfite treatment and tips for designing polymerase chain reaction (PCR) primers to amplify from bisulfite-treated DNA.     

Single-base resolution analysis of DNA epigenome via high-throughput sequencing

Epigenetic changes caused by DNA methylation and histone modifications play important roles in the regulation of various cellular processes and development. Recent discoveries of 5-methylcytosine (5mC) oxidation derivatives including 5-hydroxymethylcytosine (5hmC), 5-formylcytsine (5fC) and 5-carboxycytosine (5caC) in mammalian genome further expand our understanding of the epigenetic regulation. Analysis of DNA modification patterns relies increasingly on sequencing-based profiling methods. A number of different approaches have been established to map the DNA epigenomes with single-base resolution, as represented by the bisulfite-based methods, such as classical bisulfite sequencing (BS-seq), TAB-seq (TET-assisted bisulfite sequencing), oxBS-seq (oxidative bisulfite sequencing) and etc. These methods have been used to generate base-resolution maps of 5mC and its oxidation derivatives in genomic samples. The focus of this review will be to discuss the chemical methodologies that have been developed to detect the cytosine derivatives in the genomic DNA.     

The bisulfite genomic sequencing used in the analysis of epigenetic states, a technique in the emerging environmental genotoxicology research

Methylation at position 5 of cytosine in DNA is an important event in epigenetic changes of cells, the methylation being linked to the control of gene functions. The DNA methylation can be analyzed by bisulfite genomic sequencing, and a large body of data have now been accumulated, based on which causation of diseases, for example cancer, and many other manifestations of cellular activities have been discussed intensively. This article gives an extensive account of the chemical aspects of bisulfite modification of cytosine and 5-methylcytosine in DNA. Various factors that affect the action of bisulfite are discussed, and a recent progress from our laboratory is explained. Conventional procedures for the bisulfite treatment consist incubation of single-stranded DNA with sodium bisulfite under acidic conditions. This treatment converts cytosine into uracil, but 5-methylcytosine remains unchanged. Amplification by polymerase chain reaction (PCR) of the bisulfite-treated DNA followed by sequencing can result in revealing the positions of 5-methylcytosine in the gene. We have discovered that the whole procedure can be significantly speeded up by the use of a highly concentrated bisulfite solution, 10 M ammonium bisulfite. Another recent finding is that urea, which has been often added to the reaction mixture with the purpose of facilitating the bisulfite-mediated deamination of cytosine in DNA, may not work as anticipated: we have observed that urea does not show such promoting actions in our treatments of DNA. A laboratory protocol for quantifying bisulfite, suitably simple for routine practice to ensure valid experiments, is described.     

Purification and characterization of 5-hydroxymethyluracil-DNA glycosylase from calf thymus. Its possible role in the maintenance of methylated cytosine residues

5-Hydroxymethyluracil (HmUra) residues formed by the oxidation of thymine are removed from DNA through the action of a DNA glycosylase activity. This activity was purified over 1870-fold from calf thymus and found to be distinct from uracil (Ura)-DNA glycosylase. The HmUra-DNA glycosylase has a molecular weight of 38,000, a pH optimum of 6.7-6.8 and an apparent Km of 0.73 +/- 0.04 microM. These values are similar to those reported for other mammalian DNA glycosylases. The enzyme removed HmUra residues from single- and double-stranded DNA with almost equal efficiency. HmUra-DNA glycosylase activity was not product inhibited by free HmUra. The DNA glycosylase activity was inhibited by Mg2+, but the purest enzyme fractions contained a Mg2+-dependent apurinic/apyrimidinic endonuclease activity. HmUra-DNA glycosylase and the recently described 5-hydroxymethylcytosine (HmCyt)-DNA glycosylase (Cannon, S. V., Cummings, A. C., and Teebor, G. W. (1988) Biochem. Biophys. Res. Commun. 151, 1173-1179) are unique among known DNA glycosylases in being present in mammalian cells and absent from bacteria. These DNA glycosylase activities were shown here to reside on different proteins. We suggest that the major function of HmUra-DNA glycosylase, together with HmCyt-DNA glycosylase, is the maintenance of methylated cytosine residues in the DNA of higher organisms.     

Recognition of a cytosine base lesion by a human damage-specific DNA binding protein.

Sodium bisulfite reacts with cytosine and 5-methylcytosine, forming the 5,6-dihydrosulfonate adducts which deaminate to the uracil and thymine adducts, respectively. At alkaline pH, the sulfonate groups are then released, generating uracil and thymine. In DNA, the resulting G:U and G:T base mismatches generated are potential sites of mutagenesis. Using a human damage-specific DNA binding protein as a probe, we have found protein-recognizable lesions in bisulfite-treated DNA and poly d(I-C), but not in treated poly d(A-T) or poly d(A-U). Although this suggests that the lesion recognized is cytosine-derived, there was no correlation between the number of uracils induced and the number of binding sites, suggesting that the protein-bound damage is not a uracil-containing mismatch. Modification of the treatment protocol to reduce elimination of the bisulfite from the base adducts increased the level of binding, suggesting that the protein recognizes a base-sulfonate adduct.     

A new method for accurate assessment of DNA quality after bisulfite treatment

The covalent addition of methylgroups to cytosine has become the most intensively researched epigenetic DNA marker. The vast majority of technologies used for DNA methylation analysis rely on a chemical reaction, the so-called ‘bisulfite treatment’, which introduces methylation-dependent sequence changes through selective chemical conversion of non-methylated cytosine to uracil. After treatment, all non-methylated cytosine bases are converted to uracil but all methylated cytosine bases remain cytosine. These methylation dependent C-to-T changes can subsequently be studied using conventional DNA analysis technologies.

The bisulfite conversion protocol is susceptible to processing errors, and small deviation from the protocol can result in failure of the treatment. Several attempts have been made to simplify the procedure and increase its robustness. Although significant achievements in this area have been made, bisulfite treatment remains the main source of process variability in the analysis of DNA methylation. This variability in particular impairs assays, which strive for the quantitative assessment of DNA methylation. Here we present basic mathematical considerations, which should be taken into account when analyzing DNA methylation. We also introduce a PCR-based assay, which allows ab initio assessment of the DNA quality after bisulfite treatment and can help to prevent inaccurate quantitative measurement resulting from poor bisulfite treatment.

     

Control of carry-over contamination for PCR-based DNA methylation quantification using bisulfite treated DNA

In this study, we adapted the well known uracil DNA glycosylase (UNG) carry-over prevention system for PCR, and applied it to the analysis of DNA methylation based on sodium bisulfite conversion. As sodium bisulfite treatment converts unmethylated cytosine bases into uracil residues, bisulfite treated DNA is sensitive to UNG treatment. Therefore, UNG cannot be used for carry-over prevention of PCR using bisulfite treated template DNA, as not only contaminating products of previous PCR, but also the actual template will be degraded. We modified the bisulfite treatment procedure and generated DNA containing sulfonated uracil residues. Surprisingly, and in contrast to uracil, 6-sulfonyl uracil containing DNA (SafeBis DNA) is resistant to UNG. We showed that the new procedure removes up to 10 000 copies of contaminating PCR product in a closed PCR vessel without significant loss of analytical or clinical sensitivity of the DNA methylation analysis.     

Simple and accurate single base resolution analysis of 5-hydroxymethylcytosine by catalytic oxidative bisulfite sequencing using micelle incarcerated oxidants

Oxidation of 5-methylcytosine (5mC) is catalyzed by ten-eleven translocation (TET) enzymes to produce 5-hydroxymethylcytosine (5hmC) and following oxidative products. The oxidized nucleotides were shown to be the intermediates for DNA demethylation, as the nucleotides are removed by base excision repair system initiated by thymine DNA glycosylase. A simple and accurate method to determine initial oxidation product 5hmC at single base resolution in genomic DNA is necessary to understand demethylation mechanism. Recently, we have developed a new catalytic oxidation reaction using micelle-incarcerated oxidants to oxidize 5hmC to form 5-formylcytosine (5fC), and subsequent bisulfite sequencing can determine the positions of 5hmC in DNA. In the present study, we described the optimization of the catalytic oxidative bisulfite sequencing (coBS-seq), and its application to the analysis of 5hmC in genomic DNA at single base resolution in a quantitative manner. As the oxidation step showed quite low damage on genomic DNA, the method allows us to down scale the sample to be analyzed.     

Detection of Cytosine Methylation in Ancient DNA from Five Native American Populations Using Bisulfite Sequencing

While cytosine methylation has been widely studied in extant populations, relatively few studies have analyzed methylation in ancient DNA. Most existing studies of epigenetic marks in ancient DNA have inferred patterns of methylation in highly degraded samples using post-mortem damage to cytosines as a proxy for cytosine methylation levels. However, this approach limits the inference of methylation compared with direct bisulfite sequencing, the current gold standard for analyzing cytosine methylation at single nucleotide resolution. In this study, we used direct bisulfite sequencing to assess cytosine methylation in ancient DNA from the skeletal remains of 30 Native Americans ranging in age from approximately 230 to 4500 years before present. Unmethylated cytosines were converted to uracils by treatment with sodium bisulfite, bisulfite products of a CpG-rich retrotransposon were pyrosequenced, and C-to-T ratios were quantified for a single CpG position. We found that cytosine methylation is readily recoverable from most samples, given adequate preservation of endogenous nuclear DNA. In addition, our results indicate that the precision of cytosine methylation estimates is inversely correlated with aDNA preservation, such that samples of low DNA concentration show higher variability in measures of percent methylation than samples of high DNA concentration. In particular, samples in this study with a DNA concentration above 0.015 ng/μL generated the most consistent measures of cytosine methylation. This study presents evidence of cytosine methylation in a large collection of ancient human remains, and indicates that it is possible to analyze epigenetic patterns in ancient populations using direct bisulfite sequencing approaches.     

Microsoft Word macro for analysis of cytosine methylation by the bisulfite deamination reaction

Cytosine methylation at CpG dinucleotides is an important control mechanism in development, differentiation, and neoplasia. Bisulfite genomic sequencing and its modifications have been developed to examine methylation at these CpG dinucleotides. To use these methods, one has to (i) manually convert the sequence to that produced by bisulfite conversion and PCR amplification, taking into account that cytosine residues at CpG dinucleotides may or may not be converted depending on their methylation status, (ii) identify relevant restriction sites that may be used for methylation analysis, and (iii) conduct similar steps with the other DNA strand since the two strands of DNA are no longer complementary after bisulfite conversion. To automate these steps, we have developed a macro that can be used with Microsoft Word. This macro (i) converts genomic sequence to modified sequence that would result after bisulfite treatment facilitating primer design for bisulfite genomic sequencing and methylation-sensitive PCR assay and (ii) identifies restriction sites that are preserved in bisulfite-converted and PCR-amplified product only if cytosine residues at relevant CpG dinucleotides are methylated (and thereby not converted to uracil) in the genomic DNA.     

Role of uracil-DNA glycosylase in the repair of deaminated cytosine residues of DNA in Escherichia coli.

Uracil-DNA glycosylase, which acts specifically on uracil-containing DNA, was purified 250-fold from an extract of Escherichia coli 1100. The enzyme releases free uracil from DNA, producing alkali-labile apyrimidinic sites in the DNA. The enzyme is active on both native and heat-denatured DNA of phage PBS1, which contains uracil in place of thymine. piX174 DNA which had been treated with bisulfite and then at alkaline pH was susceptible to the action of uracil-DNA glycosylase. Since DNA treated with bisulfite alone was less susceptible to the enzyme, it is likely that the enzyme recognizes deaminated cytosine, namely uracil, but not bisulfite adducts of uracil and cytosine in the treated DNA. DNA treated with nitrite or hydroxylamine was not attacked by the enzyme. Enzyme activity acting on bisulfite-treated DNA was absent from an extract of E. coli mutant BD10 (ung). The mutant exhibited higher sensitivity to bisulfite than did the wild-type strain and was unable to reactivate phage T1 pre-exposed to bisulfite and weak alkali.     

Comparison of bisulfite modification of 5-methyldeoxycytidine and deoxycytidine residues.

Sodium bisulfite is a mutagen which can specifically deaminate more than 96% of the cytosine residues in single-stranded DNA via formation of a 5,6-dihydrocytosine-6-sulfonate intermediate. Under the same reaction conditions, only 2-3% of the 5-methylcytosine (m5Cyt) residues in single-stranded XP-12 DNA, which has 34 mole% m5Cyt, was converted to thymine (Thy) residues. In contrast, at the deoxynucleoside and free base levels, the same treatment with bisulfite and then alkali converted 51% and > 95%, respectively, of the m5Cyt to the corresponding Thy derivatives. However, the rate of reaction of m5Cyt and its deoxyribonucleoside was much slower than that of the analogous quantitative conversion of cytosine or deoxycytidine to uracil or deoxyuridine, respectively. The much lower reactivity of m5Cyt and its derivatives compared to that of the unmethylated analogs is primarily due to a decrease in the rate of formation of the sulfonate adduct.     

Bisulfite genomic sequencing: systematic investigation of critical experimental parameters

Bisulfite genomic sequencing is the method of choice for the generation of methylation maps with single-base resolution. The method is based on the selective deamination of cytosine to uracil by treatment with bisulfite and the sequencing of subsequently generated PCR products. In contrast to cytosine, 5-methylcytosine does not react with bisulfite and can therefore be distinguished. In order to investigate the potential for optimization of the method and to determine the critical experimental parameters, we determined the influence of incubation time and incubation temperature on the deamination efficiency and measured the degree of DNA degradation during the bisulfite treatment. We found that maximum conversion rates of cytosine occurred at 55°C (4–18 h) and 95°C (1 h). Under these conditions at least 84–96% of the DNA is degraded. To study the impact of primer selection, homologous DNA templates were constructed possessing cytosine-containing and cytosine-free primer binding sites, respectively. The recognition rates for cytosine (≥97%) and 5-methylcytosine (≥94%) were found to be identical for both templates.