|
|
||||||||
Molecular Diagnostics and Genetics |
1 Department of Chemistry and Bioscience, Chalmers University of Technology, Gothenburg, Sweden. 2 TATAA Biocenter, 405 30 Gothenberg, Sweden. 3 Department of Medical Biochemistry, Gothenburg University, Gothenburg, Sweden.
aAddress correspondence to this author at: TATAA Biocenter, Medicinaregatan 7B, 405 30 Gothenburg, Sweden. Fax 46-31-7733948; e-mail mikael.kubista{at}tataa.com.
| Abstract |
|---|
|
|
|---|
Methods: We used SYBR green I-based quantitative real-time PCR (QPCR) to measure the properties of reverse transcription reaction for the ß-tubulin, glyceraldehyde-3-phosphate dehydrogenase, Glut2, CaV1D, and insulin II genes, using random hexamers, oligo(dT), and gene-specific reverse transcription primers.
Results: Experimental variation in reverse transcription-QPCR (RT-QPCR) was mainly attributable to the reverse transcription step. Reverse transcription efficiency depended on priming strategy, and the dependence was different for the five genes studied. Reverse transcription yields also depended on total RNA concentration.
Conclusions: RT-QPCR gene expression measurements are comparable only when the same priming strategy and reaction conditions are used in all experiments and the samples contain the same total amount of RNA. Experimental accuracy is improved by running samples in (at least) duplicate starting with the reverse transcription reaction.
| Introduction |
|---|
|
|
|---|
Both real-time PCR and cDNA microarray measurements are highly reproducible (2)(3)(4), but before the expression of any gene can be measured, the mRNA in the sample must be copied to cDNA by reverse transcription. The reverse transcription reaction is not very well understood, and it is expected to be the uncertain step in gene expression analysis. It can introduce biases as a result of effects of the secondary and tertiary structure of mRNA, variation in priming efficiency, and properties of the reverse transcriptase. The yield of the reverse transcription reaction can also be affected by reaction inhibitors present in biological samples (5)(6)(7)(8)(9)(10). To date, no published study has considered the accuracy and precision of the reverse transcription reaction.
Our aim was to study the properties of the reverse transcription reaction, using quantitative real-time PCR (QPCR) 1 as an analytical tool. We investigated the reproducibility, yield, dynamic range, sensitivity, and specificity of the reverse transcription reaction, using the ß-tubulin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (11)(12), Glut2, CaV1D, and insulin II genes, which are expressed differently in a pancreatic ß-tumor mouse cell line (13). We also studied the effect of total RNA concentration on reverse transcription efficiency and compared priming with random hexamers, oligo(dT), and gene-specific primers.
| Materials and Methods |
|---|
|
|
|---|
reverse transcription
The SuperScript II (Invitrogen) reagent set was used for the reverse transcription reaction. We heated 13-µL samples containing total RNA, ultrapure deoxynucleotide triphosphates (dNTPs; Amersham Pharmacia Biotech), and random hexamers (Promega), oligo(dT) (Amersham Pharmacia Biotech), or gene-specific primer (MWG-Biotech) designed to anneal to target mRNA at
80 bp before the start of the PCR product at 65 °C for 5 min to denature the RNA and then chilled the samples on ice for 5 min. We then added Tris-HCl (pH 8.3), KCl, MgCl2, and dithiothreitol (Invitrogen) to a total volume of 19 µL. The random hexamer-primed samples were incubated for 10 min at 25 °C. All samples were then heated to 42 °C for 2 min, and 1 µL of SuperScript II was added to give a final volume of 20 µL containing: 500 µM dNTPs, 50 mM Tris-HCl, 75 mM KCl, 3 mM MgCl2, 5 mM dithiothreitol, 200 U of SuperScript II, and 0.1 µg/µL random hexamers, 0.05 µg/µL oligo(dT), or 1 µM gene-specific primer (Table 1 in the Data Supplement that accompanies the online version of this article at http://www.clinchem.org/content/vol50/issue3/). The reverse transcription reaction was performed at 42 °C for 50 min and was stopped by heating to 70 °C for 15 min.
qpcr
All real-time PCR assays contained 10 mM Tris (pH 8.3), 50 mM KCl, 1 U of Taq polymerase (Sigma-Aldrich), 200 ng/µL bovine serum albumin (MBI Fermentas), 3 mM MgCl2, 0.3 mM dNTPs (Sigma-Aldrich), 1:100 000x SYBR Green I (Molecular Probes), and 400 nM each PCR primer (MWG-Biotech; Table 1 in the online Data Supplement) in 20 µL. The reverse transcription and real-time PCR primers were designed using Primer3 (http://www-genome.wi.mit.edu/cgi-bin/primer/primer3_www.cgi).
Real-time PCR was performed in a LightCycler (Roche Diagnostics) starting with 3 min of preincubation at 95 °C followed by 50 amplification cycles (Table 1
in the online Data Supplement). The threshold cycle (Ct) was determined by use of the maximum-second-derivative function of the LightCycler software. Formation of expected PCR product was confirmed by agarose gel electrophoresis (2%) and melting curve analysis (14).
|
| Results |
|---|
|
|
|---|
qpcr assays
The yield and reproducibility of cDNA synthesis of the ß-tubulin, CaV1D, GAPDH, insulin II, and Glut2 genes were measured by QPCR with SYBR Green I detection (detailed protocols are given in Table 1 in the online Data Supplement). PCR efficiencies (E) and the median and SD (SDQPCR) of the Ct values were calculated for the five genes (Table 1
) (16). For the ß-tubulin, CaV1D, GAPDH, and insulin II assays, SDQPCR was <0.12 cycles (Table 1
). This corresponds to less than (1 + E)0.12
(1 + 0.85)0.12 = 1.08 (assuming an 85% PCR efficiency; Table 1
) or <8% variation in the estimated number of cDNA molecules. For Glut2, which is expressed to a lesser degree than the other genes, SDQPCR = 0.36. This corresponds to a variation of 26% (1.910.36) in the estimated number of cDNA molecules.
efficiency and reproducibility of reverse transcription
Shown in Table 2
are the mean Ct values measured by QPCR for the five genes when the different reverse transcription priming strategies were used. In general, a low Ct value corresponds to high gene expression in the biological sample. However, in our comparative study, the same starting material was used for all reverse transcription reactions; therefore, in this study, a low Ct value indicates a more efficient reverse transcription reaction. A difference of one cycle in Ct between two reverse transcription priming strategies for a particular gene corresponds to a (1 + E)-fold difference in reverse transcription yield. As the data in Table 2
indicate, no reverse transcription priming strategy was best for all five genes. For example, for ß-tubulin, the highest reverse transcription yield was obtained with oligo(dT) primer (Ct = 18.1), and random hexamers gave the lowest yield (Ct = 19.5). For CaV1D, the opposite was true: random hexamers gave the most efficient priming (Ct = 26.5), whereas oligo(dT) performed worst (Ct = 28.8). Also shown in Table 2
are the largest differences in Ct among the various priming strategies for the five genes. The difference ranged from 0.8 cycles for GAPDH to 4.4 cycles for Glut2. Expressed in terms of cDNA molecules, this corresponds to a 59% (1.780.8 = 1.59) difference in cDNA synthesis yield between the best and worst priming strategies for GAPDH and a 17-fold variation (1.914.4 = 17.2) for Glut2. Clearly, the choice of priming strategy can have profound effects on the yield of cDNA synthesis. The yields are evidently also gene dependent.
|
Reverse transcription yields when we used nonmatching or false primers were in all cases low (Table 2 in the online Data Supplement). For the highly expressed genes (ß-tubulin, GAPDH, insulin II, and CaV1D), the reverse transcription yields with false primers were always much lower than with random hexamers, oligo(dT), individual gene-specific primers, or the mixture of the five gene-specific primers, but they were not negligible. In all cases, Ct with false primers was lower than the Ct of the negative control with no primers, indicating that priming of mRNA for reverse transcription is not a particularly stringent reaction. For the Glut2 gene, reverse transcription primers designed to be specific for the other genes gave lower Ct values than the Glut2-specific reverse transcription primer (Table 2 in the online Data Supplement). Clearly, the Glut2 reverse transcription primer hybridizes poorly to Glut2 mRNA under our conditions.
The reproducibility of the reverse transcription priming strategies can be estimated from the five replicates performed for each reverse transcription reaction (Fig. 1 in the online Data Supplement). The SD of a gene expression measurement, i.e., the SD of the determination of the amount of a particular mRNA (SDmRNA), is a weighted sum of the SDs of the QPCR (SDQPCR) and the reverse transcription (SDRT) reactions (16):
![]() | (1) |
accuracy of mRNA quantification
The accuracy of the estimation of gene expression can be improved by running samples in replicate and averaging the measurement results. Because mRNA quantification is performed in two steps, reverse transcription and QPCR, repeats can be done at either one or both steps. When designing experiments, one should consider the experimental accuracy of the two steps. Assuming that samples drawn from the target population are gaussian distributed, the true mRNA concentration (µ) is:
![]() | (2) |
is the estimated mRNA concentration calculated as the mean of n measurements with standard deviation (SDmRNA), for a particular confidence (t) (16). Rearranging Eq. 2
= µ -
) gives:
![]() | (3) |
2 are both expressed as either absolute uncertainties or relative uncertainties. Because both the reverse transcription and QPCR steps contribute to experimental variation, Eq. 3
![]() | (4) |
Although PCR is a cyclic reaction that accumulates errors, its reproducibility is significantly higher than that of the single-step reverse transcription reaction (Table 1
). Using Eq. 4
, we could estimate the sampling error for different experimental setups. As an example, for GAPDH with random hexamer priming (SDQPCR = 0.02; SDRT = 0.13), the sampling error when we divided the test sample into aliquots before QPCR (DRT = 1 and DQPCR = 4) was 0.30. If instead the samples were divided into aliquots before reverse transcription (DRT = 4 and DQPCR = 1), the sampling error was only 0.17. Hence, experimental accuracy was two times higher when the test sample was split into aliquots before the reverse transcription reaction than when it was split before the QPCR. The examples in Fig. 1
show how sampling error depends on experimental design for the different cases of four QPCR aliquots. In general, experimental accuracy is higher when samples are split into aliquots before rather than after the reverse transcription step. The only circumstance in which splitting samples after reverse transcription appears to be advantageous is when SDQPCR is greater than SDRT and cost is an issue.
|
dynamic range of reverse transcription
QPCR analysis has in several studies been shown to have a large dynamic range (17)(18). For quantitative gene expression analysis, the yield of reverse transcription must also be independent of template amount. This was the focus of the present study.
The amount of mRNA as determined by reverse transcription-QPCR analyses is given by:
![]() | (5) |
is the reverse transcription efficiency defined as the fraction of mRNA molecules that are converted to cDNA in the reverse transcription reaction. The exponent in Eq. 5
must be independent of both the total RNA and target mRNA concentrations. This is generally assumed but rarely verified. We studied the dynamic range of reverse transcription using the setup shown in Fig. 2 on the online Data Supplement. A RNA stock solution (1024 ng of total RNA) was diluted in steps of 4 with either water or yeast tRNA, which kept the total RNA concentration constant. cDNA was synthesized with use of either random hexamers or oligo(dT) priming, and the samples were PCR-amplified. The Ct values of the amplification curves of the samples serially diluted with tRNA decreased linearly with dilution factor (Fig. 2
(Eq. 1
|
| Discussion |
|---|
|
|
|---|
The experimental reproducibility of QPCR, as manifested by the low SD of Ct values, is usually very high for common detection chemistries such as SYBR Green I (3), TaqMan probes(12), LightUp probes (2), and Molecular Beacons(26). Only when the number of cDNA molecules is low does SDQPCR increase as a result of interfering primerdimer formation and statistical effects (15)(27). In the present study, primerdimer formation was observed only in samples with Ct values >32. A Ct of 32 corresponds, under our conditions (SYBR Green I detection in LightCycler), to
100 cDNA molecules, which is where statistical effects start to become significant. Of the five genes studied Glut2 had the highest Ct values with the lowest reproducibility (Tables 1
and 2
). In a control experiment, larger amounts of Glut2 template were found to give Ct values as reproducible as the other QPCR assays [SDQPCR(Glut2) = 0.065]. Evidently, the poor quality of Glut2 data can be ascribed to low amounts of Glut2 cDNA. This in turn may be attributable to a low abundance of Glut2 mRNA or low reverse transcription yields for Glut2. Because all priming strategies gave rather high Ct values for Glut2 (Table 2
), it is likely that the samples contained little Glut2 mRNA.
In general, when false primers were used, reverse transcription yields were higher than for the negative controls in which no primer was used (Table 2 in the online Data Supplement). This means that the false primers nonspecifically prime the reverse transcription reaction. Interestingly, the negative controls, which contained all reverse transcription components but no primers, gave some reverse transcription products. Evidently, reverse transcription can be primed by other RNA molecules present in the sample or perhaps by dNTPs (28). The lower assay temperature for reverse transcription compared with PCR is likely to contribute to the low degree of sequence specificity in the priming event. This problem may be overcome by use of thermostable reverse transcriptases. Primer hybridization relies on access to the appropriate target site in the mRNA and may vary substantially because of mRNA folding (29)(30). In general, a higher reverse transcription annealing temperature improves reverse transcription yields by reducing formation of mRNA secondary structures (7)(8).
Our results clearly show that the Ct of a reverse transcription-QPCR assay depends not only on the amount of target mRNA but also on the total RNA concentration (Fig. 2
). We speculate this is attributable to adsorption artifacts that can be eliminated, or at least reduced, by the addition of carrier nucleic acids. We showed that yeast tRNA can be used as carrier, but other polymers, such as linear polyacrylamide, also work well (data not shown).
Gene expression measurements are usually performed as relative measurements (3)(11)(12)(31). Most experimental strategies compare the expression of target genes with the expression of nonregulated reference genes (32). In some cases it has been possible to use the relative expression of two reporter genes as an indicator (2)(33). Calculating the expression ratios, estimated as (1 + E)Ctgene 1/(1 + E)Ctgene 2 (3) for any two of the five genes studied here, we found that it depends on the priming strategy. For example, the expression of ß-tubulin relative to Glut2 is 1.9127.5/1.7919.5 = 628 when measured with use of random hexamers to prime reverse transcription and 1.9131.8/1.7918.8 = 15200 when gene-specific reverse transcription primers were used. Hence, the expression ratio of the two genes differs 15 200/628 = 24-fold when measured with use of two different priming strategies. Clearly, one cannot compare relative gene expression measurements performed under different priming conditions. In fact, it may even be uncertain to compare relative gene expression measurements performed with different batches of random hexamers or oligo(dT) because of batch-to-batch variation. Relative measurements of expression ratios, i.e., comparing the expression of genes in different samples, are possible by compensating for sample-to-sample variation in PCR efficiency by, for example, in situ calibration (2) or kinetic PCR (34)(35). Absolute measurements of expression ratios, i.e., comparing the expression of two genes in a single sample, are not meaningful because of variation in reverse transcription yield unless the experimental system is properly calibrated by use of external standards.
In conclusion, we show that experimental variation in reverse transcription-QPCR is mainly attributable to the reverse transcription step. The efficiency of the reverse transcription reaction depends on the priming strategy and also varies among different genes. The efficiency also depends on total RNA concentrations. When performing gene expression measurements by reverse transcription-QPCR we recommend (a) extensive optimization of the reverse transcription reaction, (b) running the experiment in at least duplicate starting with the reverse transcription step, (c) adjusting the total RNA concentrations to be the same in all samples by adding carrier, and (d) always using the same reverse transcription priming strategy and reaction conditions in experiments that are to be compared.
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
and
immunoglobulin light chain expression. Clin Chem 2003;49:51-59.The following articles in journals at HighWire Press have cited this article:
![]() |
M. K. Udvardi, T. Czechowski, and W.-R. Scheible Eleven Golden Rules of Quantitative RT-PCR PLANT CELL, July 1, 2008; 20(7): 1736 - 1737. [Full Text] [PDF] |
||||
![]() |
L. Metwally, D. J. Fairley, P. V. Coyle, R. J. Hay, S. Hedderwick, B. McCloskey, H. J. O'Neill, C. H. Webb, W. Elbaz, and R. McMullan Improving molecular detection of Candida DNA in whole blood: comparison of seven fungal DNA extraction protocols using real-time PCR J. Med. Microbiol., March 1, 2008; 57(3): 296 - 303. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Zheng, Y. Luo, and G. K. McMaster Sensitive and Quantitative Measurement of Gene Expression Directly from a Small Amount of Whole Blood Clin. Chem., July 1, 2006; 52(7): 1294 - 1302. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. P. Choe, A. Kato, S. Hirose, C. Plata, A. Sindic, M. F. Romero, J. B. Claiborne, and D. H. Evans NHE3 in an ancestral vertebrate: primary sequence, distribution, localization, and function in gills Am J Physiol Regulatory Integrative Comp Physiol, November 1, 2005; 289(5): R1520 - R1534. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. K.C. Brolen, N. Heins, J. Edsbagge, and H. Semb Signals From the Embryonic Mouse Pancreas Induce Differentiation of Human Embryonic Stem Cells Into Insulin-Producing {beta}-Cell-Like Cells Diabetes, October 1, 2005; 54(10): 2867 - 2874. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. D. Schirman-Hildesheim, T. Bar, N. Ben-Aroya, and Y. Koch Differential Gonadotropin-Releasing Hormone (GnRH) and GnRH Receptor Messenger Ribonucleic Acid Expression Patterns in Different Tissues of the Female Rat across the Estrous Cycle Endocrinology, August 1, 2005; 146(8): 3401 - 3408. [Abstract] [Full Text] [PDF] |
||||
![]() |
S A Bustin, V Benes, T Nolan, and M W Pfaffl Quantitative real-time RT-PCR - a perspective J. Mol. Endocrinol., June 1, 2005; 34(3): 597 - 601. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Stahlberg, M. Kubista, and M. Pfaffl Comparison of Reverse Transcriptases in Gene Expression Analysis Clin. Chem., September 1, 2004; 50(9): 1678 - 1680. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |