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.     

Quantitative DNA methylation analysis based on four-dye trace data from direct sequencing of PCR amplificates

MOTIVATION: Methylation of cytosines in DNA plays an important role in the regulation of gene expression, and the analysis of methylation patterns is fundamental for the understanding of cell differentiation, aging processes, diseases and cancer development. Such analysis has been limited, because technologies for detailed and efficient high-throughput studies have not been available. We have developed a novel quantitative methylation analysis algorithm and workflow based on direct DNA sequencing of PCR products from bisulfite-treated DNA with high-throughput sequencing machines. This technology is a prerequisite for success of the Human Epigenome Project, the first large genome-wide sequencing study for DNA methylation in many different tissues. Methylation in tissue samples which are compositions of different cells is a quantitative information represented by cytosine/thymine proportions after bisulfite conversion of unmethylated cytosines to uracil and PCR. Calculation of quantitative methylation information from base proportions represented by different dye signals in four-dye sequencing trace files needs a specific algorithm handling imbalanced and overscaled signals, incomplete conversion, quality problems and basecaller artifacts. RESULTS: The algorithm we developed has several key properties: it analyzes trace files from PCR products of bisulfite-treated DNA sequenced directly on ABI machines; it yields quantitative methylation measurements for individual cytosine positions after alignment with genomic reference sequences, signal normalization and estimation of effectiveness of bisulfite treatment; it works in a fully automated pipeline including data quality monitoring; it is efficient and avoids the usual cost of multiple sequencing runs on subclones to estimate DNA methylation. The power of our new algorithm is demonstrated with data from two test systems based on mixtures with known base compositions and defined methylation. In addition, the applicability is proven by identifying CpGs that are differentially methylated in real tissue samples.     

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.     

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.     

Profiling DNA methylation by melting analysis

The idea of modifying DNA with bisulfite has paved the way for a variety of polymerase chain reaction (PCR) methods for accurately mapping 5-methylcytosine at specific genes. Bisulfite selectively deaminates cytosine to uracil under conditions where 5-methylcytosine remains unreacted. Following conventional PCR amplification of bisulfite-treated DNA, original cytosines appear as thymine while 5-methylcytosines appear as cytosine. Because the relative thermostability of a DNA duplex increases with increasing content of G:C base pairs, PCR products originating from DNA templates with different contents of 5-methylcytosine differ in melting temperature, i.e., the temperature required to convert the double helix into random coils. We describe two methods that resolve differentially methylated DNA sequences on the basis of differences in melting temperature. The first method integrates PCR amplification of bisulfite-treated DNA and subsequent melting analysis by using a thermal cycler coupled with a fluorometer. By including in the reaction a PCR-compatible, fluorescent dye that specifically binds to double-stranded DNA, the melting properties of the PCR product can be examined directly in the PCR tube by continuous fluorescence monitoring during a temperature transition. The second method relies on resolution of alleles with different 5-methylcytosine contents by analysis of PCR products in a polyacrylamide gel containing a gradient of chemical denaturants. Optimal resolution of differences in melting temperature is achieved by a special design of PCR primers. Both methods allow resolution of "heterogeneous" methylation, i.e., the situation where the content and distribution of 5-methylcytosine in a target gene differ between different molecules in the same sample.     

Whole genome amplification of sodium bisulfite-treated DNA allows the accurate estimate of methylated cytosine density in limited DNA resources

Sodium bisulfite modification-based fine mapping of methylated cytosines represents the gold standard technique for DNA methylation studies. A major problem with this approach, however is that it results in considerable DNA degradation, and large quantities of genomic DNA material are needed if numerous genomic regions are to be profiled. In this study, we examined whether whole genome amplification (WGA) techniques can be applied to sodium bisulfite-treated DNA and whether WGA would bias DNA methylation results. Sodium bisulfite-treated DNA was amplified using a standard WGA method: optimized primer-extension preamplification (PEP) with degenerate primers. Following the PCR of bisulfite-treated DNA, the DNA methylation profiles of specific DNA fragments were assessed using three approaches: (i) direct sequencing of the overall product; (ii) the sequencing of cloned PCR products; and (iii) methylation-sensitive single nucleotide primer extension (MS-SNuPE)--and compared with those obtained from bisulfite-treated DNA not subjected to WGA. Our data indicates that the DNA methylation profiles obtained from WGA of sodium bisulfite-treated DNA are consistent with those obtained from non-WGA DNA. The average difference in methylation percentage calculated from the two sets of template using MS-SNuPE was 4%. If our results are replicated on other genomic loci, WGA may become a useful technique in DNA methylation studies.     

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.     

Whole Methylome Analysis by Ultra-Deep Sequencing Using Two-Base Encoding

Methylation, the addition of methyl groups to cytosine (C), plays an important role in the regulation of gene expression in both normal and dysfunctional cells. During bisulfite conversion and subsequent PCR amplification, unmethylated Cs are converted into thymine (T), while methylated Cs will not be converted. Sequencing of this bisulfite-treated DNA permits the detection of methylation at specific sites. Through the introduction of next-generation sequencing technologies (NGS) simultaneous analysis of methylation motifs in multiple regions provides the opportunity for hypothesis-free study of the entire methylome. Here we present a whole methylome sequencing study that compares two different bisulfite conversion methods (in solution versus in gel), utilizing the high throughput of the SOLiD™ System. Advantages and disadvantages of the two different bisulfite conversion methods for constructing sequencing libraries are discussed. Furthermore, the application of the SOLiD™ bisulfite sequencing to larger and more complex genomes is shown with preliminary in silico created bisulfite converted reads.     

DNA甲基化分析中重亚硫酸盐处理DNA转化效率的评估方法

为建立一种评估重亚硫酸盐处理DNA样本后胞嘧啶转化效率的有效方法,以两组不同的TaqMan qPCR检测梯度稀释的重亚硫酸盐处理和未处理的DNA标准品,建立转化与未转化的DNA Ct值以及对应的DNA拷贝数的标准曲线。使用相同的探针定量检测重亚硫酸盐处理后的DNA样本评估转化效率。结果显示该方法应用两组探针,根据相应的标准曲线,精确评估样本经重亚硫酸盐处理的转化效率。使用已知转化和未转化拷贝数的混合DNA作为模板,证实了该方法的可靠性。同时也对不同重亚硫酸盐试剂盒处理DNA的转化效率进行了评估,结果显示,该方法能够有效地评估DNA样品重亚硫酸盐的转化效率,为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.     

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.     

A streamlined method for analysing genome-wide DNA methylation patterns from low amounts of FFPE DNA

Background

Formalin fixed paraffin embedded (FFPE) tumor samples are a major source of DNA from patients in cancer research. However, FFPE is a challenging material to work with due to macromolecular fragmentation and nucleic acid crosslinking. FFPE tissue particularly possesses challenges for methylation analysis and for preparing sequencing-based libraries relying on bisulfite conversion. Successful bisulfite conversion is a key requirement for sequencing-based methylation analysis.

Methods

Here we describe a complete and streamlined workflow for preparing next generation sequencing libraries for methylation analysis from FFPE tissues. This includes, counting cells from FFPE blocks and extracting DNA from FFPE slides, testing bisulfite conversion efficiency with a polymerase chain reaction (PCR) based test, preparing reduced representation bisulfite sequencing libraries and massively parallel sequencing.

Results

The main features and advantages of this protocol are:

• An optimized method for extracting good quality DNA from FFPE tissues.
• An efficient bisulfite conversion and next generation sequencing library preparation protocol that uses 50 ng DNA from FFPE tissue.
• Incorporation of a PCR-based test to assess bisulfite conversion efficiency prior to sequencing.

Conclusions

We provide a complete workflow and an integrated protocol for performing DNA methylation analysis at the genome-scale and we believe this will facilitate clinical epigenetic research that involves the use of FFPE tissue.

     

High-Speed Conversion of Cytosine to Uracil in Bisulfite Genomic Sequencing Analysis of DNA Methylation

Bisulfite genomic sequencing is a widely used technique foranalyzing cytosine-methylation of DNA. By treating DNA withbisulfite, cytosine residues are deaminated to uracil, whileleaving 5-methylcytosine largely intact. Subsequent PCR andnucleotide sequence analysis permit unequivocal determinationof the methylation status at cytosine residues. A major caveatassociated with the currently practiced procedure is that ittakes 16–20 hr for completion of the conversion of cytosineto uracil. Here we report that a complete deamination of cytosineto uracil can be achieved in shorter periods by using a highlyconcentrated bisulfite solution at an elevated temperature.Time course experiments demonstrated that treating DNA with9 M bisulfite for 20 min at 90°C or 40 min at 70°C allcytosine residues in the DNA were converted to uracil. Underthese conditions, the majority of 5-methylcytosines remainedintact. When a high molecular weight DNA derived from a cellline (containing a number of genes whose methylation statuswas known) was treated with bisulfite under the above conditionsand amplified and sequenced, the results obtained were consistentwith those reported in the literature. Although some degradationof DNA occurred during this process, the amount of treated DNArequired for the amplification was nearly equal to that requiredfor the conventional bisulfite genomic sequencing procedure.The increased speed of DNA methylation analysis with this novelprocedure is expected to advance various aspects of DNA sciences.     

Conversion-specific detection of DNA methylation using real-time polymerase chain reaction (ConLight-MSP) to avoid false positives

Methylated cytosines appear as sequence variations following bisulfite treatment and polymerase chain reaction (PCR) amplification. By using methylation-specific PCR (MSP), it is possible to detect methylated sequences in a background of unmethylated DNA with a high level of sensitivity. MSP is frequently used to identify methylated alleles in carcinogenesis, and may be combined with the TaqMan real-time PCR system, which uses fluorescence-based detection of amplification products during the amplification phase of the PCR and increases the sensitivity of detection (MethyLight). Sequences that have been incompletely converted during the bisulfite treatment are frequently coamplified during MSP, resulting in an overestimation of DNA methylation. The presence of amplified sequences originating from partially unconverted material may be determined by sequencing or by restriction digests or Southern blots of MSPs. Alternately, we have developed a method where the PCR and conversion assay are combined within a single TaqMan reaction by using an additional fluorescent probe directed against unconverted DNA (ConLight-MSP). We recommend that MSP detection always should include a step to detect unconverted DNA to avoid overestimation of the frequency or level of methylated DNA in the sample.     

An optimized rapid bisulfite conversion method with high recovery of cell-free DNA

Background

Methylation analysis of cell-free DNA is a encouraging tool for tumor diagnosis, monitoring and prognosis. Sensitivity of methylation analysis is a very important matter due to the tiny amounts of cell-free DNA available in plasma. Most current methods of DNA methylation analysis are based on the difference of bisulfite-mediated deamination of cytosine between cytosine and 5-methylcytosine. However, the recovery of bisulfite-converted DNA based on current methods is very poor for the methylation analysis of cell-free DNA.

Results

We optimized a rapid method for the crucial steps of bisulfite conversion with high recovery of cell-free DNA. A rapid deamination step and alkaline desulfonation was combined with the purification of DNA on a silica column. The conversion efficiency and recovery of bisulfite-treated DNA was investigated by the droplet digital PCR. The optimization of the reaction results in complete cytosine conversion in 30 min at 70 °C and about 65% of recovery of bisulfite-treated cell-free DNA, which is higher than current methods.

Conclusions

The method allows high recovery from low levels of bisulfite-treated cell-free DNA, enhancing the analysis sensitivity of methylation detection from cell-free DNA.

     

Formamide as a denaturant for bisulfite conversion of genomic DNA: Bisulfite sequencing of the GSTPi and RARβ2 genes of 43 formalin-fixed paraffin-embedded prostate cancer specimens

Analysis of methylated DNA, which refers to 5-methycytosine (5mC) versus cytosine (C) at specific loci in genomic DNA (gDNA), has received increased attention in epigenomics, particularly in the area of cancer biomarkers. Many different methods for analysis of methylated DNA rely on initial reaction of gDNA with concentrated acidic sodium bisulfite to quantitatively convert C to uracil (U) via sulfonation of denatured, single-stranded gDNA under conditions where 5mC is resistant to analogous sulfonation leading to thymine (T). These methods typically employ polymerase chain reaction (PCR) amplification after bisulfite conversion, thereby leading to readily detectable amounts of amplicons where T and C are measured as surrogates for C and 5mC in the original unconverted gDNA. However, incomplete bisulfite conversion of C in gDNA has been reported to be a common source of error in analysis of methylated DNA. Incomplete conversion can be revealed during the course of bisulfite sequencing, which is the generally accepted “gold standard” for analysis of methylated DNA. Previous bisulfite sequencing investigations of conventional predenaturation of gDNA with NaOH followed by the use of bisulfite containing added urea to maintain denaturation and thus mitigate incomplete conversion of C have been reported to give conflicting results. The current study describes a new approach where conventional predenaturation of gDNA with NaOH is instead achieved with formamide and maintains denaturation during subsequent sample handling and sulfonation. This formamide-based method was applied to 46 formalin-fixed/paraffin-embedded (FFPE) biopsy tissue specimens from well-characterized patients with primary prostate cancer. These specimens were representative of difficult-to-analyze samples due to the chemically compromised nature of the gDNA, which was recovered by modifying the protocol for a commercially available total RNA/DNA extraction kit (RecoverALL). An additional novel aspect of this study was analysis of CpG-rich promoter regions of two prostate cancer-related genes: glutathione S-transferase pi (GSTPi) and retinoic acid receptor beta2 (RARβ2). High-quality bisulfite sequencing results were obtained for both genes in 43 of 46 (93%) specimens. Detection of methylated GSTPi and RARβ2 genes was significantly associated with primary prostate cancer as compared with the benign prostate (Fisher’s exact test, P < 0.001). The sensitivity and specificity of detection of methylated GSTPi and RARβ2 genes were 86% and 100% and 91% and 100%, respectively. Moreover, the presence of either methylated gene was detected in primary prostate cancer with sensitivity and specificity of 100% and 100%, respectively. The results demonstrated a high degree of reliability of formamide-based denaturation and bisulfite conversion that should extend, generally, to FFPE and other types of samples intended for any analytical method predicated on bisulfite conversion. This pilot study also demonstrated the efficacy of determining methylation of these two genes with high sensitivity and specificity in FFPE biopsy tissue specimens. Moreover, the results showed a highly significant association of methylated GSTPi and RARβ2 genes with primary prostate cancer. Finally, this improved procedure for determining these two methylated genes may allow the detection of prostate cancer cells in core biopsy specimens with insufficient numbers of cells and poor morphology.     

A Modified Protocol for Bisulfite Genomic Sequencing of Difficult Samples

The bisulfite genomic sequencing protocol is a widely used method for analyzing DNA methylation. It relies on the deamination of unmethylated cytosine residues to uracil; however, its high rates of DNA degradation and incomplete cytosine to uracil conversion often lead to failed experiments, uninformative results, and false positives. Here, we report the addition of a single-step multiple restriction enzyme digestion (MRED) designed to differentially digest polymerase chain reaction products amplified from unconverted DNA while leaving those of converted DNA intact. We show that for our model system, RARB2 P2 promoter, use of MRED increased informative sequencings ninefold, and MRED did not alter the clonal representation in one fully methylated cell line, H-596, treated or not with 5-azadeoxycytidine, a methylation inhibitor. We believe that this method may easily be adapted for analyzing other genes and provide guidelines for selecting the most appropriate MRED restriction enzymes.     

GPU-BSM: A GPU-Based Tool to Map Bisulfite-Treated Reads

Cytosine DNA methylation is an epigenetic mark implicated in several biological processes. Bisulfite treatment of DNA is acknowledged as the gold standard technique to study methylation. This technique introduces changes in the genomic DNA by converting cytosines to uracils while 5-methylcytosines remain nonreactive. During PCR amplification 5-methylcytosines are amplified as cytosine, whereas uracils and thymines as thymine. To detect the methylation levels, reads treated with the bisulfite must be aligned against a reference genome. Mapping these reads to a reference genome represents a significant computational challenge mainly due to the increased search space and the loss of information introduced by the treatment. To deal with this computational challenge we devised GPU-BSM, a tool based on modern Graphics Processing Units. Graphics Processing Units are hardware accelerators that are increasingly being used successfully to accelerate general-purpose scientific applications. GPU-BSM is a tool able to map bisulfite-treated reads from whole genome bisulfite sequencing and reduced representation bisulfite sequencing, and to estimate methylation levels, with the goal of detecting methylation. Due to the massive parallelization obtained by exploiting graphics cards, GPU-BSM aligns bisulfite-treated reads faster than other cutting-edge solutions, while outperforming most of them in terms of unique mapped reads.