Real-Time Relative qPCR without Reference to Control Samples and Estimation of Run-Specific PCR Parameters from Run-Internal Mini-Standard Curves

Background Real-Time quantitative PCR is an important tool in research and clinical settings. Here, we describe two new approaches that broaden the scope of real-time quantitative PCR; namely, run-internal mini standard curves (RIMS) and direct real-time relative quantitative PCR (drqPCR). RIMS are an efficient alternative to traditional standard curves and provide both run-specific and target-specific estimates of PCR parameters. The drqPCR enables direct estimation of target ratios without reference to conventional control samples. Methodology/Principal Findings In this study, we compared RIMS-based drqPCR with classical quantifications based on external standard curves and the “comparative Ct method”. Specifically, we used a raw real-time PCR dataset as the basis for more than two-and-a-half million simulated quantifications with various user-defined conditions. Compared with classical approaches, we found that RIMS-based drqPCR provided superior precision and comparable accuracy. Conclusions/Significance The obviation of referencing to control samples is attractive whenever unpaired samples are quantified. This may be in clinical and research settings; for instance, studies on chimerism, TREC quantifications, copy number variations etc. Also, lab-to-lab comparability can be greatly simplified.


Introduction
Real-time relative quantitative polymerase chain reaction (qPCR) has long been a favoured principle for relative quantifications of nucleic acid sequences. In essence, an undetectably low amount of a specific nucleic acid target sequence is expanded by PCR to a measurable level. Subsequently, the original amount of the target sequence is calculated from the parameters of the PCR. The basis of these calculations is the classical PCR equation: ð Þ Cq , modified from 1 ½ ðI:1Þ N 0 is the amount of the target sequence before PCR, N Cq is the amount of target after C q -rounds of PCR, and E is the efficiency of the PCR-amplification. Usually, E is assumed constant until the onset of PCR exhaustion. The designation ''relative quantification'' refers to the fact that the amount of the target sequence is estimated relative to that of another (or several [2]). Applying the above relationship, the ratio before PCR of sequences A and B in a given interest-sample (is) can thus be estimated as: To solve the equation above, N Cq s, C q s and Es must be accounted for. Indirect measures of N Cq s are estimated by fluorescence sampling. Various technologies exist (reviewed in [3], [4]), but the common principle of Real-Time PCR is to obtain fluorescence emissions of an intensity proportional to the amount of target at a given point of time [5]. N Cq and the corresponding C q are indirectly defined in the setting of a fluorescence-intensity threshold value. The threshold can be set by various approaches; for example, second-derivatives-maximum method, manually setting, and so on. In mainstream qPCR, E is either assumed to have a value of 1 (the ''comparative Ct method'' or ''2 DDCq '' [6], [7] or estimated target-specifically from a standard curve (SC) [8].
Despite the broad applicability of the technology, several methodological limitations have yet to be addressed. In this paper, we focus on two of these.
A general limitation is that Eq.I.2 cannot be solved trivially offhand. The problem lies in the N Cq,A /N Cq,B -term; in conventional settings, the value of this term is unknown. This may not be problematic as such, provided that the ratio of A and B in the interest-sample is of interest only relative to their ratio in a controlsample (cs): The ''double-ratio'' above can be simplified, eliminating the N Cq,A /N Cq,B -term, if both N Cq and E are preserved for each target between the samples is and cs; that is: (N Cq ) A,is /(N Cq ) B,is = (N Cq ) A,cs / (N Cq ) B,cs and E A,is = E A,cs , E B,is = E B,cs ): ) B happen to be equal, the term disappears, but this situation is unlikely to occur by chance. Simply assuming equality (i.e., ignoring the term entirely) induces proportional errors of the R is -estimate. ''Ensuring'' equality by assigning the same fluorescence-threshold-value for A and B is unreliable. The proportionality factor between fluorescence and sequence-numbers (k) differs widely between targets. For example, the k-value of the fluorophore SYBRGreen depends on the length and sequence of the amplicon as well as other factors, such as salts and temperature [9]. Use of sequence-specific probes is most likely subject to similar shortcomings, in addition to potential target-associated differences in probe-labelling efficiencies and fluorophore bleaching. The most reliable approach is therefore to obtain actual (N Cq ) A /(N Cq ) B -estimates, yet conventional approaches fall short in this respect. Another important limitation of conventional approaches concerns the use of low-capacity machinery. Conventionally, limited instrument capacity forces the investigator to estimate PCR unknowns (such as E) from standard curves analysed in separate runs or assuming a value of 1. This introduces a run-torun variability that inevitably contributes to the error of E. As such, E varies considerably between replicate runs (e.g. .5% [10]). Even tiny errors of the E-estimate are critical. These errors induce disproportionately large errors of R is /R cs or R is because E constitutes the base of the exponential PCR function (Eq. I.1). In this light, errors associated with run-to-run variability of PCR unknowns are highly undesirable.
We hypothesized that run-internal estimation of PCR unknowns (from small amounts of standard curve samples) is superior to run-external estimation. Also, by modifying the composition of standard curve samples, we hypothesized that target ratio estimates can be attained directly without reference to control samples.
Our objective was to deduce the optimal composition of standard curve samples to remedy the limitations of classical qPCR. In the process, we wanted to compare the precision and accuracy of our approach to classical approaches.

Construction of a fusion-PCR product
Blunt-ended PCR products of parts of Human Endogenous retrovirus 1 (ERV1) and TUP-like enhancer of SPLIT 1 (TUPLE1) were produced from genomic DNA by conventional PCRs (Platinum Pfx DNA Polymerase (InVitrogen)). Primer sequences were obtained from Overhauser J et al [11] and Weksberg R et al [12]. PCR products were gel electrophorized. Single bands of expected lengths were excised and PCR products purified (Illustra GFX TM PCR DNA and Gel Band Purification Kit, GE Healthcare). TUPLE1 PCR products were 59-dephosphorylated with rAPID Alkaline Phospatase in supplied buffer (Roche) and purified. Blunt-ended ERV1 PCR products-and 59dephosphorylated TUPLE1 PCR products were ligated (T4 DNA Ligase in T4 DNA Ligase Buffer, New England Biolabs) and purified. Hundred-fold diluted fused PCR products were expanded by standard PCR (ERV1 forward and TUPLE1 reverse primers, respectively). Reactions were gel electrophorized; and single bands consistent with the expected length of the fused PCR product were obtained, purified, and diluted 100-fold before an additional round of PCR, isolation, purification, and dilution. The PCR products were validated by sequencing.

Dilution series
The diluted fusion-PCR product was thoroughly mixed and stored in aliquots (220uC). On three separate days, an aliquot was thawed and used for a 10-fold dilution series in eight steps. Weight data of the pipetted volumes were sampled while the dilution series were made. Each dilution step was performed three times into the same tube (to minimize impact of stochastic errors).

Real-Time PCR
Each dilution series was analyzed by real-time PCR with the primer pairs for ERV1 and TUPLE1. Each primer pair was used in separate runs. Reactions of 20 ml were set up in LightCycler capillaries: 10 ml 26QuantiTect SYBR Green PCR mix (Qiagen), 0.5 mM primers and 8 ml template. PCRs were conducted on a LightCycler 1.0 Instrument (Roche) with the following settings: 15 min at 95uC, 45 amplification cycles (each 15 seconds at 94uC, 20 seconds annealing at 57uC, and 20 seconds at 72uC with endpoint fluorescence detection). Each of the eight concentrations of the dilution series was analysed in four replicates. Six preliminary data sets containing 32 data points each were thus generated.

Raw data sampling
We estimated C q s for the six preliminary data sets by the fitpoints-approach, which gave more linear standard curves than the second derivative maximums-method (see Appendix S1, section 3). In practice, arithmetic baseline adjustment was used, noise bands were set by default, and the lowest (eighth) concentration that yielded C q s for all replicates was excluded (explained below). The remaining data points of the six separate data sets (seven concentrations, each analysed in replicates of four) were used in the ''minimize error''-function of the LightCycler software (version 3.5) for threshold setting. The fit-point number was selected for each of the six data sets as the number providing the smallest error. Relative target concentrations (N 0 s) in the dilution-series-samples were determined from the loaded volumes unless otherwise stated.

Data handling
Data handling was done in Microsoft Excel 2007. Sampled C qand N 0 -data were used to generate standard curves. We excluded data obtained from the eighth dilution step to avoid introduction of heteroscedasticity, which would invalidate conventional linear regression analysis (see Appendix S1, section 1 and Figure S1). The data of the eighth concentrations thus merely represented a safety feature and were not used for further analysis. Standard curves based on relative target concentrations were defined as: ze, e: random error ðM:1Þ Regarding LOG (N 0 ) as the outcome variable and C q as the predictor variable is opposite of the conventional approach and may seem awkward. However, as N 0 of the-samples-to-bequantified are to be estimated as regression estimates from C q the perception makes sense. More importantly, estimating errors of R is is simpler as detailed in later sections and Appendix S1 section 2. Thus, the formula for the regression equation is: With a and b being intercept and slope, respectively. Data of the six 28-sample-standard curves generated from fusion-PCR products are presented in Table 1.

Run-Internal Mini-Standard curves (RIMS) simulation
To examine the usability of RIMS and RIMS-based drqPCR, we simulated a large number of individual quantifications from actual real-time PCR data. First, we combined the raw data (C q and N 0 ) of the six 28-sample standard curves were used to generate 8,694 different RIMS. Each RIMS was based on raw data from two different concentrations of the same 28-sample standard curve data. In other words, each concentration was perceived as a separate RIMS sample (see results). The relative target concentration of a concentration-pair was denoted C. RIMS-concentration-pairs could be chosen in 7?(721)/2 = 21 different ways from each 28-sample standard curve data set and in 6?21 = 126 ways using all six data sets. For each of the seven concentrations of a 28-sample-SC, four C q -replicates were available. To simplify our simulation, we decided to use an equal number of C q -replicates for the two concentrations in the individual RIMS. Thus, for a given RIMS-concentration-pair, 4!/(1!?(421)!)) 2 +4!/(2!?(422)!)) 2 +4!/(3!?(423)!)) 2 +4!/(4!?(424)!)) 2 = 69 different RIMS could be created, for a grand total of 8,694 ( = 126?69) RIMS using all six 28-sample-SC data sets. For each RIMS, the a (intercept) and b (slope) were estimated by conventional linear regression with C q s as the predictor variables and the relative concentrations as outcome variables.

RIMS-based drqPCR simulation
The constructed RIMS were used in simulations of 2,500,848 direct relative quantifications (i.e. quantifications without reference to a control sample). To simplify, unicate quantifications were applied (only one C q -measurement of each target per quantification). The C q -measurements included in the applied RIMS were not quantified using the particular RIMS. Also, a RIMS-set, i.e. one of each target, is necessary for quantification. Therefore, (282m)?(282n) different unicate quantifications cound be constructed from a RIMS-set (m and n being the number of C qmeasurements in the two RIMS). Because more than 18 million RIMS-sets could be constructed in this manner ( = (8,694/2) 2 ), we introduced the constraint that each RIMS-set of ERV1 and TURPLE1 must encompass C q -measurements of samples positioned similarly in the LightCycler carrousel (that is, if sample numbers 1, 2, 5, and 6 were used for the ERV1 RIMS, then the same sample numbers were used for the TURPLE1 RIMS. In this way, 4,347 ( = 8,694/2) RIMS-sets were generated. The a and b of these sets were used with Eq. R.1 and Eq. R.2 for direct relative quantification of ERV1 and TUPLE1 in the remaining samples in the run that was not included in the RIMS. The total number of quantifications were 2,500,848 ( = (7?(721)/2)?(4 2 ?(2821?2) 2 + (4?3/2) 2 ?(2822?2) 2 +4 2 ?(2823?2) 2 +1 2 ?(2824?2) 2 )?3). To facilitate Table 1. Data of the six 28-sample relative standard curves used in the study.

Target
Day evaluation of accuracy and precision, we normalized the individual ratio-estimate by the theoretical ratio. The theoretical ratio was inferred from the quantified samples relative concentration based on data from the construction of the dilution series. This parameter, R is,Norm , ideally equals 1. Any systematic discrepancy reflects inaccuracy whereas stochastic variation around the mean estimate reflects precision.
For the evaluation of double-ratio drqPCR, a combined measure of concentration difference between interest-and virtual control sample in the individual quantification was calculated:

External-SC-based quantifications
The external-SC-based quantifications were done as the 2 DDCqbased quantifications, but with efficiency corrections based on external SCs of Table 1

Statistical analysis
Regression analysis was based on least-squares methods and tdistributions. Probability testing of variance similarity (s 2 1~s 2 2 ) was based on the F-distribution: F obs = largest variance estimate/ smallest variance estimate, degrees of freedom being (f 1 ,f 2 ) and two-tailed p-values = 2?P(F$F obs ). Significance level was set to 0.001 to avoid importance of mass significance.

Direct relative quantitative PCR (drqPCR)
To estimate R is directly (without reference to a control sample), the unknown (N Cq ) A /(N Cq ) B -term presents a challenging. A simple solution is to perform real-time PCR on a sample containing A and B in equal concentrations. The (N Cq ) A /(N Cq ) B -term can then be estimated from the C q -data by rearrangement of Eq. I.2, provided that E-estimates are available (e.g. from standard curves). More elegantly, (N Cq ) A and (N Cq ) B can be inferred from the intercept of As and Bs standard curves (a = LOG (N Cq ), cf. Eq. M.1). If the underlying relation between relative and absolute scales is similar for A and B, then a meaningful estimate of the (N Cq ) A /(N Cq ) B -term can be determined from the intercepts. Similar scales can be attained simply by constructing the standard curves from a common sample containing same concentrations of A and B. If such standard curves is used then it is unnecessary to estimate the actual (N Cq ) A /(N Cq ) B -term. Instead R is can be estimated as a difference of regression estimates (from Eq. M.2): The benefit of using regression estimates is that simple statistics can be applied to determining errors of LOG (R is ) (see Appendix S1 section 2). A sample containing A and B in equal concentrations is sometimes available. If not, the sample can be constructed; for example, by cloning A and B into the same plasmid. Instead of this somewhat cumbersome cloning approach, we recommend joining the targeted PCR products in a fusion-PCR product containing A and B in equal stoichiometry. This fusion-PCR product can be expanded endlessly by PCR and allows standard curves of maximal dilution ranges. A fusion-PCR product of ERV1 and TUPLE1 was constructed as a proof of principle and used in the present study (see the Materials and Methods-section).

Run-Internal Mini Standards (RIMS)
The accuracy of Eq. R.1 hinges on use of valid estimates of slopes and intercepts. These parameters can vary significantly between targets but also between PCR runs of the same target (Tabel I). Therefore, target-specific and run-specific parameters are preferable. We hypothesized that internal standard curves based on fewer samples are preferable over larger, external standard curves. This hypothesis was confronted as follows: Initially, we sought an optimal sample composition strategy for RIMSs. The composition should minimize the errors of regression estimates of LOG ((N 0 ) A,is ) and LOG ((N 0 ) B,is ) in Eq. R.1.
Hellemans et al point out the error of the slope in linear regression is reduced by expanding the range of the dilution and including more measurements points [13]. A similar principle applies to the error of regression estimates (cf Eq. S2.1 in Appendix S1). A large number of measurement points are not desirable with RIMS. However, it is deductable that predictor variable extremes reduce the error more effectively than those close to the predictor variable mean. We therefore based our strategy solely on ''extreme concentrations'' and used only two samples, of relative concentration C, for our RIMSs. Our next step was to determine the importance of C-size and number of replicate analyses of the two RIMS-samples for the precision of aand b-estimates. The total of 8,694 different RIMS was constructed from the raw data of six 28sample standard curve data sets presented in Table 1 (see Material and Methods). The precisions of the RIMS-based aand b-estimates, compared with those of external standard curves, are presented in Figure 1. In particular, we found that precision was increased by increasing values of C. The benefit of increasing the RIMS-sample replicate number was less obvious.

RIMS-based single ratio drqPCR
Next, we investigated the quantitative precision and accuracy of our two approaches when used in combination. In total, 2,500,848 unicate quantifications were determined from the raw data of the six 28-sample data sets. The quantitative precision is summarized in Figure 2. Increases in C and RIMS-sample replicate numbers both generally conferred significant precision improvements. However, the effect of using four as opposed to three RIMS-sample replicates was insignificant. The accuracy was unaffected by C or RIMSsample replicate number and ranged between 94% and 110% of the true target ratios. Double-ratio drqPCR as argued in the introduction, the double-ratio approach (Eq. I.4) can provide R isestimates if the control sample contains A and B in equal concentrations. Using Eq. R.1, LOG (R is ) can be estimated as:

2Þ
A prerequisite is that as and bs are constant between interest and control samples for each target. Some heterogeneity is evident in a comparison of single-ratio to double-ratio based drqPCR (Eq. R.1 and Eq. R.2). Fewer different samples are required in the former approach. This confers fewer sources of errors to R is . However, erroneous b are expectedly more critical in the single-ratio approach, because b is multiplied by larger values (full C q s as opposed to DC q s in the double-ratio approach). In the double-ratio approach, errors in bestimates become less critical when DC q s decreases. A potential drawback of double-ratio drqPCR is that more capacity is required. However, this may be circumvented by using one of the RIMS samples as control sample also. The control sample should be chosen as the RIMS sample providing the smallest DC q s for the individual target in the given quantification. Naturally, one RIMS sample may be closest to A while another is closest to B. Acknowledging, a ''virtual control sample'' can be constructed for the specific interest sample-to-be-quantified comprised of C q -and N 0 -data from the closest RIMS sample (chosen target-specifically). Correction for the situation where (N 0 ) A,cs ?(N 0 ) B,cs is remedied using by the following equation:

Single-vs. double-ratio drqPCR
We compared the precision and accuracy of quantification by single-and double-ratio drqPCR in the 1,190,880 possible quantifications from the six 28-samplestandard curve data sets (C$10 3 ). We expected that increasing the DC q of interest sample and virtual control sample would decrease the precision of doubleratio drqPCR. Data were therefore split according to X, which is a measure of this distance (Eq. M.3). Results are presented in Figure 3. For low Cs (10 3 -10 4 ), the double-ratio approach was more precise when target concentrations of control and interest samples were close (X#2, see Figure 3). Single-ratio-based drqPCR was more precise for larger concentration differences. For large Cs (10 5 to 10 6 ), we observed no remarkable differences between the approaches. The accuracy of both approaches was within 68%.
Comparison of RIMS-based drqPCR, the 2 DDCq -approach, and relative quantification based on external standard curves The data of the six 28-samplestandard curves were used to generate 61,236 and 244,944 different quantifications by the 2 DDCq -approach and based on external standard curves, respectively. Quantifications were based on Eq. I.4. Control samples containing ERV1 and TUPLE1 in equal concentrations were used. Results were normalized by the theoretical ratio and compared with those obtained from RIMS based single-and   (R is,Norm ))). The 8,694 sets of RIMS-derived parameters (a and b) of similar C and RIMS sample replicate number were paired for ERV1 and TUPLE1. Each of the 4,347 paired RIMS-parameter sets were used to calculate all possible run-internal (unicate) LOG (R is ) from the remaining individual C q s not included in the specific RIMS pair. Each LOG (R is ) was normalized by subtracting the logarithmic transformed actual target ratio (determined from the sample's position in the serial dilutions). This provided a total of 2,500,848 LOG (R is, Norm)s. These were sub grouped according to C and RIMS replicate number. The SD of subgroups is illustrated. The number of R is,Norm s in each subgroup is calculable as: (72LOG (C))?(4!/ (m!?(42m)!)) 2 ?3?(282m) 2 . A ten-fold increase of C as well as introduction of an additional RIMS sample replicate provided significantly (p,0.001) better precision with the exceptions indicated by arrows in the figure. doi:10.1371/journal.pone.0011723.g002 double-ratio ( Figure 4). RIMS based drqPCR offered potential for significantly better precision, regardless of X-size. The accuracy of the conventional approaches was overall within 63% but 67% if data were split according to the used standard curves.

Discussion
Real-time relative quantitative PCR in clinical settings is hampered by a lack of accurate and precise approaches to estimate the ratio between nucleic acid sequences without reference to a control sample. Also, conventional approaches for estimating internal PCR parameters are problematic in low capacity PCR machines. This paper concerns the establishment and examination of two new approaches for real-time quantitative PCR; namely RIMS and drqPCR. RIMS concerns estimation of run-internal specific PCR parameters, such as efficiency, from a minimum of samples. The drqPCR is a universal strategy for estimation of ratios directly in the sample, alleviating the need for control samples and therefore ideal for analysis of unpaired samples. We compared RIMS and drqPCR with conventional methods on a common data set. This data set was generated from samples with known target ratios. Therefore, both the precision and accuracy of the approaches could be evaluated.
Separately, RIMS gives target-specific and run-specific estimates of the standard curve's slope and intercept (measures of N Cq and E, respectively). Estimates determined for the specific run are obviously preferable to those determined in separate runs. Runspecific differences of N Cq for the targets can be corrected by inclusion of calibrator samples. However, run-to-run differences in E are not corrected by use of calibrator samples and are critical to the use of external standard curves. It is evident from our data that run-to-run differences of E are to be anticipated ( Table 1). Testing the bs of the same target for significant differences over days discloses dissimilarity in 3 of the possible 6 comparisons (p,0.001 for ERV1 on day 1 vs. 3, TUPLE1 on day 1 vs. 3, and 2 vs. 3).
From Real-Time PCR data, we produced 8,694 individual RIMS. Not surprisingly, the number of RIMS-sample replicates and the value of C were of immense importance for precision ( Figure 1). Compared with standard curves, RIMS provided the potential for attaining estimates of the highest precision. Obviously, a potential explanation could be large run-run external standard curve-variation in our study. Rutledge and Cote [10] used a model comparable to ours, in which a PCR product was serially diluted and subjected to real-time PCR with two different primer sets five times. They reported E CVs of 2.2% and 2.1% for five repeated standard curves for each of two targets. In comparison, E CVs of our study were 2.9% (ERV1) and 2.5% (TUPLE1). The CVs of N Cq in the study of Rutledge and Cote were 19.0% and 14.7%. N Cq -detection was based on constant fluorescence threshold. We based our threshold setting on an error-minimizing strategy for the individual standard curve. Correction of our N Cq -data by fluorescence intensity at C q enables comparison. Threshold-corrected N Cq CVs of our study were 12.4% and 26.3% (data from Table 1). However, the CVs of Rutledge et al do not include the variation associated with the construction of the dilution series, insofar as all standard curves were generated from sequential analysis of only one dilution series. In this light, our raw data are at least comparable in quality.
In our examination of the precision of RIMS-based estimates, the estimates were compared to the corresponding estimate of the internal standard curves. We found a close approximation (that is, a high precision) of the large-C-based RIMS estimates to the 28sample internal standard curves' (Figure 1). We perceive this as  . Comparison of the quantitative precision of drqPCR based on RIMS, external standard curves, or the 2 DDCqapproach. Illustration of the precision of drqPCR (SD (LOG (R is,Norm ))) based on the 2 DDCq -approach (black bars), external standard curves (dark grey bars), and RIMS (exemplified by C = 10 5 and RIMS samples analyzed in duplicates). RIMS parameters were used in single-ratio drqPCR (horizontal broken, black line) and double-ratio drqPCR (white bars). External standard curve-based and 2 DDCq -based drqPCR where by the double-ratio approach only. As control sample data, we used the run-internal, identical-sample-position C q s of ERV1 and TUPLE1. The remaining run-internal combinations of ERV1 and TUPLE1 for a given control sample pair were treated as interest samples. The number of quantifications was 61,236 (27?28?3) for the 2 DDCq -approach and 244,944 (27 2 ?28?3?4) for use of external standard curves. The LOG (R is )s of the conventional approaches were normalized as those determined by RIMS-based drqPCR (Figure 2, legend). Asterisk and crosses indicates significantly (p,0.001) better precision in quantification based on single-ratios and double-ratios compared to conventional approaches, respectively. The 1 at X = 0 indicates the only X-value where external standard curves offered significantly better precision than single-ratio-based drqPCR. doi:10.1371/journal.pone.0011723.g004 indicative of RIMS potential for high accuracy. Furthermore, the result illustrates the redundancy of the intervening data points. RIMS may therefore also be considered as a cost-saving alternative in high-capacity machinery.
The following points should be considered when constructing the RIMS in practice. First, RIMS samples should be selected from the dilution series based on a preliminary standard curve to maximize C ($10 5 ) while preserving linearity. Second, very dilute samples should be avoided completely (Appendix S1, section 1). Third, the chosen samples should be aliquoted and stored. Fourth, two or more replicate C q -estimations of each RIMS sample are preferable. Fifth, the type of template (e.g. PCRproduct, cDNA, or genomic DNA) chosen for RIMS-samples should permit appropriately sized Cs. Direct estimation of nucleic acid sequences ratio is attractive in many settings, yet difficult to obtain. The problem, associated with the (N Cq ) A /(N Cq ) B -term (cf. Eq. I.2), can be dealt with in several ways. The simplest is to ignore it altogether. This is the case if ordinary relative standard curves are used. Such an approach confers systematic errors of magnitudes defined by the inverse of the (N Cq ) A /(N Cq ) B -term. In our example, we would have attained between 1.1-fold and 2.6fold errors for same-day quantifications and as much as 20-fold errors for quantifications across days (cf. N Cq -values in Table 1). Another approach is to rely on experimental conditions assumed to ensure equality of (N Cq ) A and (N Cq ) B . If equality is indeed ensured, the (N Cq ) A /(N Cq ) B -term can be safely ignored. A frequent perception in Real-Time PCR literature is that a similar detection threshold for different targets automatically ensures similar copy numbers at detection [8], [10], [14], [15]. As described earlier, this is unreliable because the proportionality factor between copy numbers and fluorescence intensity (k) can differ between targets. The magnitude of the systematic error associated with the approach can be assessed as the inverse of the targets k-ratio. With this approach, our results would have been off by factor of 1.2 to 1.6 or 2.0 to 2.7 with or without target length corrections, respectively (cf. k-values of Table 1). Although we observed a great deal of variability of the N Cq s across both days and targets, this does not prove the existence of a universal phenomenon. However, our findings do illustrate the importance of expecting variation. This definitely also applies if a common threshold setting is defined. Quantitative-methods with insufficient corrections of (N Cq ) A /(N Cq ) B -term variability can provide results of reasonable precision. But the systematic error is never spotted unless samples of known ratio are quantified. A third and more qualified approach is to use classical absolute quantifications [16]. A drawback of absolute quantification is the considerable effort required to generate samples of known absolute target concentration. Typically, the targets are cloned into plasmids, which are subsequently biologically expanded and purified. The ensuing measurements of DNA concentrations and the derivation of absolute target concentrations are both prone to errors.
Instead, we advocate approaches based on internal, relative standard curves (derived from samples containing the targets of interest in equal amounts) and use of regression estimates. This ensures both target-specific and run-specific corrections of the underlying variable (N Cq ) A /(N Cq ) B . Ligated PCR products of the targets are the ideal theoretical choice. This generally applicable method has the potential for attaining the largest Cs and absolutely defining the targets' stoichiometry. For data handling, we suggest two different approaches: single-ratio or double-ratio-based drqPCR (Eq. R.1 and Eq. R.3, respectively). In the doubleratio-based approach, a virtual control sample is constructed from the RIMS samples already used for b-estimations. Such ''recy-cling'' of RIMS-samples is not problematic, inasmuch as the samples are used to estimate two independent quantities. The known relative concentration between the RIMS-samples (e.g. C A = (N 0 ) A,RIMS1 /(N 0 ) A,RIMS2 ) is used for b-estimations, whereas the relative concentration of the targets of the virtual control sample (e.g. (N 0 ) A,RIMS1 /(N 0 ) B,RIMS1 )) is used in corrections of the (N Cq ) A / (N Cq ) B -term. Steps to reduce the error of R is -estimation are of general importance. This is obvious in clinical settings but also applies to experimental comparison of groups. The larger the variance of a group mean, the more individuals, cell cultures or such, must be included to demonstrate a given significant difference between groups. With this in mind we allocated some efforts to choose the ideal raw data sampling approaches. C qsampling by the FP-approach offered significantly better precision than the second derivatives maximum-approach, whereas estimations of relative target concentrations in standard curves based on weight as opposed to volumes was negligible (and insignificantly) better (see Appendix S1 section 3 and Figure S2, Figure S3, and Figure S4).
We evaluated drqPCR based on RIMS in quantifications of samples containing the quantified targets in known stoichiometry. This model permitted us to evaluate both quantitative precision and accuracy. Not surprisingly, we found that increasing C and the number of RIMS replicates increased the precision significantly ( Figure 2-3). However, C was the most important parameter for improving precision. More than two replicates of each RIMS-sample conferred only minimal improvements of precision. Accuracy was within 68%. Finally, we compared precision and accuracy of drqPCR based on the 2 DDCq -approach, external standard curves, and RIMS. RIMS-based drqPCR demonstrated the largest potential for precision ( Figure 4). The accuracies of the conventional approaches were comparable. The conventional approaches were conducted in a double-ratio manner where R is s were normalised by the ratio derived from analysis of a sample containing the targets in equal concentrations, to maximize accuracy. When using our two approaches in practice, we suggest inclusion of an intermediate RIMS sample. The purpose is dual: to permit evaluation of linearity within the specific RIMS, and to provide more samples available for ''virtual controls''. Also, the C q of samples to be quantified should be estimated in duplicates or more. For practical data handling, we have included an Excel-based spreadsheet in the supplementary material (Algorithm S1).
Use of drqPCR has other beneficial side-effects. The principle renders calibrator samples (or reference-control samples) superfluous. Calibrators are a necessity when interest and control samples are not in the same PCR-run [3]. Their purpose is to correct for run-to-run differences of targets N q -value. Briefly, the calibrator sample contains the targets-to-be-quantified and is PCR-expanded both in runs of controls and interest samples. Subsequently, interest and control data are made comparable by dividing each with the calibrator data of their respective runs. In drqPCR, data of control and interest samples are always immediately comparable, provided that the same RIMS samples are used in runs. Avoidance of calibrators is attractive to minimize the sources of errors of R is . The importance of drqPCR can be stretched further. Another important effect is that results obtained by drqPCR are immediately comparable between different labs. Thus, the problem of lab-to-lab comparability is avoided completely.
In summary, we suggest that RIMS and drqPCR be used separately or combined for relative quantifications of high precision and accuracyThe drqPCR allows determination of R is directly in the sample, and RIMS can replace external standard curves.

Supporting Information
Appendix S1 Found at: doi:10.1371/journal.pone.0011723.s001 (0.12 MB DOC) Figure S1 The association of replicate-C q spread and N 0 C q residuals were calculated (difference of individual C q and mean C q for the target, day and target concentration) and plotted against the sample number (1-8) of the dilution series. Levene's test for equality of C q -residuals' variances between groups where conducted. Inclusion of data of the eighth samples was associated with heteroscedasticity (p: 6?10 211 ) whereas exclusion conferred homoscedasticity (p: 0.66). NB: On this single occasion, we made two modifications to our C q -sampling approach. Firstly, noise bands were set manually for 3 out of 6 data sets. This was necessary to ensure that thresholds were not set in the lower nonlogarithmic phases. Secondly, C q s of all eight concentrations of the dilution series were used in the minimizing error strategy.  Figure S3 Precision of double-ratio-based qPCR by the 2 DDCqapproach (black and dark grey bars) or from external standard curves (light grey and white bars) for variable sizes of X (defined in main text). C q -sampling was by SDM (black and light grey bars) or FP (dark grey and white bars). Quantifications based on C qsampling by FP were more precise for both approaches of efficiency estimations for all X-values. The improvements were significant (p,0.001) except for external standard curves and X of 5 (p: 0.035) and 6 (p: 0.066) (indicated by asterisks). Found at: doi:10.1371/journal.pone.0011723.s004 (0.18 MB TIF) Figure S4 Quantitative precision of single-ratio drqPCR based on RIMS (duplicate analysis) with estimations of C based on pipetted weights (white and dark grey bars) or volumes (black and light grey). Weight as opposed to volume based C-estimation provided minute improvements of precision for all Cs regardless of C q -sampling approach (SDM: black and dark grey, FP: light grey and white). However, improvements were insignificant except where indicated by an asterisk (p: 0.0002). Analysis of RIMS samples in replicates of 1, 3, and 4 demonstrated similar findings.