Microarray-based uncovering reference genes for quantitative real time PCR in grapevine under abiotic stress
© Coito et al.; licensee BioMed Central Ltd. 2012
Received: 31 January 2012
Accepted: 18 April 2012
Published: 7 May 2012
Quantitative real time polymerase chain reaction is becoming the primary tool for detecting mRNA and transcription data analysis as it shows to have advantages over other more commonly used techniques. Nevertheless, it also presents a few shortcomings, with the most import being the need for data normalisation, usually with a reference gene. Therefore the choice of the reference gene(s) is of great importance for correct data analysis. Microarray data, when available, can be of great assistance when choosing reference genes. Grapevine was submitted to water stress and heat stress as well as a combination of both to test the stability of the possible reference genes.
Using the analysis of microarray data available for grapevine, six possible reference genes were selected for RT-qPCR validation: PADCP, ubiq, TIF, TIF-GTP, VH1-IK, aladin-related. Two additional genes that are commonly used as reference genes were included: act and L2. The stability of those genes was tested in leaves of grapevine in both field plants and in greenhouse plants under water or heat stress or a combination of both. Gene stability was analyzed with the softwares GeNorm, NormFinder and the ΔCq method resulting in several combinations of reference genes suitable for data normalisation. In order to assess the best combination, the reference genes were tested in putative stress marker genes (PCO, Galsynt, BKCoAS and HSP17) also chosen from the same microarray, in water stress, heat stress and the combination of both.
Each method selected different gene combinations (PADCP + act, TIF + TIF-GTP and ubiq + act). However, as none of the combinations diverged significantly from the others used to normalize the expression of the putative stress marker genes, then any combination is suitable for data normalisation under the conditions tested. Here we prove the accuracy of choosing grapevine reference genes for RT-qPCR through a microarray analysis.
KeywordsMicroarray Grapevine Real time RT-PCR Reference gene
Quantitative real time polymerase chain reaction (RT-qPCR) has become a mainstream research tool for the quantification of mRNA and transcription data analysis [1, 2]. This method presents many advantages when compared with the more commonly used reverse transcriptase polymerase chain reaction (RT-PCR), Northern blotting and microarrays. These advantages being its higher sensitivity, specificity, broad quantification and avoidance of post-PCR processing [2–4].
Approaches such as array-based transcription profiling technologies allow the assessment of expression levels of thousands of genes in control and stress tissues. Gene redundancy is, however, a common trait of this type of analysis and requires removal. RT-qPCR technology, due to the higher sensitivity and specificity [2, 4, 5] is excellent to confirm non-redundant gene expression obtained through microarray analysis. However RT-qPCR itself presents several shortcomings, the most important and relevant being the need for data normalisation [1, 6–9]. Normalisation and homogenization  are of the utmost importance in RT-qPCR as they allow the corrections of PCR reactions due to inaccurate quantification of RNA or problems related to RNA quality and purity [1, 11].
RT-qPCR data normalisation is usually performed using the expression of an internal control gene [4, 6, 12, 13]. Such a gene is also referred as a “reference gene” or sometimes a “housekeeping gene” i.e. a gene whose expression does not change under the different conditions or tissues under investigation [1, 14, 15]. Unfortunately no such gene exists, and the universality of such an ideal gene is not valid, since the transcript levels of all genes show some degree of variability under different experimental conditions [5–7, 16]. Initially, gene stability was assessed with the ΔCq method. In this method, the Cq value of the gene of interest (either target or reference) is related to a control/calibrator . In order to overcome these flaws, statistical algorithms such as GeNorm  and NormFinder  have been developed to evaluate the best suited reference gene or a combination of genes for normalisation of RT-qPCR data in a specific set of biological data .
The identification of suitable reference genes can be difficult. Several attempts have been made, all with different outcomes. Studies usually allocate different “model” genes to be used in data normalisation. Usually the allocated genes vary with the plant species, as well with the experimental conditions. Also the method to select the reference candidate varies with the availability of data for the plant species under study. Studies regarding reference genes often employ a variety of methods to chose possible reference genes; such as searching the bibliography for published references , using orthologs of Arabidopsis reference genes , cDNA libraries  or analysis of EST libraries . Genes that are frequently identified to be good references include elongation factors 1-α (eEF-1α) [23–26], actin (act) [8, 15, 21, 23], ubiquitin (ubiq) [8, 25, 27–29], glyceraldehyde 3-phosphate dehydrogenase (GAPDH) [7, 25, 27, 28], ribosomal proteins [7, 21, 25, 30–32], SAND family protein (SAND) [20, 27, 33] and, other less common genes have been identified to be particularly good references in very specific contexts.
Czechowski et al. , using data from Affymetrix ATH1 whole-genome GeneChip, proposed not only the typical reference genes for RT-qPCR but also new ones. After that study, microarrays, when available, have been used for identifying reference genes [34, 35], with new genes being found for RT-qPCR data normalisation.
In grapevine, several attempts have been made in identifying reference genes. Gamm et al.  indicated two genes (V-type proton ATPase 16 kDa proteolipid subunit and 60S ribosomal protein L18) as being optimal reference genes for the study of the expression of genes involved in pterostilbene synthesis in grapevine leaves infected by P. viticola and berries infected by B. cinerea. Reid et al.  also tested possible reference genes specifically suitable for use in grapevine berry development studies and suggested GAPDH, act, eEF-1α and SAND as the most stable.
Whatever the method used for choosing possible reference genes for data normalisation, stability analysis should always be performed in the optimal conditions.
Grapevine (Vitis vinifera L.) is a sessile organism and therefore cannot avoid abiotic stress. Plants have been developing mechanisms to cope with environmental changes and help to overcome them. The most relevant abiotic stresses that can affect the production of a Mediterranean crop such as grapevine are: drought, excessive light and excessive heat. In fact, climate models predict an intensification of extreme conditions, which can reduce production to below the threshold for optimal grapevine growth . This represents a serious challenge for Mediterranean agriculture.
The pattern of gene expression in response to abiotic stress has been monitored in ex vitro grapevine plants by comparing the use of the grapevine Affymetrix GeneChip with extensive RT-qPCR analysis ; and in greenhouse and field plants through genechip microarray (results under analysis).
In this paper we report several putative reference genes chosen from a grapevine microarray analysis and identify the genes to be used as references for RT-qPCR normalisation, after obtaining the “optimal combination” of reference genes using three different methods. We also apply the three “optimal combinations” obtained to quantify the expression of a set of stress-marker genes.
Results and discussion
Choice of reference genes
Possible reference genes retrieved from the microarray analysis
growth-on protein GRO10
F-box protein 7
RAB GTPase ARA3
DnaJ homolog, subfamily B, member 4
ubiquitin-like domain containing CTD phosphatase 1
peptidylprolyl isomerase PAS1 (PASTICCINO 1)
DNA polymerase eta subunit
translation initiation factor eIF-3 subunit 4
Protein kinase PKN/PRK1
translation initiation factor eIF-2B alpha subunit
RNA-binding protein Musashi
DNA repair protein RAD23
ankyrin repeat family protein
ribosomal protein L27
Plastid-specific 50S ribosomal protein 6
Primers used for the reference genes chosen after analysis of Table 1
5' ATTCATTAAAGTACCTTTCTTT 3'
5' AACACCCAAAAGATGTCGTA 3'
5' CAATTTCCTGAGTTCTACAGTT 3'
5' CCTCATTGTATGACTCCCAGT 3'
5' AAAGCAGAAGAAACCAAGATT 3'
5' TTGCCAGTGCCTGTAGTAGCC 3'
5' AGCAGCACAGAATAAGAAACT 3'
5' CCATCAGCCCCAACAAATACC 3'
5' CAGGGATTATGATAGTAGGA 3'
5' TTGTTTGGTAGAGGAGGTGG 3'
5' CCTACACTTATTCATCTTCG 3'
5' ACTTGTGGCGGTTGCTCTGC 3'
5´ TGGATTCTGATGGTGTGAGTC 3'
5' CAATTTCCCGTTCAGCAGTAGTGG 3'
5' TCTACTTCAACCGATATGC 3'
5' CAACCTGTCCGACTG 3'
Selection of the best reference genes
Ranking of reference genes using ΔCq and GeNorm
M values of Ref genes
Sum of Ref Gene values
The GeNorm application for Microsoft Excel determines the most stable reference genes from a group of genes. The application also calculates a normalisation factor of gene expression for each sample, based on the geometric mean of a user-defined number of reference genes. It is based on the assumption that the expression of two ideal reference genes will always have the same ratio among samples regardless of the experimental conditions . This average expression stability value (M) is calculated using the expression data for each gene. M is the average pairwise variation (V) of one gene compared with each of the other reference genes tested. Stepwise exclusion of the gene with the highest M value allows the ranking of the tested genes according to their expression stability, until the two most stable genes in the remaining set cannot be ranked any further. GeNorm also allows estimating the optimal number of reference genes which should be used for normalisation. It calculates the normalisation factor (NF) based on the geometric mean of the expression of more than one reference gene. GeNorm calculates the NFn of the two most stable reference genes based on the geometric mean of the expression data, and then the NFn+1 with the next most stable gene. To determine how many genes should be used for accurate normalisation, the pair-wise variation (Vn/n+1) is calculated for each two sequential normalisation factors (NFn) and NFn+1.
While GeNorm stepwise finds the two genes whose expression ratio shows the least variation in relation to the other genes; NormFinder finds the single gene with the most stable expression and the best pair of genes with the most stable combined expression . The NormFinder pair can compensate in a way that, for best performance, a gene that is slightly overexpressed in a treatment group, and slightly underexpressed in the untreated group, can be combined with a gene with the opposite bias. Furthermore, NormFinder can account for the heterogeneity in the tested samples, allowing the comparison of different treatment groups. It can thus distinguish between stability and bias  eventually being able to discard a candidate reference gene that is excellent for all treatment groups but one. This indicates that the particular treatment affects that reference gene, which is then not such a good candidate. These are the most striking differences between the two analysis softwares, and that they can account for different results when comparing them.
Ranking of reference genes using NormFinder
Greenhouse vs field
TIF + TIF-GTP
In all the GeNorm analyses L2 was the worst performing reference gene (Table 3 and Figure 1) while aladin-related was the worst performing reference gene in the ΔCq and NormFinder analyses (Table 3 and Figure 2).
Optimal number of genes for RT-qPCR data normalisation
Specific markers of abiotic stress
For the normalisation tests, we chose as best combinations of reference genes (BC) the following: GeNorm: act + ubiq; NormFinder: TIF + TIF-GTP; ΔCq: PADCP + act. L2 was also used for comparison as the most unstable gene.
Putative stress marker genes retrieved from the microarray analysis
Primers for the putative stress marker genes
5' GCGTCTCATTATCGTTGGTTC 3'
5' CAGTGTCCTCGTGGTATCG 3'
5' CCAATCCCTTCTGAATACAACC 3'
5' TTTCCCACCATTTCCTCACC 3'
5' TGCGACAAGGGCTTTCATC 3'
5' CAGGCTCCAGATCATACTCAG 3'
5' AGAAGAAGAGCCAGAAGAGAAG 3'
5' ACACACGAAGCGACCAAG 3'
GalSynt (Galactinol synthase) catalyses the first committed step in the biosynthesis of the rafinose oligosaccharide family (RFOs) and plays a key regulatory role in carbon partitioning between sucrose and RFOs . In Arabidopsis the GalSynt isoform AtGalSynt 1 is induced by drought and salinity . In the microarray analysis this experiment was based upon, GalSynt was the most up-regulated gene under WS, and was therefore chosen as an up-regulation marker of WS. Indeed GalSynt showed significant up-regulation with all the BCs and L2 under WS (Figure 5A); while in HS it was down-regulated when normalised with the TIF + TIF-GTP combination and showed a slight up-regulation when normalised with the remaining BCs and L2 (Figure 5B). The results obtained in HS were significantly lower than the expression levels verified in WS with all the BCs and L2. This suggested that not only this gene responds well to the treatment imposed, but that all the BCs are stress insensitive and therefore suitable for data normalisation.
BKCoAS (β-ketoacyl-CoA synthase) is one of the enzymes responsible for the elongation of fatty acid chains. Very-long-chain fatty acids (VLCFAs) are fatty acids that have 20 or more carbon atoms. In plants, VLCFAs are the precursors of several lipids, such as cuticular waxes, and their derivatives act as protection barriers, they also provide energy storage in the seeds and act as signalling molecules, responding to several stress stimuli . In this work, the down-regulation of BKCoAS was studied as a possible marker of HS. The expression of BKCoAS under HS was down-regulated when normalised with all BCs and L2, as expected (Figure 5B). However under WS BKCoAS also displayed down-regulation when normalised with all BCs and L2 (Figure 5A). Statistical analysis showed significant differences in the expression of BKCoAS between WS and HS (with the expression of BKCoAS being lower in HS) only when the gene expression was normalised with the TIF + TIF-GTP combination; making this combination ideal as a reference for RT-qPCR data normalisation in this case.
The accumulation of heat shock proteins (HSPs), under the control of heat stress transcription factors (HSFs), is known to play a central role in the response to heat stress and in acquired thermotolerance in plants and other organisms [44,. We have used the gene HSP17 as an up-regulation marker of HS. Under HS HSP17 was up-regulated with all the BCs and L2 (Figure 5B). Under WS HSP17 showed differential expression, being up-regulated when normalised with the TIF + TIF-GTP combination and down-regulated with the remaining BCs and L2 (Figure 5A). Statistical analysis showed that the expression of HSP17 was significantly higher in HS than in WS, when normalised with all BCs and L2, therefore making any BC and L2 suitable for RT-qPCR data normalisation.
Gene expression of the chosen markers was also tested under combined WS and HS conditions (Figure 5C) in order to evaluate the reliability of the genes chosen as stress markers in a complex environment when more than one factor is at stake. All the genes used as stress markers responded as expected, with PCO down-regulated in response to WS and BKCoAS was down-regulated due to HS. Similarly, GalSynt was significantly up-regulated in WS, whilst HSP17 was up-regulated in HS. When statistical analysis was performed between the expression values of the stress marker genes under WSHS and their expected individual stress response; we found that PCO expression, when normalised with the combinations act + ubiq, PADCP + act and L2, was significantly lower when compared with its expression under WS. In this case, the best combination for data normalisation would be TIF + TIF-GTP. Under WSHS, GalSynt expression showed significant differences in relation to WS when normalised with the three BCs. In this occasion, L2 would be the most suitable gene for data normalisation, probably due to the expression stability of GalSynt. Under WSHS no significant differences were found in BKCoAS when compared to the expression of this gene in HS; rendering all BCs suitable for data normalisation. Under WSHS HSP17 showed a significantly higher level of expression when compared to HS, only when L2 was used as reference, making it a desirable reference gene in such situation.
Most works propose different combinations of reference genes for different experimental conditions after combining the best reference genes that present a ranking consensus between different methods [15, 25, 27]. Conversely, our reference genes were selected to be used in analysing gene expression after different treatments. We used the BCs obtained and tested them individually by method. In fact, if we had used the four best ranking genes of all treatments and the three methods, we would have obtained the same four reference genes that comprise the three BCs. Therefore we propose that the BCs obtained with any of the methods are suitable reference genes for transcription studies and this option is obviously less time and resource consuming than the use of all four genes.
Regarding the most stable genes, ubiq is referred as a reference gene in scientific works [8, 21, 25, 27], as confirmed in the microarray analysis, although its presence was recorded in only one combination. act was chosen due to its common use as reference gene [8, 21, 23] and was present in two of the three combinations. Translation initiation factors have been tested by other authors but without success [17, 19], while in our study TIF and TIF-GTP showed good stability and their combination was selected by NormFinder. PADCP is, to our knowledge, referred to for the first time as a potential reference gene and not only showed stable expression in the microarray analysis, but was also selected as reference gene by the GeNorm software. Ribosomal protein genes have been used as references in several studies [7, 19, 31] while in others they were ruled out as not stable enough [24, 29]. L2, as one of those genes, was considered a good reference in one study  but did not perform so well in another . In our work L2 ranked as the least stable gene; however in some situations when normalizing stress markers, it outperformed the BCs, highlighting its major problem of being inconsistent.
This study attempts to provide the validation of reference genes in grapevine leaves under abiotic stress conditions for RT-qPCR data normalisation. The availability of a custom microarray for grapevine was of great assistance in the choice of candidates. We used plants submitted to water and heat stresses, as well as a combination of both. Six genes that presented high levels of stability in the microarray were tested in RT-qPCR, as well as two other genes commonly used in RT-qPCR data normalisation. Through the use of GeNorm, NormFinder and ΔCq methods, we obtained three possible reference gene combinations: act + ubiq; TIF + TIF-GTP and PADCP + act. With regards in obtaining the best combination of genes resulting from the analysis with the different softwares and methods, it was not possible to appoint a single optimal combination lining off from the others. All combinations of reference genes were able to normalise the putative stress markers; but not without flaw. Therefore we can conclude that any of the combinations tested is suitable to be used as reference for RT-qPCR data normalisation of grapevine leaf samples under abiotic stress; we prove that microarray analysis can be a powerful tool to obtain reference genes. When microarrays are not available some of the most commonly used references (act and ubiq) are in fact good options, but this choice must be taken with caution because some are not (e.g. L2). The careful testing of reference genes also comes out as paramount when compared to the test method. In fact, the three alternatives evaluated here performed equally well.
Greenhouse plant material and stress treatments
Cuttings from pruned wood of pre-selected plants of the variety Trincadeira were grown in pots in the greenhouse under the following controlled conditions: 200 μmol quanta m-2 s-1 irradiance, 16 h light/8 h dark photoperiod, temperature of 25°C day/ 23°C night and well watered with nutrient solution whenever necessary.
Individual stresses were applied when shoots were 50 to 60 cm high. The stresses applied were: HS – 1 hour at 42°C; WS – stop irrigation until the pre dawn leaf water potential (Ψw) was −0.9 MPa and WSHS – a combination of both. Ψw was measured with a pressure chamber, Model 600, PMS Instruments Company (Albany, OR). Samples consisted of the third, fourth and fifth totally expanded leaves and were taken immediately after the end of the stress (or, in the case of WS, after the measurement of the pre dawn leaf Ψw), frozen in liquid nitrogen and stored at −80°C until RNA extraction.
Total RNA extraction and cDNA synthesis
Samples were ground with a mortar and pestle in the presence of liquid nitrogen. Total RNA was extracted with the RNA Plant Total RNA Kit (Sigma-Aldrich, Inc) following the manufacturer’s instructions. Nucleic acid concentration of each sample was quantified by spectrophotometry using the software Gen5 1.09 (Synergy HT, Bio-Tek Instruments, Winooski, USA). Total RNA quality was assessed using the A260/A280 and A260/A230. Only RNA samples with A260/A280 between 1.8 and 2.1 and A260/A230 between 2.0 and 2.2 were used. Total RNA integrity was checked through 1% agarose gel electrophoresis under denaturing conditions.
RNA samples were treated with RQ1 RNase-Free DNase (Promega, Madison, WI). cDNA was synthesized from 2 μg of total RNA using oligo(dT)20 in a 20 μL-reaction volume using RevertAid Reverse Transcriptase (Fermentas Life Science, Helsingborg, Sweden) according to the manufacturer’s recommendations. cDNA was tested for gDNA contamination in PCRs using the intron spanning primers ubiq (Table 2) that yield a 229 bp amplicon in cDNA and a 547 amplicon in gDNA. Amplicon sizes were compared in 2% agarose gels together with the molecular weight marker 1Kb + (Invitrogen) and no gDNA contamination was detected. cDNA was stored at −20°C until further use.
Selection of reference genes and primer design
The selection of the possible reference genes for RT-qPCR was made based on a previous microarray analysis performed with two biological replicates for each experimental condition (data not shown) using a 23 096 unigene sequences array  based on the lower fold-change, function, category and presence in all replicates. The fold-change chosen was within the interval of −1.25 to 1.25, which is the smallest interval where sufficient possible reference genes could still be detected. The possible reference genes were selected from this pool, taking into account their category and function in order to obtain a mix of genes. This included members of families commonly used as reference genes in RT-qPCR data normalisation and other putative genes which do not have a clearly described function or which function remains unknown in grapevine but nevertheless displayed high stability on the microarray analysis. Six possible reference genes were selected from this group and two typically used reference genes were added to the study for comparison, act (AF369525.1) and L2 (AJ441290.2).
Primers for these eight putative reference genes were designed using the software Primer Premier 5.0 (Premier Biosoft International) using a primer length of 20 ± 2 bp, melting temperature of 60°C ± 2°C, a guanine-cytosine content of circa 50% and an expected amplicon size of 180–280 bp.
The real-time PCR was performed in 96 well white reaction plates (Bio-Rad, Hercules, CA), using an IQ5 Real Time PCR (Bio-Rad, Hercules, CA) with three biological replicates and two technical replicates. The 20 μL reaction mixture was composed of 1 μL cDNA diluted 50-fold, 0.5 μM of each gene-specific primer and 10 μL master mix (SsoFast_EvaGreen Supermix, Bio-Rad, Hercules, CA). Amplification of PCR products was monitored via intercalation of Eva-Green (included in the master mix). The following program was applied: initial polymerase activation, 95°C, 3 min; then 40 cycles at 94°C 10 s (denaturation), 60°C 20 s (annealing), 72°C 15 s (extension). The PCR products were run on 2% agarose gels to make sure that there was only one amplicon of the expected size. PCRs with each primer pair were also performed on samples lacking cDNA template, in triplicate (no template controls).
To assess amplification efficiency of the candidate genes, identical volumes of cDNA samples were diluted and used to generate five-point standard curves based on a five-fold dilution series (1;1:5;1:25;1:125;1:625), in triplicate. Amplification efficiency (E) is calculated as E = 10(−1/a)-1, “a” being the slope of the linear regression curve (y = a log(x) + b) fitted over the log-transformed data of the input cDNA dilution (y) plotted against the respective quantification cycle (Cq) values (x). E-values of the target genes were considered comparable when they did not exceed 100 ± 10%, corresponding to a standard curve slope of 3.3 ± 0.33. All cDNA samples were diluted 50 fold and were amplified in duplicate in two independent PCR runs.
To generate a baseline-subtracted plot of the logarithmic increase in fluorescence signal (ΔRn) versus cycle number, baseline data were collected between the cycles 5 and 17. All amplification plots were analysed with an R n threshold of 0.2, at the beginning of the region of exponential amplification, to obtain Cq (quantification cycle) and the data obtained were exported into a MS Excel workbook (Microsoft Inc.) for further analysis.
For the relation between the expressions of the different marker genes with the different best combination genes the relative quantity values were transformed into log2 (thus rendering them parametric) and tested through ANOVA in the program SAS 9 for Windows, SAS Institute Inc., Cary, NC, USA. When the p value of the ANOVA was lower than 0.05 a Tukey test was performed and statistically significant differences were accepted for a p value lower than 0.05.
The authors thank Prof. Leonor Morais for critical reading of the manuscript and Frederic Cobb for revising the English writing. The research was funded by Fundação para a Ciência e Tecnologia (FCT): project PTDC/AGR-GPL/099624/2008 that includes a reaserch grant to JLC; CBAA (PestOE/AGR/UI0240/2011), and the post-doc grant SFRH/BPD/64905/2009 to MR.
- Bustin SA, Benes V, Nolan T, Pfaffl MW: Quantitative real-time RT-PCR–a perspective. J Mol Endocrinol. 2005, 34: 597-601. 10.1677/jme.1.01755.PubMedView Article
- Gachon C, Mingam A, Charrier B: Real-time PCR: what relevance to plant studies?. J Exp Bot. 2004, 55: 1445-1454. 10.1093/jxb/erh181.PubMedView Article
- Ginzinger DG: Gene quantification using real-time quantitative PCR: an emerging technology hits the mainstream. Exp Hematol. 2002, 30: 503-512. 10.1016/S0301-472X(02)00806-8.PubMedView Article
- Guénin S, Mauriat M, Pelloux J, Van Wuytswinkel O, Bellini C, Gutierrez L: Normalization of qRT-PCR data: the necessity of adopting a systematic, experimental conditions-specific, validation of references. J Exp Bot. 2009, 60: 487-493. 10.1093/jxb/ern305.PubMedView Article
- Czechowski T, Stitt M, Altmann T, Udvardi MK: Genome-Wide Identification and Testing of Superior Reference Genes for Transcript Normalization. Plant Physiol. 2005, 139: 5-17. 10.1104/pp.105.063743.PubMedPubMed CentralView Article
- 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 control genes. Genome Biol. 2002, 3: research0034.1-research0034.11. 10.1186/gb-2002-3-7-research0034.View Article
- Barsalobres-Cavallari CF, Severino FE, Maluf MP, Maia IG: Identification of suitable internal control genes for expression studies in Coffea arabica under different experimental conditions. BMC Mol Biol. 2009, 10: 1-10.1186/1471-2199-10-1.PubMedPubMed CentralView Article
- Artico S, Nardeli SM, Brilhante O, Grossi-de-Sa MF, Alves-Ferreira M: Identification and evaluation of new reference genes in Gossypium hirsutum for accurate normalization of real-time quantitative RT-PCR data. BMC Plant Biol. 2010, 10: 49-10.1186/1471-2229-10-49.PubMedPubMed CentralView Article
- Derveaux S, Vandesompele J, Hellemans J: How to do successful gene expression analysis using real-time PCR. Methods. 2010, 50: 227-230. 10.1016/j.ymeth.2009.11.001.PubMedView Article
- Bustin SA, Beaulieu J-F, Huggett J, Jaggi R, Kibenge FSB, Olsvik PA, Penning LC, Toegel S: MIQE précis: Practical implementation of minimum standard guidelines for fluorescence-based quantitative real-time PCR experiments. BMC Mol Biol. 2010, 11: 74-10.1186/1471-2199-11-74.PubMedPubMed CentralView Article
- Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6: 279-284. 10.1038/sj.gene.6364190.PubMedView Article
- Bustin SA, Nolan T: Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction. J Biomol Tech. 2004, 15: 155-166.PubMedPubMed Central
- Dheda K, Huggett JF, Chang JS, Kim LU, Bustin SA, Johnson MA, Rook GAW, Zumla A: The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Anal Biochem. 2005, 344: 141-143. 10.1016/j.ab.2005.05.022.PubMedView Article
- Olsvik PA, Lie KK, Jordal A-EO, Nilsen TO, Hordvik I: Evaluation of potential reference genes in real-time RT-PCR studies of Atlantic salmon. BMC Mol Biol. 2005, 6: 21-10.1186/1471-2199-6-21.PubMedPubMed CentralView Article
- Hu R, Fan C, Li H, Zhang Q, Fu Y-F: Evaluation of putative reference genes for gene expression normalization in soybean by quantitative real-time RT-PCR. BMC Mol Biol. 2009, 10: 93-10.1186/1471-2199-10-93.PubMedPubMed CentralView Article
- Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, Grisar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits. J Biotechnol. 1999, 75: 291-295. 10.1016/S0168-1656(99)00163-7.PubMedView Article
- Brunner AM, Yakovlev IA, Strauss SH: Validating internal controls for quantitative plant gene expression studies. BMC Plant Biol. 2004, 4: 14-10.1186/1471-2229-4-14.PubMedPubMed CentralView Article
- Andersen CL, Jensen JL, Ørntoft 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-5250. 10.1158/0008-5472.CAN-04-0496.PubMedView Article
- Condori J, Nopo-Olazabal C, Medrano G, Medina-Bolivar F: Selection of reference genes for qPCR in hairy root cultures of peanut. BMC Res Notes. 2011, 4: 392-10.1186/1756-0500-4-392.PubMedPubMed CentralView Article
- Demidenko NV, Logacheva MD, Penin AA: Selection and Validation of Reference Genes for Quantitative Real-Time PCR in Buckwheat (Fagopyrum esculentum) Based on Transcriptome Sequence Data. PLoS One. 2011, 6: e19434-10.1371/journal.pone.0019434.PubMedPubMed CentralView Article
- Pellino M, Sharbel TF, Mau M, Amiteye S, Corral JM: Selection of reference genes for quantitative real-time PCR expression studies of microdissected reproductive tissues in apomictic and sexual Boechera. BMC Res Notes. 2011, 4: 303-10.1186/1756-0500-4-303.PubMedPubMed CentralView Article
- Paolacci AR, Tanzarella OA, Porceddu E, Ciaffi M: Identification and validation of reference genes for quantitative RT-PCR normalization in wheat. BMC Mol Biol. 2009, 10: 11-10.1186/1471-2199-10-11.PubMedPubMed CentralView Article
- Reid KE, Olsson N, Schlosser J, Peng F, Lund ST: An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biol. 2006, 6: 27-10.1186/1471-2229-6-27.PubMedPubMed CentralView Article
- Obrero A, Die JV, Román B, Gómez P, Nadal S, González-Verdejo CI: Selection of reference genes for gene expression studies in zucchini (Cucurbita pepo) using qPCR. J Agric Food Chem. 2011, 59: 5402-5411. 10.1021/jf200689r.PubMedView Article
- Martin RC, Hollenbeck VG, Dombrowski JE: Evaluation of Reference Genes for Quantitative RT-PCR in Lolium perenne. Crop Sci. 1881, 2008: 48-
- Nicot N, Hausman J-F, Hoffmann L, Evers D: Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. J Exp Bot. 2005, 56: 2907-2914. 10.1093/jxb/eri285.PubMedView Article
- Selim M, Legay S, Berkelmann-Löhnertz B, Langen G, Kogel K-H, Evers D: Identification of suitable reference genes for real-time RT-PCR normalization in the grapevine-downy mildew pathosystem. Plant Cell Rep. 2011, 31: 205-216.PubMedView Article
- Wan H, Yuan W, Ruan M, Ye Q, Wang R, Li Z, Zhou G, Yao Z, Zhao J, Liu S, Yang Y: Identification of reference genes for reverse transcription quantitative real-time PCR normalization in pepper (Capsicum annuum L.). Biochem Biophys Res Commun. 2011, 416: 24-30. 10.1016/j.bbrc.2011.10.105.PubMedView Article
- Tong Z, Gao Z, Wang F, Zhou J, Zhang Z: Selection of reliable reference genes for gene expression studies in peach using real-time PCR. BMC Mol Biol. 2009, 10: 71-10.1186/1471-2199-10-71.PubMedPubMed CentralView Article
- Jain M, Nijhawan A, Tyagi AK, Khurana JP: Validation of housekeeping 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 Article
- Kim B-R, Nam H-Y, Kim S-U, Kim S-I, Chang Y-J: Normalization of reverse transcription quantitative-PCR with housekeeping genes in rice. Biotechnol Lett. 2003, 25: 1869-1872.PubMedView Article
- Yan J, Yuan F, Long G, Qin L, Deng Z: Selection of reference genes for quantitative real-time RT-PCR analysis in citrus. Mol Biol Rep. 2011, 39: 1831-1838.PubMedView Article
- Remans T, Smeets K, Opdenakker K, Mathijsen D, Vangronsveld J, Cuypers A: Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations. Planta. 2008, 227: 1343-1349. 10.1007/s00425-008-0706-4.PubMedView Article
- Libault M, Thibivilliers S, Bilgin DD, Radwan O, Benitez M, Clough SJ, Stacey G: Identification of Four Soybean Reference Genes for Gene Expression Normalization. Plant Genome J. 2008, 1: 44-10.3835/plantgenome2008.02.0091.View Article
- Hruz T, Wyss M, Docquier M, Pfaffl MW, Masanetz S, Borghi L, Verbrugghe P, Kalaydjieva L, Bleuler S, Laule O, Descombes P, Gruissem W, Zimmermann P: RefGenes: identification of reliable and condition specific reference genes for RT-qPCR data normalization. BMC Genomics. 2011, 12: 156-10.1186/1471-2164-12-156.PubMedPubMed CentralView Article
- Gamm M, Héloir M-C, Kelloniemi J, Poinssot B, Wendehenne D, Adrian M: Identification of reference genes suitable for qRT-PCR in grapevine and application for the study of the expression of genes involved in pterostilbene synthesis. Mol Genet Genomics. 2011, 285: 273-285. 10.1007/s00438-011-0607-2.PubMedView Article
- Jones GV, White MA, Cooper OR, Storchmann K: Climate Change and Global Wine Quality. Clim Chang. 2005, 73: 319-343. 10.1007/s10584-005-4704-2.View Article
- Carvalho LC, Vilela BJ, Mullineaux PM, Amâncio S: Comparative Transcriptomic Profiling of Vitis vinifera Under High Light Using a Custom-Made Array and the Affymetrix GeneChip. Mol Plant. 2011, 4: 1038-1051. 10.1093/mp/ssr027. In pressPubMedView Article
- Pontin MA, Piccoli PN, Francisco R, Bottini R, Martinez-Zapater JM, Lijavetzky D: Transcriptome changes in grapevine (Vitis vinifera L.) cv. Malbec leaves induced by ultraviolet-B radiation. BMC Plant Biol. 2010, 10: 224-10.1186/1471-2229-10-224.PubMedPubMed CentralView Article
- Kumar Tewari A, Charan Tripathy B: Temperature-stress-induced impairment of chlorophyll biosynthetic reactions in cucumber and wheat. Plant Physiol. 1998, 117: 851-858. 10.1104/pp.117.3.851.PubMedPubMed CentralView Article
- Saravitz DM, Pharr DM, Carter TE: Galactinol synthase activity and soluble sugars in developing seeds of four soybean genotypes. Plant Physiol. 1987, 83: 185-189. 10.1104/pp.83.1.185.PubMedPubMed CentralView Article
- Taji T, Ohsumi C, Iuchi S, Seki M, Kasuga M, Kobayashi M, Yamaguchi-Shinozaki K, Shinozaki K: Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J. 2002, 29: 417-426. 10.1046/j.0960-7412.2001.01227.x.PubMedView Article
- Lee S-B, Jung S-J, Go Y-S, Kim H-U, Kim J-K, Cho H-J, Park OK, Suh M-C: Two Arabidopsis 3-ketoacyl CoA synthase genes, KCS20 and KCS2/DAISY, are functionally redundant in cuticular wax and root suberin biosynthesis, but differentially controlled by osmotic stress. Plant J. 2009, 60: 462-475. 10.1111/j.1365-313X.2009.03973.x.PubMedView Article
- Kim DH, Xu Z-Y, Na YJ, Yoo Y-J, Lee J, Sohn E-J, Hwang I: Small heat shock protein Hsp17.8 functions as an AKR2A cofactor in the targeting of chloroplast outer membrane proteins in Arabidopsis. Plant Physiol. 2011, 157: 132-146. 10.1104/pp.111.178681.PubMedPubMed CentralView Article
- Queitsch C, Hong SW, Vierling E, Lindquist S: Heat shock protein 101 plays a crucial role in thermotolerance in Arabidopsis. Plant Cell. 2000, 12: 479-492.PubMedPubMed CentralView Article
- Løvdal T, Lillo C: Reference gene selection for quantitative real-time PCR normalization in tomato subjected to nitrogen, cold, and light stress. Anal Biochem. 2009, 387: 238-242. 10.1016/j.ab.2009.01.024.PubMedView Article
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.