Front. Plant Sci. Frontiers in Plant Science Front. Plant Sci. 1664-462X Frontiers Media S.A. 10.3389/fpls.2016.01481 Plant Science Original Research Suitable Reference Genes for Accurate Gene Expression Analysis in Parsley (Petroselinum crispum) for Abiotic Stresses and Hormone Stimuli Li Meng-Yao Song Xiong Wang Feng Xiong Ai-Sheng * State Key Laboratory of Crop Genetics and Germplasm Enhancement, College of Horticulture, Nanjing Agricultural University Nanjing, China

Edited by: Fabio Marroni, University of Udine, Italy

Reviewed by: Julia Christine Meitz-Hopkins, Stellenbosch University, South Africa; Palakolanu Sudhakar Reddy, International Crops Research Institute for the Semi-Arid Tropics, India

*Correspondence: Ai-Sheng Xiong xiongaisheng@njau.edu.cn

This article was submitted to Technical Advances in Plant Science, a section of the journal Frontiers in Plant Science

 30 09 2016 2016 7 1481 25 07 2016 20 09 2016 Copyright © 2016 Li, Song, Wang and Xiong. 2016 Li, Song, Wang and Xiong

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Parsley, one of the most important vegetables in the Apiaceae family, is widely used in the food, medicinal, and cosmetic industries. Recent studies on parsley mainly focus on its chemical composition, and further research involving the analysis of the plant's gene functions and expressions is required. qPCR is a powerful method for detecting very low quantities of target transcript levels and is widely used to study gene expression. To ensure the accuracy of results, a suitable reference gene is necessary for expression normalization. In this study, four software, namely geNorm, NormFinder, BestKeeper, and RefFinder were used to evaluate the expression stabilities of eight candidate reference genes of parsley (GAPDH, ACTIN, eIF-4α, SAND, UBC, TIP41, EF-1α, and TUB) under various conditions, including abiotic stresses (heat, cold, salt, and drought) and hormone stimuli treatments (GA, SA, MeJA, and ABA). Results showed that EF-1α and TUB were the most stable genes for abiotic stresses, whereas EF-1α, GAPDH, and TUB were the top three choices for hormone stimuli treatments. Moreover, EF-1α and TUB were the most stable reference genes among all tested samples, and UBC was the least stable one. Expression analysis of PcDREB1 and PcDREB2 further verified that the selected stable reference genes were suitable for gene expression normalization. This study can guide the selection of suitable reference genes in gene expression in parsley.

 parsley reference gene qPCR abiotic stress hormone stimuli 31272175 BK20130027 National Natural Science Foundation of China10.13039/501100001809 Natural Science Foundation of Jiangsu Province10.13039/501100004608 Introduction

Gene expression analysis is an important tool from studying the complex biological processes, such as signal transduction, metabolic pathways, and plant development. In the qualitative or quantitative analysis of the expression of the sample or target gene, qPCR is the most effective method because of its simplicity, high sensitivity and specificity (Wong and Medrano, 2005). However, to date, universal reference gene in plants or animals has not been reported yet (Warrington et al., 2000; Schmittgen and Zakrajsek, 2000). Just borrowed the reference genes in other species might get the wrong expression model. Therefore, evaluation of the stability of reference genes is particularly important.

Housekeeping genes, such as ACTIN, EF-1α, and UBQ, are usually used as candidate reference genes because of their role in maintaining cell survival irrespective of different physiological conditions (Bustin, 2002). However, recent studies showed that the expression patterns of housekeeping genes also change under different experimental conditions (Thellin et al., 1999; Suzuki et al., 2000). In a study of carrot, GAPDH displayed the maximum stability for most of single abiotic stresses (Tian et al., 2015), whereas GAPDH appeared to be the least stable gene during the five developmental stages of the plant (Wang et al., 2016). In addition, the stably expressed reference gene in one species may not be suitable in other species (Andersen et al., 2004; Gutierrez et al., 2008). For example, TUB-B, TUB-A, and UBC are the most stable reference genes during celery development (Li et al., 2016), whereas in cherries, TUB-A is the least stable gene during fruit development (Ye et al., 2015).

Parsley (Petroselinum crispum L.), a member of the Apiaceae family, is a widely cultivated spice. Parsley is also reported as a medicinal plant because it is rich in antioxidants, essential oil, flavonoids, and vitamins (Suhaj, 2006; Zhang et al., 2006). Abiotic stresses, including extreme temperature, salinity, and drought, are the major factors affecting the growth and yield of parsley. The physiological function of the plant is affected by abiotic stresses, and the expression patterns of mRNA and protein also change accordingly (Li et al., 2014). Studies showed that most of the hormones can regulate plant physiological activities, and exogenous plant hormones can improve plant resistance to varying degrees (Eraslan et al., 2007; Bari and Jones, 2009). For example, ABA is called “stress hormone” in plants because it plays an important role in resisting salinity, drought, and low temperature (Narendra, 2007). Plants have established their own resistance strategies in response to extreme adverse conditions (Chen and Zhu, 2004; Katagiri, 2004). Research on the molecular biology of parsley can help in screening the stress resistant genes and accelerating the breeding process. One way of studying the function of stress-resistance genes is to detect their expression levels. However, no systematic study regarding the selection of suitable reference genes for qPCR normalization in parsley has been published to date.

In this study, eight candidate reference genes, including GAPDH, ACTIN, eIF-4α, SAND, UBC, TIP41, EF-1α, and TUB, were assessed by qPCR under various abiotic stresses (cold, heat, salt, and drought) and hormone stimuli treatments (GA, SA, MeJA, and ABA). The expression stabilities of these candidate genes were evaluated by four statistical tools, namely geNorm (Vandesompele et al., 2002), NormFinder (Andersen et al., 2004), Bestkeeper (Pfaffl et al., 2004), and RefFinder (Xie et al., 2012). Additionally, the target genes, PcDREB1 and PcDREB2, were used to identify the best-ranked reference genes. The stable reference genes will enable the more accurate and reliable qPCR analysis of parsley.

 Materials and methods Plant materials and treatments

The curly parsley (P. crispum L. cv. Moxi) seeds were sown in plastic pots containing a soil/vermiculite mixture (3:1) in a controlled-environment growth chamber under a photoperiod of 16 h with 300 μmol m−2s−1 light intensity at 25°C and 8 h dark condition at 16°C. After 8 weeks, the healthy seedlings were used for different treatments. For cold and heat treatments, the seedlings were kept in a growth chamber under 4°C and 38°C for 2 h. For salt and drought treatments, the soils were irrigated with 0.2 M NaCl and 20% PEG6000, respectively. For hormone treatments, the leaves were sprayed with 1.4 mM GA, or 1.4 mM SA, or 0.8 mM MeJA, or 0.1 mM ABA (Jiang et al., 2014). Plants were irrigated or sprayed only once, with untreated plants used as control. All the plants were kept in the same growth conditions for 2 h, and then the leaf samples were harvested. Three biological replicates were performed in different pots for each treatment. All the materials were immediately frozen in liquid nitrogen and then stored at −80°C.

 RNA isolation and cDNA synthesis

Total RNA was extracted from the samples using the RNA simple Total RNA Kit (Tiangen, Beijing, China), and genomic DNA was removed by RNase-free DNase I (Takara, Dalian, China). The quantity and quality of RNA samples were measured by agarose gel electrophoresis and Nanodrop ND 1000 spectrophotometer (Nanodrop Technologies Inc., Delaware, USA). RNA samples with an OD260/280 value between 1.8 and 2.2 were considered as high-purity RNA. A total of 1.0 μg RNA was reverse transcribed into cDNA using PrimeScript RT reagent Kit (Takara, Dalian, China).

 Primer design and qPCR analysis

Based on the transcriptome database of parsley built by our group (Lab of Apiaceae Plant Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, China) (Li et al., 2014), eight candidate reference genes (GAPDH, ACTIN, eIF-4α, SAND, UBC, TIP41, EF-1α, and TUB) were selected as their corresponding homologs genes performed well in other plants in previous studies (Czechowski et al., 2005; Paolacci et al., 2009) and they were cloned in parsley. We submitted the nucleotide sequences to GenBank, and the corresponding accession numbers are KX784033 (eIF-4α), KX784034 (ACTIN), KX784035 (TIP41), KX784036 (GAPDH), KX784037 (SAND), KX784038 (EF-1α), KX784039 (TUB), and KX784040 (UBC).

The reverse transcribed cDNAs were diluted to ten-fold series (10, 102, 103, 104, and 105 X dilutions) for determination of the amplification efficiency (E %) of primers and correlation coefficient (R2), and 18-fold dilution was conducted for qPCR analysis. qPCR was performed using the MyiQ Single Color Real-Time PCR Detection System (Bio-rad, Hercules, CA, USA). The reaction mixture with a total volume of 20 μL contained 2 μL diluted cDNA, 0.4 μL of each 10 μM primer, 7.2 μL ddH2O, and 10 μL SYBR Green I mix (Takara, Dalian, China). The amplification program was 95°C for 30 s, followed by 40 cycles at 94°C for 5 s, and 60°C for 30 s. After amplification, a melting curve (65–95°C with at increments of 0.5°C) was generated for each reaction to verify the specific amplification. A no-template control (cDNAs were replaced with sterilized water) was included in each run for each gene to confirm the absence of nonspecific products. Each PCR reaction was repeated thrice with three biological replicates. Primers for the eight genes were designed by Primer Premier 6 with the following parameters: 57–63°C annealing temperature, 18–22 bp primer length, 40–60% GC contents, and 90–300 bp amplicon length. The primer sequences are listed in Table 1.

Details of candidate reference genes and primers used in qPCR.

Gene Primer sequence (5′–3′) forward/reverse Amplicon length (bp) PCR efficiency (E %) Correlation coefficient (R2) Tm (°C)
ACTIN CTGGATTCTGGTGATGGTGTGA/CTCAGCAGTGGTGGTGAACAT 159 106.0 0.997 82.5
EF-1α AGGCTCTTCAGGAGGCTCTTC/CAATGTGACAGGTGTGGCAATC 216 95.2 0.992 83.5
eIF-4α CACGGAGACATGGATCAGAACAC/GAGACCTGCTGGACATCAATACC 119 95.2 0.993 81.5
GAPDH TCGGACGCATTGGCAGGAA/GCTGGAGTGGATCTCTGTTGGA 228 96.5 0.990 83.5
TIP41 GGAGGACTGTGAGGAACGAATTG/TAAGCACGCCATCAACTCTAAGC 194 108.6 0.992 81.5
TUB ATGGTTCTTGACAATGAGGCACTA/GCTTCCGCAGATCCGAGTTG 160 101.3 0.994 84.0
SAND CAGATGCCGTGTTCTCTTCTCT/CACCAACCAGACTGATGACCTT 192 105.7 0.993 83.0
UBC ATGGCGAATAGCAGCAATCTCC/TGTTATCTTCAGACGGCGATGC 103 101.2 0.993 83.5
Data analysis

The Cq value (quantification cycle) of each reference gene under different conditions was recorded using the qPCR system. Four statistical tools were used to rank to expression stability of the reference genes: geNorm (Vandesompele et al., 2002), NormFinder (Andersen et al., 2004), Bestkeeper (Pfaffl et al., 2004), and RefFinder (Xie et al., 2012). For geNorm and NormFinder, the raw Cq values were converted into the relative quantity values via the formula, 2−ΔCq. geNorm ranks the gene expression stability (M-value) based on the pairwise variation (V) for the gene compared with all other reference genes. The default limit of M-value is 1.5, and the lowest M-value indicates the most stable gene. geNorm also calculates the pairwise variation (Vn/n+1) between normalization factors (NFn and NFn+1, n ≥ 2) to determine the optimal number of reference genes in qPCR normalization. With a Vn/n+1 below 0.15, introducing an additional reference gene is not necessary. NormFinder generates a similar measure through estimated intragroup and intergroup expression variations and combines the variation into a stability value for each gene. BestKeeper ranks the stability based on the calculation of the standard deviation (SD) and the coefficient of variance (CV) with the Cq values. The reference gene with the lowest SD- and CV-values is identified as the most stable gene. RefFinder is a web-based tool that integrates the four computational programs (geNorm, Normfinder, BestKeeper, and the comparative delta-Ct method) to comprehensively rank the tested candidate reference genes.

 Assessment of normalization by qPCR

qPCR was also performed to analyze the expression levels of the target genes, PcDREB1and PcDREB2. The relative expression levels were calculated with the 2−ΔΔCT method (Pfaffl, 2001). The qPCR primer pairs for PcDREB1 were 5′-ACTCGGATCTGGCTAGGCACAT-3′ and 5′-TGGAGAAGCGGTCACCTCATC A-3′, and for PcDREB2, the primer paris were 5′-TTCCAGAACTCGCTCACCACT T-3′ and 5′- GACTCGCTCCATTCGTAACTCG-3′.

 Results Primer specificity and efficiency of candidate reference genes

The eight candidate reference genes (GAPDH, ACTIN, eIF-4α, SAND, UBC, TIP41, EF-1α, and TUB) were selected from the transcriptome database, and they were chosen as their corresponding homologs genes performed well in qPCR normalization. The performances of all primer pairs were tested by standard curve analysis. The amplification efficiencies (E%) ranged from 95.2% to 108.6% and correlation coefficient (R2) ranged from 0.990 to 0.997 (Table 1, Figure S1), the results agreed with the standard (Ramakers et al., 2003). A single peak from all samples in the melting curve analysis showed that the primers were displayed without non-specific amplification (Figure 1).

Melting curves generated for eight candidate reference genes by qPCR in parsley. (A) eIF-4α; (B) ACTIN; (C) TIP41; (D) GAPDH; (E) SAND; (F) EF-1α; (G) TUB; (H) UBC.

 Cq value analysis

The expression levels of the eight candidate reference genes (GAPDH, ACTIN, eIF-4α, SAND, UBC, TIP41, EF-1α, and TUB) were presented as Cq values. To further understand of the Cq values, a distribution diagram was drawn (Figure 2). The Cq values of eight reference genes presented a relatively wide range from 18.94 for EF-1α to 29.92 for UBC in all tested samples. A low Cq value represents high expression level, which indicates that the gene has low expression. In our study, EF-1α was the most abundant gene with the lowest mean Cq value (20.26), whereas TIP41 was the least expressed gene with the highest mean Cq value (29.92). In addition, among all the samples, the SD of TIP41 was the smallest (0.58), whereas UBC showed the maximum variability (SD = 2.28).

Data statistics of Cq values of candidate reference genes in parsley.

 Expression stability of candidate reference genes

In our study, each reference gene was subjected to eight experiment treatments. The eight treatment sets were analyzed individually. To achieve a more comprehensive analysis, these eight sets were divided into three groups for further analysis: “Abiotic stress” (heat, cold, salt, and drought), “Hormone stimuli” (SA, GA, ABA, and MeJA), and “Total” (composed of all the treatment sets). Four different software, namely geNorm, NormFinder, Bestkeeper, and RefFinder, were used to analyze the stability of candidate reference genes.

geNorm analysis

geNorm ranks reference genes by calculating the M-value. A candidate gene with the M-value below 1.5 is considered as a stably expressed gene, and a lower M-value indicates greater stability. The expression stability rankings based on the M-values are displayed in Tables 2, 3. The M-values for the tested genes in all samples and groups were all lower than 1.5. For each individual treatment (Table 2), TUB was the most stably expressed gene with the least M-value under heat, drought, and GA conditions, whereas EF-1α and GAPDH were the most stably expressed under salt, ABA, MeJA, and SA conditions. For cold treatment, GAPDH and eIF-4α were the most stable candidate genes. For the three groups (Table 3), EF-1α and GAPDH were the most stably expressed genes under “Hormone stimuli,” and TUB and EF-1α were the most stable genes under “Abiotic stress” and “Total.” Among all the groups, eIF-4α and UBC were the least stable genes.

Gene expression stability under individual stress conditions ranked by geNorm, NormFinder, and BestKeeper.

Treatments Rank geNorm NormFinder BestKeeper RefFinder
Gene Stability Gene Stability Gene SD CV Gene Comprehensive ranking values
Heat 1 TUB 0.25 TIP41 0.02 GAPDH 0.14 0.57 EF-1α 1.68
2 ACTIN 0.25 EF-1α 0.02 ACTIN 0.69 3.01 TUB 2.63
3 eIF-4α 0.48 ACTIN 0.07 TUB 0.71 2.99 GAPDH 3.22
3 EF-1α 0.49 SAND 0.10 eIF-4α 1.22 5.29 TIP41 5.00
5 SAND 0.58 TUB 0.17 EF-1α 1.25 5.88 eIF-4α 6.24
6 TIP41 0.66 eIF-4α 0.17 SAND 1.42 5.02 SAND 5.96
7 GAPDH 0.79 GAPDH 0.21 TIP41 1.80 6.32 ACTIN 6.74
8 UBC 1.42 UBC 0.26 UBC 4.23 15.49 UBC 8.00
Cold 1 GAPDH 0.15 TUB 0.01 eIF-4α 0.30 1.26 EF-1α 1.57
2 eIF-4α 0.15 eIF-4α 0.01 GAPDH 0.32 1.27 GAPDH 1.73
3 EF-1α 0.42 GAPDH 0.02 ACTIN 0.35 1.48 eIF-4α 2.11
3 TUB 0.44 EF-1α 0.02 EF-1α 0.80 3.69 ACTIN 3.94
5 UBC 0.52 SAND 0.06 TUB 0.81 3.44 TUB 4.47
6 SAND 0.65 UBC 0.06 UBC 1.06 3.49 TIP41 6.00
7 TIP41 0.76 ACTIN 0.12 SAND 1.37 4.80 UBC 7.00
8 ACTIN 0.90 TIP41 0.17 TIP41 1.69 5.92 SAND 8.00
Drought 1 TUB 0.21 EF-1α 0.02 ACTIN 0.17 0.74 ACTIN 1.86
2 GAPDH 0.21 GAPDH 0.03 eIF-4α 0.34 1.41 TUB 1.97
3 EF-1α 0.24 TUB 0.06 EF-1α 1.00 4.69 GAPDH 2.91
3 SAND 0.53 SAND 0.09 TUB 1.06 4.50 EF-1α 3.31
5 TIP41 0.65 eIF-4α 0.12 GAPDH 1.16 4.82 eIF-4α 3.31
6 eIF-4α 0.79 TIP41 0.13 SAND 1.53 5.42 TIP41 5.18
7 ACTIN 0.91 ACTIN 0.18 TIP41 1.89 6.66 SAND 7.00
8 UBC 1.48 UBC 0.22 UBC 4.07 14.83 UBC 8.00
Salt 1 EF-1α 0.22 TUB 0.03 ACTIN 0.21 0.91 EF-1α 1.19
2 GAPDH 0.22 SAND 0.04 eIF-4α 0.29 1.18 ACTIN 1.41
3 TUB 0.31 UBC 0.06 TUB 1.24 5.33 UBC 3.22
3 SAND 0.48 EF-1α 0.11 EF-1α 1.45 6.90 GAPDH 3.94
5 TIP41 0.59 GAPDH 0.20 GAPDH 1.52 6.41 TUB 4.73
6 eIF-4α 1.02 eIF-4α 0.22 SAND 1.83 6.55 SAND 6.00
7 ACTIN 1.20 ACTIN 0.24 TIP41 2.18 7.74 TIP41 7.24
8 UBC 1.45 TIP41 0.26 UBC 3.26 11.54 eIF-4α 7.74
GA 1 TUB 0.18 EF-1α 0.02 eIF-4α 0.72 2.88 EF-1α 1.19
2 GAPDH 0.18 GAPDH 0.02 ACTIN 0.82 3.37 TUB 2.45
3 EF-1α 0.21 TUB 0.03 TUB 1.00 4.24 SAND 3.13
3 TIP41 0.47 TIP41 0.07 EF-1α 1.03 4.79 GAPDH 3.46
5 SAND 0.72 UBC 0.10 GAPDH 1.05 4.34 ACTIN 3.66
6 eIF-4α 1.21 eIF-4α 0.19 TIP41 1.60 5.57 eIF-4α 5.23
7 ACTIN 1.42 SAND 0.23 SAND 2.14 7.73 TIP41 7.00
8 UBC 1.84 ACTIN 0.24 UBC 3.76 13.56 UBC 8.00
ABA 1 EF-1α 0.21 UBC 0.01 GAPDH 0.58 2.36 EF-1α 1.00
2 GAPDH 0.21 EF-1α 0.01 EF-1α 0.62 2.82 GAPDH 1.68
3 TUB 0.32 GAPDH 0.03 eIF-4α 0.64 2.56 TUB 3.00
3 TIP41 0.51 TUB 0.03 ACTIN 0.80 3.30 ACTIN 4.60
5 SAND 0.62 TIP41 0.09 TUB 0.90 3.81 eIF-4α 5.23
6 UBC 0.75 eIF-4α 0.09 SAND 1.21 4.23 UBC 5.63
7 eIF-4α 1.04 SAND 0.11 TIP41 1.30 4.50 TIP41 6.19
8 ACTIN 1.21 ACTIN 0.14 UBC 1.78 5.99 SAND 8.00
MeJA 1 EF-1α 0.21 UBC 0.04 ACTIN 0.27 1.15 EF-1α 1.68
2 GAPDH 0.40 TIP41 0.07 eIF-4α 0.70 2.98 GAPDH 2.06
3 TUB 0.44 EF-1α 0.08 GAPDH 1.57 6.64 eIF-4α 2.45
3 TIP41 0.44 GAPDH 0.09 EF-1α 1.60 7.68 TUB 3.83
5 SAND 0.54 SAND 0.11 TUB 2.01 8.96 UBC 4.16
6 eIF-4α 0.78 TUB 0.15 TIP41 2.08 7.37 SAND 4.86
7 ACTIN 1.02 eIF-4α 0.20 SAND 2.29 8.32 TIP41 6.44
8 TUB 1.41 ACTIN 0.31 UBC 3.99 14.49 ACTIN 8.00
SA 1 EF-1α 0.29 EF-1α 0.02 ACTIN 0.20 0.85 EF-1α 1.19
2 GAPDH 0.29 TUB 0.02 eIF-4α 0.35 1.47 ACTIN 1.41
3 TUB 0.35 GAPDH 0.03 GAPDH 0.74 3.02 TIP41 3.41
3 SAND 0.46 SAND 0.07 EF-1α 0.92 4.26 eIF-4α 3.94
5 eIF-4α 0.56 eIF-4α 0.09 TUB 1.12 4.78 TUB 4.47
6 ACTIN 0.66 ACTIN 0.14 SAND 1.14 3.97 GAPDH 6.24
7 TIP41 0.83 TIP41 0.15 TIP41 1.89 6.67 SAND 6.96
8 UBC 1.47 UBC 0.30 UBC 4.18 15.28 UBC 7.74

Gene expression stability under multiple stress conditions ranked by geNorm,NormFinder, and BestKeeper.

Treatments Rank geNorm NormFinder BestKeeper RefFinder
Gene Stability Gene Stability Gene SD CV Gene Comprehensive ranking values
Abiotic stress 1 TUB 0.48 EF-1α 0.07 ACTIN 0.62 2.65 TUB 1.68
2 EF-1α 0.48 TUB 0.10 TUB 0.66 2.84 EF-1α 1.73
3 SAND 0.60 SAND 0.10 eIF-4α 0.76 3.24 TIP41 2.63
3 TIP41 0.64 TIP41 0.17 GAPDH 0.77 3.14 SAND 2.78
5 GAPDH 0.78 ACTIN 0.18 EF-1α 0.77 3.68 ACTIN 4.73
6 ACTIN 1.00 eIF-4α 0.25 SAND 1.08 3.91 eIF-4α 6.24
7 eIF-4α 1.10 GAPDH 0.26 TIP41 1.35 4.88 GAPDH 6.74
8 UBC 1.53 UBC 0.28 UBC 3.62 13.48 UBC 8.00
Hormone stimuli 1 EF-1α 0.21 EF-1α 0.05 EF-1α 0.84 4.04 EF-1α 1.00
2 GAPDH 0.21 GAPDH 0.09 GAPDH 0.85 3.61 GAPDH 2.28
3 TUB 0.35 TIP41 0.10 ACTIN 0.86 3.57 TUB 2.63
3 TIP41 0.54 TUB 0.13 TUB 0.92 4.07 TIP41 3.36
5 SAND 0.70 eIF-4α 0.14 eIF-4α 1.00 4.08 ACTIN 5.00
6 eIF-4α 0.96 ACTIN 0.20 TIP41 1.21 4.40 eIF-4α 6.24
7 ACTIN 1.08 SAND 0.21 SAND 1.44 5.31 SAND 6.74
8 UBC 1.45 UBC 0.23 UBC 2.96 11.38 UBC 8.00
Total 1 TUB 0.56 EF-1α 0.06 TUB 0.68 3.03 EF-1α 1.19
2 EF-1α 0.56 TIP41 0.14 EF-1α 0.68 3.31 TUB 2.11
3 GAPDH 0.64 TUB 0.15 ACTIN 0.72 3.05 TIP41 2.28
3 TIP41 0.72 ACTIN 0.19 TIP41 0.89 3.27 ACTIN 4.43
5 SAND 0.76 SAND 0.19 eIF-4α 0.96 4.02 SAND 4.73
6 ACTIN 0.98 GAPDH 0.21 GAPDH 0.97 4.10 GAPDH 5.05
7 eIF-4α 1.09 eIF-4α 0.22 SAND 1.02 3.80 eIF-4α 7.00
8 UBC 1.44 UBC 0.27 UBC 2.64 10.28 UBC 8.00

The pairwise variation (Vn/n+1) was calculated to determine the optimal number of reference genes for various groups. A Vn/n+1 below 0.15 indicates that introducing an additional reference gene for normalization is not necessary. As shown in Figure 3, except for the cold and heat treatments, the V2/3 values were lower than 0.15 under the other six conditions, indicating that the two reference genes were sufficient for normalization; three reference genes (V3/4 < 0.15) were proposed to be used in cold and heat conditions. For group “Hormone stimuli,” two reference genes were sufficient for gene normalization with the V2/3 value below the threshold value of 0.15. For groups “Abiotic stress” and “Total,” three and four reference genes were required, respectively. Using more number of reference genes may help in reducing system deviation, but it does not result in more reliable results.

Determination of the optimal numbers of reference genes in parsley.

 NormFinder analysis

The calculation principle of NormFinder differed slightly from geNorm. NormFinder ranks the candidate reference genes based on intragroup and intergroup expression variations. A smaller M-value indicates that the gene is more stable. As shown in Table 2, TIP41 and EF-1α were the two most stable genes under heat condition, eIF-4α and TUB under cold and salt conditions, EF-1α under salt, SA, and GA treatments, and UBC under ABA and MeJA treatments. For the combination analysis (Table 3), EF-1α ranked the first in all three groups. Similar to the results generated by geNorm, EF-1α and TUB displayed good performance, whereas UBC was the least stable reference gene.

 BestKeeper analysis

BestKeeper ranks the reference genes by calculating the SD and the CV of their Cq values, and small SD- and CV-values indicate that the gene is more stable. For the eight individual treatments (Table 2), ACTIN had the lowest SD- and CV-values under drought, salt, MeJA, and SA treatments. eIF-4α was the most stable gene in cold and GA treatments. For the three groups (Table 3), both ACTIN and TUB were the most stable genes in “Abiotic stress,” and EF-1α and GAPDH were the most stable in “Hormone stimuli.” In “Total,” TUB and EF-1α were more stable than other genes. In addition, we also discovered that UBC was the least stable gene in most sets.

 RefFinder

RefFinder was used to integrate a comprehensive ranking of the most stable candidate genes. However, the ranking orders generated by geNorm, NormFinder, BestKeeper, and RefFinder were not entirely consistent for the eight individual treatments (Table 2). Comprehensive ranking for the three groups revealed that EF-1α and TUB were the most stable genes in “Abiotic stresses,” while EF-1α, GAPDH, and TUB had good performance in “Hormone stimuli” (Table 3). In “Total,” EF-1α and TUB can be recommended as the optimal reference genes for normalization in various conditions.

 Validation of the best and least ranked reference genes

The DREB transcription factors are important members of the AP2/ERF family, which play important roles in the regulation of plant stress response (Liu et al., 2000; Agarwal et al., 2006; Zhuang et al., 2008; Wu et al., 2015). In Arabidopsis, a number of DREB genes have been induced to be expressed under drought and cold stresses (Seki et al., 2001). Overexpression of a cotton DREB gene, GhDREB, can increase drought, salt, and freezing tolerance of transgenic wheat (Gao et al., 2009). To validate the selected reference genes, the expression levels of PcDREB1 and PcDREB2 genes were evaluated by qPCR under cold and drought conditions. Four reference genes, including the two most stable genes, TUB and EF-1α, one less stable gene eIF-4α, and the least stable gene UBC, were selected to normalize the expression of target genes.

As shown in Figure 4A, under cold condition, PcDREB1 showed similar response pattern when normalized by the four reference genes (TUB, EF-1α, eIF-4α, and UBC): expression level initially increased, reached the highest level at 8 h, and finally decreased afterward. Meanwhile, the expression levels during normalization with TUB and EF-1α were maintained at high levels at 24 h. In response to drought treatment, the change in the expression level of PcDREB1 was different (Figure 4B). Normalization with TUB and EF-1α indicated a 2- or 3-fold increase after 1 h of drought treatment and maintenance of high levels at 2 h, whereas normalization with UBC indicated little change during these stages. Differences in results were also observed in PcDREB2. The expression patterns of PcDREB2 were similar when the most stable genes TUB and EF-1α were used for normalization, but varied greatly during normalization with the least stable gene UBC (Figure 4C). Under drought stress, PcDREB2 responded quickly to drought, and high expression level was maintained at early stages (1–4 h) when using TUB and EF-1α as reference genes. On the other hand, expression levels of PcDREB2 at 4 and 8 h was higher when using eIF-4α and UBC as reference genes (Figure 4D). Overall, the expression patterns of PcDREB1 or PcDREB2 were consistent when using the TUB and EF-1α, whereas the expression levels demonstrated some deviations when the less stable genes were used. In view of these results, TUB and EF-1α were considered as the suitable reference genes for normalizing qPCR data.

Relative expression patterns of PcDREB1 (A,B) and PcDREB2 (C,D) under cold and drought conditions. TUB, EF-1α, eIF-4α, and UBC were used as reference genes for expression normalization. Different letters on the vertical bars indicate significant difference at 0.05 levels.

 Discussion

qPCR is one of the most sensitive methods that can detect the low expression of target genes (Bustin, 2000). In qPCR quantification, reference genes are required for eliminating the errors during mRNA extraction, amplification efficiency, qPCR procedure, and so on (Vanguilder et al., 2008; Derveaux et al., 2010). An ideal reference gene should satisfy the following requirements: stable expression in different tissue or organ, expression is not affected by any experimental conditions, gene is expressed to a certain degree, and the expression level is similar to that of the target gene (Thellin et al., 1999; Suzuki et al., 2000; Udvardi et al., 2008). Some genes, which are involved in cytoskeleton structure (ACTIN and TUB), protein synthesis (EF-1α and eIF-4α), and biological metabolic processes (GAPDH and UBQ), are called housekeeping genes and are usually used as reference genes (Rebouças et al., 2013). To date, with the obtained genome and transcriptome data, some new generation reference genes such as Fb15 (fiber protein15), ABCT (ATP-binding cassette transporter), and CAC (clathrin adaptor complexes medium) were shown to have stable expression similar to the traditional reference genes (Xu et al., 2015; Reddy D. S. et al., 2016).

In this study, using parsley, the stability of eight reference genes (GAPDH, ACTIN, eIF-4α, SAND, UBC, TIP41, EF-1α, and TUB) under different abiotic stresses and hormone stimuli was analyzed. Three software, geNorm (Vandesompele et al., 2002), NormFinder (Andersen et al., 2004), BestKeeper (Pfaffl et al., 2004), were used to evaluate the expression stability of candidate reference genes. Results obtained from geNorm, NormFinder, and BestKeeper were not consistent, especially under specific individual conditions. For example, in heat treatment, GAPDH was the most stable gene in BestKeeper analysis, yet GAPDH preformed unsatisfactorily in geNorm and NormFinder. In cold treatment, geNorm and BestKeeper both recommended GAPDH and eIF-4α as the suitable reference genes, whereas TUB and eIF-4α were identified as suitable genes by NormFinder. However, no significant difference was observed in the expression stability between TUB and GAPDH under cold conditions. Therefore, we inferred that these three genes can be used for the correction of relative gene expression. The divergences among the three software were possibly due to the differences in algorithms. However, synthesizing several multiple results from different software can minimize the errors on the selection of reference genes. Two or more statistical software were used for reference gene selection in previous literature, geNorm, NormFinder, and BestKeeper were the three most commonly used software (Qi et al., 2016; Reddy P. S. et al., 2016). RefFinder (Xie et al., 2012) was used to generate a comprehensive ranking of candidate reference genes. Although ACTIN ranked the first in BestKeeper analysis under abiotic stress condition, by comprehensive analysis, RefFinder recommended TUB and EF-1α as the most stable reference genes. Results from geNorm, NormFinder, BestKeeper, and RefFinder in “Hormone stimuli” were consistent, which showed that EF-1α and GAPDH were the most stably expressed genes. We also found that TUB had good performance in “Hormone stimuli.” By comprehensive comparison, EF-1α and TUB can be used as the most suitable reference genes in various conditions.

In most previous studies, reference genes were not always stable in different tissues, organs, genotypes, and experimental conditions (Schmittgen and Zakrajsek, 2000). In rice, UBQ5 and eEF-1α are the most stable genes in all the tissue samples, and 18S and 25S rRNAs are the most stable genes under various treatment conditions (Jain et al., 2006). The other study showed that eIF-4α and ACT1 performed well at different developmental stages and different varieties of rice, while 18S and 25S rRNAs had the least stable expression (Li et al., 2010). Moreover, the selections of reference genes are not the same in various species. Carrot is another important vegetable of the Apiaceae family. Tian et al. (2015) discovered that ACTIN and TUB in carrot are the most suitable reference genes among “Abiotic stress” and “Hormone stimuli,” but ACTIN was not the suitable choice in our study. GAPDH showed good performance in “Hormone stimuli” in parsley but was the least stable in carrot. However, TUB performed satisfactorily in these two plants. To accurately analyze the expression patterns of the target genes, selecting the proper reference genes according to the species and experimental requirement is necessary.

Plant growth and development are generally affected by abiotic stresses, such as extreme temperature, drought or high-salinity conditions. Increasing evidence have demonstrated that AP2/ERF transcription factors play important roles in plant development and stress response (Xu et al., 2011; Licausi et al., 2013; Zhuang et al., 2014). DREB proteins are members of a class of AP2/ERF transcription factor family that can directly regulate various stress related genes by binding the DNA with the DRE/CRT element, which can increase plant resistance to abiotic stress (Sakuma et al., 2002). The functions of DREB genes have been characterized in many plants. For example, overexpression of ZmDREB1A in transgenic plants can increase tolerance of plants to drought and freezing stresses (Qin et al., 2004). Similarly, overexpression of GmDREB gene from soybean resulted in enhanced drought and cold stress tolerance in transgenic wheat plants (Gao et al., 2005). In this study, the expression levels of PcDREB1 and PcDREB2 were assessed under cold and drought conditions to validate the selected reference genes. When using the two most stable genes TUB and EF-1α, the expression of the target gene was consistent. PcDREB1 and PcDREB2 genes were up-regulated under cold and drought stresses, which were in agreement with other DREB genes responding to cold and drought stresses (Shinwari et al., 1998; Ito et al., 2006). Normalization with the less stable gene eIF-4α showed little change in the expression of the gene. In contrast, some divergences were observed in the expression patterns, which normalized by the least stable reference gene UBC. The observations were as follows: the expression of PcDREB1 had no significant change under drought condition, and the expression level of PcDREB2 under cold stress increased initially, decreased at 2 h, increased again, and finally declined. Results demonstrated that using an unstable reference gene for normalization may contribute to inaccurate results.

In conclusion, we evaluated the stabilities of eight candidate reference genes during treatments with various abiotic stresses (heat, cold, salt, and drought) and hormone stimuli treatments (GA, SA, MeJA, and ABA). Our study proposed that EF-1α and TUB are the suitable genes under abiotic stresses; EF-1α, GAPDH, and TUB are suitable for normalization during hormone stimuli treatments. Overall, EF-1α and TUB can be used under all above conditions. In addition, to verify the selected reference genes, the expression patterns of PcDREB1 and PcDREB2 were analyzed. The reference genes selected in this study provide more choices in genes expression analysis and functional studies in parsley.

 Author contributions

AX and ML initiated and designed the research. ML, XS, and FW performed the experiments. ML, XS, FW, and AX analyzed the data. AX contributed reagents/materials/analysis tools. ML wrote the paper. ML and AX revised the paper.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The research was supported by the National Natural Science Foundation of China (31272175), Natural Science Foundation of Jiangsu Province (BK20130027) and Priority Academic Program Development of Jiangsu Higher Education Institutions.

 Supplementary material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.01481

Standard curves of eight candidate reference genes in parsley. (A) eIF-4α; (B) ACTIN; (C) TIP41; (D) GAPDH; (E) SAND; (F) EF-1α; (G) TUB; (H) UBC.

 References Agarwal P. K. Agarwal P. Reddy M. K. Sopory S. K. (2006). Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Rep. 25, 1263–1274. 10.1007/s00299-006-0204-816858552 Andersen C. L. Jensen J. L. Ørntoft T. F. (2004). 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. 64, 5245–5250. 10.1158/0008-5472.CAN-04-049615289330 Bari R. Jones J. D. G. (2009). Role of plant hormones in plant defence responses. Plant Mol. Biol. 69, 473–488. 10.1007/s11103-008-9435-019083153 Bustin S. A. (2000). Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J. Mol. Endocrinol. 25, 169–193. 10.1677/jme.0.025016911013345 Bustin S. A. (2002). Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J. Mol. Endocrinol. 29, 23–39. 10.1677/jme.0.029002312200227 Chen W. J. Zhu T. (2004). Networks of transcription factors with roles in environmental stress response. Trends Plant Sci. 9, 591–596. 10.1016/j.tplants.2004.10.00715564126 Czechowski T. Stitt M. Altmann T. Udvardi M. K. Scheible W. R. (2005). Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 139, 5–17. 10.1104/pp.105.06374316166256 Derveaux S. Vandesompele J. Hellemans J. (2010). How to do successful gene expression analysis using real-time PCR. Methods 50, 227–230. 10.1016/j.ymeth.2009.11.00119969088 Eraslan F. Inal A. Gunes A. Alpaslan M. (2007). Impact of exogenous salicylic acid on the growth, antioxidant activity and physiology of carrot plants subjected to combined salinity and boron toxicity. Sci. Hortic. 113, 120–128. 10.1016/j.scienta.2007.03.012 Gao S. Q. Chen M. Xia L. Q. Xiu H. J. Xu Z. S. Li L. C. . (2009). A cotton (Gossypium hirsutum) DRE-binding transcription factor gene, GhDREB, confers enhanced tolerance to drought, high salt, and freezing stresses in transgenic wheat. Plant Cell Rep. 28, 301–311. 10.1007/s00299-008-0623-919005655 Gao S. Xu H. Cheng X. Chen M. Xu Z. Li L. . (2005). Improvement of wheat drought and salt tolerance by expression of a stress-inducible transcription factor GmDREB of soybean (Glycine max). Chin. Sci. Bull. 50, 2714–2723. 10.1360/982005-1234 Gutierrez L. Mauriat M. Guénin S. Pelloux J. Lefebvre J. F. Louvet R. . (2008). The lack of a systematic validation of reference genes: a serious pitfall undervalued in reverse transcription-polymerase chain reaction (RT-PCR) analysis in plants. Plant Biotechnol. J. 6, 609–618. 10.1111/j.1467-7652.2008.00346.x18433420 Ito Y. Katsura K. Maruyama K. Taji T. Kobayashi M. Seki M. . (2006). Functional analysis of rice DREB1/CBF-type transcription factors involved in cold-responsive gene expression in transgenic rice. Plant Cell Physiol. 47, 141–153. 10.1093/pcp/pci23016284406 Jain M. Nijhawan A. Tyagi A. K. Khurana J. P. (2006). Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochem. Biophys. Res. Commun. 345, 646–651. 10.1016/j.bbrc.2006.04.14016690022 Jiang Q. Wang F. Li M. Y. Ma J. Tan G. F. Xiong A. S. (2014). Selection of suitable reference genes for qPCR normalization under abiotic stresses in Oenanthe javanica (BI.) DC. PLoS ONE 9:e92262. 10.1371/journal.pone.009226224651080 Katagiri F. (2004). A global view of defense gene expression regulation – a highly interconnected signaling network. Curr. Opin. Plant Biol. 7, 506–511. 10.1016/j.pbi.2004.07.01315337092 Li M. Y. Tan H. W. Wang F. Jiang Q. Xu Z. S. Tian C. . (2014). De novo transcriptome sequence assembly and identification of AP2/ERF transcription factor related to abiotic stress in parsley (Petroselinum crispum). PLoS ONE 9:e108977. 10.1371/journal.pone.010897725268141 Li M. Y. Wang F. Jiang Q. Wang G. L. Tian C. Xiong A. S. (2016). Validation and comparison of reference genes for qPCR normalization of celery (Apium graveolens) at different development stages. Front. Plant Sci. 7:313. 10.3389/fpls.2016.0031327014330 Li Q. F. Sun S. S. M. Yuan D. Y. Yu H. X. Gu M. H. Liu Q. Q. (2010). Validation of candidate reference genes for the accurate normalization of real-time quantitative RT-PCR data in rice during seed development. Plant Mol. Biol. Rep. 28, 49–57. 10.1007/s11105-009-0124-1 Licausi F. Ohme-Takagi M. Perata P. (2013). APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors: mediators of stress responses and developmental programs. New Phytol. 199, 639–649. 10.1111/nph.1229124010138 Liu Q. Zhao N. Yamaguch-Shinozaki K. Shinozaki K. (2000). Regulatory role of DREB transcription factors in plant drought, salt and cold tolerance. Chin. Sci. Bull. 45, 970–975. 10.1007/BF02884972 Narendra T. (2007). Abscisic acid and abiotic stress signaling. Plant Signal. Behav. 2, 135–138. 10.4161/psb.2.3.415619516981 Paolacci A. R. Tanzarella O. A. Porceddu E. Ciaffi M. (2009). Identification and validation of reference genes for quantitative RT-PCR normalization in wheat. BMC Mol. Biol. 10:11. 10.1186/1471-2199-10-1119232096 Pfaffl M. W. (2001). A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29:e45. 10.1093/nar/29.9.e4511328886 Pfaffl M. W. Tichopad A. Prgomet C. Neuvians T. P. (2004). Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol. Lett. 26, 509–515. 10.1023/B:BILE.0000019559.84305.4715127793 Qi S. Yang L. Wen X. Hong Y. Song X. Zhang M. . (2016). Reference gene selection for RT-qPCR analysis of flower development in Chrysanthemum morifolium and Chrysanthemum lavandulifolium. Front. Plant Sci. 7:287. 10.3389/fpls.2016.0028727014310 Qin F. Sakuma Y. Li J. Liu Q. Li Y. Q. Shinozaki K. . (2004). Cloning and functional analysis of a novel DREB1/CBF transcription factor involved in cold-responsive gene expression in Zea mays L. Plant Cell Physiol. 45, 1042–1052. 10.1093/pcp/pch11815356330 Ramakers C. Ruijter J. M. Deprez R. H. L. Moorman A. F. M. (2003). Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci. Lett. 339, 62–66. 10.1016/S0304-3940(02)01423-412618301 Rebouças E. D. L. Passos M. J. Passos J. R. D. S. Hurk R. V. D. Silva J. R. V. (2013). Real time PCR and importance of housekeepings genes for normalization and quantification of mRNA expression in different tissues. Braz. Arch. Biol. Technol. 56, 143–154. 10.1590/S1516-89132013000100019 Reddy D. S. Bhatnagar-Mathur P. Reddy P. S. Sri C. K. Sivaji G. A. Sharma K. K. (2016). Identification and validation of reference genes and their impact on normalized gene expression studies across cultivated and wild cicer species. PLoS ONE 11:e0148451. 10.1371/journal.pone.014845126863232 Reddy P. S. Reddy D. S. Sivasakthi K. Bhatnagar-Mathur P. Vadez V. Sharma K. K. (2016). Evaluation of sorghum [Sorghum bicolor (L.)] reference genes in various tissues and under abiotic stress conditions for quantitative real-time PCR data normalization. Front. Plant Sci. 7:529. 10.3389/fpls.2016.0052927200008 Sakuma Y. Liu Q. Dubouzet J. Abe H. Shinozaki K. YamaguchiShinozaki K. (2002). DNA-binding specificity of the ERF/AP2 domain of Arabidopsis DREBs, transcription factors involved in dehydration- and cold-inducible gene expression. Biochem. Biophys. Res. Commun. 290, 998–1009. 10.1006/bbrc.2001.629911798174 Schmittgen T. D. Zakrajsek B. A. (2000). Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. J. Biochem. Biophys. Methods 46, 69–81. 10.1016/S0165-022X(00)00129-911086195 Seki M. Narusaka M. Abe H. Kasuga M. Yamaguchi-Shinozaki K. Carninci P. . (2001). Monitoring the expression pattern of 1300 Arabidopsis genes under drought and cold stresses by using a full-length cDNA microarray. Plant Cell 13, 61–72. 10.1105/tpc.13.1.6111158529 Shinwari Z. K. Nakashima K. Miura S. Kasuga M. Seki M. Yamaguchi-Shinozaki K. . (1998). An Arabidopsis gene family encoding DRE/CRT binding proteins involved in low-temperature-responsive gene expression. Biochem. Biophys. Res. Commun. 250, 161–170. 10.1006/bbrc.1998.92679735350 Suhaj M. (2006). Spice antioxidants isolation and their antiradical activity: a review. J. Food Compost. Anal. 19, 531–537. 10.1016/j.jfca.2004.11.005 Suzuki T. Higgins P. J. Crawford D. R. (2000). Control selection for RNA quantitation. BioTechniques 29, 332–337. 10948434 Thellin O. Zorzi W. Lakaye B. Borman B. D. Coumans B. Hennen G. . (1999). Housekeeping genes as internal standards: use and limits. J. Biotechnol. 75, 291–295. 10.1016/S0168-1656(99)00163-710617337 Tian C. Jiang Q. Wang F. Wang G. L. Xu Z. S. Xiong A. S. (2015). Selection of suitable reference genes for qPCR normalization under abiotic stresses and hormone stimuli in carrot leaves. PLoS ONE 10:e0117569. 10.1371/journal.pone.011756925658122 Udvardi M. K. Czechowski T. Scheible W. (2008). Eleven golden rules of quantitative RT-PCR. Plant Cell 20, 1736–1737. 10.1105/tpc.108.06114318664613 Vandesompele J. Preter K. D. Pattyn F. Poppe B. Roy N. V. Paepe A. D. . (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3:research0034. 10.1186/gb-2002-3-7-research003412184808 Vanguilder H. D. Vrana K. E. Freeman W. M. (2008). Twenty-five years of quantitative PCR for gene expression analysis. BioTechniques 44, 619–626. 10.2144/00011277618474036 Wang G. L. Tian C. Jiang Q. Xu Z. S. Wang F. Xiong A. S. (2016). Comparison of nine reference genes for real-time quantitative PCR in roots and leaves during five developmental stages in carrot (Daucus carota L.). J. Hortic. Sci. Biotechnol. 91, 264–270. 10.1080/14620316.2016.1148372 Warrington J. A. Nair A. Mahadevappa M. Tsyganskaya M. (2000). Comparison of human adult and fetal expression and identification of 535 housekeeping/maintenance genes. Physiol. Genomics 2, 143–147. 11015593 Wong M. L. Medrano J. F. (2005). Real-time PCR for mRNA quantitation. BioTechniques 39, 75–85. 10.2144/05391RV0116060372 Wu Z. J. Li X. H. Liu Z. W. Li H. Wang Y. X. Zhuang J. (2015). Transcriptome-based discovery of AP2/ERF transcription factors related to temperature stress in tea plant (Camellia sinensis). Funct. Integr. Genomics 15, 741–752. 10.1007/s10142-015-0457-926233577 Xie F. Peng X. Chen D. Xu L. Zhang B. (2012). miRDeepfinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol. Biol. 80, 75–84. 10.1007/s11103-012-9885-2 Xu H. Bao J. D. Dai J. S. Li Y. Zhu Y. (2015). Genome-wide identification of new reference genes for qRT-PCR normalization under high temperature stress in rice endosperm. PLoS ONE 10:e0142015. 10.1371/journal.pone.014201526555942 Xu Z. S. Chen M. Li L. C. Ma Y. Z. (2011). Functions and application of the AP2/ERF transcription factor family in crop improvement. J. Integr. Plant Biol. 53, 570–585. 10.1111/j.1744-7909.2011.01062.x21676172 Ye X. Zhang F. Tao Y. Song S. Fang J. (2015). Reference gene selection for quantitative real-time PCR normalization in different cherry genotypes, developmental stages and organs. Sci. Hortic. 181, 182–188. 10.1016/j.scienta.2014.10.027 Zhang H. Chen F. Wang X. Yao H. Y. (2006). Evaluation of antioxidant activity of parsley (Petroselinum crispum) essential oil and identification of its antioxidant constituents. Food Res. Int. 39, 833–839. 10.1016/j.foodres.2006.03.007 Zhuang J. Cai B. Peng R. H. Zhu B. Jin X. F. Xue Y. . (2008). Genome-wide analysis of the AP2/ERF gene family in Populus trichocarpa. Biochem. Biophys. Res. Commun. 371, 468–474. 10.1016/j.bbrc.2008.04.08718442469 Zhuang J. Zhang J. Hou X. L. Wang F. Xiong A. S. (2014). Transcriptomic, proteomic, metabolomic and functional genomic approaches for the study of abiotic stress in vegetable crops. Crit. Rev. Plant Sci. 33, 225–237. 10.1080/07352689.2014.870420 Abbreviations ABA

abscisic acid

 eIF-4α

Eukaryotic translation initiation factor 4α-1 gene

 EF-1α

Elongation factor-1α gene

 GA

gibberellins

 GAPDH

Glyceraldehyde-3-phosphate dehydrogenase gene

 MeJA

methyl jasmonate

 qPCR

quantitative real-time reverse transcription PCR

 SA

salicylic acid

 TIP41

Tap42-interacting protein of 41 kDa gene.