High recovery of cell-free methylated DNA based on a rapid bisulfite-treatment protocol

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.


Background
Bisulfite induced modification of nucleic acids was originally described in the 1970s [1][2][3] and since the emergence of PCR and sequencing based detection methods [4,5] bisulfite treatment has played a pivotal role in the analysis of DNA methylation. In its original form it is a time consuming and labour intensive procedure involving numerous experimental steps: DNA denaturation, treatment with bisulfite for 12-16 hr, desalting and desulfonation with NaOH, and finally neutralisation and desalting. Recently published improvements include an accelerated deamination step, cutting down incubation time from 12-16 hr to 40 min, achieved by the use of a more concentrated bisulfite solution at higher temperatures [6,7]. The accelerated and the conventional methods have been explicitly compared by Genereux et al. [8]. The accelerated method leads to a more homogenous conversion of cytosine both across sites and molecules, conceivably due to facilitation of DNA denaturation in concentrated bisulfite solution at 70°C. Hence, inappropriate conversion of 5-methylcytosine, a result of prolonged bisulfite exposure of molecules with complete conversion of cytosine, is more controllable.
In addition to deamination of cytosine, bisulfite also induces chain breakage of DNA [9]. The DNA degradation caused by bisulfite treatment results in DNA fragments of an average length of approximately 600 nucleotides [10]. Real-time PCR based methods rely on amplification of short DNA fragments of 60-150 nucleotides. Hence, the use of real-time PCR limits the direct influence of fragmentation on the detection step. However, fragmentation affects recovery of DNA after bisulfite treatment. If the starting material is < 200 ng DNA, more than 95% of bisulfite-treated DNA is lost during desulfonation and purification with standard procedures [10,11]. This is a serious problem, especially when analysing material with very small amounts of DNA available, such as plasma with a median level of 10 ng cellfree DNA/ml in normal controls [11]. Although, the level of cell-free DNA is slightly increased in cancer patients [12], the combination of minute amounts of cell-free DNA in plasma and poor recovery after bisulfite treatment may lead to stochastic sampling issues. Improved recovery can be achieved by incorporation of DNA into agarose prior to bisulfite treatment. Denaturation, deamination and desulfonation are subsequently carried out on DNA embedded in the agarose beads [13]. This method has successfully been used to analyse DNA from microdissected cells [14]. Embedding of DNA in agarose is, however, a labour intensive method not amneable to automation, limiting its suitability in a clinical setting.
The ability to easily analyse sparse amounts of methylated DNA is a prerequisite for the usefulness of cellfree methylated DNA in plasma as a diagnostic or prognostic marker for cancer. An enormous amount of work has been put into identification of putative methylated DNA biomarkers and optimisation of the final analytical step: the detection of bisulfite-treated DNA. Hence, several reliable methods enabling the detection of minute amounts of bisulfite-treated DNA have been published [15][16][17][18][19][20]. Most recently the qMAMBA protocol (quantitative Methylation Analysis of Minute DNA amounts after whole Bisulfitome Amplification) has been developed, elegantly addressing several of the issues related to methylation analysis of samples with very low amounts of starting material [19,20]. However, as stated by Paliwal et al. [19] the most critical determinant of successful application of qMAMBA is the quality and quantity of starting material, emphasising the importance of the initial steps of the analysis: DNA isolation and bisulfitetreatment. Very few studies have addressed the loss of analytical sensitivity associated with the bisulfite treatment itself. Here we present a fast and reliable method, optimised in order to achieve high recovery, for the detection of methylated DNA from biospecimens with sparse amounts of DNA.

Bisulfite treatment
Cell-free DNA was isolated from 1 ml EDTA plasma on the EasyMAG nucleic acid purification platform (Biomeriéux), using the recommended protocol for plasma.
Bisulfite treatment was based on previously published accelerated methods [6,7], with some modifications. The optimised protocol is as follows: 50 μl of 10 M (NH 4 ) HSO 3 -NaHSO 3 solution were added to 25 μl of purified DNA in PCR strips. The strips were placed in a PCR machine, heated for 10 min at 90°C and subsequently cooled to room temperature. The solution containing the DNA-bisulfite adducts were purified on the Easy-MAG nucleic acid purification platform (Biomeriéux) using the off-board procedure according to manufacture's instructions, except for changes made to lysis buffer, extraction buffers A and B, and elution buffer. To ensure high recovery, 2 ml EasyMAG lysis buffer were replaced by a mix of 0.5 ml ethanol and 0.5 ml H 2 O, whereas extraction buffers A and B (Biomeriéux) both were replaced by 33% ethanol in H 2 O. For this extraction 25 μl magnetic beads were used. Finally, DNA was eluted in 25 μl 10 mM KOH. The alkaline solution in combination with the heating occurring during elution leads to desulfonation of DNA-bisulfite adducts.

DNA quantitation
To investigate recovery of bisulfite-treated DNA, 3 different primersets, all sharing the same detection probe, were designed: MLH1UF and MLH1 R detect DNA regardless of deamination [11]. MLH1 DF and MLH1 R detect deaminated DNA, whereas MLH1 UDF and MLH1 UDR detect undeaminated DNA (Table 1 and Figure 1B).

Measurement of analytical sensitivity
The sensitivity of methylation detection was assessed by analysing a dilution series of Universal Methylated DNA (Chemicon) in purified plasma DNA. After deamination the recovered DNA was quantified in triplicate for methylated RASFF1A and unmethylated RASSF1A using primers and probes listed in Table 1.
The efficiency of the deamination of unmethylated cytosines and the recovery of both methylated and unmethylated DNA fragments following bisulfite-treatment were elucidated using the hemimethylated MEST promoter as a model. Plasma DNA was diluted to concentrations of 10, 5, 2.5 and 1.25 genome equivalents pr. 25 μl of methylated and unmethylated MEST, respectively. Twentyfour samples of each dilution were deaminated and used in methylation-specific real-time PCR with MEST UM primers and probe and MEST M primers and probe, respectively. Primer and probe sequences are listed in Table 1. The detection frequencies were calculated as the percentage of positive reactions out of the 24 replicates.

Measurement of recovery
Recovery of deaminated plasma DNA was investigated using DNA from 24 plasma samples. DNA from 1 ml of each sample was extracted and eluted in 75 ul. Twentyfive ul were deaminated and extracted, another 25 ul were mock-deaminated (no heating) and extracted, whereas the residual was used for quantitation using MLH1 primers described in Table 1.

Results and discussion
Extensive optimisation of a previously published bisulfite treatment method [6,7] has led to the development of a fast bisulfite treatment with high recovery and potential for semi-automation. Optimisation involved the purification procedure after deamination in order to facilitate recovery of the fragmented DNA resulting from bisulfite treatment. This part of the procedure has previously been identified as a critical step [10]. Optimisation of the purification method had a pronounced effect on recovery, especially replacement of lysis and extraction buffers with ethanol increased recovery. In addition, adjustments of the reaction time were made. During optimisation the dynamics of the reaction had to be measured in order to ensure optimal conversion of cytosine and limited conversion of methyl-cytosine.
In the discovery phase bisulfite treatment, cloning and sequencing of individual clones is the method of choice, because of the detailed information of methylation status of individual CpG sites achieved. This procedure, from primer design to presentation of results has been described in detail [22]. When sequencing individual clones, complete deamination of unmethylated cytosine is extremely important. Incomplete deamination of unmethylated cytosine could lead to a false positive result. Typically more than one methylated CpG site is investigated when using methylated DNA as a tumour marker. Hence, this set up is less likely to yield a false positive result due to incomplete conversion of cytosine. However, if the reaction is poorly optimised and conversion of cytosine is far from complete, real-time detection will fail, leading to apparently low recovery. Therefore, while optimising recovery, dynamics of the reaction have been extensively monitored using both a real-time PCR based and a HPLC based method. The results of optimisation of the reaction time are shown in Figures 1  and 2. After 5 min of deamination, no undeaminated Table 1 Primer and beacon sequences   Primer/beacon  Sequence   MLH1 UF  TGT GAI AAA AAA TGT GAA GGG   MLH1 DF  GAA GAT ATT AGA TTT TAT GGG TTA TTT   MLH1 R  CAA CTI AAT TTT AAC AAA ATA ATC T   MLH1 UDF  ACC AGA TTT TAT GGG TCA TCC   MLH1 UDR  TTC TAT TAA    time PCR, where both primers and probes generally are specific for either the methylated or the unmethylated sequence. However, no cross-reactivity is observed in our experiment. The primers used to discriminate between undeaminated and deaminated product cover 5 cytosines (Figure 1). Primer binding mainly depends on the conversion of the 2 cytosines in the 3' end of the primer binding site. Therefore complete conversion of cytosine appears to be achieved faster when measured by real-time PCR compared with HPLC, where a limited amount of cytosines (< 1%) is detectable after 10 min deamination, and none after 15 min (Figure 2). Longer deamination time than 10 min results in reduced detection of both deaminated and undeaminated DNA, since the purpose is to optimise this protocol in order to achieve the best possible recovery of bisulfite treated DNA, the deamination time was set to 10 min.
In addition to incomplete conversion of cytosine, inappropriate conversion of methylated cytosine could lead to reduced sensitivity when detecting methylated DNA. The real-time PCR based experiment ( Figure 1) does not address the problem of over-conversion, since the primers and the probe are designed to an unmethylated part of MLH1. Methylated cytosine is monitored by the HPLC procedure ( Figure 2). However, methylated cytosine constitutes only a small fraction of the nucleotides. Therefore, to further elucidate the reaction dynamics and recovery of both methylated and unmethylated DNA, the hemimethylated MEST promoter has been used as a naturally occurring model system. Comparable detection of methylated and un-methylated MEST was observed, when calculating detection frequency as the percentage of positive reactions out of the 24 replicates ( Figure 3).
Similar experiments monitoring the dynamics of the reaction were carried out when changing other parameters, such as reaction temperature and denaturation of DNA prior to bisulfite treatment (data not shown). Denaturation of DNA did not affect recovery or conversion in the current study, even though this step has been pointed out to be critical, and different denaturants such as NaOH, urea and formamide have been used. Formamide was shown to have a great effect when working on DNA from formalin-fixed paraffinembedded tissue [23]. The lack of need of denaturants in the protocol presented here possibly reflects the denaturating effect of the high concentration of bisulfite and high temperature used for deamination.
Further investigation of the analytical sensitivity of the newly optimised procedure was carried out using a dilution of universally methylated DNA in purified plasma DNA. RASSF1A was used as a model system. RASSF1A is fully methylated in the universally methylated DNA and unmethylated in the plasma DNA sample. Table 2 illustrates that low copy numbers of methylated RASSF1A promotor region was reliably detected in a background of approximately 2000 copies of unmethylated DNA.
Recovery of deaminated DNA using the optimised method was compared with recovery of DNA in a matrix of deamination reagents. Re-extraction of DNA in a matrix of deamination reagents results in a recovery of 83% whereas recovery of deaminated DNA is 60% ( Table 3). The likely explanation of the difference is the unavoidable degradation of DNA during deamination.
Markers such as septin 9 hold great promise of clinical usefulness. The improved recovery rate of the bisulfite reaction achieved by the presented protocol may further improve its potential in a clinical setting, and avoid some of the expensive and labour intensive steps such as the need to purify multiple plasma samples in parallel and subsequent pooling of DNA [16,24]. Indeed, the initial promising reports on the potential of septin 9 have been strengthened by a recent report employing silica based magnetic DNA purification [25]. This work by deVos et al. address the problem of poor recovery, and the reported recovery rates are comparable to our results. However, the quantitation of DNA prior to and after bisulfite treatment is based on two different real-time PCR assays. We have employed a real-time PCR assay capable of detecting both treated and untreated DNA [11] allowing the same assay to be used for measuring DNA levels both pre and post bisulfite treatment resulting in improved accuracy

Conclusions
In this report we present a fast and reliable bisulfite treatment with high recovery from samples with minute amounts of DNA (< 0.1 ng/ml), opening the possibility of improving diagnostic sensitivity of the wide-range of potential markers already identified. The major improvement in recovery is achieved by alterations in the purification of bisulfite treated DNA. In addition, the method presented here allows preparation of deaminated DNA from patient samples in less than 2 hours. The development of this method amenable to semiautomation should considerably enhance the usefulness of methylation analysis in a clinical setting.