Methylation-dependent fragment separation: direct detection of DNA methylation by capillary electrophoresis of PCR products from bisulfite-converted genomic DNA

Fundamental to understanding the role of cytosine (C) methylation in genomic DNA (gDNA) is the need for robust analysis methods to determine the location and degree of this modification. We report a novel method for methylation detection by denaturing capillary electrophoresis (CE) using standard fragment analysis conditions. Bisulfite treatment of gDNA will selectively deaminate C but not 5-methylcytosine (5mC). Amplicons generated from bisulfite-converted gDNA are analyzed immediately after PCR using a 6-carboxy fluorescein (6-FAM) dye-labeled primer. The amplicons from methylated and unmethylated gDNA separate based solely on base composition due to the presence of multiple C versus thymine (T) differences. By direct detection of PCR amplicons following PCR using primers that anneal independent of methylation status, the overall workflow from gDNA sample input to data analysis is relatively simple. Furthermore, the same PCR product is suitable for additional analyses such as direct sequencing, cloning and sequencing, single-base extension, and post-PCR incorporation of a modified dCTP, the latter of which allows resolution of amplicons with as little as a single C/T difference. We show the utility of this novel CE detection assay by analyzing the hypermethylated region of the fragile-X FMR1 locus.     

Capillary electrophoretic analysis of methylation status in CpG-rich regions by single-base extension of primers modified with N6-methoxy-2,6-diaminopurine

Polymerase-mediated single-base extension (SBE) of primers using a fluorescently labeled 2′,3′-dideoxynucleotide triphosphate terminator was originally commercialized as SNaPshot for analysis of single-nucleotide polymorphisms by capillary electrophoresis (CE). Application of this general method to bisulfite-converted/PCR-amplified genomic DNA (gDNA) to analytically infer polymorphic methylation status (i.e., 5-methylcytosine [5mC] vs. cytosine [C]) in CpG-rich regions of gDNA has been noted previously by others to be limited by CE mobility-shifted peaks for SBE products derived from guanine (G)/adenine (A)-mixed-base primers used to hybridize to possible polymorphic sites (i.e., C vs. thymine [T] resulting from 5mC vs. C, respectively). This limitation is precluded in the current study by using novel SNaPshot primers modified with N⁶-methoxy-2,6-diaminopurine (K), which was originally described in 1991 by Brown and Lin as a unique adenine–guanine analog capable of participating in three H-bonds with C or T in DNA. Oligonucleotides modified by K as a bispecific complement for C/T are commercially available or can be readily synthesized, and they may have utility in other assay formats used to analyze CpG methylation status.     

Methylcytosine and Normal Cytosine Deamination by the Foreign DNA Restriction Enzyme APOBEC3A

Multiple studies have indicated that the TET oxidases and, more controversially, the activation-induced cytidine deaminase/APOBEC deaminases have the capacity to convert genomic DNA 5-methylcytosine (MeC) into altered nucleobases that provoke excision repair and culminate in the replacement of the original MeC with a normal cytosine (C). We show that human APOBEC3A (A3A) efficiently deaminates both MeC to thymine (T) and normal C to uracil (U) in single-stranded DNA substrates. In comparison, the related enzyme APOBEC3G (A3G) has undetectable MeC to T activity and 10-fold less C to U activity. Upon 100-fold induction of endogenous A3A by interferon, the MeC status of bulk chromosomal DNA is unaltered, whereas both MeC and C nucleobases in transfected plasmid DNA substrates are highly susceptible to editing. Knockdown experiments show that endogenous A3A is the source of both of these cellular DNA deaminase activities. This is the first evidence for nonchromosomal DNA MeC to T editing in human cells. These biochemical and cellular data combine to suggest a model in which the expanded substrate versatility of A3A may be an evolutionary adaptation that occurred to fortify its innate immune function in foreign DNA clearance by myeloid lineage cell types.     

116 Deamination of both methyl- and normal-cytosine by the foreign DNA restriction enzyme APOBEC3A

Multiple studies have indicated that the TET oxidases and, more controversially, the AID/APOBEC deaminases have the capacity to convert genomic DNA 5-methyl-cytosine (MeC) into altered nucleobases that provoke excision repair and culminate in the replacement of the original MeC with a normal cytosine (C). We show that human APOBEC3A (A3A) efficiently deaminates both MeC to thymine (T) and normal C to uracil (U) in single-stranded DNA substrates. In comparison, the related enzyme APOBEC3G (A3G) has undetectable MeC-to-T activity and 10-fold less C-to-U activity. Upon 100-fold induction of endogenous A3A by interferon, the MeC status of bulk chromosomal DNA is unaltered whereas both MeC and C nucleobases in transfected plasmid DNA substrates are highly susceptible to editing. Knockdown experiments show that endogenous A3A is the source of both of these cellular DNA deaminase activities. This is the first evidence for non-chromosomal DNA MeC-to-T editing in human cells. These biochemical and cellular data combine to suggest a model in which the expanded substrate versatility of A3A may be an evolutionary adaptation that occurred to fortify its innate immune function in foreign DNA clearance by myeloid lineage cell types.     

Characterization of the level,target sites and inheritance of cytosine methylation in tomato nuclear DNA

The tomato nuclear genome was determined to have a G+C content of 37% which is among the lowest reported for any plant species. Non-coding regions have a G+C content even lower (32% average) whereas coding regions are considerably richer in G+C (46%).5-methyl cytosine was the only modified base detected and on average 23% of the cytosine residues are methylated. Immature tissues and protoplasts have significantly lower levels of cytosine methylation (average 20%) than mature tissues (average 25%). Mature pollen has an intermediate level of methylation (22%). Seeds gave the highest value (27%), suggesting de novo methylation after pollination and during seed development.Based on isoschizomer studies we estimate 55% of the CpG target sites (detected by Msp I/Hpa II) and 85% of the CpNpG target sites (detected by Bst NI/Eco RI)are methylated. Unmethylated target sites (both CpG and CpNpG) are not randomly distributed throughout the genome, but frequently occur in clusters. These clusters resemble CpG islands recently reported in maize and tobacco.The low G+C content and high levels of cytosine methylation in tomato may be due to previous transitions of 5mCT. This is supported by the fact that G+C levels are lowest in non-coding portions of the genome in which selection is relaxed and thus transitions are more likely to be tolerated. This hypothesis is also supported by the general deficiency of methylation target sites in the tomato genome, especially in non-coding regions.Using methylation isoschizomers and RFLP analysis we have also determined that polymorphism between plants, for cytosine methylation at allelic sites, is common in tomato. Comparing DNA from two tomato species, 20% of the polymorphisms detected by Bst NI/Eco RII could be attributed to differential methylation at the CpNpG target sites. With Msp I/Hpa II, 50% of the polymorphisms were attributable to methylation (CpG and CpNpG sites). Moreover, these polymorphisms were demonstrated to be inherited in a mendelian fashion and to co-segregate with the methylation target site and thus do not represent variation for transacting factors that might be involved in methylation of DNA. The potential role of heritable methylation polymorphism in evolution of gene regulation and in RFLP studies is discussed.     

Measurement of DNA methylation using stable isotope dilution and gas chromatography-mass spectrometry

A simple, highly selective, and sensitive method has been developed to quantify methylation of DNA extracted from human peripheral blood mononuclear cells. Assay has been performed at nucleobases level. Cytosine and 5-methylcytosine DNA content has been detected by gas chromatography-mass spectrometry using [2-(13)C]cytosine and [2-(13)C]5-methylcytosine as internal standards. The methylation level has been calculated as 5-methylcytosine/total cytosine ratio. The working range selected on calibration curve, obtained by evaluation of standards and matrix-added standards measurements, is suitable for 5 microg DNA analysis. In this range, healthy human DNA methylation percentage is within 5-6%.     

甲基化特异性PCR检测FMR1 和XIST基因甲基化实验方法的建立

建立一种快速、灵敏的检测脆性X智障基因（Fragile X mental retardation, FMR1）和X染色体失活基因（X chromosome inactivation,XIST）甲基化的方法,用亚硫酸氢钠和对苯二酚对基因组DNA进行脱氨基修饰。以修饰后的DNA为模板,用两套不同的引物对：1对甲基化特异性引物和1对非甲基化特异性引物扩增FMR1基因（CGG）n重复序列区、FMR1 和XIST 基因的启动子区。PCR产物进一步克隆、测序。以亚硫酸氢钠和对苯二酚脱氨基修饰后的DNA为模板,进行PCR扩增后的产物与预期基因目的基因片段大小相符合,无非特异性扩增产物。测序结果表明,FMR1、XIST基因中的非甲基化的C碱基转变为U碱基,而CpG岛被甲基化的C碱基不改变。成功地建立了检测FMR1、XIST甲基化的方法,为实验室诊断脆性X综合征提供了新的方法。     

DNA Methylation: Bisulphite Modification and Analysis

Epigenetics describes the heritable changes in gene function that occur independently to the DNA sequence. The molecular basis of epigenetic gene regulation is complex, but essentially involves modifications to the DNA itself or the proteins with which DNA associates. The predominant epigenetic modification of DNA in mammalian genomes is methylation of cytosine nucleotides (5-MeC). DNA methylation provides instruction to gene expression machinery as to where and when the gene should be expressed. The primary target sequence for DNA methylation in mammals is 5''-CpG-3'' dinucleotides (Figure 1). CpG dinucleotides are not uniformly distributed throughout the genome, but are concentrated in regions of repetitive genomic sequences and CpG "islands" commonly associated with gene promoters (Figure 1). DNA methylation patterns are established early in development, modulated during tissue specific differentiation and disrupted in many disease states including cancer. To understand the biological role of DNA methylation and its role in human disease, precise, efficient and reproducible methods are required to detect and quantify individual 5-MeCs.This protocol for bisulphite conversion is the "gold standard" for DNA methylation analysis and facilitates identification and quantification of DNA methylation at single nucleotide resolution. The chemistry of cytosine deamination by sodium bisulphite involves three steps (Figure 2). (1) Sulphonation: The addition of bisulphite to the 5-6 double bond of cytosine (2) Hydrolic Deamination: hydrolytic deamination of the resulting cytosine-bisulphite derivative to give a uracil-bisulphite derivative (3) Alkali Desulphonation: Removal of the sulphonate group by an alkali treatment, to give uracil. Bisulphite preferentially deaminates cytosine to uracil in single stranded DNA, whereas 5-MeC, is refractory to bisulphite-mediated deamination. Upon PCR amplification, uracil is amplified as thymine while 5-MeC residues remain as cytosines, allowing methylated CpGs to be distinguished from unmethylated CpGs by presence of a cytosine "C" versus thymine "T" residue during sequencing.DNA modification by bisulphite conversion is a well-established protocol that can be exploited for many methods of DNA methylation analysis. Since the detection of 5-MeC by bisulphite conversion was first demonstrated by Frommer et al.¹ and Clark et al.², methods based around bisulphite conversion of genomic DNA account for the majority of new data on DNA methylation. Different methods of post PCR analysis may be utilized, depending on the degree of specificity and resolution of methylation required. Cloning and sequencing is still the most readily available method that can give single nucleotide resolution for methylation across the DNA molecule.     

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.     

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.     

Disclosing Bias in Bisulfite Assay: MethPrimers Underestimate High DNA Methylation

Discordant results obtained in bisulfite assays using MethPrimers (PCR primers designed using MethPrimer software or assuming that non-CpGs cytosines are non methylated) versus primers insensitive to cytosine methylation lead us to hypothesize a technical bias. We therefore used the two kinds of primers to study different experimental models and methylation statuses. We demonstrated that MethPrimers negatively select hypermethylated DNA sequences in the PCR step of the bisulfite assay, resulting in CpG methylation underestimation and non-CpG methylation masking, failing to evidence differential methylation statuses. We also describe the characteristics of “Methylation-Insensitive Primers” (MIPs), having degenerated bases (G/A) to cope with the uncertain C/U conversion. As CpG and non-CpG DNA methylation patterns are largely variable depending on the species, developmental stage, tissue and cell type, a variable extent of the bias is expected. The more the methylome is methylated, the greater is the extent of the bias, with a prevalent effect of non-CpG methylation. These findings suggest a revision of several DNA methylation patterns so far documented and also point out the necessity of applying unbiased analyses to the increasing number of epigenomic studies.     

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.     

Preparation of genome-wide DNA fragment libraries using bisulfite in polyacrylamide gel electrophoresis slices with formamide denaturation and quality control for massively parallel sequencing by oligonucleotide ligation and detection

Bisulfite sequencing is widely used for analysis of DNA methylation status (i.e., 5-methylcytosine [5mC] vs. cytosine [C]) in CpG-rich or other loci in genomic DNA (gDNA). Such methods typically involve reaction of gDNA with bisulfite followed by polymerase chain reaction (PCR) amplification of specific regions of interest that, overall, converts C→T (thymine) and 5mC→C and then capillary sequencing to measure C versus T composition at CpG sites. Massively parallel sequencing by oligonucleotide ligation and detection (SOLiD) has recently enabled relatively low-cost whole genome sequencing, and it would be highly desirable to apply such massively parallel sequencing to bisulfite-converted whole genomes to determine DNA methylation status of an entire genome, which has heretofore not been reported. As an initial step toward achieving this goal, we have extended our ongoing interest in improving bisulfite conversion sample preparation to include a human genome-wide fragment library for SOliD. The current article features novel use of formamide denaturant during bisulfite conversion of a suitably constructed library directly in a band slice from polyacryamide gel electrophoresis (PAGE). To validate this new protocol for 5mC-protected fragment library conversion, which we refer to as Bis-PAGE, capillary-based size analysis and Sanger sequencing were carried out for individual amplicons derived from single-molecule PCR (smPCR) of randomly selected library fragments. smPCR/Capillary Sanger sequencing of approximately 200 amplicons unambiguously demonstrated greater than 99% C→T conversion. All of these approximately 200 Sanger sequences were analyzed with a previously published web-accessible bioinformatics tool (methBLAST) for mapping to human chromosomes, the results of which indicated random distribution of analyzed fragments across all chromosomes. Although these particular Bis-PAGE conversion and quality control methods were exemplified in the context of a fragment library for SOLiD, the concepts can be generalized to include other genome-wide library constructions intended for DNA methylation analysis by alternative high-throughput or massively parallelized methods that are currently available.     

Roles,and establishment,maintenance and erasing of the epigenetic cytosine methylation marks in plants

Heritable information in plants consists of genomic information in DNA sequence and epigenetic information superimposed on DNA sequence. The latter is in the form of cytosine methylation at CG, CHG and CHH elements (where H = A, T or C) and a variety of histone modifications in nucleosomes. The epialleles arising from cytosine methylation marks on the nuclear genomic loci have better heritability than the epiallelic variation due to chromatin marks. Phenotypic variation is increased manifold by epiallele comprised methylomes. Plants (angiosperms) have highly conserved genetic mechanisms to establish, maintain or erase cytosine methylation from epialleles. The methylation marks in plants fluctuate according to the cell/tissue/organ in the vegetative and reproductive phases of plant life cycle. They also change according to environment. Epialleles arise by gain or loss of cytosine methylation marks on genes. The changes occur due to the imperfection of the processes that establish and maintain the marks and on account of spontaneous and stress imposed removal of marks. Cytosine methylation pattern acquired in response to abiotic or biotic stress is often inherited over one to several subsequent generations. Cytosine methylation marks affect physiological functions of plants via their effect(s) on gene expression levels. They also repress transposable elements that are abundantly present in plant genomes. The density of their distribution along chromosome lengths affects meiotic recombination rate, while their removal increases mutation rate. Transposon activation due to loss of methylation causes rearrangements such that new gene regulatory networks arise and genes for microRNAs may originate. Cytosine methylation dynamics contribute to evolutionary changes. This review presents and discusses the available evidence on origin, removal and roles of cytosine methylation and on related processes, such as RNA directed DNA methylation, imprinting, paramutation and transgenerational memory in plants.     

Cytosine methylated DNA synthesized by Taq polymerase used to assay methylation sensitivity of restriction endonuclease HinfI.

We have studied the resistance of cytosine methylated DNA to digestion by the restriction endonuclease HinfI, using a simple PCR procedure to synthesize DNA of known sequence in which every cytosine is methylated at the 5 position. We find that HinfI cannot digest cytosine methylated DNA at the concentrations normally used in restriction digests. Complete digestion is possible using a vast excess of enzyme; under these conditions, the rate of HinfI digestion for cytosine methylated DNA is at least 1440-fold slower than for unmethylated DNA. The presence of an additional methylated cytosine at the degenerate position internal to the recognition sequence does not appear to increase the resistance to HinfI digestion. We also tested HhaII, an isoschizomer of HinfI, and found that it is completely inactive on cytosine methylated DNA. The procedure we have used should be of general applicability in determination of the methylation sensitivities of other restiction enzymes, as well as studies of the effects of methylation on gene expression in direct DNA transfer experiments.     

Profiling DNA methylation by melting analysis

The idea of modifying DNA with bisulfite has paved the way for a variety of polymerase chain reaction (PCR) methods for accurately mapping 5-methylcytosine at specific genes. Bisulfite selectively deaminates cytosine to uracil under conditions where 5-methylcytosine remains unreacted. Following conventional PCR amplification of bisulfite-treated DNA, original cytosines appear as thymine while 5-methylcytosines appear as cytosine. Because the relative thermostability of a DNA duplex increases with increasing content of G:C base pairs, PCR products originating from DNA templates with different contents of 5-methylcytosine differ in melting temperature, i.e., the temperature required to convert the double helix into random coils. We describe two methods that resolve differentially methylated DNA sequences on the basis of differences in melting temperature. The first method integrates PCR amplification of bisulfite-treated DNA and subsequent melting analysis by using a thermal cycler coupled with a fluorometer. By including in the reaction a PCR-compatible, fluorescent dye that specifically binds to double-stranded DNA, the melting properties of the PCR product can be examined directly in the PCR tube by continuous fluorescence monitoring during a temperature transition. The second method relies on resolution of alleles with different 5-methylcytosine contents by analysis of PCR products in a polyacrylamide gel containing a gradient of chemical denaturants. Optimal resolution of differences in melting temperature is achieved by a special design of PCR primers. Both methods allow resolution of "heterogeneous" methylation, i.e., the situation where the content and distribution of 5-methylcytosine in a target gene differ between different molecules in the same sample.     

DNA methylation in states of cell physiology and pathology

DNA methylation is one of epigenetic mechanisms regulating gene expression. The methylation pattern is determined during embryogenesis and passed over to differentiating cells and tissues. In a normal cell, a significant degree of methylation is characteristic for extragenic DNA (cytosine within the CG dinucleotide) while CpG islands located in gene promoters are unmethylated, except for inactive genes of the X chromosome and the genes subjected to genomic imprinting. The changes in the methylation pattern, which may appear as the organism age and in early stages of cancerogenesis, may lead to the silencing of over ninety endogenic genes. It has been found, that these disorders consist not only of the methylation of CpG islands, which are normally unmethylated, but also of the methylation of other dinucleotides, e.g. CpA. Such methylation has been observed in non-small cell lung cancer, in three regions of the exon 5 of the p53 gene (so-called "non-CpG" methylation). The knowledge of a normal methylation process and its aberrations appeared to be useful while searching for new markers enabling an early detection of cancer. With the application of the Real-Time PCR technique (using primers for methylated and unmethylated sequences) five new genes which are potential biomarkers of lung cancer have been presented.     

Sequence analysis and editing for bisulphite genomic sequencing projects

Bisulphite genomic sequencing is a widely used technique for detailed analysis of the methylation status of a region of DNA. It relies upon the selective deamination of unmethylated cytosine to uracil after treatment with sodium bisulphite, usually followed by PCR amplification of the chosen target region. Since this two-step procedure replaces all unmethylated cytosine bases with thymine, PCR products derived from unmethylated templates contain only three types of nucleotide, in unequal proportions. This can create a number of technical difficulties (e.g. for some base-calling methods) and impedes manual analysis of sequencing results (since the long runs of T or A residues are difficult to align visually with the parent sequence). To facilitate the detailed analysis of bisulphite PCR products (particularly using multiple cloned templates), we have developed a visually intuitive program that identifies the methylation status of CpG dinucleotides by analysis of raw sequence data files produced by MegaBace or ABI sequencers as well as Staden SCF trace files and plain text files. The program then also collates and presents data derived from independent templates (e.g. separate clones). This results in a considerable reduction in the time required for completion of a detailed genomic methylation project.