Suitable Reference Genes for Accurate Gene Expression Analysis in Parsley (Petroselinum crispum) for Abiotic Stresses and Hormone Stimuli

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.

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⁻²s⁻¹ 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 OD_260/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, 10², 10³, 10⁴, and 10⁵ X dilutions) for determination of the amplification efficiency (E %) of primers and correlation coefficient (R²), 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.

TABLE 1

Table 1. Details of candidate reference genes and primers used in qPCR.

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 (V_n/n+1) between normalization factors (NFn and NFn+1, n ≥ 2) to determine the optimal number of reference genes in qPCR normalization. With a V_n/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 (R²) 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).

FIGURE 1

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

FIGURE 2

Figure 2. 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.

TABLE 2

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

TABLE 3

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

The pairwise variation (V_n/n+1) was calculated to determine the optimal number of reference genes for various groups. A V_n/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 V_2/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 (V_3/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 V_2/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.

FIGURE 3

Figure 3. 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.

FIGURE 4

Figure 4. 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.

Acknowledgments

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

Figure S1. 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.

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.

References

Agarwal, P. K., Agarwal, P., Reddy, M. K., and Sopory, S. K. (2006). Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Rep. 25, 1263–1274. doi: 10.1007/s00299-006-0204-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersen, C. L., Jensen, J. L., and Ø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. doi: 10.1158/0008-5472.CAN-04-0496

PubMed Abstract | CrossRef Full Text | Google Scholar

Bari, R., and Jones, J. D. G. (2009). Role of plant hormones in plant defence responses. Plant Mol. Biol. 69, 473–488. doi: 10.1007/s11103-008-9435-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Bustin, S. A. (2000). Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J. Mol. Endocrinol. 25, 169–193. doi: 10.1677/jme.0.0250169

PubMed Abstract | CrossRef Full Text | Google Scholar

Bustin, S. A. (2002). Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J. Mol. Endocrinol. 29, 23–39. doi: 10.1677/jme.0.0290023

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W. J., and Zhu, T. (2004). Networks of transcription factors with roles in environmental stress response. Trends Plant Sci. 9, 591–596. doi: 10.1016/j.tplants.2004.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Czechowski, T., Stitt, M., Altmann, T., Udvardi, M. K., and Scheible, W. R. (2005). Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 139, 5–17. doi: 10.1104/pp.105.063743

PubMed Abstract | CrossRef Full Text | Google Scholar

Derveaux, S., Vandesompele, J., and Hellemans, J. (2010). How to do successful gene expression analysis using real-time PCR. Methods 50, 227–230. doi: 10.1016/j.ymeth.2009.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Eraslan, F., Inal, A., Gunes, A., and 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. doi: 10.1016/j.scienta.2007.03.012

CrossRef Full Text | Google Scholar

Gao, S. Q., Chen, M., Xia, L. Q., Xiu, H. J., Xu, Z. S., Li, L. C., et al. (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. doi: 10.1007/s00299-008-0623-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, S., Xu, H., Cheng, X., Chen, M., Xu, Z., Li, L., et al. (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. doi: 10.1360/982005-1234

CrossRef Full Text | Google Scholar

Gutierrez, L., Mauriat, M., Guénin, S., Pelloux, J., Lefebvre, J. F., Louvet, R., et al. (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. doi: 10.1111/j.1467-7652.2008.00346.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ito, Y., Katsura, K., Maruyama, K., Taji, T., Kobayashi, M., Seki, M., et al. (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. doi: 10.1093/pcp/pci230

PubMed Abstract | CrossRef Full Text | Google Scholar

Jain, M., Nijhawan, A., Tyagi, A. K., and 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. doi: 10.1016/j.bbrc.2006.04.140

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Q., Wang, F., Li, M. Y., Ma, J., Tan, G. F., and Xiong, A. S. (2014). Selection of suitable reference genes for qPCR normalization under abiotic stresses in Oenanthe javanica (BI.) DC. PLoS ONE 9:e92262. doi: 10.1371/journal.pone.0092262

PubMed Abstract | CrossRef Full Text | Google Scholar

Katagiri, F. (2004). A global view of defense gene expression regulation – a highly interconnected signaling network. Curr. Opin. Plant Biol. 7, 506–511. doi: 10.1016/j.pbi.2004.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M. Y., Tan, H. W., Wang, F., Jiang, Q., Xu, Z. S., Tian, C., et al. (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. doi: 10.1371/journal.pone.0108977

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M. Y., Wang, F., Jiang, Q., Wang, G. L., Tian, C., and 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. doi: 10.3389/fpls.2016.00313

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q. F., Sun, S. S. M., Yuan, D. Y., Yu, H. X., Gu, M. H., and 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. doi: 10.1007/s11105-009-0124-1

CrossRef Full Text | Google Scholar

Licausi, F., Ohme-Takagi, M., and Perata, P. (2013). APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors: mediators of stress responses and developmental programs. New Phytol. 199, 639–649. doi: 10.1111/nph.12291

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Zhao, N., Yamaguch-Shinozaki, K., and Shinozaki, K. (2000). Regulatory role of DREB transcription factors in plant drought, salt and cold tolerance. Chin. Sci. Bull. 45, 970–975. doi: 10.1007/BF02884972

CrossRef Full Text | Google Scholar

Narendra, T. (2007). Abscisic acid and abiotic stress signaling. Plant Signal. Behav. 2, 135–138. doi: 10.4161/psb.2.3.4156

PubMed Abstract | CrossRef Full Text | Google Scholar

Paolacci, A. R., Tanzarella, O. A., Porceddu, E., and Ciaffi, M. (2009). Identification and validation of reference genes for quantitative RT-PCR normalization in wheat. BMC Mol. Biol. 10:11. doi: 10.1186/1471-2199-10-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29:e45. doi: 10.1093/nar/29.9.e45

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaffl, M. W., Tichopad, A., Prgomet, C., and 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. doi: 10.1023/B:BILE.0000019559.84305.47

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, S., Yang, L., Wen, X., Hong, Y., Song, X., Zhang, M., et al. (2016). Reference gene selection for RT-qPCR analysis of flower development in Chrysanthemum morifolium and Chrysanthemum lavandulifolium. Front. Plant Sci. 7:287. doi: 10.3389/fpls.2016.00287

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, F., Sakuma, Y., Li, J., Liu, Q., Li, Y. Q., Shinozaki, K., et al. (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. doi: 10.1093/pcp/pch118

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramakers, C., Ruijter, J. M., Deprez, R. H. L., and Moorman, A. F. M. (2003). Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci. Lett. 339, 62–66. doi: 10.1016/S0304-3940(02)01423-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Rebouças, E. D. L., Passos, M. J., Passos, J. R. D. S., Hurk, R. V. D., and 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. doi: 10.1590/S1516-89132013000100019

CrossRef Full Text | Google Scholar

Reddy, D. S., Bhatnagar-Mathur, P., Reddy, P. S., Sri, C. K., Sivaji, G. A., and 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. doi: 10.1371/journal.pone.0148451

PubMed Abstract | CrossRef Full Text | Google Scholar

Reddy, P. S., Reddy, D. S., Sivasakthi, K., Bhatnagar-Mathur, P., Vadez, V., and 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. doi: 10.3389/fpls.2016.00529

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakuma, Y., Liu, Q., Dubouzet, J., Abe, H., Shinozaki, K., and 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. doi: 10.1006/bbrc.2001.6299

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmittgen, T. D., and 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. doi: 10.1016/S0165-022X(00)00129-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Seki, M., Narusaka, M., Abe, H., Kasuga, M., Yamaguchi-Shinozaki, K., Carninci, P., et al. (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. doi: 10.1105/tpc.13.1.61

PubMed Abstract | CrossRef Full Text | Google Scholar

Shinwari, Z. K., Nakashima, K., Miura, S., Kasuga, M., Seki, M., Yamaguchi-Shinozaki, K., et al. (1998). An Arabidopsis gene family encoding DRE/CRT binding proteins involved in low-temperature-responsive gene expression. Biochem. Biophys. Res. Commun. 250, 161–170. doi: 10.1006/bbrc.1998.9267

PubMed Abstract | CrossRef Full Text | Google Scholar

Suhaj, M. (2006). Spice antioxidants isolation and their antiradical activity: a review. J. Food Compost. Anal. 19, 531–537. doi: 10.1016/j.jfca.2004.11.005

CrossRef Full Text | Google Scholar

Suzuki, T., Higgins, P. J., and Crawford, D. R. (2000). Control selection for RNA quantitation. BioTechniques 29, 332–337.

PubMed Abstract | Google Scholar

Thellin, O., Zorzi, W., Lakaye, B., Borman, B. D., Coumans, B., Hennen, G., et al. (1999). Housekeeping genes as internal standards: use and limits. J. Biotechnol. 75, 291–295. doi: 10.1016/S0168-1656(99)00163-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, C., Jiang, Q., Wang, F., Wang, G. L., Xu, Z. S., and 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. doi: 10.1371/journal.pone.0117569

PubMed Abstract | CrossRef Full Text | Google Scholar

Udvardi, M. K., Czechowski, T., and Scheible, W. (2008). Eleven golden rules of quantitative RT-PCR. Plant Cell 20, 1736–1737. doi: 10.1105/tpc.108.061143

PubMed Abstract | CrossRef Full Text | Google Scholar

Vandesompele, J., Preter, K. D., Pattyn, F., Poppe, B., Roy, N. V., Paepe, A. D., et al. (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3:research0034. doi: 10.1186/gb-2002-3-7-research0034

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanguilder, H. D., Vrana, K. E., and Freeman, W. M. (2008). Twenty-five years of quantitative PCR for gene expression analysis. BioTechniques 44, 619–626. doi: 10.2144/000112776

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, G. L., Tian, C., Jiang, Q., Xu, Z. S., Wang, F., and 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. doi: 10.1080/14620316.2016.1148372

CrossRef Full Text | Google Scholar

Warrington, J. A., Nair, A., Mahadevappa, M., and Tsyganskaya, M. (2000). Comparison of human adult and fetal expression and identification of 535 housekeeping/maintenance genes. Physiol. Genomics 2, 143–147.

PubMed Abstract | Google Scholar

Wong, M. L., and Medrano, J. F. (2005). Real-time PCR for mRNA quantitation. BioTechniques 39, 75–85. doi: 10.2144/05391RV01

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Z. J., Li, X. H., Liu, Z. W., Li, H., Wang, Y. X., and 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. doi: 10.1007/s10142-015-0457-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, F., Peng, X., Chen, D., Xu, L., and Zhang, B. (2012). miRDeepfinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol. Biol. 80, 75–84. doi: 10.1007/s11103-012-9885-2

CrossRef Full Text | Google Scholar

Xu, H., Bao, J. D., Dai, J. S., Li, Y., and 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. doi: 10.1371/journal.pone.0142015

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Z. S., Chen, M., Li, L. C., and Ma, Y. Z. (2011). Functions and application of the AP2/ERF transcription factor family in crop improvement. J. Integr. Plant Biol. 53, 570–585. doi: 10.1111/j.1744-7909.2011.01062.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, X., Zhang, F., Tao, Y., Song, S., and 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. doi: 10.1016/j.scienta.2014.10.027

CrossRef Full Text | Google Scholar

Zhang, H., Chen, F., Wang, X., and 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. doi: 10.1016/j.foodres.2006.03.007

CrossRef Full Text | Google Scholar

Zhuang, J., Cai, B., Peng, R. H., Zhu, B., Jin, X. F., Xue, Y., et al. (2008). Genome-wide analysis of the AP2/ERF gene family in Populus trichocarpa. Biochem. Biophys. Res. Commun. 371, 468–474. doi: 10.1016/j.bbrc.2008.04.087

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhuang, J., Zhang, J., Hou, X. L., Wang, F., and 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. doi: 10.1080/07352689.2014.870420

CrossRef Full Text | Google Scholar