 Methodology Article
 Open access
 Published:
A fixedpoint algorithm for estimating amplification efficiency from a polymerase chain reaction dilution series
BMC Bioinformatics volumeÂ 15, ArticleÂ number:Â 372 (2014)
Abstract
Background
The polymerase chain reaction amplifies and quantifies small amounts of DNA. It is a cyclic process, during each cycle of which each strand of template DNA is copied with probability approaching one: the amount of DNA approximately doubles and this amount can be estimated fluorimetrically each cycle, producing a set of fluorescence values hereafter referred to as the amplification curve. Commonly the biological question of relevance is one of the ratio of DNA concentrations in two samples: a ratio that is deduced by comparing the two amplification curves, usually by way of a plot of fluorescence against cycle number. Central to this analysis is measuring the extent to which one amplification curve is shifted relative to the other, a measurement often accomplished by defining a threshold or quantification cycle, C _{ q }, for each curve: the fractional cycle number at which fluorescence reaches some threshold or at which some other criterion (maximum slope, maximum rate of change of slope) is satisfied.
We propose an alternative where position is measured relative to a reference curve; position equates to the cycle shift which maximizes the correlation between the reference and the observed fluorescence sequence. A key parameter of the reference curve is obtained by fixedpoint convergence.
Results
We consider the analysis of dilution series constructed for the estimation of qPCR amplification efficiency. The estimate of amplification efficiency is based on the slope of the regression line when the C _{ q } is plotted against the logarithm of dilution. We compare the approach to three commonly used methods for determining C _{ q }; each is applied to publicly accessible calibration data sets, and to ten from our own laboratory. As in the established literature we judge their relative merits both from the standard deviation of the slope of the calibration curve, and from the variance in C _{ q } for replicate fluorescence curves.
Conclusions
The approach does not require modification of experimental protocols, and can be applied retrospectively to existing data. We recommend that it be added to the methodological toolkit with which laboratories interpret their realtime PCR data.
Background
Since its introduction by Mullis et al. [1], the polymerase chain reaction (PCR) has been widely used to amplify and quantify small amounts of DNA. Briefly it constitutes a few dozen cycles in each of which there are three stages: denaturation, annealing and extension.
At each cycle each DNA strand either doubles (with probability p) or fails to double. Defining amplification efficiency, E, as E=1+p the expected number N _{ c } of strands after C cycles, given constant p, is
where N _{0} is the initial number. We use this definition of E to be consistent with Ruitjer et al. [2] and in that we deviate from MIQE guidelines [3].
The amount of DNA can be estimated fluorimetrically using, for instance, Sybr Green 1 dye which binds to the minor groove in doublestranded DNA [4]. Commonly the relevant question biologically is one of relative quantification: the ratio of initial DNA in two samples. If samples A and B, initially containing N _{ A } and N _{ B } strands, exhibit the same fluorescence after C _{ a } and C _{ b } cycles respectively, then assuming constant p,
from which we have
where âˆ† C=C _{ b }C _{ a }. Accordingly, the cycle difference, âˆ† C, and E (or equivalently, p), are the key to estimation of the ratio N _{ A }/N _{ B }.
In common use there are two fundamentally different approaches to the estimation of E. One involves the generation of curves from a series of dilutions: an eightfold dilution, for instance, would delay the fluorescence by three cycles if E=2 because 2^{3}=8. If E=1.8, however, it would delay the fluorescence curve by 3.53 cycles, because 1.8^{3.53}=8. Alternatively (Gentle et al., 2001. [5], Zhao and Fernald, 2005 [6]), E can be estimated from the realtime fluorescence data. At each cycle early in the sequence, fluorescence above background will increase by a factor E, which can therefore be estimated from the data in this `loglinearâ€™ phase. The early phase of exponential increase is shortlived. As resources become limiting, the fluorescence curve flattens out, and Liu and Saint [7] have used the sigmoid or logistic function,
where F _{ b } is background, and F _{ max } is maximum contribution of the reaction to fluorescence, (the asymptote, rather than maximum observed experimentally), to describe the data. Rutledge and Stewart [8] introduced an analysis which takes into account the linear decrease in amplification under this model, simplifying the estimation of the initial amplification efficiency from the curve itself. MIQE guidelines [3] recommend the former approach: `PCR amplification efficiency must be established by means of calibration curves...â€™ but we acknowledge ongoing debate on this issue.
Strictly speaking the data from a tube are discontinuous; fluorescence is measured at the end of each cycle, and there is no such thing as a fluorescence after a fractional number of cycles as implied by the continuous functions above. We use the term reference curve to imply an abstraction; a smooth continuous curve of fluorescence as a function of x, which we observe at cycles C which are integer values of x. The observed fluorescence is the fluorescence at these integer values, but with the addition of error or noise.
A key to analyzing PCR, therefore, is, given two fluorescence curves, to measure âˆ† C, the extent to which one curve is shifted laterally relative to the other. There are two very different circumstances under which one may need to do this. If E is to be estimated from the cycletocycle increase in fluorescence of a single assay tube, then quantifying some aspect of the fluorescence is important. Conversely, if one is using dilution to estimate amplification then the shape of the curve is of less import so long as the data and the reference curve have in common that they are Sshaped (sigmoid): interest lies only in the extent to which dilution has shifted the curve, of whatever shape, to the right. Whatever the method of estimating E, that estimate is commonly used subsequently to derive, in concert with a measured cycle difference between two tubes, âˆ† C, the initial concentration ratio implied by Eq. 1.
It is the estimation of cycle shift in these scenarios which we address; to what extent is one fluorescence curve shifted relative to another? There is, of course, a significant literature detailing several algorithms to do just that, and we should justify any attempt to add another. Ruitjer et al. [2] have examined the performance of nine estimators of E and have proposed several measures of their relative merits. In using the publicly available data sets comprising dilution series for establishing amplification efficiency, two measures are of central importance. One is the withinreplicate variance; most data sets have three or more replicates at each dilution, and for a good estimator we expect values of C _{ q } from these replicates to be close. The second measure is the standard deviation of the estimate of the slope when C _{ q } is regressed against the logarithm of dilution; the smaller the standard deviation of the slope, the smaller will be that of the estimated efficiency. Following Ruitjer et al., we use both of these, and compare approaches using Friedman's nonparametric rank sum.
The three algorithms which we examine in detail, and which performed very well in the review by Ruitjer et al. are C y0, Standard C _{ q }, and PCRMiner. The latter algorithm includes both an estimate of C _{ q }, and an estimate of efficiency derived from each curve, and we should emphasize that we are implementing only the C _{ q }estimating component of PCRMiner. To avoid confusion with the full PCRMiner algorithm we will refer to it as the SDMl5 method (second derivative maximum of the model designated l5 in the qpcR package [9] associated with the R statistical software [10]).
Notwithstanding their established utility we have concerns about each of these approaches. For C y0, the derived C _{ q } depends on the baseline. We regard baseline fluorescence as a `nuisanceâ€™ parameter, as do several algorithms that attempt to eliminate it. Our bias (and we accept that it is personal bias) is to use an estimator independent of baseline. Standard C _{ q } finds the fluorescence (F _{ q }) at the second derivative maximum (SDM) for the (mean) undiluted sample, and for subsequent samples C _{ q } is defined as the (interpolated) cycle at which fluorescence achieves F _{ q }. Again this is influenced by sampletosample variation in baseline, and for the subsequent diluted samples it takes information from only two readings out of the entire curve. The SDMl4 approach overcomes the above reservations by fitting a fourparameter sigmoid curve and calculating the cycle of SDM as implemented, for instance in some commercially available software [11]. This approach is independent both of baseline and of scale, but it raises a more subtle problem. Each reference curve is of a different shape. The distance between the curves in a dilution series is not well defined if each curve is of a different shape; they are not the same curve translated laterally along the dilution axis. The distance between the second derivative maxima, for instance, is different from the distance between the first derivative maxima.
Methods
Biological methods
Ethics approval for the use of peripheral blood leucocytes was obtained from the Human Researach Ethics Committee of The Queen Elizabeth Hospital (South Australia), and the use of samples followed the protocol approved by that committee, as documented in BiancoMiotto et al. [12].
RNA extraction and reverse transcription
RNA was extracted from cells grown in tissue culture using Trizol (Invitrogen, USA) according to the manufacturer's protocol. The concentration of RNA was determined using a Biophotometer (Eppendorf, North America Inc, Westbury, USA). DNAse treatment of total RNA was performed prior to reverse transcription in order to minimize PCR signal arising from carryover genomic DNA (Ambion DNAfree kit). RNA was reverse transcribed using Superscript III RT (Invitrogen, USA). cDNA was diluted 20 fold in ultra pure water (Fischer Biotech) prior to real time PCR.
Preparation of genomic DNA
Mononuclear cells were isolated from the peripheral blood of healthy donors using Lymphoprep (AxisShield, Oslo, Norway) according to the manufacturer's instructions. Genomic DNA (gDNA) was purified from the mononuclear cells using Trizol (Invitrogen Life Technologies, NY, USA) according to the manufacturer's instructions.
Preparation of dilution series
50 Î¼l of ultra pure water was aliquoted into a series of 0.5 ml PCR tubes, and either 50 Î¼l of gDNA or 50 Î¼l of cDNA was added to the first tube and mixed by pipetting up and down 10 times. 50 Î¼l of this mixture was then pipetted to the next tube and mixed, and the process repeated across the tubes, to produce a twofold serial dilution.
Realtime polymerase chain reaction
PCR amplification was performed in 20 Î¼L final volumes containing 6 Î¼L of cDNA or gDNA template, 2 Î¼L of each forward and reverse primer (5 Î¼M), and 10 Î¼L of 2 Ã— Quantitect Sybr Green Master Mix (Qiagen, Germany). Thermocycling was performed in a Rotorgene 6000 thermocycler (Corbett, Australia) with an initial activation/denaturation (hot start) at 95Â°C for 15 min; followed by 45 cycles of 20 sec at 95Â°C, 30 sec at the annealing temperature, and 30 sec extension at 72Â°C. After the cycling there was a final extension at 72Â°C for 4 min. Triple replicates of twelve (sometimes eleven) 2fold dilutions reactions were performed on all samples. Products were then melted in the Rotorgene 6000 thermocycler from 60Â°C to 99Â°C at 0.5Â°C for 5 sec per step to determine if the PCR products melted at the same temperature as PCR products that had been fractionated through 1% agarose gel to confirm that the product was of the predicted size.
Details of amplicons and primers appear under Additional file 1: Table S1.
Numerical methods
Data analysis was carried out under GNU/Linux Ubuntu 14.04 LTS using the R programming language [10] and the associated packages qpcR [9] and ggplot2 [13]. The fixedpoint estimator is as documented below. The methods Standard C _{ q }, SDMl4 nd C y0 were implemented as follows.
Standard C _{ q }
The essence of standard C _{ q } is to locate the fractional cycle corresponding to the SDM of the (averaged) undiluted sample, and to define F _{ q } as the (interpolated) fluorescence at that fractional cycle. The C _{ q } of each cycle, diluted or undiluted, is the fractional cycle at which F _{ q } is achieved.
If F _{ i } denotes the fluorescence at the i ^{th} cycle of the averaged undiluted samples we find i for which the second derivative, (F _{i1}2F _{ i }+F _{i+1}) is maximal, and then assuming that the second derivative of fluorescence, if continuous, would be adequately approximated by a quadratic around the i ^{th} cycle we now have as fractional cycle maximizing that quadratic as the location of SDM. The (mean) fluorescence of the undiluted sample is then found by interpolating the cubic through the adjacent four fluorescence values. This defines F _{ q }. For each sample we then find k such that F _{k1}<F _{ q }<F _{ k } implying that C _{ q } for that sample lies between (k1) and k. Again the fractional cycle at which F _{ q } occurs is found by cubic interpolation of the observed fluorescence at F _{k2}â€¦F _{k+1}.
SDML4
The function pcrfit() from the qpcR statistical package finds for each dilution curve, the parameters of best fit for the fourparameter model defined (using the nomenclature from Zhao and Ferdinand [6])
from which the location of the SDM is given by
Cy0
The function pcrfit() as above finds the parameters of best fit for a fiveparameter sigmoidal curve. As introduced by Guiscini [14], the function used was
although the example in the qpcR package uses the closely related
and we have implemented both for comparison. We denote the former Cy0b5, the latter Cy0l5, referring to the fiveparameter functions b5 and l5 of the qpcR package.
The function Cy0 from the qpcR package takes the five parameter function and returns C y0 as the point of intersection with the abscissa of the tangent through the maximum first derivative.
Theoretical development of fixedpoint approach
In estimating âˆ† C we are quantifying the extent to which one curve needs to be shifted horizontally (on the cycle axis) in order that it might overlie the other. That aim requires three qualifications: first, that there may need to be some vertical shift to accommodate different baselines; second, that the same applies to scale; third that we have equallyspaced points rather than a continuous curve.
If, as in some standard analyses, the `positionâ€™ of a fluorescence curve is taken to be the fractional cycle at which fluorescence attains some arbitrary threshold, then the tubetotube variation in the baseline and scale of fluorescence becomes a problem; scale is particularly so where fluorescence has not reached a terminal plateau. The appeal of using position of maximum first or second derivative (as in PCRMiner software) is that these are not influenced by changes in baseline or scale.
We can ask how much one fluorescence curve needs to be shifted such that it overlies another, but because we have points, rather than continuous curves we will usually find that, at best, one set of points lies close to, but between, the points of the other.
A reference curve
The strategy commonly adopted, and which we adopt here, is to fit our observed fluorescence data to some continuous function. The functions discussed below have in common that they are Sshaped (sigmoid). The common definition of a good fit is one which minimizes the sum of squares of differences between the observed data and the continuous function being fitted, and again we follow that practice. The qpcR package [9] uses the MarquardtLevenberg algorithm to accomplish this. Commonly used are the fiveparameter functions adopted in the variants of the C y0 estimator introduced earlier simplified version of which are obtained with the constraint f=1, giving the fourparameter curve used in SDMl4 or the fourparameter curve of Eq. 2 as in Liu and Saint [7]. In explaining the fixedpoint approach is convenient to note that we can recast Eq. 2 as
where A=e ^{Î²}. We still have four parameters, but instead of varying Î² to obtain best fit we are varying A=e ^{Î²}. This makes no difference to the fit; it just makes the physical significance of the parameters more obvious to the reader.
Looking at the four parameters in turn we have
This is the background fluorescence of an assay which we are assuming to be a nuisance variable. We want our estimate, C _{ q } to be independent of Fb.
This is the difference in fluorescence between Fb and the asymptote which fluorescence is approaching. We take tubetotube variation in F _{ max } to result from differences in such factors as the opacity of the assay tubes; it is a nuisance variable, and C _{ q } should be invariant with F _{ max }.
This determines shift along the abscissa (cycle axis). In the fourparameter models mentioned, it is the fractional cycle at which the fluorescence representing reaction product is 50% of F _{ max }.
This determines the shape of the curve and we note that the increase in fluorescence due to polymerase reaction during the first cycle is
which, for large {A}^{{x}_{0}} tends to A. That is, the parameter, A is the amplification efficiency during the early, exponential, part of the chain reaction.
If we now want to fit the fluorescence in each tube of a dilution series to the `sameâ€™ continuous function, it seems preferable to use the function which models the appropriate amplification efficiency. By the `sameâ€™ continuous function we mean functions for which A is the same. We accept that Fb, F _{ max } and x _{0} will vary from tube to tube, because baseline, scale and C _{ q } will vary from tube to tube. But the idea of running a dilution series to determine amplification efficiency is predicated on the assumption that amplification efficiency, and hence A, is the same for every tube. We might seek, therefore, to find for each tube the values for Fb, F _{ max } and x _{0} that best fit the observed fluorescence in that tube, keeping A fixed at the amplification efficiency as derived in the usual way from a regression of C _{ q } against logarithm of dilution.
The impasse is now obvious; the point of dilution assay is to determine the amplification efficiency. Until we know the amplification efficiency we do not know the appropriate value of A to use in Eq. 5 to determine the C _{ q } for each tube. The Bauer fixedpoint theorem resolves that impasse. Using fixedpoint convergence (see, for instance, [15]), we begin with an initial guess, A _{0}. This leads to estimates C _{ q } for each tube on the basis of which the slope of regression against logarithm of dilution gives a first estimate E _{1} of amplification efficiency. We now replace our initial guesstimate A _{0} with A _{1}=E _{1} and repeat the process giving a second estimate E _{2} and so on, until subsequent estimates are unchanged and convergence has been achieved.
For the process to converge, the requirement of the fixedpoint theorem is that a plot of E _{1} against A _{0} (which is, of course also that of E _{2} against A _{1}and so on), should have an absolute slope less than one. Providing this condition is satisfied (and for the data sets considered here it is), the theorem asserts that
the process will converge
the smaller the slope, the faster it will converge
the value to which it converges is independent of the starting estimate A _{0}
We note that, in fitting Eq. 5 to an observed fluorescence curve, if A is fixed, there are only three parameters, Fb, F _{ max } and x _{0} of which both Fb and F _{ max } enter the equation linearly. Since we are really only interested in the parameter x _{0} we can, under these circumstances, avoid using the MarquardtLevenberg algorithm and find the value of x _{0} which maximizes the Pearson correlation coefficient between observed fluorescence and the fitted curve. That correlation will be independent of Fb and F _{ max }, and will give the same `best fitâ€™ x _{0} as a MarquardtLevenberg approach minimizing the sum of squares of differences. If we denote by x _{ i } the value of x _{0} maximizing correlation in the i ^{th} tube, each tube now has a reference curve of the form
The reference curves for all tubes are now the same shape, apart from their shift along the abscissa determined by x _{ i }. Because shift along the abscissa is the only difference between any two curves the concept of âˆ† C _{ q } as that difference is now unambiguous; we could define it as the difference between first derivative maximum, second derivative maximum, or difference between cycles at which some fraction of the increase has been achieved and all these definitions would result in the same âˆ† C _{ q }. The simplest definition is to use x _{ i } as the C _{ q } for the corresponding tube (the cycle at which the reaction is 50% of its maximum fluorescence).
In summary the steps in the fixedpoint algorithm are as follows. We will use the publicly available dataset batsch1 from the qpcR package as an illustration, including the numerical values which we obtain.

1.
Using an initial guesstimate of amplification efficiency (we have used A=1.5) define
f(A,{x}_{0},x)=f(1.5,{x}_{0},x)=\frac{1}{1+1.{5}^{({x}_{0}x)}}. 
2.
Use algorithm of choice to find, for each tube, the value of x _{0} which maximizes the correlation between f(A,x _{0},x) and the fluorescence data from that tube. For that tube let C _{ q }=x _{0}. The first fluorescence data from batsch1 is shown in Figure 1. By inspection we can see that half the generated fluorescence occurs by about cycle 29, so the reference curve showing greatest correlation will have a value of x _{0} of about 29. Figure 2 shows a plot of correlation against x _{0} for values of x _{0} from 1 to 45 and we can see the correlation maximizing at cycle 29. A good firstestimate of the fractional cycle at which correlation maximizes is obtained by quadratic interpolation using the three correlations at this and the adjacent cycles. Given the correlations are 0.99497, 0.99766, and 0.99408 at cycles 28, 29 and 30 respectively we imply a maximum at cycle
{x}_{0}=\frac{0.994970.99408}{2(0.994972\xc3\u20140.99766+0.99408)}=28.929.
Of course, any iterative approach can be used, and as we are only looking for a maximum in one dimension the R procedure optimize() is appropriate. In practice, for these data there is no improvement; a maximum correlation is found at cycle fractional cycle 28.930. Set C _{ q } for this tube as 28.93. Repeat for each tube.

3.
Regressing C _{ q } against logarithm of dilution determine estimated amplification efficiency, and return to step 1 replacing A with this estimate. We prefer to use logarithms to base 2 as in Figure 3 because the implication of the regression slope is clear from inspection; a doubling at each cycle (E=2) would imply that a twofold dilution shifts the fluorescence curve by exactly one cycle. The regression slope in Figure 3 is 1.185, implying that it takes 1.185 cycles to compensate for a twofold dilution. If there is an Efold increase each cycle then
{E}^{1.185}=2
from which it is immediate that

4.
Return to step 1, replacing the initial A=1.5 with A=1.795. Iterate until estimated efficiency is unchanged. This is the fixedpoint iteration, and is illustrated in Figure 4. In the above we started with guesstimate A _{0}=1.5 (deliberately far from what we expect, so as to illustrate the method) and the first cycle returns an estimated efficiency 1.795. This corresponds to the vertical blue line on Figure 4 in which the red curve shows for this process the output E for a range on input A from 1Â·4 to 2Â·2. Bauer's fixed point theorem guarantees convergence if the slope of this line is absolutely less than one in the region of interest (as shown here for these data). We then replace our initial A _{0}=1.5 with the revised A _{1}=1.795 (horizontal blue line, completing the first iteration of fixedpoint convergence. The second cycle is shown in green, giving the second estimate E _{2}=1.848 and so on, converging rapidly to an amplification efficiency of 1.8517.
Results and discussion
Fixed point convergence
For all data sets with an initial estimate A _{0}=1.5 used in the logistic function, deliberately far from the expected E, the fixedpoint iteration had converged to within.001 of its asymptote in four iterations; a convergence well beyond any biological relevance. The rapid convergence results from the relative insensitivity of estimated E to the A used in the logistic reference curve. Figure 5 shows estimated E _{1} as a function of input A _{0} for all ten data sets from our laboratory. The slope of all curves in the interval of interest is well within the (1,1) required for convergence by the Banach fixedpoint theorem. For nine of the ten data sets, however, the parameter of the logistic function does influence the corresponding E estimate and it is important that the parameter be determined objectively. Additional files 2 and 3 show the corresponding figures for the publicly available data sets.
There are 23 data sets, giving 23 standard deviations for estimates of slope by each of the five analysis methods. The rank sums appear in the first line of Table 1. Using Friedman's rank sum to compare the standard deviation of the estimate of slope the null hypothesis of equal distributions is rejected (Ï‡ _{4}^{2}=25.7,p<0.0001). The best two estimators are SDMl4 and Fixedpoint. Proceeding to a direct comparison of the two, the Fixedpoint estimator performs better, having a smaller standard deviation in 15 of 23 comparison. The difference, however, is not of statistical significance.
Within the 23 data sets are 181 sets of replicate C _{ q }. The rank sums of replicate standard deviations appear as the second line of Table 1. Again using Friedman's rank sum test on the null hypothesis of equal distributions of replicate standard deviations the null hypothesis is rejected (Ï‡ _{4}^{2}=149,p<0.0001). On direct comparison of the better two estimators, SDMl4 and Fixedpoint, the Fixedpoint is the better estimator in 96 of 181 comparisons but again this is not statistically significant.
Caveats
We would caution, however, against too literal an interpretation of the standard deviation of the slope estimate. The linear regression of position against log(d i l u t i o n) on which it is based assumes normallydistributed homoscedastic data. Data from qPCR dilution series are almost never homoscedastic because higher dilutions lead to more variable fluorescence sequences than do lower dilutions. Even if errors at low dilutions are Gaussian, those at higher dilutions result from Poisson effects and will not be Gaussian. Figure 6 is characteristic of the regression line and associated data, and demonstrates the increased variance at higher dilutions. In addition, the linear regression assumes that the dilution itself is errorfree, and we doubt that this is true of our own data. Although these concerns lead us to doubt the absolute values of the standard deviations, the concerns apply equally to the five estimators, and are unlikely to draw into question the ordering of these for any one data set.
When C _{ q } is regressed against log dilution the standard deviation of the slope depends in part on whether C _{ q } increases linearly with log dilution, and this in turn depends on amplification efficiency being constant at all dilutions. In using the variance of slope as a measure of the merit of a method, we assume that amplification efficiency is invariant with dilution. In practice the data sets we have analyzed show a remarkable linearity; the reason for assaying at intermediate dilutions is to confirm that linearity, without which a dilution series would be difficult to interpret.
The fixedpoint method assumes that the fluorescence data approach the plateau. If, in a dilution series, the higher dilutions result in only the very early part of the fluorescence curve emerging, then the estimated C _{ q } at these dilutions will be unreliable.
Finally we have presented a comparison of the five methods discussed, and use Friedman's nonparametric rank sum as a test of the null hypothesis that the methods are equivalent. Our data, however, are not randomly selected from the population of dilution series in general, and the Friedman's test should be interpreted with caution in this context. We have examined two `meritsâ€™ of the methods: replicate standard deviation and slope standard deviation. These are not independent: the standard deviation of the slope estimate takes into account that of the replicates.
Conclusion
The use of a reference curve, (in this case logistic) relative to which the position of fluorescence data can be measured, avoids subjective decisions as to baseline and scale and threshold. Using data from the whole curve, rather than just a few points, it offers an approach to the estimation of amplification efficiency from a dilution series. The logistic function represents a family of curves, however, and the specific curve appropriate to a given dilution series can be defined by fixedpoint iteration. Convergence is rapid and for the illustrative data used here the method is often, but not always, an improvement on existing estimators.
Authors' contributions
DW and DH conceived and oversaw the research. DH, TW and GM constructed the primers and carried out the experimental work. MJ and GM devised, carried out, and documented the data analysis. All authors read and approved the final manuscript.
Additional files
References
Mullis K, Faloona F, Scharf S, Saiki R, Horn G, Erlich H: Specific enzymatic amplification of dna in vitro; the polymerase chain reaction . Cold Spring Harbor Symp. Quant. Biol . 1986, Cold Spring Harbor Laboratory Press, New York, 263273.
Ruitjer J, Pfaffl M, Zhao S, Speiss A, Boggy G, Blom J, Rutledge R, Sisti D, Lievens A, De Preter K, Derveaux S, Hellemans J: JVandesompele: Evaluation of qpcr curve analysis methods for reliable biomarker discovery: Bias, resolution, precision and implications . Methods. 2013, 59: 3246. 10.1016/j.ymeth.2012.08.011.
Bustin S, Benes V, Garson J, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl M, Shipley G, Vandesompele J, Wittwer C: The miqe guidelines: Minimum information for publication of quantitative realtime pcr experiements . Clin Chem. 2009, 55: 611622. 10.1373/clinchem.2008.112797.
Morrison T, Weis J, Wittwer C: Quantification of lowcopy transcripts by continuous sybr green i monitoring during amplification . Biotechniques. 1998, 24: 954962.
Gentle A, Anastasopoulos F, McBrien N: Highresolution semiquantitative realtime pcr without the use of a standard curve . Biotechniques. 2001, 31: 502508.
Zhao S, Fernald R: Comprehensive algorithm for quantitative realtime polymerase chain reaction . J Comput Biol. 2005, 12 (8): 10471064. 10.1089/cmb.2005.12.1047.
Liu W, Saint D: Validation of a quantitative method for real time pcr kinetics . Biochem Biophys Res Comm. 2002, 249: 347353. 10.1016/S0006291X(02)004783.
Rutledge R, Stewart D: A kineticbased sigmoidal model for the polymerase chain reaction and its application to highcapacity absolute quantative realtime pcr . BMC Biotechnol. 2008, 8: 4710.1186/14726750847.
AndrejNikolai Spiess: qpcR: Modelling and Analysis of Realtime PCR Data2014. http://cran.rproject.org/web/packages/qpcR/index.html.
R Foundation for Statistical Computing: R: A Language and Environment for Statistical Computing . 2014, R Foundation for Statistical Computing, Vienna, Austria
Roche Molecular Biochemical: LightCycler Software Version 3.5 . 2001, Indiana, Roche Molecular Biochemical
BiancoMiotto T, Hussey D, Kay T, Oâ€™Keefe D, Dobrovic A: Dna methylation of the ABO promoter underlies loss of ABO allelic expresson in a significant proportion of leukaemic patients . PLoS ONE. 2009, 4 (3): e478810.1371/journal.pone.0004788.
Wickham H: Ggplot2: Elegant Graphics for Data Analysis . 2009, Springer, New York
Guiscini M, Sisti D, Rocchi M, Stocchi L, Stocchi V: A new realtime pcr method to overcome significant quantitative inaccuracy due to slight amplification inhibition . BMC Bioinformatics. 2008, 9: 32610.1186/147121059326.
Chapra S, Canale R: Numerical Methods for Engineers . 2006, McGrawHill, Boston
Acknowledgments
This work was funded by grant APP1008337 from the National Health and Medical Research Council of Australia. We gratefullly acknowledge the recommendations of two anonymous reviewers contributed substantially to the final form of this manuscript.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Electronic supplementary material
12859_2014_372_MOESM2_ESM.pdf
Additional file 2: Convergence of Fixedpoint Iteration. This figure corresponds to Figure 5 of the main text, but looks at the publiclyavailable data sets batsch1 to batsch5 (blue, brown, red, purple, black) and reps (turquoise), reps2 (green), reps3 (yellow). (PDF 74 KB)
12859_2014_372_MOESM3_ESM.pdf
Additional file 3: Convergence of Fixedpoint Iteration. As for Additional file 2, but with the publiclyavailable data sets boggy (blue), guiscini (brown), lievens1 (red), rutledge (purple) and sisti (green). (PDF 70 KB)
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Jones, M.E., Mayne, G.C., Wang, T. et al. A fixedpoint algorithm for estimating amplification efficiency from a polymerase chain reaction dilution series. BMC Bioinformatics 15, 372 (2014). https://doi.org/10.1186/s1285901403724
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1285901403724