Assessment of Mathematical Approaches for the Estimation and Comparison of Efficiency in qPCR Assays for a Prokaryotic Model

: Quantitative PCR is a molecular technique for DNA quantification that depends on reaction efficiency and the Ct value (“cycle threshold”). However, the results are dependent on laboratory conditions and mathematical approaches. Thus, the data of 16 genes from Pseudomonas aeruginosa strain AG1 were generated using qPCR to assess the effect of DNA concentration and three mathematical methods (a standard curve and two individual-curve-based approaches called exponential and sigmoidal models) on efficiency and DNA quantification. Differences in efficiency were revealed depending on the mathematical method used; the values were 100% in three out of the four standard curves, but estimations of the expected fold change in DNA serial dilutions were not achieved, indicating the possible overestimation of efficiency. Moreover, when efficiency was compared to DNA concentration, a decreasing trend in efficiency as DNA concentration increased in the reaction was observed in most cases, which is probably related to PCR inhibitors. For all 16 genes at a single DNA concentration, the efficiencies for the exponential model were found in the range of 1.5–2.79 (50–79%), and for the sigmoidal approach, the range was 1.52–1.75 (52–75%), with similar impact on normalized expression values, as indicated by the genes for standard curves. Jointly, DNA concentration and mathematical model choice were demonstrated to impact the estimation of reaction efficiency and, subsequently, DNA quantification when using qPCR.


Introduction
Quantitative polymerase chain reaction (qPCR) is a technique that amplifies and estimates the concentration of a DNA amplicon after each amplification cycle [1].qPCR is considered the gold standard in the field of relative nucleic acid measurements [2].This highly sensitive method is frequently used as an analytical tool with biological and clinical applications.Despite the use of qPCR for more than 30 years since it was invented [3], there are different mathematical approaches to model amplification curves, efficiency, and DNA concentration estimation [2,4].Many methods exist and have already been benchmarked in a previous work [5], but classic approaches based on standard curve, exponential, and sigmoidal models are frequently used in several applications.
As shown in Figure 1A, a classical model of DNA amplification by PCR is represented by a curve after plotting the number of reaction cycles (x-axis) versus fluorescence signal (R n ) accumulation (y-axis) [6].During the first amplification cycles, the low amounts DNA 2024, 4 190 of produced amplicon cannot be detected (the fluorescence signal is below the limit of detection).This is called the baseline phase.Then, an exponential increase is observed, in which a maximum production rate is achieved.Finally, a post-exponential (also known as the plateau) phase is established when the product is no longer produced [4].
As shown in Figure 1A, a classical model of DNA amplification by PCR is represented by a curve after plotting the number of reaction cycles (x-axis) versus fluorescence signal (Rn) accumulation (y-axis) [6].During the first amplification cycles, the low amounts of produced amplicon cannot be detected (the fluorescence signal is below the limit of detection).This is called the baseline phase.Then, an exponential increase is observed, in which a maximum production rate is achieved.Finally, a post-exponential (also known as the plateau) phase is established when the product is no longer produced [4].The efficiency (E) of a PCR is defined as the fraction of target molecules that are copied in one PCR cycle [7].After the assay is finished, calculations of the Ct values ("cycle threshold" is the number of cycles required to reach the arbitrary threshold of the fluorescence signal) are performed at the beginning of the exponential phase and are subsequently used for quantitation purposes, either absolute or relative [6].
Classically, an efficiency value (theoretically ranged from 0-100% or equivalently ranged from 1-2 for form 1 + E) can be estimated from a standard curve, in which serial dilutions of known concentrations of DNA are used to estimate Ct values [2].
Then, a plot of the logarithm of concentrations (x-axis) and Ct values (y-axis) is established to perform a linear regression (Figure 1B), and by using the slope of the regression, the efficiency is then estimated as in [8]: Alternatively, a specific efficiency for each amplification curve can be obtained individually (from now on, we refer to this as "individual-curve-based" approaches) by considering exponential or sigmoidal methods [9].For the first case, the exponential phase of the amplification curve (signal Rn) can be represented as an exponential function of the initial DNA quantity (baseline signal R0), amplification cycle, n, and the efficiency, E, to produce an Rn signal: The efficiency (E) of a PCR is defined as the fraction of target molecules that are copied in one PCR cycle [7].After the assay is finished, calculations of the Ct values ("cycle threshold" is the number of cycles required to reach the arbitrary threshold of the fluorescence signal) are performed at the beginning of the exponential phase and are subsequently used for quantitation purposes, either absolute or relative [6].
Classically, an efficiency value (theoretically ranged from 0-100% or equivalently ranged from 1-2 for form 1 + E) can be estimated from a standard curve, in which serial dilutions of known concentrations of DNA are used to estimate Ct values [2].
Then, a plot of the logarithm of concentrations (x-axis) and Ct values (y-axis) is established to perform a linear regression (Figure 1B), and by using the slope of the regression, the efficiency is then estimated as in [8]: Alternatively, a specific efficiency for each amplification curve can be obtained individually (from now on, we refer to this as "individual-curve-based" approaches) by considering exponential or sigmoidal methods [9].For the first case, the exponential phase of the amplification curve (signal R n ) can be represented as an exponential function of the initial DNA quantity (baseline signal R 0 ), amplification cycle, n, and the efficiency, E, to produce an R n signal: In contrast, under the sigmoidal method, the entire amplification curve (all phases) can be modeled.In this case, the function requires maximum, R max , and minimum, R min , signals, a cycle, n 1/2 , in which half of R max is achieved, and the slope, k, in a quasi-linear part of the amplification [9,10]: For individual curves, modelling can be carried out by optimizing parameters to reduce the differences between the observed values (R n signal in each cycle) and the modeled values.Then, for specific cycles a and b, efficiency can also be calculated [11]: For normalization, a fold change or ratio is established between two experimental conditions, as detailed mathematically in several studies by Pfaffl et al. [8,12,13].In this context, the use of a reference gene with "stable expression" (biologically known as a housekeeping gene) is critical to normalizing the values of a gene of interest (target).In this case, the initial values for target and reference genes are divided to estimate the ratio, r = R 0,T /R 0,R , or normalized values.After the isolation of value R 0 in Equation ( 2), and considering that the threshold value is obtained for R n when n = Ct with specific efficiencies for each gene, the normalization is This normalization can be applied to two different conditions, such as experimental and control, and subsequently, the fold change (FC) as the ratio between normalized values, FC = r experiment /r control , can be obtained.
From Equation (5), when an efficiency of 100% is assumed, the normalization corresponds to the classical 2-delta-Ct for differential expression analysis [14].Theoretically, in this later case of E = 1, the DNA quantity will double in each cycle.Nonetheless, PCR efficiency usually reaches 65-90% due to, for example, differences in reaction inhibitors, enzymes, primers, and probes [2,6,15].Thus, this scenario shows that efficiency calculation is not straightforward, and different parameters and mathematical models can impact the results of the gene expression estimations.In this sense, this is particularly critical due to the exponential nature of PCR, in which the reaction efficiency can have dramatic effects on quantitative determinations [2].
Here, we benchmarked different conditions to estimate efficiency for qPCR, including DNA concentration, amplicon length and mathematical models based on the standard curve and exponential and sigmoidal methods, as well as the ideal case 2-delta-Ct for normalizations.For this, qPCR data from a biological model were used.A total of 16 genes were selected from Pseudomonas aeruginosa AG1 [16].These genes are related to the response to perturbations and are part of previous and further studies to gain biological insights into their regulation at the transcriptomic level [17][18][19].
Jointly, the study aimed to assess the effect of mathematical methods (in silico conditions) and experimental conditions on the quantification of qPCR efficiency using data from a prokaryotic model.

Biological Model and Experimental Analyses
The prokaryotic model Pseudomonas aeruginosa AG1 was used to generate qPCR data.A total of 16 genes found in the assembled genome [20] were considered for primer design and DNA quantification.The primers were designed using Primer-BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/,accessed on 14 January 2020) and were synthesized by Macrogen Inc. (Seoul, Republic of Korea).Details of the primers and amplicon size are presented in Table 1.The standardization of qPCR assays for three genes corresponding to two metallo-βlactamases alleles (VIM-2 and IMP-18) and rpoD as a reference gene [21] are reported in a previous study [19].For the other 13 genes (proC, gcdH, dhcA, braZ, PrtN, pyoS5, prtR, rpoS, lexA, recA, recN, sulA, and lpxA), a qPCR protocol was standardized at the experimental level.DNA was isolated from a bacterial culture using the QIAGEN DNeasy Kit (QIAGEN, UK).A single DNA extract was considered for all assays, including dilutions.The reactions were prepared with 12.5 µL of SYBR green Master Mix (Thermo Scientific™ Inc., Carlsbad, CA, USA), 10 µL of PCR-grade water, 0.25 µL of each primer (concentrations in Table 1), and 2 µL of DNA.A StepOnePlus Real-Time PCR System (Thermo Scientific™, Inc.) was used for thermocycling.The conditions were denaturation at 95 • C (5 min) and 35 amplification cycles of 95 • C (20 s), 60 • C (20 s), and 72 • C (30 s), with data acquisition at 72 • C. For the melting curve, the range of 60-95 • C (increment: 0.03%) was established to detect nonspecific amplicons.After establishing the appropriate experimental conditions, definitive assays with all genes were performed.For all genes but VIM-2 and IMP-18, quantification was performed with a single DNA concentration of 10 ng/µL, including proC as a reference gene [21,22].For genes VIM-2, IMP-18, rpoD, and proC, five DNA concentrations, with serial dilutions (diluent: PCR-grade water) from the highest concentration to obtain 0.01, 0.1, 1, 10, and 100 ng/µL, were used to generate calibration curves (Table 2).Three replicates were considered for each concentration.DNA quantification was performed using NanoDrop equipment (Thermo Scientific™ Inc.).By using the Ct values obtained in each cycle for each assay, the efficiency values were obtained using the exponential and sigmoidal methods (Figure 1C,D).Equations ( 2) and (3) were used as mathematical models to fit the experimental curves using GraphPad Prism (version 5.03, GraphPad Software, San Diego, CA, USA), in which parameter optimization was implemented to minimize the error between the function and the actual data of the qPCR amplification curve (i.e., until the curve fit the experimental data).For example, Figure 1C,D show the actual parameter fitting for the proC gene when [DNA] = 0.1 ng/µL for the exponential and sigmoidal methods.Similarly, a total of 91 qPCR assays were individually fitted by using both individual-curve-based approaches.After parameter optimization, curves were used to estimate efficiency with the cycles around the Ct value of Equation (4) (for example, proC gene with [DNA] = 0.1 ng/µL had Ct = 21.84,so cycles a = 21 and b = 22 were used to estimate efficiency; see Figure 1C,D).

Comparison of Mathematical Approaches for the Estimation of qPCR Efficiency
Replicates used in the standard curves for each gene (IMP-18, VIM-2, proC, and rpoD) were used to assess the effect of DNA concentration on PCR efficiency.A Kruskal-Wallis test (95% confidence, after validation of assumptions) was implemented to determine if there was a significant difference between the template concentration and the efficiency among the different mathematical methods.In a second analysis, data for all the 16 genes (same DNA concentration) were used to compare amplicon size and PCR efficiency estimated by each individual-curve-based approach (95% confidence, using linear regression test).

Fold Change between Points of the Standard Curve
Finally, by using the Ct values obtained from the standard curves and efficiency values per DNA concentration, the fold was calculated against the lower concentration, 0.01 ng/µL.Due to the nature of the standard curve, in which decimal dilutions were made, expected folds of 1, 10, 100, 1000, and 10,000 were compared to the calculated values from the different approaches for estimating qPCR efficiency.The ideal theoretical case with 100% (E = 2) was included as part of the 2-delta-Ct approach.The calculations were derived from Equation ( 5), including a normalization using the reference gene proC at the same DNA concentration, as per Equation (5).Because rpoD was also used for normalization [21][22][23], a preliminary analysis demonstrated lower variability for proC compared to rpoD.

Results
In order to assess the effect of mathematical methods and other parameters on the quantification of the efficiency of qPCR reactions, data from P. aeruginosa AG1 were used.
For the calculation of qPCR efficiency using the standard curve method, standard curves with five 10-fold dilutions (triplicates) were obtained for each of the following genes: IMP-18, VIM-2, proC, and rpoD.Table 2 shows that a consensus E value was obtained for each gene (using Equation (1)).For three out of the four genes, the efficiency value was 2 (i.e., 100%, with R-square values of >98.7%), whereas for VIM-2, the value was 1.9 (90%; R-square: 99.5%).In the case of the individual-curve-based approaches using the same data as above, diverse E values were obtained.For the exponential method, the values ranged between 1.36 and 2.57, while a range of E = 1.30-2.24was found for all cases using the sigmoidal method, with a relevant number of cases higher than the expected maximum value of 2 for both approaches.
For each assay (any row in Table 2), the efficiency and normalized data are dependent on the mathematical model.For example, for VIM-2 with [DNA] = 100 ng/µL, the efficiency values were 1.90 under the standard curve, 1.50 for the exponential method, and 1.43 for the sigmoidal method.When the expression is calculated and normalized according to proC at [DNA] = 100 ng/µL, the values resulted in a high variation, with a coefficient of variation (CV) of CV = 87.4obtained from normalized expression values of 0.8606, 1.6678, 0.2515, and 0.2705 for the 2-delta-Ct, standard curve, exponential, and sigmoidal methods, respectively.Similar patterns were found for all other genes.
When efficiency is compared to DNA concentration (Figure 2), the pattern suggests that reaction efficiency is influenced by the concentration despite maintaining the same amplification conditions for the same gene.A decreasing trend in efficiency can be observed in most cases as DNA concentration increased in the reaction.
using the sigmoidal method, with a relevant number of cases higher than the expected maximum value of 2 for both approaches.
For each assay (any row in Table 2), the efficiency and normalized data are dependent on the mathematical model.For example, for VIM-2 with [DNA] = 100 ng/µL, the efficiency values were 1.90 under the standard curve, 1.50 for the exponential method, and 1.43 for the sigmoidal method.When the expression is calculated and normalized according to proC at [DNA] = 100 ng/µL, the values resulted in a high variation, with a coefficient of variation (CV) of CV = 87.4obtained from normalized expression values of 0.8606, 1.6678, 0.2515, and 0.2705 for the 2-delta-Ct, standard curve, exponential, and sigmoidal methods, respectively.Similar patterns were found for all other genes.
When efficiency is compared to DNA concentration (Figure 2), the pattern suggests that reaction efficiency is influenced by the concentration despite maintaining the same amplification conditions for the same gene.A decreasing trend in efficiency can be observed in most cases as DNA concentration increased in the reaction.This phenomenon is observed for the four genes using specific-curve-based approaches when using either an exponential or sigmoidal method, with statistical support (p < 0.05) in all cases except for the rpoD gene in the exponential method.For the other qPCR data at a single DNA concentration, efficiencies for the exponential model were found in the range of 1.50-1.79,while the range was 1.52-1.75for the sigmoidal approach.For these cases, we also compared efficiency and the amplicon size (Figure 3), in which a stability of efficiency is observed when size varies in the range of 150-200 bp (p > 0.05).This phenomenon is observed for the four genes using specific-curve-based approaches when using either an exponential or sigmoidal method, with statistical support (p < 0.05) in all cases except for the rpoD gene in the exponential method.For the other qPCR data at a single DNA concentration, efficiencies for the exponential model were found in the range of 1.50-1.79,while the range was 1.52-1.75for the sigmoidal approach.For these cases, we also compared efficiency and the amplicon size (Figure 3), in which a stability of efficiency is observed when size varies in the range of 150-200 bp (p > 0.05).
The estimated fold-change among samples with different DNA concentrations was computed with respect to the lower concentration (0.01 ng/µL).The expected values were compared to the predicted values according to the Ct and efficiency values for each mathematical approach.As shown in Table 3, in the theoretical case of 100% efficiency being assumed (2-delta-Ct) or predicted by the standard curve (for all genes but VIM-2), the calculations do not achieve the expected value of the fold and are critical to higher concentrations.For example, for the case of 100 ng/µL (10,000 times the concentration of 0.01 ng/µL), an efficiency of 100% predicts a 5000-fold change after calculations for genes IMP-18, proC, and rpoD.In contrast, when efficiency was 90% (1.9 for VIM-2), the same case provides a fold change of 9878×.The estimated fold-change among samples with different DNA concentrations was computed with respect to the lower concentration (0.01 ng/µL).The expected values were compared to the predicted values according to the Ct and efficiency values for each mathematical approach.As shown in Table 3, in the theoretical case of 100% efficiency being assumed (2-delta-Ct) or predicted by the standard curve (for all genes but VIM-2), the calculations do not achieve the expected value of the fold and are critical to higher concentrations.For example, for the case of 100 ng/µL (10,000 times the concentration of 0.01 ng/µL), an efficiency of 100% predicts a 5000-fold change after calculations for genes IMP-18, proC, and rpoD.In contrast, when efficiency was 90% (1.9 for VIM-2), the same case provides a fold change of 9878×.For the individual-curve-based approaches, both the exponential and sigmoidal methods are greatly impacted by calculations using non-ideal efficiencies.Very extreme values out of the expected folds are found mainly for samples in which the concentration is increased, including log variations.For example, for the case of 1 ng/µL (100 times the concentration of 0.01 ng/µL) for IMP-18, using specific efficiencies predicts a fold change of 400× under an exponential method, which is in contrast to 1343× using the sigmoidal approach.More extreme values are found for the rpoD gene at the same point, in which the exponential approach predicted ~151× in the exponential method and ~17,000× using the sigmoidal approach when 100× is expected.

Discussion
Amplification curves in qPCR can be modeled using mathematical approaches to estimate reaction efficiency [24].In this study, the data of 16 genes from the strain P. aeruginosa AG1 were quantified using qPCR at specific DNA concentrations.Reaction efficiency was estimated using the standard curve and two individual-curve-based approaches: exponential and sigmoidal.The models showed differences in accuracy and reproducibility in efficiency estimation.This represents a critical issue due to a key point in qPCR: accurate quantitation requires that all samples have equal amplification effi-ciency [15] or when considering the limitations of specific assay and statistical parameters when interpreting results.
Notwithstanding the fact that an "ideal value" of 100% efficiency (E = 2) was achieved by almost all standard curves, the mathematical calculations demonstrated that it is not possible to obtain the expected folds of DNA concentrations (with respect to the lower concentration), as evidenced in Table 3.This indicates a possible overestimation of the efficiency using serial dilutions.This comparison was inspired by previous work [11], in which the theoretical foundations demonstrated the reliability of dilutions to predict folds.In our model, this was not achieved, but we used specific Ct data for each condition of the reference gene, which is different from previous work, which used a unique Ct for all dilutions.
For the case of the dilutions in the current study, due to the calculations being based only on efficiency and Ct values, where efficiency is assumed to be perfect, the critical parameter is the Ct value.More precisely, the delta-Ct value between samples (serial dilutions) alters the calculated fold.Thus, the real delta-Ct is affected by the kinetics of the qPCR assay [11].This influences subsequent calculations of the fold and, more importantly, normalizations due to inaccurate Ct values when assuming perfect efficiency.These discordances are like those of other studies, in which up to 42.5% of uncertainty could be found for efficiency estimations using standard curves (in this case, using only one replicate) [7].
Regarding efficiency under individual-curve-based approaches, the efficiency values resulted in variable and lower-than-ideal values when using the sigmoidal method, with an average value of 1.67 (with scarce values with E > 2) for all genes.However, with an average efficiency of 1.85, several values with E > 2 were reported for the exponential model.In both cases, the mean values do not achieve the recommended values of 90% (E = 1.9) [7].In practice, as in our study with the individual methods, the efficiencies obtained from experimental qPCR assays are often considerably less than 100% and are often not equal between individual samples [24].In the case of the expected efficiency, dramatic values were obtained, with predictions that can be considered not theoretically reasonable.As studied before in [2], this is explained by the nature of an exponential model, in which small changes in the input imply a great impact on the calculations of the output.In another study, it was found that a difference of 4% in efficiency could translate into a 400% error in DNA quantification [25].If assumptions of perfect efficiency are not satisfied, the reliability of quantitation is lost [24].Some of the sources of variation in efficiency estimation are related to the number of replicates in standard curves (at least three replicates are recommended for concentration), the instruments, and the sample volumes [7].
In addition, because efficiency calculations are based on different mathematical models in particular contexts, interpretations can be quite different.For example, in the context of a standard curve, efficiency reflects an average value across various dilution levels.In contrast, within a specific curve-fitting model, efficiency represents a model-based estimate of the amplification rate over the PCR cycles.Another consideration concerns any differences in product size among the amplification curve and its effect on the mathematical approach, as was previously reported [26].Long products dominate the first PCR cycles (where Ct is calculated), which is unlike short products that are found at the exponential phase.Thus, efficiency values for standard curves and individual approaches deal with experimental differences that affect their approximations.
Based on these results and mainly considering DNA concentration folds, it seems reasonable to follow the current instructions of the Minimum Information for publication of Quantitative real-time PCR Experiments (MIQE) guidelines, which has recommended the standard curve as the method of choice for the estimation of qPCR efficiency since 2009 [27].Other methods, such as the individual-curve-based approaches, are suggested as complementary options to quantify DNA, although different authors have several opinions.The standard curve remains the most reliable and robust approach to estimating PCR assay efficiency and is broadly accepted by the scientific community according to current agreements [7,27].However, the standard curve has the requirement of several replicates in each concentration to report confident results, including the subsequent drawbacks of being time-consuming and highly expensive [28].Methods using individual amplification curves are reported to give more confidence in estimating efficiency in each reaction and require reduced workloads and costs; however, as was found in this study, the results can be greatly impacted by this method.
In addition, based on the appropriate selection of a specific mathematical method depending on the experimental design for particular contexts [29], some scenarios make sigmoidal and exponential methods candidates for implementation.For example, when an RNA source is scarce (by degradation or during the intracellular cycle of a pathogen, for example), transcript quantification can benefit from the use of the individual-curve methods without the need for large amounts of copy-DNA for building standard curves.
On the other hand, regarding DNA concentration, it was evidenced that efficiency is negatively related to DNA concentration.This phenomenon was also found in another study, in which efficiencies systematically tended to increase as DNA serial dilutions were made [25].A plausible explanation of the dilutional effect is related to PCR inhibitors, which are also diluted with the DNA.PCR inhibitors may affect qPCR assays via different mechanisms, such as disturbing the annealing of primers, affecting DNA polymerase activity, or impairing fluorescence detection [30].Inhibitors from source samples or contamination with chemical compounds (salts, phenol, chloroform, heparin, and/or ethanol) in previous steps, including DNA extraction in our case, result in a lower-than-expected efficiency [25].The application of strategies to purify DNA and eliminate interferences is suggested (but not considered in the protocols presented here).In addition, other explanations of the occurrence of this pattern can be a product of enzyme kinetics-in which substrate concentration can have effects on qPCR efficiency-or the effect of reduced priming accuracy events from an excess of genomic DNA (for example, primer annealing at similar but not identical primer-binding sites along the genome) [31].In the case of amplicon length, which ranged 150-200 bp for all 16 genes, the association was discarded between size and efficiency using individual-curve methods.This is in line with the previous recommendations in which amplicons length <200 bp are suggested [31].
Other experimental conditions can impact qPCR reactions and subsequent analyses, including primer design [32], template quality [25], enzymes (DNA polymerase types and conditions) [33], other reaction components [24], cycling conditions [34], sample handling and pipetting errors, reference genes for normalization [35], and others.Jointly, these parameters will be considered as part of further assays using bacterial models in a more biological context, including the quantification of determinants that we have previously found in the perturbome [17], the SOS response [36], or the response to ciprofloxacin [18] in P. aeruginosa AG1.

Conclusions
In conclusion, qPCR performance is dependent on experimental conditions and mathematical models, as we demonstrated by using 16 different genes for a prokaryotic model, comparing DNA concentration, amplicon length, and three mathematical approaches to assess reaction efficiency.This is critical because it can impact the estimation of DNA concentration and the confidence of qPCR results in different scientific and professional settings.Funding: This work was funded by projects "C1163 pro-NGS 2.0: Protocolos operativos estandarizados de análisis de datos moleculares obtenidos por NGS o afines y de algoritmos de inteligencia artificial en modelos biológicos", Vicerrectoría de Investigación, Universidad de Costa Rica (period 2021-2023) and "C4604 iPAT: Plataforma genómica, bioinformática y de inteligencia artificial

Figure 1 .
Figure 1.Experimental qPCR data and mathematical methods to model amplification and reaction efficiency.Gene proC was used as an example including five DNA concentrations as templates (a single replicate is shown) for amplification (A), to establish calibration curve (B).Model fitting for a specific concentration by exponential (C) and sigmoidal (D) methods are shown, including the points (a) and (b) to determine efficiency (see text).Rn: fluorescence signal; Ct: cycle threshold.

Figure 1 .
Figure 1.Experimental qPCR data and mathematical methods to model amplification and reaction efficiency.Gene proC was used as an example including five DNA concentrations as templates (a single replicate is shown) for amplification (A), to establish calibration curve (B).Model fitting for a specific concentration by exponential (C) and sigmoidal (D) methods are shown, including the points (a) and (b) to determine efficiency (see text).Rn: fluorescence signal; Ct: cycle threshold.

Figure 2 .
Figure 2. Assessment of qPCR reaction efficiency with DNA concentration for four genes and two individual-curve-based mathematical models, for genes IMP-18, VIM-2, proC, and rpoD.

Figure 2 .
Figure 2. Assessment of qPCR reaction efficiency with DNA concentration for four genes and two individual-curve-based mathematical models, for genes IMP-18, VIM-2, proC, and rpoD.

Figure 3 .
Figure 3. Assessment of efficiency of qPCR reaction based on amplicon size for 16 genes.Curvespecific models (exponential and sigmoidal methods) were used to estimate efficiency among all genes with the same DNA concentration as template.

Figure 3 .
Figure 3. Assessment of efficiency of qPCR reaction based on amplicon size for 16 genes.Curvespecific models (exponential and sigmoidal methods) were used to estimate efficiency among all genes with the same DNA concentration as template.

Author
Contributions: J.A.M.-M., C.G. and F.G. participated in the conception and design of the study.M.S.-A., J.A.M.-M., L.R.-M.and D.C.-A.performed experimental assays.M.S.-A.and J.A.M.-M.participated in data analysis.J.A.M.-M.drafted the manuscript.All authors have read and agreed to the published version of the manuscript.

Table 1 .
List of primers for qPCR assays of 16 genes in the bacterial model.

Table 2 .
Comparison of different mathematical methods for efficiency estimation and normalization in qPCR assays.Estimation of qPCR Efficiency Using Different Mathematical Models 2.2.1.Calibration or Standard Curves Four calibration curves were established for gene IMP-18, VIM-2, proC, and rpoD.Efficiency was estimated using Equation (1) (Figure 1B).2.2.2.Individual-Curve-Based Approaches (Exponential and Sigmoidal Methods) * Normalization was undertaken using gene proC by considering the same experiment, in which target genes were quantified at the same DNA concentration.Two independent assays for proC with [DNA] = 10 ng/µL were performed, one inside the standard curves and another with the other genes at a single concentration.** Values with "-" indicate that no standard curve was available.*** For standard curves, a single replicate (median) is reported in the table, but all replicates were considered for the parameter estimation.Abbreviations: [DNA]: DNA concentration; CV: coefficient of variation (among values using different approaches).

Table 3 .
Comparison of the expected and predicted fold change among qPCR assays based on serial dilutions.

Table 3 .
Comparison of the expected and predicted fold change among qPCR assays based on serial dilutions.