Unaccounted uncertainty from qPCR efficiency estimates entails uncontrolled false positive rates

被引:3
作者
Bilgrau, Anders E. [1 ,2 ]
Falgreen, Steffen [1 ]
Petersen, Anders [1 ]
Kjeldsen, Malene K. [1 ]
Bodker, Julie S. [1 ]
Johnsen, Hans E. [1 ,3 ]
Dybkaer, Karen [1 ,3 ]
Bogsted, Martin [1 ,3 ]
机构
[1] Aalborg Univ Hosp, Dept Haematol, Sdr Skovvej 15, DK-9000 Aalborg, Denmark
[2] Aalborg Univ, Dept Math Sci, Fredrik Bajers Vej 7G, DK-9220 Aalborg O, Denmark
[3] Aalborg Univ Hosp, Dept Clin Med, Sdr Skovvej 15, DK-9000 Aalborg, Denmark
关键词
qPCR; Amplification efficiency; Delta-delta Cq; Delta Delta C-q; Error propagation; Efficiency adjusted; POLYMERASE-CHAIN-REACTION; BRANCHING-PROCESS; GENE-EXPRESSION; PCR; QUANTIFICATION; MODEL; PUBLICATION; PROPAGATION; PRECISION; ACCURACY;
D O I
10.1186/s12859-016-0997-6
中图分类号
Q5 [生物化学];
学科分类号
071010 ; 081704 ;
摘要
Background: Accurate adjustment for the amplification efficiency (AE) is an important part of real-time quantitative polymerase chain reaction (qPCR) experiments. The most commonly used correction strategy is to estimate the AE by dilution experiments and use this as a plug-in when efficiency correcting the Delta Delta C-q. Currently, it is recommended to determine the AE with high precision as this plug-in approach does not account for the AE uncertainty, implicitly assuming an infinitely precise AE estimate. Determining the AE with such precision, however, requires tedious laboratory work and vast amounts of biological material. Violation of the assumption leads to overly optimistic standard errors of the Delta Delta C-q, confidence intervals, and p-values which ultimately increase the type I error rate beyond the expected significance level. As qPCR is often used for validation it should be a high priority to account for the uncertainty of the AE estimate and thereby properly bounding the type I error rate and achieve the desired significance level. Results: We suggest and benchmark different methods to obtain the standard error of the efficiency adjusted Delta Delta C-q using the statistical delta method, Monte Carlo integration, or bootstrapping. Our suggested methods are founded in a linear mixed effects model (LMM) framework, but the problem and ideas apply in all qPCR experiments. The methods and impact of the AE uncertainty are illustrated in three qPCR applications and a simulation study. In addition, we validate findings suggesting that MGST1 is differentially expressed between high and low abundance culture initiating cells in multiple myeloma and that microRNA-127 is differentially expressed between testicular and nodal lymphomas. Conclusions: We conclude, that the commonly used efficiency corrected quantities disregard the uncertainty of the AE, which can drastically impact the standard error and lead to increased false positive rates. Our suggestions show that it is possible to easily perform statistical inference of Delta Delta C-q, whilst properly accounting for the AE uncertainty and better controlling the false positive rate.
引用
收藏
页数:14
相关论文
共 42 条
[21]   Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method [J].
Livak, KJ ;
Schmittgen, TD .
METHODS, 2001, 25 (04) :402-408
[22]   No Control Genes Required: Bayesian Analysis of qRT-PCR Data [J].
Matz, Mikhail V. ;
Wright, Rachel M. ;
Scott, James G. .
PLOS ONE, 2013, 8 (08)
[23]   Error propagation in relative real-time reverse transcription polymerase chain reaction quantification models:: The balance between accuracy and precision [J].
Nordgard, Oddmund ;
Kvaloy, Jan Terje ;
Farmen, Ragne Kristin ;
Heikkila, Reino .
ANALYTICAL BIOCHEMISTRY, 2006, 356 (02) :182-193
[24]   Theoretical uncertainty of measurements using quantitative polymerase chain reaction [J].
Peccoud, J ;
Jacob, C .
BIOPHYSICAL JOURNAL, 1996, 71 (01) :101-108
[25]   Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR -: art. no. e36 [J].
Pfaffl, MW ;
Horgan, GW ;
Dempfle, L .
NUCLEIC ACIDS RESEARCH, 2002, 30 (09)
[26]   A new mathematical model for relative quantification in real-time RT-PCR [J].
Pfaffl, MW .
NUCLEIC ACIDS RESEARCH, 2001, 29 (09) :E45
[27]  
Rabinow Paul., 1996, Making PCR: A Story of Biotechnology
[28]  
Rao Xiayu, 2013, Biostat Bioinforma Biomath, V3, P71
[29]   A New Method for Quantitative Real-Time Polymerase Chain Reaction Data Analysis [J].
Rao, Xiayu ;
Lai, Dejian ;
Huang, Xuelin .
JOURNAL OF COMPUTATIONAL BIOLOGY, 2013, 20 (09) :703-711
[30]   Specific expression of miR-17-5p and miR-127 in testicular and central nervous system diffuse large B-cell lymphoma [J].
Robertus, Jan-Lukas ;
Harms, Geert ;
Blokzijl, Tjasso ;
Booman, Marije ;
de Jong, Daphne ;
van Imhoff, Gustaaf ;
Rosati, Stefano ;
Schuuring, Ed ;
Kluin, Philip ;
van den Berg, Anke .
MODERN PATHOLOGY, 2009, 22 (04) :547-555