- Methodology article
- Open Access
- Published:

# Critical evaluation of methods used to determine amplification efficiency refutes the exponential character of real-time PCR

*BMC Molecular Biology*
**volume 9**, Article number: 96 (2008)

## Abstract

### Background

The challenge of determining amplification efficiency has long been a predominant aspect of implementing real-time qPCR, playing a critical role in the accuracy and reliability that can be achieved. Based upon analysis of amplification profile position, standard curves are currently the gold standard for amplification efficiency determination. However, in addition to being highly resource intensive, the efficacy of this approach is limited by the necessary assumption that all samples are amplified with the same efficiency as predicted by a standard curve. These limitations have driven efforts to develop methods for determining amplification efficiency by analyzing the fluorescence readings from individual amplification reactions. The most prominent approach is based on analysis of the "log-linear region", founded upon the presumption that amplification efficiency is constant within this region. Nevertheless, a recently developed sigmoidal model has provided new insights that challenge such historically held views, dictating that amplification efficiency is not only dynamic, but is linearly coupled to amplicon DNA quantity. Called "linear regression of efficiency" or LRE, this kinetic-based approach redefines amplification efficiency as the maximal efficiency (E_{max}) generated at the onset of thermocycling.

### Results

This study presents a critical evaluation of amplification efficiency determination, which reveals that potentially large underestimations occur when exponential mathematics is applied to the log-linear region. This discrepancy was found to stem from misinterpreting the origin of the log-linear region, which is derived not from an invariant amplification efficiency, but rather from an exponential loss in amplification rate. In contrast, LRE analysis generated E_{max} estimates that correlated closely to that derived from a standard curve, despite the fact that standard curve analysis is founded upon exponential mathematics. This paradoxical result implies that the quantitative efficacy of positional-based analysis relies not upon the exponential character of real-time PCR, but instead on the ability to precisely define the relative position of an amplification profile.

### Conclusion

In addition to presenting insights into the sigmoidal character of the polymerase chain reaction, LRE analysis provides a viable alternative to standard curves for amplification efficiency determination, based on analysis of high-quality fluorescence readings within the central region of SYBR Green I generated amplification profiles.

## Background

All commercial real-time quantitative PCR platforms currently rely on defining the relative position of amplification profiles. As such, they are reliant on amplification of a serially diluted target to provide an estimate of amplification efficiency, which is essential to accurate and reliable quantification [1, 2]. However, a major caveat of this approach is that sample-specific inhibitors can compromise both the reliability and accuracy of an assay. This can be a major concern, particularly for samples derived from sources known to contain PCR inhibitors, in that any loss of amplification efficiency will generate unidentified and potentially large quantitative errors [3, 4].

Driven by the highly desirable ability to assess amplification kinetics within individual PCR reactions, a large number of studies have attempted to utilize the kinetics of amplicon DNA accumulation as the basis for amplification efficiency determination, through the application of exponential mathematics [5–15]. Derived from the intuitive notion that PCR amplification is inherently exponential in nature, this approach attempts to exploit the presence of a "log-linear region" within the lower region of real-time amplification profiles. Founded on the presumption that log-linearity reflects constant amplification efficiency, amplification efficiency is calculated from the slope of the log-linear region, similar to that conducted for standard curves.

An alternative approach to fluorescence analysis is based upon the presumption that PCR amplification is inherently sigmoidal in nature, allowing amplification efficiency to be estimated by fitting fluorescence readings to the classic Boltzmann sigmoid function using nonlinear regression analysis [16–26]. Unfortunately, effective implementation of this approach has been impeded by errors produced by distortions within the plateau phase [18, 23–25]. Recent development of a simplified approach to sigmoidal analysis based on recognition that amplification rate is linearly coupled to amplicon DNA quantity, circumvents such plateau phase anomalies. Called "linear regression of efficiency" or LRE analysis, amplification efficiency is determined by applying linear regression analysis to the fluorescence readings within the central region of an amplification profile [27].

Under the LRE model, amplification efficiency is maximal at the onset of thermocycling, with amplification rate progressively slowing as amplicon DNA accumulates, such that each cycle produces a unique amplification efficiency, with entry into the plateau phase occurring as amplification efficiency approaches zero. As such, amplification efficiency is defined as the maximal efficiency (E_{max}) produced in the absence of amplicon DNA. Although this clearly conflicts with the exponential model of real-time qPCR, which dictates that amplification efficiency is constant, it is unclear as to what extent these opposing interpretations impact the efficacy of amplification efficiency determination.

A central objective of this study was to critically evaluate exponential- and sigmoidal-based fluorescence analysis for determining amplification efficiency, with the expectation that standard curves would provide a gold standard from which to base the comparison. Notwithstanding the extensive resources required for their construction, the positional analysis upon which standard curves are based proved to be an effective platform from which to conduct the analysis.

## Results and discussion

### C_{t}-based standard curve analysis

The current gold standard for conducting amplification efficiency determinations is based on analysis of a serially diluted target. An example of this approach is presented in Figure 1 in which lambda gDNA is diluted in 10-fold increments to cover a quantitative range of five magnitudes. Examination of the resulting amplification profiles illustrates the central principle underpinning real-time qPCR, which is that profile position is precisely related to target quantity, as reflected by the regular spacing of these amplification profiles (Figure 1A). A prominent approach to positional analysis, called the threshold method, exploits this principle by defining profile position in terms of the fraction cycle (called the threshold cycle or C_{t}) at which reaction fluorescence reaches an arbitrary quantity (called the fluorescence threshold or F_{t}). Target quantification can then be achieved by analyzing the position of the amplification profile generated by a sample, in relation to that produced by a serially diluted, target-specific quantitative standard.

The mathematics of the threshold method is relatively simple, originating from the premise that PCR amplification is exponential. As such, target quantity can be expressed in terms of the threshold and amplification efficiency as:

where *N*_{
0
}is the number of target molecules, *N*_{
t
}is the number of amplicon molecules at threshold and *E* is the amplification efficiency. Target quantity can thus be calculated once values for *N*_{
t
}and *E* have been obtained, which can be accomplished by constructing a standard curve in which Log(N_{0}) is plotted against C_{t}, such that amplification efficiency is defined by the slope:

*E*
*slope*
_{= 10}
*slope*
^{- 1}

where *E**slope*is the slope-derived estimate of amplification efficiency [2]. Figure 1B illustrates this approach, which produced an *E*_{
slope
}of 92.0% with a linear correlation coefficient (r^{2}) of 0.9997.

Such a high level of linearity is consistent with the contention that PCR amplification is exponential in nature, which in turn supports the presumption that amplification efficiency is invariant. Having said this, however, it is also evident from the sigmoidal shape of a typical SYBR Green I real-time profile that amplification efficiency is not completely invariant, as reflected by the eventual cessation of amplification that defines the end point of the amplification process, commonly referred to as the "plateau phase" (Figure 1A). This potential contradiction has led many prominent studies to surmise that an "exponential region" exists within the lower domain of amplification profiles [6, 9, 10, 13–15, 20, 28–30], with the inference that the upper boundary of this exponential region is defined as the point at which the amplification efficiency begins to decrease.

The general validity of this presumption can be tested empirically by compiling a series of standard curves generated by progressively increasing F_{t}. This also allows testing of another general presumption, which is that the integrity of a standard curve relies on placing F_{t} within this putative exponential region [14, 15, 28, 29]. As illustrated in Figure 2, this appears not to be the case, in that E_{slope} was unaffected until the fluorescence threshold was placed into the extreme upper region of the profiles.

Not only does this contest the existence of an exponential region, it also challenges the perception that the efficacy of the threshold method is reliant on the exponential character of real-time amplification profiles. Partial resolution to this apparent paradox can be gained by considering PCR amplification to be an inherently sigmoidal process in which amplification efficiency is dynamic [18, 27]. Under such a scenario, it could be surmised that the efficacy of the threshold method is derived from the application of a purely positional-based analysis. A slope-based estimate of amplification efficiency (E_{slope}) would thus reflect the rate of amplification, not at C_{t}, but rather the maximal amplification efficiency (E_{max}) as defined under a sigmoidal paradigm (see below). Indeed, sigmoidal analysis provides insights that support such an interpretation.

### Kinetic-based sigmoidal analysis

Recognition that amplification rate is linearly correlated with amplicon quantity led to the proposal that the dynamics of PCR amplification can be described by the linear equation:

*E*_{
C
}= *ΔE* × *F*_{
C
}+*E*_{max}

where *E*_{
C
}is the amplification efficiency produced during cycle C (referred to as "cycle efficiency"), *F*_{
C
}is the reaction fluorescence at cycle C and is proportional to the mass of amplicon DNA in the reaction, *ΔE* is the rate of loss in amplification efficiency and *E*_{
max
}is the maximal amplification efficiency, that is, when *F*_{
C
}= 0 [27]. Called "linear regression of efficiency" or LRE analysis, this approach is implemented by first estimating cycle efficiency from the relative increase in reaction fluorescence over each individual cycle:

where *F*_{
C-1
}is the fluorescence reading produced by the previous cycle. Plotting E_{C} against F_{C} produces a line as described by equation 3, from which *E*_{
max
}is determined from the intercept. Application of LRE analysis to the amplification profiles presented in Figure 1 produced an average E_{max} of 95.5% ± 0.9% (Figure 3), as compared with the 92.0–92.6% E_{slope} estimates produced from the standard curve analysis (Figures 1B and 2B).

In addition to greatly simplifying amplification efficiency determination, development of LRE analysis was also instrumental to the derivation of two equations that allow real-time PCR amplification to be mathematically modeled. As described previously [27], the LRE model was developed by using *ΔE* and *E*_{
max
}to adapt the classic Boltzmann sigmoid function to PCR. Although the primary application of LRE modeling is target quantification, it also allows PCR amplification to be predicted with a very high degree of precision, as is illustrated in Figure 4.

The LRE model can also be used to further evaluate the positional analysis upon which the threshold method is based, through the ability to describe profile position via the fraction cycle, called C_{1/2}, that defines the mid-point at which reaction fluorescence reaches precisely half of its maxima [27]. C_{1/2} can be calculated for an individual amplification profile once values for *ΔE*, *E*_{
max
}and average *F*_{
0
}have been obtained, using a derivative generated during development of the LRE method (Figure 5A).

Although C_{1/2} may appear to be similar to C_{t}, a central attribute of C_{1/2} is that it is purely sigmoidal in origin. This then provides the opportunity to examine the efficacy of positional analysis by constructing a standard curve using the C_{1/2} values from each of the amplification profiles presented in Figure 1A. As presented in Figure 5B, this produces a slope-derived amplification efficiency estimate of 92.6%, a value that falls within the range of estimates of 92.0–92.6% produced by the C_{t}-based standard curve analyses (Figures 1B and 2B). This further supports the contention that the efficacy of the threshold method is dependent on accurately defining profile position, rather than on the exponential nature of real-time PCR.

### Exponential analysis of the log-linear region

First introduced during the early history of quantitative PCR, a number of studies have attempted to determine amplification efficiency by exploiting the existence of a "log-linear region" within the lower region of an amplification profile [5–9]. Differing primarily in the number and position of the cycles included in the analysis, several implementations of this approach have been reported for real-time qPCR [9–15], which have subsequently been utilized for large-scale gene expression profiling in which standard curve construction is impractical [20, 31, 32]. Similar to the exponential mathematics of the threshold method, amplification efficiency is calculated from the slope of the log-linear region:

*E*_{
loglin
}= 10^{slope} - 1

where *E*_{
loglin
}is the log-linear estimate of amplification efficiency. It is important to note, however, that the validity of this approach rests on the supposition that amplification efficiency remains constant throughout the log-linear region. Indeed, the very existence of a log-linear region lends itself to the compelling, albeit implicit, contention that amplification efficiency is invariant within this region. Nevertheless, critical examination reveals that linearity alone is insufficient to validate this contention. Figure 6 provides an example of this approach, based upon analysis of the five amplification profiles presented in Figure 1.

This generated an average r^{2} of 0.9991, corroborating the presence of a substantive log-linear region within the lower region of all five amplification profiles. Nevertheless, the first indication of an anomaly is provided by the E_{loglin} values produced by the five profiles, which averaged 80.3% ± 2.5% (Figure 6). This is substantially lower than the C_{t}-E_{slope} of 92.0–92.6%, E_{max} of 95.5%, and C_{1/2}-E_{slope} of 92.6% generated from the same five amplification profiles (Figures 1B and 2B, 3, and 5B, respectively).

An essential insight into the nature of the log-linear region is provided by LRE analysis, in that cycle efficiency can be calculated from the reaction fluorescence using equation 3. This predicts that a 15–30% loss in amplification efficiency occurs across the six cycles included in this analysis, despite the high level of linearity (Figure 6). In addition to challenging the contention that the log-linear region is representative of an exponential region, this predicted reduction in amplification efficiency would be expected to reduce the slope as compared with an invariant amplification efficiency, which in turn would reduce the resulting E_{loglin} values. Based on comparison with the "gold standard" E_{slope} values, this indeed appears to be the case. Ultimately, however, conflicting interpretations raise the question as to what is the true origin of the log-linear region.

### Resolving the origin of the log-linear region

The sigmoidal modeling presented in Figure 4 provides the opportunity to examine the characteristics of the log-linear region under a purely sigmoidal paradigm. Importantly, this also allows the dynamics of PCR amplification to be examined within the lowest regions of an amplification profile, where low reaction fluorescence intensity precludes effective E_{C} determination using equation 4. Figure 7A presents a classic Log(F_{C}) vs. cycle number plot of a sigmoidal F_{C} dataset taken from Figure 4.

This confirms the presence of a log-linear region that generates an E_{loglin} estimate of 83.5%, similar to the 79.6% E_{loglin} estimate generated from analysis of the actual F_{C} readings (Figure 6). LRE modeling predicts a large loss in amplification efficiency within this log-linear region in close agreement with that predicted from the actual F_{C} dataset (Figure 6). Another key attribute of this sigmoid-derived log-linear region is illustrated in Figure 7B, which is that the loss in amplification efficiency is exponential in nature, suggesting that it is an exponential loss in amplification efficiency that underpins the linear character of the log-linear region.

Additional perspectives on the kinetics of real-time PCR are provided by Figure 7C, which demonstrates that loss in amplification efficiency can be estimated by expressing reaction fluorescence as a percentage of F_{max}. For example, a reaction fluorescence that is 1% of F_{max} corresponds to an approximate 1% loss in amplification efficiency. Importantly, this provides a simple method for illustrating at what point loss in amplification efficiency becomes nontrivial. In relation to SYBR Green I detection under the optical capabilities of the instrumentation used in this study, 1% of F_{max} corresponds to about 100 fluorescence units, which is below the lowest fluorescence intensity that can be measured with acceptable accuracy. Accordingly, at 1,000 fluorescence units, which roughly corresponds to the middle of the log-linear region, the loss in amplification efficiency is estimated to be 10%. It would thus be difficult to obtain accurate fluorescence readings in which loss in amplification efficiency can be considered to be trivial, without a substantive increase in detection sensitivity. Of greater significance, however, is that a measurable loss in amplification efficiency occurs within 99% of a real-time amplification profile, irrespective of detection sensitivity, refuting the exponential character that has historically been ascribed to real-time PCR.

### Comparison of automated data processing packages

Although detailed analysis of a single standard curve may be sufficient to demonstrate lack of exponential character, it only provides a limited perspective as to the potential impact on currently employed data processing packages. The evaluation was therefore expanded to include four identical standard curves in which the reaction mixes were supplemented with increasing quantities of SYBR Green I. Initially based on concerns that SYBR Green I quantity could impact the quantitative accuracy of LRE modeling [27], this approach also provides some perspective as to how each method responds to the inhibition of amplification, produced by increasing SYBR Green I quantity. It should also be noted that DyNAmo (formulated with an engineered *T. brockianus* DNA polymerase fused to a non-specific DNA-binding region) was used for this analysis, demonstrating that the trends described here are not unique to enzyme formulations containing *T. aquaticus* DNA polymerase.

The evaluation was also extended to include three publicly available software packages that provide automated amplification efficiency determination, two of which are based on analysis of the log-linear region. A major challenge for both of these exponential-based packages is selecting which cycles to include in the analysis. The first, called "LinReg", uses an iterative method that searches for the greatest linear correlation coefficient based on a minimum four cycle window, in which amplification efficiency is determined from the slope [11], as is used in this study for E_{loglin} determination. The second, called "Miner", uses a complex, multi-step approach in which cycle selection is based on a series of calculations using nonlinear regression to a four-parameter logistic model to define the upper boundary of the analysis. This is followed by fitting the selected fluorescence readings to an exponential equation using an iterative nonlinear regression analysis approach, from which amplification efficiency is determined using a weighted average [15]. The third is a prototypic Java program that automates LRE analysis, which is provided as supplementary material by Rutledge and Stewart 2008 [27]. The results of the evaluation are summarized in Table 1.

Similar to that observed previously, positional analysis produced a high level of correlation between C_{t}- and C_{1/2}-based E_{slope} values. Both present a progressive loss in E_{max} as SYBR Green I quantity was increased, resulting in an 8% reduction at the highest quantity examined. Increasing SYBR Green I quantity had no apparent impact on the r^{2} of each standard curve, indicating that increasing reaction fluorescence did not improve the precision of the analysis.

As predicted by sigmoidal modeling (Figure 7) all three methods that use exponential analysis produced lower efficiency values, although the magnitude of this difference varied. Potentially more significant is that exponential analysis generated no discernible trend when SYBR Green I quantity was progressively increased, with the possible exception of Miner, which predicted a 6% loss at 2.0X SYBR Green I. It is important to also note that exponential analysis is highly influenced by the number and position of the cycles included in the analysis. This is particularly evident for LinReg in which cycle selection is based on searching for a maximal r^{2} using a variable cycle window size and position. This can produce large variances and frequently requires manual adjustment. However, regardless of implementation, the invalidity of applying exponential mathematics to the log-linear region brings into question both the reliability and accuracy generated by LinReg and Miner.

Consistent with that seen in Figure 3, automated LRE analysis produced E_{max} values similar to the E_{slope} values, except at the highest SYBR Green I quantity, which produced a ~4% lower efficiency estimate. Notwithstanding this discrepancy, this dataset does provide substantive supporting evidence that LRE analysis generates amplification efficiency estimates that correlate more closely to E_{slope}, and generate less variance than methods that rely on exponential analysis of the log-linear region.

## Conclusion

Kinetic analysis of PCR amplification based on the LRE model reveals fundamental flaws in the current interpretation of amplification efficiency dynamics, demonstrating that SYBR Green I-derived amplification profiles lack the exponential character that has historically been ascribed to real-time PCR. Consequently, methods for amplification efficiency determination that rely on exponential analysis of the log-linear region generate systematic underestimations, differing only in the extent of the bias and variability in the resulting amplification efficiency estimates. Although this apparent lack of exponential character could also be expected to compromise the efficacy of the threshold method, empirical testing demonstrates that positional analysis does not rely on the exponential character of a real-time amplification profile. Furthermore, it was demonstrated that the slope-derived efficiency estimate produced by C_{t}-based standard curves reflect the amplification efficiency not at C_{t}, but rather the maximal amplification efficiency as defined under a sigmoidal paradigm.

This study further corroborates the efficacy of LRE analysis for amplification efficiency determination, as well as providing additional insights into the linear coupling between amplification efficiency and amplicon DNA quantity. LRE analysis thus not only provides a gateway to sigmoid-based quantification, but also provides a simple methodology for analyzing amplification kinetics within individual amplification reactions. Based on analysis of the high-quality fluorescence readings within the central region of an amplification profile, LRE analysis avoids errors associated with both low reaction fluorescence and distortions associated with the plateau phase.

## Methods

PCR amplification was conducted as previously described [27] in which 3–4 replicate reactions were run for each quantity of lambda gDNA, and the F_{C} datasets averaged to generate a single amplification profile for analysis. Briefly, replicate amplification sets consisting of 5.0 μl reactions containing lambda gDNA (New England BioLabs) at the specified quantity and 500 μM of the lambda primers K7B (CTGCTGGCCGGAACTAATGAATTTATTGGT) and K12 (ATGCCACGATGCCTCATCACTGTTG). The standard curve presented in Figure 1 employed QuantiTect (Qiagen) enzyme formulation, whereas DyNAmo (Finnzymes, distributed by New England BioLabs) was used for the standard curves containing increasing quantities of SYBR Green I (Table 1). SYBR Green I was diluted to the appropriate quantity using ddH_{2}0 before addition to the PCR master mix just prior to amplification reaction preparation and is expressed in units designated by the manufacturer (Invitrogen).

All amplifications were conducted with a Mx3000P spectrofluorometric thermal cycler (Stratagene) using a two temperature cycling regime initiated with a 15 min activation at 95°C, followed by 50 cycles of 120 s annealing and elongation at 65°C and a 10 s denaturation at 95°C. To increase optical precision, three fluorescent reads were taken at the end of the annealing and elongation step and the average used as an estimate of reaction fluorescence. Specificity of amplification was confirmed by melting curve analysis conducted at the end of each run.

An extensive description of the development and implementation of the LRE method is provided by Rutledge and Stewart [27]. Automated LRE analysis was conducted using the prototypic Java program provided as supplementary materials in this earlier study using default values.

## References

Pfaffl MW:

**A new mathematical model for relative quantification in real-time RT-PCR.***Nucleic Acids Res*2001,**29**(9):e45. 10.1093/nar/29.9.e45.Rutledge RG, Côté C:

**Mathematics of quantitative kinetic PCR and the application of standard curves.***Nucleic Acids Res*2003,**31**(16):e93. 10.1093/nar/gng093.Meijerink J, Mandigers C, Locht L, Tönnissen E, Goodsaid F, Raemaekers J:

**A novel method to compensate for different amplification efficiencies between patient DNA samples in quantitative real-time PCR.***J Mol Diagn*2001,**3**(2):55-61.Ståhlberg A, Åman P, Ridell B, Mostad P, Kubista M:

**Quantitative real-time PCR method for detection of B-lymphocyte monoclonality by comparison of kappa and lambda immunoglobulin light chain expression.***Clin Chem*2003,**49**(1):51-59. 10.1373/49.1.51.Chelly J, Kaplan JC, Maire P, Gautron S, Kahn A:

**Transcription of the dystrophin gene in human muscle and non-muscle tissue.***Nature*1988,**333**(6176):858-860. 10.1038/333858a0.Wiesner RJ, Rüegg JC, Morano I:

**Counting target molecules by exponential polymerase chain reaction: copy number of mitochondrial DNA in rat tissues.***Biochem Biophys Res Commun*1992,**183**(2):553-559. 10.1016/0006-291X(92)90517-O.Wiesner RJ:

**Direct quantification of picomolar concentrations of mRNAs by mathematical analysis of a reverse transcription/exponential polymerase chain reaction assay.***Nucleic Acids Res*1992,**20**(21):5863-5864. 10.1093/nar/20.21.5863.Harting I, Wiesner RJ:

**Quantification of transcript-to-template ratios as a measure of gene expression using RT-PCR.***BioTechniques*1997,**23**(3):450-455.Gentle A, Anastasopoulos F, McBrien NA:

**High-resolution semi-quantitative real-time PCR without the use of a standard curve.***Biotechniques*2001,**31**(3):502.Liu W, Saint DA:

**A new quantitative method of real time reverse transcription polymerase chain reaction assay based on simulation of polymerase chain reaction kinetics.***Anal Biochem*2002,**302**(1):52-59. 10.1006/abio.2001.5530.Ramakers C, Ruijter JM, Deprez RH, Moorman AF:

**Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data.***Neurosci Lett*2003,**339**(1):62-66. 10.1016/S0304-3940(02)01423-4.Peirson SN, Butler JN, Foster RG:

**Experimental validation of novel and conventional approaches to quantitative real-time PCR data analysis.***Nucleic Acids Res*2003,**31**(14):e73. 10.1093/nar/gng073.Bar T, Stahlberg A, Muszta A, Kubista M:

**Kinetic Outlier Detection (KOD) in real-time PCR.***Nucleic Acids Res*2003,**31**(17):e105. 10.1093/nar/gng106.Wilhelm J, Pingoud A, Hahn M:

**SoFAR: software for fully automatic evaluation of real-time PCR data.***BioTechniques*2003,**34**(2):324-332.Zhao S, Fernald RD:

**Comprehensive algorithm for quantitative real-time polymerase chain reaction.***J Comput Biol*2005,**12**(8):1047-1064. 10.1089/cmb.2005.12.1047.Liu W, Saint DA:

**Validation of a quantitative method for real time PCR kinetics.***Biochem Biophys Res Commun*2002,**294**(2):347-353. 10.1016/S0006-291X(02)00478-3.Tichopad A, Dzidic A, Pfaffl MW:

**Improving quantitative real-time RT-PCR reproducibility by boosting primer-linked amplification efficiency.***Biotechnol Lett*2002,**24:**2053-2056. 10.1023/A:1021319421153.Rutledge RG:

**Sigmoidal curve-fitting redefines quantitative real-time PCR with the prospective of developing automated high-throughput applications.***Nucleic Acids Res*2004,**32**(22):e178. 10.1093/nar/gnh177.Goll R, Olsen T, Cui G, Florholmen J:

**Evaluation of absolute quantitation by nonlinear regression in probe-based real-time PCR.***BMC Bioinformatics*2006,**7:**107. 10.1186/1471-2105-7-107.Karlen Y, McNair A, Perseguers S, Mazza C, Mermod N:

**Statistical significance of quantitative PCR.***BMC Bioinformatics*2007,**8:**131. 10.1186/1471-2105-8-131.Roussel Y, Harris A, Lee MH, Wilks M:

**Novel methods of quantitative real-time PCR data analysis in a murine Helicobacter pylori vaccine model.***Vaccine*2007,**25**(15):2919-2929. 10.1016/j.vaccine.2006.07.013.Qiu H, Durand K, Rabinovitch-Chable H, Rigaud M, Gazaille V, Clavere P, Sturtz FG:

**Gene expression of HIF-1alpha and XRCC4 measured in human samples by real-time RT-PCR using the sigmoidal curve-fitting method.***Biotechniques*2007,**42**(3):355-362. 10.2144/000112331.Cikos S, Bukovska A, Koppel J:

**Relative quantification of mRNA: comparison of methods currently used for real-time PCR data analysis.***BMC Mol Biol*2007,**8:**113. 10.1186/1471-2199-8-113.Durtschi JD, Stevenson J, Hymas W, Voelkerding KV:

**Evaluation of quantification methods for real-time PCR minor groove binding hybridization probe assays.***Anal Biochem*2007,**361**(1):55-64. 10.1016/j.ab.2006.11.023.Swillens S, Dessars B, Housni HE:

**Revisiting the sigmoidal curve fitting applied to quantitative real-time PCR data.***Anal Biochem*2008,**373**(2):370-376. 10.1016/j.ab.2007.10.019.Feng J, Zeng R, Chen J:

**Accurate and efficient data processing for quantitative real-time PCR using a tripartite plant virus as a model.***Biotechniques*2008,**44**(7):901-912. 10.2144/000112750.Rutledge RG, Stewart D:

**A kinetic-based sigmoidal model for the polymerase chain reaction and its application to high-capacity absolute quantitative real-time PCR.***BMC Biotechnol*2008,**8:**47. 10.1186/1472-6750-8-47.Tichopad A, Dilger M, Schwarz G, Pfaffl MW:

**Standardized determination of real-time PCR efficiency from a single reaction set-up.***Nucleic Acids Res*2003,**31**(20):e122. 10.1093/nar/gng122.Schefe JH, Lehmann KE, Buschmann IR, Unger T, Funke-Kaiser H:

**Quantitative real-time RT-PCR data analysis: current concepts and the novel "gene expression's C**_{ T }**difference" formula.***J Mol Med*2006,**84**(11):901-910. 10.1007/s00109-006-0097-6 10.1007/s00109-006-0097-6Chervoneva I, Li Y, Iglewicz B, Waldman S, Hyslop T:

**Relative quantification based on logistic models for individual polymerase chain reactions.***Stat Med*2007,**26**(30):5596-5611. 10.1002/sim.3127.Czechowski T, Bari RP, Stitt M, Scheible WR, Udvardi MK:

**Real-time RT-PCR profiling of over 1400 Arabidopsis transcription factors: unprecedented sensitivity reveals novel root- and shoot-specific genes.***Plant J*2004,**38**(2):366-379. 10.1111/j.1365-313X.2004.02051.x.Caldana C, Scheible WR, Mueller-Roeber B, Ruzicic S:

**A quantitative RT-PCR platform for high-throughput expression profiling of 2500 rice transcription factors.***Plant Methods*2007,**3:**7. 10.1186/1746-4811-3-7.

## Acknowledgements

We are indebted to Frédéric Vigneault and Brian Boyle for critically reviewing the manuscript, and thank Pamela Cheers for editorial assistance. The research was supported in part by a grant to RGR from the Canadian Biotechnology Strategy.

## Author information

### Authors and Affiliations

### Corresponding authors

## Additional information

### Authors' contributions

RGR conceived and led the project, in addition to drafting the manuscript. DS conducted all of the real-time qPCR and contributed to the data analysis.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Rutledge, R.G., Stewart, D. Critical evaluation of methods used to determine amplification efficiency refutes the exponential character of real-time PCR.
*BMC Molecular Biol* **9**, 96 (2008). https://doi.org/10.1186/1471-2199-9-96

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1471-2199-9-96

### Keywords

- Amplification Efficiency
- Threshold Method
- Fluorescence Reading
- Amplification Profile
- Positional Analysis