A systematic evaluation of whole genome amplification of bisulfite-modified DNA
© Bundo et al.; licensee BioMed Central Ltd. 2012
Received: 3 September 2012
Accepted: 5 November 2012
Published: 22 November 2012
Skip to main content
© Bundo et al.; licensee BioMed Central Ltd. 2012
Received: 3 September 2012
Accepted: 5 November 2012
Published: 22 November 2012
Studying DNA methylation profiles in detail should be the first step in epigenetic research. Although sodium bisulfite modification of genomic DNA is the gold standard method for DNA methylation analysis, this method results in the loss of the majority of the DNA material. Whole genome amplification (WGA) of bisulfite-modified DNA is expected to provide a rich source of materials, but its validity has not been thoroughly evaluated. In this study, we evaluated the extent of biased amplification in the WGA of bisulfite-modified DNA and the reproducibility of independent WGA reactions. We performed the multiple displacement amplification-based WGA separately three times. Each experiment included two reactions using 10 or 50 ng of bisulfite-modified DNA as template. DNA methylation levels were compared between WGA products and original bisulfite-modified DNA at about 450,000 CpG sites.
Using a sufficient amount of bisulfite-modified DNA for WGA was critical for downstream application. The considerable deviations from original bisulfite-modified DNA were found in the middle range of DNA methylation levels. Distribution of hyper- and hypomethylation were equal, which suggested that the deviation at each CpG site occurred randomly. Averaging the data from independently amplified WGA products dramatically improved the overall quality.
WGA of bisulfite-modified DNA could be a valuable tool for epigenetic research, but careful experimental design and data interpretation are required.
In mammals, DNA methylation is mainly observed at the cytosine residues of CpG dinucleotides. The methyl group is transferred to the fifth position of cytosine by DNA methyltransferases. This modification plays important roles in the regulation of gene expression . In the promoter region, where a CpG-rich region known as the CpG island is often situated, DNA methylation is generally involved in gene silencing. In the intragenic regions, where most of the methylcytosine is enriched, DNA methylation is associated with highly expressed genes and alternative splicing, although its precise role remains unclear [2–5].
DNA methylation is involved in genomic imprinting, X chromosome inactivation, and tissue-specific gene expression. Alteration of the DNA methylation results in developmental deficits and diseases [1, 6]. Studies have shown that, in addition to cancer, various kinds of diseases, such as autoimmune diseases, diabetes, and neuropsychiatric diseases, are associated with altered DNA methylation [6–8]. Detailed qualitative and quantitative analyses of DNA methylation profiles should be the first step of epigenetic research in clinical medicine.
Sodium bisulfite modification of genomic DNA, which converts non-methyl cytosine to uracil, has been the gold standard method for DNA methylation analysis for decades. However, this method causes the degradation of genomic DNA, resulting in a loss of the majority of DNA material. Therefore, the amount of DNA required for such analyses is often on the order of micrograms.
Whole genome amplification (WGA) has been used to amplify genomic DNA for sequencing and genotyping analyses [9–11]. This method may also be valuable for epigenetic research, as WGA could provide a large amount of DNA. Because WGA products lose DNA methylation during amplification, bisulfite-modified DNA must be used as the template for WGA. One expected challenge is the unbiased amplification of bisulfite-modified genome DNA, which shows reduced genome complexity. In addition, unlike conventional genotyping, which is basically represented as categorical values (two homozygous alleles and one heterozygous allele), methylation levels are represented as continuous values ranging from 0 to 100%. It will be important, therefore, to develop a precise quantitative assay using WGA products, as well as a strategy for appropriate data interpretation.
Previous studies have reported on the validity and limitations of the application of WGA to bisulfite-modified genomic DNA [12–15], but only a handful of CpG sites were evaluated. Therefore, the extent of biased amplification in the WGA products and reproducibility of independent WGA reactions at the genome-wide level remain largely unclear, as does the effect of the quantity of starting materials. In this study, we systematically examined the validity and limitations of the WGA of bisulfite-modified genomic DNA.
Yield and overall quality of whole genome amplification (WGA) products of bisulfite-modified DNA
Numberof the detected CpGa
Average methylation signal
Average unmethylation signal
no amplification 1
484,857 (99.81 %)
no amplification 2
484,841 (99.81 %)
no amplification 3
484,831 (99.81 %)
WGA 50 ng 1
480,869 (98.99 %)
WGA 50 ng 2
480,159 (98.85 %)
WGA 50 ng 3
479,465 (98.70 %)
WGA 10 ng 1
465,730 (95.88 %)
WGA 10 ng 2
463,572 (95.43 %)
WGA 10 ng 3
467,036 (96.14 %)
We performed an Infinium HumanMethylation450 assay, which evaluates DNA methylation levels of about 450,000 CpG sites, using WGA products of bisulfite-modified DNA and original bisulfite-modified DNA. The total number of reliably detected CpG sites was significantly decreased after WGA, and depended on the quantity of bisulfite-modified DNA used for WGA reaction (Table 1). Average signal intensities of both methylation and unmethylation probes were also concordantly decreased in the WGA products (Table 1).
We calculated pairwise correlations of the methylation profile between all possible pairs (Figure 2c). Whereas correlations among original bisulfite-modified DNAs were very high (R = 0.991 to 0.992), Those correlations between original bisulfite-modified DNAs and WGA products decreased progressively according to the quantity of bisulfite-modified DNA, ranging from R = 0.923 to 0.935 for 50 ng and from R = 0.855 to 0.865 for 10 ng. Likewise, the average correlation among triplicates also decreased and was dependent upon the quantity of bisulfite-modified DNA (average R = 0.911 and 0.826 for 50 and 10 ng, respectively).
We examined the characteristics of MDA-based WGA of bisulfite-modified genomic DNA in detail. First, our analyses revealed that when performing WGA, using a sufficient amount of bisulfite-modified genome DNA is critical. Comparison between WGA products using 50 or 10 ng of bisulfite-modified DNA clearly showed more deviations in the 10 ng reactions. Second, although methylation levels were relatively conserved and showed little deviation in the hypomethylated (beta value <0.3) and hypermethylated (>0.7) regions, considerable deviations were found in the middle range of DNA methylation levels. These findings were not platform-dependent results, as the independent pyrosequencing analysis confirmed the extent of deviations (Figure 6). This result is consistent with previous examinations of several representative CpG sites [12–14]. Third, given the nature of random amplification deviations, averaging the multiple WGA products considerably reduces the deviation. Averaging three WGA products of 10 ng of bisulfite-modified DNA showed better results with regard to beta value difference in the middle range of methylation levels as compared with WGA product of 50 ng of bisulfite-modified DNA (Figure 7).
In terms of clinical settings, researchers must be cautious when using WGA-based products for the detection of subtle methylation differences, which are often required for case–control studies of common diseases. There are several requirements for such an application: (1) the expected methylation difference should be relatively large between groups; (2) the use of a sufficient amount of bisulfite-modified DNA for WGA; and (3) preparation of independent WGA replicates whenever possible. Practically, pooling the independent WGA reactions is expected to be effective to reduce the deviations. One alternative approach would be to treat the methylation level as a categorical variable. For example, when we categorized methylation levels of the reference sample as hypomethylation (beta value <0.3) and hypermethylation (beta value >0.7), 94% and 92% of the hypomethylated and 92% and 83% of the hypermethylated probes were correctly detected in WGA product of 50 and 10 ng of bisulfite-modified DNA, respectively (data not shown).
In this study, we employed the MDA-based WGA method. Alternative WGA would be the primer extension preamplification (PEP)-based method . As PEP involves PCR reaction by DNA polymerase, MDA is believed to produce more unbiased amplified products . However, a previous study reported that both methods provided comparable results when bisulfite-modified DNA was used as template . One major drawback of MDA-based method would be the insufficient amplification from severely degraded DNA template . As bisulfite modification causes DNA degradation, low amount of input DNA will be resulted in the failure of unbiased amplification. Therefore, improving the bisulfite modification method, which prevents high degradation of DNA template, would be one of useful steps for further reduction of input DNA.
WGA of bisulfite-modified DNA may serve as a valuable tool for epigenetics research in the clinical medicine, although careful experimental design and data interpretation will be required.
This study was approved by the ethics committee at the University of Tokyo Hospital. Genomic DNA was extracted from peripheral blood of an adult Japanese female by using a standard phenol-chloroform extraction. Quality and integrity was assessed by optical density (OD) measurement and gel electrophoresis. Quantity of genomic DNA was measured using a Qubit dsDNA BR assay kit (Life Technologies, Carlsbad, Califolnia, USA) with a Qubit fluorometer (Life Technologies).
Sodium bisulfite modification of genomic DNA was performed with an EpiTect Bisulfite kit (Qiagen, Venlo, Netherlands), according to the manufacturer’s instructions. The quantity of bisulfite-modified DNA was measured using Qubit ssDNA assay kit (Life Technologies) with a Qubit fluorometer.
WGA of bisulfite-modified DNA was performed using EpiTect Whole Bisulfitome kit (Qiagen), according to the manufacturer’s instructions. In brief, bisulfite-modified DNA was amplified with a reaction buffer containing phi29 DNA polymerase at 28 °C for 8 h. WGA included three independent experiments. In each experiment, the sample set contained 50 or 10 ng of bisulfite-modified DNA and DDW as the negative control. In performing WGA, UV-irradiation of all the equipment successfully suppressed nonspecific amplification from negative controls  (data not shown).
An Illumina Infinium HumanMethylation450 assay (Illumina) was performed according to the manufacturer’s protocol. The DNA methylation level was represented as a beta value, which was calculated as the ratio of fluorescent signal intensity of the methylated probe to those of total (methylated and unmethylated) probes by the GenomeStudio software (Illumina) with the default settings. All data are publicly available (GSE39565).
We used the detection P value, which estimates the confidence that the signals from the target CpG probe are above those from the negative control probes, and average signal intensities for methylation and unmethylation probes were used to compensate for the quality of WGA products. P values below 0.01 were considered to represent specific detection of the target CpG probes. Spearman’s correlation and principal component analysis were performed with R (ver. 2.13.0; http://www.R-project.org/). Hierarchical clustering analysis was performed by another multidimensional analysis package (amap) implemented in the R software package .
DNA methylation levels of the selected CpG sites were measured by PSQ 96MA (Qiagen), according to the manufacturer’s instructions. Briefly, bisulfite-PCR product using a biotin labeled primer was mixed with a binding buffer containing streptavidin-sepharose beads. The reaction mixture was placed onto a MultiScreen-HV, Clear Plate (Millipore, Billerica, Massachusetts, USA). After applying the vacuum, the beads were treated with a denaturation solution, and were suspended with an annealing buffer containing a sequencing primer. The mixture was transferred onto a PyroMark Q96 Plate Low (Qiagen). Sequencing reaction was performed with a PyroMark Gold Q96 Reagents Kit (Qiagen). The percentage of methylation was calculated using the allele quantification algorithm of the PyroMark Q96 ID software 188.8.131.52 (Qiagen). Primers were listed in Table S1 in Additional file 1. Detailed methods are available upon request.
Whole genome amplification
Multiple displacement amplification
Deionized distilled water.
This work was partly supported by Grant-in-Aid for Scientific Research on Innovative Areas (23118002; Adolescent Mind and Self-Regulation) from the Ministry of Education, Culture, Sports, Science and Technology of Japan (MEXT). A part of this study is the result of “Development of biomarker candidates for social behavior” carried out under the Strategic Research Program for Brain Sciences by MEXT.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.