- Research article
- Open Access
Validation of reference genes as internal control for studying viral infections in cereals by quantitative real-time RT-PCR
- Jana Jarošová1 and
- Jiban K Kundu1Email author
https://doi.org/10.1186/1471-2229-10-146
© Jarošová and Kundu; licensee BioMed Central Ltd. 2010
- Received: 16 September 2009
- Accepted: 15 July 2010
- Published: 15 July 2010
Abstract
Background
Reference genes are commonly used as the endogenous normalisation measure for the relative quantification of target genes. The appropriate application of quantitative real-time PCR (RT-qPCR), however, requires the use of reference genes whose level of expression is not affected by the test, by general physiological conditions or by inter-individual variability. For this purpose, seven reference genes were investigated in tissues of the most important cereals (wheat, barley and oats). Titre of Barley yellow dwarf virus (BYDV) was determined in oats using relative quantification with different reference genes and absolute quantification, and the results were compared.
Results
The expression of seven potential reference genes was evaluated in tissues of 180 healthy, physiologically stressed and virus-infected cereal plants. These genes were tested by RT-qPCR and ranked according to the stability of their expression using three different methods (two-way ANOVA, GeNorm and NormFinder tools). In most cases, the expression of all genes did not depend on abiotic stress conditions or virus infections. All the genes showed significant differences in expression among plant species. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), beta-tubulin (TUBB) and 18S ribosomal RNA (18S rRNA) always ranked as the three most stable genes. On the other hand, elongation factor-1 alpha (EF1A), eukaryotic initiation factor 4a (EIF4A), and 28S ribosomal RNA (28S rRNA) for barley and oat samples; and alpha-tubulin (TUBA) for wheat samples were consistently ranked as the less reliable controls.
The BYDV titre was determined in two oat varieties by RT-qPCR using three different quantification approaches. There were no significant differences between the absolute and relative quantifications, or between quantification using GAPDH + TUBB + TUBA +18S rRNA and EF1A + EIF4A + 28S rRNA. However, there were discrepancies between the results of individual assays.
Conclusions
The geometric average of GAPDH, 18S rRNA and TUBB is suitable for normalisation of BYDV quantification in barley tissues. For wheat and oat samples, a combination of four genes is necessary: GAPDH, 18S rRNA, TUBB and EIF4A for wheat; and GAPDH, 18S rRNA, TUBB and TUBA for oat is recommended.
Keywords
- Reference Gene
- Absolute Quantification
- Wheat Sample
- Pairwise Variation
- Barley Yellow Dwarf Virus
Background
Barley yellow dwarf virus (BYDV) is one of the economically most important viral diseases of cereals worldwide. It can cause significant yield losses in major cereal crops like wheat, barley, rice, maize, oat and ryegrass [1]. Virtually all species of the family Poaceae (Graminae) can be infected, providing > 150 species as putative sources of these viruses [2].
In many cases, the virus titre in plants and its ability to multiply need not correspond with symptomatic manifestation of the infection. For example, this can be the case with tolerant plants [3], plants with high nitrogen uptake [4], and can also be dependent on environment conditions [5]. Furthermore, it is often a combination of different viral genes' expression that determines the severity of symptoms [6]. Quantification by quantitative real-time PCR (RT-qPCR) makes it possible to monitor the titre of all the different viral genes as well as to quantify the virus particles. The determination of individual viral genes' expression can greatly enhance our knowledge and reveal many aspects of disease aetiology and virus ecology.
However, for such methods, normalisation is required to correct for any variation in RNA integrity, reverse transcription efficiency, and initial sample amount among different samples [7]. Different approaches have been proposed to normalise measurements of expression levels [8], but this is generally done using an internal 'reference gene', under the assumption that this has a constant level of expression in the chosen tissue, is not affected by the treatment and has no inter-individual variability [9]. Therefore, the most prominent problem in quantitative RNA expression analysis is the selection of an appropriate reference gene. Misinterpretation of data occurs when expression measures are erroneously normalised to a subset of mRNAs that are subject to strong regulation [10, 11]. While it seems unreasonable that the transcription of any gene in a living cell is absolutely resistant to cell cycle fluctuations or nutrient status, it is important to identify candidate genes that are at least minimally regulated during the individual experiment to allow the accuracy of RNA transcription analysis that real-time PCR offers. The correct reference genes can be selected by evaluating data from RT-qPCR with statistical algorithms such as GeNorm [12], BestKeeper [13] or NormFinder [14]. Commonly used reference genes for normalisation of RT-qPCR data in plants include ACTβ, TUBA, TUBB, EF-1-α, EIF4A, GAPDH, and 18S, 25S and 28S rRNAs [15–19].
In this paper, commonly used reference genes' expressions in wheat, barley and oats under abiotic and biotic [Barley yellow dwarf virus (BYDV) infection] stress were compared. Furthermore, an assay for BYDV relative quantification in these cereal species was developed. Different approaches of RT-qPCR were tested in a case study that quantified the titre of BYDV in tolerant and susceptible oat varieties.
Results and discussion
C q values distribution of candidate reference genes in different species. The values are given as real-time PCR quantification cycle (Cq) values. The boxes represent the lower and upper quartiles with medians; the whiskers illustrate the 10-90 percentiles of the samples. All Cq values significantly differed between species (one-way ANOVA). (W) = wheat samples (60 samples), (B) = barley samples (60 samples), (O) = oat samples (60 samples).
Furthermore, the titre of BYDV in two oat varieties (BYDV susceptible and BYDV tolerant) was quantified by absolute quantification with cloned standards and no reference gene, and by relative quantification with the reference genes that performed best and worst in the NormFinder and GeNorm analyses.
ANOVA results
Two-way ANOVA results
Gene | Fitness | Species | Interaction |
---|---|---|---|
28S rRNA | ns | *** | ns |
18S rRNA | ** | *** | ** |
TUBA | * | *** | ** |
TUBB | * | *** | ** |
EIF4A | ns | *** | * |
EF1A | ns | *** | ns |
GAPDH | * | *** | ** |
Results of Bonferroni tests
Gene | Healthy vs. Abiotic stress | Healthy vs. Virus Infected | Abiotic stress vs. Virus infected | ||||||
---|---|---|---|---|---|---|---|---|---|
barley | oat | wheat | barley | oat | wheat | barley | oat | wheat | |
28S rRNA | ns | ns | ns | ns | ns | ns | ns | ns | ns |
18S rRNA | *** | *** | *** | *** | ns | ns | ns | * | *** |
TUBA | ns | ns | *** | ns | ns | ** | ns | ns | *** |
TUBB | ns | ns | ns | ns | ns | ** | ns | ns | *** |
EIF4A | ns | ns | ns | ns | ns | ns | ns | ns | ns |
EF1A | ns | ns | ns | ns | ns | ns | ns | ns | ns |
GAPDH | ns | ns | *** | ns | ns | *** | ns | ns | *** |
GeNorm Results
Ranking of candidate reference genes and choice of best pair of reference genes by NormFinder and GeNorm tools
gene | GeNorm stability value (M) | NormFinder stability value | |||||||
---|---|---|---|---|---|---|---|---|---|
wheat | barley | oat | all | wheat | barley | oat | all | ||
GAPDH | 1,09 | 1,06 | 1,08 | 2,0 | 0,088 | 0,112 | 0,035 | 0,074 | |
TUBB | 1,14 | 1,27 | 1,12 | 2,0 | 0,041 | 0,085 | 0,108 | 0,098 | |
TUBA | 2,13 | 1,54 | 1,59 | 4,43 | 0,785 | 0,201 | 0,205 | 0,289 | |
18S rRNA | 1,48 | 1,27 | 1,14 | 2,15 | 0,076 | 0,158 | 0,125 | 0,248 | |
28S rRNA | 1,82 | 1,71 | 1,90 | 2,78 | 0,378 | 0,354 | 0,258 | 0,487 | |
EIF4a | 1,48 | 2,81 | 2,59 | 2,82 | 0,352 | 0,489 | 0,301 | 0,499 | |
EEF1A | 2,08 | 2,48 | 2,49 | 2,87 | 0,181 | 0,187 | 0,355 | 0,241 |
Stability indices (A), and pairwise variation to determine the optimal number of normalisation genes (B) calculated with GeNorm, and stability indices calculated with NormFinder (C). (A) and (C): The stability indices are shown for all species combined, and also barley, wheat and oats individually. The stability of gene expression is inversely proportional to the stability index, so least stable genes have the lowest achieved values and vice versa. (B): The optimal number of genes was determined separately for barley, wheat and oats as well as for all samples combined (all). The recommended cut-off value under which there is no need for inclusion of another gene is 0.15.
Theoretically, the high expression stability of a gene indicates that the use of a single reference gene is appropriate. However, for many studies, no single gene may be adequate and may require normalisation with two or more stable reference genes. Therefore, pairwise variations were calculated using GeNorm for each data set to determine the optimal number of reference genes for normalisation. For this, first the normalisation factors (NF) were calculated for the most stable reference genes (with lowest M-value) and then for other genes by stepwise inclusion of the reference gene that remains most stable. Subsequently, pairwise variations of NFn and NFn+1 were calculated, reflecting the effect of including additional (n + 1) gene.
For wheat and oats, four genes with the most stable expression were optimal for reliable normalisation with a pairwise variation cut-off value of 0.15, whereas for barley three genes were necessary (Fig. 2B). Below this cut-off value there is no need for inclusion of an additional reference gene. It was apparent that the inclusion of a fourth/fifth gene had no significant effect on pairwise variation. However, an attempt to calculate the optimal number of reference genes for all samples combined failed; even if all the reference genes were taken into account, the pairwise variation exceeded the cut-off value of 0.15 (Fig. 2B). The quantification might not be accurately normalised even if all the tested genes were used in the assay. This was additional evidence of the infeasibility of comparing BYDV between different cereal species by relative quantification, and also further confirmed the suspicion raised by analysis of raw Cq values that the expression of the tested genes differed among species.
NormFinder results
NormFinder is another approach to assess the stability of expression of reference genes [14]. In the present study, NormFinder indicated that the genes with the most stable levels of transcript abundance were GAPDH and TUBB for all the samples taken as one set and for barley and oat samples, and TUBB and 18S rRNA for wheat samples (Table 3, Fig. 2C). These were almost the same genes as those identified by GeNorm, with the exception that 18S rRNA performed better than GAPDH for wheat samples in NormFinder. The least stable genes, according to NormFinder, were EIF4A for all the samples analysed as one set and also for barley, TUBA for wheat samples, and EF1A for oat samples. Therefore, there was a very strong correlation between the results obtained from GeNorm and NormFinder, despite the fact that the methods of calculation are fundamentally different. NormFinder ranks reference genes according to the least estimated intra- and inter-group variation, which is more effective to control the influence of co-regulation of reference genes. NormFinder can account for heterogeneity in the tested samples, such as different treatment groups, and so distinguishes between stability and bias.
The wheat samples differed more from barley and oat samples than these two species differed from each other. This was especially the case for TUBA which performed well for barley and oats but very poorly for wheat. However, there were significant differences in expression of most of the genes, even between oats and barley, and therefore no intra-species comparison by real-time relative RT-qPCR can be recommended.
Choice of reference genes
The results of NormFinder and GeNorm analyses corresponded well. The two tools have not always agreed on the particular order of the individual genes' expressions' stability values; however, the final choice of the best reference genes was almost uniform. The calculation of V-values by GeNorm for the proposed genes (Fig. 2B) is useful for deciding the optimal number to be used in an expression study; pairwise variation between samples is reduced by the inclusion of additional reference genes and therefore indicates the number of genes required to achieve an arbitrarily selected threshold of reference gene stability; a recommended cut-off value is 0.15. If we took into account the results of the pairwise variations, then four genes were necessary for wheat and oats, and three genes for barley. For barley, the three genes with the most stable expression across the samples as determined by GeNorm and NormFinder were GAPDH, 18S rRNA and TUBB. For wheat and oat, the combination of GAPDH, 18S rRNA, TUBB and TUBA (oat) or EIF4A (wheat) should provide safe normalisation. For oats and barley, the poorly performing genes were 28S rRNA, EIF4A and EF1A; these could not be recommended as reference genes for these species. The use of TUBA as a reference gene is not recommended as appropriate for wheat tissues, as with an M-value of 2.13 it appears to be regulated.
Virus titre comparison results
BYDV symptoms manifest in two oat varieties compared to a healthy control plant. On the left (A), is a healthy control; in the middle (B) the tolerant cv. Atego, and on the right (C), susceptible cv. Saul.
Titre of BYDV in two oat varieties as determined by three quantification approaches. (A): The relative BYDV titre levels in all oat samples determined by relative quantification using GAPDH, TUBB, TUBA and 18S rRNA (GAPDH group); 28S rRNA, EIF4A and EF1A (EIF4A group); and by absolute quantification with cloned standards (absolute). The absolute numbers were transformed into relative ratios for comparison. The boxes represent the lower and upper quartiles with medians; the whiskers illustrate the 10-90 percentiles of the samples. (B): The relative BYDV titre levels in the two oat cvs. Atego (tolerant) and Saul (susceptible) determined by relative quantification using GAPDH, TUBB, TUBA and 18S rRNA (GAPDH group); 28S rRNA, EIF4A and EF1A (EIF4A group); and by absolute quantification with cloned standards (absolute). The absolute numbers were transformed into relative ratios for comparison. No significant differences were recorded.
There were no significant differences between the viral titre in cv. Saul (susceptible) and cv. Atego (tolerant) according to relative quantification, independently of the reference genes used in the assays. However, the trends of the two reference genes groups were slightly different (Fig. 4B). The GAPDH, TUBB, TUBA and 18S rRNA assays resulted in almost identical data from the two groups; however, the 28S rRNA, EIF4A and EF1A results might lead to the assumption that the tolerant cv. Atego had a lower virus load than the susceptible variety. Furthermore, the BYDV titre was determined to be statistically higher in the tolerant compared to the susceptible variety, by absolute quantification. This seems quite unlikely, and contrary to the plant symptoms, indicating that the absolute quantification with no normalisation is inadequate to quantify viruses in plants. These results support the necessity of the correct choice of reference genes for valid experimental data.
A similar study was conducted by Balaji et al. [20], in which the titre of BYDV-PAV was determined in susceptible and resistant wheat lines, using 18S rRNA as a reference gene. In their study, the titre of the virus was measured 1, 2, 4, 6, 8, 10, 12, and 14 d post-infection (dpi); the greatest differences were in the early stages of infection (4 dpi), when the titre varied as much as tenfold, the most susceptible wheat line being that with higher titre. However, later in infection (≥ 12 dpi), the titre in the susceptible line decreased to almost the same as in the resistant line. As the oat samples in this study were collected at 30 dpi, when the symptom differences appeared, the results of the relative quantification seem to be accordant. Nevertheless, for the comparison of such subtle differences, missing or inappropriately chosen normalisation could easily lead to erroneous interpretation of data.
Conclusions
Our results indicate GAPDH, TUBB and 18S rRNA were the most stable reference genes for virus-infected cereals, giving very good statistical reliability according to the two software packages employed. Moreover, the use of four reference genes (GAPDH, 18S rRNA, TUBB and TUBA-oat/EIF4A-wheat) was sufficient for a reliable normalisation of the viral genes in wheat and oats, and three genes (GAPDH, TUBB and 18S rRNA) were necessary for barley. The absolute quantification with no reference gene is not recommended, since it might lead to inaccurate and misleading conclusions, due to lack of normalising factors. Studies aimed at the relative comparison of the BYDV titre among different cereal species also seem to be unfeasible because of different gene expression between species.
Methods
Plant, virus and vector material
Laboratory isolates of BYDV-PAV and BYDV-PAS maintained on barley seedlings and transmitted by aphids Rhopalosiphum padi and Sitobion avenae were used for the study. For each of wheat (Triticum aestivum L.), barley (Hordeum vulgare L.) and oats (Avena sativa L.), 60 plants were sown into individual pots. For each species, 30 plants were of the BYDV-susceptible and 30 of the BYDV-tolerant variety: wheat cvs. Anza (tolerant) and Vlada (susceptible); oat cvs. Atego (tolerant) and Saul (susceptible); and barley cvs. Wysor (tolerant) and Finesse (susceptible). They were maintained in greenhouses with controlled conditions of 16-h-light period and 20°C. When the plants were at the growth stage where two leaves were unfolded, 20 plants per species (10 per variety) were separated and infected with BYDV-PAS strain by vector transmission. For each plant, approximately five viruliferous R. padi aphids were transferred onto and allowed to feed on the plant for 2 d. After that, the plants were treated with insecticide and kept in the greenhouses for one more week. On the last day of the week period, 20 additional plants of each species were separated and abiotic stress conditions were created for them. The plants were moved into a chamber with 4°C constant temperature and darkness and were kept there for 24 h. Then, all 180 plants were harvested and the whole aboveground biomass of individual plants was used for further analysis.
For the case study, oat plants from a field experiment testing the resistance of different varieties to BYDV were used: cvs. Atego and Saul. Atego did not manifest many symptoms, while Saul symptomatically reacted as very sensitive to the infection. The leaf samples from 20 randomly chosen plants were collected one month after the inoculation with BYDV-PAV and BYDV-PAS. For each plant, the most symptomatic and the least symptomatic leaves were chosen.
RNA purification and cDNA synthesis
The whole plants were ground with a mortar and a pestle in liquid nitrogen, and 100 mg of the material was then used for RNA isolation. The RNA was isolated with a Spectrum™ Plant Total RNA Kit (Sigma Aldrich) according to the manufacturer's instructions. The concentration and purity of the isolated RNA was then spectrophotometrically measured. Of total RNA, 5 μg was diluted to a total volume of 25 μL and digested with DNase I (DNA-free™; Ambion) according to the manufacturer's instructions. cDNA was produced using the M-MLV Reverse Transcriptase (Promega) according to the manufacturer's recommendations for oligo(dT)20 or random hexamers primed cDNA-synthesis. cDNA synthesis was performed using 1 μg of RNA, at 40°C. Finally, cDNA was diluted 1:5 before use in qPCR.
Standards RNA preparation
A specific BYDV nucleotide sequence (294 bp) amplified by RT-PCR was inserted into the vector pGem-T (Promega) and cloned into E. coli JM-109. The plasmid was linearised at the RsaI site and used as a target in an in vitro transcription reaction performed with Megascript T7 kit (Ambion) followed by DNase I (Ambion). The amount of RNA was quantified by spectrophotometry (Nanophotometer, Implen). The μg of single stranded RNA was converted to ρmol using the average molecular weight of a ribonucleotide (340 Da) and the number of bases of the transcript (Nb). The following mathematical formula was applied: ρmol of ssRNA = μg (of ssRNA) × (106 ρg/1 μg) × (1 ρmol/340 ρg) × (1/Nb). Avogadro's constant was used to estimate the number of transcripts (6.023 × 1023 molecules/mol). Thereafter, tenfold serial dilutions of the transcripts were prepared.
Real-Time PCR conditions
Characteristics of gene specific real-time RT-PCR assays
Gene symbol | Gene name | Accession No. | Primer sequence (5'→3') | Amplicon size | PCR efficiency |
---|---|---|---|---|---|
GAPDH | Glyceraldehyde-3-phosphate dehydrogenase | AK251456 | Forward: TGTCCATGCCATGACTGCAA | 105 | 101% |
Reverse: CCAGTGCTGCTTGGAATGATG | |||||
TUBA | Alpha tubulin-2B | AK250165 | Forward: TTCGCCCGTGGTCATTACA | 113 | 100% |
Reverse: GCATTGAAGACAAGGAAGCCC | |||||
TUBB | Beta-tubulin | U76897 | Forward: CAAGGAGGTGGACGAGCAGATG | 84 | 97% |
Reverse: GACTTGACGTTGTTGGGGATCCA | |||||
ELF1A | Elongation factor-1 alpha | AF479046 | Forward: CAGTGCTGGACTGCCACA | 164 | 91% |
Reverse: CTCCACCACCATGGGCTT | |||||
EIF4A | Eukaryotic initiation factor 4a | EU850433 | Forward: TTGTGCTGGATGAAGCTGATG | 76 | 99% |
Reverse: ACACCAACAGCCACAGTTTGC | |||||
18S rRNA | 18S ribosomal RNA | M82356 | Forward: GTGACGGGTGACGGAGAATT | 151 | 85%, 99%* |
Reverse: GACACTAATGCGCCCGGTAT | |||||
28s rRNA | 28S ribosomal RNA | M82206 | Forward: CCTGATCTTCTGTGAAGGGTTCGA | 172 | 94% |
Reverse: GGTTCGATTAGTCTTTCGCCCCTA | |||||
BYDV cp | coat protein of BYDV | EF521849 | Reverse:TGTTGAGGAGTCTACCTATTTG | 294 | 99% |
Forward:GTTGAGTTTAAGTCACACGC |
Results analysis
The relative gene expression ratios were calculated by a mathematical model, which includes an efficiency correction for real-time PCR efficiency of the individual transcripts [8]. The amplification efficiency was established for each of the targets from serial dilutions of cereal leaves within range 0.80-1.0. The absolute viral gene's quantification values were transformed into relative values by simple proportion. Two-way ANOVA was performed on data, followed by the Bonferroni post-tests, where appropriate. Data in graphs or tables are presented as means with their standard errors of the means. P > 0.05 was considered as non-significant, P < 0.01 very significant and P < 0.001 extremely significant. Reference gene selective analysis was performed with the NormFinder [14] and GeNorm [12] tools.
Declarations
Acknowledgements
We are deeply grateful to Prof. Ulrich K. Melcher for his critical reading of the manuscript and valuable comments. This work was supported by grants from the Ministry of Agriculture, Czech Republic for the project Nos. MZE 0002700604 (50%) and QH81269 (50%).
Authors’ Affiliations
References
- Rastgou M, Khatabi B, Kvarnheden A, Izadpanah K: Relationships of Barley yellow dwarf virus-PAV and Cereal yellow dwarf virus-RPV from Iran with viruses of the family Luteoviridae. Eur J Plant Pathology. 2005, 113: 321-326. 10.1007/s10658-005-1231-y. 4View ArticleGoogle Scholar
- D'Arcy CJ: Symptomatology and host range of barley yellow dwarf. Barley yellow dwarf: 40 years of progress. Edited by: D'Arcy CJ, Burnett PA, APS Press, St. Paul MN, 2005. 9-28.Google Scholar
- Cooper JI, Jones AT: Responses of plants to viruses: proposals for the use of terms. Phytopathology. 1983, 73 (2): 127-128. 10.1094/Phyto-73-127.View ArticleGoogle Scholar
- Osmond CB: Photorespiration and photoinhibition: some implications for the energetics of photosynthesis. Biochim Biophys Acta. 1981, 639: 77-98.View ArticleGoogle Scholar
- Kassanis B: Effects of changing temperature on virus diseases. Adv Virus Res. 2003, IV: 221-241.Google Scholar
- Hull R: Matthew's plant virology, Academic Press, 2002:1056-FourthGoogle Scholar
- Chen BE, Kondo M, Garnier A, Watson FL, Puettmann-Holgado R, Lamar DR, Schmucker D: The molecular diversity of Dscam is functionally required for neuronal wiring specificity in Drosophila. Cell. 2006, 125 (3): 607-620. 10.1016/j.cell.2006.03.034.PubMedView ArticleGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.PubMedPubMed CentralView ArticleGoogle Scholar
- Cappelli K, Felicetti M, Capomaccio S, Spinsanti G, Silvestrelli M, Supplizi AV: Exercise induced stress in horses: selection of the most stable reference genes for quantitative RT-PCR normalization. BMC Mol Biol. 2008, 9: 49-10.1186/1471-2199-9-49.PubMedPubMed CentralView ArticleGoogle Scholar
- Bustin SA, Nolan T: Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction. J Biomol Tech. 2004, 15: 155-66.PubMedPubMed CentralGoogle Scholar
- Pfaffl MW, Windisch W: Influence of zinc deficiency on the mRNA expression of zinc transporters in adult rats. J Trace Elem Med Biol. 2003, 17: 97-106. 10.1016/S0946-672X(03)80005-6.PubMedView ArticleGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal reference genes. Genome Biology. 2002, 3: research0034.1-0034.11. 10.1186/gb-2002-3-7-research0034.View ArticleGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP: Determination of stable reference genes differentially regulated target genes and sample integrity: BestKeeper--Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26: 509-15. 10.1023/B:BILE.0000019559.84305.47.PubMedView ArticleGoogle Scholar
- Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization applied to bladder and colon cancer data sets. Cancer Res. 2004, 64: 5245-50. 10.1158/0008-5472.CAN-04-0496.PubMedView ArticleGoogle Scholar
- Brunner A, Yakovlev I, Straus S: Validating internal controls for quantitative plant gene expression studies. BMC Plant Biology. 2004, 4: 14-10.1186/1471-2229-4-14.PubMedPubMed CentralView ArticleGoogle Scholar
- Faccioli P, Ciceri GP, Provero P, Stanca AM, Morcia C, Terzi V: A combined strategy of "in silico" transcriptome analysis and web search engine optimization allows an agile identification of reference genes suitable for normalization in gene expression studies. Plant Mol Biol. 2007, 63: 679-688. 10.1007/s11103-006-9116-9.PubMedView ArticleGoogle Scholar
- Jain M, Nijhawan A, Tyagi AK, Khurana JP: Validation of reference genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochem Biophys Res Commun. 2006, 345: 646-651. 10.1016/j.bbrc.2006.04.140.PubMedView ArticleGoogle Scholar
- Nicot N, Hausman J-F, Hoffmann L, Evers D: Reference gene selection for RT-qPCRnormalization in potato during biotic and abiotic stress. J Exp Bot. 2005, 56: 2907-2914. 10.1093/jxb/eri285.PubMedView ArticleGoogle Scholar
- Xu WF, Shi WM: Expression Profiling of the 14-3-3 Gene Family in Response to Salt Stress and Potassium and Iron Deficiencies in Young Tomato (Solanum lycopersicum) Roots: Analysis by Real-time RT-PCR. Ann Bot. 2006, 98: 965-974. 10.1093/aob/mcl189.PubMedPubMed CentralView ArticleGoogle Scholar
- Balaji B, Bucholtz DB, Anderson JM: Barley yellow dwarf virus and Cereal yellow dwarf virus quantification by real-time polymerase chain reaction in resistant and susceptible plants. Phytopathology. 2003, 93: 1386-1392. 10.1094/PHYTO.2003.93.11.1386.PubMedView ArticleGoogle Scholar
Copyright
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.