High-resolution melting (HRM) shows great promise for high-throughput, rapid genotyping of individual polymorphic loci. We have developed HRM assays for genotyping single nucleotide polymorphisms (SNP) in several key genes that are involved in methyl metabolism and may directly or indirectly affect the methylation status of the DNA. The SNPs are in the 5,10-methylenetetrahydrofolate reductase (MTHFR; C677T and A1298C), methionine synthetase (MTR; 5-methyltetrahydrofolate-homocysteine methyltransferase; A2756G), and DNA methyltransferase 3b (DNMT3b; C46359T and C31721T) loci. The choice of short amplicons led to greater melting temperature (Tm) differences between the two homozygous genotypes, which allowed accurate genotyping without the use of probes or spiking with control DNA. In the case of MTHFR, there is a second rarer SNP (rs4846051) close to the A1298C SNP that may result in inaccurate genotyping. We masked this second SNP by placing the primer over it and choosing a base at the polymorphic position that was equally mismatched to both alleles. The HRM assays were done on HRM capable real-time PCR machines rather than stand-alone HRM machines. Monitoring the amplification allows ready identification of samples that may give rise to aberrant melting curves because of PCR abnormalities. We show that samples amplifying markedly late can give rise to shifted melting curves without alteration of shapes and potentially lead to misclassification of genotypes. In conclusion, rapid and high-throughput SNP analysis can be done with probe-free HRM if sufficient attention is paid to amplicon design and quality control to omit aberrantly amplifying samples. (Cancer Epidemiol Biomarkers Prev 2008;17(5):1240–7)

Abnormal methylation patterns are one of the hallmarks of cancer. Methylation of CpG islands in the promoter region of many genes, including tumor suppressor genes such as the cell cycle inhibitor p16INK4a, the DNA repair genes BRCA1, MLH1, and MGMT, and the p53 regulator p14ARF, has been shown to shut down their expression (1, 2). It is still incompletely understood what underlies this alteration of methylation patterns and which susceptibility factors are involved. Much effort has been put into solving these questions, but they still remain largely unanswered.

Common variants in genes involved in the metabolism of the methyl group are likely candidates for the variation underlying propensity to methylation in normal tissues as well as in tumors (3, 4). A sufficient supply of the methyl group donor S-adenosyl methionine (SAM) is important to maintain a normal methylation pattern (5, 6). Because the polymorphisms of the 5,10-methylenetetrahydrofolate reductase (MTHFR), 5-methyltetrahydrofolate-homocysteine methyltransferase (MTR), and DNA methyltransferase 3b (DNMT3b) genes studied here either influence or are influenced by the levels of SAM, they are of particular interest, especially because they have been reported to modify the risk of getting different types of malignancies. However, further investigations are needed. Not all forms of cancers have been investigated in this regard, and some results need validation. Interestingly, some of the variants have been shown to be associated with an increased risk of getting some cancers and a decreased risk of getting others (7-9).

SAM is synthesized using dietary methionine or methionine generated from homocysteine. MTR methylates homocysteine to generate methionine and thus influences the cellular levels of SAM. MTHFR catalyses the reduction of 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate, which is the carbon donor for the de novo synthesis of methionine (Fig. 1). The common germ-line variants of MTHFR studied here (677T and 1298C) are less active (10, 11), and this can lead to higher levels of homocysteine and a deficiency of methyl group donors. The same is likely to be true for the MTR 2756G allele, but so far no one has been able to express human MTR in active form at sufficient levels to evaluate the biochemical effects of this polymorphism (12, 13). DNMT3b uses SAM to transfer methyl groups to DNA and is both responsible for de novo and maintenance DNA methylation. Overactivity of this gene has been linked to methylation of tumor suppressor genes and cancer (14, 15). Folate is an important player in these pathways as well, as it is required for the synthesis of methionine and SAM (Fig. 1). It has been shown that a deficiency of this nutrient results in DNA hypomethylation. The role of folate and the above germ-line variants in relation to cancer has been reviewed (16, 17).

Figure 1.

Folate and methionine metabolism. Enzymes are indicated in bold. THF, tetrahydrofolate; TS, thymidylate synthase; SAH, S-adenosyl homocysteine.

Figure 1.

Folate and methionine metabolism. Enzymes are indicated in bold. THF, tetrahydrofolate; TS, thymidylate synthase; SAH, S-adenosyl homocysteine.

Close modal

MTHFR also plays a role in the synthesis of dTMP from dUMP as MTHFR and the enzyme catalyzing this reaction (thymidylate synthase) both utilize 5,10-methylenetetrahydrofolate (Fig. 1). Reduced MTHFR activity leads to elevated levels of 5,10-methylenetetrahydrofolate, which in turn leads to lower levels of dUMP and higher levels of dTMP. This can help prevent cancer because deficiency in dTMP increases the rate of misincorporation of dUMP into DNA, which has been shown to increase the rate of DNA strand breaks and other chromosomal damage (18). Thus, reduced MTHFR activity has been shown to correlate with both increased and decreased risk of getting different forms of malignancies (8, 9).

The DNMT3b promoter C46359T single nucleotide polymorphism (SNP) is located 149 bp upstream from the transcription start site. The C-to-T change significantly increases the promoter activity of the DNMT3b gene (19). This up-regulates the expression of the gene and may in turn lead to aberrant methylation of CpG islands in some tumor suppressor genes. Nevertheless, both the C and the T alleles have been associated with increased risk of getting different malignancies (7). For this reason, further investigations are needed to clarify the function of this SNP in relation to cancer.

The DNMT3b C31721T SNP is intronic, but when the T allele is present the sequence motif MATWAAT (where M is A or C and where W is A or T) is missing (20). This motif is recognized by the transcription factor N-Oct-3 encoded by the POU3F2 gene. This gene has been shown to regulate the tumorigenic potential of human melanoma cells (21). Strong evidence for an association between this SNP and breast cancer risk has been reported (20).

High-resolution melting (HRM) analysis is a recently developed methodology made possible by the use of a new generation of fluorescent dyes, which do not inhibit PCR when they intercalate into double-stranded DNA at saturating levels (22). As the temperature increases, and the DNA “melts,” the dye is released and stops fluorescing. This change in fluorescence is sequence specific and can be monitored by appropriately designed instrumentation. Heterozygotes are particularly easy to identify because the heteroduplexes formed before the melting step have a characteristic melting profile. Homozygotes are more difficult to differentiate because their melting curves are usually very similar with often marginal differences in their melting temperature (Tm).

We used HRM to develop reliable and low-cost SNP genotyping assays that did not require additional probes or spiking with control DNA. This involved developing short amplicons that had Tm differences above 0.4°C to 0.5°C for the homozygotes. These assays were tested on the Rotor-gene 6000 (Corbett Research) and on the LC480 (Roche). Although we did not test the stand-alone HR-1 and LightScanner HRM machines (Idaho), the assays we developed should translate well to those machines (23).

Samples and DNA Extraction

The investigations were done after approval by the Peter MacCallum Cancer Centre Ethics of Human Research Committee (Project 02/70). We obtained 94 peripheral blood samples from the Melbourne Branch of the Red Cross Blood Bank. DNA was extracted from the buffy-coat fraction of blood using the X-tractor Gene (Corbett Research) according to the manufacturer's protocol. We did the incubation step after loading the Liquid Sample Digest Buffer for 20 min at 37°C instead of 10 min at room temperature. The DNA concentrations of the samples were measured using a NanoDrop ND-1000 UV-Vis Spectrophotometer (NanoDrop Technologies).

Primer Design, PCR, and HRM Conditions

One set of primers was designed for each of the MTHFR and DNMT3b (C46359T) assays and two for each of the MTR and DNMT3b (C31721T) assays. Primers were designed to have an optimal annealing temperature of 60°C using Primer Express 1.5 (Applied Biosystems). All amplicons were chosen so that they contained a single melting domain using the Poland program (http://www.biophys.uni-duesseldorf.de/local/POLAND/poland.html). The primer sequences and amplicon information can be found in Table 1. PCR cycling and HRM analysis were done on the Rotor-gene 6000 (Corbett Research) or on the LC480 (Roche).

Table 1.

Primers and amplicon information (University of California-Santa Cruz Genome Browser, September 2007)

PolymorphismDirectionSequenceAmplicon size (bp)Region spanned
MTHFR Forward 5′-GCACTTGAAGGAGAAGGTGTCTG 50 11778944-11778994 of  
C677T Reverse 5′-AGCTGCGTGATGATGAAATCG  chromosome 1 
MTHFR Forward 5′-GGGGAGGAGCTGACCAGTGA 54 11777032-11777085 of  
A1298C Reverse 5′-GAGGTAAAGAACCAAGACTTCAAAGACAC  chromosome 1 
MTR A2756G Forward 5′-TGTTCCCAGCTGTTAGATGAAAATC 129 235115044-235115172 of  
 Reverse 5′-AACAAGCAAAAATCTGTTTCTACCACTT  chromosome 1 
MTR A2756G Forward 5′-ATCATGGAAGAATATGAAGATATTAGACAGG 57 235115092-235115148 of 
 Reverse 5′-ACTTACCTTGAGAGACTCATAATGG  chromosome 1 
DNMT3b Forward 5′-TTCTGGCCCCGCCAGAC 58 30837901-30837958 of  
C46359T Reverse 5′-GAGAGGACACTCACTGGGCCTT  chromosome 20 
DNMT3b Forward 5′-ATGAGATGGCATTATTGTAAAAATGAGAC 70 30863274-30863343 of  
C31721T Reverse 5′-CATCCAGTCTTCATTGTCATTAGCAC  chromosome 20 
DNMT3b Forward 5′-CCTGGCAGGCAATAAACTTCACA 100 30863248-30863347 of  
C31721T Reverse 5′-CGTGCATCCAGTCTTCATTGTCA  chromosome 20 
PolymorphismDirectionSequenceAmplicon size (bp)Region spanned
MTHFR Forward 5′-GCACTTGAAGGAGAAGGTGTCTG 50 11778944-11778994 of  
C677T Reverse 5′-AGCTGCGTGATGATGAAATCG  chromosome 1 
MTHFR Forward 5′-GGGGAGGAGCTGACCAGTGA 54 11777032-11777085 of  
A1298C Reverse 5′-GAGGTAAAGAACCAAGACTTCAAAGACAC  chromosome 1 
MTR A2756G Forward 5′-TGTTCCCAGCTGTTAGATGAAAATC 129 235115044-235115172 of  
 Reverse 5′-AACAAGCAAAAATCTGTTTCTACCACTT  chromosome 1 
MTR A2756G Forward 5′-ATCATGGAAGAATATGAAGATATTAGACAGG 57 235115092-235115148 of 
 Reverse 5′-ACTTACCTTGAGAGACTCATAATGG  chromosome 1 
DNMT3b Forward 5′-TTCTGGCCCCGCCAGAC 58 30837901-30837958 of  
C46359T Reverse 5′-GAGAGGACACTCACTGGGCCTT  chromosome 20 
DNMT3b Forward 5′-ATGAGATGGCATTATTGTAAAAATGAGAC 70 30863274-30863343 of  
C31721T Reverse 5′-CATCCAGTCTTCATTGTCATTAGCAC  chromosome 20 
DNMT3b Forward 5′-CCTGGCAGGCAATAAACTTCACA 100 30863248-30863347 of  
C31721T Reverse 5′-CGTGCATCCAGTCTTCATTGTCA  chromosome 20 

NOTE: The position of the interfering SNP is shown in bold.

For the Rotor-gene 6000, the intercalating dye used was SYTO 9 (Invitrogen), and the reaction mixtures were made up using HotStarTaq (Qiagen) and consisted of 2 to 20 ng genomic DNA, 1× PCR buffer, 2.5 mmol/L MgCl2 total, 200 nmol/L of each primer [400 nmol/L were used in the DNMT3b (C46359T) assay], 200 μmol/L deoxynucleotide triphosphates, 5 μmol/L SYTO 9, 0.5 units HotStarTaq polymerase, and PCR-grade water in a total volume of 20 μL. For the LC480, the LC480 HRM Scanning Master (Roche) was used with the same concentrations of primers, MgCl2, and genomic DNA in the same total volume as the HotStarTaq/SYTO 9-containing MasterMix. The cycling protocol when using HotStarTaq started with one cycle of 95°C for 15 min, whereas the cycling protocol when using the LC480 HRM Scanning Master started with one cycle of 95°C for 10 min as less time is needed to activate the enzyme used in the latter MasterMix.

The cycling protocol for both machines was 45 cycles of 95°C for 10 s, 60°C for 10 s, and 72°C for 20 s, 1 cycle of 97°C for 1 min, and a melt from 65°C to 90°C for all assays, except for the DNMT3b C46359T assay. In this assay, 50 cycles were used, and the annealing step was done as a touchdown PCR with 6 cycles decreasing 1°C/cycle from 70°C to 64°C. For the melt on the Rotor-gene 6000, the temperature was increased at the rate of 0.2°C/s, and for the LC480, the temperature was increased at the rate of 1°C/s with the acquisitions set at 30 per °C. All reactions were done in duplicate.

The MTHFR C677T and the DNMT3b C31721T assays also worked well when using a more rapid cycling protocol of 45 cycles of 95°C for 5 s, 60°C for 5 s, and 72°C for 10 s. This was not the case for the MTR A2756G assay as the melting curves of all the replicates showed considerable variation with these conditions (see Results and Fig. 2C and D).

Figure 2.

MTR A2756G. A. Rotor-gene 6000 HRM normalized graph. B. LC480 normalized graph. A and B,red, heterozygous (AG); blue, homozygous for the A allele (AA); green, homozygous for the G allele (GG). Y axis, normalized fluorescence; X axis, temperature (°C). C. Amplification (raw channel). Uneven amplification was observed when using shorter PCR steps. D. HRM normalized graph of the samples shown in C. Black, homozygous for the A allele. Right-shifted melting curves for these samples were observed.

Figure 2.

MTR A2756G. A. Rotor-gene 6000 HRM normalized graph. B. LC480 normalized graph. A and B,red, heterozygous (AG); blue, homozygous for the A allele (AA); green, homozygous for the G allele (GG). Y axis, normalized fluorescence; X axis, temperature (°C). C. Amplification (raw channel). Uneven amplification was observed when using shorter PCR steps. D. HRM normalized graph of the samples shown in C. Black, homozygous for the A allele. Right-shifted melting curves for these samples were observed.

Close modal

HRM Analysis

The analysis was done with the software supplied with each machine. For the LC480, there is the possibility of temperature shifting the melting curves. This software feature makes it easier to distinguish different samples giving rise to differently shaped melting curves. However, because the two different homozygotes have identically shaped melting curves, temperature shifting was not done in this study.

The samples are from individuals of a variety of ethnicities. Because of this, the actual genotyping results may not be comparable with results derived from a single population. However, the data are shown as it has been proven useful as a quality-control measure to ensure that expected proportions of all the genotypes were observed. In all cases, heterozygotes were readily distinguished by their characteristic biphasic melting patterns. However, in some cases, the two homozygote genotypes were less readily distinguished from each other.

MTHFR C677T (rs1801133)

The homozygotes had Tm values that differed by ∼0.6°C independent of the MasterMix used. These were easily distinguished by HRM analysis on both machines when the data were normalized for the fluorescence before and after the melting transitions (Fig. 3A and B). We typed 94 samples and found 35 CC, 48 CT, and 11 TT.

Figure 3.

MTHFR C677T. A. Rotor-gene 6000 HRM normalized graph. B. LC480 normalized graph. LC480 normalized graph. Genotypes were readily distinguised and are indicated by arrows for both.

Figure 3.

MTHFR C677T. A. Rotor-gene 6000 HRM normalized graph. B. LC480 normalized graph. LC480 normalized graph. Genotypes were readily distinguised and are indicated by arrows for both.

Close modal

MTHFR A1298C (rs1801131)

It was difficult to design primers in this region because of the proximity of other SNPs. In particular, the position of the rs4846051 SNP is close to the MTHFR A1298C SNP. To avoid having two SNPs in between the primers, we designed the reverse primer to span the interfering SNP (see Table 1). The SNP is a C/T SNP, [listed in the SNP database (dbSNP BUILD 126) as an A/G SNP]. The reverse primer was designed to match neither the C nor the T expected in the antisense strand of this SNP by having a C at that position.

The size of the amplicon in this assay resulted in homozygous samples that differed by ∼0.8°C in Tm when using the SYTO 9 MasterMix and by ∼1.2°C when using the LC480 HRM Scanning MasterMix. These were readily distinguished by HRM analysis on both machines when the data were normalized for the fluorescence before and after the melting transitions (Fig. 4A and B). We typed 94 samples and found 9 CC, 46 AC, and 39 AA.

Figure 4.

MTHFR A1298C. A. Rotor-gene 6000 HRM normalized graph. B. LC480 normalized graph. Genotypes were readily distinguished and are indicated by arrows for both. C. Amplification (raw channel). Two samples stayed in the exponential phase, denoted as outliers. D. HRM normalized graph of the samples shown in C. Right-shifted melting curves of the two samples staying in the exponential phase were observed. These are denoted as outliers.

Figure 4.

MTHFR A1298C. A. Rotor-gene 6000 HRM normalized graph. B. LC480 normalized graph. Genotypes were readily distinguished and are indicated by arrows for both. C. Amplification (raw channel). Two samples stayed in the exponential phase, denoted as outliers. D. HRM normalized graph of the samples shown in C. Right-shifted melting curves of the two samples staying in the exponential phase were observed. These are denoted as outliers.

Close modal

As with other studies, we found partial linkage disequilibrium with the 667C>T polymorphism; in the 19 samples that were homozygous for the C allele, 18 were found to be homozygous for the C allele at the MTHFR C677T SNP as well, and all samples that were homozygous for the T allele at the MTHFR C677T SNP were homozygous for the A allele at the MTHFR A1298C SNP.

In one run on the Rotor-gene 6000, one replicate of two of the samples remained in the exponential phase, whereas all the other sample replicates went into the plateau phase (Fig. 4C). This resulted in a melting curve that was shifted to the right for both (Fig. 4D). One of these two samples is actually homozygous for the A allele and the other for the C allele. This could have resulted in an AA sample being interpreted as a CC sample if samples had not been run in duplicate and if normal PCR had been used instead of real-time PCR before the melting analysis.

DNMT3b C46359T (rs2424913)

The homozygotes had Tm values that differed by ∼0.7°C. These were readily distinguished by HRM analysis on both machines when the data were normalized for the fluorescence before and after the melting transitions (see Fig. 5A and B). We typed 94 samples and found 33 CC, 38 CT, and 23 TT.

Figure 5.

DNMT3b C46359T and C31721T. A.DNMT3b C46359T SNP: Rotor gene 6000 HRM normalized graph. B.DNMT3b C46359T SNP: LC480 normalized graph. C.DNMT3b C31721T SNP: Rotor gene 6000 HRM normalized graph. D.C31721T SNP: LC480 HRM normalized graph. Genotypes were readily distinguised and are indicated by arrows for all.

Figure 5.

DNMT3b C46359T and C31721T. A.DNMT3b C46359T SNP: Rotor gene 6000 HRM normalized graph. B.DNMT3b C46359T SNP: LC480 normalized graph. C.DNMT3b C31721T SNP: Rotor gene 6000 HRM normalized graph. D.C31721T SNP: LC480 HRM normalized graph. Genotypes were readily distinguised and are indicated by arrows for all.

Close modal

DNMT3b C31721T (rs406193)

For the large amplicon, the homozygotes had Tm values that differed by ∼0.4°C (Fig. 5C). The size of the amplicon in this assay resulted in poor separation of the homozygotes on the LC480 (data not shown). The smaller amplicon resulted in homozygotes that differed by ∼0.7°C in Tm values (Fig. 5D). All genotypes could then easily be determined on both machines. We typed 94 samples and found 72 CC, 21 CT, and 1 TT.

MTR A2756G (rs1805087)

For the large amplicon, the homozygotes had Tm values that differed by ∼0.3°C (Fig. 2A), once again leading to poor separation of the homozygotes on the LC480 (data not shown). The shorter amplicon resulted in homozygotes that differed by ∼0.8°C in Tm values (Fig. 2B). All genotypes could easily be determined on both machines. We typed 94 samples and found 6 GG, 24 AG, and 64 AA.

In one run on the Rotor-gene 6000 using shorter denaturation, annealing, and elongation steps (see Materials and Methods), we found samples that amplified much later than the others (Fig. 2C). This resulted in melting curves that were shifted to the right (Fig. 2D). These samples are actually homozygous for the A allele not the C allele. Again, this could have resulted in AA samples being interpreted as CC samples if samples had not been run in duplicate and if normal PCR was used instead of real-time PCR. The melting curves of all the replicates varied in Fig. 2D. This is caused by uneven amplification of the samples in general under these conditions (Fig. 2C).

Many different methods for detecting SNPs have been developed. All have their inherent strengths and weaknesses (24). All of these techniques, with the exception of fluorescent probe-based SNP genotyping, require removal of the PCR product for further analysis, making them more laborious and prone to human error as well as causing potential problems with PCR contamination.

Closed-tube assays, in which the genomic DNA and all the reagents required for amplification and genotyping are added at the same time, have become increasingly more attractive due to the continuous development of this approach. Originally, labeled primers or TaqMan-type probes were required. However, due to the introduction of saturating fluorescent DNA dyes, closed-tube SNP genotyping by amplicon melting has become simple, rapid, and inexpensive. These dyes, such as LCGreen, SYTO 9, and the dye in the LC480 HRM Scanning Master (this dye has no commercial name yet), bind to double-stranded DNA and fluoresce. Because the dyes do not inhibit PCR, it is possible to use them at saturating concentrations and thus to accurately monitor the melting of double-stranded DNA. Again, this has eliminated the need for any fluorescently labeled oligonucleotides and the produced amplicons can be directly analyzed after PCR. Generally, the method is called closed-tube SNP genotyping by amplicon melting (25).

The assays developed here are based on this method. We used the Rotor-gene 6000 and the LC480 because they are HRM-enabled real-time machines. The real-time capability allows us to record the amplification profiles and then to readily identify which aberrant melting patterns may have been due to irregular amplification. Real-time data are convenient for assay optimization as well.

Additional probes (26, 27) or control DNA spiking (28) have been used previously to separate homozygotes, but genotyping of nearly all SNPs can be done without the use of probes, and only G/C and A/T SNPs may require spiking with control DNA if the nearest neighboring bases are complementary (29). We obtained identical results for all 94 samples when using two different dyes on the two different HRM instruments. Furthermore, two different primer pairs for the MTR and DNMT3b SNPs gave identical results for all samples as well. Thus, SNP genotyping by HRM is a very reliable method that does not require confirmation by sequencing or other methods.

Homozygotes for all amplicons were readily distinguished by the Rotor-gene 6000 although their Tm values differed by only 0.3°C and 0.4°C in the MTR A2756G and DNMT3b C31721T assays respectively. The smaller size of the amplicon in the MTHFR C677T, MTHFR A1298C, and DNMT3b C46359T assays resulted in the two different homozygotes having Tm values that differed by 0.6°C, 0.8°C/1.2°C (dependent on the MasterMix; see Results), and 0.7°C respectively.

We found a bigger spread of the melting curves for each genotype when using the LC480. This is not surprising because a plate-based system has intrinsically more temperature variation than a rotor-based system. This did not matter as long as the amplicon size was kept small. However, we could not reliably type some larger amplicons, e.g., the MTR A2756G 129 bp and DNMT3b C31721T 100 bp (data not shown), but had no problems typing the same SNPs using the primers giving rise to shorter amplicons.

We recommend designing SNP assays so that the Tm difference is greater than 0.4°C for the LC480. In general, this means that the amplicon size at least should be less than 100 bp, but the Tm difference is sequence dependent as well. However, amplicon size seems to be the most important factor in determining Tm differences. The MTHFR A1298C assay showed that Tm differences can be reagent and/or equipment dependent as well. It has been shown previously that ionic strength significantly affects Tm (30). The advantage of optimizing the reaction for the LC480 is that higher throughput is possible compared to the rotor-based system of the Rotor-gene 6000.

Efficient SNP genotyping requires the presence of only one SNP per amplicon. The primers should not overlie a SNP as this may affect the differential amplification. The MTHFR C1298A showed that SNPs in primers cannot always be avoided but that the problem can be worked around by introducing a mismatch at the SNP position in the primer. In spite of the mismatch in our primer, we were able to get good amplification. This might not have been the case if the position of the mismatch were closer to the 3′ end, because this end needs to be stable for the polymerase to extend.

It has been considered that real-time amplification data is unnecessary for SNP genotyping by HRM analysis (27, 29). However, when we occasionally found outliers, these were almost always the result of individual reactions remaining in the exponential phase or amplifying very late. It has been noted previously by us that late amplifying samples should be treated with caution (31). For this reason, observation of real-time amplification data should be included in HRM result interpretations, and we recommend not scoring samples that amplified particularly late.

The results for the MTHFR A1298C and MTR A2756G SNPs showed that individual reactions staying in the exponential phase or amplifying markedly later than the majority of reactions can show right-shifted melting curves without alteration of shapes. This can lead to misinterpretation of results if certain precautions are not taken. Standardizing for DNA input helps prevent these problems, but even if samples are standardized for DNA input, some replicates occasionally amplify later than the majority. We were able to genotype samples that differed by ∼10-fold in input DNA and found no direct association between outliers and DNA concentration within this window. However, when trying to genotype samples that differed in the range from 100-fold to 1,000-fold, we had major problems and found many more outliers (data not shown). Running samples in duplicate or triplicate is preferable for this reason as it can identify outliers. We have found that choosing longer denaturation, annealing, and elongation steps can in some cases also facilitate even amplification of samples. This was the case for the MTR A2756G assay, in which the uneven amplification (Fig. 2C) was due to shorter denaturation, annealing, and elongation steps. Thus, the rare event that samples amplify markedly late in spite of standardized DNA concentrations can sometimes be explained by less optimized PCR conditions. Apart from that, we suggest that random effects in the PCR as well as pipetting errors may also be part of the explanation. The phenomena that individual replicates stay in the exponential phase and thus amplify to much higher levels are rarely observed and hard to explain. After optimization of assays, and running samples giving rise to outliers again on both machines, we were able to score all samples for all the SNPs. Thus, a very high percentage of success can be expected, and generally minimal optimization is needed when using the HRM approach.

Compared with TaqMan-based SNP genotyping, the HRM approach is much more cost-effective. Less optimization of assays is needed when no probes are used, and higher percentage of success can be expected. Pyrosequencing is another convenient method for SNP genotyping. However, this technique requires expensive and dedicated equipment and the percentage of success is highly dependent on DNA quality and quantity. The use of SNP chips has also become more widely used, and would be preferred relative to HRM, when conducting studies involving very large numbers of SNPs. However, in studies involving small numbers of SNPs, HRM will be the method of choice.

Germ-line variants in the MTHFR, MTR, and DNMT3b genes have undergone many association studies in relation to different malignancies (3, 4, 7-9, 32-37). Some of those are contradictory and some have been difficult to replicate. We predict that more correlations between different malignancies and germ-line variants in these genes are to be discovered, and for this reason, reliable, low-cost, and high-throughput methodologies for detecting these variants should be welcomed in many laboratories.

In conclusion, we showed that HRM is one such method by developing genotyping assays specific for the MTHFR C677T, MTHFR A1298C, DNMT3b C46359T, DNMT3b C31721T, and MTR A2756G SNPs. This was done without the use of additional probes or spiking with control DNA. These assays are reliable, of low cost, and are capable of high throughput analysis. Being able to review the real-time PCR amplification is a valuable tool for checking whether aberrant HRM profiles may be due to abnormal amplification.

No potential conflicts of interest were disclosed.

Grant support: U.S. Department of Defense Breast Cancer Research Program and Cancer Council of Victoria (A. Dobrovic).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

We thank Angela Tan, Ida Candiloro, Thomas Mikeska, and Michael Krypuy for critical comments on the manuscript drafts.

1
Esteller M. Aberrant DNA methylation as a cancer-inducing mechanism.
Annu Rev Pharmacol Toxicol
2005
;
45
:
629
–56.
2
Herman JG, Baylin SB. Gene silencing in cancer in association with promoter hypermethylation.
N Engl J Med
2003
;
349
:
2042
–54.
3
Friso S, Girelli D, Trabetti E, et al. The MTHFR 1298A>C polymorphism and genomic DNA methylation in human lymphocytes.
Cancer Epidemiol Biomarkers Prev
2005
;
14
:
938
–43.
4
Paz MF, Avila S, Fraga MF, et al. Germ-line variants in methyl-group metabolism genes and susceptibility to DNA methylation in normal tissues and human primary tumors.
Cancer Res
2002
;
62
:
4519
–24.
5
Jacob RA, Gretz DM, Taylor PC, et al. Moderate folate depletion increases plasma homocysteine and decreases lymphocyte DNA methylation in postmenopausal women.
J Nutr
1998
;
128
:
1204
–12.
6
Pogribny IP, Basnakian AG, Miller BJ, Lopatina NG, Poirier LA, James SJ. Breaks in genomic DNA and within the p53 gene are associated with hypomethylation in livers of folate/methyl-deficient rats.
Cancer Res
1995
;
55
:
1894
–901.
7
Montgomery KG, Liu MC, Eccles DM, Campbell IG. The DNMT3B C->T promoter polymorphism and risk of breast cancer in a British population: a case-control study.
Breast Cancer Res
2004
;
6
:
R390
–4.
8
Piyathilake CJ, Macaluso M, Johanning GL, Whiteside M, Heimburger DC, Giuliano A. Methylenetetrahydrofolate reductase (MTHFR) polymorphism increases the risk of cervical intraepithelial neoplasia.
Anticancer Res
2000
;
20
:
1751
–7.
9
Slattery ML, Potter JD, Samowitz W, Schaffer D, Leppert M. Methylenetetrahydrofolate reductase, diet, and risk of colon cancer.
Cancer Epidemiol Biomarkers Prev
1999
;
8
:
513
–8.
10
Frosst P, Blom HJ, Milos R, et al. A candidate genetic risk factor for vascular disease: a common mutation in methylenetetrahydrofolate reductase.
Nat Genet
1995
;
10
:
111
–3.
11
Weisberg I, Tran P, Christensen B, Sibani S, Rozen R. A second genetic polymorphism in methylenetetrahydrofolate reductase (MTHFR) associated with decreased enzyme activity.
Mol Genet Metab
1998
;
64
:
169
–72.
12
Harmon DL, Shields DC, Woodside JV, et al. Methionine synthase D919G polymorphism is a significant but modest determinant of circulating homocysteine concentrations.
Genet Epidemiol
1999
;
17
:
298
–309.
13
Skibola CF, Forrest MS, Coppede F, et al. Polymorphisms and haplotypes in folate-metabolizing genes and risk of non-Hodgkin lymphoma.
Blood
2004
;
104
:
2155
–62.
14
Kimura F, Seifert HH, Florl AR, et al. Decrease of DNA methyltransferase 1 expression relative to cell proliferation in transitional cell carcinoma.
Int J Cancer
2003
;
104
:
568
–78.
15
Robertson KD, Uzvolgyi E, Liang G, et al. The human DNA methyltransferases (DNMTs) 1, 3a and 3b: coordinate mRNA expression in normal tissues and overexpression in tumors.
Nucleic Acids Res
1999
;
27
:
2291
–8.
16
Fang JY, Xiao SD. Folic acid, polymorphism of methyl-group metabolism genes, and DNA methylation in relation to GI carcinogenesis.
J Gastroenterol
2003
;
38
:
821
–9.
17
Friso S, Choi SW. Gene-nutrient interactions in one-carbon metabolism.
Curr Drug Metab
2005
;
6
:
37
–46.
18
Blount BC, Mack MM, Wehr CM, et al. Folate deficiency causes uracil misincorporation into human DNA and chromosome breakage: implications for cancer and neuronal damage.
Proc Natl Acad Sci U S A
1997
;
94
:
3290
–5.
19
Wang GLZ, Mao L, Spitz MR, Wei Q. Functional relevance of C46359T in the promoter region of human DNMT3B6 [Abstract].
Proc Am Assoc Cancer Res
2004
;
45
:
2913
.
20
Cebrian A, Pharoah PD, Ahmed S, et al. Genetic variants in epigenetic genes and breast cancer risk.
Carcinogenesis
2006
;
27
:
1661
–9.
21
Thomson JA, Murphy K, Baker E, et al. The brn-2 gene regulates the melanocytic phenotype and tumorigenic potential of human melanoma cells.
Oncogene
1995
;
11
:
691
–700.
22
Gudnason H, Dufva M, Bang DD, Wolff A. Comparison of multiple DNA dyes for real-time PCR: effects of dye concentration and sequence composition on DNA amplification and melting temperature.
Nucleic Acids Res
2007
;
35
:
e127
.
23
Herrmann MG, Durtschi JD, Wittwer CT, Voelkerding KV. Expanded instrument comparison of amplicon DNA melting analysis for mutation scanning and genotyping.
Clin Chem
2007
;
53
:
1544
–8.
24
Mamotte CD. Genotyping of single nucleotide substitutions.
Clin Biochem Rev
2006
;
27
:
63
–75.
25
Wittwer CT, Reed GH, Gundry CN, Vandersteen JG, Pryor RJ. High-resolution genotyping by amplicon melting analysis using LCGreen.
Clin Chem
2003
;
49
:
853
–60.
26
Liew M, Seipp M, Durtschi J, et al. Closed-tube SNP genotyping without labeled probes/a comparison between unlabeled probe and amplicon melting.
Am J Clin Pathol
2007
;
127
:
341
–8.
27
Vandersteen JG, Bayrak-Toydemir P, Palais RA, Wittwer CT. Identifying common genetic variants by high-resolution melting.
Clin Chem
2007
;
53
:
1191
–8.
28
Palais RA, Liew MA, Wittwer CT. Quantitative heteroduplex analysis for single nucleotide polymorphism genotyping.
Anal Biochem
2005
;
346
:
167
–75.
29
Liew M, Pryor R, Palais R, et al. Genotyping of single-nucleotide polymorphisms by high-resolution melting of small amplicons.
Clin Chem
2004
;
50
:
1156
–64.
30
Owczarzy R, You Y, Moreira BG, et al. Effects of sodium ions on DNA duplex oligomers: improved predictions of melting temperatures.
Biochemistry
2004
;
43
:
3537
–54.
31
Krypuy M, Newnham GM, Thomas DM, Conron M, Dobrovic A. High resolution melting analysis for the rapid and sensitive detection of mutations in clinical samples: KRAS codon 12 and 13 mutations in non-small cell lung cancer.
BMC Cancer
2006
;
6
:
295
.
32
Campbell IG, Baxter SW, Eccles DM, Choong DY. Methylenetetrahydrofolate reductase polymorphism and susceptibility to breast cancer.
Breast Cancer Res
2002
;
4
:
R14
.
33
Li SY, Rong M, Iacopetta B. Germ-line variants in methyl-group metabolism genes and susceptibility to DNA methylation in human breast cancer.
Oncol Rep
2006
;
15
:
221
–5.
34
Matsuo K, Suzuki R, Hamajima N, et al. Association between polymorphisms of folate- and methionine-metabolizing enzymes and susceptibility to malignant lymphoma.
Blood
2001
;
97
:
3205
–9.
35
Shannon B, Gnanasampanthan S, Beilby J, Iacopetta B. A polymorphism in the methylenetetrahydrofolate reductase gene predisposes to colorectal cancers with microsatellite instability.
Gut
2002
;
50
:
520
–4.
36
Sharp L, Little J, Schofield AC, et al. Folate and breast cancer: the role of polymorphisms in methylenetetrahydrofolate reductase (MTHFR).
Cancer Lett
2002
;
181
:
65
–71.
37
Curtin K, Bigler J, Slattery ML, Caan B, Potter JD, Ulrich CM. MTHFR C677T and A1298C polymorphisms: diet, estrogen, and risk of colon cancer.
Cancer Epidemiol Biomarkers Prev
2004
;
13
:
285
–92.