- Research article
- Open access
- Published:
Relative transcript quantification by Quantitative PCR: Roughly right or precisely wrong?
BMC Molecular Biology volume 6, Article number: 10 (2005)
Abstract
Background
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.
Results
Analyses of Q-PCR results from a salmon louse starvation experiment show that, even with apparently good raw data, different analytical approaches [1, 2] may lead to opposing biological conclusions.
Conclusion
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.
Communication
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.
We have analysed results from a salmon lice (Lepeophtheirus salmonis) starvation experiment using the 2-ΔΔCT method [1] and the "DART method" adjusting for PCR efficiency differences [2]. When analysed using the 2-ΔΔCT method, our results show that LsTryp1 transcript levels decrease following starvation and return to normal adult level when the louse subsequently gets access to food (Fig. 1). The inclinations obtained when plotting ΔCT or CT against log RNA concentration do not indicate significant differences in PCR efficiencies between LsTryp1 and eEF1α (Table 1). When analysed using the "DART method" the results indicate that LsTryp1 transcript levels decrease 2–3 fold when lice are starved and remain low when lice subsequently get access to food (Fig. 1). The PCR efficiencies, calculated from at least 3 points for each reaction [2], indicated significant differences in PCR efficiency between eEF1α and LsTryp1 in starved and refed lice but not in unstarved lice (Table 1).
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 [10]. 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.
Methods
Salmon lice (Lepeophtheirus salmonis) were reared as earlier described [10]. 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 [10] and the reference gene eEF1α [11] 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 [10]. The RNA purification protocol is previously described [10] 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 [10] and a method adjusting for PCR efficiency differences described by Peirson et al. [2]. The latter analysis was performed partially in the DART-PCR Excel spreadsheet [2]. 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.[10]. 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.
References
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.1262
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/gng073
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.0290023
Freeman WM, Walker SJ, Vrana KE: Quantitative RT-PCR: pitfalls and potential. Biotechniques. 1999, 26: 112-125.
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.e45
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.e36
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-3
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.5530
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.004
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.020
Acknowledgements
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.
Author information
Authors and Affiliations
Corresponding author
Additional information
Authors' contributions
RS and PF conceived the study and designed the Q-PCR assay. RS carried out the assay and analyses, and wrote the first draft of the communication. PF contributed to development of the communication. FN provided expert input for the writing and supervised the study.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
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.
About this article
Cite this article
Skern, R., Frost, P. & Nilsen, F. Relative transcript quantification by Quantitative PCR: Roughly right or precisely wrong?. BMC Molecular Biol 6, 10 (2005). https://doi.org/10.1186/1471-2199-6-10
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1471-2199-6-10