|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: JMBJ LØ BKM. Performed the experiments: JMBJ MS. Analyzed the data: JMBJ MSP. Contributed reagents/materials/analysis tools: JMBJ LØ BKM. Wrote the paper: JMBJ MSP MS LØ BKM. Deduced the principles: JMBJ.
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.
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.
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.
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:
N0 is the amount of the target sequence before PCR, NCq is the amount of target after Cq-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 ). 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, NCqs, Cqs and Es must be accounted for. Indirect measures of NCqs are estimated by fluorescence sampling. Various technologies exist (reviewed in , ), 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 . NCq and the corresponding Cq 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ΔΔCq” ,  or estimated target-specifically from a standard curve (SC) .
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 NCq,A/NCq,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 control-sample (cs):
The “double-ratio” above can be simplified, eliminating the NCq,A/NCq,B-term, if both NCq and E are preserved for each target between the samples is and cs; that is: (NCq)A,is/(NCq)B,is=(NCq)A,cs/(NCq)B,cs and EA,is=EA,cs, EB,is=EB,cs):
Although appropriate for paired samples, Eq. I.4 is generally unsuitable if samples are unpaired: the ratio of A and B in a sample can be of immediate interest, and any reference to a control sample can be inconvenient or even meaningless. In these situations, Eq. I.4 can still be useful if the control sample contains A and B in equal numbers. Nonetheless, two problems remain. First, such control samples are not necessarily available. Second, and more importantly, the sources of errors increases by the doubling of sample numbers to be determined with Eq. I.4 compared with Eq. I.2. Ultimately, increased error of the final ratio-estimate is very likely. It may therefore be attractive to actually use Eq. I.2 directly. If so, the inherent (NCq)A/(NCq)B-term must be accounted for. Of course, if (NCq)A and (NCq)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 Ris-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 . 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 (NCq)A/(NCq)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-to-run variability that inevitably contributes to the error of E. As such, E varies considerably between replicate runs (e.g. >5% ). Even tiny errors of the E-estimate are critical. These errors induce disproportionately large errors of Ris/Rcs or Ris 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.
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  and Weksberg R et al . PCR products were gel electrophorized. Single bands of expected lengths were excised and PCR products purified (Illustra GFX™ PCR DNA and Gel Band Purification Kit, GE Healthcare). TUPLE1 PCR products were 5′-dephosphorylated with rAPID Alkaline Phospatase in supplied buffer (Roche) and purified. Blunt-ended ERV1 PCR products- and 5′-dephosphorylated 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.
The diluted fusion-PCR product was thoroughly mixed and stored in aliquots (−20°C). 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).
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 µl were set up in LightCycler capillaries: 10 µl 2× QuantiTect SYBR Green PCR mix (Qiagen), 0.5 µM primers and 8 µl template. PCRs were conducted on a LightCycler 1.0 Instrument (Roche) with the following settings: 15 min at 95°C, 45 amplification cycles (each 15 seconds at 94°C, 20 seconds annealing at 57°C, and 20 seconds at 72°C 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.
We estimated Cqs for the six preliminary data sets by the fit-points-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 Cqs 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 (N0s) in the dilution-series-samples were determined from the loaded volumes unless otherwise stated.
Data handling was done in Microsoft Excel 2007. Sampled Cq- and N0-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:
Regarding LOG (N0) as the outcome variable and Cq as the predictor variable is opposite of the conventional approach and may seem awkward. However, as N0 of the-samples-to-be-quantified are to be estimated as regression estimates from Cq the perception makes sense. More importantly, estimating errors of Ris is simpler as detailed in later sections and Appendix S1 section 2. Thus, the formula for the regression equation is:
With α and β being intercept and slope, respectively. Data of the six 28-sample- standard curves generated from fusion-PCR products are presented in Table 1.
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 (Cq and N0) 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·(7−1)/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 Cq-replicates were available. To simplify our simulation, we decided to use an equal number of Cq-replicates for the two concentrations in the individual RIMS. Thus, for a given RIMS-concentration-pair, 4!/(1!·(4−1)!))2+4!/(2!·(4−2)!))2+4!/(3!·(4−3)!))2+4!/(4!·(4−4)!))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 α (intercept) and β (slope) were estimated by conventional linear regression with Cqs as the predictor variables and the relative concentrations as outcome variables.
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 Cq-measurement of each target per quantification). The Cq-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, (28−m)·(28−n) different unicate quantifications cound be constructed from a RIMS-set (m and n being the number of Cq-measurements 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 Cq-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 α and β 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·(7−1)/2)·(42·(28−1·2)2+(4·3/2)2·(28−2·2)2+42·(28−3·2)2+12·(28−4·2)2)·3). To facilitate 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, Ris,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:
Double-ratio-based quantification was used (Eq. I.4, EA=EB=1). Samples of identical position of ERV1 and TUPLE1 for each serial dilution were used as control samples, and the remaining combinations of samples used as interest samples. This provided 61,236 (272·28·3) double-ratios to be estimated. These were normalized by the theoretical ratio as above. X-values were determined as above.
The external-SC-based quantifications were done as the 2ΔΔCq-based quantifications, but with efficiency corrections based on external SCs of Table 1. Thus, 244,944 (272·28·3·4) normalized double-ratios were calculated.
Regression analysis was based on least-squares methods and t-distributions. Probability testing of variance similarity () was based on the F-distribution: Fobs=largest variance estimate/smallest variance estimate, degrees of freedom being (f1,f2) and two-tailed p-values=2·P(F≥Fobs). Significance level was set to 0.001 to avoid importance of mass significance.
To estimate Ris directly (without reference to a control sample), the unknown (NCq)A/(NCq)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 (NCq)A/(NCq)B-term can then be estimated from the Cq-data by rearrangement of Eq. I.2, provided that E-estimates are available (e.g. from standard curves). More elegantly, (NCq)A and (NCq)B can be inferred from the intercept of As and Bs standard curves (α=LOG (NCq), 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 (NCq)A/(NCq)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 (NCq)A/(NCq)B-term. Instead Ris 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 (Ris) (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).
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 ((N0)A,is) and LOG ((N0)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 . 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 α- and β-estimates. The total of 8,694 different RIMS was constructed from the raw data of six 28-sample standard curve data sets presented in Table 1 (see Material and Methods). The precisions of the RIMS-based α- and β-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.
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 RIMS-sample 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 Ris-estimates if the control sample contains A and B in equal concentrations. Using Eq. R.1, LOG (Ris) can be estimated as:
A prerequisite is that αs and βs 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 Ris. However, erroneous β are expectedly more critical in the single-ratio approach, because β is multiplied by larger values (full Cqs as opposed to ΔCqs in the double-ratio approach). In the double-ratio approach, errors in β-estimates become less critical when ΔCqs 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 ΔCqs 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 Cq- and N0-data from the closest RIMS sample (chosen target-specifically). Correction for the situation where (N0)A,cs≠(N0)B,cs is remedied using by the following equation:
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≥103). We expected that increasing the ΔCq of interest sample and virtual control sample would decrease the precision of double-ratio 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 (103–104), 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 (105 to 106), we observed no remarkable differences between the approaches. The accuracy of both approaches was within ±8%.
The data of the six 28-samplestandard curves were used to generate 61,236 and 244,944 different quantifications by the 2ΔΔCq-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 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 ±3% but ±7% if data were split according to the used standard curves.
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 NCq and E, respectively). Estimates determined for the specific run are obviously preferable to those determined in separate runs. Run-specific differences of NCq 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 βs 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  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 NCq in the study of Rutledge and Cote were 19.0% and 14.7%. NCq-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 NCq -data by fluorescence intensity at Cq enables comparison. Threshold-corrected NCq 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 28-sample internal standard curves' (Figure 1). We perceive this as 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 (≥105) 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 Cq-estimations of each RIMS sample are preferable. Fifth, the type of template (e.g. PCR-product, 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 (NCq)A/(NCq)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 (NCq)A/(NCq)B-term. In our example, we would have attained between 1.1-fold and 2.6-fold errors for same-day quantifications and as much as 20-fold errors for quantifications across days (cf. NCq-values in Table 1). Another approach is to rely on experimental conditions assumed to ensure equality of (NCq)A and (NCq)B. If equality is indeed ensured, the (NCq)A/(NCq)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 , , , . 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 NCqs 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 (NCq)A/(NCq)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 . 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 (NCq)A/(NCq)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 double-ratio-based approach, a virtual control sample is constructed from the RIMS samples already used for β-estimations. Such “recycling” 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. CA=(N0)A,RIMS1/(N0)A,RIMS2) is used for β-estimations, whereas the relative concentration of the targets of the virtual control sample (e.g. (N0)A,RIMS1/(N0)B,RIMS1)) is used in corrections of the (NCq)A/(NCq)B-term. Steps to reduce the error of Ris-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. Cq-sampling 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 ±8%. Finally, we compared precision and accuracy of drqPCR based on the 2ΔΔCq-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 Riss 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 Cq 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 . Their purpose is to correct for run-to-run differences of targets Nq-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 Ris. 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 Ris directly in the sample, and RIMS can replace external standard curves.
(0.12 MB DOC)
The association of replicate-Cq spread and N0Cq residuals were calculated (difference of individual Cq and mean Cq 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 Cq-residuals' variances between groups where conducted. Inclusion of data of the eighth samples was associated with heteroscedasticity (p: 6·10−11) whereas exclusion conferred homoscedasticity (p: 0.66). NB: On this single occasion, we made two modifications to our Cq-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 non-logarithmic phases. Secondly, Cqs of all eight concentrations of the dilution series were used in the minimizing error strategy.
(0.21 MB TIF)
Precision of RIMS based single-ratio drqPCR with Cq-sampling by SDM and FP (black and white bars, respectively) for different values of C (RIMS samples analysed in duplicates). Cq-sampling by FP prompted significantly (p<0.001) better quantitative precision for all values of C.
(0.11 MB TIF)
Precision of double-ratio-based qPCR by the 2ΔΔCq-approach (black and dark grey bars) or from external standard curves (light grey and white bars) for variable sizes of X (defined in main text). Cq-sampling was by SDM (black and light grey bars) or FP (dark grey and white bars). Quantifications based on Cq-sampling 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).
(0.18 MB TIF)
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 Cq-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.
(0.13 MB TIF)
Algorithm for practical use of RIMS based drqPCR.
(1.26 MB XLS)
The authors thank laboratory technician Erik Hagen for performing the sequencing of the TUPLE1-ERV1 fusion PCR product and PhD Martin Roelsgaard for technical discussions.
Competing Interests: The authors have declared that no competing interests exist.
Funding: The authors have no support or funding to report.