Published online 15 December 2004
Nucleic Acids Research, Vol. 32 No. 22 © Oxford University Press 2004; all rights reserved
Sigmoidal curve-fitting redefines quantitative real-time PCR with the prospective of developing automated high-throughput applications
Natural Resources Canada, 1055 du P.E.P.S, Sainte-Foy, Quebec, Canada G1V 4C7
* Tel: +1 418 648 2582; Fax: +1 418 648 5849; Email: Bob.Rutledge{at}NRCan.gc.ca
Received May 25, 2004; Revised October 15, 2004; Accepted November 24, 2004
| ABSTRACT |
|---|
|
|
|---|
Quantitative real-time PCR has revolutionized many aspects of genetic research, biomedical diagnostics and pathogen detection. Nevertheless, the full potential of this technology has yet to be realized, primarily due to the limitations of the threshold-based methodologies that are currently used for quantitative analysis. Prone to errors caused by variations in reaction preparation and amplification conditions, these approaches necessitate construction of standard curves for each target sequence, significantly limiting the development of high-throughput applications that demand substantive levels of reliability and automation. In this study, an alternative approach based upon fitting of fluorescence data to a four-parametric sigmoid function is shown to dramatically increase both the utility and reliability of quantitative real-time PCR. By mathematically modeling individual amplification reactions, quantification can be achieved without the use of standard curves and without prior knowledge of amplification efficiency. Combined with provision of quantitative scale via optical calibration, sigmoidal curve-fitting could confer the capability for fully automated quantification of nucleic acids with unparalleled accuracy and reliability.
| INTRODUCTION |
|---|
|
|
|---|
First introduced commercially in 1996, fluorescence-based detection of amplicon DNA allowed the kinetics of PCR amplification to be monitored in real time, providing the ability to quantify nucleic acids with extraordinary ease and accuracy (13). With a large dynamic range (78 magnitudes) and a high degree of sensitivity (15 molecules), quantitative real-time PCR has greatly impacted many aspects of genetic research, in addition to facilitating the development of new applications in biomedical diagnostics for pathogen detection, viral load, and minimal residual disease, and for gene expression analysis in relation to disease, pathogenesis and oncogenesis (49).
Despite the significance of this technology, its full potential has yet to be realized, primarily due to limitations associated with the threshold-based methodologies that currently predominate. Developed upon comparing amplification reactions at a point in which they have identical amounts of amplicon DNA, a theoretical threshold cycle (Ct) is the primary quantitative output (1012). Much of the technical limitation of the threshold approach is associated with conversion of Ct values into the number of target molecules that requires construction of a standard curve for each target sequence (1214). Although standard curves can be effective, their construction requires extensive effort and is prone to errors in DNA standard preparation, which is further exacerbated by the difficulties of assessing accuracy-of-scale (5,15,16).
An alternative approach was recently proposed by Liu and Saint (17), in which a sigmoid function was shown to model PCR amplification more effectively than the exponential model upon which the threshold method is based. In the study presented here, the remarkable implications and quantitative effectiveness of sigmoidal modeling were examined in detail, which demonstrates that many of the limitations of the threshold method can be circumvented, including abrogating the need for standard curves. Most significant, however, is the potential to develop automated high-throughput applications for quantitative real-time PCR, with unprecedented capabilities for quality control over both PCR amplification and establishment of quantitative scale.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Fluorescence datasets
All the fluorescence readings used in this study were obtained from a previous study, in which five replicate standard curves were constructed for each of two amplicons (K3/K2, 218 bp and K1/K2, 102 bp), with the smaller amplicon nested within the larger (18). In brief, PCR amplifications were conducted using QuantiTectTM SYBR® Green PCR Kit (Qiagen Inc.) according to the manufacturer's instructions, and an Opticon2 DNA Engine (MJ Research Inc.) with thermocycling initiated by a 15 min incubation at 94°C, followed by 45 cycles (90°C for 1 s; 62°C for 120 s) with a single fluorescent reading taken at the end of each cycle.
Specificity of amplification and absence of primer dimers was confirmed by melting curve analysis at the end of each run. To maintain consistency with the previously generated Ct data, fluorescence readings were exported after background subtraction via baseline averaging of the 5 cycles immediately preceding the cycles in which fluorescence was first detected, and data for each amplicon were exported to an Excel workbook for analysis (provided as Supplementary Datasets 1 and 2).
Curve-fitting
The nonlinear regression function of SigmaPlot (Version 8) was used to fit fluorescence readings to Equation 1. Following extensive analyses, it became evident that cycles within the plateau phase diverged significantly from that of predicted by sigmoidal modeling, an anomaly that impacted the effectiveness of the curve-fitting process. It was thus necessary to exclude plateau cycles, a process that was based upon selection of a cut-off cycle, beyond which cycles were excluded from the regression analysis (Figure 1; Supplementary Datasets 1 and 2). The criterion used for the selection of the cut-off cycle was developed upon subjecting each amplification curve to repetitive regression analyses, in which the last cycle was excluded and the regression analysis repeated. As best illustrated in the Excel workbooks provided as Supplementary Datasets 1 and 2, this generated a series of r2, Fmax, C1/2 and k values, from which a F0 value was calculated for each sequential data-subset using Equation 2. The impact of this sequential exclusion of cycles was then evaluated, which revealed a highly regular trend in which the calculated F0 value reached a minimum, followed by a small but progressive increase (Figure 2). Selection of the cut-off cycle for all amplification curves was based upon the subset that produced the minimum-calculated F0 value (for additional details see Supplementary Datasets 1 and 2).
|
|
| RESULTS |
|---|
|
|
|---|
Sigmoidal curve-fitting (SCF)
As previously described by Liu and Saint (17), real-time PCR can be effectively modeled using the four-parametric sigmoid function:
![]() | (1) |
The general efficacy of SCF for modeling of real-time PCR amplification can be demonstrated by the correlation of experimentally derived FC datasets to that predicted by Equation 1. Indeed, extensive analyses confirmed that SCF is effective, producing correlation coefficients (r2) that reach averages of about 0.99998 for highly replicated amplifications (n = 20), down to an average of about 0.99992 for individual amplification reactions (Supplementary Datasets 1 and 2; see below). An important exception, however, is evident in Figure 1, in that the reaction fluorescence produced within the plateau phase was found to diverge from that predicted by Equation 1. Although this anomaly is much less prevalent under high amplification efficiencies (data not shown), it was found necessary to exclude plateau cycles from the curve-fitting process through the selection of a cut-off cycle (Figure 1). Nonetheless, extensive application of SCF demonstrated a capacity to model real-time PCR with a truly remarkable degree of precision.
Of all the challenges encountered during development of the SCF method, developing a criterion for the selection of the cut-off cycle was the most difficult. An effective solution was based upon the empirical observation that as cycles from the plateau were sequentially removed from the curve-fitting process, the resulting changes in F0 value were found to be highly reproducible (Figure 2). Selection of the cut-off cycle for each target concentration was subsequently based upon that which produced a minimum F0 value, as illustrated in Figure 2 (see Materials and Methods, and Supplementary Datasets 1 and 2 for additional details).
Quantitative precision of SCF
Previous quantitative methods have been based upon the direct relationship between the target quantity and the number of cycles needed to produce quantifiable amounts of amplicon DNA. An alternative approach is provided by a simple derivative of Equation 1, when C = 0:
![]() | (2) |
The SCF method presents two aspects of practical importance. First is the potential to conduct quantification without the knowledge of amplification efficiency, and second, that quantification and establishment of quantitative scale are separable processes. Thus, if sufficient levels of quantitative precision can be demonstrated, and if establishment of quantitative scale (i.e. fluorescence calibration) can be effectively achieved, SCF could provide capabilities either difficult to achieve or unattainable using a threshold approach.
A previous study in which the quantitative accuracy of the Ct method was examined provided an ideal dataset for evaluating the SCF method (18). Encompassing 10 standard curves constructed from replicate amplifications (n = 4) of six target concentrations, this earlier study generated a total of 120 amplification reactions for each of the two amplicons (K1/K2 = 102 bp and K3/K2 = 218 bp). In addition to providing a large number of amplification reactions from which fluorescence datasets could be obtained, this provided the ability to directly compare quantitative determinations generated by SCF and Ct methodologies. Owing to the large size of the datasets and extensive scope of the work, it is only possible to present summaries of the analyses. The reader is, however, referred to the Excel workbooks provided as supplementary data online (Supplementary Datasets 1 and 2), which contain the fluorescence readings along with the calculations and data compilations used to prepare the summaries presented in this report.
As a starting point for assessing quantitative precision, the relative differences in F0 values derived from each of the six target quantities were examined for each of the two amplicons (Table 1). This indicated that a high level of precision can be achieved with SCF quantification, as reflected by the low variance from that of predicted (Table 1, ±3.0% and ±7.7% for K3/K2 and K1/K2, respectively), and that these are very similar to the variances produced by Ct-based quantification (±4.8% and ±10.0%, respectively, Supplementary Datasets 1 and 2).
|
The quantitative precision of SCF was further assessed by examining the relationship of F0 to amplicon size, as based on the linear relationship between SYBR® Green I fluorescence and DNA mass; i.e. for an identical target quantity, the respective F0 values for each amplicon should differ in direct relation to their size (218 bp versus 102 bp; ratio = 2.14). Indeed, the average F0 values adjusted to the highest target quantity differ from theoretical value by only 6.8% (5.91 x 105 versus 2.95 x 105, Table 1; ratio = 2.00). In addition to supporting the precision of SCF quantification, this further confirms that the fluorescence characteristics of these two amplicons are similar, as was previously concluded based upon their respective DNA mass at threshold (Mt, see below) that differ by 7.3% (18).
Establishment of quantitative scale via optical calibration
Conversion of F0 to the number of target molecules would appear to be straightforward, accomplished by first correlating fluorescence to DNA mass such that:
![]() | (3) |
![]() | (4) |
Amplification of a quantified DNA standard provides a straightforward method for CF determination, based upon rearrangement of Equation 3:
![]() | (5) |
![]() | (6) |
is the Ct-based calibration factor, Mt is amplicon DNA mass at threshold and Ft is the fluorescence threshold, in which Mt is derived from the number of amplicon molecules at threshold (Nt), which in turn is derived from the intercept of a Ct-based standard curve (18).
The relative accuracies of SCF- and Ct-based quantification were assessed by comparing CFSCF and
generated for each of the two amplicons (Table 2). This revealed that each of the four quantitative scales agree within ±5%, supporting the effectiveness of these two diverse quantitative methods in relation to a shared DNA standard. The effectiveness of a common quantitative scale was further examined by converting F0 and Ct values into the number of target molecules using the average CF derived in Table 2. This demonstrated that the combined variance in N0 determination for both amplicons averaged ±6.4% from that of predicted over the six magnitudes of target quantity examined (Table 3).
|
|
Dynamics of amplification efficiency
A fundamental difference between exponential and sigmoidal models relates to the amplification efficiency, such that for threshold-based quantification, amplification efficiency is presumed to be constant (14,18), whereas under SCF quantification, amplification efficiency is dynamic. This can be illustrated with an equation described by Liu and Saint (17):
![]() | (7) |
![]() | (8) |
Relating exponential and sigmoidal models
Quantification under the threshold method is based upon the exponential equation:
![]() | (9) |
![]() | (10) |
Analogy of Ct with C1/2, can be demonstrated by plotting Log(N0) against C1/2, which also produces a line (Figure 3) such that:
![]() | (11) |
![]() | (12) |
![]() | (13) |
|
This provides mathematical support for the analogy of the initial amplification efficiency (E0) as is defined by Equation 8, to the slope-derived efficiency (ES) as defined by a Ct-based standard curve (18). Indeed, this contention is borne out by the close correlation observed between E0 and ES values derived for each amplicon (Table 4; Supplementary Datasets 1 and 2).
|
Analogy with a Ct-based quantification can be extended further by substituting Nmax for Nt in Equation 9, and converting amplicon number to reaction fluorescence:
![]() | (14) |
![]() | (15) |
![]() | (16) |
| DISCUSSION |
|---|
|
|
|---|
Real-time PCR is well suited for high-throughput quantification of nucleic acids that could greatly benefit a variety of applications in medical diagnostics and genomics research, to name only two prominent fields. Large-capacity microtiter plates, combined with software automation of data collection and processing, provide the fundamental capabilities required for high-volume analysis. Nevertheless, high-throughput quantification is still technically difficult to achieve, primarily due to deficiencies of the threshold-based methodologies that have prevailed since the introduction of real-time PCR (1).
Most of these limitations center on the requirement for a reliable estimation of amplification efficiency that is generally acquired through construction of a standard curve, a process that is repeated for each target sequence. The technical challenges are not only exacerbated by difficulties of preparing quantified standards for each target, but also by the necessary assumption that the amplification efficiency of samples is identical, or at least similar, to that predicted by the standard curve. Such an assumption has been reported to be patently invalid for many cases in which amplification efficiency of samples has been determined (21,22). Not only can small differences in amplification efficiency produce large quantitative errors, the frequency and magnitude of these errors are virtually impossible to ascertain using a threshold approach.
The study described here extends the pioneering work of Liu and Saint (17), demonstrating that SCF can circumvent many of these deficiencies. The outstanding capabilities of SCF center on mathematic modeling of individual amplifications, which requires no prior assumptions other than dependency on the optical precision of the instrumentation. From this work, three aspects of practical significance are evident: quantification without prior knowledge of amplification efficiency, assessment of amplification efficacy and establishment of quantitative scale.
The effectiveness of SCF-based quantification is fundamentally linked to curve-fitting of experimental data, such that variations unique to each amplification reaction are incorporated into the analysis. Thus, instead of deriving a single quantitative entity, as is the case for threshold methodologies, SCF generates values for three kinetic parameters from which target quantity is calculated. Although lacking scale, the resulting target quantity could be directly used for relative quantification (16,23,24), with the advantage that determination of amplification efficiency is unnecessary.
This is not to say that amplification efficiency is unimportant, nor is it necessarily difficult to determine via SCF; indeed, amplification efficiency and target quantity are simply different expressions of the same kinetics parameters. Consequently, SCF allows the routine evaluation of amplification efficacy within individual reactions via monitoring of the initial amplification efficiency. In addition to identification of aberrant reaction preparation, this could allow detection of enzymatic inhibitors within the sample. Although effective quality control is clearly desirable for high-throughput applications in research, it is of even greater importance to biomedical diagnostics where reliability of analysis is paramount. It is equally evident that accuracy-of-scale is critical to many diagnostic applications, particularly for those involving pathogenesis and residual disease (5).
Of the many issues encompassing quantitative real-time PCR, provision of quantitative scale has some of the most profound implications, and yet it is one of the least understood. The simplicity of concept provided by standard curves employed under the threshold approach has likely contributed to the general oversight that quantitative scale is directly linked to amplicon fluorescence through its relationship to DNA mass. Based upon this principle, it was previously proposed that a common quantitative scale could be established for Ct-based analyses via determination of amplicon mass at threshold (18), a concept that was utilized further in this study. The relationship of quantitative scale to amplicon fluorescence is more evident under SCF in that target quantification is derived in fluorescence units. Determination of target copy number thus simply requires a calibration factor that relates fluorescence to DNA mass. It is the method for establishing this calibration factor, and the general effectiveness of applying a single calibration factor to multiple amplicons, that have major implications for the accuracy and reliability that can be achieved.
In this regard, it is important to note that optical calibration determines the absolute accuracy or exactness of quantitative scale, whereas the curve-fitting process determines the precision of the assay (25). Although subtle, this is an important distinction in that acquisition of quantitative scale may be linked to PCR amplification (and thus susceptible to the factors impacting the precision of PCR quantification), but not necessarily so. Alternative methods for optical calibration can be envisaged that are independent of PCR, similar, for example, to the standard curves generated for DNA quantification using a dedicated fluorometer.
Equally significant is that once a calibration factor is established, it can provide a common quantitative scale for any target sequence, if it is assumed that the base pair composition and amplicon size do not significantly impact the fluorescence characteristics of SYBR® Green I. Under such an approach, the need for preparing a quantified DNA standard for each individual target sequence would be circumvented, and would allow large numbers of different targets to be quantified simultaneously, based upon a single, pre-established quantitative scale. Ultimately, it may be possible to introduce a gold standard for the establishment of quantitative scale, which would further increase the utility and accuracy of quantitative real-time PCR.
In conclusion, SCF provides many of the fundamental capabilities required for fully automated high-throughput quantification that are lacking under currently employed methodologies. These include routine assessment of amplification efficacy within individual amplification reactions and elimination of PCR-generated standards curves, which when combined with provision of quantitative scale via optical calibration could allow reliable high-throughput quantification with minimal intervention. Based upon redefining the mathematics of PCR, sigmoidal modeling also elaborates more clearly on the theoretical aspects of PCR amplification upon which to base further investigations into the factors impacting amplification kinetics, and to more effectively expand the application of quantitative real-time PCR.
| SUPPLEMENTARY MATERIAL |
|---|
|
|
|---|
Supplementary Material is available at NAR Online. Supplementary Dataset 1: K1/K2 amplicon Excel workbook containing the fluorescence reading, SCF analysis and data summaries of the K1/K2 amplifications. Supplementary Dataset 2: K3/K2 amplicon Excel workbook containing the fluorescence reading, SCF analysis and data summaries of the K3/K2 amplifications.
| ACKNOWLEDGEMENTS |
|---|
The author thanks Brian Boyle and Richard Hamelin for their critical review of the manuscript, and Isabelle Lamarre for editorial assistance. This research was supported by a grant from the National Biotechnology Strategy of Canada.
| REFERENCES |
|---|
|
|
|---|
- Higuchi,R., Fockler,C., Dollinger,G. and Watson,R. ( (1993) ) Kinetic PCR analysis: real-time monitoring of DNA amplification reactions. Biotechnology, , 11, , 10261030.[CrossRef][Medline]
- 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, , 954962.[Web of Science][Medline]
- 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.
[Abstract/Free Full Text] - Bustin,S.A. ( (2000) ) Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J. Mol. Endocrinol., , 25, , 169193.[Abstract]
- Niesters,H.G.M. ( (2001) ) Quantitation of viral load using real-time amplification techniques. Methods, , 25, , 419429.[CrossRef][Web of Science][Medline]
- Mackay,I.M., Arden,K.E. and Nitsche,A. ( (2002) ) Real-time PCR in virology. Nucleic Acids Res., , 30, , 12921305.
[Abstract/Free Full Text] - Ivnitski,D., O'Neil,D.J., Gattuso,A., Schlicht,R., Calidonna,M. and Fisher,R. ( (2003) ) Nucleic acid approaches for detection and identification of biological warfare and infectious disease agents. BioTechniques, , 35, , 862869.[Web of Science][Medline]
- Klein,D. ( (2002) ) Quantification using real-time PCR technology: applications and limitations. Trends Mol. Med., , 8, , 257260.[CrossRef][Web of Science][Medline]
- Bustin,S.A. ( (2002) ) Quantification of mRNA using real-time reverse transcription PCR (RTPCR): trends and problems. J. Mol. Endocrinol., , 29, , 2339.[Abstract]
- Wittwer,C.T., Herrmann,M.G., Moss,A.A. and Rasmussen,R.P. ( (1997) ) Continuous fluorescence monitoring of rapid cycle DNA amplification. BioTechniques, , 22, , 130138.[Web of Science][Medline]
- 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. 129144.
- Freeman,W.M., Walker,S.J. and Vrana,K.E. ( (1999) ) Quantitative RTPCR: pitfalls and potential. BioTechniques, , 26, , 112125.[Web of Science][Medline]
- Pfaffl,M.W. and Hageleit,M. ( (2001) ) Validities of mRNA quantification using recombinant RNA and recombinant DNA external calibration curves in real-time RTPCR. Biotechnol. Lett., , 23, , 275282.[CrossRef]
- Rasmussen,R. ( (2001) ) Quantification on the LightCycler. In Meuer,S., Wittwer,C. and Nakagawara,K. (eds), Rapid Cycle RealTime PCR: Methods and Applications. Springer Press, Heidelberg, pp. 2134.
- Lehmann,U. and Kreipe,H. ( (2001) ) Real-time PCR analysis of DNA and RNA extracted from formalin-fixed and paraffin-embedded biopsies. Methods, , 25, , 409418.[CrossRef][Web of Science][Medline]
- Peirson,S.N., Butler,J.N. and Foster,R.G. ( (2003) ) Experimental validation of novel and conventional approaches to quantitative real-time PCR data analysis. Nucleic Acids Res., , 31, , e73.
[Abstract/Free Full Text] - Liu,W. and Saint,D.A. ( (2002) ) Validation of a quantitative method for real time PCR kinetics. Biochem. Biophys. Res. Commun., , 294, , 347353.[CrossRef][Web of Science][Medline]
- Rutledge,R.G. and Côté,C. ( (2003) ) Mathematics of quantitative kinetic PCR and the application of standard curves. Nucleic Acids Res., , 31, , e93.
[Abstract/Free Full Text] - Tichopad,A., Dzidic,A. and Pfaffl,M.W. ( (2002) ) Improving quantitative real-time RTPCR reproducibility by boosting primer-linked amplification efficiency. Biotechnol. Lett., , 24, , 20532056.[CrossRef]
- Wilhelm,J., Pingoud,A. and Hahn,M. ( (2003) ) Validation of an algorithm for automatic quantification of nucleic acid copy numbers by real-time polymerase chain reaction. Anal. Biochem., , 317, , 218225.[CrossRef][Web of Science][Medline]
- Meijerink,J., Mandigers,C., van de Locht,L., Tönnissen,E., Goodsaid,F. and Raemaekers,J. ( (2001) ) A novel method to compensate for different amplification efficiencies between patient DNA samples in quantitative real-time PCR. J. Mol. Diagn., , 3, , 5561.
[Abstract/Free Full Text] - Ståhlberg,A., Åman,P., Ridell,B., Mostad,P. and Kubista,M. ( (2003) ) Quantitative real-time PCR method for detection of B-lymphocyte monoclonality by comparison of kappa and lambda immunoglobulin light chain expression. Clin. Chem., , 49, , 5159.
[Abstract/Free Full Text] - Livak,K.J. and Schmittgen,T.D. ( (2001) ) Analysis of relative gene expression data using real-time quantitative PCR and the 2

Ct method. Methods, , 25, , 402408.[CrossRef][Web of Science][Medline]
- Pfaffl,M.W. ( (2001) ) A new mathematical model for relative quantification in real-time RTPCR. Nucleic Acids Res., , 29, , e45.
[Abstract/Free Full Text] - Ferré,F. ( (1998) ) Key issues, challenges, and future opportunities in gene quantification. In Ferré,F. (ed.), Gene Quantification. Birkhäuser, Boston, MA, pp. 116.
This article has been cited by other articles:
![]() |
J. M. Ruijter, C. Ramakers, W. M. H. Hoogaars, Y. Karlen, O. Bakker, M. J. B. van den Hoff, and A. F. M. Moorman Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data Nucleic Acids Res., April 1, 2009; 37(6): e45 - e45. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Ritz and A.-N. Spiess qpcR: an R package for sigmoidal model selection in quantitative real-time polymerase chain reaction analysis Bioinformatics, July 1, 2008; 24(13): 1549 - 1551. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Neuzil, C. Zhang, J. Pipper, S. Oh, and L. Zhuo Ultra fast miniaturized real-time PCR: 40 cycles in less than six minutes Nucleic Acids Res., June 28, 2006; 34(11): e77 - e77. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




















