Skip to main content

Selection and validation of reference genes for RT-qPCR analysis in potato under abiotic stress



Real-time quantitative PCR (RT-qPCR) is the most commonly used method for accurately detecting gene expression patterns. As part of RT-qPCR analysis, normalization of the data requires internal control gene(s) that display uniform expression under different biological conditions. However, no invariable internal control gene exists, and therefore more than one reference gene is needed to normalize RT-qPCR results. Identification of stable reference genes in potato will improve assay accuracy for selecting stress-tolerance genes and identifying molecular mechanisms conferring stress tolerance in this species.


In the experiment, we assessed the expression of eight candidate internal control genes, namely elongation factor-1alpha (EF1α), actin, tubulin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), adenine phosphoribosyl transferase (APRT), 60S ribosomal protein L8 (L8), Cullin 3A (CUL3A), and exocyst complex component sec3 (sec3), in a diverse set of potato samples representing drought stress and osmotic stress challenges, and using geNorm, NormFinder, BestKeeper and RefFinder softwares.


The results indicated that EF1α and sec3 were the most stably expressed genes in the potato under drought and osmotic stress conditions. This work will facilitate future work on gene expression studies in potato and also benefit other species of the Solanaceae, such as tomato.


Drought and osmotic stress can seriously limit plant growth and productivity. Plants have developed multiple strategies to cope with drought and osmotic stress. These normally involve a mixture of stress avoidance and tolerance adaptations, which produce a range of changes at the morphological, physiological, cellular, and molecular levels [1]. Modern potato (Solanum tuberosum L.) varieties are considered to be sensitive to drought, but they have different morphological and physiological responses to water deficit compared to other crops [2]. Drought in field conditions results in significant losses in the yield and/or quality of potato tubers. With the increased global warming, there is a need to identify genotypes of the potato showing high tolerance to drought stress [3, 4].

Real-time quantitative PCR (RT-qPCR) is the most commonly used method for accurately detecting gene expression patterns [5]. It is necessary to use the normalization method for the variation of the control sample, caused by changes in RNA samples, reverse transcription efficiency and the quality of PCR efficiency. The use of one or more stable reference genes is the most common method for normalization of RT-qPCR data. As part of RT-qPCR analysis, normalization of the data requires internal control gene(s) that display uniform expression under different biological conditions [6]. However, no invariable internal control gene exists, and therefore more than one reference gene is needed to normalize RT-qPCR results [7]. Identification of stable reference genes in potato will improve assay accuracy for selecting stress-tolerance genes and identifying molecular mechanisms conferring stress tolerance in this species.

A reference gene is usually a housekeeping gene, ubiquitously expressed in all cells, and its product is essential for cell structure or metabolism, such as actin, ribosomal protein, EF1-α, and GAPDH [8]. However, these traditional reference genes are not always stably expressed under different circumstances. Therefore, it is necessary to select corresponding reference genes that are expressed at a constant level in certain cases. In potato, EF1-α was confirmed to be the most suitable reference gene under salt stress and late blight, and EF1-α and APRT were considered to be the most stable under cold stress [9, 10]. In addition, sec3, CUL3A, and C2 were found to be most suitable for screening of reference genes in edible tubers [11]. However, the stability of these reference genes under other abiotic stresses has not been confirmed. To the best of our knowledge, there is no previous report on the selection of suitable reference genes for potato under drought and osmotic stress. Therefore, the selection of stable reference genes for potato is helpful for future molecular studies using RT-qPCR.

RNA-seq is a powerful technique that can be used to provide estimates of gene and/or transcript expression, and RPKM (Reads Per Kilobase per Million reads) is widely used to represent the relative abundance of mRNAs for a gene [12, 13]. In our previous study, we used RNA-sequenced potato samples to gain insight into the molecular basis of drought adaptation by comparing digital expression profiles between drought treatment and the control treatment [14]. To gain insight into the transcriptome dynamics that are associated with drought stress, genome-wide gene expression profile was conducted by Solexa sequencing to generate a large dataset and a comprehensive transcriptome profiling for potato. Finally, we identified a number of genes that were stably expressed under drought stress, including many recognized housekeeping genes. An effective method for selecting a stable expression of a reference gene candidate from the RNA-Seq data can be done by using the coefficient of variation (CV) [15]. However, the transcriptome with corresponding treatment is inherently limited; it cannot be used to determine a suitable reference gene(s) for other experimental conditions. Therefore, for a systematic selection of reference genes, RT-qPCR is still the primary approach [16].

In the present study, we selected eight stable genes of different metabolic pathways according to the previous RNAseq analysis and screened them as reference genes under drought and osmotic stress. The RT-qPCR data generated were analyzed using the three most widely used algorithms, namely geNorm, NormFinder and Bestkeeper, to determine sets of reference genes suitable for gene expression studies in different experimental conditions [17]. Additionally, a comprehensive reference gene stability analysis tool RefFinder was used to confirm the ranking results obtained from geNorm, NormFinder and Bestkeeper. This work will facilitate future work on gene expression studies in potato and also benefit other species of the Solanaceae, such as tomato.


Plant material, growth conditions and stress treatments

The experiments were carried out with the potato (Solanum tubersosum L.) tetraploid cultivar ‘Atlantic’. Potatoes were grown in pots of Gansu Agricultural University (Lanzhou, China) greenhouse facilities, maintained at 18–23 °C and 70% relative humidity under natural light conditions. Three replicate pots per treatment were arranged in a randomized block design. Drought stresses were applied to the 6-week-old plants by stopping irrigation in the treatment plots; the control plots were irrigated continuously. Fresh plants were collected from drought-treated potato plants every 5 days; plants were collected 6 times.

To prevent microbial contamination, potato shoots were removed from tubers, surface sterilized (ethanol, HgCl2 and sodium hypochlorite) and placed in 30% (w/v) sucrose Murashige-Skoog (MS) medium. The shoots were incubated at 20 °C for 20 days. The stem sections of the test-tube seedlings were cut and transplanted into a special culture tube with an opening at the bottom containing the new MS medium. When the plants were 6 cm long, replacement of MS medium containing 20% polyethylene glycol 6000 (C6 M), while the other three plants were grown under MS medium. Fresh leaves were collected after 0, 2, 4, 8, 16 and 32 h.

We also carry out simulated drought treatment as follows: some of the test tube seedlings were planted on quartz sand containing MS medium, and when the plants grew to 6 cm, the MS medium was removed from the bottom of the special culture tube.

Fresh leaves were collected at 0, 2, 4, 8, 12 and 32 days. All treated plant materials were immediately frozen in liquid nitrogen and stored at − 80 °C until further RNA extraction and target gene expression analysis.

Potential reference gene selection

We selected 8 (CUL3A, EF1α, GAPDH, sec3, tubulin, L8, APRT and actin) commonly used reference genes for RT-qPCR analysis based on our previous publication [14] and our unpublished second RNAseq data set; these genes contain more than one exon. The sequences of these eight potato reference genes were obtained from the GenBank database and from the Potato Genomics Resource ( Primer pairs were designed from these sequences with the NCBI Primer-BLAST program (Table 1) ( [18]. Primers were designed across exon boundaries to avoid genomic DNA contamination and exon analysis using Splign ( [11]. Before RT-qPCR analysis, PCR was performed using primers as shown in Table 1 to determine the size specificity of the amplicons, then electrophoresed on ethoxylated gels and ethidium bromide, and the target amplicons were sequenced to confirm the identity of the PCR product.

Table 1 Candidate reference genes and primer sequences

RNA extraction and first strand cDNA synthesis

The sampled plants were grounded to fine powder with mortar and pestle in liquid nitrogen, and 100 mg of the material was used for RNA isolation. Total RNA was extracted using Trizol TIANGEG (TIANGEN, Beijing, China), according to the manufacturer’s instructions. Isolated RNA was treated with RNase-free DNase (TaKaRa, Japan) to remove genomic DNA contamination. Purity and concentration of RNA samples were measured using a micro-volume UV spectrophotometer (Quawell Q5000, Quawell, USA) and integrity was checked on agarose gel electrophoresis. RNA samples with 260/280 ratio between 1.9 and 2.1 were used for subsequent experimentation. First strand cDNA synthesis using cDNA synthesis kit (Sangon, Shanghai, China) was performed according to manufacturer’s instructions in a total volume of 20 μl containing 2 μg total RNA. RNA extraction and cDNA synthesis from all samples were performed for three biological replicates. The cDNA solution was 10 times diluted with nuclease-free water and aliquots were stored at − 20 °C until use in RT-qPCR [19].

Real-time quantitative PCR analysis

RT-qPCR analysis was carried out in 96-well plates with a Mx3005 multiplex quantitative RT-PCR system (Stratagene, USA) using the EvaGreen-based PCR assay. Each reaction (final volume of 20 µl) contained 1.0 μl template, 10.0 μl 2 ×  SuperReal Premix (TIANGEN, Beijing, China), 0.4  μl 50× ROX Reference Dye, 0.6  μl each primer (10 μmol/μl) and 6.2  μl double distilled water was carried out. The RT-qPCR conditions were pre-denaturation at 94  °C for 40 s, followed by 38 cycles of denaturation at 94  °C for 10 s, annealing at 57.2  °C for 30 s and with a final extension step at 72  °C for 10 min. Three technical replicates were set for each cDNA [20].

Gene expression stability analysis

A standard curve, repeated in three independent plates using a tenfold serial dilution of the mixed cDNAs was obtained from all tested samples as templates. The correlation coefficients (R2) and slope values were acquired from the standard curve. Then, we calculated the gene-specific PCR amplification efficiency of each gene. The corresponding real-time PCR efficiencies were calculated according to the equation: \({\text{E}} = \left[ {10^{{\left( {{ - }1/slope} \right)}} - 1} \right] \times 100\) [21].

Simultaneously, the amplicon characteristics, including Tm, length, amplification efficiency with standard deviation, and correlation coefficient, of the eight candidate reference genes are listed in Table 1.

To compare stability of expression among the candidate reference genes, the computational methods, geNorm [22], Normfinder [23], and BestKeeper [24] were applied to quantification cycle (Cq) for each gene’s expression data. These tools are based on different models and assumptions and each produced different results for the same gene’s expression data [25]. RefFinder was used to calculate a recommended comprehensive ranking based on the results of computational analysis, which in turn allowed us to identify the best reference genes for RT-qPCR data normalization in potato samples [26].

For geNorm and NormFinder analysis, the raw Cq values under different experimental designs were transformed into relative quantities using the formula 2−ΔCq (ΔCq = each corresponding Cq value-lowest Cq value) and then imported to geNorm to analyze gene expression stability value (M1). Similar to geNorm, NormFinder was further used to investigate the expression stability values (M2) for each gene and the pairwise variation of that gene with other reference genes. The reference gene with the highest M (M1 or M2) value is considered as the most unstable gene while the lowest M (M1 or M2) value indicated the most stable gene [22]. BestKeeper analysis was based on the untransformed Cq values and using coefficient of variance (CV) and the standard deviation (SD) of the Cq to evaluate the stability of reference genes. All three of the software programs were run based on the software manuals to select suitable reference genes [24]. By the combination of the three kinds of RefFinder ( software, we could easily rank the expression stability of reference genes in different experimental sets [27].


Assessment of primer specificity and PCR amplification efficiency

The sequences of CUL3A, EF1α, GAPDH, sec3, tubulin, L8, APRT and actin in potato were targeted by using specific primers on a template cDNA. The specificity of the designed primers was identified by gel electrophoresis and target amplicons were sequenced. The results showed a single band with the expected size by gel electrophoresis (Fig. 1), and the sequencing results of the target amplicon are also consistent with those of the target amplicons. The amplification efficiencies and correlation coefficients (R2) of 8 candidate reference genes in potato were calculated by slopes of the standard curves. The RT-qPCR amplification efficiencies for the 8 reference genes ranged from 95.56 to 103.40, and correlation coefficients ranged from 0.995 to 0.999 (Table 1). Thus, these primers can be used in RT-qPCR analysis.

Fig. 1
figure 1

Specificity of PCR and amplicon length. Amplified fragments shown by agarose gel electrophoresis with ethidium bromide staining. M marker 2000. Lanes 1, 2, 3, 4, 5, 6, 7 and 8 were the genes of EF1α, tubulin, GAPDH, sec3, CUL3A, L8, APRT and actin from potato

Cq Values of candidate reference genes

The quantification cycle (Cq) values of the eight potential reference genes were assayed using RT-qPCR, with lower Cq values reflecting higher mRNA transcript levels. The Cq values of all the potential reference genes were between 17 and 30 under the three treatments (Fig. 2). The expression level of the potential reference gene was different under different treatments. CUL3A had the highest expression level under drought treatment (mean Cq of 21.2), while APRT had the lowest expression level (mean Cq of 26.3). Under the osmotic stress, tubulin had the highest expression level (mean Cq of 21.0), while EF1a had the lowest expression level (mean Cq of 22.8). CUL3A had the highest expression level under simulated drought (mean Cq of 22.7), while sec3 had the lowest expression level (mean Cq of 26.1).

Fig. 2
figure 2

Expression levels of candidate reference genes in experimental samples. Expression data are displayed as Cq values for each reference gene in all samples. The line across the box indicates the median. The box indicates the 25 and 75th percentiles. Whiskers represent the maximum and minimum values. Points represent the average. a Drought, b osmotic, c simulated drought

Candidate reference genes expression stability

geNorm analysis

The gene expression stability of all the 8 candidate reference genes was evaluated using geNorm statistical algorithm. This software determines the normalization value based on the geometric mean of various candidate reference genes and mean pair wise variation of each gene from all the reference genes in a given set of samples. According to geNorm analysis the cut-off range of stability value (M1) is < 1.5, so the gene with lowest M1 value is considered to be the most stable reference gene in terms of gene expression, and vice versa. We analyzed our data for all three experimental sets and found that under drought and simulated drought all the 8 candidate reference genes exhibited high expression stability with low (< 0.8) M1 values, which were much lower than the default threshold of 1.5. However, only 5 genes are stable under osmotic stress, they are EF1α, sec3, tubulin, L8, and actin. Among all treatments, EF1α has the lowest M1 value followed by sec3, however CUL3A exhibited highest M1 value indicating that EF1α and sec3 were most stable in expression and CUL3A the least (Table 2).

Table 2 Ranking of candidate reference genes based on stability values calculated by four softwares for three treatments

NormFinder analysis

The expression stability of 8 candidate genes was further analyzed using NormFinder software. NormFinder measures gene expression stability by comparing the variation within and between user-defined sample groups. Candidate control genes with lowest stability values have the minimum intra and intergroup variation and thus are top ranked. For each candidate gene, NormFinder provides a stability value (M2) that is a direct measurement of expression variation. Hence, it could easily be seen that EF1α, actin and sec3 are the most stable reference genes for the three treatments. Among the most stable reference genes, EF1α had the lowest value which may be considered as the most important reference genes. More interestingly, EF1α has the lowest value under the three treatments consistent with the use of GeNorm analysis. However, there are also slight differences between the results of geNorm and NormFinder analysis. For instance, actin was considered as the most stable reference gene by geNorm under osmotic stress, while it was ranked third by NormFinder (Table 2).

BestKeeper analysis

BestKeeper calculates the BestKeeper Index from the geometric mean of the reference genes and performs Pearson correlation of each of the reference genes to the BestKeeper Index to indicate the correlation of that gene with the Index [28]. BestKeeper also calculates the standard deviation (SD) and the coefficient of variation (CV) based on the Cq values of all candidate reference genes [24]. Genes with SD greater than 1 are considered unacceptable. Reference genes are identified as the most stable genes, i.e. those that exhibit the lowest coefficient of variance and standard deviation (CV ± SD) [29]. The results of BestKeeper analysis are shown in Table 2. The best correlations were obtained for APRT (0.52), actin (0.67), sec3 (0.75) and L8 (0.85) in the drought stress, and for sec3 (0.55), actin (0.80), tubulin (0.93) and EF1α (0.91) in the osmotic treatment, and for tubulin (0.39), L8 (0.53), sec3 (0.58) and EF1α (0.58) in the simulated drought treatment. Under the three treatments, the most unstable genes are identical to CUL3A.

RefFinder analysis

The results obtained from geNorm, NormFinder and BestKeeper were further confirmed using the comprehensive ranking platform RefFinder. RefFinder is a web-based tool which integrates three current computing programs to compare and re-rank the tested reference genes based on the geometric mean of the weights of every single gene calculated by each program. The final ranking results are shown in Table 2. For all the tested samples, a similar ranking order was obtained using RefFinder as compared to geNorm and NormFinder. Under drought treatment, the ranking order (from the most stable to the least stable) was: EF1α > sec3 > actin > GAPDH > APRT > L8 > tubulin > CUL3A. Under osmotic stress, the ranking order was: sec3 > EF1α > actin > L8 > tubulin > GAPDH > APRT > CUL3A. Under simulated drought treatment, the ranking order was: EF1α > sec3 > L8 > actin > tubulin > APRT > GAPDH > CUL3A.

Our combined results show, for normalization, EF1α and sec3 exhibited the best stability under the three treatments, although their ranking was different. On the other hand CUL3A is the most unstable.


According to previous studies on the selection of plant reference genes for RT-qPCR, the expression level of a reference gene might not be constant across various species. Even the stable reference genes of tomato that are homologous to potato are also different [30, 31]. In addition, the expression of a reference gene can be different in the same species in response to various treatments or different plant tissues. For example, EF1α is considered to be the most stable reference gene under potato biotic (late blight) and abiotic stresses (salt stress) [9], but EF1α and APRT are the most stable genes under cold stress. In the potato tubers, the most suitable reference genes are C2, sec3 and CUL3A [11].

In this study, 8 genes that have been commonly used as the candidate reference gene in many species were evaluated. Interestingly, EF1a and sec3 exhibited good stability in potato under drought, osmotic and simulated drought treatments in this study. This indicates that there is an inherent link between the three processes. However, there are some differences in the rankings of potential reference gene stability. For example, the three most stable reference genes are EF1α, sec3 and actin under drought stress, and the three most stable reference genes are sec3, EF1α and actin, while the three most stable reference gene of the simulated drought are EF1α, sec3 and L8.

GeNorm, NormFinder, and BestKeeper are three programs based on statistical analysis. They are commonly used by researchers to assess the robustness of a gene that is used as a reference gene in RT-qPCR analysis. The operating principle of NormFinder is similar to that of the GeNorm program, but the latter can select suitable reference gene combinations and the optimal number of reference genes. In contrast to GeNorm and NormFinder, BestKeeper software directly makes calculations using Cq values [24]. In our analysis, the rankings created by GeNorm and NormFinder were similar, while the ranking obtained by the BestKeeper program was almost always different. A previous report revealed a similar difference between BestKeeper and other methods [32, 33]. The final rank of reference genes was determined with the RefFinder program, a web-based user-friendly comprehensive tool that integrates geNorm, Normfinder and BestKeeper [34]. Therefore, the use of multiple softwares in a comprehensive analysis will help to obtain a more accurate reference gene.

Finally, the ability to control water stress in plants is desirable for drought resistance studies. When plants are growing in soil, studies of the effects of water stress on physiology are usually either very short so that soil and plant water stress change little, or they involve repeated soil drying cycles. For many studies, however, water stress of soil-grown plants cannot be manipulated well enough for careful experiments. Therefore, the addition of osmotic agents, such as PEGs, sorbitol, or mannitol, to liquid nutrient media is very useful and has received considerable attention in the literature [35]. However, the toxicity of osmotic agents has also received increasing attention. PEG toxicity may come from toxic contaminants, and the other may be the accumulation of salt toxic levels. Therefore, we decided to look at additional drought simulation in this study; we planted the plant on quartz sand, and when the plant grows to 6 cm high, the culture medium is dried to avoid the use of osmotic regulators. It is worth noting that the stable reference gene under simulated drought conditions is similar to the other two treatments. In summary, in this study, we analyzed the stability of reference genes for RT-qPCR in potatoes under drought and osmotic stress conditions. The combination of EF1α with sec3 was suitable for gene quantification in potato under drought and osmotic stress. This study also proved that the drought formed by the quartz sand culture has a high degree of similarity in the internal reference gene.


This study represents the first attempt to select a set of commonly used candidate reference genes in potato under drought and osmotic stress for the normalization of gene expression data using RT-qPCR. We showed that the most suitable reference gene is EF1α in drought and simulated drought environments, and the most suitable gene under osmotic stress is sec3. While different pairs were found to be the most appropriate for these biological contexts, we observed that the EF1α and sec3 showed the most stable under all three treatments.


  1. Zhu JK. Abiotic stress signaling and responses in plants. Cell. 2016;167:313–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Sprenger H, Kurowsky C, Horn R, Erban A, Seddig S, Rudack K, et al. The drought response of potato reference cultivars with contrasting tolerance. Plant, Cell Environ. 2016;39:2370–89.

    Article  CAS  Google Scholar 

  3. Penuelas J, Sardans J, Estiarte M, Ogaya R, Carnicer J, Coll M, et al. Evidence of current impact of climate change on life: a walk from genes to the biosphere. Glob Chang Biol. 2013;19:2303–38.

    Article  PubMed  Google Scholar 

  4. Vasquez-Robinet C, Mane SP, Ulanov AV, Watkinson JI, Stromberg VK, De Koeyer D, et al. Physiological and molecular adaptations to drought in Andean potato genotypes. J Exp Bot. 2008;59:2109–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. 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–8.

    Article  CAS  PubMed  Google Scholar 

  6. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55:611–22.

    Article  CAS  PubMed  Google Scholar 

  7. Udvardi MK, Czechowski T, Scheible WR. Eleven golden rules of quantitative RT-PCR. Plant Cell. 2008;20:1736–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 2005;139:5–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Nicot N, Hausman JF, 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–14.

    Article  CAS  PubMed  Google Scholar 

  10. Lopez-Pardo R. Selection of housekeeping genes for qRT-PCR analysis in potato tubers under cold stress. Mol Breed. 2013;31:39–45.

    Article  CAS  Google Scholar 

  11. Mariot RF, de Oliveira LA, Voorhuijzen MM, Staats M, Hutten RC, Van Dijk JP, et al. Selection of reference genes for transcriptional analysis of edible tubers of potato (Solanum tuberosum L.). PLoS ONE. 2015;10:e120854.

    Article  Google Scholar 

  12. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.

    Article  CAS  PubMed  Google Scholar 

  13. Fonseca NA, Marioni J, Brazma A. RNA-Seq gene profiling–a systematic empirical comparison. PLoS ONE. 2014;9:e107026.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Zhang N, Liu B, Ma C, Zhang G, Chang J, Si H, et al. Transcriptome characterization and sequencing-based identification of drought-responsive genes in potato. Mol Biol Rep. 2014;41:505–17.

    Article  CAS  PubMed  Google Scholar 

  15. Yim AK, Wong JW, Ku YS, Qin H, Chan TF, Lam HM. Using RNA-Seq data to evaluate reference genes suitable for gene expression studies in soybean. PLoS ONE. 2015;10:e136343.

    Google Scholar 

  16. Gantasala NP, Papolu PK, Thakur PK, Kamaraju D, Sreevathsa R, Rao U. Selection and validation of reference genes for quantitative gene expression studies by real-time PCR in eggplant (Solanum melongena L.). BMC Res Notes. 2013;6:312.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wang Q, Ishikawa T, Michiue T, Zhu BL, Guan DW, Maeda H. Stability of endogenous reference genes in postmortem human brains for normalization of quantitative real-time PCR data: comprehensive evaluation using geNorm, NormFinder, and BestKeeper. Int J Legal Med. 2012;126:943–52.

    Article  PubMed  Google Scholar 

  18. Sivakumar T, Lan DT, Long PT, Yoshinari T, Tattiyapong M, Guswanto A, et al. PCR detection and genetic diversity of bovine hemoprotozoan parasites in Vietnam. J Vet Med Sci. 2013;75:1455–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Li X, Zhang D, Li H, Gao B, Yang H, Zhang Y, et al. Characterization of reference genes for RT-qPCR in the desert moss Syntrichia caninervis in response to abiotic stress and desiccation/rehydration. Front Plant Sci. 2015;6:38.

    PubMed  PubMed Central  Google Scholar 

  20. Kundu A, Patel A, Pal A. Defining reference genes for qPCR normalization to study biotic and abiotic stress responses in Vigna mungo. Plant Cell Rep. 2013;32:1647–58.

    Article  CAS  PubMed  Google Scholar 

  21. Kong Q, Yuan J, Gao L, Zhao S, Jiang W, Huang Y, et al. Identification of suitable reference genes for gene expression normalization in qRT-PCR analysis in watermelon. PLoS ONE. 2014;9:e90612.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:H34.

    Article  Google Scholar 

  23. 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.

    Article  CAS  PubMed  Google Scholar 

  24. Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004;26:509–15.

    Article  CAS  PubMed  Google Scholar 

  25. Long X, He B, Gao X, Qin Y, Yang J, Fang Y, et al. Validation of reference genes for quantitative real-time PCR during latex regeneration in rubber tree. Gene. 2015;563:190–5.

    Article  CAS  PubMed  Google Scholar 

  26. Huang L, Yan H, Jiang X, Yin G, Zhang X, Qi X, et al. Identification of candidate reference genes in perennial ryegrass for quantitative RT-PCR under various abiotic stress conditions. PLoS ONE. 2014;9:e93724.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Chen Y, Hu B, Tan Z, Liu J, Yang Z, Li Z, et al. Selection of reference genes for quantitative real-time PCR normalization in creeping bentgrass involved in four abiotic stresses. Plant Cell Rep. 2015;34:1825–34.

    Article  CAS  PubMed  Google Scholar 

  28. De Spiegelaere W, Dern-Wieloch J, Weigel R, Schumacher V, Schorle H, Nettersheim D, et al. Reference gene validation for RT-qPCR, a note on different available software packages. PLoS ONE. 2015;10:e122515.

    Google Scholar 

  29. Guo J, Ling H, Wu Q, Xu L, Que Y. The choice of reference genes for assessing gene expression in sugarcane under salinity and drought stresses. Sci Rep. 2014;4:7042.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Lovdal 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–42.

    Article  CAS  PubMed  Google Scholar 

  31. Muller OA, Grau J, Thieme S, Prochaska H, Adlung N, Sorgatz A, et al. Genome-wide identification and validation of reference genes in infected tomato leaves for quantitative RT-PCR analyses. PLoS ONE. 2015;10:e136499.

    Google Scholar 

  32. Mallona I, Lischewski S, Weiss J, Hause B, Egea-Cortines M. Validation of reference genes for quantitative real-time PCR during leaf and flower development in Petunia hybrida. BMC Plant Biol. 2010;10:4.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Rapacz M, Pień AS, Skorupa K. Internal standards for quantitative RT-PCR studies of gene expression under drought treatment in barley (Hordeum vulgare L.): the effects of developmental stage and leaf age. Acta Physiol Plant. 2012;34:1723–33.

    Article  CAS  Google Scholar 

  34. Janska A, Hodek J, Svoboda P, Zamecnik J, Prasil IT, Vlasakova E, et al. The choice of reference gene set for assessing gene expression in barley (Hordeum vulgare L.) under low temperature and drought stress. Mol Genet Genomics. 2013;288:639–49.

    Article  CAS  PubMed  Google Scholar 

  35. Lagerwerff JV, Ogata G, Eagle HE. Control of osmotic pressure of culture solutions with polyethylene glycol. Science. 1961;133:1486–7.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

XT, NZ, and HS designed the experiments. XT performed the experiments and analyzed the data. XT, HS, and AC wrote the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets supporting the conclusions and description of a complete protocol are included within the article.

Consent for publication

All authors have consented for publication.

Ethics approval and consent to participate

Not applicable.


This work was supported by the National Natural Science Foundation of China (31460370, 31660416) and the “Fuxi Talent” Plan of Gansu Agricultural University (FXRC20130102).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Huaijun Si.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tang, X., Zhang, N., Si, H. et al. Selection and validation of reference genes for RT-qPCR analysis in potato under abiotic stress. Plant Methods 13, 85 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: