A novel method to quantify local CpG methylation density by regional methylation elongation assay on microarray
© Zhang et al. 2008
Received: 08 May 2007
Accepted: 31 January 2008
Published: 31 January 2008
Skip to main content
© Zhang et al. 2008
Received: 08 May 2007
Accepted: 31 January 2008
Published: 31 January 2008
DNA methylation based techniques are important tools in both clinical diagnostics and therapeutics. But most of these methods only analyze a few CpG sites in a target region. Indeed, difference of site-specific methylation may also lead to a change of methylation density in many cases, and it has been found that the density of methylation is more important than methylation of single CpG site for gene silencing.
We have developed a novel approach for quantitative analysis of CpG methylation density on the basis of microarray-based hybridization and incorporation of Cy5-dCTP into the Cy3 labeled target DNA by using Taq DNA Polymerase on microarray. The quantification is achieved by measuring Cy5/Cy3 signal ratio which is proportional to methylation density. This methylation-sensitive technique, termed RMEAM (regional methylation elongation assay on microarray), provides several advantages over existing methods used for methylation analysis. It can determine an exact methylation density of the given region, and has potential of high throughput. We demonstrate a use of this method in determining the methylation density of the promoter region of the tumor-related gene MLH1, TERT and MGMT in colorectal carcinoma patients.
This technique allows for quantitative analysis of regional methylation density, which is the representative of all allelic methylation patterns in the sample. The results show that this technique has the characteristics of simplicity, rapidness, specificity and high-throughput.
In the human genome, GC-rich DNA sequences are found frequently within the promoter and first exon of ~50% of all genes . These sequences, also known as CpG islands, can be targets of DNA methylation. An epigenetic phenomenon is known to be associated with genomic imprinting and X-chromosome inactivation, and essential for normal mammalian development . Both global hypomethylation and regional hypermethylation have been described in human tumor cell lines and a wide spectrum of cancers . Global hypomethylation has been associated with instability of chromosomal or microsatellite, while regional hypermethylation of CpG islands within promoter region of tumor suppressor genes is associated with transcriptional inactivation and represents an important mechanism of gene silencing in the pathogenesis of neoplasia [4, 5]. There is emerging evidence that each tumor may harbor multiple genes susceptible to promoter hypermethylation. Methylation profiles of multiple genes for each cancer type might have important prognostic implications for clinical monitoring, risk assessment, and even therapeutic considerations [6–8].
The present techniques commonly used for the methylation analysis are based on the bisulfite modification of the genomic DNA. Since sodium bisulfite treatment exclusively converts unmethylated cytosine to uracil under appropriate conditions , subsequent analysis to differentiate unconverted cytosine from converted uracil enables us to know the primary methylation status. It is the basis for methylation-specific PCR (MSP) , bisulfite DNA sequencing , enzymatic regional methylation assay (ERMA) , pyrosequencing , and mass spectrometry . MSP is widely used to analyze promoter methylation, although only a limited number of CpG sites can be analyzed by this method. Bisulfite DNA sequencing provides precise methylation status over an amplified region, but it requires large-scale sequencing of multiple plasmid clones. The major advantage of the pyrosequencing method compared to MSP is that the data are actual sequences rather than fluorescence data from PCR-based amplification. But highly repetitive thymines and the limitation of ~75 bp extension length  can affect assay reproducibility and accuracy. Mass spectrometry permits the high-throughput identification of methylation sites and the semiquantitative measurement at single or multiple CpG positions, but this method needs high quality and a large amount of samples, and sometimes false positive results are inevitable. Even though ERMA can determine overall methylation level within a region containing a number of CpG sites, but it requires radioactive labeling of DNA samples and subsequent cumbersome purification steps of the radiolabeled products. Recently, several groups have shown that the methylation status can be achieved by 5mC-antibodies or MBD-proteins [16–18]. However, the affinity and specificity would be further improved for methylation density detection.
Microarrays provide the powerful tools for mapping the epigenome and detecting patterns of DNA methylation in the genome level [19–24]. Methods Coupling CpG-recognizing restriction enzymes with microarray technique have been reported, such as differential methylation hybridization (DMH) [19, 23], methylation amplification DNA chip (MAD) , and microarray-based integrated analysis of methylation by isoschizomers (MIAMI) . The above methods are only able to analyze methylation at the restriction sites of the enzymes, which are potentially biased by mutations or polymorphisms at the sites. Another microarray approach for methylation pattern analysis is methylation-specific oligonucleotide (MSO) microarray [19, 20]. MSO uses bisulfite-treated DNA as template for non-discriminatory PCR amplification, followed by hybridization of the PCR products to glass slides with oligonucleotides that could discriminate between methylated and unmethylated cytosines at specific CpG positions. One of the great potentials of MSO is that multiple genes can be analyzed on the same array. A potential limitation of this method is that closely spaced CpGs may not be amenable to analyze if the gene in question is heterogeneously methylated. Recently, genomic tiling microarray has been used to profile DNA methylation patterns . But this method can not obtain exact informations of methylation level within the given region of specific caner-related genes.
Here, we describe a novel technique for the quantification of CpG methylation density of a given DNA region, which combines microarray-based hybridization and enzymatic elongation on microarray. The quantification is achieved by measuring Cy5/Cy3 signal ratio which is proportional to methylation density. In contrast to the other protocols, this method determines the methylation density of the entire amplified region, not only of a few CpG sites or those CpG dinucleotides that are covered by PCR primers/probes. The results show that this technique allows for the quantitative analysis of regional methylation density with simplicity, rapidness, specificity and high-throughput.
To evaluate the sensitivity and quantitative accuracy of our assay as well as determining any possible PCR bias , we then carried out mixing experiments using methylation positive clones and methylation negative clones. The mixtures were prepared prior to PCR in order to assess the possibility of bias during the PCR amplification . The mixtures of Cy3-labeled fully-methylated allele and non-methylated allele with different proportion merged the Cy3-labeled positive control. And then a series of microarray hybridization were performed. After hybridization, elongation was carried out. Optimal elongation conditions on microarray were determined in order to allow elongation reactions to go to completion. Figure 3 shows that there is a linear correlation between the defined methylation density of the standard clone dilutions and the (Cy5sam*Cy3pos)/(Cy5pos*Cy3sam) signal ratios. By including mixtures of clones as standards (fully-methylated only, 75% fully-methylated/25% non-methylated, 50% fully-methylated/50% non-methylated, 25% fully-methylated/75% non-methylated and non-methylated only, corresponding to 100%, 75%, 50%, 25% and 0% methylation density) in every assay , it is possible to create a standard curve in order to control for differential specific activity and incorporation efficiency of the Cy5-dCTP. This standard curve can then be used to calculate the methylation density of unknown samples from the ratios (Cy5sam*Cy3pos)/(Cy5pos*Cy3sam). All standards and unknown samples were analyzed with eight replications.
We developed a novel approach for quantitative analysis of CpG methylation density on the basis of microarray-based hybridization and incorporation of Cy5-dCTP into the Cy3 labeled target DNAs by using Taq DNA Polymerase on microarray. This method provides several advantages over existing methods for quantitative methylation analysis. Firstly, it has potential for high-throughput analysis of DNA methylation profiles. CpG island hypermethylation has been reported to be linked to the silencing of many cancer-related genes. A DNA microarray can be designed and generated to contain a large amount of oligonucleotide probes capturing specifically target genes. Cancer-related genes can be parallelly amplified from investigated samples in a 96-well format or with multiple PCR, which can generate multiple target genes for hybridization and detection. Secondly, RMEAM can determine the methylation density of all CpG dinucleotides within the entire amplified region other than individual CpG sites. In the case of the MGMT promoter region, the detected methylation density includes the methylation status of 28 CpG dinucleotides. RMEAM produced a linear response when used for quantitative methylation analysis in mixing experiments by use of DNA from fully-methylated and non-methylated clones with a defined methylation status. The third feature of RMEAM is its simplicity and precision. Unlike ERMA method , RMEAM doesn't require radioactive labeling of DNA samples and subsequent cumbersome purification steps of the radiolabeled products in tubes. Hybridization, elongation and washing on microarray can be quickly and easily achieved within 4 hours. The use of an internal control (Cy3) allows for correction of potential substrate loss and for normalization of the results to the DNA amount that is finally analyzed by microarray scaner. The negative control in Figure 3 can remove the influence of cross-hybridization. For example, the net Cy3 intensities of row 1, 3, 4 and 5 of Figure 3A can be derived by subtraction of the row 2 negative control background. In row 2, the non-specific hybridization can't bring incorporation of Cy5-dCTP into the DNA due to the specificity of Taq DNA Polymerase. Dual-color fluorescence can be standardized by positive control, just like house-keeping gene in cDNA microarray. The accuracy of our new technique is shown by the high concordance with bisulfite sequencing data. However, in this study, only 7 clones or less per gene in each patient sample were sequenced, which might account for the slight discrepancies in methylation percentages (Figure 5). As cloning and final sequencing of a small number of alleles does not necessarily give representative data for the initial sample, the approach used in our current assay can provide more informations regarding the overall methylation density of the region from the sum of all alleles. Furthermore, RMEAM approach is less costly and labor intensive than the extensive bisulfite sequencing analysis needed for precise quantitation.
One possible limitation of our method is that it is not suitable for a detailed characterization of the pattern of DNA or the determination of methylation status of individual CpG sites. However, since difference of site-specific methylation may also lead to a change of methylation density in many tumor-related genes, It has been indicated that the methylation density of promoter CpG islands is more important than methylation occurred at a specific single CpG site for gene silencing [28–32]. Nevertheless, it has been very difficult to measure the level of methylation in a target area containing multiple CpG sites. Results with our new assay might be more biologically relevant, especially in patient samples with heterogeneous methylation patterns. Furthermore, it was suggested in an in vitro model that de novo CpG island methylation is not a single event, but rather a progressive process, and appears to be region specific . Therefore, it is reasonable to apply a method that accurately determines the methylation density of a CpG island region in order to investigate the dynamics of methylation changes in human malignancies.
One of the potential applications of RMEAM approach is to examine the changes in methylation patterns following treatment with demethylating drugs . Galm et al.  had demonstrated the application of ERMA technique for the quantitative analysis of methylaiton density, in which KG1a cells were treated in vitro with the increasing 5'aza-2'-deoxycytidine (DAC) dose within a given time course. The findings in that paper indicated that the changes of methylation density could help to optimize the dosage of demethylating drugs for clinical use by determining the extent of demethylation at a dose range with little cytotoxic side effects. Large-scale screenings of tumor-related genes under various conditions are needed to elucidate the functional role of DNA methylation and to screen for new antineoplastic drugs operative on that mechanism. Our new method could help elucidate the consequences of treating patients with such drugs and eventually help clarify their therapeutic value in future antineoplastic therapy regimens. Recent evidence has indicated that increased density of methylation within a susceptible CpG island is associated with more advanced stages of tumors. Our RMEAM method coupled with quantitative real-time PCR approaches to detect gene expression levels will contribute to more accurately assess the dynamics of methylation-mediated transcriptional silencing of tumor-related genes in carcinogenesis. Such assays would preferably be genome-wide to monitor both desired and undesired effects on methylation and gene expression.
We described a novel approach for quantitative analysis of CpG methylation density on the basis of microarray-based hybridization and incorporation of Cy5-dCTP into the Cy3 labeled target DNA by using Taq DNA Polymerase on microarray. Our results showed that this RMEAM technique allows for the quantitative analysis of regional methylation density with simplicity, rapidness, specificity and high-throughput.
Informed written consent was obtained from all patients and donors, and tissue collection was approved by each Institutional Review Board. The samples of primary tumor tissues were collected from colorectal carcinoma patients (n = 18) during surgery at Gulou Hospital and Zhongda Hospital (Nanjing, China). All tissue samples were fresh-frozen and stored at -80°C until further processing. Genomic DNA was extracted from tissues by digesting with proteinase K (0.5 mg/ml) in 10 mM TE (pH 8.0) buffer and 0.5% SDS at 48°C overnight, followed by a standard phenol/chloroform (1:1) extraction and precipitated with ethanol in a standard fashion.
Bisulfite treatment was carried out following the modified  procedure of Frommer et al. . Briefly, 5 μg of genomic DNA were digested with EcoRI (New England Biolabs) and denatured with 0.3 M NaOH for 15 min at 37°C. A freshly prepared solution of sodium bisulfite (2.5 M, pH 5.0) and hydroquinone (100 mM) was added to the denatured DNA, and the mixture was incubated at 55°C for 5 h. After desalting (Wizard Clean Up System; Promega), the DNA was desulphonated with 0.3 M NaOH for 15 min at 37°C. The solution was neutralized with 75 μl of 5 M ammonium acetate (pH 7.0), and the DNA was ethanol precipitated and resuspended in Tris-EDTA (pH 7.5).
Placental DNA treated in vitro with SssI methyltransferase (New England Biolabs) was used as fully-methylated alleles for APC, MLH1, TERT and MGMT gene. Meanwhile, DNA from normal lymphocytes was used as unmethylated alleles for MLH1, TERT and MGMT gene. Sodium bisulfite modification of the unmethylated DNA and M. SssI-treated DNA were then performed. The targeted CpG islands (see Fig. 2) were amplified from bisulfite-treated genomic DNA by triplex PCR using the following primer pairs [26, 36, 37]. The primers used were MLH1-Forward, 5'-TTTTTTAGGAGTGAAGGAGG-3'; MLH1-Reverse, 5'-ATAAAACCCTATACCTAATCTATC-3'; MGMT-Forword, 5'-GGATATGTTGGGATAGTT-3'; MGMT-Reverse, 5'-CCAAAAACCCCAAACCC-3';TERT-Forward, 5'-GGGTTATTTTATAGTTTAGGT-3'; TERT-Reverse, 5'-AATCCCCAATCCCTC-3'. These primers are specific for bisulfite-modified DNA without CpG dinucleotides. The PCRs were performed in 25-μl reactions containing PCR buffer with 1.8 mM MgCl2, 5 pmol primer, 1.25 units Hot-start DNA polymerase (TaKaRa), and 1 μl bisulfite-modified DNA (equivalent to 50 ng genomic DNA). After an initial preheating step of 5 min at 95°C, PCR was performed 40 cycles of 95°C for 30 s, 54°C for 30 s, and 72°C for 30 s, with a final extension at 72°C for 7 min. APC fully-methylated allele was used as positive control which was amplified by forward primer 5'-GGGGTTAGGGTTAGGTAGG-3' and reverse primer 5'-AACTACACCAATACAACCACATA-3' . All PCR products were gel purified and cloned into the pMD18-T vector according to the manufacturer's instructions (TaKaRa) and sequenced using the ABI sequencing system. Quantitation of the PCR product was performed with a UV spectrometer; 1 OD260 unit was calculated as 50 ng/μL.
Five probes (18–20 nucleotides in length, Fig. 2) were synthesized with an amino-linked C6 [NH2 (CH2)6] linker attached to its 5'end. Additional 10 nucleotides T were used to reduce the space influence. Each oligonucleotide was printed on the aldehyde-coated glass slides using a PixSys5500 microarrayer (Cartesian Technology Inc). After printed, the glass slides were incubated in a humid chamber at room temperature overnight, and then at 37°C for 2 h. The slides were washed thoroughly in 1% SDS solution to remove unbound oligonucleotides. After further treatment with a NaBH4 solution for 20 min, the slides were ready for hybridization. For target DNA labeling, PCR products of bisulfite-treated DNA were labeled at the 5'end of reverse primer with Cy3. The labeled products were resuspended in hybridization solution (1:3 dilution v/v). Then the mixture was denatured at 95°C for 5 min, immediately cooled on ice for 10 min and subsequently applied to the DNA microarray slides. Microarray hybridization was conducted in a moist hybridization chamber under a cover slip at 42°C for 3 h recommended by Schumacher et al .
After hybridization, the slide was rinsed and washed at room temperature with 2 × SSC-0.1% SDS for 10 min, 0.1 × SSC-0.1% SDS for 5 min, water for 5 min, and then dried by flowing nitrogen. Then 15 μL elongation system was applied to the hybridization region with a cover slip. The elongation system contained 1 × PCR buffer, 2.5 mM Mg2+, dNTPs (20 μM dGTP, 20 μM dATP, 20 μM dTTP, 0.2 μM dCTP, 0.2 μM Cy5-dCTP), and 1.5 IU Taq DNA Polymerase (TaKaRa). The incubation procedure was 42°C 10 min, 55°C 5 min, and 72°C 3 min. Subsequently, the slide was washed with 2 × SSC-0.5% SDS and dried, and it was ready to be scanned.
These microarray slides were scanned with ScanArray Lite microarray analysis systems (A Packard BioScience Company, USA). The fluorescence images were analyzed with GenePix Pro3.0 software. Each spot was defined by the positioning of a grid of circles over the array image. For each fluorescent image, the average pixel intensity within each circle was determined and a local background using mean pixel intensity was computed for each spot. Net signal was determined by subtraction of this local background from the mean average intensity for each spot and the nospecific influence of cross-hybridization from negative control. Signal intensities of individual spots were obtained and exported to Excel spreadsheets for further analysis. The intensity ratio of Cy5/Cy3 for each probe set was then obtained.
This work was funded by the National Natural Science Foundation of China (Project 60121101, 60701008), the Hi-Tech Research and Development Program of China (Project 2006AA020702), and the National Basic Research and Development Program (973 Program) of China (Grant No. 2006CB705602-9).
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.