BMC Bioinformatics

official impact factor 3.03

Open Access Methodology article

Simultaneous fitting of real-time PCR data with efficiency of amplification modeled as Gaussian function of target fluorescence

Anke Batsch1, Andrea Noetel1, Christian Fork1, Anita Urban1, Daliborka Lazic1, Tina Lucas1, Julia Pietsch1, Andreas Lazar1, Edgar Schömig1,2 and Dirk Gründemann1,2*

Author Affiliations

1 Department of Pharmacology, University of Cologne, Gleueler Straße 24, 50931 Cologne, Germany

2 Center for Molecular Medicine (CMMC), University of Cologne, Joseph-Stelzmann-Straße 52, 50931 Cologne, Germany

For all author emails, please log on.

BMC Bioinformatics 2008, 9:95 doi:10.1186/1471-2105-9-95

Published: 12 February 2008

Abstract

Background

In real-time PCR, it is necessary to consider the efficiency of amplification (EA) of amplicons in order to determine initial target levels properly. EAs can be deduced from standard curves, but these involve extra effort and cost and may yield invalid EAs. Alternatively, EA can be extracted from individual fluorescence curves. Unfortunately, this is not reliable enough.

Results

Here we introduce simultaneous non-linear fitting to determine – without standard curves – an optimal common EA for all samples of a group. In order to adjust EA as a function of target fluorescence, and still to describe fluorescence as a function of cycle number, we use an iterative algorithm that increases fluorescence cycle by cycle and thus simulates the PCR process. A Gauss peak function is used to model the decrease of EA with increasing amplicon accumulation. Our approach was validated experimentally with hydrolysis probe or SYBR green detection with dilution series of 5 different targets. It performed distinctly better in terms of accuracy than standard curve, DART-PCR, and LinRegPCR approaches. Based on reliable EAs, it was possible to detect that for some amplicons, extraordinary fluorescence (EA > 2.00) was generated with locked nucleic acid hydrolysis probes, but not with SYBR green.

Conclusion

In comparison to previously reported approaches that are based on the separate analysis of each curve and on modelling EA as a function of cycle number, our approach yields more accurate and precise estimates of relative initial target levels.