MBD-isolated Genome Sequencing provides a high-throughput and comprehensive survey of DNA methylation in the human genome

DNA methylation is an epigenetic modification involved in both normal developmental processes and disease states through the modulation of gene expression and the maintenance of genomic organization. Conventional methods of DNA methylation analysis, such as bisulfite sequencing, methylation sensitive restriction enzyme digestion and array-based detection techniques, have major limitations that impede high-throughput genome-wide analysis. We describe a novel technique, MBD-isolated Genome Sequencing (MiGS), which combines precipitation of methylated DNA by recombinant methyl-CpG binding domain of MBD2 protein and sequencing of the isolated DNA by a massively parallel sequencer. We utilized MiGS to study three isogenic cancer cell lines with varying degrees of DNA methylation. We successfully detected previously known methylated regions in these cells and identified hundreds of novel methylated regions. This technique is highly specific and sensitive and can be applied to any biological settings to identify differentially methylated regions at the genomic scale.     

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区域提供了重要的甲基化检测平台,对表观遗传研究具有重要意义。     

Melting curve analysis of SNPs (McSNP): a gel-free and inexpensive approach for SNP genotyping

High-throughput methods for assaying DNA variation require two important steps: (i) discriminating the variation and (ii) detecting the signal. In this report, we describe a novel SNP genotyping method that we refer to as melting curve analysis of SNPs (McSNP). McSNP combines a classic approach for discriminating alleles, restriction enzyme digestion, with a more recent method for detecting DNA fragments, melting curve analysis. Melting curve analysis is performed by slowly heating DNA fragments in the presence of the dsDNA-specific fluorescent dye SYBR Green I. As the sample is heated, fluorescence rapidly decreases when the melting temperature of a particular fragment is reached. We show that it is possible to determine the composition of simple mixtures of DNA fragments, such as those that result from restriction enzyme digestions of short PCR products. McSNP is well suited for high-throughput genotyping because 96 samples can be analyzed and automatically scored in 20 min. Our results clearly demonstrate that McSNP is a simple, inexpensive, and accurate means of genotyping SNP variation.     

Detection and analysis of enzymatic DNA methylation of oligonucleotide substrates by matrix-assisted laser desorption ionization time-of-flight mass spectrometry

Matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF) mass spectrometry was employed to analyze DNA methylation carried out by the Escherichia coli dam DNA methyltransferase using oligonucleotide substrates with molecular masses of 5000-10,000 Da per strand. The mass spectrometry assay offers several advantages: (i) it directly shows the methylation as the increase in the mass of the substrate DNA, (ii) it is nonradioactive, (iii) it is quantitative, and (iv) it can be automated for high-throughput applications. Since unmethylated and methylated DNA are detected, the ratio of methylation can be determined directly and accurately. Furthermore, the assay allows detection individually of the methylation of several substrates in competition, offering an ideal setup to analyze the specificity of DNA interacting with enzymes. We could not identify methylation at any noncanonical site, indicating that the dam MTase is a very specific enzyme. Finally, MALDI-TOF mass spectrometry permitted assessment of the number of methyl groups incorporated into each DNA strand, thereby, allowing study of mechanistic details such as the processivity of the methylation reaction. We provide evidence that the dam MTase modifies DNA in a processive reaction, confirming earlier findings.     

Region of interest methylation analysis: a comparison of MSP with MS-HRM and direct BSP

The aim of this study was to compare and contrast three DNA methylation methods of a specific region of interest (ROI): methylation-specific PCR (MSP), methylation-sensitive high resolution melting (MS-HRM) and direct bisulfite sequencing (BSP). The methylation of a CpG area in the promoter region of Estrogen receptor alpha (ESR1) was evaluated by these three methods with samples and standards of different methylation percentages. MSP data were neither reproducible nor sensitive, and the assay was not specific due to non-specific binding of primers. MS-HRM was highly reproducible and a step forward into categorizing the methylation status of the samples as percent ranges. Direct BSP was the most informative method regarding methylation percentage of each CpG site. Though not perfect, it was reproducible and sensitive. We recommend the use of either method depending on the research question and target amplicon, and provided that the designed primers and expected amplicons are within recommendations. If the research question targets a limited number of CpG sites and simple yes/no results are enough, MSP may be attempted. For short amplicons that are crowded with CpG sites and of single melting domain, MS-HRM may be the method of choice though it only indicates the overall methylation percentage of the entire amplicon. Although the assay is highly reproducible, being semi-quantitative makes it of lesser interest to study ROI methylation of samples with little methylation differences. Direct BSP is a step forward as it gives information about the methylation percentage at each CpG site.     

Heterogeneity and randomness of DNA methylation patterns in human embryonic stem cells

DNA methylation has been proposed to be important in many biological processes and is the subject of intense study. Traditional bisulfite genomic sequencing allows detailed high-resolution methylation pattern analysis of each molecule with haplotype information across a few hundred bases at each locus, but lacks the capacity to gather voluminous data. Although recent technological developments are aimed at assessing DNA methylation patterns in a high-throughput manner across the genome, the haplotype information cannot be accurately assembled when the sequencing reads are short or when each hybridization target only includes one or two cytosine-phosphate-guanine (CpG) sites. Whether a distinct and nonrandom DNA methylation pattern is present at a given locus is difficult to discern without the haplotype information, and the DNA methylation patterns are much less apparent because the data are often obtained only as methylation frequencies at each CpG site with some of these methods. It would facilitate the interpretation of data obtained from high-throughput bisulfite sequencing if the loci with nonrandom DNA methylation patterns could be distinguished from those that are randomly methylated. In this study, we carried out traditional genomic bisulfite sequencing using the normal diploid human embryonic stem (hES) cell lines, and utilized Hamming distance analysis to evaluate the existence of a distinct and nonrandom DNA methylation pattern at each locus studied. Our findings suggest that Hamming distance is a simple, quick, and useful tool to identify loci with nonrandom DNA methylation patterns and may be utilized to discern links between biological changes and DNA methylation patterns in the high-throughput bisulfite sequencing data sets.     

Rapid quantitation of methylation differences at specific sites using methylation-sensitive single nucleotide primer extension (Ms-SNuPE).

We have developed a rapid quantitative method (Ms-SNuPE) for assessing methylation differences at specific CpG sites based on bisulfite treatment of DNA followed by single nucleotide primer extension. Genomic DNA was first reacted with sodium bisulfite to convert unmethylated cytosine to uracil while leaving 5-methylcytosine unchanged. Amplification of the desired target sequence was then performed using PCR primers specific for bisulfite-converted DNA and the resulting product isolated and used as a template for methylation analysis at the CpG site(s) of interest. This methylation-sensitive technique has several advantages over existing methods used for detection of methylation changes because small amounts of DNA can be analyzed including microdissected pathology sections and it avoids utilization of restriction enzymes for determining the methylation status at CpG sites.     

Methylation-sensitive high resolution melting (MS-HRM): a new approach for sensitive and high-throughput assessment of methylation

In this article, we show that high resolution melting analysis (HRM) is a sensitive and specific method for the detection of methylation. Methylated DNA and unmethylated DNA acquire different sequences after bisulphite treatment resulting in PCR products with markedly different melting profiles. We used PCR to amplify both methylated and unmethylated sequences and assessed HRM for the determination of the methylation status of the MGMT promoter region. Reconstruction experiments showed that MGMT methylation could be detected at levels as low as 0.1%. Moreover, MS-HRM allows for estimation of the methylation level by comparing the melting profiles of unknown PCR products to the melting profiles of PCR products derived from standards with a known unmethylated to methylated template ratio. We used MS-HRM for the analysis of eight cell lines of known methylation status and a panel of colorectal cancer specimens. The simplicity and high reproducibility of the MS-HRM protocol makes MS-HRM the method of choice for methylation assessment in many diagnostic and research applications.     

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.     

Microarray-based Ms-SNuPE: near-quantitative analysis for a high-throughput DNA methylation

Aberrant DNA methylation of CpG site in the gene promoter region has been confirmed to be closely associated with carcinogenesis. In the present study, a microarray-based methylation-sensitive single-nucleotide primer extension (Ms-SNuPE) for parallel detecting changes of DNA methylation in cancer was developed. After modification by sodium sulfite, the unmethylated cytosine in the genomic DNA is converted to uracil while leaving the 5-methylcytosine unchanged, which can be detected by bifunctional primer carrying a unique sequence tag in addition to a locus-specific sequence. Because each locus has a distinct tag, the detecting reactions can be performed in a highly multiplexed fashion and the resulting product then be hybridized to the reverse complements of the sequence tags arrayed on a glass slide for methylation analysis. The calibration curves with the correlation coefficient >0.97 were established, which suggested that the method could be used in near-quantitative DNA methylation analysis. Two breast tumor-related genes (E-cad and p16) are successfully analyzed by two group primers (22 primers total), and the results are compatible with that of methylation-specific PCR (MSP). Our research proved that the method is simple and inexpensive, and could be applied as a high-throughput tool to quantitatively determine methylation status of the investigated genes.     

High-throughput assay of DNA methylation based on methylation-specific primer and SAGE

Mapping of genomic DNA methylation is a dispensable part of functional genome. We have developed a novel method based on methylation-specific primer and serial analysis of gene expression, called MSP-SAGE, with potential of high-throughput quantification of genomic DNA methylation. We used a 6-mer methylation-specific primer to extend the methylated CpG sequences other than non-methylated CpG sequences. The 17 bp tags contained methylated CpG sequence, which were obtained from extended methylation sequence by digestion of restriction endonuclease, and then the tags were concatenated and cloned for sequencing. We can identify the locations of methylation according to the sequences of tags and quantify the methylation status from the frequency of the tags. MSP-SAGE has a good linearity in a broad methylation range from 5% to 100% with good accuracy and high precision. The proof-of-principle study shows that MSP-SAGE is a reliable high-throughput assay for quantification of DNA methylation.     

Explorative screening of complex microbial communities by real-time 16S rDNA restriction fragment melting curve analyses

We have developed restriction fragment melting curve analyses (RFMCA), which is a novel method for the real-time analysis of microbial communities. The major advantage of RFMCA compared to, for example, terminal restriction fragment length polymorphism (T-RFLP) or temperature/denaturing gradient gel electrophoresis (TGGE/DGGE) is that the physical separation of DNA fragments is avoided. The RFMCA detection is done by melting point analyses in closed tube systems, which enables high-throughput applications. The robustness of RFMCA was demonstrated by analyzing both mixtures of known samples and the microbial communities in the cecal content of poultry. Our conclusions are that RFMCA is robust, gives a relatively high resolution, and has the potential for high-throughput explorative screenings of microbial communities and large clone libraries.     

Methylation-sensitive single-nucleotide primer extension (Ms-SNuPE) for quantitative measurement of DNA methylation

Methylation-sensitive single-nucleotide primer extension (Ms-SNuPE) is a technique that can be used for rapid quantitation of methylation at individual CpG sites. Treatment of genomic DNA with sodium bisulfite is used to convert unmethylated Cytosine to Uracil while leaving 5-methylcytosine unaltered. Strand-specific PCR is performed to generate a DNA template for quantitative methylation analysis using Ms-SNuPE. SNuPE is then performed with oligonucleotide(s) designed to hybridize immediately upstream of the CpG site(s) being interrogated. Reaction products are electrophoresed on polyacrylamide gels for visualization and quantitation by phosphorimage analysis. The Ms-SNuPE technique is similar to other quantitative assays that use bisulfite treatment of genomic DNA to discriminate unmethylated from methylated Cytosines (i.e., COBRA, pyrosequencing). Ms-SNuPE can be used for high-throughput methylation analysis and rapid quantitation of Cytosine methylation suitable for a wide range of biological investigations, such as checking aberrant methylation changes during tumorigenesis, monitoring methylation changes induced by DNA methylation inhibitors or for measuring hemimethylation. Approximately two to four CpG sites can be interrogated in up to 40 samples by Ms-SNuPE in less than 5 h, after PCR amplification of the desired target sequence and preparation of PCR amplicons.     

DNA methylation profiling in cancer: from single nucleotides towards methylome

Genomic DNA methylation pattern (methylome) represents epigenetic program of a cell. It controls expression of genetic information. In tumor cells, significant alterations in DNA methylation take place, which can be identified as one of the earliest and most consistent features of tumorigenesis. Detailed survey of methylcytosines' distribution in genome is extremely important for understanding of real tumor etiology and early diagnostics. Progress in the field has been hampered by the unavailability of methods for large-scale determination of methylation patterns. Nowadays, variety of techniques is in development that allow for highly parallel regime of samples analysis (high-throughput analysis) or large loci DNA profiling (large-scale analysis). Aim of the work is to consider the main trends in the field of new methods development. The principles of the most frequently used approaches to DNA methylation studies are reviewed as well as their application and results. Most attention is paid to DNA microarrays as a technology of choice for epigenetic tumor analysis (oligonucleotide microarrays, BAC-arrays etc.). Alternative DNA sequencing based techniques are discussed, which can soon take on the leadership. Results of a large-scale analysis can be used for identification of new epigenetic markers and epigenetic classification of neoplasia.     

Differentiating between monozygotic twins through DNA methylation-specific high-resolution melt curve analysis

Although short tandem repeat profiling is extremely powerful in identifying individuals from crime scene stains, it is unable to differentiate between monozygotic (MZ) twins. Efforts to address this include mutation analysis through whole genome sequencing and through DNA methylation studies. Methylation of DNA is affected by environmental factors; thus, as MZ twins age, their DNA methylation patterns change. This can be characterized by bisulfite treatment followed by pyrosequencing. However, this can be time-consuming and expensive; thus, it is unlikely to be widely used by investigators. If the sequences are different, then in theory the melting temperature should be different. Thus, the aim of this study was to assess whether high-resolution melt curve analysis can be used to differentiate between MZ twins. Five sets of MZ twins provided buccal swabs that underwent extraction, quantification, bisulfite treatment, polymerase chain reaction amplification and high-resolution melting curve analysis targeting two markers, Alu-E2F3 and Alu-SP. Significant differences were observed between all MZ twins targeting Alu-E2F3 and in four of five MZ twins targeting Alu-SP (P < 0.05). Thus, it has been demonstrated that bisulfite treatment followed by high-resolution melting curve analysis could be used to differentiate between MZ twins.     

DNA sequencing and helix–coil transition. III. DNA sequencing

We show that the fine oscillatory structure of the DNA melting curve can be used to determine explicitly the nucleotide composition and the order of certain domains within the DNA. If DNA is specifically fragmented, the order of fragments can be learned directly from a comparison of the differential melting curves of the nonfragmented and fragmented DNA. The indicated information may complement exact methods of DNA sequencing. The proposed analysis is applied to bacteriophage ?X-174, whose melting curve is known. Compared to the known ?X-174 DNA sequence, the results of the analysis are found to be very accurate.