Primer design versus PCR bias in methylation independent PCR amplifications

Many protocols in methylation studies utilize one primer set to generate a PCR product from bisulfite modified template regardless of its methylation status (methylation independent amplification MIP). However, proportional amplification of methylated and unmethylated alleles is hard to achieve due to PCR bias favoring amplification of unmethylated relatively GC poor sequence. Two primer design systems have been proposed to overcome PCR bias in methylation independent amplifications. The first advises against including any CpG dinucleoteides into the primer sequence (CpG-free primers) and the second, recently published by us, is based on inclusion of a limited number of CpG sites into the primer sequence. Here we used the Methylation Sensitive High Resolution Melting (MS-HRM) technology to investigate the ability of primers designed according to both of the above mentioned primer design systems to proportionally amplify methylated and unmethylated templates. Ten “CpG-free” primer pairs and twenty primers containing limited number of CpGs were tested. In reconstruction experiments the “CpG-free” primers showed primer specific sensitivity and allowed us to detect methylation levels in the range from 5 to 50%. Whereas while using primers containing limited number of CpG sites we were able to consistently detect 1–0.1% methylation levels and effectively control PCR amplification bias. In conclusion, the primers with limited number of CpG sites are able to effectively reverse PCR bias and therefore detect methylated templates with significantly higher sensitivity than CpG free primers.     

A rapid, quantitative, non-radioactive bisulfite-SNuPE- IP RP HPLC assay for methylation analysis at specific CpG sites

The precise mapping and quantification of DNA methylation as an epigenetic parameter during development and in diseased tissues is of great importance for functional genomics. Here we describe a rapid, quantitative method to assess methylation levels at specific CpG sites using PCR products of bisulfite-treated genomic DNA. Using single nucleotide primer extension (SNuPE) assays in combination with ion pair reverse phase high performance liquid chromatography (IP RP HPLC) separation techniques, methylated and unmethylated CpGs can be discriminated and quantified based on the different masses and hydrophobicities of the extended primer products. The assay is linear, highly reproducible and several sites can be measured simultaneously in one reaction. It can be semi-automated and eliminates the need for cloning and sequencing of individual bisulfite PCR products.     

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.     

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.     

Microarray-based methylation analysis using dual-color fluorescence hybridization

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

Characterization of Genome-Methylome Interactions in 22 Nuclear Pedigrees

Genetic polymorphisms can shape the global landscape of DNA methylation, by either changing substrates for DNA methyltransferases or altering the DNA binding affinity of cis-regulatory proteins. The interactions between CpG methylation and genetic polymorphisms have been previously investigated by methylation quantitative trait loci (mQTL) and allele-specific methylation (ASM) analysis. However, it remains unclear whether these approaches can effectively and comprehensively identify all genetic variants that contribute to the inter-individual variation of DNA methylation levels. Here we used three independent approaches to systematically investigate the influence of genetic polymorphisms on variability in DNA methylation by characterizing the methylation state of 96 whole blood samples in 52 parent-child trios from 22 nuclear pedigrees. We performed targeted bisulfite sequencing with padlock probes to quantify the absolute DNA methylation levels at a set of 411,800 CpG sites in the human genome. With mid-parent offspring analysis (MPO), we identified 10,593 CpG sites that exhibited heritable methylation patterns, among which 70.1% were SNPs directly present in methylated CpG dinucleotides. We determined the mQTL analysis identified 49.9% of heritable CpG sites for which regulation occurred in a distal cis-regulatory manner, and that ASM analysis was only able to identify 5%. Finally, we identified hundreds of clusters in the human genome for which the degree of variation of CpG methylation, as opposed to whether or not CpG sites were methylated, was associated with genetic polymorphisms, supporting a recent hypothesis on the genetic influence of phenotypic plasticity. These results show that cis-regulatory SNPs identified by mQTL do not comprise the full extent of heritable CpG methylation, and that ASM appears overall unreliable. Overall, the extent of genome-methylome interactions is well beyond what is detectible with the commonly used mQTL and ASM approaches, and is likely to include effects on plasticity.     

A simple algorithm for quantifying DNA methylation levels on multiple independent CpG sites in bisulfite genomic sequencing electropherograms

DNA methylation at cytosines is a widely studied epigenetic modification. Methylation is commonly detected using bisulfite modification of DNA followed by PCR and additional techniques such as restriction digestion or sequencing. These additional techniques are either laborious, require specialized equipment, or are not quantitative. Here we describe a simple algorithm that yields quantitative results from analysis of conventional four-dye-trace sequencing. We call this method Mquant and we compare it with the established laboratory method of combined bisulfite restriction assay (COBRA). This analysis of sequencing electropherograms provides a simple, easily applied method to quantify DNA methylation at specific CpG sites.     

Restriction endonuclease isoschizomers ItaI, BsoFI and Fsp4HI are characterised by differences in their sensitivities to CpG methylation.

BsoFI , ItaI and Fsp4HI are isoshizomers of Fnu4HI (5'-GC NGC-3'). Both Fnu4HI and BsoFI have previously been shown to be inhibited by cytosine-specific methylation within the recognition sequence. Fnu4HI is inhibited if either the internal cytosine at position 2 or the external cytosine at position 5 of the restriction sequence is methylated, but the precise nature of the methylation sensitivity of BsoFI is unclear from the literature. The methylation sensitivities of ItaI and Fsp4HI have not previously been reported. By methylating the plasmid pUC18 with M.SssI (a DNA cytosine-5'-methyltransferase with a specificity for CpG), we have determined that ItaI is sensitive only to methylation of internal CpG sites within the restriction sequence. The methylation sensitivity of Fsp4HI is identical to that of Fnu4HI, being inhibited by methylation of either internal CpG sites or overlapping CpG sites. BsoFI , like the other isoschizomers tested, is sensitive to a combination of internal and overlapping CpG methylation. BsoFI is also sensitive to overlapping CpG methylation (in the absence of internal CpG methylation) if CpG overlap with both sides of the recognition sequence. Sites containing one overlapping CpG (in the absence of internal CpG) are cut when methylated but show marked individual variation in their rates of cleavage. Considerable variation in the rate of cleavage by BsoFI is also observed at sites containing only internal methylated CpG. Some sites are cut slowly, whilst others fail to cut even after prolonged incubation with excess of enzyme.     

The survival decrease in gastric cancer is associated with the methylation of B-cell CLL/lymphoma 6 member B promoter

The methylation of B-cell CLL/lymphoma 6 member B (BCL6B) DNA promoter was detected in several malignancies. Here, we quantitatively detect the methylated status of CpG sites of BCL6B DNA promoter of 459 patients with gastric cancer (GC) by using bisulfite gene sequencing. We show that patients with three or more methylated CpG sites in the BCL6B promoter were significantly associated with poor survival. Furthermore, by using the Akaike information criterion value calculation, we show that the methylated count of BCL6B promoter was identified to be the optimal prognostic predictor of GC patients.     

Regionally Specific and Genome-Wide Analyses Conclusively Demonstrate the Absence of CpG Methylation in Human Mitochondrial DNA

Although CpG methylation clearly distributes genome-wide in vertebrate nuclear DNA, the state of methylation in the vertebrate mitochondrial genome has been unclear. Several recent reports using immunoprecipitation, mass spectrometry, and enzyme-linked immunosorbent assay methods concluded that human mitochondrial DNA (mtDNA) has much more than the 2 to 5% CpG methylation previously estimated. However, these methods do not provide information as to the sites or frequency of methylation at each CpG site. Here, we have used the more definitive bisulfite genomic sequencing method to examine CpG methylation in HCT116 human cells and primary human cells to independently answer these two questions. We found no evidence of CpG methylation at a biologically significant level in these regions of the human mitochondrial genome. Furthermore, unbiased next-generation sequencing of sodium bisulfite treated total DNA from HCT116 cells and analysis of genome-wide sodium bisulfite sequencing data sets from several other DNA sources confirmed this absence of CpG methylation in mtDNA. Based on our findings using regionally specific and genome-wide approaches with multiple human cell sources, we can definitively conclude that CpG methylation is absent in mtDNA. It is highly unlikely that CpG methylation plays any role in direct control of mitochondrial function.     

Demethylation of somatic and testis-specific histone H2A and H2B genes in F9 embryonal carcinoma cells.

In contrast to many other genes containing a CpG island, the testis-specific H2B (TH2B) histone gene exhibits tissue-specific methylation patterns in correlation with gene activity. Characterization of the methylation patterns within a 20-kb segment containing the TH2A and TH2B genes in comparison with that in a somatic histone cluster revealed that: (i) the germ cell-specific unmethylated domain of the TH2A and TH2B genes is defined as a small region surrounding the CpG islands of the TH2A and TH2B genes and (ii) somatic histone genes are unmethylated in both liver and germ cells, like other genes containing CpG islands, whereas flanking sequences are methylated. Transfection of in vitro-methylated TH2B, somatic H2B, and mouse metallothionein I constructs into F9 embryonal carcinoma cells revealed that the CpG islands of the TH2A and TH2B genes were demethylated like those of the somatic H2A and H2B genes and the metallothionein I gene. The demethylation of those CpG islands became significantly inefficient at a high number of integrated copies and a high density of methylated CpG dinucleotides. In contrast, three sites in the somatic histone cluster, of which two sites are located in the long terminal repeat of an endogenous retrovirus-like sequence, were efficiently demethylated even at a high copy number and a high density of methylated CpG dinucleotides. These results suggest two possible mechanisms for demethylation in F9 cells and methylation of CpG islands of the TH2A and TH2B genes at the postblastula stage during embryogenesis.     

A Bayesian hierarchical model to detect differentially methylated loci from single nucleotide resolution sequencing data

DNA methylation is an important epigenetic modification that has essential roles in cellular processes including gene regulation, development and disease and is widely dysregulated in most types of cancer. Recent advances in sequencing technology have enabled the measurement of DNA methylation at single nucleotide resolution through methods such as whole-genome bisulfite sequencing and reduced representation bisulfite sequencing. In DNA methylation studies, a key task is to identify differences under distinct biological contexts, for example, between tumor and normal tissue. A challenge in sequencing studies is that the number of biological replicates is often limited by the costs of sequencing. The small number of replicates leads to unstable variance estimation, which can reduce accuracy to detect differentially methylated loci (DML). Here we propose a novel statistical method to detect DML when comparing two treatment groups. The sequencing counts are described by a lognormal-beta-binomial hierarchical model, which provides a basis for information sharing across different CpG sites. A Wald test is developed for hypothesis testing at each CpG site. Simulation results show that the proposed method yields improved DML detection compared to existing methods, particularly when the number of replicates is low. The proposed method is implemented in the Bioconductor package DSS.     

Active mammalian replication origins are associated with a high-density cluster of mCpG dinucleotides.

ori-beta is a well-characterized origin of bidirectional replication (OBR) located approximately 17 kb downstream of the dihydrofolate reductase gene in hamster cell chromosomes. The approximately 2-kb region of ori-beta that exhibits greatest replication initiation activity also contains 12 potential methylation sites in the form of CpG dinucleotides. To ascertain whether DNA methylation might play a role at mammalian replication origins, the methylation status of these sites was examined with bisulfite to chemically distinguish cytosine (C) from 5-methylcytosine (mC). All of the CpGs were methylated, and nine of them were located within 356 bp flanking the minimal OBR, creating a high-density cluster of mCpGs that was approximately 10 times greater than average for human DNA. However, the previously reported densely methylated island in which all cytosines were methylated regardless of their dinucleotide composition was not detected and appeared to be an experimental artifact. A second OBR, located at the 5' end of the RPS14 gene, exhibited a strikingly similar methylation pattern, and the organization of CpG dinucleotides at other mammalian origins revealed the potential for high-density CpG methylation. Moreover, analysis of bromodeoxyuridine-labeled nascent DNA confirmed that active replication origins were methylated. These results suggest that a high-density cluster of mCpG dinucleotides may play a role in either the establishment or the regulation of mammalian replication origins.     

Assessing combined methylation–sensitive high resolution melting and pyrosequencing for the analysis of heterogeneous DNA methylation

Heterogeneous DNA methylation leads to difficulties in accurate detection and quantification of methylation. Methylation-sensitive high resolution melting (MS-HRM) is unique among regularly used methods for DNA methylation analysis in that heterogeneous methylation can be readily identified, although not quantified, by inspection of the melting curves. Bisulfite pyrosequencing has been used to estimate the level of heterogeneous methylation by quantifying methylation levels present at individual CpG dinucleotides. Sequentially combining the two methodologies using MS-HRM to screen the amplification products prior to bisulfite pyrosequencing would be advantageous. This would not only replace the quality control step using agarose gel analysis prior to the pyrosequencing step but would also provide important qualitative information in its own right. We chose to analyze DAPK1 as it is an important tumor suppressor gene frequently heterogeneously methylated in a number of malignancies, including chronic lymphocytic leukemia (CLL). A region of the DAPK1 promoter was analyzed in ten CLL samples by MS-HRM. By using a biotinylated primer, bisulfite pyrosequencing could be used to directly analyze the samples. MS-HRM revealed the presence of various extents of heterogeneous DAPK1 methylation in all CLL samples. Further analysis of the biotinylated MS-HRM products by bisulfite pyrosequencing provided quantitative information for each CpG dinucleotide analyzed, and confirmed the presence of heterogeneous DNA methylation. Whereas each method could be used individually, MS-HRM and bisulfite pyrosequencing provided complementary information for the assessment of heterogeneous methylation.     

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.