Microarray-based methylation analysis using dual-color fluorescence hybridization

BACKGROUND: Aberrant DNA methylation of CpG sites is among the earliest and most frequent alterations in cancer. It is of great importance to develop simple, high-throughput and quantitative methods for methylation detection. METHODS: A high-throughput methylation analysis method has been developed based on microarray and dual-color fluorescence hybridization. The genomic DNA was treated with bisulfite, resulting in conversion of non-methylated cytosine, but not methylated cytosine, into uracil within CpG islands of interest. PCR products of the treated genomic templates were spotted and immobilized onto a poly-l-lysine coated glass slide to fabricate a microarray and then interrogated by hybridization with dual-color probes to determine the methylation status. The hybridized signals were obtained with a scanner and the results were analyzed with the software Genepix Pro 3.0. RESULTS: The methylation status of the CpG islands of IGFBP7 gene has been successfully evaluated by the microarray method for twenty-seven samples. All the investigated samples, including twenty human breast tumor tissues, six corresponding normal human breast tissues and one liver cell line, all CpG sites were found completely methylated. CONCLUSIONS: The microarray technology has been proven to have potential for high-throughput detection of the methylation status for a given gene in multi-genomic samples, which could be a novel approach for rapidly screening DNA methylation marker for early stage cancer diagnosis.     

Methylation assay by nucleotide incorporation: a quantitative assay for regional CpG methylation density

Although the aberrant methylation in CpG islands is of great interest as a causative role in human malignancies, it has been very difficult to accurately determine methylation density. Here we report a novel microplate-based quantitative methylation assay, designated MANIC, for a region containing a number of CpG sites based on incorporation of hapten-labeled dCTP at cytosine sites where the methylated cytosines have not been converted to uracil by the bisulfite treatment. Validation using control DNAs revealed that the method was sensitive enough to detect < 1.25% methylated DNA and that calibration curve was linear. With this approach, we determined relative methylation density of O6-methylguanine-DNA methyltransferase gene promoter containing 12 CpG sites among the 12 colorectal cancers and corresponding normal mucosal tissues. Consequently, MANIC showed a high concordance with results by a quantitative method, bisulfite PCR single-stranded conformational polymorphism (BiPS). MANIC is a technique that avoids cumbersome procedures such as electrophoresis or the use of radiolabeling and is applicable to any sequence regardless of the total number of CpG sites or heterogeneity in methylation status.     

High-throughput method for detecting DNA methylation

Aberrant DNA methylation of CpG site is among the earliest and most frequent alterations in cancer. Detection of promoter hypermethylation of cancer-related gene may be useful for cancer diagnosis or the detection of recurrence. However, most of the studies have focused on a single gene only and gave little information about the concurrent methylation status of multiple genes. In this study, we attempted to develop a microarray method coupled with linker-PCR for detecting methylation status of multiple genes in the tumor tissue. A series of synthesized oligonucleotides were synthesised and purified to completely match with 16 investigated targets. Then they were immobilized on the aldehyde-coated glass slide to fabricate a DNA microarray for detecting methylation status of these genes. The results indicated that these genes were all methylated in the positive control. However, no methylated was found in these genes for the negative control. Only p16 and p15 genes were methylated in investigated genes for the gastric tumor tissue, whereas others were not methylated. The above results were validated by bisulfite DNA sequencing. Our experiments successfully demonstrated that the DNA microarray could be applied as a high-throughput tool to determine methylation status of the investigated genes.     

Detection method for methylation density on microarray using methyl-CpG-binding domain protein

A method for determining methylation density of target CpG islands has been established. In the method, DNA microarray was prepared by spotting a set of PCR products amplified from bisulfite-converted sample DNAs. The PCR products on the microarray were treated by SssI methyltransferase and labeled with TAMRA fluorescence. A recombinant, antibody-like methyl-CpG-binding protein labeled with Cy5 fluorescence was used to identify symmetrical methyl-CpG dinucleotide of the PCR products on the microarray. By use of a standard curve with control mixtures, the ratio of two fluorescence signals can be converted into percentage values to assess methylation density of targeted fragments. We obtained the methylation density of six CpG islands on the two tumor suppressor genes of CDK2A and CDK2B from seven cancer cell line samples and two normal blood samples. The validity of this method was tested by bisulfite sequencing. This method not only allows the quantitative analysis of regional methylation density of a set of given genes but also could provide information of methylation density for a large amount of clinical samples.     

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

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

Methylation profiling using degenerated oligonucleotide primer-PCR specific for genome-wide amplification of bisulfite-modified DNA

DNA methylation is one of the essential epigenetic processes that play a role in regulating gene expression. Aberrant methylation of CpG-rich promoter regions has been associated with many forms of human cancers. The current method for determining the methylation status relies mainly on bisulfite treatment of genomic DNA, followed by methylation-specific PCR (MSP). The difficulty in acquiring a methylation profiling often is limited by the amount of genomic DNA that can be recovered from a given sample, whereas complex procedures of bisulfite treatment further compromise the effective template for PCR analysis. To circumvent these obstacles, we developed degenerated oligonucleotide primer (DOP)-PCR to enable amplification of bisulfite-modified genomic DNA at a genome-wide scale. A DOP pair was specially designed as follows: first 3' DOP, CTCGAGCTGHHHHHAACTAC, where H is a mixture of base consisting of 50% A, 25% T, and 25% C; and second 5' DOP, CTCGAGCTGDDDDDGTTTAG, where D is a mixture of base consisting of 50% T, 25% G, and 25% A. Our results showed that bisulfite-modified DNAs from a cell line, cord blood cells, or cells obtained by laser capture microdissection were amplified by up to 1000-fold using this method. Subsequent MSP analysis using these amplified DNAs on nine randomly selected cancer-related genes revealed that the methylation status of these genes remained identical to that derived from the original unamplified template.     

Accurate quantification of DNA methylation using combined bisulfite restriction analysis coupled with the Agilent 2100 Bioanalyzer platform

DNA methylation is the best-studied epigenetic modification and describes the conversion of cytosine to 5-methylcytosine. The importance of this phenomenon is that aberrant promoter hypermethylation is a common occurrence in cancer and is frequently associated with gene silencing. Various techniques are currently available for the analysis of DNA methylation. However, accurate and reproducible quantification of DNA methylation remains challenging. In this report, we describe Bio-COBRA (combined bisulfite restriction analysis coupled with the Agilent 2100 Bioanalyzer platform), as a novel approach to quantitative DNA methylation analysis. The combination of a well-established method, COBRA, which interrogates DNA methylation via the restriction enzyme analysis of PCR-amplified bisulfite treated DNAs, with the Bioanalyzer platform allows for the rapid and quantitative assessment of DNA methylation patterns in large sample sets. The sensitivity and reproducibility of Bio-COBRA make it a valuable tool for the analysis of DNA methylation in clinical samples, which could aid in the development of diagnostic and prognostic parameters with respect to disease detection and management.     

High density DNA methylation array with single CpG site resolution

We have developed a new generation of genome-wide DNA methylation BeadChip which allows high-throughput methylation profiling of the human genome. The new high density BeadChip can assay over 480K CpG sites and analyze twelve samples in parallel. The innovative content includes coverage of 99% of RefSeq genes with multiple probes per gene, 96% of CpG islands from the UCSC database, CpG island shores and additional content selected from whole-genome bisulfite sequencing data and input from DNA methylation experts. The well-characterized Infinium® Assay is used for analysis of CpG methylation using bisulfite-converted genomic DNA. We applied this technology to analyze DNA methylation in normal and tumor DNA samples and compared results with whole-genome bisulfite sequencing (WGBS) data obtained for the same samples. Highly comparable DNA methylation profiles were generated by the array and sequencing methods (average R² of 0.95). The ability to determine genome-wide methylation patterns will rapidly advance methylation research.     

A systematic study of normalization methods for Infinium 450K methylation data using whole-genome bisulfite sequencing data

DNA methylation plays an important role in disease etiology. The Illumina Infinium HumanMethylation450 (450K) BeadChip is a widely used platform in large-scale epidemiologic studies. This platform can efficiently and simultaneously measure methylation levels at ∼480,000 CpG sites in the human genome in multiple study samples. Due to the intrinsic chip design of 2 types of chemistry probes, data normalization or preprocessing is a critical step to consider before data analysis. To date, numerous methods and pipelines have been developed for this purpose, and some studies have been conducted to evaluate different methods. However, validation studies have often been limited to a small number of CpG sites to reduce the variability in technical replicates. In this study, we measured methylation on a set of samples using both whole-genome bisulfite sequencing (WGBS) and 450K chips. We used WGBS data as a gold standard of true methylation states in cells to compare the performances of 8 normalization methods for 450K data on a genome-wide scale. Analyses on our dataset indicate that the most effective methods are peak-based correction (PBC) and quantile normalization plus β-mixture quantile normalization (QN.BMIQ). To our knowledge, this is the first study to systematically compare existing normalization methods for Illumina 450K data using novel WGBS data. Our results provide a benchmark reference for the analysis of DNA methylation chip data, particularly in white blood cells.     

巢式甲基化特异性PCR检测不同类型组织标本中WIF-1基因启动子异常甲基化

目的:采用一种高灵敏度的DNA甲基化分析方法,即巢式甲基化特异性PCR法(nested-MSP,nMSP),检测外科手术切除新鲜组织、石蜡包埋组织及纤维支气管镜活检组织中WIF-1基因启动子的异常甲基化状态。方法:将基因组DNA变性成为单链,用亚硫酸氢盐修饰单链DNA,所有未甲基化的胞嘧啶被转变为尿嘧啶,而甲基化的胞嘧啶则不变。设计针对甲基化和非甲基化等位基因的特异引物,进行巢式PCR扩增,最后经凝胶电泳检测目的片段。结果:在3种类型的原发性非小细胞肺癌标本中都检测出了WIF-1基因启动子的异常甲基化。结论:巢式甲基化特异性PCR是一种灵敏度高、特异性强的甲基化检测方法,可广泛应用于不同类型标本基因启动子甲基化的分析。     

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.     

Droplet digital PCR is an accurate method to assess methylation status on FFPE samples

Most tissue samples available for cancer research are archived as formalin-fixed paraffin-embedded (FFPE) samples. However, the fixation process and the long storage duration lead to DNA fragmentation and hinder epigenome analysis. The use of droplet digital PCR (ddPCR) to detect DNA methylation has recently emerged. In this study, we compare an optimized ddPCR assay with a conventional qPCR assay by targeting a dilution series of control DNA. In addition, we compare the ddPCR technology with results from Infinium arrays targeting two separate CpG sites on a set of colon adenoma FFPE samples. Our data demonstrate that qPCR and ddPCR assess methylation status equally well on dilution controls with a high DNA input. However, the methylation detection on low-input samples is more accurate using ddPCR. The proposed primer design (methylation-independent primers with amplification of solely the converted DNA target) will allow for methylation detection, independent of bisulfite conversion efficiency. Those data show that ddPCR can be used for methylation analysis on FFPE samples with a wide range of DNA input and that the precision of the assay depends largely on the total amount of amplifiable DNA fragments. Due to accessibility of the ddPCR technology and its accuracy on high- as well as low-DNA input samples, we propose the use of this approach for studies involving degraded FFPE samples.     

Reducing the risk of false discovery enabling identification of biologically significant genome-wide methylation status using the HumanMethylation450 array

Background

The Illumina HumanMethylation450 BeadChip (HM450K) measures the DNA methylation of 485,512 CpGs in the human genome. The technology relies on hybridization of genomic fragments to probes on the chip. However, certain genomic factors may compromise the ability to measure methylation using the array such as single nucleotide polymorphisms (SNPs), small insertions and deletions (INDELs), repetitive DNA, and regions with reduced genomic complexity. Currently, there is no clear method or pipeline for determining which of the probes on the HM450K bead array should be retained for subsequent analysis in light of these issues.

Results

We comprehensively assessed the effects of SNPs, INDELs, repeats and bisulfite induced reduced genomic complexity by comparing HM450K bead array results with whole genome bisulfite sequencing. We determined which CpG probes provided accurate or noisy signals. From this, we derived a set of high-quality probes that provide unadulterated measurements of DNA methylation.

Conclusions

Our method significantly reduces the risk of false discoveries when using the HM450K bead array, while maximising the power of the array to detect methylation status genome-wide. Additionally, we demonstrate the utility of our method through extraction of biologically relevant epigenetic changes in prostate cancer.

Electronic supplementary material

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

Library-free methylation sequencing with bisulfite padlock probes

Targeted quantification of DNA methylation allows for interrogation of the most informative loci across many samples quickly and cost-effectively. Here we report improved bisulfite padlock probes (BSPPs) with a design algorithm to generate efficient padlock probes, a library-free protocol that dramatically reduces sample-preparation cost and time and is compatible with automation, and an efficient bioinformatics pipeline to accurately obtain both methylation levels and genotypes from sequencing of bisulfite-converted DNA.     

Bisulfite-modified target DNA array for aberrant methylation analysis

Aberrant DNA methylation of CpG islands is among the earliest and most frequent alterations in cancer. It is of great importance to develop simple and high-throughput methods of methylation analysis for earlier cancer diagnosis or the detection of recurrence. In this study, bisulfite-modified target DNA arrays were prepared on positively charged nylon membrane with two different procedures: fixing PCR products and fixing genomic DNA. First, a bisulfite PCR product array was prepared through fixing PCR products amplified in bisulfite sequencing primers from the bisulfite-modified genomic DNA of different clinical samples on membrane. Furthermore, bisulfite-modified genomic DNA of the different samples was directly fixed on membrane to fabricate bisulfite genomic DNA arrays. The two kinds of arrays were hybridized by probes labeled with digoxigenin, and the hybridization signals were obtained through chemiluminescent detection. The methylation statuses of the IGFBP7 gene for breast tumor and normal tissue samples and for normal human blood cell samples were detected successfully by the two procedures. It was shown that the methods are reliable and sensitive and that they have high potential in screening molecular methylation markers from a large number of clinical samples.     

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.     

Sensitive measurement of unmethylated repeat DNA sequences by end-specific PCR

We describe a new method that is well-suited for the determination of the methylation level of repetitive sequences such as human Alu elements. We have applied the method to the analysis of cell and tissue DNAs and expect it to have wide utility in studies of DNA methylation in cancer and other disease states, in monitoring response to epigenetic cancer therapies and in epidemiological studies. Only 1 ng DNA is needed for a duplex, one-tube real-time PCR in which methylation level and the amount of input DNA are concurrently measured. The relative cutting by the methylation-sensitive enzyme BstUI is compared with that of the methylation-insensitive enzyme DraI to give a measure of DNA methylation. The method depends upon the use of 5'-tailed, 3'-blocked oligonucleotides called facilitator oligonucleotides (Foligos). Only cut DNAs with specific matching sequences at their 3' ends can copy the tails of the Foligos and thus become tagged and available for subsequent PCR. Both the tagging and PCR are carried out by the same enzyme, Taq DNA polymerase. Because amplification only occurs if suitable ends have been generated in the target DNA, we have called this method end-specific PCR (ESPCR). ESPCR avoids the bisulfite treatment step that is usually required to measure methylation.     

Microsoft Word macro for analysis of cytosine methylation by the bisulfite deamination reaction

Cytosine methylation at CpG dinucleotides is an important control mechanism in development, differentiation, and neoplasia. Bisulfite genomic sequencing and its modifications have been developed to examine methylation at these CpG dinucleotides. To use these methods, one has to (i) manually convert the sequence to that produced by bisulfite conversion and PCR amplification, taking into account that cytosine residues at CpG dinucleotides may or may not be converted depending on their methylation status, (ii) identify relevant restriction sites that may be used for methylation analysis, and (iii) conduct similar steps with the other DNA strand since the two strands of DNA are no longer complementary after bisulfite conversion. To automate these steps, we have developed a macro that can be used with Microsoft Word. This macro (i) converts genomic sequence to modified sequence that would result after bisulfite treatment facilitating primer design for bisulfite genomic sequencing and methylation-sensitive PCR assay and (ii) identifies restriction sites that are preserved in bisulfite-converted and PCR-amplified product only if cytosine residues at relevant CpG dinucleotides are methylated (and thereby not converted to uracil) in the genomic DNA.     

Advanced methylome analysis after bisulfite deep sequencing: an example in Arabidopsis

Deep sequencing after bisulfite conversion (BS-Seq) is the method of choice to generate whole genome maps of cytosine methylation at single base-pair resolution. Its application to genomic DNA of Arabidopsis flower bud tissue resulted in the first complete methylome, determining a methylation rate of 6.7% in this tissue. BS-Seq reads were mapped onto an in silico converted reference genome, applying the so-called 3-letter genome method. Here, we present BiSS (Bisufite Sequencing Scorer), a new method applying Smith-Waterman alignment to map bisulfite-converted reads to a reference genome. In addition, we introduce a comprehensive adaptive error estimate that accounts for sequencing errors, erroneous bisulfite conversion and also wrongly mapped reads. The re-analysis of the Arabidopsis methylome data with BiSS mapped substantially more reads to the genome. As a result, it determines the methylation status of an extra 10% of cytosines and estimates the methylation rate to be 7.7%. We validated the results by individual traditional bisulfite sequencing for selected genomic regions. In addition to predicting the methylation status of each cytosine, BiSS also provides an estimate of the methylation degree at each genomic site. Thus, BiSS explores BS-Seq data more extensively and provides more information for downstream analysis.