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.     

Control of carry-over contamination for PCR-based DNA methylation quantification using bisulfite treated DNA

In this study, we adapted the well known uracil DNA glycosylase (UNG) carry-over prevention system for PCR, and applied it to the analysis of DNA methylation based on sodium bisulfite conversion. As sodium bisulfite treatment converts unmethylated cytosine bases into uracil residues, bisulfite treated DNA is sensitive to UNG treatment. Therefore, UNG cannot be used for carry-over prevention of PCR using bisulfite treated template DNA, as not only contaminating products of previous PCR, but also the actual template will be degraded. We modified the bisulfite treatment procedure and generated DNA containing sulfonated uracil residues. Surprisingly, and in contrast to uracil, 6-sulfonyl uracil containing DNA (SafeBis DNA) is resistant to UNG. We showed that the new procedure removes up to 10 000 copies of contaminating PCR product in a closed PCR vessel without significant loss of analytical or clinical sensitivity of the DNA methylation analysis.     

A Modified Protocol for Bisulfite Genomic Sequencing of Difficult Samples

The bisulfite genomic sequencing protocol is a widely used method for analyzing DNA methylation. It relies on the deamination of unmethylated cytosine residues to uracil; however, its high rates of DNA degradation and incomplete cytosine to uracil conversion often lead to failed experiments, uninformative results, and false positives. Here, we report the addition of a single-step multiple restriction enzyme digestion (MRED) designed to differentially digest polymerase chain reaction products amplified from unconverted DNA while leaving those of converted DNA intact. We show that for our model system, RARB2 P2 promoter, use of MRED increased informative sequencings ninefold, and MRED did not alter the clonal representation in one fully methylated cell line, H-596, treated or not with 5-azadeoxycytidine, a methylation inhibitor. We believe that this method may easily be adapted for analyzing other genes and provide guidelines for selecting the most appropriate MRED restriction enzymes.     

一种用于DNA甲基化分析的改进型亚硫酸氢盐转化法

DNA甲基化分析是认识生理、病理条件下基因表达变化的重要途径.亚硫酸氢盐转化是DNA甲基化分析的瓶颈.本文旨在改进琼脂糖 亚硫酸氢盐DNA处理方案(agarose bisulfite method),建立一种简便稳定、适合常规甲基化分析的亚硫酸氢盐转化法.把DNA包入普通琼脂糖,以饱和亚硫酸氢盐在较高的温度下快速处理,然后用离心柱型琼脂糖凝胶DNA回收试剂盒,集DNA凝胶回收、脱盐、脱磺基和纯化于一体,完成整个转化过程.Bisulfite-PCR、克隆测序和酶切法分析转化率、转化特异性和转化物的质量.用该方案处理的HeLa细胞DNA,多个片段的转化率均大于98%,甲基化片段96.2%的CpG保持不变,可以扩增605 bp的较大片段,灵敏度介于普通法和琼脂糖亚硫酸氢盐法之间,而重复性较二者都好.改良后的方案简化了操作流程,快速稳定,易学易用,可实现高效特异转化,适合于一般实验者对常规检材进行DNA甲基化分析.     

The bisulfite genomic sequencing used in the analysis of epigenetic states, a technique in the emerging environmental genotoxicology research

Methylation at position 5 of cytosine in DNA is an important event in epigenetic changes of cells, the methylation being linked to the control of gene functions. The DNA methylation can be analyzed by bisulfite genomic sequencing, and a large body of data have now been accumulated, based on which causation of diseases, for example cancer, and many other manifestations of cellular activities have been discussed intensively. This article gives an extensive account of the chemical aspects of bisulfite modification of cytosine and 5-methylcytosine in DNA. Various factors that affect the action of bisulfite are discussed, and a recent progress from our laboratory is explained. Conventional procedures for the bisulfite treatment consist incubation of single-stranded DNA with sodium bisulfite under acidic conditions. This treatment converts cytosine into uracil, but 5-methylcytosine remains unchanged. Amplification by polymerase chain reaction (PCR) of the bisulfite-treated DNA followed by sequencing can result in revealing the positions of 5-methylcytosine in the gene. We have discovered that the whole procedure can be significantly speeded up by the use of a highly concentrated bisulfite solution, 10 M ammonium bisulfite. Another recent finding is that urea, which has been often added to the reaction mixture with the purpose of facilitating the bisulfite-mediated deamination of cytosine in DNA, may not work as anticipated: we have observed that urea does not show such promoting actions in our treatments of DNA. A laboratory protocol for quantifying bisulfite, suitably simple for routine practice to ensure valid experiments, is described.     

Enhanced Reduced Representation Bisulfite Sequencing for Assessment of DNA Methylation at Base Pair Resolution

DNA methylation pattern mapping is heavily studied in normal and diseased tissues. A variety of methods have been established to interrogate the cytosine methylation patterns in cells. Reduced representation of whole genome bisulfite sequencing was developed to detect quantitative base pair resolution cytosine methylation patterns at GC-rich genomic loci. This is accomplished by combining the use of a restriction enzyme followed by bisulfite conversion. Enhanced Reduced Representation Bisulfite Sequencing (ERRBS) increases the biologically relevant genomic loci covered and has been used to profile cytosine methylation in DNA from human, mouse and other organisms. ERRBS initiates with restriction enzyme digestion of DNA to generate low molecular weight fragments for use in library preparation. These fragments are subjected to standard library construction for next generation sequencing. Bisulfite conversion of unmethylated cytosines prior to the final amplification step allows for quantitative base resolution of cytosine methylation levels in covered genomic loci. The protocol can be completed within four days. Despite low complexity in the first three bases sequenced, ERRBS libraries yield high quality data when using a designated sequencing control lane. Mapping and bioinformatics analysis is then performed and yields data that can be easily integrated with a variety of genome-wide platforms. ERRBS can utilize small input material quantities making it feasible to process human clinical samples and applicable in a range of research applications. The video produced demonstrates critical steps of the ERRBS protocol.     

Kismeth: Analyzer of plant methylation states through bisulfite sequencing

Background  

There is great interest in probing the temporal and spatial patterns of cytosine methylation states in genomes of a variety of organisms. It is hoped that this will shed light on the biological roles of DNA methylation in the epigenetic control of gene expression. Bisulfite sequencing refers to the treatment of isolated DNA with sodium bisulfite to convert unmethylated cytosine to uracil, with PCR converting the uracil to thymidine followed by sequencing of the resultant DNA to detect DNA methylation. For the study of DNA methylation, plants provide an excellent model system, since they can tolerate major changes in their DNA methylation patterns and have long been studied for the effects of DNA methylation on transposons and epimutations. However, in contrast to the situation in animals, there aren't many tools that analyze bisulfite data in plants, which can exhibit methylation of cytosines in a variety of sequence contexts (CG, CHG, and CHH).     

Bisulfite genomic sequencing: systematic investigation of critical experimental parameters

Bisulfite genomic sequencing is the method of choice for the generation of methylation maps with single-base resolution. The method is based on the selective deamination of cytosine to uracil by treatment with bisulfite and the sequencing of subsequently generated PCR products. In contrast to cytosine, 5-methylcytosine does not react with bisulfite and can therefore be distinguished. In order to investigate the potential for optimization of the method and to determine the critical experimental parameters, we determined the influence of incubation time and incubation temperature on the deamination efficiency and measured the degree of DNA degradation during the bisulfite treatment. We found that maximum conversion rates of cytosine occurred at 55°C (4–18 h) and 95°C (1 h). Under these conditions at least 84–96% of the DNA is degraded. To study the impact of primer selection, homologous DNA templates were constructed possessing cytosine-containing and cytosine-free primer binding sites, respectively. The recognition rates for cytosine (≥97%) and 5-methylcytosine (≥94%) were found to be identical for both templates.     

应用亚硫酸氢盐修饰后PCR联合TA克隆测序对E-cadherin启动子区甲基化水平的检测

E-cadherin是一种细胞粘附因子,通过增强细胞之间的粘附而起到抑制肿瘤转移的作用.Ecadherin基因启动子区的高甲基化是导致其在众多肿瘤细胞中表达下调甚至缺失的主要原因之一.本实验首先抽提SGC-7901细胞(胃腺癌细胞)、A549细胞(肺腺癌细胞)、MCF-7细胞(乳腺癌细胞)等3个肿瘤细胞株的全基因组DNA,然后对抽提的DNA进行亚硫酸氢盐修饰和纯化回收,根据修饰后的DNA序列设计引物并对其进行PCR扩增.然后将PCR扩增产物与pUC-T TA载体连接并转化入感受态大肠杆菌DH5α中进行培养,对筛选出的含有阳性重组子的菌落进行测序.测序结果显示,3个肿瘤细胞株的E-cadherin基因启动子区的CpG岛都呈现了高度的甲基化,亚硫酸氢盐的修饰效率达到了99.2%.综上研究表明,亚硫酸氢盐修饰后PCR(BSP)联合TA克隆测序可以对肿瘤细胞某基因启动子区CpG岛的甲基化水平进行精确量化,研究所使用的3个肿瘤细胞株均可作为研究肿瘤细胞E-cadherin基因甲基化的细胞模型.     

A new method for accurate assessment of DNA quality after bisulfite treatment

The covalent addition of methylgroups to cytosine has become the most intensively researched epigenetic DNA marker. The vast majority of technologies used for DNA methylation analysis rely on a chemical reaction, the so-called ‘bisulfite treatment’, which introduces methylation-dependent sequence changes through selective chemical conversion of non-methylated cytosine to uracil. After treatment, all non-methylated cytosine bases are converted to uracil but all methylated cytosine bases remain cytosine. These methylation dependent C-to-T changes can subsequently be studied using conventional DNA analysis technologies.

The bisulfite conversion protocol is susceptible to processing errors, and small deviation from the protocol can result in failure of the treatment. Several attempts have been made to simplify the procedure and increase its robustness. Although significant achievements in this area have been made, bisulfite treatment remains the main source of process variability in the analysis of DNA methylation. This variability in particular impairs assays, which strive for the quantitative assessment of DNA methylation. Here we present basic mathematical considerations, which should be taken into account when analyzing DNA methylation. We also introduce a PCR-based assay, which allows ab initio assessment of the DNA quality after bisulfite treatment and can help to prevent inaccurate quantitative measurement resulting from poor bisulfite treatment.

     

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.     

Methylation status of Sillago japonica satellite DNA examined by bisulfite modification

A member of Sillago japonica satellite DNA contained internal subrepeats in its 174 bp unit. S. Japonica genomic DNA isolated from liver tissue was subjected to bisulfite modification, and the DNA sequences of about 40 bp flanked by both subrepeats were amplified by polymerase chain reaction (PCR). This protocol, combination of bisulfite reaction and PCR, converts cytosines in the genomic DNA to thymines in the amplified DNA, whereas 5-methylcytosines in the genomic DNA remain as cytosines. Sequence analysis of the amplified DNA fragments revealed that most of the cytosine residues at CpG were methylated in this region.     

Systematic assessment of reduced representation bisulfite sequencing to human blood samples: A promising method for large-sample-scale epigenomic studies

Complementary to the time- and cost-intensive direct bisulfite sequencing, we applied reduced representation bisulfite sequencing (RRBS) to the human peripheral blood mononuclear cells (PBMC) from YH, the Asian individual whose genome and epigenome has been deciphered in the YH project and systematically assessed the genomic coverage, coverage depth and reproducibility of this technology as well as the concordance of DNA methylation levels measured by RRBS and direct bisulfite sequencing for the detected CpG sites. Our result suggests that RRBS can cover more than half of CpG islands and promoter regions with a good coverage depth and the proportion of the CpG sites covered by the biological replicates reaches 80-90%, indicating good reproducibility. Given a smaller data quantity, RRBS enjoys much better coverage depth than direct bisulfite sequencing and the concordance of DNA methylation levels between the two methods is high. It can be concluded that RRBS is a time and cost-effective sequencing method for unbiased DNA methylation profiling of CpG islands and promoter regions in a genome-wide scale and it is the method of choice to assay certain genomic regions for multiple samples in a rapid way.     

MethylQuant: a sensitive method for quantifying methylation of specific cytosines within the genome

Here we present MethylQuant, a novel method that allows accurate quantification of the methylation level of a specific cytosine within a complex genome. This method relies on the well-established treatment of genomic DNA with sodium bisulfite, which converts cytosine into uracil without modifying 5-methyl cytosine. The region of interest is then PCR-amplified and quantification of the methylation status of a specific cytosine is performed by methylation-specific real-time PCR with SYBR Green I using one of the primers whose 3′ end discriminates between the methylation states of this cytosine. The presence of a locked nucleic acid at the 3′ end of the discriminative primer provides the specificity necessary for accurate and sensitive quantification, even when one of the methylation states is present at a level as low as 1% of the overall population. We demonstrate that accurate quantification of the methylation status of specific cytosines can be achieved in biological samples. The method is high-throughput, cost-effective, relatively simple and does not require any specific equipment other than a real-time PCR instrument.     

Rapid and quantitative method of allele-specific DNA methylation analysis

Several biological phenomena depend on differential methylation of chromosomal strands. While understanding the role of these processes requires information on allele-specific methylation, the available methodologies are not quantitative or labor-intensive. We describe a novel, rapid method to quantitate allele-specific DNA methylation based on the combination of bisulfite PCR and Pyrosequencing. In this method, DNA is first treated with sodium bisulfite, which converts cytosine but not 5-methylcytosine to uracil. Genes of interest are subsequently amplified using PCR. Allele-specific methylation can then be determined by pyrosequencing each allele individually using sequencing primers that incorporate single nucleotide polymorphisms (SNPs) that allow differentiation between the two parental alleles. This allele-specific methylation methodology can potentially afford quantitative analyses relevant to the regulation of X chromosome inactivation, allele-specific expression of genes in the immune system, repetitive elements, and genomic imprinting. As an illustration of our new method, we quantitated allele-specific methylation of the differentially methylated region of the H19 gene, which is imprinted. Although we could reliably determine allele-specific methylation with our technique, additional studies will be required to confirm the ability of our assay to measure loss of imprinting.     

CpG_MPs: identification of CpG methylation patterns of genomic regions from high-throughput bisulfite sequencing data

High-throughput bisulfite sequencing is widely used to measure cytosine methylation at single-base resolution in eukaryotes. It permits systems-level analysis of genomic methylation patterns associated with gene expression and chromatin structure. However, methods for large-scale identification of methylation patterns from bisulfite sequencing are lacking. We developed a comprehensive tool, CpG_MPs, for identification and analysis of the methylation patterns of genomic regions from bisulfite sequencing data. CpG_MPs first normalizes bisulfite sequencing reads into methylation level of CpGs. Then it identifies unmethylated and methylated regions using the methylation status of neighboring CpGs by hotspot extension algorithm without knowledge of pre-defined regions. Furthermore, the conservatively and differentially methylated regions across paired or multiple samples (cells or tissues) are identified by combining a combinatorial algorithm with Shannon entropy. CpG_MPs identified large amounts of genomic regions with different methylation patterns across five human bisulfite sequencing data during cellular differentiation. Different sequence features and significantly cell-specific methylation patterns were observed. These potentially functional regions form candidate regions for functional analysis of DNA methylation during cellular differentiation. CpG_MPs is the first user-friendly tool for identifying methylation patterns of genomic regions from bisulfite sequencing data, permitting further investigation of the biological functions of genome-scale methylation patterns.     

Quantitative DNA methylation analysis based on four-dye trace data from direct sequencing of PCR amplificates

MOTIVATION: Methylation of cytosines in DNA plays an important role in the regulation of gene expression, and the analysis of methylation patterns is fundamental for the understanding of cell differentiation, aging processes, diseases and cancer development. Such analysis has been limited, because technologies for detailed and efficient high-throughput studies have not been available. We have developed a novel quantitative methylation analysis algorithm and workflow based on direct DNA sequencing of PCR products from bisulfite-treated DNA with high-throughput sequencing machines. This technology is a prerequisite for success of the Human Epigenome Project, the first large genome-wide sequencing study for DNA methylation in many different tissues. Methylation in tissue samples which are compositions of different cells is a quantitative information represented by cytosine/thymine proportions after bisulfite conversion of unmethylated cytosines to uracil and PCR. Calculation of quantitative methylation information from base proportions represented by different dye signals in four-dye sequencing trace files needs a specific algorithm handling imbalanced and overscaled signals, incomplete conversion, quality problems and basecaller artifacts. RESULTS: The algorithm we developed has several key properties: it analyzes trace files from PCR products of bisulfite-treated DNA sequenced directly on ABI machines; it yields quantitative methylation measurements for individual cytosine positions after alignment with genomic reference sequences, signal normalization and estimation of effectiveness of bisulfite treatment; it works in a fully automated pipeline including data quality monitoring; it is efficient and avoids the usual cost of multiple sequencing runs on subclones to estimate DNA methylation. The power of our new algorithm is demonstrated with data from two test systems based on mixtures with known base compositions and defined methylation. In addition, the applicability is proven by identifying CpGs that are differentially methylated in real tissue samples.     

Rapid analysis of DNA methylation using new restriction enzyme sites created by bisulfite modification.

Bisulfite converts non-methylated cytosine in DNA to uracil leaving 5-methylcytosine unaltered. Here, predicted changes in restriction enzyme sites following reaction of genomic DNA with bisulfite and amplification of the product by the polymerase chain reaction (PCR) were used to assess the methylation of CpG sites. This procedure differs from conventional DNA methylation analysis by methylation-sensitive restriction enzymes because it does not rely on an absence of cleavage to detect methylated sites, the two strands of DNA produce different restriction enzyme sites and may be differentially analyzed, and closely related sequences may be separately analyzed by using specific PCR primers.