Cubic regression-based degree of correction predicts the performance of whole bisulfitome amplified DNA methylation analysis

Epigenetic mechanisms, including DNA methylation, are important determinants in development and disease. There is a need for technologies capable of detecting small variations in methylation levels in an accurate and reproducible manner, even if only limited amounts of DNA are available (which is the case in many studies in humans). Quantitative methylation analysis of minute DNA amounts after whole bisulfitome amplification (qMAMBA) has been proposed as an alternative, but this technique has not been adequately standardized and no comparative study against conventional methods has been performed, that includes a wide range of methylation percentages and different target assays. We designed an experiment to compare the performance of qMAMBA and bisulfite-treated genomic (non-amplified) DNA pyrosequencing. Reactions were performed in duplicate for each technique in eight different target genes, using nine artificially constructed DNA samples with methylation levels ranging between 0% and 100% with intervals of 12.5%. Cubic polynomial curves were plotted from the experimental results and the real methylation values and the resulting equation was used to estimate new corrected data points. The use of the cubic regression-based correction benefits the accuracy and the power of discrimination in methylation studies. Additionally, dispersion of the new estimated data around a y = x line (R²) served to fix a cutoff that can discriminate, with a single 9-point curve experiment, whether whole bisulfitome amplification and subsequent qMAMBA can produce accurate methylation results. Finally, even with an optimized reagent kit, DNA samples subjected to whole bisulfitome amplification enhance the preferential amplification of unmethylated alleles, and subtle changes in methylation levels cannot be detected confidently.     

Quantitative analysis of DNA methylation after whole bisulfitome amplification of a minute amount of DNA from body fluids

Cell-free circulating DNA isolated from the plasma of individuals with cancer has been shown to harbor cancer-associated changes in DNA methylation, and thus it represents an attractive target for biomarker discovery. However, the reliable detection of DNA methylation changes in body fluids has proven to be technically challenging. Here we describe a novel combination of methods that allows quantitative and sensitive detection of DNA methylation in minute amounts of DNA present in body fluids (quantitative Methylation Analysis of Minute DNA amounts after whole Bisulfitome Amplification, qMAMBA). This method involves genome-wide amplification of bisulphite-modified DNA template followed by quantitative methylation detection using pyrosequencing and allows analysis of multiple genes from a small amount of starting DNA. To validate our method we used qMAMBA assays for four genes and LINE1 repetitive sequences combined with plasma DNA samples as a model system. qMAMBA offered high efficacy in the analysis of methylation levels and patterns in plasma samples with extremely small amounts of DNA and low concentrations of methylated alleles. Therefore, qMAMBA will facilitate methylation studies aiming to discover epigenetic biomarkers, and should prove particularly valuable in profiling a large sample series of body fluids from molecular epidemiology studies as well as in tracking disease in early diagnostics.     

Assessment of RainDrop BS-seq as a method for large-scale,targeted bisulfite sequencing

We present a systematic assessment of RainDrop BS-seq, a novel method for large-scale, targeted bisulfite sequencing using microdroplet-based PCR amplification coupled with next-generation sequencing. We compared DNA methylation levels at 498 target loci (1001 PCR amplicons) in human whole blood, osteosarcoma cells and an archived tumor tissue sample. We assessed the ability of RainDrop BS-seq to accurately measure DNA methylation over a range of DNA quantities (from 10 to 1500 ng), both with and without whole-genome amplification (WGA) following bisulfite conversion. DNA methylation profiles generated using at least 100 ng correlated well (median R = 0.92) with those generated on Illumina Infinium HumanMethylation450 BeadChips, currently the platform of choice for epigenome-wide association studies (EWAS). WGA allowed for testing of samples with a starting DNA amount of 10 and 50 ng, although a reduced correlation was observed (median R = 0.79). We conclude that RainDrop BS-seq is suitable for measuring DNA methylation levels using nanogram quantities of DNA, and can be used to study candidate epigenetic biomarker loci in an accurate and high-throughput manner, paving the way for its application to routine clinical diagnostics.     

Assessment of RainDrop BS-seq as a method for large-scale,targeted bisulfite sequencing

We present a systematic assessment of RainDrop BS-seq, a novel method for large-scale, targeted bisulfite sequencing using microdroplet-based PCR amplification coupled with next-generation sequencing. We compared DNA methylation levels at 498 target loci (1001 PCR amplicons) in human whole blood, osteosarcoma cells and an archived tumor tissue sample. We assessed the ability of RainDrop BS-seq to accurately measure DNA methylation over a range of DNA quantities (from 10 to 1500 ng), both with and without whole-genome amplification (WGA) following bisulfite conversion. DNA methylation profiles generated using at least 100 ng correlated well (median R = 0.92) with those generated on Illumina Infinium HumanMethylation450 BeadChips, currently the platform of choice for epigenome-wide association studies (EWAS). WGA allowed for testing of samples with a starting DNA amount of 10 and 50 ng, although a reduced correlation was observed (median R = 0.79). We conclude that RainDrop BS-seq is suitable for measuring DNA methylation levels using nanogram quantities of DNA, and can be used to study candidate epigenetic biomarker loci in an accurate and high-throughput manner, paving the way for its application to routine clinical diagnostics.     

Synthesis of universal unmethylated control DNA by nested whole genome amplification with phi29 DNA polymerase

Optimization of highly sensitive methods to detect methylation of CpG islands in gene promoter regions requires adequate methylated and unmethylated control DNA. Whereas universal methylated control DNA is available, universal unmethylated control (UUC) DNA has not been made because demethylase is not available to remove methyl groups from all methylated cytosines. On the basis that DNA synthesized by DNA polymerase does not contain methylated cytosines, we developed a method to create UUC DNA by nested whole genome amplification (WGA) with phi29 DNA polymerase. Contamination of the template genomic DNA in UUC was only 3.1 x 10(-7), below the detection limit of sensitive methods used for methylation studies such as methylation-specific PCR. Assessment of microsatellite markers demonstrated that even nested phi29 WGA achieves highly accurate and homogeneous amplification with very low amounts of genomic DNA as an initial template. The UUC DNA created by nested phi29 WGA is practically very useful for methylation analysis.     

Quantitative measurement of genome-wide DNA methylation by a reliable and cost-efficient enzyme-linked immunosorbent assay technique

DNA methylation, the conversion of cytosine to 5-methylcytosine, is an important epigenetic modification involved in gene regulation. DNA methylation is essential for normal development whereas abnormal methylation has been implicated in pathological conditions including cancer. To evaluate the extent and variation of genome-wide DNA methylation and its changes during cellular differentiation and tumorgenesis as well as the interplay with histone modifications, accurate and reproducible quantification of the genomic DNA methylation level is required. These measurements have so far been achieved only by sophisticated and costly techniques. Here we report the generation of an enzyme-linked immunosorbent assay (methDNA-ELISA) for the accurate quantification of global DNA methylation levels. The linear region of this methDNA-ELISA ranges from 1 to 10%, making it highly suitable for the typical ranges from 2 to 6% in mammalian genomes. This method requires 10 ng of isolated DNA per sample, thus permitting investigation with minimal amounts of DNA previously not applicable for global DNA methylation analysis, e.g., clinical biopsies or cells collected by microdissection.     

High sensitivity mapping of methylated cytosines.

An understanding of DNA methylation and its potential role in gene control during development, aging and cancer has been hampered by a lack of sensitive methods which can resolve exact methylation patterns from only small quantities of DNA. We have now developed a genomic sequencing technique which is capable of detecting every methylated cytosine on both strands of any target sequence, using DNA isolated from fewer than 100 cells. In this method, sodium bisulphite is used to convert cytosine residues to uracil residues in single-stranded DNA, under conditions whereby 5-methylcytosine remains non-reactive. The converted DNA is amplified with specific primers and sequenced. All the cytosine residues remaining in the sequence represent previously methylated cytosines in the genome. The work described has defined procedures that maximise the efficiency of denaturation, bisulphite conversion and amplification, to permit methylation mapping of single genes from small amounts of genomic DNA, readily available from germ cells and early developmental stages.     

An isothermal method for whole genome amplification of fresh and degraded DNA for comparative genomic hybridization, genotyping and mutation detection.

Molecular genotyping has important biomedical and forensic applications. However, limiting amounts of human biological material often yield genomic DNA (gDNA) in insufficient quantity and of poor quality for a reliable analysis. This motivated the development of an efficient whole genome amplification method with quantitatively unbiased representation usable on fresh and degraded gDNA. Amplification of fresh frozen, formalin-fixed paraffin-embedded (FFPE) and DNase-degraded DNA using degenerate oligonucleotide-primed PCR or primer extension amplification using a short primer sequence bioinformatically optimized for coverage of the human genome was compared with amplification using current primers by chromosome-based and BAC-array comparative genomic hybridization (CGH), genotyping at short tandem repeats (STRs) and single base mutation detection. Compared with current primers, genome amplification using the bioinformatically optimized primer was significantly less biased on CGH in self-self hybridizations, and replicated tumour genome copy number aberrations, even from FFPE tissue. STR genotyping could be performed on degraded gDNA amplified using our technique but failed with multiple displacement amplification. Of the 18 different single base mutations 16 (89.5%) were correctly identified by sequencing gDNA amplified from clinical samples using our technique. This simple and efficient isothermal method should be helpful for genetic research and clinical and forensic applications.     

A comparison of existing global DNA methylation assays to low-coverage whole-genome bisulfite sequencing for epidemiological studies

DNA methylation is an epigenetic mark at the interface of genetic and environmental factors relevant to human disease. Quantitative assessments of global DNA methylation levels have therefore become important tools in epidemiology research, particularly for understanding effects of environmental exposures in complex diseases. Among the available methods of quantitative DNA methylation measurements, bisulfite sequencing is considered the gold standard, but whole-genome bisulfite sequencing (WGBS) has previously been considered too costly for epidemiology studies with high sample numbers. Pyrosequencing of repetitive sequences within bisulfite-treated DNA has been routinely used as a surrogate for global DNA methylation, but a comparison of pyrosequencing to WGBS for accuracy and reproducibility of methylation levels has not been performed. This study compared the global methylation levels measured from uniquely mappable (non-repetitive) WGBS sequences to pyrosequencing assays of several repeat sequences and repeat assay-matched WGBS data and determined uniquely mappable WGBS data to be the most reproducible and accurate measurement of global DNA methylation levels. We determined sources of variation in repetitive pyrosequencing assays to be PCR amplification bias, PCR primer selection bias in methylation levels of targeted sequences, and inherent variability in methylation levels of repeat sequences. Low-coverage, uniquely mappable WGBS showed the strongest correlation between replicates of all assays. By using multiplexing by indexed bar codes, the cost of WGBS can be lowered significantly to improve the accuracy of global DNA methylation assessments for human studies.     

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.     

Hi-Meth：特定位点DNA甲基化高通量检测平台

胞嘧啶甲基化是DNA表观遗传修饰的主要类型之一,在维持正常细胞功能和调控基因表达中具有重要作用。重亚硫酸盐测序法(bisulfite sequencing PCR,BSP)是特异性位点DNA甲基化检测的通用方法,能明确目的片段中每一个CpG位点的甲基化状态,但此方法需要大量的单克隆测序,操作过程较繁琐、成本昂贵。因此,开发准确、高效、便捷的DNA甲基化检测技术对提升表观遗传研究效率具有重要意义。基于本课题组开发的高通量突变类型检测平台Hi-TOM (high-throughput tracking of mutations),我们进一步建立了特定位点DNA甲基化高通量检测平台Hi-Meth (high-throughput detection of DNA methylation)。DNA样品通过重亚硫酸盐处理之后,仅需一轮PCR扩增即可通过Hi-Meth平台获得特定位点DNA甲基化分析结果。利用Hi-Meth平台,对水稻不同基因启动子区域进行了DNA甲基化检测分析,并与基于BSP方法获得的结果进行了比较。结果表明,Hi-Meth策略与BSP策略检测结果基本一致。而且通过Hi-Meth平台可以更准确、便捷地获得特异性位点DNA甲基化分析结果。综上所述,Hi-Meth为特定DNA区域提供了重要的甲基化检测平台,对表观遗传研究具有重要意义。     

Genome-Wide Analysis of DNA Methylation in Five Tissues of Zhikong Scallop,Chlamys farreri

DNA methylation plays a vital role in tissue development and differentiation in eukaryotes. Epigenetic studies have been seldom conducted in the extremely diverse and evolutionarily highly successful bilaterian lineage Mollusca. In the present study, we conducted the genome-wide profiling of DNA methylation for five tissues of a bivalve mollusc, Chlamys farreri using the methylation-sensitive amplification polymorphism (MSAP) technique. The methylation levels were quite similar among tissues, ranging from 20.9% to 21.7%. CG methylation was the dominant type (14.9%–16.5%) in the C. farreri genome, but CHG methylation also accounted for a substantial fraction of total methylation (5.1%–6.3%). Relatively high methylation diversity was observed within tissues. Methylation differentiation between tissues was evaluated and 460 tissue-specific epiloci were identified. Kidney differs from the other tissues in DNA methylation profiles. Our study presents the first look at the tissue-specific DNA methylation patterns in a bivalve mollusc and represents an initial step towards understanding of epigenetic regulatory mechanism underlying tissue development and differentiation in bivalves.     

Evaluation of a Quantitative DNA Methylation Analysis Technique using Methylation-Sensitive/Dependent Restriction Enzymes and Real-Time PCR

DNA methylation in mammals has been shown to play many important roles in diverse biological phenomena. Several methods have been developed for the measurement of region-specific levels of DNA methylation. We sought a technique that could be used to quantitatively evaluate multiple independent loci in several tissues in a quick and cost-effective manner. Recently, a few quantitative techniques have been developed by employing the use of real-time PCR, though they require the additional step of sodium bisulfite conversion. Here we evaluate a technique that involves the digestion of non-sodium bisulfite-treated genomic DNA using methylation-sensitive and methylation-dependent restriction enzymes followed by real-time PCR. The utility of this method is tested by analyzing seventeen genomic regions of known tissue-specific levels of DNA methylation including three imprinted genes. We find that this approach generates rapid, reproducible and accurate results (range= ±5%) without the additional time required for bisulfite conversion. This approach is also adaptable for use with smaller amounts of starting material. We propose this method as a rapid, quantitative method for the analysis of DNA methylation at single sites or within small regions of DNA.      

A novel method for whole blood PCR without pretreatment

PCR is usually performed on purified DNA. However, the extraction of DNA from whole blood is time consuming and involves the risk of contamination at every step. Hence, it is desirable to amplify DNA directly from whole blood. Earlier, investigators tried to achieve this target by either pretreatment of whole blood samples with different agents or by altering the conventional thermal cyclic conditions. This would make the technique cumbersome and time consuming. Here, we describe a simple protocol to amplify DNA directly from whole blood without the need of pretreatment. PCR buffer system was optimized in the laboratory and Apolipoprotein B gene was used as a model for this experiment. 480 bp was the target site for amplification. Fresh whole blood samples were used both from healthy and diseased individuals (coronary artery disease patients). Successful amplification was achieved with 1 μl volume of whole blood and it was comparable to that of genomic DNA. No pretreatment of whole blood samples was required with the optimized buffer system. 3mM concentration of MgCl(2) was observed to be optimal and hence used in the reaction mixture. Amplification was relatively better with this buffer system as compared to that of commercially available PCR buffer. With the present technique, amplicon detection did not require the centrifugation/dilution of the PCR products which further saves time. Successful amplification was achieved in both the healthy and diseased blood samples, indicating the robustness of the technique as changed blood composition and presence of increased inhibitory molecules in the diseased state did not seem to affect the efficacy of the present technique. In conclusion, as compared to the existing protocols for whole blood PCR, the present technique is relatively novel, simple, requires minimal steps and eliminates the need for additional standardizations.     

Assessing combined methylation-sensitive high resolution melting and pyrosequencing for the analysis of heterogeneous DNA methylation

Heterogeneous DNA methylation leads to difficulties in accurate detection and quantification of methylation. Methylation-sensitive high resolution melting (MS-HRM) is unique among regularly used methods for DNA methylation analysis in that heterogeneous methylation can be readily identified, although not quantified, by inspection of the melting curves. Bisulfite pyrosequencing has been used to estimate the level of heterogeneous methylation by quantifying methylation levels present at individual CpG dinucleotides. Sequentially combining the two methodologies using MS-HRM to screen the amplification products prior to bisulfite pyrosequencing would be advantageous. This would not only replace the quality control step using agarose gel analysis prior to the pyrosequencing step but would also provide important qualitative information in its own right. We chose to analyze DAPK1 as it is an important tumor suppressor gene frequently heterogeneously methylated in a number of malignancies, including chronic lymphocytic leukemia (CLL). A region of the DAPK1 promoter was analyzed in ten CLL samples by MS-HRM. By using a biotinylated primer, bisulfite pyrosequencing could be used to directly analyze the samples. MS-HRM revealed the presence of various extents of heterogeneous DAPK1 methylation in all CLL samples. Further analysis of the biotinylated MS-HRM products by bisulfite pyrosequencing provided quantitative information for each CpG dinucleotide analyzed, and confirmed the presence of heterogeneous DNA methylation. Whereas each method could be used individually, MS-HRM and bisulfite pyrosequencing provided complementary information for the assessment of heterogeneous methylation.Key words: MS-HRM, pyrosequencing, digital PCR, heterogeneous DNA methylation, DAPK1, chronic lymphocytic leukemia     

Capillary electrophoretic analysis of genomic DNA methylation levels

Changes in DNA methylation have been found in the large majority of tumors. This phenomenon includes both genome-wide hypomethylation and gene- specific hypermethylation. However, the clinical relevance of either mechanism has remained contentious. In order to determine DNA methylation levels from a large number of clinical samples, we have established a method for accurate high-throughput quantification of 5-methylcytosine in genomic DNA. Our protocol requires a small amount (<1 µg) of DNA that is enzymatically hydrolyzed to single nucleotides. Single nucleotides are then derivatized with a fluorescent marker and separated by capillary electrophoresis. After calibration of the method, we have determined cytosine methylation levels from tumor samples of 81 patients that had been diagnosed with chronic lymphocytic leukemia (CLL). These patients showed a high variability in their methylation levels with a general trend towards hypomethylation. Because of its high accuracy and throughput our method will be useful in determining the role of genomic DNA methylation levels in tumorigenesis.     

Correction of PCR-bias in quantitative DNA methylation studies by means of cubic polynomial regression

DNA methylation profiling has become an important aspect of biomedical molecular analysis. Polymerase chain reaction (PCR) amplification of bisulphite-treated DNA is a processing step that is common to many currently used methods of quantitative methylation analysis. Preferential amplification of unmethylated alleles-known as PCR-bias-may significantly affect the accuracy of quantification. To date, no universal experimental approach has been reported to overcome the problem. This study presents an effective method of correcting biased methylation data. The procedure includes a calibration performed in parallel to the analysis of the samples under investigation. DNA samples with defined degrees of methylation are analysed. The observed deviation of the experimental results from the expected values is used for calculating a regression curve. The equation of the best-fitting curve is then used for correction of the data obtained from the samples of interest. The process can be applied irrespective of the locus interrogated and the number of sites analysed, avoiding an optimization of the amplification conditions for each individual locus.     

Accurate detection of subclonal single nucleotide variants in whole genome amplified and pooled cancer samples using HaloPlex target enrichment

Background

Target enrichment and resequencing is a widely used approach for identification of cancer genes and genetic variants associated with diseases. Although cost effective compared to whole genome sequencing, analysis of many samples constitutes a significant cost, which could be reduced by pooling samples before capture. Another limitation to the number of cancer samples that can be analyzed is often the amount of available tumor DNA. We evaluated the performance of whole genome amplified DNA and the power to detect subclonal somatic single nucleotide variants in non-indexed pools of cancer samples using the HaloPlex technology for target enrichment and next generation sequencing.

Results

We captured a set of 1528 putative somatic single nucleotide variants and germline SNPs, which were identified by whole genome sequencing, with the HaloPlex technology and sequenced to a depth of 792–1752. We found that the allele fractions of the analyzed variants are well preserved during whole genome amplification and that capture specificity or variant calling is not affected. We detected a large majority of the known single nucleotide variants present uniquely in one sample with allele fractions as low as 0.1 in non-indexed pools of up to ten samples. We also identified and experimentally validated six novel variants in the samples included in the pools.

Conclusion

Our work demonstrates that whole genome amplified DNA can be used for target enrichment equally well as genomic DNA and that accurate variant detection is possible in non-indexed pools of cancer samples. These findings show that analysis of a large number of samples is feasible at low cost, even when only small amounts of DNA is available, and thereby significantly increases the chances of indentifying recurrent mutations in cancer samples.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-856) contains supplementary material, which is available to authorized users.     

Reference Materials for Calibration of Analytical Biases in Quantification of DNA Methylation

Most contemporary methods for the quantification of DNA methylation employ bisulfite conversion and PCR amplification. However, many reports have indicated that bisulfite-mediated PCR methodologies can result in inaccurate measurements of DNA methylation owing to amplification biases. To calibrate analytical biases in quantification of gene methylation, especially those that arise during PCR, we utilized reference materials that represent exact bisulfite-converted sequences with 0% and 100% methylation status of specific genes. After determining relative quantities using qPCR, pairs of plasmids were gravimetrically mixed to generate working standards with predefined DNA methylation levels at 10% intervals in terms of mole fractions. The working standards were used as controls to optimize the experimental conditions and also as calibration standards in melting-based and sequencing-based analyses of DNA methylation. Use of the reference materials enabled precise characterization and proper calibration of various biases during PCR and subsequent methylation measurement processes, resulting in accurate measurements.