| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |






||**








From the Departments of Pathology,
* Medicine,
and Radiation Oncology,
|| Brigham and Womens Hospital, Boston, Massachusetts; the Departments of Medical Oncology
and Radiation Oncology,
** Dana-Farber Cancer Institute, Boston, Massachusetts; Harvard Medical School,
Boston, Massachusetts; the Departments of Epidemiology and Biostatistics,
¶ Harvard School of Public Health, Boston, Massachusetts; and the Departments of Surgery and of Biochemistry and Molecular Biology,

University of Southern California Keck School of Medicine, Norris Comprehensive Cancer Center, Los Angeles, California
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
A number of methods to determine DNA methylation status in tumor tissues have been developed. For formalin-fixed, paraffin-embedded tissue, methylation-specific polymerase chain reaction (MSP) after sodium bisulfite conversion13
is widely used. However, because of the qualitative nature of the assay, MSP cannot reliably distinguish low levels of methylation from high levels of methylation. Furthermore, MSP does not appear to be highly reproducible for some samples that may have low levels of methylation, and it is difficult to assess performance characteristics. Thus, clinical use of qualitative MSP is a major concern in terms of quality control and assurance. Quantitative measurement of methylation is important because low levels of methylation (below the threshold of transcriptional silencing) may not be biologically important. Hence, a variety of quantitative assays to measure DNA methylation have been developed, including combined bisulfite restriction analysis,14
restriction ligation-mediated polymerase chain reaction (PCR),15
methylation-sensitive single nucleotide primer extension (Ms-SNuPE),16
ion-pair reverse-phase high performance liquid chromatography,17
denaturing high performance liquid chromatography,18
pyrosequencing,19, 20, 21
MALDI-TOF,22
and real-time PCR.5, 23, 24, 25, 26, 27, 28, 29, 30, 31
Most of these methods use DNA treated with sodium bisulfite. Therefore, sodium bisulfite treatment is a critical step in the measurement of DNA methylation. However, reproducibility of the sodium bisulfite conversion step has not been well investigated. Sodium bisulfite treatment presents a harsh environment for DNA molecules, and previous data suggest that
84 to 96% of DNA is degraded during 4 hours of bisulfite treatment at 55°C.32
Assurance of reproducibility in this critical bisulfite treatment step is essential in applying quantitative methylation assays to clinical practice.
In this study, we evaluated the precision and performance characteristics of sodium bisulfite treatment and MethyLight, a real-time PCR assay.5, 26, 28, 33 Our data indicate that both sodium bisulfite conversion and MethyLight assays have good reproducibility and precision and can be effectively used to quantify DNA methylation in paraffin-embedded tumor tissue.
| Materials and Methods |
|---|
|
|
|---|
Sodium Bisulfite Treatment and DNA Extraction
Hematoxylin and eosin (H&E)-stained slides of the tumors were reviewed, and areas of tumor were marked, to exclude pure normal tissue and enrich tumor DNA. H&E-stained tumor sections (10 µm x 14 sections) were scraped off slides from each of case 1 through case 4 and suspended in 140 µl of tissue lysate solution (100 mmol/L Tris, pH 8, 10 mmol/L ethylenediaminetetraacetic acid, pH 8, 1 mg/ml proteinase K, and 0.05 mg/ml tRNA). The lysate solutions were incubated overnight at 50°C. The lysate was aliquoted into seven tubes, each containing 18 µl of tissue lysate, and stored at 20°C until sodium bisulfite modification was performed. For methylation and immunohistochemistry study (using cohort colorectal cancer cases), H&E-stained tissue sections (depending on tissue and tumor size, in average, large tumor tissue 10 µm x 1 section) from each case were scraped off slides, suspended in 20 µl of the tissue lysate, and incubated at 50°C overnight. The tissue lysate was then stored at 20°C until sodium bisulfite modification was performed.
Sodium metabisulfite (1.9 g) was dissolved in mixture of 3.2 ml of 0.44 mol/L NaOH at 50°C. Then, 0.5 ml of 1 mol/L hydroquinone was added to the dissolved sodium bisulfite mixture. An 18-µl aliquot of the tissue lysate was denatured at 100°C for 10 minutes and chilled on ice. Then, 2 µl of 3 mol/L NaOH was added and incubated at 42°C for 20 minutes. The bisulfite solution (120 µl) was added (total volume of 140 µl) and incubated at 50°C for 15 hours in the dark. The bisulfite-converted DNA was recovered using a Qiagen QIAamp viral RNA mini kit (Qiagen, Valencia, CA) according to the manufacturers instructions with some modifications. Buffer AVL/carrier RNA (560 µl) was added to the 140 µl of bisulfite-converted DNA sample and incubated at room temperature for 10 minutes. Ethanol (560 µl) was then added, and after extensive mixing, the mixture was loaded onto the provided spin columns in consecutive 630-µl aliquots. After each loading, the columns were centrifuged at full speed (21,000 x g) for 1 minute. Both the filtrate and spin column were saved, and both filtrates were passed through the column a second time in the same manner to increase the yield of recovery. The spin column was then washed with 500 µl of buffer AW1, followed by centrifugation at 21,000 x g for 1 minute. Buffer AW2 (500 µl) was then added to the column, and the column was centrifuged at 21,000 x g for 4 minutes to wash the column and eliminate possible buffer AW2 carry over. DNA in the spin column was eluted by the addition of 40 µl of buffer AVE, followed by a 1-minute incubation at room temperature and centrifugation at 7600 x g for 1 minute. This elution step was repeated with a second 40-µl volume of buffer AVE. Fifty µl of 0.2 mol/L NaOH was added to the 80-µl pooled eluate for 15 minutes at room temperature to desulphonate the sample, and then 10 µl of 1 mol/L HCl was added to for neutralization. Buffer AVL/carrier RNA (560 µl) was then added to the 140-µl sample mixture, and the recovery procedure was repeated with a new spin column. The eluted DNA (80-µl volume) was then used for MethyLight analysis.
Quantitative Real-Time PCR (MethyLight)
Real-time PCR assays to measure DNA methylation (MethyLight) have been described.5, 26, 28
Briefly, three sets of primers and probes designed specifically for bisulfite-converted DNA were used: a set for the gene of interest and two sets for ACTB and COL2A1 to normalize for the amount of input DNA. We used an ABI 7300 real-time PCR instrument (Applied Biosystems, Foster City, CA). Primers and probes for CDKN2A, MGMT, ACTB, and COL2A1 were previously described.5
The MLH1 forward primer is 5'-AGG AAG AGC GGA TAG CGA TTT-3', the MLH1 reverse primer is 5'-TCT TCG TCC CTC CCT AAA ACG-3', and the MLH1 probe is 6FAM-5'-CCC GCT ACC TAA AAA AAT ATA CGC TTA CGC G-3'-BHQ1. The amount of methylated DNA (PMR, percentage of methylated reference36
) at a specific locus was calculated by dividing the GENE:ACTB (or COL2A1) ratio of a sample by the GENE:ACTB (or COL2A1) ratio of SssI-treated human genomic DNA (presumably fully methylated) and multiplying by 100. Reactions using SssI-treated DNA were used to normalize for any difference in amplification efficiencies between GENE and ACTB (or COL2A1).
Measuring the Precision of Sodium Bisulfite Conversion and MethyLight Assay
Figure 1
shows our strategy to measure precision and reproducibility of bisulfite conversion of DNA and quantitative MethyLight assay. We performed bisulfite conversion on seven different aliquots from each of the four cell lysate samples (case 1 through case 4), starting on 7 different days. Thus, from each cell lysate, we have seven separate bisulfite-converted DNA samples (designated as B1 through B7), which ideally should give at least similar values for methylation at a specific locus. We measured PMR for each locus by the MethyLight assay on each of these seven separate bisulfite-DNA samples. For each case, we measured the coefficient of variation (CV) of these seven PMR values on B1 through B7, which would primarily depend on variations in bisulfite conversion. We repeated the MethyLight assays five times (designated as M1 through M5) and determined PMR values on the exact same set of bisulfite-converted DNA samples (B1 through B7) on 5 different days. We calculated the CV of these five PMR values (M1 through M5) on each one of the seven bisulfite-converted DNA samples (B1 through B7). Thus, these seven CVs would primarily depend on day-to-day variations of each MethyLight assay.
|
DNA Mixing Study to Assess Assay Linearity
We performed a mixing study to assess linearity of sodium bisulfite treatment and quantitative MethyLight assays (Figure 2)
. We first treated human peripheral blood leukocyte DNA (Promega, Madison, WI) with SssI to methylate all CpG sites. Assuming near complete methylation and no loss of DNA during SssI treatment, we prepared DNA mixtures in triplicates (each containing 1 µg of template DNA) for SssI-treated DNA:untreated DNA as follows: 50:50, 30:70, 20:80, 10:90, and 5:95 (Figure 2)
. We used 100% SssI-treated DNA as control reference DNA to calculate PMR as described above. We treated DNA mixtures in triplicates with sodium bisulfite and duplicated MethyLight assays, resulting in six data points of PMR values for each dilution. Linear regression analysis was performed to assess assay linearity.
|
Statistical Analysis
For statistical analysis, paired t-test to compare means, F-test to assess the equality of variances38
and linear regression analysis to assess assay linearity were performed using the Microsoft Excel 2000 Analysis ToolPak. For F-test, raw PMR values were log-transformed to make distributions close to normal. Interclass correlation coefficients were calculated to assess assay reproducibility using the SAS program (version 9.1; SAS Institute, Cary, NC). To assess the association between promoter methylation and loss of protein expression, Fishers exact test was performed using the SAS program. All P values were two-sided.
| Results |
|---|
|
|
|---|
|
|
4 indicated methylation-negative.33
In our analysis, negative markers typically showed PMR
1 and were most frequently undetected (Ct >45 and PMR = 0). There were a total of five case-negative marker combinations (case 1, CDKN2A; case 1, MGMT; case 2, MLH1; case 3, MLH1; and case 4, MGMT). Within a total of 175 PMR results (five case-marker combinations times seven bisulfite-treated samples times five MethyLight runs) for these negative reactions, when using COL2A1 as a control, 160 runs (91%) showed PMR values of 0, 14 runs (8%) showed a PMR <1, and only one run (0.6%) showed PMR >1 (PMR 1.01). However, when using ACTB as a control, 4 of 175 negative reactions (2.3%) showed PMR >1 (PMR ranging from 1.10 to 3.59). Notably with either reference, none of these 175 presumed negative reactions showed PMR >4.
Mean and variance of PMR values in MethyLight using ACTB as a control were compared to those using COL2A1 as a control (Table 4)
. For each case-marker combination, we obtained PMR values 35 times (seven different bisulfite-treated samples times five MethyLight runs) using each one of ACTB and COL2A1 as a control. Distributions of PMR values became closer to normal distributions after log transformation of raw PMR values. Thereafter, we assessed the equality of variances in the 35 repeated runs. PMR variances using ACTB were significantly larger than those using COL2A1 in four of seven case-marker combinations (P < 0.003) (Table 4)
, indicating COL2A1 as a superior control. CVs of PMR values in MethyLight using COL2A1 as a control are shown in Table 5
. Both bisulfite-to-bisulfite variation [CV(B)] and run-to-run variation [CV(M)] were similar, ranging from 0.05 to 0.6. Therefore, MethyLight assay showed fair reproducibility in terms of both bisulfite-to-bisulfite variations [CV(B)] and run-to-run variations [CV(M)]. To assess within-sample variation of PMR values, interclass correlation coefficients were calculated as 0.84, 0.74, and 0.79 for CDKN2A, MLH1, and MGMT, respectively, indicating fair reproducibility of MethyLight assays.
|
|
Assay Linearity
We performed repeated measurements of PMR values on DNA mixtures containing 50%, 30%, 20%, 10%, and 5% of methylated (M-SssI-treated) DNA (Figures 2
and 3)
. Correlation coefficients and P values for the test of a linear association were as follows: r = 0.962, P = 3 x 1017 for CDKN2A; r = 0.913, P = 2 x 1012 for MLH1; and r = 0.904, P = 109 for MGMT. Thus, our MethyLight assay showed fair to good linearity for all three genes.
|
4.33
The bimodal distributions of PMR values in Table 6
|
4 for the promoters of CDKN2A, MLH1, and MGMT, respectively. Loss of protein expression was observed in 62% (45 of 73), 96% (43 of 45), and 66% (69 of 104) of samples with PMR >4 for CDKN2A, MLH1, and MGMT, respectively. Using a PMR cutoff of 4, promoter methylation was significantly associated with loss of respective protein expression (P = 1 x 1021 for CDKN2A, P = 6 x 1044 for MLH1, and P = 1 x 1017 for MGMT). | Discussion |
|---|
|
|
|---|
In this study, we evaluated precision of sodium bisulfite conversion as well as subsequent real-time PCR (MethyLight) assays for the CDKN2A (p16), MLH1, and MGMT promoters. Advantages of MethyLight technology include its quantitative and high-throughput nature and relatively simple assay procedures, which do not require opening of tubes containing PCR products. All of these features are extremely beneficial in clinical molecular diagnostics. Any of these properties is not a feature of qualitative MSP, which has been widely used for methylation detection. As for disadvantages of MethyLight, it cannot quantify methylation at the individual nucleotide level; rather, it can assess methylation levels at the primer sites and/or a probe site as a whole. On the other hand, bisulfite-sequencing, such as bisulfite-pyrosequencing,19, 20, 21 can achieve resolution at the individual nucleotide level. However, in contrast to MethyLight, pyrosequencing involves end-point analysis of PCR products and requires procedures that need opening of PCR tubes.
To measure the relative amount of methylated DNA, we divided the measured amount of a methylated gene of interest by the amount of COL2A1 (or ACTB) for input DNA. One should also be aware of PCR bias in quantitative molecular assays,42 so we further normalized for a potential difference in PCR efficiencies between the gene of interest and COL2A1 by dividing the gene:COL2A1 ratio in a sample of interest by that ratio in SssI-treated DNA (presumably fully methylated). For all three genes tested, we showed good precision in threshold cycle (Ct) values in terms of both bisulfite-to-bisulfite variation (among seven bisulfite-treated aliquots B1 to B7) in the same run and MethyLight-to-MethyLight run-to-run variation (among five repeated MethyLight runs M1 to M5) in the same bisulfite-treated DNA sample. COL2A1 reactions, as well as PMR values based on COL2A1 reactions, showed much smaller variances than ACTB. Thus, COL2A1 is a superior control gene to ACTB in terms of precision of MethyLight assays.
We also demonstrated that run-to-run CVs of PMR values were larger than bisulfite-to-bisulfite CVs of PMR values. This likely reflects the introduction of other sources of variations when we measured run-to-run CVs of PMR, including different standard curves. Thus, we also examined variations of standard curves in a single plate as well as between different MethyLight runs. We demonstrated that variations of standard curves were small and acceptable. Other sources of run-to-run variations are independent control reactions (COL2A1 or ACTB) in different MethyLight runs. A gene that showed small variation in repeated amplification reactions, such as COL2A1, would be suitable for a control. Normalizing methylation measurement for input DNA by a control reaction is necessary to compare quantitative methylation data across samples. Quality and quantity of DNA and the degree of bisulfite conversion can vary from sample to sample.
To assess acceptability of the precision of sodium bisulfite conversion and subsequent MethyLight assays, we measured PMR values (degrees of methylation) for the CDKN2A, MLH1, and MGMT promoters in 272 cases of colon cancer. We demonstrated that overall distributions of PMR values for all three loci were bimodal (either PMR
1 or >10) with only rare cases showing PMR values between 3 and 5. For most cases, CVs of PMR were sufficiently low to definitively differentiate methylation-positive tumor from methylation-negative tumor, with a PMR cutoff of 4. Thus, our data indicate that precision of sodium bisulfite conversion and quantitative Methy-Light assays is acceptable, and we can reliably determine the status of methylation for most cases by a single MethyLight run.
In addition, for each of the three genes tested, promoter methylation was strongly associated with loss of respective protein expression. Low PMR values between 0 and 4 were principally unrelated to loss of expression and, therefore, seemed to represent biologically insignificant low levels of methylation. There were infrequent cases with loss of expression without promoter methylation and cases with intact protein expression despite PMR values >10. The former may reflect different mechanisms of gene silencing, including gene deletion or mutation, and the latter may reflect partial methylation or mono-allelic methylation. Difficulty in interpreting immunohistochemistry, particularly for MGMT and CDKN2A (p16), may also contribute to some of the discrepant results.
In conclusion, sodium bisulfite conversion is reproducible, and subsequent quantitative real-time PCR methylation assays have acceptable precision. Quantitative promoter methylation data are highly correlated with loss of protein expression. Carefully validated quantitative MethyLight assays will be useful in both research and clinical molecular diagnostics.
| Note Added in Proof |
|---|
|
|
|---|
|
| Acknowledgments |
|---|
| Footnotes |
|---|
Supported by the National Institutes of Health (grants P01 CA87969-03 and P01 CA55075-13).
P.W.L. is shareholder and Scientific Advisory Board Member of Epigenomics, AG, which has a commercial interest in the development of DNA methylation markers for disease detection and diagnosis. None of the work described in this manuscript was supported by Epigenomics, AG.
Accepted for publication December 6, 2005.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. C. Christensen, J. J. Godleski, C. J. Marsit, E. A. Houseman, C. Y. Lopez-Fagundo, J. L. Longacker, R. Bueno, D. J. Sugarbaker, H. H. Nelson, and K. T. Kelsey Asbestos exposure predicts cell cycle control gene promoter methylation in pleural mesothelioma Carcinogenesis, August 1, 2008; 29(8): 1555 - 1559. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. A.T. Rodriguez, A. S.L. Cheng, P. S. Yan, D. Potter, F. J. Agosto-Perez, C. L. Shapiro, and T. H.-M. Huang Epigenetic repression of the estrogen-regulated Homeobox B13 gene in breast cancer Carcinogenesis, July 1, 2008; 29(7): 1459 - 1465. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. Iafrate and D. N. Louis "MGMT for pt mgmt": Is Methylguanine-DNA Methyltransferase Testing Ready for Patient Management? J. Mol. Diagn., July 1, 2008; 10(4): 308 - 310. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Vlassenbroeck, S. Califice, A.-C. Diserens, E. Migliavacca, J. Straub, I. Di Stefano, F. Moreau, M.-F. Hamou, I. Renard, M. Delorenzi, et al. Validation of Real-Time Methylation-Specific PCR to Determine O6-Methylguanine-DNA Methyltransferase Gene Promoter Methylation in Glioma J. Mol. Diagn., July 1, 2008; 10(4): 332 - 337. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. S. Kristensen, T. Mikeska, M. Krypuy, and A. Dobrovic Sensitive Melting Analysis after Real Time- Methylation Specific PCR (SMART-MSP): high-throughput and probe-free quantitative DNA methylation detection Nucleic Acids Res., April 1, 2008; 36(7): e42 - e42. [Abstract] [Full Text] [PDF] |
||||
![]() |
Q. Feng, S. E. Hawes, J. E. Stern, L. Wiens, H. Lu, Z. M. Dong, C. D. Jordan, N. B. Kiviat, and H. Vesselle DNA Methylation in Tumor and Matched Normal Tissues from Non-Small Cell Lung Cancer Patients Cancer Epidemiol. Biomarkers Prev., March 1, 2008; 17(3): 645 - 654. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ogino and A. Goel Molecular Classification and Correlates in Colorectal Cancer J. Mol. Diagn., January 1, 2008; 10(1): 13 - 27. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Bonanno, E. Shehi, D. Adlerstein, and G. M. Makrigiorgos MS-FLAG, a Novel Real-Time Signal Generation Method for Methylation-Specific PCR Clin. Chem., December 1, 2007; 53(12): 2119 - 2127. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ogino, T. Kawasaki, G. J Kirkner, Y. Suemoto, J. A Meyerhardt, and C. S Fuchs Molecular correlates with MGMT promoter methylation and silencing support CpG island methylator phenotype-low (CIMP-low) in colorectal cancer Gut, November 1, 2007; 56(11): 1564 - 1571. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. W. Coolen, A. L. Statham, M. Gardiner-Garden, and S. J. Clark Genomic profiling of CpG methylation and allelic specificity using quantitative high-throughput mass spectrometry: critical evaluation and improvements Nucleic Acids Res., September 25, 2007; 35(18): e119 - e119. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Dejeux, V. Audard, C. Cavard, I. G. Gut, B. Terris, and J. Tost Rapid Identification of Promoter Hypermethylation in Hepatocellular Carcinoma by Pyrosequencing of Etiologically Homogeneous Sample Pools J. Mol. Diagn., September 1, 2007; 9(4): 510 - 520. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ogino, A. Hazra, G. J. Tranah, G. J. Kirkner, T. Kawasaki, K. Nosho, M. Ohnishi, Y. Suemoto, J. A. Meyerhardt, D. J. Hunter, et al. MGMT germline polymorphism is associated with somatic MGMT promoter methylation and gene silencing in colorectal cancer Carcinogenesis, September 1, 2007; 28(9): 1985 - 1990. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ogino, T. Kawasaki, G. J. Kirkner, P. Kraft, M. Loda, and C. S. Fuchs Evaluation of Markers for CpG Island Methylator Phenotype (CIMP) in Colorectal Cancer by a Large Population-Based Sample J. Mol. Diagn., July 1, 2007; 9(3): 305 - 314. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Mikeska, C. Bock, O. El-Maarri, A. Hubner, D. Ehrentraut, J. Schramm, J. Felsberg, P. Kahl, R. Buttner, T. Pietsch, et al. Optimization of Quantitative MGMT Promoter Methylation Analysis Using Pyrosequencing and Combined Bisulfite Restriction Analysis J. Mol. Diagn., July 1, 2007; 9(3): 368 - 381. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Nosho, H. Yamamoto, T. Takahashi, M. Mikami, H. Taniguchi, N. Miyamoto, Y. Adachi, Y. Arimura, F. Itoh, K. Imai, et al. Genetic and epigenetic profiling in early colorectal tumors and prediction of invasive potential in pT1 (early invasive) colorectal cancers Carcinogenesis, June 1, 2007; 28(6): 1364 - 1370. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ogino, T. Kawasaki, G. J. Kirkner, M. Loda, and C. S. Fuchs CpG Island Methylator Phenotype-Low (CIMP-Low) in Colorectal Cancer: Possible Associations with Male Sex and KRAS Mutations J. Mol. Diagn., November 1, 2006; 8(5): 582 - 588. [Abstract] [Full Text] [PDF] |
||||
![]() |
S Ogino, M Cantor, T Kawasaki, M Brahmandam, G J Kirkner, D J Weisenberger, M Campan, P W Laird, M Loda, and C S Fuchs CpG island methylator phenotype (CIMP) of colorectal cancer is best characterised by quantitative DNA methylation analysis and prospective cohort studies Gut, July 1, 2006; 55(7): 1000 - 1006. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. B. Coleman and A. G. Rivenbark Quantitative DNA Methylation Analysis: The Promise of High-Throughput Epigenomic Diagnostic Testing in Human Neoplastic Disease J. Mol. Diagn., May 1, 2006; 8(2): 152 - 156. [Full Text] [PDF] |
||||
|
|