Home | About | Journals | Submit | Contact Us | Français |

**|**Nucleic Acids Res**|**v.31(16); 2003 August 15**|**PMC169985

Formats

Article sections

Authors

Related links

Nucleic Acids Res. 2003 August 15; 31(16): e93.

PMCID: PMC169985

Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, 1055 du P.E.P.S., PO Box 3800, Sainte-Foy, Quebec G1V 4C7, Canada

Received 2003 May 15; Accepted 2003 June 18.

Copyright © 2003 Oxford University Press

This article has been cited by other articles in PMC.

Fluorescent monitoring of DNA amplification is the basis of real-time PCR, from which target DNA concentration can be determined from the fractional cycle at which a threshold amount of amplicon DNA is produced. Absolute quantification can be achieved using a standard curve constructed by amplifying known amounts of target DNA. In this study, the mathematics of quantitative PCR are examined in detail, from which several fundamental aspects of the threshold method and the application of standard curves are illustrated. The construction of five replicate standard curves for two pairs of nested primers was used to examine the reproducibility and degree of quantitative variation using SYBER® Green I fluorescence. Based upon this analysis the application of a single, well- constructed standard curve could provide an estimated precision of ±6–21%, depending on the number of cycles required to reach threshold. A simplified method for absolute quantification is also proposed, in which quantitative scale is determined by DNA mass at threshold.

Kinetic PCR (kPCR) allows quantification of a target DNA within a sample, with the advantage that sensitivity is independent of copy number (1–4). The key aspect differentiating kPCR from previous quantitative PCR methodologies is that target copy number is determined from the fractional cycle at which a threshold amount of amplicon DNA is reached (threshold cycle or C_{t}), set at a point where amplicon DNA just becomes detectable, but is still within the exponential phase of the amplification (5–7). This approach ensures that interfering factors associated with the late stages of amplification are minimized, and provides the potential for unprecedented precision for quantitative determinations.

Although several methods have been developed to measure C_{t}, all are based upon fluorescent monitoring of amplicon DNA generation (8–10). Absolute quantification can be achieved using a standard curve, constructed by amplifying known amounts of target DNA in a parallel group of reactions run under identical conditions to that of the sample (7,11). Standard curve preparation is both labour intensive and error prone, with quantitative accuracy being dependent on both the accuracy of DNA standard quantification and the quality of standard curve construction (1,12,13).

In this study, a detailed examination of the mathematics governing PCR yielded insights into the process of quantitative kPCR and into some of the fundamental aspects of the threshold method. This provided a foundation from which to examine the reproducibility of standard curve construction. Assessment of intra- and inter-run variation indicates that a substantial degree of precision can be achieved, even with the application of a single standard curve to multiple runs. An alternative approach is also proposed in which quantitative scale is determined by DNA mass at threshold, such that absolute quantification would only require determination of amplification efficiency.

The DNA standard consisted of a 218 bp amplicon produced by the K3/K2 primer pair (forward K3: GGCACCTC AGGAATGGGCTATTACAA and reverse K2: AGAATA ACACAGAAATCTGTAGGTGGAATTGAA) that was purified by chloroform extraction followed by isopropanol precipitation, and quantified by averaging three replicate A_{260} absorbance determinations conducted on two spectrophotometers. A second 102 bp amplicon was produced by pairing of K2 with another primer (forward K1: TCCTATGAGATTATGACGCATTTCTCCAAA) located near the center of the K3/K2 amplicon. The primer pair combinations of K3/K2 and the nested K1/K2 thus allowed the production of two different-sized amplicons (218 and 102 bp, respectively) using the same DNA standard dilution series.

PCR amplifications were conducted using QuantiTect™ Syber® Green PCR Kit (Qiagen Inc.) according to the manufacturer’s instructions, with 0.25 µM primers and a variable amount of DNA standard in a 35 µl final reaction volume. Thermocycling was conducted using an Opticon2 DNA Engine (MJ Research Inc.) initiated by a 15 min incubation at 94°C, followed by 45 cycles (90°C, 1 s; 62°C, 120 s) with a single fluorescent reading taken at the end of each cycle. Each run was completed with a melting curve analysis to confirm the specificity of amplification and lack of primer dimers. C_{t} values were determined by the Opticon2 software using a fluorescence threshold manually set to 0.0160 for all runs and exported into a MS Excel workbook (Microsoft Inc.) for analysis (available as Supplementary Material).

The basic equation describing PCR amplification is:

N_{C} = N_{0}·(E + 1)^{C}**1**

where C is the number of thermocycles, E is amplification efficiency (also expressed as %E = E ×100%), N_{C} is the number of amplicon molecules and N_{0} is the initial number of target molecules.

In simple terms, each thermocycle produces an increase in N_{C} in proportion to amplification efficiency, such that 100% efficiency produces a doubling in the number of amplicon molecules. Additionally, the quantity of N_{C} present after any specific number of thermocycles is dependent on N_{0}. Rearrangement of equation 1 provides the mathematical relationship upon which quantitative kPCR is based:

N_{0} = N_{C}/(E + 1)^{C}**2**

Quantification of N_{C} thus allows N_{0} to be calculated if amplification efficiency is known. A major breakthrough for quantitative PCR came with the use of DNA fluorescence to monitor amplicon accumulation (3,5). Based upon this technique Higuchi *et al*. (5) developed an elegant method that simplifies N_{C} determination, such that individual amplification reactions are compared at the point at which they contain identical amounts of amplicon DNA. This is accomplished by selecting a fluorescent threshold (F_{t}) from which the fractional thermocycle (C_{t}) is calculated that defines the theoretical point at which each amplification reaction reaches fluorescence threshold.

Under this ‘threshold’ method, N_{C} becomes a constant such that equation **2** becomes:

N_{0} = N_{t}/(E + 1)^{Ct}**3**

where C_{t} is the threshold cycle and N_{t} is the number of amplicon molecules at fluorescent threshold.

Absolute quantification can be achieved using a standard curve constructed by amplification of known amounts of target DNA and plotting the resulting C_{t} values against target DNA concentration. The mathematical basis of a standard curve can be derived by taking the logarithm of equation **3**:

Log(N_{0}) = Log(N_{t}) – Log[(E + 1)^{Ct}]

Log(N_{0}) = Log (N_{t}) – Log(E + 1)·C_{t}

Log(N_{0}) = –Log(E + 1)·C_{t} + Log (N_{t})**4**

Assuming E and N_{t} are constants, equation **4** has the general structure of a line (y = mx + b) such that plotting Log(N_{0}) versus C_{t} produces a line with:

Slope = –Log(E + 1)

E_{S} = 10^{–Slope} – 1**5**

and

Intercept = Log(N_{t})

N_{t} = 10^{Intercept}**6**

where E_{S} is the slope-derived estimate of amplification efficiency.

Although the ability to derive amplification efficiency from the slope of a standard curve has been widely reported, it has not been generally recognized that the number of amplicon molecules at threshold can be directly determined from the intercept. It must also be stressed that these derivations are valid only if all PCR reactions have identical amplification efficiencies, and only if amplification efficiency is invariant over the number of thermocycles required to reach C_{t}.

Another important but often overlooked aspect of the threshold method is the interdependency of C_{t} and N_{t} on F_{t}, which has two important implications. First, C_{t} values generated from different amplification runs can be directly compared only if an identical F_{t} is used for each run. Second, the relationship between N_{t} and F_{t} is dependent on amplicon size. This is due to the fact that the underlying determinant of F_{t} is DNA fluorescence, which in turn has a linear relationship with DNA mass. As such F_{t} directly reflects DNA mass at threshold, which is related to N_{t} as described by:

M_{t} = (N_{t}·A_{S})/9.1 × 10^{11}**7**

where M_{t} is the DNA mass at threshold in nanograms, A_{S} is the amplicon size in base pairs and 9.1 × 10^{11} is the number of single base pair molecules per nanogram.

A less obvious but potentially significant extension of this is that if M_{t} is known, N_{t} can be predicted for any amplicon of known size, if it is assumed that amplicon size and base pair composition do not significantly influence DNA fluorescence. To test the general utility of PCR mathematics for standard curve evaluation and to examine the effectiveness of M_{t} for predicting N_{t}, a series of replicate standard curves was constructed for two amplicons that differ significantly in size.

Figure Figure11 is an example of the two types of graphic output generated by the instrument used in this study, and illustrates the two basic steps in quantitative kPCR using the threshold method, i.e. the selection of a fluorescent threshold from which C_{t} values are generated (Fig. (Fig.1A),1A), followed by linear regression analysis of a Log(N_{0}) versus C_{t} plot, from which E_{S} and N_{t} are estimated (Fig. (Fig.11B).

Output of a typical amplification run used for construction of the standard curves in this study. Four replicate amplification reactions were conducted for each concentration of DNA standard, covering six magnitudes of target DNA concentration. (**A**) Plot **...**

The major consideration for F_{t} selection is that it falls within the exponential phase of the amplification reaction, best illustrated by plotting log fluorescence versus cycle number (Fig. (Fig.1A).1A). As long as F_{t} is within this log-linear region, the absolute value of F_{t} was found to have only a modest impact on the slope-derived estimate of amplification efficiency (data not shown). However, as outlined above, F_{t} does have a direct impact on both C_{t} and N_{t} such that F_{t} must be fixed if data from multiple runs are to be directly compared.

To evaluate the reproducibility and quantitative variation of the threshold method, five replicate standard curves were generated from two pairs of nested primers (K3/K2 and K1/K2, see Materials and Methods for additional details) using a DNA standard dilution series covering six magnitudes of target DNA concentration. The use of nested primers allowed two different-sized amplicons (218 and 102 bp, respectively) to be amplified side-by-side within the same run, using the same DNA standard dilution series. Intra- and inter-run variation could then be examined for each of the two amplicons, free of errors caused by variations in the DNA standard. Using an identical F_{t} for all runs, the average C_{t} of four replicate amplifications for each DNA concentration were used in the analysis (Table (Table1).1). A spreadsheet containing the individual C_{t} values and the calculations used for their analysis is provided as Supplementary Material.

As an initial step for evaluation of quantitative precision, the reproducibility of amplification under our experimental conditions was estimated, based upon the standard deviation in C_{t} values generated from replicate amplifications. Moreover, due to the exponential scale of C_{t}, the impact of its variation can be difficult to assess, and thus the standard deviations in C_{t} were also used to estimate the variation in percent molecules based upon the equation:

±%Molecules = [(E + 1)^{SD} – 1] ×100%**8**

where SD is the standard deviation in C_{t} generated from replicate amplifications.

Overall the standard deviation in C_{t} of replicate amplifications ranged from 0.036 to 0.367 cycles, with an average of 0.183 cycles (Table (Table11 and Supplementary Material). This corresponds to an estimated variation in molecules that ranges from ±2.3 to ±26.6% with an average of ±12.4%, using an amplification efficiency of 90% taken from the slope-based estimate of amplification efficiency determined below.

Based upon the average standard deviation produced from each individual run, estimates of the intra-run variation were similar for both amplicons, ranging from ±9.6% to 14.9% of molecules (runs 1–5, Table Table1).1). When combined with inter-run variation, this increased to ±17.4 and ±21.3% of molecules for each amplicon, respectively, based upon averaging the standard deviation in C_{t} for each DNA concentration from all runs (‘Combined’, Table Table1).1). These variations, although significant, indicate that C_{t} values have an acceptable level of reproduciblity over the six magnitudes of target DNA concentration that were examined.

Evaluation of the quantitative variation between replicate standard curves was conducted by generating N_{t} and %E_{S} values for each amplification run listed in Table Table1.1. This was done by exporting the C_{t} values into a spreadsheet, and calculating the slope and intercept for each run using linear regression analysis of log(N_{0}) versus C_{t}. Two methods were then used to assess the quantitative variation between the five replicate standard curves constructed for each of the two amplicons (Table (Table22).

Examination of the absolute values of N_{t} and %E_{S} revealed similar trends for both amplicons, with an inter-curve variation in %E_{S} of ±2.2 and ±2.1%, and variation in N_{t} of ±19.0 and ±14.7%, respectively, as based upon their standard deviations (Table (Table2).2). Taken individually, the magnitude of variation in %E_{S} and N_{t} suggests the resulting variation in N_{0} determination could be large. For example, for a C_{t} of 25 cycles, a ±2.2% variance in the estimate of amplification efficiency would produce an approximate ±33% variation in N_{0} that, when combined with the apparent ±19% variance in N_{t}, could produce an overall variation of about ±52% for N_{0}. It must be noted, however, that further examination suggests that these estimates of variance are most certainly erroneous, due to an apparent intra-curve correlation between slope and intercept.

Comparing the %E_{S} and N_{t} values generated from each individual standard curve reveals that for both amplicons, the curve that produced the highest %E_{S} also produced the highest N_{t} (Table (Table2,2, K1/K2, run 2 and K3/K2, run 1). Similarly, the standard curves producing the lowest %E_{S} also had the lowest N_{t} (Table (Table2,2, K1/K2, run 1 and K3/K2, run 5). Taken together, these trends suggest that variations in intercept and slope are not solely caused by inter-run variation in instrumentation and/or amplification, but also reflect an innate characteristic of linear regression in which variations in slope can be compensated for to some degree by a corresponding variation in intercept.

This can be best illustrated through an alternative approach to evaluating quantitative differences produced by replicate standard curves. As illustrated in Table Table2,2, inter-curve variation can be estimated by comparing the calculated N_{0} for a series of simulated C_{t} values using equation **3**. Thus, for the five standard curves constructed from the K1/K2 amplicon, the calculated N_{0} for C_{t}=10 cycles ranges from 3.67 × 10^{7} to 4.64 × 10^{7} molecules, with an average of 4.13 × 10^{7} molecules and a standard deviation corresponding to ±8.7% of molecules (Table (Table2).2). Furthermore, a general increase in variation is observed with increasing C_{t} such that for C_{t} = 30 cycles, a variation of ±18.1% of molecules is produced. Very similar results were produced by the larger K3/K2 amplicon (Table (Table22).

Overall, this analysis demonstrates that quantitative variations produced by replicate standard curves can be relatively small, ranging in this study from a low of about ±6% to a high of about ±21% depending on the number of cycles needed to reach threshold. The observed inter-curve variation in the absolute values of slope and intercept also suggests that curve-based estimates of amplification efficiency and N_{t} require a larger data set than would normally be used for construction of a single standard curve. Indeed, the relative accuracy of the N_{t} estimates for each of the two amplicons can be tested through the correlation of their respective M_{t} values, as described by equation **7**. Based upon the N_{t} values derived from each respective ‘combined’ data set, the estimated M_{t} values differ by 7.3% (Table (Table2).2). This provides support for both the optical precision of the instrumentation and similarity in the SYBER® Green I fluorescent characteristics of these two amplicons.

Despite the extensive use of the threshold method for absolute quantification, there exists a paucity of studies that have examined the utility of the underlying mathematics. Furthermore, the general simplicity and widespread use of standard curves has led to the automation of quantitative determinations, which can obscure the mathematical principles upon which the analysis is based. Familiarity with the fundamentals of PCR mathematics cannot only yield important insights, but as well provide a foundation from which to address some of the major limitations of quantitative kPCR.

At the most basic level, the threshold method does not generally provide an effective indication of quantitative precision or accuracy. Although the standard deviation in C_{t} produced from replicate amplifications can provide an estimate of reproducibility, there is a general deficiency in reporting the errors associated with standard curve construction. This makes it difficult to evaluate the effectiveness of any specific quantitative determination, or of comparing results produced by different studies. As was demonstrated in this study, a basic assessment of standard curve construction can be conducted, if it is understood that slope and intercept are directly correlated to amplification efficiency and the number of amplicon molecules at threshold (N_{t}), respectively.

In this study, comparison of replicate standard curves revealed potentially large inter-curve variations, based upon the absolute values of slope and intercept. This initially led to the conclusion that this was caused by substantial inter-run variations in amplification and/or instrumentation. However, upon closer examination, an intra-curve correlation between slope and intercept became apparent, such that differences in slope are compensated for to a significant degree by corresponding differences in intercept.

This can be demonstrated through simple mathematical modeling, in which the initial number of target molecules (N_{0}) is calculated for a series of simulated C_{t} values (equation **3**, Table Table2).2). This showed that despite the differences in the absolute values of amplification efficiency and N_{t}, the resulting N_{0} values generated by each standard curve were unexpectedly similar. Based upon this analysis the application of a single, well-constructed standard curve could provide an estimated precision of ±6–21% of molecules, depending on the number of cycles required to reach threshold.

Notwithstanding the interrelationship of slope and intercept, it must be stressed that the mathematics of PCR dictate that amplification efficiency and N_{t} are independent entities. In reality N_{t} is determined solely by the fluorescent threshold (F_{t}), and as such its value is independent of the parameters impacting PCR amplification. Indeed, this interrelationship between N_{t} and F_{t} has important practical implications, based on the principle that F_{t} does not directly reflect the number of amplicon molecules, but rather DNA mass at fluorescent threshold (M_{t}). This in turn dictates that M_{t} could be used to predict N_{t} for any amplicon of known size, if it is assumed that amplicon size and base composition do not significantly impact DNA fluorescence. Support for the validity of this assumption was provided by the N_{t} estimates generated from the two amplicons used in this study, for which the predicted M_{t} values differ by 7.3% (equation **7**, Table Table22).

The practical significance of this becomes apparent if it is noted that N_{t} is the sole determinant of scale (equation **3**), the accuracy of which is dependent on the quantitative accuracy of the DNA standard used for standard curve construction. If, however, M_{t} can be used to predict N_{t} with sufficient precision, a common quantitative scale could be applied to all amplicons. In addition to circumventing the necessity of preparing a quantified DNA standard for each individual amplicon, the major source of variation in quantitative scale would become the optical precision of the instrument. Equally significant is that absolute quantification would be simplified, requiring only determination of amplification efficiency once M_{t} has been established.

The authors thank Richard Hamelin, Krystyna Klimaszewska and Brian Boyle for helpful comments, and Pamela Cheers for editorial assistance. This research was supported by a grant from the National Biotechnology Strategy of Canada.

1. Freeman W.M., Walker,S.J. and Vrana,K.E. (1999) Quantitative RT-PCR: pitfalls and potential. Biotechniques, 26, 112–125. [PubMed]

2. Bustin S.A. (2000) Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J. Mol. Endocrinol., 25, 169–193. [PubMed]

3. Kang J.J., Watson,R.M., Fisher,M.E., Higuchi,R., Gelfand,D.H. and Holland,M.J. (2000) Transcript quantitation in total yeast cellular RNA using kinetic PCR. Nucleic Acids Res., 28, e2. [PMC free article] [PubMed]

4. Schmittgen T.D. (2001) Real-time quantitative PCR. Methods, 25, 383–385. [PubMed]

5. Higuchi R., Fockler,C., Dollinger,G. and Watson,R. (1993) Kinetic PCR analysis: real-time monitoring of DNA amplification reactions. Biotechnology, 11, 1026–1030. [PubMed]

6. Morrison T.B., Weis,J.J. and Wittwer,C.T. (1998) Quantification of low-copy transcripts by continuous SYBR® Green I monitoring during amplification. Biotechniques, 24, 954–962. [PubMed]

7. Rasmussen R. (2001) Quantification on the LightCycler. In Meuer,S., Wittwer,C. and Nakagawara,K. (eds), Rapid Cycle Real-time PCR, Methods and Applications. Springer Press, Heidelberg, pp. 21–34.

8. Wittwer C.T., Herrmann,M.G., Moss,A.A. and Rasmussen,R.P. (1997) Continuous fluorescence monitoring of rapid cycle DNA amplification. Biotechniques, 22, 130–138. [PubMed]

9. Wittwer C., Ririe,K. and Rasmussen,R. (1998) Fluorescence monitoring of rapid cycle PCR for quantification. In Ferré,F. (ed.), Gene Quantification. Birkhäuser, Boston, MA, pp. 129–144.

10. Giulietti A., Overbergh,L., Valckx,D., Decallonne,B., Bouillon,R. and Mathieu,C. (2001) An overview of real-time quantitative PCR: applications to quantify cytokine gene expression. Methods, 25, 386–401. [PubMed]

11. Pfaffl M.W. and Hageleit,M. (2001) Validities of mRNA quantification using recombinant RNA and recombinant DNA external calibration curves in real-time RT-PCR. Biotechnol. Lett., 23, 275–282.

12. Niesters H.G.M. (2001) Quantitation of viral load using real-time amplification techniques. Methods, 25, 419–429. [PubMed]

13. Lehmann U. and Kreipe,H. (2001) Real-time PCR analysis of DNA and RNA extracted from formalin-fixed and paraffin-embedded biopsies. Methods, 25, 409–418. [PubMed]

Articles from Nucleic Acids Research are provided here courtesy of **Oxford University Press**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |