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


From the Department of Pathology,*
University of Utah
School of Medicine, and the Associated and Regional University
Pathologists Institute for Clinical and Experimental
Pathology,
Salt Lake City, Utah
| Abstract |
|---|
|
|
|---|
20 bp). We describe a solution-based
fluorescence method that discriminates between wild-type and mutant
sequences using a dsDNA binding dye, and interrogates a region
of >200 nucleotides. This method is based on melting theory and
entails fluorescence monitoring of the melting temperatures of
GC-clamped amplicons subjected to gradual and progressive thermal
denaturation in the presence of a constant concentration of urea.
Heterozygous samples are easily identified by the lower melting
temperatures of the less thermodynamically stable heteroduplex
mismatches from the wild-type:mutant DNA hybrids as compared to the
more stable wild-type Watson-Crick duplexes. All of the four possible
sets of mismatches (A·G/T·C, T·G/A·C,
G·G/C·C, and T·T/A·A) represented in 17 heterozygous
mutations distributed throughout the length of 20 different amplicons
(104 to 212 bp), were distinguished from the wild-type by their
altered melting profiles. This methodology is advantageous in that it
obviates gel electrophoresis or labeled oligonucleotide probes.
Significantly, it expands the region of interrogation for
detection of single-base changes using fluorescence-based methods in
solution, and is amenable for automation and adaptation to
high-throughput systems.
| Introduction |
|---|
|
|
|---|
Recent trends favor fluorescence-based technologies for the detection of specific nucleic acid sequences.12 These assays have typically exploited any one of several fluorescence chemistries.13,14 The nonspecific methods incorporate a double-stranded DNA dsDNA binding dye such as SYBR Green I into the amplification reaction, and the fidelity of product detection is dependent on the inherent specificity of the PCR amplification conditions.14 The sequence-specific probe-based methods incorporate oligonucleotides that hybridize to a sequence within the amplicon, thereby providing an additional parameter for verification of product identity. The specificity of this hybridization interaction has been further exploited for the identification of single-base changes by virtue of the fact that the single-base mismatches within the hybridization probe to DNA target hybrids exhibit lower Tms than perfectly complementary strands.14 However, the detection of single-base changes using probe-based methods is limited to very short segments of DNA (<20 bp) and would require several probes in multiple separate reactions for screening larger regions.
It is desirable to detect single-base changes in long DNA fragments in solution because such methods dramatically reduce assay times and are easier to automate. In this study, we describe a strategy, designated "mutational scanning by temperature ramping and chemical denaturation" (MUSTERD) that permits the detection of single-base alterations in solution using SYBR Green I and fluorescence analysis of PCR product Tm. After amplification, the sample to be investigated is rapidly cooled to effect annealing of DNA strands and subjected to a gradual and uniform temperature gradient augmented by a constant concentration of urea. These denaturing conditions convert the dsDNA into single strands as the Tms of the respective domains are reached. Attachment of a high-melting domain (GC clamp) to the DNA fragment of interest converts the fragment into a lower melting domain and renders its entire length accessible for detection of point mutations. Using synthetic mutant, genomic WT, and heterozygous cell line DNAs each containing 1 of 17 single-base substitutions in a total of 20 different amplicons, we show that this method detects single-base alterations in longer PCR products than has been feasible with solution-based fluorescence methods.
| Materials and Methods |
|---|
|
|
|---|
WT Sequences
DNA was extracted from peripheral blood lymphocytes of healthy individuals and human placenta using standard methods. All samples were confirmed by direct sequencing using the ABI PRISM 377 (Perkin Elmer Applied Biosystems, Emeryville, CA) to contain the WT N-RAS and Factor V genes.
Mutant Sequences
We synthesized nine 104-bp DNA fragments each containing one of
the possible base substitutions in codon 61 of the N-RAS
gene (GenBank Accession no. L00041; nucleotides 38 to 141) by PCR
site-directed mutagenesis as previously described.15,16
Three additional PCR products with single-base substitutions at
nucleotide 18 (CTG to ATG transversion at codon
56), nucleotide 61 (CAA to CTA transversion at
codon 70), and nucleotide 89 (CTC to CTA
transversion at codon 79). Segments (160 and 212 bp) of the Factor V
gene (GenBank Accession no. Z99572, nucleotides 62,832 to 62,991 and
62,783 to 62,994, respectively) were also synthesized containing base
substitutions at positions 1689, 1691, and 169217
(Tables 1 and 2)
.
The heterozygous state for each mutation was simulated by mixing
equivalent quantities of synthesized homozygous mutant DNA and WT
genomic DNA. Four cell lines harboring mutations in codon 61 of the
N-RAS gene were also investigated (HL-60, SK-N-SH, SW1271,
and HT1080 that harbor heterozygous CAA
CTA,
CAA
AAA, CAA
CGA, and
CAA
AAA mutations,
respectively).16
Patient samples previously determined to
contain heterozygous mutations of the Factor V Leiden locus were also
evaluated using MUSTERD. All samples were confirmed by direct
sequencing using the ABI PRISM 377 (Perkin Elmer Applied Biosystems) to
contain the WT and the desired mutant sequences of the N-RAS
and Factor V genes.
|
|
Computer algorithms that predict the melting behavior of a DNA fragment based on its nucleotide sequence have been previously described.18 We assessed the 104-bp N-RAS and 160-bp and 212-bp Factor V gene sequences with and without the attached 40-bp GC clamp using the WINMELT DNA analysis program (MedProbe AS, Oslo, Norway), which predicts the thermodynamic stability of DNA duplexes based on the DNA sequence.
PCR Amplification of Target Sequences
Rapid-cycle PCR amplification was performed in a microvolume
fluorimeter (LightCycler; Roche Molecular Biochemicals, Indianapolis,
IN)19
using primers directed at the N-RAS and
Factor V genes as summarized in Tables 1 and 2
. Briefly, a 100-ng
aliquot of template DNA was amplified in a 20-µl reaction containing
1x PCR buffer (50 mmol/L Tris, pH 8.5, 3.0 mmol/L
MgCl2, 500 µg/ml bovine serum albumin), four
deoxynucleotide triphosphates at 200 µmol/L each, 0.2 µmol/L of
each primer, 10% dimethyl sulfoxide (by volume), and 0.8 U of
Promega DNA polymerase (Promega Corporation, Madison, WI) with
176 ng of TaqStart antibody (ClonTech, Palo Alto, CA). PCR entailed a
denaturation step at 94°C (0 seconds), annealing at 50°C (0
seconds), and extension at 72°C (10 seconds) for 45 cycles. Each
reaction included the dsDNA-binding dye SYBR Green I (1:30,000
dilution; Molecular Probes, Eugene, OR). In the LightCycler, filtered
excitation light (450 to 490 nm) from a blue light-emitting diode is
reflected from a 505-nm dichroic filter and focused on the capillary
tip. The samples are interrogated by paraxial epi-illumination. A
portion of the excitation light is passed up the capillary tube by
total internal reflection at the glass/air interface. Similarly,
emitted light is passed down the capillary tube and exits out at the
tip. The light is then reflected off a 505-nm dichroic filter, passed
through a 520- to 560-nm interference filter, and focused onto silicon
photodiodes for detection. Fluorescence signals were obtained once in
each cycle by sequential fluorescence monitoring of each tube at the
end of extension.
Standard Fluorescence-Melting Curve Analysis
Standard fluorescence-melting curve analysis entailing postamplification denaturation of non-GC-clamped amplified DNA segments in the absence of urea was performed as previously described.20,21
MUSTERD
The optimal conditions for MUSTERD were empirically determined by
varying the length of the GC clamp, SYBR Green I concentration,
temperature ramping rate, and urea concentration, both individually and
in combination. Once established, no further manipulation of the
conditions was necessary as different amplicons were evaluated. The
optimal conditions were applied to all amplicons as detailed below.
After amplification, the amplicon was removed from the capillary tubes
by reverse centrifugation. Ten µl of the PCR products were added to
10 µl of 26.6 mol/L urea and mixed vigorously to produce a final urea
concentration of 13.3 mol/L. Ten µl of this amplicon/urea mixture
were subjected to MUSTERD analysis. The procedure entailed sample
heating to 95°C (10 seconds), cooling to 33°C (2 minutes), followed
by continuous acquisition of fluorescence during gradual temperature
elevation from 33 to 85°C at 0.02°C/second in the presence of a
constant concentration of urea (13.3 mol/L). The ramp rate of
0.02°C/second was achieved using custom software kindly provided by
Dr. Carl T. Wittwer (University of Utah School of Medicine).
SYBR Green I fluorescence signals were continuously monitored during
the postamplification denaturation process. For easier visualization of
the Tm, the fluorescence versus temperature curves
(F versus T) were converted
into fluorescence-melting peaks by plotting the negative derivative of
fluorescence over temperature versus temperature
(-dF/dT versus T). Up to six
replicates per sample were assessed by MUSTERD analysis, and the
melting profiles were analyzed in a blinded manner and scored for the
presence of mutations (Table 3)
.
|
The within-run standard deviations and P values (paired two-tailed Students t-test) for the experimentally derived Tm values from replicate samples assessed by MUSTERD were determined using the Microsoft Excel 2000 statistical tool (Microsoft, Redmond, WA).
Dilutional Analyses
Dilutional assays were performed using serial dilutions of HT1080 cell line DNA into WT genomic DNA.
| Results |
|---|
|
|
|---|
The 104-bp N-RAS amplicon showed two distinct melting
domains (Figure 1A)
. Similar analysis of
the 160-bp Factor V DNA fragment showed three distinct melting domains
with a different distribution from that observed for the
N-RAS amplicon (data not shown). The 212-bp Factor V
fragment yielded a profile with three melting domains showing Tms
differing by only 7°C (data not shown). Importantly, the GC-clamped
counterparts for all amplicons showed only two melting domains
comprising a higher melting domain corresponding to the 40-bp GC clamp
and a lower melting domain corresponding to the entire length of the
DNA fragment of interest. Figure 1A
illustrates the predicted melting
profiles of the 104-bp N-RAS product, with and without the
40-bp clamp.
|
Standard postamplification-melting curve analysis was performed in
the absence of urea using non-GC-clamped amplicons created from WT DNA
sequences, heterozygous cell line DNA, and artificial heterozygotes
(50/50 mixtures of homozygous mutant and WT DNA mixtures). The three
targets (104-bp N-RAS, 160-bp and 212-bp Factor V) showed an
abrupt decline in fluorescence corresponding to their experimental Tms.
Figure 1B
, inset, shows the WT 104-bp N-RAS amplicon with an
abrupt decline in fluorescence at 80°C on the F
versus T curves. A corresponding melting peak was
evident at 80°C on the -dF/dT versus
T plots (Figure 1B)
. In our experiments, standard
fluorescence-melting curve analysis failed to distinguish between WT
and heterozygous mutant samples in products of
104 bp in length.
Figure 1B
illustrates identical fluorescence-melting profiles
(same Tm) for a 104-bp segment of N-RAS exon 2 for both the
WT sample and the HT1080 cell line, which harbors a heterozygous
CAA
AAA point mutation at codon 61 of the
N-RAS gene.
MUSTERD
The incorporation of a GC clamp to a primer generates amplicons with a high-melting domain (GC clamp) attached to a lower melting domain (fragment of interest).22 To create GC-clamped amplicons for MUSTERD, 40-bp clamps were attached to one primer to generate PCR amplicons with a high-melting domain juxtaposed to the DNA segment of interest. The amplified fragment maintains a double-stranded configuration until it reaches its Tm.11,22 Because the disparity in Tms is most evident in the lower melting domain of the fragment, single-base changes are distinguished in the lower melting domain whereas the higher melting domains of both WT and mutant samples show identical Tms.22
Wild-Type Sequences
MUSTERD analysis performed on GC-clamped amplicons derived from
the 104-bp N-RAS and 160- and 212-bp Factor V amplicons
showed two distinct melting domains for all products (Figures 2 to 4)
. Acquisition of fluorescence was
performed during temperature ramping from 33 to 85°C at
0.02°C/second. F versus T plots
revealed similar profiles consisting of a lower melting domain
corresponding to the amplified DNA segment of interest, and a higher
melting domain corresponding to the 40-bp GC clamp. The Tms of these
two melting domains were better visualized by assessing the
fluorescence-melting peaks derived from the negative derivative of
fluorescence over temperature (-dF/dT) versus
temperature (T) plots (Figure 2A
, peaks I and II). WT
samples assessed in quintuplicate yielded consistent melting profiles
with identical Tms.
|
|
|
We assessed a total of 17 different heterozygous single-base
substitutions in 20 amplicons using MUSTERD. As detailed in Table 3
,
the mutations included all nine possible single-base substitutions in
codon 61 of a 104-bp N-RAS exon 2 amplicon and three
single-base changes distributed throughout the length of the 104-bp
N-RAS amplicon (Figure 3)
. In addition, we studied five
single-base substitutions in a 160-bp Factor V gene segment
encompassing the Leiden locus (Table 3)
. Three of these mutations were
also evaluated in the context of a 212-bp amplicon (Table 3)
.
Using the 104-bp N-RAS amplicon, all nine heterozygous
mutations in N-RAS codon 61 were distinguishable from the WT
sequence by virtue of their variant melting profiles. In contrast to
the WT sequence that displayed only one melting transition in the lower
melting domain, the heterozygous mutants exhibited either one or two
melting subtransitions with lower Tms than the WT in the lower melting
domain (Figure 2, B and C)
.
Analysis of artificially synthesized single-base substitutions
distributed throughout the length of the 104-bp N-RAS DNA
sequence (within the lower melting domain of the GC-clamped amplicon)
was also performed using artificial heterozygotes composed of a 50/50
mixtures of WT genomic DNA and homozygous mutants. These mutants
consisted of a C to A substitution at nucleotide 18 (codon 56), an A to
T substitution at nucleotide 61 (codon 70), and a C to A substitution
at nucleotide 89 (codon 79). All of these heterozygous mutations were
easily distinguished from the WT sequence by their variant melting
profiles (Figure 3)
. Although melting
curve analysis of non-GC-clamped amplicons in the absence of urea
failed to distinguish mutant from WT sequences (Figure 1B)
, melting
curve analysis of non-GC-clamped amplicons in a 13.3 mol/L urea
solution discriminated between WT and heterozygous mutant sequences,
provided the mutation was located in the lower melting domain (data not
shown).
To assess the functionality of MUSTERD in a different sequence context,
we evaluated five different single-base substitutions in a 160-bp
Factor V gene amplicon (Table 3)
. The heterozygous G1689A, G1691A,
G1691T, and A1692G substitutions were readily distinguished from the WT
sequence by the variant melting peaks with lower Tms distinct from that
of the WT sample (Figure 4)
. However, the
heterozygous G1689T transversion was essentially indistinguishable from
the WT sequence (Table 3)
. In like manner, the G1689A transition and
the G1691T transversion in the 212-bp Factor V product were
distinguished from the WT sequence (Figure 4)
.
Dilutional Analysis
Serial dilutions were performed using DNA extracted from the N-RAS mutation-bearing cell line (HT1080) into WT genomic DNA. The assay sensitivity was 25% using our MUSTERD method (data not shown).
| Discussion |
|---|
|
|
|---|
Fluorescence detection of specific PCR products may use a nonspecific DNA dye or sequence-specific fluorescently labeled probes.20,21 dsDNA binding dyes such as ethidium bromide and SYBR Green I are relatively inexpensive and easy to use. For our experiments, we have used SYBR Green I that binds to the minor groove of dsDNA and exhibits an excitation maximum (497 nm) similar to that of fluorescein. SYBR Green I fluorescence correlates with the quantity of dsDNA in a sample such that continuous monitoring of the dye fluorescence during PCR has been used for quantitative PCR analysis in real-time.27 Because dsDNA denaturation is accompanied by loss of binding of the dsDNA dye, the Tm of the DNA duplex is observed as a dramatic decrease in SYBR Green I fluorescence in the F versus T curves, or alternatively as fluorescence melting peaks in the derivative (-dF/dT versus T) plots. Because the Tm of a DNA fragment is determined by its length, GC content, and specific DNA sequence, important information can be deduced from the melting characteristics of long PCR products. In this regard, we have previously shown that products of identical size, but with dissimilar GC content and sequence can be distinguished by their Tms.20,21 Other investigators have demonstrated the feasibility of single nucleotide polymorphism genotyping by Tm analysis, but these assays have required previous knowledge of the precise site of the base substitution for effective design of the appropriate allele-specific PCR primer.28 More recently, Lipsky and colleagues29 have also described a solution-based approach for the detection of single nucleotide polymorphisms, but the majority of the base changes were located within the middle third of the fragment of interest. Thus suitability for detection of mutations throughout the length of >100-bp fragments was not systematically addressed.
Our current study demonstrates another approach that takes advantage of the melting characteristics of dsDNA for the detection of single-base substitutions. Amplified mutant and WT DNA from the genomic region of interest are subjected to a gradual and uniform linear temperature gradient augmented by a constant concentration of urea. These denaturing conditions convert the double-stranded molecule into single strands as the Tms of the respective domains are reached. Fluorescence monitoring of the melting profiles during application of the temperature gradient yielded melting transitions corresponding to the GC clamp and the DNA fragment of interest, confirming the existence of the two distinct DNA-melting domains in solution.
The cyclical reiteration of denaturation and annealing during PCR
generates heteroduplex molecules composed of WT and mutant DNA
sequences in heterozygous samples. As predicted by melting theory and
experimental data,11
our method identified single-base
substitutions by the presence of one or two subtransitions within the
lower melting domain with lower Tms than the WT. These results can be
explained by the fact that heterozygous samples contain both the WT and
mutant sequences that when annealed give rise to homoduplexes
containing WT fragments and heteroduplexes consisting of mismatched
mutant:WT hybrids. Because the Watson-Crick homoduplexes are more
thermodynamically stable than the heteroduplexes, they exhibit higher
Tms than the less stable heteroduplex mismatches. Thus the heterozygous
mutants were identifiable by the presence of an additional
subtransition with a lower Tm than was observed for the WT sample. The
different heterozygous single-base substitutions yielded Tm shifts of
varying magnitudes depending on the destabilizing ability of the
resulting heteroduplex mismatch (Table 3)
. Nevertheless, no
false-positives were recorded and 18 of the heterozygous mutations in
20 amplicons analyzed were detected (sensitivity,
90%) (Table 3)
.
Our inability to detect all mutations may result from the relatively
stable mismatches present in these samples and the sequence contexts in
which they appear. Overall GC content as well as nearest neighbor
sequences may affect the applicability of this method to some targets.
However, we also observed that there was a diminution in the magnitude
of Tm shifts associated with specific mutations as the length of the
interrogated DNA segment increased (Table 3)
, suggesting that this
method may be most applicable to amplicons of the order of 200 bp. In
addition, our dilutional studies suggest that this method is better
suited to the detection of polymorphisms or germline mutations, rather
than somatic/acquired mutations.
In summary, we have shown here that screening for heterozygous mutations is feasible in solution using SYBR Green I generated fluorescence-melting peak analysis during gradual thermal denaturation of GC-clamped amplicons. This methodology is advantageous in that it is simple and rapid requiring only 40 minutes of postamplification analysis for mutational screening of single nucleotide polymorphisms or point mutations. Furthermore, gel electrophoresis is not required, neither is there a need for radioactive, biotinylated, or sequence-specific fluorescently labeled probes. Although other solution-based fluorescent methodologies used for detection of point mutations are limited to short segments, our current method is particularly applicable when the precise site of the mutation or polymorphism is unknown in advance, or spans a length of more than a few nucleotides. We anticipate that this method will find utility in the screening for single nucleotide polymorphisms or point mutations that may subsequently be confirmed by DNA sequencing or other methods of genotyping.
| Acknowledgements |
|---|
| Footnotes |
|---|
Supported by a National Institutes of Health grant (no. CA 83984 to K. S. J. E.-J.) and by the ARUP Institute for Clinical and Experimental Pathology.
Accepted for publication May 25, 2001.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. P. Vaughn and K. S. J. Elenitoba-Johnson High-Resolution Melting Analysis for Detection of Internal Tandem Duplications J. Mol. Diagn., August 1, 2004; 6(3): 211 - 216. [Abstract] [Full Text] |
||||
![]() |
C. T. Wittwer, G. H. Reed, C. N. Gundry, J. G. Vandersteen, and R. J. Pryor High-Resolution Genotyping by Amplicon Melting Analysis Using LCGreen Clin. Chem., June 1, 2003; 49(6): 853 - 860. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. N. Gundry, J. G. Vandersteen, G. H. Reed, R. J. Pryor, J. Chen, and C. T. Wittwer Amplicon Melting Analysis with Labeled Primers: A Closed-Tube Method for Differentiating Homozygotes and Heterozygotes Clin. Chem., March 1, 2003; 49(3): 396 - 406. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |