Relative transcript quantification by Quantitative PCR: Roughly right or precisely wrong?
© Skern et al; licensee BioMed Central Ltd. 2005
Received: 01 December 2004
Accepted: 26 April 2005
Published: 26 April 2005
When estimating relative transcript abundances by quantitative real-time PCR (Q-PCR) we found that the results can vary dramatically depending on the method chosen for data analysis.
The results emphasise the importance of being cautious when analysing Q-PCR data and indicate that uncritical routine application of an analytical method will eventually result in incorrect conclusions. We do not know the extent of, or have a universal solution to this problem. However, we strongly recommend caution when analysing Q-PCR results e.g. by using two or more analytical approaches to validate conclusions. In our view a common effort should be made to standardise methods for analysis and validation of Q-PCR results.
Reverse transcription (RT) followed by quantitative polymerase chain reaction (Q-PCR) is at present the most sensitive method for transcript abundance measurement. However, there are many sources of errors, both when purifying RNA, performing the RT reaction and during the PCR setup [3, 4]. Q-PCR utilises optical measurement of generated amplicons to survey PCR amplifications. It is common to derive the initial template concentration from the number of amplification cycles required for a signal to reach a threshold chosen by the investigator [1, 2, 5]. In relative quantification the expression of a target gene is stated relative to a standard gene, which is assumed to be constitutively and uniformly expressed. One popular approach, the 2-ΔΔCT method, assumes ≈100% efficient target and standard gene PCR reactions given that the results conform to certain criteria [1, 5]. In recognition of the fact that PCR efficiencies may vary between runs or between target and standard genes, other numerous methods have emerged that calculate template concentrations using amplification simulations or PCR efficiencies derived from CT values or fluorescence data [2, 6–9]. We here present the results of a case study showing that the interpretation of results may vary dramatically with the chosen method for data analysis.
Quality assessment of the retrieved data. Quality assessment of the retrieved data. For the 2-ΔΔCT results the table shows inclination and R2 for ΔCT plotted against log RNA concentration and inclinations and R2 for CT plotted against log RNA concentration for eEF1α and LsTryp1. For the DART-PCR results the table shows PCR-efficiencies for eEF1α and LsTryp1 calculated by DART-PCR and the p-value (one-way ANOVA) for the hypothesis that there is no difference between the efficiencies.
ΔCT vs. log [RNA] (R2)
eEF1α CT vs. log [RNA] (R2)
LsTryp1 CT vs. log [RNA] (R2)
By intuition it appears that surveying PCR efficiencies using several measured fluorescence points from each PCR reaction, as done using the "DART method", is superior to using one point from each reaction, as done when comparing ΔCT values using the 2-ΔΔCT method. However, since PCR efficiencies calculated using the "DART method" exceed 100% in some instances, it is clear that this approach also has weaknesses. In the present example (Fig. 1) we would not have more confidence in one method than the other unless we had data from supplementary methods (e.g. microarrays) to support this. Consequently these data indicate that LsTryp1 transcript levels decrease when lice are starved, which is in accordance with the alleged digestive function of the encoded protein . However, since the result varies between the "DART-method" and the 2-ΔΔCT method, we are unable to determine how transcription is regulated after lice resume feeding. Thus, despite the fact that both the 2-ΔΔCT method and the "DART-method" are theoretically sound given a number of assumptions [1, 2], we may be mislead when these assumptions are not fulfilled.
All strategies for analysing Q-PCR data are based on a number of assumptions, and due to experimental errors none or few of these assumptions will be fulfilled entirely. Unfortunately, it is not always obvious when assumptions are broken to a degree that invalidates the conclusions. Since the sources of potential problems are diverse, no simple solution is available. Therefore we do not offer a universal analytical approach that can be applied to any given set of data and ensure a correct conclusion. Rather, we suggest investigators to urge caution when analysing results and hope that future discussions will lead to a more unified approach to Q-PCR data analysis and improved reliability of published results.
Salmon lice (Lepeophtheirus salmonis) were reared as earlier described . After development to the adult stage, 15 lice were removed with forceps from their anaesthetised (80 μg/ml benzocaine) salmon hosts (Salmo salar) and 3 lice were stored in RNA later (Ambion). The remaining 12 lice were starved in incubators with flowing seawater for 14 days. After starvation, 3 lice were sampled and stored as described above, and the remaining 9 lice were put in a tank with uninfected salmon where they could settle on their salmon hosts and resume feeding. After 15 days on their new hosts 3 lice were sampled and stored as described above. The experimental procedures were carried out in accordance with national regulations for use of animals in scientific research.
The transcript levels of LsTryp1  and the reference gene eEF1α  in 1 selected unstarved, starved and refed lice were determined by quantitative real time PCR carried out with 3 parallels at 5 sequential 2-fold dilutions as previously described . The RNA purification protocol is previously described  and cDNA syntheses were performed using MultiScribe™ according to the manufacturers recommendations (Applied Biosystems). The Q-PCR results were analysed by the 2-ΔΔCT method as earlier described  and a method adjusting for PCR efficiency differences described by Peirson et al. . The latter analysis was performed partially in the DART-PCR Excel spreadsheet . When using the 2-ΔΔCT method, at least 2 parallels were required at each dilution. Parallels were removed when the CT value differed more than 0.3 (CT<32) or 0.4 (CT = 32) from the most similar parallel at the same dilution. At least 4 dilutions were required for each stage. The resulting data were calibrated to unstarved lice and analysed as described by Kvamme et al.. When using the "DART-method", dilutions were removed when PCR efficiency differed significantly (one way ANOVA, α = 0.05) from the other dilutions. The signal corresponding to the initial template concentration (R0) was derived using the average PCR efficiency for LsTryp1 and eEF1α when the PCR efficiencies were not significantly different (one way ANOVA, α = 0.05). When the PCR efficiency differed significantly, R0 was calculated using individual gene specific mean efficiencies. The mean R0 for each dilution of LsTryp1 was normalised to corresponding eEF1α values. The normalised R0 values were calibrated to the values for unstarved lice. 95% confidence intervals (CI) were derived from normalised R0 values.
We would like to thank Dr. Mette Mauritzen, Dr. Kevin Glover and Bjørn Olav Kvamme for constructive criticism of the communication. This work has been supported by Intervet International bv and the Norwegian Research Council.
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262View ArticlePubMed
- Peirson SN, Butler JN, Foster RG: Experimental validation of novel and conventional approaches to quantitative real-time PCR data analysis. Nucleic Acids Res. 2003, 31 (14): e73- 10.1093/nar/gng073PubMed CentralView ArticlePubMed
- Bustin SA: Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol. 2002, 29: 23-39. 10.1677/jme.0.0290023View ArticlePubMed
- Freeman WM, Walker SJ, Vrana KE: Quantitative RT-PCR: pitfalls and potential. Biotechniques. 1999, 26: 112-125.PubMed
- Livak KJ: ABI Prims 7700 Sequencer detection system. User bulletin 2. PE Applied Biosystems. 1997, 1-36.
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: 2002-2007. 10.1093/nar/29.9.e45. 10.1093/nar/29.9.e45View Article
- Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30: 1-10. 10.1093/nar/30.9.e36View Article
- Liu W, Saint DA: Validation of a quantitative method for real time PCR kinetics. Biochem Biophys Res Commun. 2002, 294: 347-353. 10.1016/S0006-291X(02)00478-3View ArticlePubMed
- Liu W, Saint DA: A new quantitative method of real time reverse transcription polymerase chain reaction assay based on simulation of polymerase chain reaction kinetics. Anal Biochem. 2002, 302: 52-59. 10.1006/abio.2001.5530View ArticlePubMed
- Kvamme BO, Skern R, Frost P, Nilsen F: Molecular characterisation of five trypsin-like peptidase transcripts from the salmon louse (Lepeophtheirus salmonis) intestine. Int J Parasitol. 2004, 34: 823-832. 10.1016/j.ijpara.2004.02.004View ArticlePubMed
- Frost P, Nilsen F: Validation of Reference Genes for Transcription Profiling in the Salmon Louse, Lepeophtheirus Salmonis, by Quantitative Real-Time Pcr. Veterinary Parasitology. 2003, 118: 169-174. 10.1016/j.vetpar.2003.09.020View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.