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.     

Reaction of bisulfite with the 5-hydroxymethyl group in pyrimidines and in phage DNAs.

5-Hydroxymethylcytosine reacted with bisulfite and, instead of undergoing usual deamination process, gave cytosine 5-methylenesulfonate as the product. The conversion was rapid and quantitative, and the optimum pH was 4.5. The product was isolated as crystals and characterized. Cytosine 5-methylenesulfonate was only very slowly deaminated by treatment with bisulfite. 5-Hydroxymethyl-2'-deoxycytidine 5'-phosphate reacted with bisulfite in the same way as 5-hydroxymethylcytosine. Residues of 5-hydroxymethylcytosine in native as well as denatured T2 DNA were convertible to those of cytosine 5-methylenesulfonate by treatment of the DNA with bisulfite. While it is known that the 5-hydroxy-methyl groups of T-even bacteriophage DNA can be enzymatically glucosylated, this observation offers chemical evidence that the 5-hydrozymethyl groups in DNA are situated in such a way that they can readily react with external agents. 5-Hydroxymethyluracil gave uracil 5-methylenesulfonate on treatment with bisulfite. This reaction was much slower than that of 5-hydroxymethylcytosine, and the optimum pH was between 6 and 7.     

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.

     

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.     

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.     

A polymerase engineered for bisulfite sequencing

Bisulfite sequencing is a key methodology in epigenetics. However, the standard workflow of bisulfite sequencing involves heat and strongly basic conditions to convert the intermediary product 5,6-dihydrouridine-6-sulfonate (dhU6S) (generated by reaction of bisulfite with deoxycytidine (dC)) to uracil (dU). These harsh conditions generally lead to sample loss and DNA damage while milder conditions may result in incomplete conversion of intermediates to uracil. Both can lead to poor recovery of bisulfite-treated DNA by the polymerase chain reaction (PCR) as either damaged DNA and/or intermediates of bisulfite treatment are poor substrate for standard DNA polymerases. Here we describe an engineered DNA polymerase (5D4) with an enhanced ability to replicate and PCR amplify bisulfite-treated DNA due to an ability to bypass both DNA lesions and bisulfite intermediates, allowing significantly milder conversion conditions and increased sensitivity in the PCR amplification of bisulfite-treated DNA. Incorporation of the 5D4 DNA polymerase into the bisulfite sequencing workflow thus promises significant sensitivity and efficiency gains.     

Bisulfite modification of immobilized DNAs for methylation detection

We have developed a novel method for detecting DNA methylation status of multiple samples, in which the DNA samples were firstly immobilized on the slide and treated with bisulfite directly on the chip. In this experiment, DNAs of pUC19 plasmid were restricted by the enzymes, and ligated with a linker bearing 5'-terminal acrylamide group at the sticky ends. Using universal acrylamide gel polymerization technique, a large amount of DNAs could be immobilized on the slide. The immobilized DNAs were converted by soaking the chip in bisulfite reaction mixtures for 16 h. The probes for detection of the methylation patterns of CpG sites hybridized with the converted DNAs on the microarray, and non-specifically bound probes were cleaned by electrophoresis. We have optimized the experimental conditions of both bisulfite treatment and electrophoresis to increase sensitivity and specificity. The results were further validated by bisulfite DNA sequencing. The experiments show that the method can simplify the experimental processes and increase the efficiency of the bisulfite treatment. This novel method could be used as a convenient tool to detect the methylation status of the multiple genes for a large amount of samples in the future.     

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.     

Tritium labelling of nucleic acids accompanied by conversion of cytosine to uracil.

A new method of incorporation of tritium into nucleic acids with an accompanying conversion of cytosine to uracil is proposed. The method is based on the reaction of nucleic acids with bisulfite in the presence of 3H2O. Under certain conditions poly(C) is quantitatively converted to a radioactive poly(U), whereas similar bisulfite treatment of poly(U) does not result in any tritium incorporation. Specificity of the reaction is confirmed by the results of analysis of modified tRNA and rRNA. Incubation of tRNA with bisulfite and 3H2O does not lead to cleavage of the polynucleotide chain. Similar treatment of the denatured DNA results in tritium incorporation into DNA which is accompanied by a conversion of cytosine to uracil. There is virtually no reaction between native DNA and bisulfite. Only certain cytosone residues in yeast tRNAVal/2a interact with bisulfate providing that reaction is carried out under sufficiently mild conditions.     

Archived Guthrie blood spots as a novel source for quantitative DNA methylation analysis

Sodium bisulfite treatment followed by PCR and DNA sequencing is widely considered the gold standard for the analysis of DNA methylation patterns. However, this technique generally requires substantial quantities of genomic DNA as starting material and is often associated with degradation of DNA. Here, we assess the feasibility of performing bisulfite sequencing on DNA isolated from 3-mm diameter punches of dried blood Guthrie spots. We demonstrate that it is possible to perform bisulfite sequencing from both freshly prepared and archived dried blood spots, using a combination of high purity DNA extraction and efficient bisulfite conversion. With the number of new technologies available for DNA methylation studies, we have extended this analysis and have successfully used a high-throughput mass spectrometry method for DNA methylation analysis on these samples. This provides a new source of material for epigenetic analysis of birth samples and provides an invaluable reference point to track temporal change in epigenetic profiles possibly linked with health and disease.     

Methylation sequencing from limiting DNA: embryonic,fixed, and microdissected cells

It is frequently useful to determine the methylation state of samples containing limited amounts of DNA such as from embryos, or from fixed tissue samples in which DNA is degraded or difficult to isolate. By modification of the standard protocols for DNA preparation and bisulfite treatment, it is possible to obtain DNA methylation sequence data for such samples. We present methods for bisulfite treatment of embryos, fixed sections, and samples obtained by laser capture microdissection, and discuss the additional experimental considerations required when working with small numbers of cells or degraded DNA samples.     

Enhanced Methylation Analysis by Recovery of Unsequenceable Fragments

Bisulfite sequencing is a valuable tool for mapping the position of 5-methylcytosine in the genome at single base resolution. However, the associated chemical treatment causes strand scission, which depletes the number of sequenceable DNA fragments in a library and thus necessitates PCR amplification. The AT-rich nature of the library generated from bisulfite treatment adversely affects this amplification, resulting in the introduction of major biases that can confound methylation analysis. Here, we report a method that enables more accurate methylation analysis, by rebuilding bisulfite-damaged components of a DNA library. This recovery after bisulfite treatment (ReBuilT) approach enables PCR-free bisulfite sequencing from low nanogram quantities of genomic DNA. We apply the ReBuilT method for the first whole methylome analysis of the highly AT-rich genome of Plasmodium berghei. Side-by-side comparison to a commercial protocol involving amplification demonstrates a substantial improvement in uniformity of coverage and reduction of sequence context bias. Our method will be widely applicable for quantitative methylation analysis, even for technically challenging genomes, and where limited sample DNA is available.     

Effects of bisulfite (sulfur dioxide) on DNA synthesis and fidelity by DNA polymerase I

In an attempt to explain the mechanism of comutagenesis by bisulfite, the extent and accuracy of DNA synthesis by E. coli DNA polymerase I was examined in the presence of sodium bisulfite. Bisulfite concentration of 100 mM caused nearly complete inhibition of dNTP incorporation into activated calf thymus DNA. Other salts (NaCL, Na2SO4) at the same concentration had no effect on enzyme activity. Preincubation of the various DNA synthesis assay components in 100 mM bisulfite showed that only preincubation of DNA polymerase I caused inhibition of DNA synthesis. Exonuclease functions of DNA polymerase I were unaffected by up to 100 mM bisulfite. Accuracy of DNA synthesis in the presence of bisulfite was determined using poly (dA-dT) as a template-primer. Concentrations of bisulfite greater than 50 mM caused a progressive decrease in enzyme accuracy. At 100 mM bisulfite there was an approximate 7.5-fold decrease in the fidelity of DNA synthesis, compared to control values, as measured by the ratio of noncomplementary (dGTP) to complementary (dTTP) nucleotide incorporated. Based on the known chemistry of bisulfite, it is hypothesized that sulfitolysis of the one disulfide group in DNA polymerase I by bisulfite might be responsible for the reduced polymerase activity and accuracy. The exonuclease functions of DNA polymerase I do not seem to require the disulfide linkage. These results suggest that the effects of bisulfite on mutation frequency might be mediated by effects on the fidelity of DNA repair systems.     

Anisotropic staining of apurinic acid with toluidine blue.

Using normal rat liver imprints, studies were carried out on the effects of histone extraction and the formation of aldehyde groups from deoxyribose on anisotropic toluidine blue staining of depurinized DNA after sodium bisulfite treatment. The anisotropic effect of bisulfite was found to be determined by binding of bisulfite ions to the aldehyde groups of apurinic acid which, together with free phosphate groups of DNA ensure coparallel attachment of the dye molecules. It was also shown that at pH 5.0 toluidine blue binds with both the phosphate and aldehyde groups of apurinic acid, to give anisotropic staining.     

Anisotropic staining of apurinic acid with toluidine blue

Summary Using normal rat liver imprints, studies were carried out on the effects of histone extraction and the formation of aldehyde groups from deoxyribose on anisotropic toluidine blue staining of depurinized DNA after sodium bisulfite treatment. The anisotropic effect of bisulfite was found to be determined by binding of bisulfite ions to the aldehyde groups of apurinic acid which, together with free phosphate groups of DNA ensure coparallel attachment of the dye molecules. It was also shown that at pH 5.0 toluidine blue binds with both the phosphate and aldehyde groups of apurinic acid, to give anisotropic staining.     

Directional telomeric silencing and lack of canonical B1 elements in two silencer Autonomously Replicating Sequences in S. cerevisiae

Background

Detection of cell-free methylated DNA in plasma is a promising tool for tumour diagnosis and monitoring. Due to the very low amounts of cell-free DNA in plasma, analytical sensitivity is of utmost importance. The vast majority of currently available methods for analysing DNA methylation are based on bisulfite-mediated deamination of cytosine. Cytosine is rapidly converted to uracil during bisulfite treatment, whereas 5-methylcytosine is only slowly converted. Hence, bisulfite treatment converts an epigenetic modification into a difference in sequence, amenable to analysis either by sequencing or PCR based methods. However, the recovery of bisulfite-converted DNA is very poor.

Results

Here we introduce an alternative method for the crucial steps of bisulfite treatment with high recovery. The method is based on an accelerated deamination step and alkaline desulfonation in combination with magnetic silica purification of DNA, allowing preparation of deaminated DNA from patient samples in less than 2 hours.

Conclusions

The method presented here allows low levels of DNA to be easily and reliably analysed, a prerequisite for the clinical usefulness of cell-free methylated DNA detection in plasma.     

Enzymatic degradation of uracil-containing deoxyribonucleic acid. V. Survival of Escherichia coli and coliphages treated with sodium bisulfite.

A number of mutants of Escherichia coli defective in the ung gene (structural gene for uracil-deoxyribonucleic acid [ura-DNA] glycosylase) are shown to be abnormally sensitive to treatment with sodium bisulfite when compared with congenic ung+ strains. These results provide further evidence that sodium bisulfite causes the deamination of cytosine to uracil in DNA and that ura-DNA glycosylase is required for the repair of U-G mispairs. The effect of the chemical is apparently selective with respect to base damage; coliphages containing cytosine in their DNA are inactivated by treatment with sodium bisulfite, whereas those containing hydroxymethylcytosine are not. ura-DNA glycosylase and the major apurinic-apyrimidinic endonuclease of E. coli may function in the same repair pathway, since the extent of inactivation of a congenic set of strains which are ung xth (structural gene for the major apurinic-apyrimidinic endonuclease of E. coli) or ung xth+ is the same.     

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.     

Editor's choice: Highly sensitive targeted methylome sequencing by post-bisulfite adaptor tagging

The current gold standard method for methylome analysis is whole-genome bisulfite sequencing (WGBS), but its cost is substantial, especially for the purpose of multi-sample comparison of large methylomes. Shotgun bisulfite sequencing of target-enriched DNA, or targeted methylome sequencing (TMS), can be a flexible, cost-effective alternative to WGBS. However, the current TMS protocol requires a considerable amount of input DNA and hence is hardly applicable to samples of limited quantity. Here we report a method to overcome this limitation by using post-bisulfite adaptor tagging (PBAT), in which adaptor tagging is conducted after bisulfite treatment to circumvent bisulfite-induced loss of intact sequencing templates, thereby enabling TMS of a 100-fold smaller amount of input DNA with far fewer cycles of polymerase chain reaction than in the current protocol. We thus expect that the PBAT-mediated TMS will serve as an invaluable method in epigenomics.     

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.