Assessment for Melting Temperature Measurement of Nucleic Acid by HRM
High resolution melting (HRM), with a high sensitivity to distinguish the nucleic acid species with small variations, has been widely applied in the mutation scanning, methylation analysis, and genotyping. For the aim of extending HRM for the evaluation of thermal stability of nucleic acid secondary structures on sequence dependence, we investigated effects of the dye of EvaGreen, metal ions, and impurities (such as dNTPs) on melting temperature () measurement by HRM. The accuracy of HRM was assessed as compared with UV melting method, and little difference between the two methods was found when the DNA was higher than 40°C. Both insufficiency and excessiveness of EvaGreen were found to give rise to a little bit higher , showing that the proportion of dye should be considered for precise measurement of nucleic acids. Finally, HRM method was also successfully used to measure s of DNA triplex, hairpin, and RNA duplex. In conclusion, HRM can be applied in the evaluation of thermal stability of nucleic acid (DNA or RNA) or secondary structural elements (even when dNTPs are present).
High resolution melting (HRM) is a robust technique for mutation scanning, methylation, and genotyping, which possesses high sensitivity to distinguish DNA species from tiny differences in their melting temperatures based on sequence composition, length, GC content, or strand complementarity [1–4]. To obtain precise , data are required to be collected at narrow temperature increments of 0.02~0.2°C. Due to the high sensitivity to small variations, HRM can distinguish small differences in nucleic acid sequences, even for single nucleotide polymorphism (SNP) [5, 6]. Meanwhile, rapid measurement and analysis become easy due to the high-throughput work platforms (96 or 384 wells) coupled with HRM-compatible software. Therefore, HRM can be speculated as a potential approach for fast, accurate, and high-throughput evaluation of the thermal stability of nucleic acid, in the case that the conventional UV melting method does not work .
However, at present, HRM has been mainly applied to distinguish two DNA species longer than 60 bp in PCRs. It is unknown whether it can be applied to the high-throughput evaluation of thermal stability of short DNA or RNA sequences, which provide significant information for investigating sequence dependence or predication of secondary structures. From the viewpoint of basic data for thermodynamics, there has been no specific report on measuring of short DNA sequences by HRM, not to mention RNA . In this study, we checked whether HRM could measure s of short DNA duplex, triplex, or hairpin and even RNA duplex. We found that s of DNA duplexes with various GC contents and lengths (even shorter than 10 bp) could be well measured by HRM. Effects of the dye and metal ion as well as other composites (such as deoxyribonucleoside triphosphates and free single strands) were also investigated.
2. Materials and Methods
2.1. Materials and Reagents
2.2. Measurement Using HRM Method
DNA duplex, DNA hairpin, or RNA duplex was prepared in a solution containing 10 mM phosphate buffer (pH 7.4), EvaGreen (0.25x, 0.5x, 1x, 2x, or 5x), and NaCl (0-1 M). For experiments involving Mg2+, the buffer of 10 mM Tris-HCl (pH 7.1), 1x EvaGreen, and Mg2+ (1–100 mM) was used. For DNA triplex, DNA oligomers were prepared in a solution containing 1x EvaGreen, 2.5 mM KCl, 10 mM Na2HPO4, 2 mM KH2PO4, 10 mM MgCl2, and 0.5 mM spermine (pH 5.5).
For fluorescence melting profiles, the mixed oligomer solution (10 μL) was pipetted into 96-well microtiter plates; then, the microtiter plates were transferred to a PikoReal Real-Time PCR instrument (Thermo Scientific, Finland). Annealing was performed with a cooling rate of 0.1°C/s from 95°C to 10°C; then, fluorescence data were collected over a temperature range of 10–95°C in 0.1°C increments (the holding time was 2 seconds). All fluorescence data were acquired from at least three parallel tests in one plate and repeated twice.
2.3. Measurement Using UV Melting Method
UV melting assays were operated by a UV-Vis Spectrophotometer (UV-1800, SHIMADZU, Tokyo, Japan). For measurement of DNA duplexes, each DNA duplex (1 μM) was dissolved in a buffer as the cases of HRM, but with no EvaGreen. The solution was added into a cuvette with 1 mm or 10 mm path lengths. Melting profiles were obtained at 260 nm over the temperature ranges from 95°C to 10°C and then from 10°C to 95°C at a ramp rate of 1°C/min. All UV melting assays were performed for three times at least, and s were calculated by the first derivatives of melting curves.
3.1. Evaluation of Measurement by HRM
In this study, EvaGreen was adopted as the fluorescent dye in measurement by HRM, which was one of the most popular saturating dyes for HRM analysis. The detailed effects of DNA sequences, metal ions, and other factors on measurement by HRM were investigated.
3.1.1. The Limit of DNA Duplex Length for HRM
In order to investigate effects of DNA duplex length on measurement by HRM, we used a series of DNA duplexes with lengths from 6 bp to 20 bp (L6–L10, L15, and L20 in Table 1). The melting profiles of DNA duplexes are shown in Figure 1(a). The initial relative fluorescence intensity of 8–20 bp duplexes was at the range of 8000–11000, the intensity for the 7 bp duplex declined to about 5500, and the intensity of the 6 bp one was only about 1300, indicating that the 6 bp duplex was hard to be stained. The differential curves shown in Figure 1(b) also confirmed that only the melting peak of the 6 bp duplex was not detected and calculated. Therefore, EvaGreen can be used to measure of 7 bp or longer DNA duplex by HRM.
3.1.2. Effect of EvaGreen Concentration on Measurement by HRM
The fluorescent dye concentration is an important factor for HRM analysis [9, 10]. When the dye is unsaturated, the dissociated dye from melted duplexes may reincorporate into other DNA duplexes, resulting in less fluorescent signal change and obtaining a higher . In order to evaluate the effect of EvaGreen concentration on measurement by HRM, two sequences of GC3/15 (20% GC content) and GC12/15 (80% GC content) were used (Table 1). The obtained s for various EvaGreen concentrations (0.25x, 0.5x, 1x, 2x, or 5x) were shown in Figure 2. Interestingly, GC3/15 (1 μM) showed a lowest when the EvaGreen was 0.5x (Figure 2(a)). For higher concentrations, 2x or 5x, became higher. In the case of 5x, was even 2-3°C higher than that of 0.5x. At 0.25x concentration, became a little bit higher. This indicated that 0.5x EvaGreen was sufficient to bind with 1 μM 15 bp DNA duplex. Similar results were also obtained for duplex GC12/15 (Figure 2(b)).
Meanwhile, the DNA concentration effect on HRM was also investigated. We measured the DNA duplexes GC3/15 and GC12/15 at various strand concentrations from 0.1 to 10 μM in 1x EvaGreen by HRM. s by HRM were compared with that predicted by mfold, which is a popular software for calculating (see Supplementary Figure S1 in Supplementary Material available online at http://dx.doi.org/10.1155/2016/5318935). When the DNA concentration was less than 1 μM, the difference between two methods was less than 2.5°C. The difference became greater as DNA concentration increased to 5 or 10 μM. Obviously, for a higher concentration of DNA, more EvaGreen is required, and a reasonable proportion of dye to DNA is essential for measurement by HRM.
3.1.3. Evaluation of Accuracy by HRM
To apply HRM to measure of nucleic acid, its accuracy is a necessary aspect to be investigated by comparison with the traditional UV melting method. Two series of 15 bp or 20 bp DNA duplexes with various GC contents ( range of 30–90°C, Supplementary Tables S1 and S2) were measured by HRM and UV melting method, respectively. Some single stranded DNA may form secondary structures (e.g., the dG of F or R strand homodimer of GC14/15 were lower than −22 kcal/mol); however, we found that they did not affect the measurement of the duplex by HRM (Supplementary Figure S2). As shown in Figure 3(a) (15 bp DNA duplexes) and Figure 3(b) (20 bp DNA duplexes), there was little difference between the two methods when was higher than 40°C, although the difference became greater at the lower range. In detail, for DNA duplexes with s higher than 70°C, their s by HRM were almost consistent with those by UV melting method. For s of 40–65°C, the differences were 1-2°C. For s lower than 40°C, the difference was within 4°C. We also found that there was a good linear relationship between s by the two methods (Figure 3(c)) which was used as the standard curve to correct HRM data. Then, this linear formula was adopted to correct the HRM s of 20-bp DNA duplexes, and the results showed that the error of most corrected data was within 0.9°C as compared with UV data (Supplementary Table S2) and that the root mean-squared deviation (RMSD) of corrected was 0.5°C, less than that of mfold to UV (1.7°C).
3.1.4. Effect of Na+ or Mg2+ on Measurement by HRM
In terms of the metal ion effect, we measured DNA duplexes GC3/15 and GC12/15 at various concentrations of Na+ or Mg2+ by HRM (Figure 4). HRM s of both GC3/15 and GC12/15 increased greatly with Na+ concentration from 0 to 100 mM, and the increasement became slightly at the range from 100 to 1000 mM (Figure 4(a)). The same trend was also observed for s obtained by UV, reflecting the effect of Na+ on stabilization of DNA duplex. When the Na+ was lower than 400 mM, almost the same values were obtained by UV and HRM. When the Na+ was higher than 400 mM, s by UV became higher than those by HRM, and the difference increased with ion strength. According to their fluorescence melting curves (data not shown), the initial fluorescence intensity was lower when Na+ was at 0.6–1.0 M. The lower by HRM indicated that Na+ at a high concentration affected the binding of EvaGreen to DNA.
The effect of Mg2+ concentration was shown in Figure 4(b). Here, the Tris-HCl buffer was used instead of phosphate buffer because magnesium phosphate has a low solubility. As compared with NaCl (Figure 4(a)), the difference between HRM and UV became much greater, especially at higher concentrations of MgCl2. Even at 10 mM MgCl2, of GC3/15 by HRM was 2.7°C lower than that by UV. Similar results were also obtained for GC12/15, and the difference between two methods became even greater. These results indicated that Mg2+ competed for DNA binding with EvaGreen containing only one positive charge; then, EvaGreen dissociated at a lower temperature.
3.1.5. Effect of dNTPs on HRM
In some biological reactions, dNTPs were present in the solution. Because dNTPs have UV absorption at 260 nm, measurement of DNA duplexes by UV melting method becomes difficult. Whether dNTPs can affect measurement by HRM was determined. In the presence of 0.02–2.0 mM dNTPs, s of GC0/15~GC15/15 with various GC contents were measured, and almost no change was observed at various concentration of dNTPs (Figure 5). However, high concentration of dNTPs (≥0.2 mM) resulted in excessively high absorbance (Supplementary Figure S3), beyond the accurate scale of UV spectrophotometers.
3.2. HRM Applied in DNA Triplex, Hairpin, and RNA Duplex
Except for DNA duplex, measurement of DNA triplex and hairpin was also investigated by HRM. We constituted a DNA triplex (Tri-a/a′/b in Table 2) and identified it by UV melting method and HRM. As shown in Figure 6(a), two melting regions occurred in the fluorescence melting curve, which indicated that the s of duplex and triplex were 73.9 and 64.5°C. This was in line with the UV melting curves (Figure 6(b)) indicating that s of duplex and triplex were 73.7 and 63.9°C, respectively. Thus, HRM could detect the formation of DNA triplex.
In terms of DNA hairpin, we investigated a series of triloop hairpins possessing the same loop sequence, GAA, but differing in their length of stems (Supplementary Figure S4). We found that when the stem was 5 bp or longer, s could be obtained and the full width at half maximum of peak heights was as narrow as 10–12°C. It can be concluded that HRM can measure the of a DNA hairpin as short as 5 bp. In order to further investigate whether HRM can detect the effect of the loop sequence on the of a hairpin, we measured the GAA-loop hairpin (HP-S6, Table 2) and its mutants HP-M1~M4 with loops of “IAA” “GIA” “GAI” or “IAI” (where I, inosine, served to replace the certain base of GAA) by HRM. The obtained s of IAA-, GAI-, and IAI-loop hairpins (, , and °C, resp.) were much lower than those of GAA- and GIA-loop hairpins (, °C). This was in line with Moody and Bevilacqua’s findings by UV melting method indicating that the first and third base of loop dominated the stability of DNA triloop hairpin [11, 12]. Thereby, this method can exactly evaluate the stability of short DNA hairpin on sequence dependence.
Considering the significance of RNA SSEs, we measured the of a 21-bp RNA duplex (RR, Table 1) by HRM. It was found that RR presented good fluorescence melting curves at 1 μM and 10 μM, and the full width at half maximum of peak heights was only 5–5.5°C (Supplementary Figure S5). It was interesting that EvaGreen could bind to A-form conformation of RNA duplex. Therefore, HRM could be applied in analysis of short RNA duplex, or other RNA secondary structures.
In this study, whether HRM could measure s of small nucleic secondary structures was investigated. Surprisingly, HRM can measure as short as 7-bp DNA duplex with of about 20°C (Figure 1). For the first time, HRM was proved to be capable of measuring the of DNA triplex and short DNA hairpin (Figure 6). For hairpin structures, a minihairpin with 5 bp stem and 3 nt loop could also be analyzed by HRM. It can be concluded that HRM is suitable for evaluation of thermal stability of most secondary structural elements. Moreover, HRM measurement performs better stability and reproducibility compared with UV melting method. The standard deviation of the average by HRM was within 0.5°C, lower than that by UV melting method (within 1.0°C) (Supplementary Table S1). Because the work platform of 96- or 384-well plate could be used for HRM and the data could be analyzed by well-edited processing software, HRM has been used as a high-throughput approach for determining SNP, gene methylation, and mutant analysis on a large scale [13–15]. Obviously, based on our study, HRM shows the potential for high-throughput evaluation of nucleic acid secondary structural elements, which is significant to the thermodynamic studies of nucleic acid. s of thousands of DNA sequences can be obtained in hours, even using 96-well plates (Supplementary Figure S6).
Another advantage of measurement by HRM is that it is less affected by the presence of dNTPs (Figure 5). In routine PCRs or other approaches, the concentration of dNTPs is normally more than 0.2 mM [16–19]. However, UV melting method is hard to be used if the solution has compositions with strong UV absorption at 260 nm, such as enzymes and other additives in the buffer [20, 21]. Taken together, it can be concluded that HRM is suitable for measurement of practical values in most biological reactions.
In terms of measurement by HRM, some points should be noted including the proportion of dye to nucleic acid, the metal ion concentration, and other compositions. In some cases, the s obtained by HRM were higher than that of UV melting method, especially for low s (Figures 3(a) and 3(b)). The difference may be explained from two aspects. On the one hand, the binding of dye with positive charge can enhance the stability of DNA duplex, resulting in the fact that the dissociation of EvaGreen from the DNA duplex may occur at a relatively higher temperature. On the other hand, EvaGreen can bind to ssDNA, especially at a temperature below 40°C, and EvaGreen may remain binding with dissociated DNA, resulting in the fluorescence change delay and higher . Fortunately, there was a good linear relationship (the slope = 1.07, ) between HRM and UV (Figure 3(c)), which can be used as the standard curve to correct HRM s (Supplementary Table S2). Therefore, HRM can be utilized as an accurate approach to evaluate the thermal stability of nucleic acids.
Another important point to be considered is the proportion of dye to nucleic acid when highly accurate values are required [8, 9]. Our results showed that both the insufficiency and excessiveness of dye gave rise to a higher (Figure 2 and Supplementary Figure S1). In addition, the metal ion at high concentrations may affect the accuracy of HRM. When the Na+ concentration was higher than 400 mM or Mg2+ higher than 10 mM, HRM s remained almost constant (Figure 4), which did not conform to the fact that nucleic acid with high concentration of metal ions shows higher [22, 23]. This may be attributed to the competition between Na+ (or Mg2+) and EvaGreen for binding with DNA. Accordingly, HRM using EvaGreen is not suitable for the studies of salt dependence of nucleic acid stability [24, 25]. However, in most biological reactions, the Na+ concentration is lower than 400 mM, and Mg2+ concentration is not higher than 10 mM; thus, the accuracy of HRM is acceptable for most cases.
In conclusion, HRM can be used as a high-throughput approach for measurement of small nucleic acid secondary structures, including minihairpins, triplex, or RNA structures. The high throughput of this approach makes it possible for obtaining large amounts of data, for example, the sequence dependence of nucleic acid secondary structure. In fact, the work of thermodynamic evaluation of DNA and RNA hairpins by HRM is being carried out by our group and some valuable characteristics of sequence dependence have been obtained. It can be expected that HRM will play an important role in the thermodynamic study of nucleic acid.
The authors declare that there is no conflict of interests regarding the publication of this manuscript.
This work was supported by the National Natural Science Foundation of China [31571937 to Xingguo Liang] and by the National Natural Science Foundation of China (NSFC) Shandong Joint Fund for Marine Science Research Centre/NSFC-Shandong Joint Fund for Marine Science Research Centre [U1406402] to Xingguo Liang.
The supplementary material contains additional data used for assessing the melting temperature measurement of nucleic acid by HRM. Supplementary Tables S1 and S2 supply the detailed melting temperature (Tm) data of 15 bp and 20 bp DNA duplexes, respectively. The Tm difference between UV and HRM in Table S1 indicated that HRM measurement performed better stability and reproducibility compared with UV melting method. Moreover, the linear formula obtained from Table S1 can be used to correct the HRM Tms of 20-bp DNA duplexes (Table S2), and the error of most corrected data was within 0.9oC as compared with UV data. The HRM Tm of DNA duplex GC3/15 and GC12/15 at various strand concentrations from 0.1 to 10 μM in 1× EvaGreen was investigated (Figure S1).The Tm difference between HRM and prediction became greater as DNA concentration increased to 5 or 10 μM. The results demonstrated that a reasonable proportion of dye to DNA is essential for Tm measurement by HRM. The HRM melting curves of GC14/15 duplex, GC14/15 duplex with additional F or R strand, and only F strand or R strand were investigated and the results showed that the single strands did not affect the Tm measurement of the duplex by HRM (Figure S2). For the UV melting method, it was found that dNTPs (≥0.2 mM) resulted in a high absorbance beyond the accurate scale of UV spectrophotometers (Figure S3). Melting curves of a series of triloop hairpins with 3~8-bp stem and a 21-bp RNA duplex were detected by HRM (Figures S4 and S5), which demonstrated that HRM could be applied in Tm analysis of DNA and RNA secondary structural elements. For the melting curves of 96 samples on a 96-well plate by HRM (Figure S6), all the data could be obtained in 1 to 2 h. The results showed its high efficiency and high throughput.
A. Lehmensiek, M. W. Sutherland, and R. B. McNamara, “The use of high resolution melting (HRM) to map single nucleotide polymorphism markers linked to a covered smut resistance gene in barley,” Theoretical and Applied Genetics, vol. 117, no. 5, pp. 721–728, 2008.View at: Publisher Site | Google Scholar
M. Nakano, E. M. Moody, J. Liang, and P. C. Bevilacqua, “Selection for thermodynamically stable DNA tetraloops using temperature gradient gel electrophoresis reveals four motifs: d(cGNNAg), d(cGNABg), d(cCNNGg), and d(gCNNGc),” Biochemistry, vol. 41, no. 48, pp. 14281–14292, 2002.View at: Publisher Site | Google Scholar
H. Gudnason, M. Dufva, D. D. Bang, and A. Wolff, “Comparison of multiple DNA dyes for real-time PCR: effects of dye concentration and sequence composition on DNA amplification and melting temperature,” Nucleic Acids Research, vol. 35, no. 19, article e127, 2007.View at: Publisher Site | Google Scholar
M. Krypuy, G. M. Newnham, D. M. Thomas, M. Conron, and A. Dobrovic, “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, vol. 6, no. 1, article 295, 2006.View at: Publisher Site | Google Scholar