MethyScan：一种甲基化特异性PCR引物设计及评估工具

DNA甲基化是重要的表观遗传现象,对基因表达发挥重要调控功能.大量研究表明,基因DNA甲基化是重要的临床诊断生物标志物.在临床上,实施快速、准确的DNA甲基化状态检测是诊断应用的前提和关键.甲基化特异性PCR(methylation specific PCR,MSP)通过将两种引物与甲基化、非甲基化模板各自特异性结合和扩增,实现基因甲基化状态的区分,是切实可行、简单便捷的临床诊断实验技术.但是,不同于常规PCR,MSP主要存在如何强化引物-甲基化/非甲基化模板特异性结合、降低引物序列Tm值差异、去除假阳性扩增及提高敏感性等四大难点.尽管大多数MSP引物设计软件对上述难题都提出了各自解决办法,但在引物设计影响因素考虑、设计与评估并行处理及特异性扩增预测等方面仍然存在较大缺陷.为此,本研究通过对MethPrimer、MSPPrimer、MethBlast、BiSearch等现有MSP引物设计软件原理的深入探究,以及对Bowtie、SAMtools和BEDTools等工具的有效综合整合,基于图形库Matplotlib和第三方Python功能库BioPython与Primer3-py实现了具有系列优点的甲基化特异性PCR引物设计与评估可视化工具MethyScan.它具有引物设计、基因组索引、引物评估等三大完整功能模块,不仅可快速进行MSP引物设计,实现巢式(Nested)引物适配,还可基于4种基因组碱基转换模板分析引物结合信息,图形化展示非特异性扩增与目的片段差异,从而综合评估引物特异性-非特异性扩增.同时,对食管癌、结直肠癌等多种恶性肿瘤中6个潜在生物标志物TFPI-2、NDRG4、CDKN2A、CD44、CASP8和SDHD的甲基化引物设计对比结果表明,MethyScan不仅可获得更多CpG位点的检测引物,而且所获得MSP引物位置与其他软件结果相同或相近,且引物间Tm值差值更小.总之,作为首个图形化展示特异性-非特异性扩增差异MSP引物设计工具,MethyScan可有效提高甲基化引物设计准确性,为临床DNA甲基化检测项目开展、检测试验实施及诊断试剂盒研发提供有力支撑.MethyScan工具下载地址:https://github.com/bioinfo-ibms-pumc/MethyScan.     

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.     

Single-molecule polymerase chain reaction reduces bias: application to DNA methylation analysis by bisulfite sequencing

The treatment of DNA with bisulfite, which converts C to U but leaves 5-methyl-C unchanged, forms the basis of many analytical techniques for DNA methylation analysis. Many techniques exist for measuring the methylation state of a single CpG but, for analysis of an entire region, cloning and sequencing remains the gold standard. However, biases in polymerase chain reaction (PCR) amplification and in cloning can skew the results. We hypothesized that single-molecule PCR (smPCR) amplification would eliminate the PCR amplification bias because competition between templates that amplify at different efficiencies no longer exists. The amplified products can be sequenced directly, thus eliminating cloning bias. We demonstrated this accurate and unbiased approach by analyzing a sample that was expected to contain a 50:50 ratio of methylated to unmethylated molecules: a region of the X-linked FMR1 gene from a human female cell line. We compared traditional cloning and sequencing to smPCR and sequencing. Sequencing smPCR products gave an expected methylated to unmethylated ratio of 48:52, whereas conventional cloning and sequencing gave a biased ratio of 72:28. Our results show that smPCR sequencing can eliminate both PCR and cloning bias and represents an attractive approach to bisulfite sequencing.     

Detection and measurement of PCR bias in quantitative methylation analysis of bisulphite-treated DNA.

Methylation analysis of individual cytosines in genomic DNA can be determined quantitatively by bisulphite treatment and PCR amplification of the target DNA sequence, followed by restriction enzyme digestion or sequencing. Methylated and unmethylated molecules, however, have different sequences after bisulphite conversion. For some sequences this can result in bias during the PCR amplification leading to an inaccurate estimate of methylation. PCR bias is sequence dependent and often strand-specific. This study presents a simple method for detection and measurement of PCR bias for any set of primers, and investigates parameters for overcoming PCR bias.     

Headloop suppression PCR and its application to selective amplification of methylated DNA sequences

Selective amplification in PCR is principally determined by the sequence of the primers and the temperature of the annealing step. We have developed a new PCR technique for distinguishing related sequences in which additional selectivity is dependent on sequences within the amplicon. A 5′ extension is included in one (or both) primer(s) that corresponds to sequences within one of the related amplicons. After copying and incorporation into the PCR product this sequence is then able to loop back, anneal to the internal sequences and prime to form a hairpin structure—this structure is then refractory to further amplification. Thus, amplification of sequences containing a perfect match to the 5′ extension is suppressed while amplification of sequences containing mismatches or lacking the sequence is unaffected. We have applied Headloop PCR to DNA that had been bisulphite-treated for the selective amplification of methylated sequences of the human GSTP1 gene in the presence of up to a 10⁵-fold excess of unmethylated sequences. Headloop PCR has a potential for clinical application in the detection of differently methylated DNAs following bisulphite treatment as well as for selective amplification of sequence variants or mutants in the presence of an excess of closely related DNA sequences.     

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.     

Bisulphite Differential Denaturation PCR for Analysis of DNA Methylation

Differential denaturation during PCR can be used to selectively amplify unmethylated DNA from a methylated DNA background. The use of differential denaturation in PCR is particularly suited to amplification of undermethylated sequences following treatment with bisulphite, since bisulphite selectively converts cytosines to uracil while methylated cytosines remain unreactive. Thus amplicons derived from unmethylated DNA retain less cytosines and their lower G + C content allows for their amplification at the lower melting temperatures, while limiting amplification of the corresponding methylated amplicons (Bisulphite Differential Denaturation PCR, BDD-PCR). Selective amplification of unmethylated DNA of four human genomic regions from three genes, GSTP1, BRCA1 and MAGE-A1, is demonstrated with selectivity observed at a ratio of down to one unmethylated molecule in 105 methylated molecules. BDD-PCR has the potential to be used to selectively amplify and detect aberrantly demethylated genes, such as oncogenes, in cancers. Additionally BDD-PCR can be effectively utilised in improving the specificity of methylation specific PCR (MSP) by limiting amplification of DNA that is not fully converted, thus preventing misinterpretation of the methylation versus non-conversion.      

A novel real-time PCR assay for quantitative analysis of methylated alleles (QAMA): analysis of the retinoblastoma locus

Altered methylation patterns have been found to play a role in developmental disorders, cancer and aging. Increasingly, changes in DNA methylation are used as molecular markers of disease. Therefore, there is a need for reliable and easy to use techniques to detect and measure DNA methylation in research and routine diagnostics. We have established a novel quantitative analysis of methylated alleles (QAMA) which is essentially a major improvement over a previous method based on real-time PCR (MethyLight). This method is based on real-time PCR on bisulfite-treated DNA. A significant advantage over conventional MethyLight is gained by the use of TaqMan probes based on minor groove binder (MGB) technology. Their improved sequence specificity facilitates relative quantification of methylated and unmethylated alleles that are simultaneously amplified in single tube. This improvement allows precise measurement of the ratio of methylated versus unmethylated alleles and cuts down potential sources of inter-assay variation. Therefore, fewer control assays are required. We have used this novel technical approach to identify hypermethylation of the CpG island located in the promoter region of the retinoblastoma (RB1) gene and found that QAMA facilitates reliable and fast measurement of the relative quantity of methylated alleles and improves handling of diagnostic methylation analysis. Moreover, the simplified reaction setup and robustness inherent to the single tube assay facilitates high-throughput methylation analysis. Because the high sequence specificity inherent to the MGB technology is widely used to discriminate single nucleotide polymorphisms, QAMA potentially can be used to discriminate the methylation status of single CpG dinucleotides.     

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.     

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.     

甲基化特异性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综合征提供了新的方法。     

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.     

Bias in Template-to-Product Ratios in Multitemplate PCR

Bias introduced by the simultaneous amplification of specific genes from complex mixtures of templates remains poorly understood. To explore potential causes and the extent of bias in PCR amplification of 16S ribosomal DNAs (rDNAs), genomic DNAs of two closely and one distantly related bacterial species were mixed and amplified with universal, degenerate primers. Quantification and comparison of template and product ratios showed that there was considerable and reproducible overamplification of specific templates. Variability between replicates also contributed to the observed bias but in a comparatively minor way. Based on these initial observations, template dosage and differences in binding energies of permutations of the degenerate, universal primers were tested as two likely causes of this template-specific bias by using 16S rDNA templates modified by site-directed mutagenesis. When mixtures of mutagenized templates containing AT- and GC-rich priming sites were used, templates containing the GC-rich permutation amplified with higher efficiency, indicating that different primer binding energies may to a large extent be responsible for overamplification. In contrast, gene copy number was found to be an unlikely cause of the observed bias. Similarly, amplification from DNA extracted from a natural community to which different amounts of genomic DNA of a single bacterial species were added did not affect relative product ratios. Bias was reduced considerably by using high template concentrations, by performing fewer cycles, and by mixing replicate reaction preparations.     

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.     

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

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

BSMAP: whole genome bisulfite sequence MAPping program

Background  

Bisulfite sequencing is a powerful technique to study DNA cytosine methylation. Bisulfite treatment followed by PCR amplification specifically converts unmethylated cytosines to thymine. Coupled with next generation sequencing technology, it is able to detect the methylation status of every cytosine in the genome. However, mapping high-throughput bisulfite reads to the reference genome remains a great challenge due to the increased searching space, reduced complexity of bisulfite sequence, asymmetric cytosine to thymine alignments, and multiple CpG heterogeneous methylation.     

Directional cloning of native PCR products with preformed sticky ends (Autosticky PCR)

A novel method for the directional cloning of native PCR products was developed. Abasic sites in DNA templates make DNA polymerases stall at the site during synthesis of the complementary strand. Since the 5′ ends of PCR product strands contain built-in amplification primers, abasic sites within the primers result in the formation of 5′ single-stranded overhangs at the ends of the PCR product, enabling its direct ligation to a suitably cleaved cloning vector without any further modification. This “autosticky PCR” (AS-PCR) overcomes the problems caused by end sensitivity of restriction enzymes, or internal restriction sites within the amplified sequences, and enables the generation of essentially any desired 5′ overhang.