Analysis of repetitive element DNA methylation by MethyLight

Repetitive elements represent a large portion of the human genome and contain much of the CpG methylation found in normal human postnatal somatic tissues. Loss of DNA methylation in these sequences might account for most of the global hypomethylation that characterizes a large percentage of human cancers that have been studied. There is widespread interest in correlating the genomic 5-methylcytosine content with clinical outcome, dietary history, lifestyle, etc. However, a high-throughput, accurate and easily accessible technique that can be applied even to paraffin-embedded tissue DNA is not yet available. Here, we report the development of quantitative MethyLight assays to determine the levels of methylated and unmethylated repeats, namely, Alu and LINE-1 sequences and the centromeric satellite alpha (Satα) and juxtacentromeric satellite 2 (Sat2) DNA sequences. Methylation levels of Alu, Sat2 and LINE-1 repeats were significantly associated with global DNA methylation, as measured by high performance liquid chromatography, and the combined measurements of Alu and Sat2 methylation were highly correlative with global DNA methylation measurements. These MethyLight assays rely only on real-time PCR and provide surrogate markers for global DNA methylation analysis. We also describe a novel design strategy for the development of methylation-independent MethyLight control reactions based on Alu sequences depleted of CpG dinucleotides by evolutionary deamination on one strand. We show that one such Alu-based reaction provides a greatly improved detection of DNA for normalization in MethyLight applications and is less susceptible to normalization errors caused by cancer-associated aneuploidy and copy number changes.     

MethyLight droplet digital PCR for detection and absolute quantification of infrequently methylated alleles

Aberrant DNA methylation is a common epigenetic alteration found in colorectal adenomas and cancers and plays a role in cancer initiation and progression. Aberrantly methylated DNA loci can also be found infrequently present in normal colon tissue, where they seem to have potential to be used as colorectal cancer (CRC) risk biomarkers. However, detection and precise quantification of the infrequent methylation events seen in normal colon is likely beyond the capability of commonly used PCR technologies. To determine the potential for methylated DNA loci as CRC risk biomarkers, we developed MethyLight droplet digital PCR (ddPCR) assays and compared their performance to the widely used conventional MethyLight PCR. Our analyses demonstrated the capacity of MethyLight ddPCR to detect a single methylated NTRK3 allele from among more than 3125 unmethylated alleles, 25-fold more sensitive than conventional MethyLight PCR. The MethyLight ddPCR assay detected as little as 19 and 38 haploid genome equivalents of methylated EVL and methylated NTRK3, respectively, which far exceeded conventional MethyLight PCR (379 haploid genome equivalents for both genes). When assessing methylated EVL levels in CRC tissue samples, MethyLight ddPCR reduced coefficients of variation (CV) to 6–65% of CVs seen with conventional MethyLight PCR. Importantly, we showed the ability of MethyLight ddPCR to detect infrequently methylated EVL alleles in normal colon mucosa samples that could not be detected by conventional MethyLight PCR. This study suggests that the sensitivity and precision of methylation detection by MethyLight ddPCR enhances the potential of methylated alleles for use as CRC risk biomarkers.     

DNA methylation analysis by digital bisulfite genomic sequencing and digital MethyLight

Alterations in cytosine-5 DNA methylation are frequently observed in most types of human cancer. Although assays utilizing PCR amplification of bisulfite-converted DNA are widely employed to analyze these DNA methylation alterations, they are generally limited in throughput capacity, detection sensitivity, and or resolution. Digital PCR, in which a DNA sample is analyzed in distributive fashion over multiple reaction chambers, allows for enumeration of discrete template DNA molecules, as well as sequestration of non-specific primer annealing templates into negative chambers, thereby increasing the signal-to-noise ratio in positive chambers. Here, we have applied digital PCR technology to bisulfite-converted DNA for single-molecule high-resolution DNA methylation analysis and for increased sensitivity DNA methylation detection. We developed digital bisulfite genomic DNA sequencing to efficiently determine single-basepair DNA methylation patterns on single-molecule DNA templates without an interim cloning step. We also developed digital MethyLight, which surpasses traditional MethyLight in detection sensitivity and quantitative accuracy for low quantities of DNA. Using digital MethyLight, we identified single-molecule, cancer-specific DNA hypermethylation events in the CpG islands of RUNX3, CLDN5 and FOXE1 present in plasma samples from breast cancer patients.     

Cluster analysis for DNA methylation profiles having a detection threshold

Background  

DNA methylation, a molecular feature used to investigate tumor heterogeneity, can be measured on many genomic regions using the MethyLight technology. Due to the combination of the underlying biology of DNA methylation and the MethyLight technology, the measurements, while being generated on a continuous scale, have a large number of 0 values. This suggests that conventional clustering methodology may not perform well on this data.     

Quantification of epigenetic biomarkers: an evaluation of established and emerging methods for DNA methylation analysis

Background

DNA methylation is an important epigenetic mechanism in several human diseases, most notably cancer. The quantitative analysis of DNA methylation patterns has the potential to serve as diagnostic and prognostic biomarkers, however, there is currently a lack of consensus regarding the optimal methodologies to quantify methylation status. To address this issue we compared five analytical methods: (i) MethyLight qPCR, (ii) MethyLight digital PCR (dPCR), methylation-sensitive and -dependent restriction enzyme (MSRE/MDRE) digestion followed by (iii) qPCR or (iv) dPCR, and (v) bisulfite amplicon next generation sequencing (NGS). The techniques were evaluated for linearity, accuracy and precision.

Results

MethyLight qPCR displayed the best linearity across the range of tested samples. Observed methylation measured by MethyLight- and MSRE/MDRE-qPCR and -dPCR were not significantly different to expected values whilst bisulfite amplicon NGS analysis over-estimated methylation content. Bisulfite amplicon NGS showed good precision, whilst the lower precision of qPCR and dPCR analysis precluded discrimination of differences of < 25% in methylation status. A novel dPCR MethyLight assay is also described as a potential method for absolute quantification that simultaneously measures both sense and antisense DNA strands following bisulfite treatment.

Conclusions

Our findings comprise a comprehensive benchmark for the quantitative accuracy of key methods for methylation analysis and demonstrate their applicability to the quantification of circulating tumour DNA biomarkers by using sample concentrations that are representative of typical clinical isolates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1174) contains supplementary material, which is available to authorized users.     

A comparison of cluster analysis methods using DNA methylation data

MOTIVATION: Aberrant DNA methylation is common in cancer. DNA methylation profiles differ between tumor types and subtypes and provide a powerful diagnostic tool for identifying clusters of samples and/or genes. DNA methylation data obtained with the quantitative, highly sensitive MethyLight technology is not normally distributed; it frequently contains an excess of zeros. Established tools to analyze this type of data do not exist. Here, we evaluate a variety of methods for cluster analysis to determine which is most reliable. RESULTS: We introduce a Bernoulli-lognormal mixture model for clustering DNA methylation data obtained using MethyLight. We model the outcomes using a two-part distribution having discrete and continuous components. It is compared with standard cluster analysis approaches for continuous data and for discrete data. In a simulation study, we find that the two-part model has the lowest classification error rate for mixture outcome data compared with other approaches. The methods are illustrated using DNA methylation data from a study of lung cancer cell lines. Compared with competing hierarchical clustering methods, the mixture model approaches have the lowest cross-validation error for detecting lung cancer subtype (non-small versus small cell). The Bernoulli-lognormal mixture assigns observations to subgroups with the lowest uncertainty. AVAILABILITY: Software is available upon request from the authors. SUPPLEMENTARY INFORMATION: http://www-rcf.usc.edu/~kims/SupplementaryInfo.html     

Correlations in global DNA methylation measures in peripheral blood mononuclear cells and granulocytes

Alterations in global DNA methylation levels have been associated with chronic diseases. Despite the increase in the number of studies measuring markers of global methylation, few have adequately examined within-individual differences by source of DNA and whether within-individual differences by source of DNA differ by age, race and other lifestyle factors. We examined correlations between peripheral mononuclear cell (PBMC) and granulocyte DNA methylation levels measured by the luminometric methylation assay (LUMA), and in LINE-1, Sat2, and Alu by MethyLight and pyrosequencing, in the same individual in 112 women participating in The New York City Multiethnic Breast Cancer Project. Levels of DNA methylation of Sat2 by MethyLight (r = 0.57; P < 0.01) and LINE-1 by pyrosequencing (r = 0.30; P < 0.01) were correlated between PBMC and granulocyte DNA of the same individuals, but LUMA and Alu levels were not. The magnitude of the correlations for Sat2 and LINE-1 varied when stratified by selected demographic and lifestyle factors, although the study sample size limited our comparisons across subgroups. These results lend further support to the importance of considering the source of DNA in epidemiologic studies of white blood cell DNA methylation. Results from studies that combine individuals with different available DNA sources need to be interpreted with caution.     

Correlations in global DNA methylation measures in peripheral blood mononuclear cells and granulocytes

Alterations in global DNA methylation levels have been associated with chronic diseases. Despite the increase in the number of studies measuring markers of global methylation, few have adequately examined within-individual differences by source of DNA and whether within-individual differences by source of DNA differ by age, race and other lifestyle factors. We examined correlations between peripheral mononuclear cell (PBMC) and granulocyte DNA methylation levels measured by the luminometric methylation assay (LUMA), and in LINE-1, Sat2, and Alu by MethyLight and pyrosequencing, in the same individual in 112 women participating in The New York City Multiethnic Breast Cancer Project. Levels of DNA methylation of Sat2 by MethyLight (r = 0.57; P < 0.01) and LINE-1 by pyrosequencing (r = 0.30; P < 0.01) were correlated between PBMC and granulocyte DNA of the same individuals, but LUMA and Alu levels were not. The magnitude of the correlations for Sat2 and LINE-1 varied when stratified by selected demographic and lifestyle factors, although the study sample size limited our comparisons across subgroups. These results lend further support to the importance of considering the source of DNA in epidemiologic studies of white blood cell DNA methylation. Results from studies that combine individuals with different available DNA sources need to be interpreted with caution.     

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.     

Physical activity and global genomic DNA methylation in a cancer-free population

Changes in DNA methylation may represent an intermediate step between the environment and human diseases. Little is known on whether behavioral risk factors may modify gene expression through DNA methylation. To assess whether DNA methylation is associated with different levels of physical activity, we measured global genomic DNA methylation using bisulfite converted DNA and real time PCR (MethyLight) for LINE-1 in peripheral blood of 161 participants aged 45-75 years enrolled in the North Texas Healthy Heart Study and levels of physical activity using an accelerometer (Actigraph GT1M Monitor). We found that individuals with physical activity 26-30 min/day had a significantly higher level of global genomic DNA methylation compared to those with physical activity ≤ 10 min/day (β=2.52, 95%CI: 0.70, 4.35) However, the association was attenuated and became statistically insignificant after multivariate adjustment (β=1.24, 95%CI:-0.93, 3.40). There were some suggestions of a positive association between physical activity and global genomic DNA methylation in non-Hispanics (β=1.50, 95%CI: -0.08, 3.08) that warrants further investigation.     

DNA methylation: a profile of methods and applications

Ever since methylcytosine was found in genomic DNA, this epigenetic alteration has become a center of scientific attraction, especially because of its relation to gene silencing in disease. There is currently a wide range of methods designed to yield quantitative and qualitative information on genomic DNA methylation. The earliest approaches were concentrated on the study of overall levels of methylcytosine, but more recent efforts havefocused on the study ofthe methylation status of specific DNA sequences. Particularly, optimization of the methods based on bisulfite modification of DNA permits the analysis of limited CpGs in restriction enzyme sites (e.g., combined bisulfite restriction analyses and methylation-sensitive single nucleotide primer extension) and the overall characterization based on differential methylation states (e.g., methylation-specific PCR, MethyLight, and methylation-sensitive single-stranded conformational polymorphism) and allows very specific patterns of methylation to be revealed (bisulfite DNA sequencing). In addition, novel methods designed to search for new methylcytosine hot spots have yielded further data without requiring prior knowledge of the DNA sequence. We hope this review will be a valuable tool in selecting the best techniques to address particular questions concerning the cytosine methylation status of genomic DNA.     

Comparative analysis of methods to determine DNA methylation levels of a tumor-related microRNA gene

Quantifying levels of DNA methylation in tumors is a useful approach for the identification of potential tumor suppressors and to find biomarkers that can be used as prognostic or therapeutic indicators. In the current study, we compared three methods commonly used for quantifying DNA methylation—bisulfite pyrosequencing, quantitative methylation-specific PCR (Q-MSP), and MethyLight—by focusing on the CpG island of the gene encoding the microRNA-34b and microRNA-34c (miR-34b/c); aberrant regulation of this miR is associated with various human malignancies, including gastric cancer. Standard curve analysis using control DNA samples demonstrated the highest quantitative accuracy in Q-MSP analysis. We also carried out methylation analysis using gastric mucosa specimens obtained from gastric cancer patients. We found a high correlation between methylation levels determined by Q-MSP and those determined by MethyLight (R² = 0.952), whereas the results of bisulfite pyrosequencing and the other two methods were less well correlated (R² = 0.864 and R² = 0.804 for Q-MSP and MethyLight, respectively). This may reflect possible PCR bias in the pyrosequencing technique, which we show can be corrected for by applying a cubic approximate equation to the original data. Thus, although results obtained by the different DNA methylation analysis techniques are largely comparable, an appropriate correction may be necessary for stringent comparison.     

Altered LINE-1 Methylation in Mothers of Children with Down Syndrome

Down syndrome (DS, also known as trisomy 21) most often results from chromosomal nondisjunction during oogenesis. Numerous studies sustain a causal link between global DNA hypomethylation and genetic instability. It has been suggested that DNA hypomethylation might affect the structure and dynamics of chromatin regions that are critical for chromosome stability and segregation, thus favouring chromosomal nondisjunction during meiosis. Maternal global DNA hypomethylation has not yet been analyzed as a potential risk factor for chromosome 21 nondisjunction. This study aimed to asses the risk for DS in association with maternal global DNA methylation and the impact of endogenous and exogenous factors that reportedly influence DNA methylation status. Global DNA methylation was analyzed in peripheral blood lymphocytes by quantifying LINE-1 methylation using the MethyLight method. Levels of global DNA methylation were significantly lower among mothers of children with maternally derived trisomy 21 than among control mothers (P = 0.000). The combination of MTHFR C677T genotype and diet significantly influenced global DNA methylation (R² = 4.5%, P = 0.046). The lowest values of global DNA methylation were observed in mothers with MTHFR 677 CT+TT genotype and low dietary folate. Although our findings revealed an association between maternal global DNA hypomethylation and trisomy 21 of maternal origin, further progress and final conclusions regarding the role of global DNA methylation and the occurrence of trisomy 21 are facing major challenges.     

Physical activity and global genomic DNA methylation in a cancer-free population

Changes in DNA methylation may represent an intermediate step between the environment and human diseases. Little is known on whether behavioral risk factors may modify gene expression through DNA methylation. To assess whether DNA methylation is associated with different levels of physical activity, we measured global genomic DNA methylation using bisulfite-converted DNA and real-time PCR (MethyLight) for LINE-1 in peripheral blood of 161 participants aged 45–75 years enrolled in the North Texas Healthy Heart Study and levels of physical activity using an accelerometer (Actigraph GT1M Monitor). We found that individuals with physical activity 26–30 min/day had a significantly higher level of global genomic DNA methylation compared to those with physical activity ≤10 min/day (β = 2.52, 95% CI: 0.70, 4.35). However, the association was attenuated and became statistically insignificant after multivariate adjustment (β = 1.24, 95% CI: −0.93, 3.40). There were some suggestions of a positive association between physical activity and global genomic DNA methylation in non-Hispanics (β = 1.50, 95% CI: −0.08, 3.08) that warrants further investigation.Key words: DNA methylation, physical activity, peripheral blood     

Global DNA methylation levels in girls with and without a family history of breast cancer

Lower levels of global DNA methylation in white blood cell (WBC) DNA have been associated with adult cancers. It is unknown whether individuals with a family history of cancer also have lower levels of global DNA methylation early in life. We examined global DNA methylation in WBC (measured in three repetitive elements, LINE1, Sat2 and Alu, by MethyLight and in LINE1 by pyrosequencing) in 51 girls aged 6–17 years. Compared to girls without a family history of breast cancer, methylation levels were lower for all assays in girls with a family history of breast cancer and statistically significantly lower for Alu and LINE1 pyrosequencing. After adjusting for age, body mass index (BMI) and Tanner stage, only methylation in Alu was associated with family history of breast cancer. If these findings are replicated in larger studies, they suggest that lower levels of global WBC DNA methylation observed later in life in adults with cancer may also be present early in life in children with a family history of cancer.Key words: Alu, DNA global methylation, early life exposure, epigenetics, LINE1, methylight, pyrosequencing, Sat2     

Sensitive measurement of unmethylated repeat DNA sequences by end-specific PCR

We describe a new method that is well-suited for the determination of the methylation level of repetitive sequences such as human Alu elements. We have applied the method to the analysis of cell and tissue DNAs and expect it to have wide utility in studies of DNA methylation in cancer and other disease states, in monitoring response to epigenetic cancer therapies and in epidemiological studies. Only 1 ng DNA is needed for a duplex, one-tube real-time PCR in which methylation level and the amount of input DNA are concurrently measured. The relative cutting by the methylation-sensitive enzyme BstUI is compared with that of the methylation-insensitive enzyme DraI to give a measure of DNA methylation. The method depends upon the use of 5'-tailed, 3'-blocked oligonucleotides called facilitator oligonucleotides (Foligos). Only cut DNAs with specific matching sequences at their 3' ends can copy the tails of the Foligos and thus become tagged and available for subsequent PCR. Both the tagging and PCR are carried out by the same enzyme, Taq DNA polymerase. Because amplification only occurs if suitable ends have been generated in the target DNA, we have called this method end-specific PCR (ESPCR). ESPCR avoids the bisulfite treatment step that is usually required to measure methylation.     

Correlation of DNA methylation levels in blood and saliva DNA in young girls of the LEGACY Girls study

Many epidemiologic studies of environmental exposures and disease susceptibility measure DNA methylation in white blood cells (WBC). Some studies are also starting to use saliva DNA as it is usually more readily available in large epidemiologic studies. However, little is known about the correlation of methylation between WBC and saliva DNA. We examined DNA methylation in three repetitive elements, Sat2, Alu, and LINE-1, and in four CpG sites, including AHRR (cg23576855, cg05575921), cg05951221 at 2q37.1, and cg11924019 at CYP1A1, in 57 girls aged 6–15 years with blood and saliva collected on the same day. We measured all DNA methylation markers by bisulfite-pyrosequencing, except for Sat2 and Alu, which were measured by the MethyLight assay. Methylation levels measured in saliva DNA were lower than those in WBC DNA, with differences ranging from 2.8% for Alu to 14.1% for cg05575921. Methylation levels for the three repetitive elements measured in saliva DNA were all positively correlated with those in WBC DNA. However, there was a wide range in the Spearman correlations, with the smallest correlation found for Alu (0.24) and the strongest correlation found for LINE-1 (0.73). Spearman correlations for cg05575921, cg05951221, and cg11924019 were 0.33, 0.42, and 0.79, respectively. If these findings are replicated in larger studies, they suggest that, for selected methylation markers (e.g., LINE-1), methylation levels may be highly correlated between blood and saliva, while for others methylation markers, the levels may be more tissue specific. Thus, in studies that differ by DNA source, each interrogated site should be separately examined in order to evaluate the correlation in DNA methylation levels across DNA sources.