Mengmeng Duan1,2
Jinglei Wang1
Xiaohui Zhang1Haohui Yang1Haiping Wang1Yang Qiu1Jiangping Song1
Yangdong Guo2
Xixiang Li1*- 1Key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Ministry of Agriculture, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, China
- 2Department of Vegetable Sciences, College of Horticulture, China Agricultural University, Beijing, China
Radish (Raphanus sativus) is an important cruciferous root crop with a close relationship to Chinese cabbage (Brassica rapa). RT-qPCR is used extensively to evaluate the expression levels of target genes, and accurate measurement of target gene expression with this method is determined by the valid reference genes used for data nomalization in different experimental conditions. Screening for appropriate reference genes with stable expression based on RT-qPCR data is important for gene expression and functional analysis research in radish and its relatives. However, many researches have thought that almost no single reference gene is widely suitable for all experimental conditions, and few researchers have paid attention to the validation of reference genes in radish gene expression analysis. In the present study, 12 candidate reference genes were selected for analysis. Their expression in 28 samples, including 20 radish samples from different organs and conditions, four Chinese cabbage organs and four organs of their distant hybrid, was assessed by RT-qPCR and then five software tools—ΔCt, geNorm, NormFinder, BestKeeper and RefFinder—were used to compare their expression stability. The results showed that the most suitable reference genes were different in different organs and conditions. GAPDH, DSS1, and UP2 were optimal reference genes for gene expression analysis in all organs and conditions in radish. UPR, GSNOR1, and ACTIN2/7 were the most stable reference genes in different radish organs. UP2 and GAPDH were suitable reference genes for radish pistil development studies. RPII, UBC9, and GAPDH had the most stable expression in radish under various stresses. DSS1, UP2, and TEF2 were the optimal reference genes for Chinese cabbage organs, whereas TUA was optimal for the distant hybrid. UP2, and TEF2 were appropriate reference genes for all of the samples together. The optimal reference genes we identified, UP2, GAPDH, UPR, and GSNOR1 were verified by normalizing the expression patterns of YAB3, RPL, and FUL. These results will provide important information for selecting target reference genes in different research contexts and improve the accuracy and precision of gene expression analysis for radish, Chinese cabbage and their distant hybrid.
Introduction
RT-qPCR (Quantitative Real-time PCR) is an important and effective method to evaluate the expression of target genes in different tissues, organs and conditions (Bustin, 2000, 2002; Gachon et al., 2004). Although RT-qPCR has many advantages including rapidity, sensitivity, and specificity, the expression data are affected by experimental conditions or inherent technical variations as well as true biological variation (Bustin, 2002; Derveaux et al., 2010). Stably expressed reference genes are usually used to normalize transcriptome quantification through exposure to the same preparation processes as the target genes. Therefore, appropriate reference genes for data normalization are critical to obtain accurate expression data by RT-qPCR (Ginzinger, 2002; Gachon et al., 2004; Dheda et al., 2005; Guenin et al., 2009; Schmidt and Delaney, 2010). Otherwise, inappropriate reference genes will lead to errors in the expression data for the target gene (Dheda et al., 2005). An appropriate reference gene should be expressed stably in different types of cells, tissues and organs, and at different developmental stages. At the same time, its expression should be high or moderate and assumed to be unaffected under different experimental conditions (Bustin, 2002; Brunner et al., 2004; Czechowski et al., 2005; Gutierrez et al., 2008).
The most commonly used reference genes for the normalization of RT-qPCR data in plants include glyceraldehyde-3-phosphate dehydrogenase (GAPDH), actin (ACT), ubiquitin (UBI), 18S ribosomal RNA (18S), α-tubulin and β-tubulin (TUA and TUB, respectively), and the ubiquitin-conjugating enzyme (UBC), which function in maintaining cell survival irrespective of physiological conditions (Bustin, 2002; Brunner et al., 2004; Radonić et al., 2004; Czechowski et al., 2005; Jain et al., 2006; Gutierrez et al., 2008). These reference genes were identified by Northern-blotting in the pre-genomic era and were assumed to have stable expression at various developmental stages and under various experimental conditions (Jain et al., 2006). Unfortunately, several studies have proved that the transcript levels of these traditional reference genes vary under different experimental conditions (Thellin et al., 1999; Suzuki et al., 2000; Czechowski et al., 2005; Zhu et al., 2011). Thus, many validation methods have been used to confirm the variability of conventional reference genes though systematic research in various plants, such as Arabidopsis thaliana (Han et al., 2013), potato (Castro-Quezada et al., 2013), tomato (Expósitorodríguez et al., 2008), cabbage (Chen et al., 2010), soybean (Bo et al., 2008), tobacco (Schmidt and Delaney, 2010), watermelon (Kong et al., 2014), melon (Kong et al., 2016), pearl millet (Shivhare and Lata, 2016), and celery (Li et al., 2016b). These experiments show that no gene has universal expression in all tissues types or under all experiment conditions. Thus, one or more different reference genes need to be selected according to the specific set of biological samples being studied.
With the development of molecular biological technologies such as the Affymetrix GeneChip, microarrays and high-throughput sequencing technologies in recent years, a large number of novel reference genes have been developed in A. thaliana (Czechowski et al., 2005) and crops such as Brassica juncea (Qi et al., 2010), Brassica napus (Yang et al., 2014), papaya (Zhu et al., 2011) and rice (Jain, 2009). In one study, an F-box protein (F-box), a SAND family protein and the mitosis protein YLS8 were expressed more stably than the traditional reference genes ACTIN-2, elongation-factor-1-a (EF1-α) and ubiquitin-conjugating enzyme 10 (UBC10) in A. thaliana (Remans et al., 2008). UBC9 and UP2 were chosen as reference genes on the basis of microarray data described by Schmid et al. in tobacco, and have also shown more stable expression than some traditional reference genes in other crops (Czechowski et al., 2005; Schmid et al., 2005; Kwon et al., 2009; Kumar et al., 2011). Twelve novel reference genes obtained from genomic and transcriptomic data were identified to perform better than the traditional reference genes ACTIN7 and GAPDH in B. napus (Yang et al., 2014). Additionally, the most commonly used tools to assess the stability of reference genes in a set of samples are the delta cycle threshold (ΔCt) method (Silver et al., 2006) and the software tools geNorm (Vandesompele et al., 2002), NormFinder (Andersen et al., 2004), BestKeeper (Pfaffl et al., 2004), and RefFinder (Xie et al., 2012). The integrated application of these software tools has improved the accuracy of candidate reference gene identification (Hao et al., 2014; Niu L. et al., 2015).
Radish (Raphanus sativus L.), a member of the Cruciferae family, is a widely cultivated root vegetable across the world. It is also a donor of elite genes for the genetic improvement of other cruciferous crops. Many researchers have worked on the selection of appropriate reference genes in Cruciferae plants, such as Brassica rapa in different tissues and under different abiotic and biotic stresses (Xiao et al., 2012; Wang et al., 2016), B. juncea in different developmental stages and under hormone treatments and drought stress (Chandna et al., 2012), and B. napus in different tissues (Chen et al., 2010). Few researchers have paid attention to the validity of reference genes in radish gene expression analysis, although many studies on gene expression in radish have been carried out in recent years. At present, the most commonly used reference gene in radish is still the traditional gene ACTIN2/7 (Xu et al., 2013; Wang et al., 2014; Yu et al., 2015) though microarray and next-generation sequencing (NGS) have been applied in radish. Therefore, it is urgent to identify optimal reference genes for RT-qPCR normalization in radish, and even its distant hybrids (such as RRAA, n = 38, a hybrid of Chinese cabbage (B. rapa; AA, n = 10) and radish (RR, n = 9), which are important materials for research related to genome evolution, genetic variability, gene exchange and germplasm enhancement in Cruciferae (Lange et al., 1989; Peterka et al., 2004; Lee et al., 2011).
The silique is an important organ for oilseed production and the main reproductive organ in radish. Some research in Arabidopsis suggests that MADS-box genes play an important role in pistil development and molecular interactions (Alvarezbuylla et al., 2010). The establishment of floral organ polarity leads to the expression of YABBY3 (YAB3) on one side of an organ, and YAB3 influences the formation of valve margin tissue. FRUITFUL (FUL) encodes a MADS-domain transcription factor that is expressed in the valves, and has multiple functions during cell differentiation, promoting cell expansion and inhibiting cell division in specific cell types (Gu et al., 1998), and REPLUMLESS (RPL) encodes a BELL-family homeodomain transcription factor. The expression of the valve margin identity genes is limited to the valve margin through negative regulation by FUL in the valves and RPL in the replum (Roeder et al., 2003).
In order to determine suitable reference genes, or a combination of the most stable reference genes for accurate quantification of target genes in radish and its relatives, we selected traditional candidate reference genes from previous studies and the new ones from our preliminary analysis of gene stability based on different sets of transcriptome data for further analysis in the present research. The traditional nine genes widely used in radish and other crops included GAPDH (glyceraldehyde-3-phosphate-dehydrogenase), TEF2 (translation elongation factor 2), ACTIN2/7, TUA (tubulin alpha-5), TUB (tubulin beta-1), RPII (RNA polymerase-II transcription factor), 18S rRNA (18S ribosomal RNA), UBC9 (ubiquitin-conjugating enzyme 9) and UP2 (uncharacterized conserved protein), and three new genes screened out by radish transcriptome analysis were GSNOR (GroES-like zinc-binding dehydrogenase family protein), UPR (uncharacterized protein family), and DSS1 (deletion of SUV3 suppressor 1(I)). RT-qPCR was used to validate the applicability of these genes in different tissues and organs, at different development stages, and under different stress conditions in radish, Chinese cabbage and their distant hybrid. Furthermore, the expression patterns of three key genes, YAB3, RPL, and FUL, which are thought to be related to silique development in Arabidopsis and Brassica plants, were analyzed as a case study to investigate the efficiency of the reference genes. These results will be helpful for the selection of target reference genes to ensure accuracy and precision in gene expression analysis in different research contexts for radish, Chinese cabbage and their distant hybrids.
Materials and Methods
Plant Materials and Treatments
The whole genome-sequenced radish inbred line “36-2,” the Chinese cabbage inbred line “chiifu” and their distant hybrid (RRAA) were used as study materials. Geminating seeds were vernalized at 4°C for about 30 days, sown in plastic pots containing a soil/vermiculite mixture (3:1) and grown in a greenhouse at the research station of the Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, China. Radish seedlings with two true leaves were exposed to the following biotic and abiotic treatments. For TuMV infection treatment, a virus inoculum was prepared by homogenizing infected fresh leaves in a homogenizer with four volumes (w/v) of 0.05 M phosphate buffer (pH 7.0). The inoculation of TuMV into radish was performed as described by Zhang et al. (2007). Under a 12 h diurnal light cycle, the temperature was maintained at 25–28°C in the daytime and 20–22°C at night in an incubator. Leaf samples were obtained on 5th, 10th, and 15th days after inoculation. For pest treatment, diamondback moths were maintained on cabbage (Brassica oleracea) plants in a climate-controlled room at 25°C with a 12:12 h photoperiod and 50–60% relative humidity, and four diamondback moths (second or third instars) were placed on each leaf of the radish seedlings. Leaves exposed to diamondback moths were harvested at 0, 4, 24, and 48 h after infestation (Wei et al., 2013). For cold treatment, radish seedlings were stored at 2–4°C in an incubator under a 12 h diurnal light cycle for 24 h and leaf samples were collected on the 2nd, 4th, and 6th days after exposure to low temperature. For the biotic and abiotic stress treatments, leaf samples collected from two-true-leaf seedlings under normal conditions were used as a control. To examine expression during radish pistil development, flower buds of 5 mm in length and siliques at the 0th, 5th, 15th, and 30th days after pollination were collected at different radish reproductive stages under normal growth conditions. To examine expression in different organs, root, stem and leaf samples of radish plants were collected at the vegetative phase and calyx, petal, stamen and pistil samples of radish, Chinese cabbage and their distant hybrid were collected during the flowering period.
All 28 samples had three biological replicates and each replicate was obtained from three plants. The samples were immediately frozen in liquid nitrogen and stored at −80°C for further use.
Candidate Reference Gene Selection and Primer Design
Seven traditional candidate reference genes (GAPDH, RPII, ACTIN2/7, TEF2, 18S, TUA, and TUB) were selected from previous studies on radish (Xu et al., 2012) and other crops (Jain et al., 2006; Remans et al., 2008; Castro-Quezada et al., 2013; Han et al., 2013). The reference genes UBC9 and UP2, which have been used in some plants (Czechowski et al., 2005; Schmid et al., 2005; Kwon et al., 2009; Kumar et al., 2011) but not in radish, were also selected.
For the selection of new candidate genes, we analyzed publically transcriptomic data from the following 21 different radish tissues and organs at different developmental stages: 7, 14, 20 days root and leaf, 40, 60, 90 days cortical, cambium, xylem, root tip and leaf (Mitsui et al., 2015; http://www.nodai-genome-d.org/download.html). A candidate reference gene was defined as a gene with the most constant expression level, i.e., a gene with a small coefficients of variation (CVs) (De Jonge et al., 2007). Therefore, the raw RNA-seq data were used to calculate the mean expression value, standard deviation, and CVs according to the following formula firstly: CVs = standard deviation of RPKM/average of RPKM. Besides the stability of gene expression level, the expression intensity of candidate reference genes is also significant. Genes with high or lowly expression abundance are not appropriate for being the reference genes (Xu et al., 2015). Hence, these selected new genes had to meet the following requirements: CV ≤ 30%, 100 ≤ RPKM ≤ 500.
The coding DNA sequences (CDS) and DNA sequences of GSNOR1, UPR, DSS1, UBC9, and UP2 were obtained from previously reported radish genomic data (Mitsui et al., 2015) by homology analysis. Primers for these genes were designed using Primer 5.0 for RT-qPCR. The product sizes were set in the range of 80–200 bp. At least one primer in each pair spanned the exon-intron junction to avoid amplification of gDNA in possibly contaminated samples. A single band of expected size in 1% agarose gel electrophoresis and a single peak in RT-qPCR melting curve were used as criteria to ensure the specificity of amplification for every candidate reference gene. In addition, primers for the seven traditional reference genes (GAPDH, RPII, ACTIN2/7, TEF2, 18S, TUA, and TUB) were used as previously reported (Xu et al., 2012).
Total RNA Extraction and cDNA Synthesis
Total RNA was extracted from all samples using an RNAprep pure Plant Kit (TransGen, Beijing, China) according to the manufacturer's instructions. Genomic DNA was eliminated from the total RNA using RNase-free DNase I. The integrity of the RNA was checked by electrophoresis in a 2% agarose gel. The quantity and purity of the RNA were evaluated using a NanoDrop™ 2000 spectrophotometer (Thermo Scientific). Samples with A260/A280 > 1.8 and A260/A230 < 2.0 were used for subsequent cDNA synthesis. First-strand cDNA was synthesized with TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (TransGen, Beijing, China) according to the manual, and oligo dT was used for the cDNA synthesis. For each sample, 1 μg of total RNA was used for every 20 μL of the reverse transcription reaction system.
Real-Time Quantitative RT-qPCR
RT-qPCR was carried out on an ABI StepOne Real-Time PCR System using TranStart Top Green PCR Super Mix (TransGen). RT-qPCR was carried out on a LightCycler480 System. Reactions were performed in a total volume of 20 μL, which contained 2 μL cDNA template, 0.4 μL each primer, 10 μL 2× Top Green RT-qPCR SuperMix and 6.8 μL ddH2O. The PCR cycling conditions were as follows: 94°C for 30 s and 40 cycles of 95°C for 5 s, 55°C for 15 s and 72°C for 10 s. Melting curve analysis was performed after 40 cycles to test the primer specificity by heating from 65° to 95°C with a stepwise increase of 0.5°C every 10 s. Three technical replicates were used for each sample. Controls without a template were also included. For each gene, the full sample set in each replication was run on the same plate to exclude any technical variation. The amplification efficiencies for all primer pairs were evaluated using fivefold dilutions of the pooled cDNA (1/5, 1/25, 1/125, 1/625, and 1/3125). EASY dilution solution (Takara, Japan) was used for primer dilution.
Data Analysis
The Ct value of each reference gene was used to evaluate its expression level. The amplification efficiencies of the candidate reference genes were calculated according to the following formula: E (%) = (10−1/slope − 1) × 100; the slope was the standard curve of Ct values of the five gradients for each reference gene. The expression stability was evaluated using the ΔCt method (Silver et al., 2006), geNorm (Vandesompele et al., 2002), NormFinder (Andersen et al., 2004), and BestKeeper (Pfaffl et al., 2004), and then comprehensively analyzed using RefFinder (available online: http://omictools.com/reffinders2857.html) (Xie et al., 2012). The ΔCt method was employed to rank the stability of the candidate reference genes by calculating the standard deviation (SD). The gene with the lowest SD was identified as the most stable reference gene (Silver et al., 2006). geNorm makes decisions based on the principle that the expression levels of two appropriate reference genes have a stable ratio across the investigated samples. The expression stability (M) of each candidate gene in geNorm is the variation of the given gene compared with the all other candidate reference genes. The candidate reference genes were ranked based on their M values, and the least stable gene with the highest M value was excluded stepwise. Finally, two genes remained, which were the most stable reference genes. To determine the optimal number of reference genes required for normalization, normalization factors (NFs) were calculated using the geNorm software, with a cut-off value for Vn/n + 1 of 0.15, below which no additional reference gene is required for normalization (Vandesompele et al., 2002). NormFinder ranks candidate reference genes using an ANOVA-based model, which takes inter-group and intra-group relationships into consideration (Andersen et al., 2004). BestKeeper determines the expression stability of candidate reference genes based on SD and CV calculations for Ct values. The most stable reference gene has the lowest SD and CV values (Pfaffl et al., 2004). RefFinder was used to generate a comprehensive ranking based on the geometric mean of the three programs mentioned above (Xie et al., 2012).
Normalization of Gene Expression in Pistil Development and Verification of the Screened-Out Reference Genes
To validate the selected reference genes, the relative expression levels of three transcription factors, YAB3, FUL, and RPL, were analyzed during pistil development in radish. The gene sequences were obtained from the genome sequence data of radish (http://www.nodai-genome-d.org/download.html), and primers for the three genes were designed according to the aforementioned methods (primer sequences are listed in Table 1). The specificity of the primers was confirmed by electrophoresis in 2% agarose gels and melting curve analysis. Samples of siliques at different developmental stages were collected as described above. The relative expression level was evaluated using 2−ΔΔCt method. The best reference genes identified by RefFinder (UP2, GAPDH, and UPR) were used for normalization. The least stable reference gene identified by RefFinder (GSNOR1) was also used for normalization. Three biological and three technical replicates were used for the measurements at each sampling point.
Table 1. Candidate reference genes and primers used for RT-qPCR in radish, Chinese cabbage, and their distant hybrid.
Results
Verification of Primer Efficiency for the Candidate Reference Genes
Based on these selection procedures for the transcriptome sequencing data, 3 genes that had a minor variation in expression were selected (Table S1). Nine traditional reference genes (GAPDH, RPII, ACTIN2/7, TEF2, 18S, TUA, TUB, UBC9, and UP2) from previous studies and three new reference genes (GSNOR1, UPR, and DSS1) from our analysis were selected as candidate reference genes in our study. Using primers designed for each gene, the PCR amplification specificities of the 12 candidate genes were checked by 2% agarose gel electrophoresis with pistil cDNA samples from radish, Chinese cabbage and their distant hybrid. The results showed that all 12 candidate genes had specific amplification and the product lengths were consistent with the expected lengths in radish, Chinese cabbage and their distant hybrid (Figure S1). Additionally, melting curve analysis of the candidate genes showed a single peak for each gene (Figures S2, S3), which confirmed that the primers for these genes were specific in radish, Chinese cabbage and the distant hybrid. The amplification efficiencies of the 12 candidate genes in different samples varied from 0.86 (UP2 in Chinese cabbage) to 1.12 (RPII in radish) (Table 1, Figures S4, S5).
Expression Stability Analysis Based on the Expression Profiles of the Candidate Reference Genes
The expression levels of the 12 candidate reference genes were evaluated in 28 samples collected from different organs in radish, Chinese cabbage and their distant hybrid under abiotic and biotic stress conditions using Ct values. As shown in the boxplot (Figure 1), the Ct values of the 12 candidate reference genes varied from 9.4 (18S) to 31.4 (UPR). UP2 had the highest mean Ct value (25.49) with lowest expression abundance among these genes. By contrast, 18S had the lowest mean value (15.6) with the highest expression abundance. In addition, the expression variation among the 28 samples for each candidate reference gene ranged from 6.61 (TUA) to 12.12 (18S). These results showed that no candidate reference gene had stable expression under all conditions, and that it is necessary to identify appropriate reference genes for precise normalized expression under specific conditions in R. sativus.
Figure 1. Boxplot analysis of the expression profiles of 12 candidate reference genes across all 28 samples. The line across the box represents the median. The boxes represent the 25/75 percentiles. The whiskers show the maximum and minimum values. The circles indicate outliers.
The 12 candidate reference genes were then subjected to further analysis based on seven sets of samples: radish organs under normal conditions (RsOs), radish pistils at different developmental stages (RsPs), radish seedlings under biotic (diamondback moth) and abiotic (cold, TuMV) stresses (RsSs), all radish samples (RsAll), Chinese cabbage organs (BrOs), organs of the distant hybrid (HOs) and all samples (All). The ΔCt method, geNorm, NormFinder, BestKeeper, and RefFinder were used to evaluate the expression stability of the candidate reference genes.
In radish organs under normal conditions, GSNOR1 and UPR were the most stable reference genes with M values of 0.34, and ACTIN2/7 was the next most stable gene with an M value of 0.473 in geNorm analysis. In addition, the recommended number of reference genes for expression normalization was three, with a V3/4 value of 0.119. Therefore, GSNOR1, UPR, and ACTIN2/7 were identified as the optimal reference genes in geNorm analysis. UPR was also the most stable reference gene in NormFinder, ΔCt and BestKeeper analysis. The comprehensive ranking order suggested that UPR, GSNOR1, and ACTIN2/7 were the optimal reference genes and 18S was the least stable reference gene in radish organs (Table 2, Table S2).
Table 2. Stability rankings of candidate reference genes in seven sets by ΔCt, BestKeeper, NormFinder, geNorm, and RefFinder.
During radish pistil development, DSS1 and GAPDH were identified as the most stable reference genes with M values of 0.39, and three reference genes were recommended for expression normalization according to the V4/5 value (0.111) by geNorm analysis. UP2 was the most stable reference gene in NormFinder and ΔCt analysis. In addition, TUA was the most stable gene in BestKeeper analysis. The comprehensive ranking indicated that UP2 and GAPDH were the most appropriate reference genes and GSNOR1 was the least stable reference gene during radish pistil development (Table 2, Table S2).
Under biotic and abiotic stress conditions, UP2 and UBC9 were the most stable reference genes with M values of 0.69, and four reference genes were recommended for expression normalization in geNorm analysis. The most stable reference gene was RPII according to NormFinder and ΔCt analysis. In addition, GSNOR1 was the most stable reference gene in BestKeeper analysis. According to the calculations performed by RefFinder, RPII, UBC9, and GAPDH were the three most stable reference genes in radish samples under biotic and abiotic stresses, whereas 18S was the least stable reference gene in all of the above evaluation systems (Table 2, Table S2).
Across all of the radish samples, DSS1 and GAPDH were the most stable reference genes, and five reference genes were recommended for expression normalization by geNorm analysis. GAPDH was also the best reference gene in ΔCt analysis. However, TEF2 was the best reference gene in NormFinder analysis and UPR was recommended by BestKeeper. The comprehensive ranking by RefFinder showed that GAPDH, DSS1, and UP2 were the most stably expressed reference genes across all radish samples, while 18S was the least stable reference gene (Table 2, Table S2).
In Chinese cabbage organs, UP2 and GAPDH were considered the best reference genes by geNorm analysis, and the appropriate number of genes for normalization was four, with a V4/5 value of 0.126. In contrast, NormFinder recognized ACTIN2/7 as the most stable reference gene. DSS1 and UPR were the most stable reference genes in ΔCt and BestKeeper analysis, respectively. DSS1, UP2, and TEF2 were ranked highly and 18S was ranked last in the comprehensive analysis by RefFinder (Table 2, Table S2).
In organs of the distant hybrid, GSNOR1 and ACTIN2/7 were found to be the best reference genes with M values of 0.26, and four reference genes were recommended for normalization because the pairwise value of V6/7 was 0.142. Conversely, TUA, DSS1, and RPII were the most stable reference genes in ΔCt, BestKeeper, and NormFinder analysis, respectively. In the comprehensive analysis, TUA was determined to be the most stable reference gene in organs of the distant hybrid, and 18S was identified as the least stable gene (Table 2, Table S2).
Across all 28 samples, UP2 and GAPDH were identified as the best reference genes with the lowest M values (0.75), whereas TUB was identified as the worst reference gene with the highest M value (1.53) in the geNorm analysis (Table 2, Table S2). The results in Figure 2 show that the V8/9 and V9/10 values were 0.152 and 0.14, respectively, which suggests that nine genes are required for reliable normalization when all samples are considered. GAPDH was also the best reference gene with the lowest stability value (0.35) according to NormFinder analysis, while TUB had the highest stability value (1.43) and was also the least stable reference gene. UP2 was the most stable reference gene in ΔCt analysis. UPR was the most stable reference gene in BestKeeper analysis. In the comprehensive analysis, UP2, GAPDH, and TEF2 were the recommended reference genes, and 18S was the least stable gene in all five ranking lists.
Figure 2. Determination of the optimal number of reference genes. Pairwise variation (Vn/n + 1) was calculated by geNorm to determine the number of reference genes required for accurate normalization in different sample sets. The dashed lines indicate 0.15 was the cut-off value to determine the optimal number of reference genes for gene normalization. RsOs, radish organs under normal conditions; RsPs, radish pistils at different developmental stages; RsSs, radish seedlings under biotic (diamondback moth) and abiotic (cold, TuMV) stresses; RsAll, all radish samples; BrOs, Chinese cabbage organs; HOs, organs of the distant hybrid; All, all samples.
Validation of Candidate Reference Genes
YAB3, FUL, and RPL, three important transcription factors during pistil development, were used as examples for expression analysis with UP2, GAPDH, UPR, and GSNOR1 as reference genes for normalization. UP2, GAPDH, and UPR were the top-ranked genes in RefFinder analysis, and were suggested for accurate expression normalization during pistil development in radish. GSNOR1 was ranked at the bottom by RefFinder analysis in radish during pistil development (Table 2, Table S2).
The expression levels of YABB3, FUL, and RPL showed similar change patterns when the stable reference genes UP2, GAPDH, and UPR were used for normalization. In contrast, the expression profiles of YAB3, FUL, and RPL were distorted when GSNOR1 was used for normalization. The relative expression levels of YAB3 and FUL at 0 and 15 days were overestimated, and the expression profile of RPL at the first stage was abnormally upregulated when GSNOR1 was used as the reference gene (Figure 3).
Figure 3. Relative expression levels of YAB3 (A), FUL (B), and RPL (C) in radish pistil development. The three top-ranked genes (UP2, GAPDH, and UPR) and the last-ranked genes GSNOR1 by RefFinder were used for expression normalization. Error bars indicate the standard deviation of three biological replicates.
Discussion
RT-qPCR is an essential method for obtaining expression profiles and understanding the biological functions of target genes (Bustin, 2000). The accuracy of RT-qPCR results is mainly dependent on a suitable normalization strategy that selects appropriate reference genes (Vanguilder et al., 2008; Derveaux et al., 2010). Ideally, reference genes should have stable expression in different developmental stages and different tissues or organs, as well as under different experimental conditions, and the expression level should be similar to that of the target gene (Thellin et al., 1999; Suzuki et al., 2000). Reference genes involved in cytoskeleton structure (ACTIN2/7, TUA, TUB, and 18S), protein synthesis (RPII and TEF2) and biological metabolic processes (GAPDH) are often used as reference genes (Chen et al., 2010; Rebouças et al., 2013). Using genomic and transcriptomic data, novel reference genes such as UBC9 and UP2 have also been selected for gene normalization (Schmid et al., 2005; Kwon et al., 2009). In our research, three new genes (GSNOR, UPR, and DSS1) were selected from radish transcriptome data as candidate reference genes in addition to the seven classical reference genes mentioned above.
However, no candidate reference gene had invariable expression across all samples in our research, which highlights the importance of employing suitable reference genes for particular crops under specific conditions using statistical approaches. In previous studies, the expression stability of reference genes also varied in different species, genotypes, developmental stages, organs, tissues and experimental conditions (Livak and Schmittgen, 2001; Jin et al., 2013). In Chinese cabbage, Qi et al. (2010) found that EF-1-α was the best reference gene in a given set of tissues. However, another study showed that UP1 and UBC9 were the best choices for vegetative tissues of Chinese cabbage (Chen et al., 2010). In radish, TEF2 and RPII performed well in a range of different tissue types (Xu et al., 2012). However, in our research, UPR, UP2, and GAPDH were ranked first in different organs, pistils at different developmental stages and all radish samples, respectively, while RPII was the most suitable reference gene only under biotic and abiotic stresses. Moreover, 18S rRNA, which is a component of the small subunit of eukaryotic ribosomes (40S), has been considered the most stable gene under various treatment conditions (Jain et al., 2006; Niu X. et al., 2015). However, 18S was the most unstable reference gene in all samples in our research, which was similar to research by Xu et al. (2012) and in Pisum sativum (Die et al., 2010). The TUB gene, which participates in cell structural maintenance, was the optimum reference gene under excess salinity and drought conditions in kenaf (Niu X. et al., 2015), but had poor performance in our research, which was similar to results in potato (Reid et al., 2006) and soybean (Hu et al., 2009). Therefore, it is necessary to select appropriate reference genes for normalization according to different species and experimental conditions.
The four computational methods (ΔCt, BestKeeper, NormFinder and geNorm) (Niu L. et al., 2015; He et al., 2016; Li et al., 2016a) for expression stability are based on different algorithms and analytical procedures (Vandesompele et al., 2002; Andersen et al., 2004; Pfaffl et al., 2004; Silver et al., 2006). In our research, we found that the most unstable genes identified by the four algorithms were mostly the same, but the most stable genes were not consistent. This was possibly due to differences among the algorithms. Integrated analysis using different programs can minimize errors in the stability evaluation of candidate reference genes (Li et al., 2016a; Shivhare and Lata, 2016). Here, RefFinder (Xie et al., 2012) was used to combine the results of the four computational methods and generate a comprehensive ranking list for the candidate reference genes. In the comprehensive analysis with RefFinder, the seven traditional reference genes, GAPDH, ACTIN2/7, UP2, RPII, UBC9, TEF2, and TUA, and three new reference genes, GSNOR1, UPR, and DSS1, were found to be the optimal reference genes for gene expression analysis of different species, organs and conditions in the present study.
GAPDH is an abundant glycolytic enzyme present in most cell types (Giulietti et al., 2001) and has been extensively used as a reference gene in RT-qPCR experiment. For example, it was recommended for measuring the expression of genes of interest in diverse tissues and genotypes of sugarcane (Iskandar et al., 2004). Furthermore, Mamo indicated that GAPDH had consistent changes at different stages of embryonic period (Mamo et al., 2007). Medrano found that GAPDH was one of the best genes in combination for embryonic samples (Medrano et al., 2017). Our results were consistent with these previous studies suggesting that GAPDH can also be universally used as a reference gene in radish and their relatives, although several reports in animal and plant systems have suggested that GAPDH has some limitations as an internal control gene (Sirover, 1999; Nazari et al., 2015). The reference gene UP2, encoding an uncharacterized conserved protein, is related to cell structure. Its stability has been verified in many studies (Czechowski et al., 2005; Schmid et al., 2005; Pollier et al., 2014), and in the present research, it also had widely expression stability in all samples. RPII was the most stable reference gene under biotic and abiotic stresses and also performed as well as TUA in the distant hybrid. In addition, DSS1 had the least variability across all radish samples and in Chinese cabbage organs. DSS1 participates in the deletion of SUV3 suppressor 1 and is a 26S proteasome ubiquitin receptor that binds ubiquitin chains (Paraskevopoulos et al., 2014), which may explain why it had stable expression in various conditions. GSNOR1 is a member of the GroES-like zinc-binding dehydrogenase family and was a suitable reference gene across radish organs. By comparison, UPR is an uncharacterized protein and was the best fit for research on radish organs. Accordingly, these reference genes can be selected for gene normalization in different species, tissues and experimental conditions.
Intergenetic and interspecific hybridization are important ways of enriching genetic backgrounds and developing new crop cultivars with novel traits for Cruciferous vegetables (Gueritaine et al., 2002; Darmency et al., 2010; Tonosaki et al., 2013). Radish, Chinese cabbage and their distant hybrid were used as samples for the first time to evaluate the stability of candidate reference genes in this study. Our identification of optimal reference genes will make an important contribution to comparative analyses of target gene expression in these species and their interspecific hybrids.
Author Contributions
Conceived and designed the experiments: XL and MD; Performed the experiments: MD and XZ; Analyzed the data: MD, JW, and YG; Prepared materials: MD, HY, HW, YQ, and JS; Wrote the paper: MD and XL. All of the authors read and approved the final manuscript.
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
This work was supported by grants from the national key research and development plan of the Ministry of Science and Technology of China (2016YFD0100204-02, 2013BAD01B04-1) and the Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences (CAAS-XTCX2016017, CAAS-XTCX2016016-4-4, CAAS-XTCX2016001-5-3).
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.01605/full#supplementary-material
References
Alvarezbuylla, E. R., Benítez, M., Corverapoiré, A., Chaos Cador, Á., De Folter, S., Alicia, G. D. B., et al. (2010). Flower development. Arabidopsis Book 8:e0127. doi: 10.1199/tab.0127
Andersen, C. L., Jensen, J. L., and Orntoft, 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
Bo, J., Liu, B., Bi, Y., Hou, W., Wu, C., and Han, T. (2008). Validation of internal control for gene expression study in soybean by quantitative real-time PCR. BMC Mol. Biol. 9:59. doi: 10.1186/1471-2199-9-59
Brunner, A. M., Yakovlev, I. A., and Strauss, S. H. (2004). Validating internal controls for quantitative plant gene expression studies. BMC Plant Biol. 4:14. doi: 10.1186/1471-2229-4-14
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
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
Castro-Quezada, P., Aarrouf, J., Claverie, M., Favery, B., Mugniéry, D., Lefebvre, V., et al. (2013). Identification of reference genes for normalizing RNA expression in potato roots infected with cyst nematodes. Plant Mol. Biol. Rep. 31, 936–945. doi: 10.1007/s11105-013-0566-3
Chandna, R., Augustine, R., and Bisht, N. C. (2012). Evaluation of candidate reference genes for gene expression normalization in Brassica juncea using real time quantitative RT-PCR. PLoS ONE 7:e36918. doi: 10.1371/journal.pone.0036918
Chen, X., Truksa, M., Shah, S., and Weselake, R. J. (2010). A survey of quantitative real-time polymerase chain reaction internal reference genes for expression studies in Brassica napus. Anal. Biochem. 405, 138–140. doi: 10.1016/j.ab.2010.05.032
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
Darmency, H., Lefol, E., and Fleury, A. (2010). Spontaneous hybridizations between oilseed rape and wild radish. Mol. Ecol. 7, 1467–1473. doi: 10.1046/j.1365-294x.1998.00464.x
De Jonge, H. J. M., Fehrmann, R. S.N., De Bont, E. S. J. M., Hofstra, R. M. W., Frans, G., Kamps, W. A., et al. (2007). Evidence based selection of housekeeping genes. PLoS ONE 2:e898. doi: 10.1371/journal.pone.0000898
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
Dheda, K., Huggett, J. F., Chang, J. S., Kim, L. U., Bustin, S. A., Johnson, M. A., et al. (2005). The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Anal. Biochem. 344, 141–143. doi: 10.1016/j.ab.2005.05.022
Die, J. V., Román, B., Nadal, S., and González-Verdejo, C. I. (2010). Evaluation of candidate reference genes for expression studies in Pisum sativum under different experimental conditions. Planta 232, 145–153. doi: 10.1007/s00425-010-1158-1
Expósitorodríguez, M., Borges, A. A., Borgespérez, A., and Pérez, J. A. (2008). Selection of internal control genes for quantitative real-time RT-PCR studies during tomato development process. BMC Plant Biol. 8:131. doi: 10.1186/1471-2229-8-131
Gachon, C., Mingam, A., and Charrier, B. (2004). Real-time PCR: what relevance to plant studies? J. Exp. Bot. 55, 1445–1454. doi: 10.1093/jxb/erh181
Ginzinger, D. G. (2002). Gene quantification using real-time quantitative PCR: an emerging technology hits the mainstream. Exp. Hematol. 30, 503–512. doi: 10.1016/S0301-472X(02)00806-8
Giulietti, A., Overbergh, L., Valckx, D., Decallonne, B., Bouillon, R., and Mathieu, C. (2001). An overview of real-time quantitative PCR: applications to quantify cytokine gene expression. Methods 25:386. doi: 10.1006/meth.2001.1261
Gu, Q., Ferrandiz, C., Yanofsky, M. F., and Martienssen, R. (1998). The FRUITFULL MADS-box gene mediates cell differentiation during Arabidopsis fruit development. Development 125, 1509–1517.
Guenin, S., Mauriat, M., Pelloux, J., Van Wuytswinkel, O., Bellini, C., and Gutierrez, L. (2009). Normalization of qRT-PCR data: the necessity of adopting a systematic, experimental conditions-specific, validation of references. J. Exp. Bot. 60, 487–493. doi: 10.1093/jxb/ern305
Gueritaine, G., Sester, M., Eber, F., Chevre, A. M., and Darmency, H. (2002). Fitness of backcross six of hybrids between transgenic oilseed rape (Brassica napus) and wild radish (Raphanus raphanistrum). Mol. Ecol. 11, 1419–1426. doi: 10.1046/j.1365-294X.2002.01540.x
Gutierrez, L., Mauriat, M., Pelloux, J., Bellini, C., and Van Wuytswinkel, O. (2008). Towards a systematic validation of references in real-time rt-PCR. Plant Cell 20, 1734–1735. doi: 10.1105/tpc.108.059774
Han, B., Yang, Z., Samma, M. K., Wang, R., and Shen, W. (2013). Systematic validation of candidate reference genes for qRT-PCR normalization under iron deficiency in Arabidopsis. Biometals 26, 403–413. doi: 10.1007/s10534-013-9623-5
Hao, X., Horvath, D. P., Chao, W. S., Yang, Y., Wang, X., and Xiao, B. (2014). Identification and evaluation of reliable reference genes for quantitative real-time PCR analysis in tea plant (Camellia sinensis (L.) O. Kuntze). Int. J. Mol. Sci. 15:22155. doi: 10.3390/ijms151222155
He, Y., Yan, H., Hua, W., Huang, Y., and Wang, Z. (2016). Selection and validation of reference genes for quantitative real-time PCR in Gentiana macrophylla. Front. Plant Sci. 7:945. doi: 10.3389/fpls.2016.00945
Hu, R., Fan, C., Li, H., Zhang, Q., and Fu, Y. F. (2009). Evaluation of putative reference genes for gene expression normalization in soybean by quantitative real-time RT-PCR. BMC Mol. Biol. 10:93. doi: 10.1186/1471-2199-10-93
Iskandar, H. M., Simpson, R. S., Casu, R. E., Bonnett, G. D., Maclean, D. J., and Manners, J. M. (2004). Comparison of reference genes for quantitative real-time polymerase chain reaction analysis of gene expression in sugarcane. Plant Mol. Biol. Rep. 22, 325–337. doi: 10.1007/BF02772676
Jain, M. (2009). Genome-wide identification of novel internal control genes for normalization of gene expression during various stages of development in rice. Plant Sci. 176, 702–706. doi: 10.1016/j.plantsci.2009.02.001
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. Biophy. Res. Commun. 345, 646–651. doi: 10.1016/j.bbrc.2006.04.140
Jin, X., Fu, J., Dai, S., Sun, Y., and Hong, Y. (2013). Reference gene selection for qPCR analysis in cineraria developing flowers. Sci. Hortic-amsterd. 153, 64–70. doi: 10.1016/j.scienta.2013.01.023
Kong, Q., Gao, L., Cao, L., Liu, Y., Saba, H., Huang, Y., et al. (2016). Assessment of suitable reference genes for quantitative gene expression studies in melon fruits. Front. Plant Sci. 7:1178. doi: 10.3389/fpls.2016.01178
Kong, Q., Yuan, J., Gao, L., Zhao, S., Jiang, W., Huang, Y., et al. (2014). Identification of suitable reference genes for gene expression normalization in qRT-PCR analysis in watermelon. PLoS ONE 9:e90612. doi: 10.1371/journal.pone.0090612
Kumar, V., Sharma, R., Trivedi, P. C., Vyas, G. K., and Khandelwal, V. (2011). Traditional and novel references towards systematic normalization of qRT-PCR data in plants. Aust. J. Crop. Sci. 5, 1455–1468.
Kwon, M. J., Oh, E., Lee, S., Roh, M. R., Kim, S. E., Lee, Y., et al. (2009). Identification of novel reference genes using multiplatform expression data and their validation for quantitative gene expression analysis. PLoS ONE 4:e6162. doi: 10.1371/annotation/695436c7-3329-4bdc-9832-f427ecc33698
Lange, W., Toxopeus, H., Lubberts, J. H., Dolstra, O., and Harrewijn, J. L. (1989). The development of Raparadish (x Brassicoraphanus, 2n = 38), a new crop in agriculture. Euphytica 40, 1–14.
Lee, S. S., Lee, S. A., Yang, J., and Kim, J. (2011). Developing stable progenies of xBrassicoraphanus, an intergeneric allopolyploid between Brassica rapa and Raphanus sativus, through induced mutation using microspore culture. Theor. Appl. Genet. 122, 885–891. doi: 10.1007/s00122-010-1494-3
Li, M. Y., Song, X., Wang, F., and Xiong, A. S. (2016a). Suitable Reference Genes for Accurate Gene Expression Analysis in Parsley (Petroselinum crispum) for abiotic stresses and hormone stimuli. Front. Plant Sci. 7:1481. doi: 10.3389/fpls.2016.01481
Li, M. Y., Wang, F., Jiang, Q., Wang, G. L., Tian, C., and Xiong, A. S. (2016b). 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
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Mamo, S., Gal, A. B., Bodo, S., and Dinnyes, A. (2007). Quantitative evaluation and selection of reference genes in mouse oocytes and embryos cultured in vivo and in vitro. BMC Dev. Biol. 7, 14. doi: 10.1186/1471-213X-7-14
Medrano, G., Guan, P., Barlow-Anacker, A. J., and Gosain, A. (2017). Comprehensive selection of reference genes for quantitative RT-PCR analysis of murine extramedullary hematopoiesis during development. PLoS ONE 12:e0181881. doi: 10.1371/journal.pone.0181881
Mitsui, Y., Shimomura, M., Komatsu, K., Namiki, N., Shibatahatta, M., Imai, M., et al. (2015). The radish genome and comprehensive gene expression profile of tuberous root formation and development. Sci. Rep. 5:10835. doi: 10.1038/srep10835
Nazari, F., Parham, A., and Maleki, A. F. (2015). GAPDH, β-actin and β2-microglobulin, as three common reference genes, are not reliable for gene expression studies in equine adipose- and marrow-derived mesenchymal stem cells. J. Anim. Sci. 57, 1–8. doi: 10.1186/s40781-015-0050-8
Niu, L., Tao, Y. B., Chen, M. S., Fu, Q., Li, C., Dong, Y., et al. (2015). Selection of reliable reference genes for gene expression studies of a promising oilseed crop, Plukenetia volubilis, by real-time quantitative PCR. Int. J. Mol. Sci. 16, 12513–12530. doi: 10.3390/ijms160612513
Niu, X., Qi, J., Chen, M., Zhang, G., Tao, A., Fang, P., et al. (2015). Reference genes selection for transcript normalization in kenaf (Hibiscus cannabinus L.) under salinity and drought stress. PeerJ. 3:e1347. doi: 10.7717/peerj.1347
Paraskevopoulos, K., Kriegenburg, F., Tatham, M. H., Rosner, H. I., Medina, B., Larsen, I. B., et al. (2014). Dss1 is a 26S proteasome ubiquitin receptor. Mol. Cell 56, 453–461. doi: 10.1016/j.molcel.2014.09.008
Peterka, H., Budahn, H., Schrader, O., Ahne, R., and Schutze, W. (2004). Transfer of resistance against the beet cyst nematode from radish (Raphanus sativus) to rape (Brassica napus) by monosomic chromosome addition. Theor. Appl. Genet. 109, 30–41. doi: 10.1007/s00122-004-1611-2
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
Pollier, J., Vanden Bossche, R., Rischer, H., and Goossens, A. (2014). Selection and validation of reference genes for transcript normalization in gene expression studies in Catharanthus roseus. Plant Physiol. Biochem. 83, 20–25. doi: 10.1016/j.plaphy.2014.07.004
Qi, J., Yu, S., Zhang, F., Shen, X., Zhao, X., Yu, Y., et al. (2010). Reference gene selection for real-time quantitative polymerase chain reaction of mrna transcript levels in chinese cabbage (Brassica rapa L. ssp. pekinensis). Plant Mol. Biol. Rep. 28, 597–604. doi: 10.1007/s11105-010-0185-1
Radonić, A., Thulke, S., Mackay, I. M., Landt, O., Siegert, W., and Nitsche, A. (2004). Guideline to reference gene selection for quantitative real-time PCR. Biochem. Biophys. Res. Commun. 313, 856–862. doi: 10.1016/j.bbrc.2003.11.177
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
Reid, K. E., Olsson, N., Schlosser, J., Peng, F., and Lund, S. T. (2006). An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biol. 6:27. doi: 10.1186/1471-2229-6-27
Remans, T., Smeets, K. K., Mathijsen, D., Vangronsveld, J., and Cuypers, A. (2008). Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations. Planta 227, 1343–1349. doi: 10.1007/s00425-008-0706-4
Roeder, A. H. K., Ferrándiz, C., and Yanofsky, M. F. (2003). The role of the REPLUMLESS homeodomain protein in patterning the arabidopsis fruit. Curr. Biol. 13, 1630–1635. doi: 10.1016/j.cub.2003.08.027
Schmid, M., Davison, T. S., Henz, S. R., Pape, U. J., Demar, M., Vingron, M., et al. (2005). A gene expression map of Arabidopsis thaliana development. Nat. Genet. 37, 501–506. doi: 10.1038/ng1543
Schmidt, G. W., and Delaney, S. K. (2010). Stable internal reference genes for normalization of real-time RT-PCR in tobacco (Nicotiana tabacum) during development and abiotic stress. Mol. Genet. Genomics 283, 233–241. doi: 10.1007/s00438-010-0511-1
Shivhare, R., and Lata, C. (2016). Selection of suitable reference genes for assessing gene expression in pearl millet under different abiotic stresses and their combinations. Sci. Rep. 6:23036. doi: 10.1038/srep23036
Silver, N., Best, S., Jiang, J., and Thein, S. L. (2006). Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol. Biol. 7:33. doi: 10.1186/1471-2199-7-33
Sirover, M. A. (1999). New insights into an old protein: the functional diversity of mammalian glyceraldehyde-3-phosphate dehydrogenase. Biochim. Biophys. Acta 1432, 159–184. doi: 10.1016/S0167-4838(99)00119-3
Suzuki, T., Higgins, P. J., and Crawford, D. R. (2000). Control selection for RNA quantitation. BioTechniques 29, 332–337.
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
Tonosaki, K., Michiba, K., Bang, S. W., Kitashiba, H., Kaneko, Y., and Nishio, T. (2013). Genetic analysis of hybrid seed formation ability of Brassica rapa in intergeneric crossings with Raphanus sativus. Theor. Appl. Genet. 126, 837–846. doi: 10.1007/s00122-012-2021-5
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
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
Wang, C., Cui, H. M., Huang, T. H., Liu, T. K., Hou, X. L., and Li, Y. (2016). Identification and validation of reference genes for RT-qPCR analysis in non-heading chinese cabbage flowers. Front. Plant Sci. 7:811. doi: 10.3389/fpls.2016.00811
Wang, H., Wang, J., Jiang, J., Chen, S., Guan, Z., Liao, Y., et al. (2014). Reference genes for normalizing transcription in diploid and tetraploid Arabidopsis. Sci. Rep. 4:6781. doi: 10.1038/srep06781
Wei, X., Zhang, X., Shen, D., Wang, H., Wu, Q., Lu, P., et al. (2013). Correction: transcriptome analysis of Barbarea vulgaris infested with diamondback moth (Plutella xylostella) Larvae. PLoS ONE 8:e64481. doi: 10.1371/annotation/83c5e6a5-2633-46d8-b138-2ac1f0f2706a
Xiao, D., Zhang, N.-W., Zhao, J.-J., Bonnema, G., and Hou, X.-L. (2012). Validation of reference genes for real-time quantitative PCR normalisation in non-heading Chinese cabbage. Funct. Plant Biol. 39, 342. doi: 10.1071/FP11246
Xie, F., Xiao, P., 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
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
Xu, Y., Zhu, X., Chen, Y., Gong, Y., and Liu, L. (2013). Expression profiling of genes involved in ascorbate biosynthesis and recycling during fleshy root development in radish. Plant Physiol. Biochem. 70, 269–277. doi: 10.1016/j.plaphy.2013.05.041
Xu, Y., Zhu, X., Gong, Y., Xu, L., Wang, Y., and Liu, L. (2012). Evaluation of reference genes for gene expression studies in radish (Raphanus sativus L.) using quantitative real-time PCR. Biochem. Biophys. Res. Commun. 424, 398–403. doi: 10.1016/j.bbrc.2012.06.119
Yang, H., Liu, J., Huang, S., Guo, T., Deng, L., and Hua, W. (2014). Selection and evaluation of novel reference genes for quantitative reverse transcription PCR (qRT-PCR) based on genome and transcriptome data in Brassica napus L. Gene 538, 113–122. doi: 10.1016/j.gene.2013.12.057
Yu, R., Wang, Y., Xu, L., Zhu, X., Zhang, W., Wang, R., et al. (2015). Transcriptome profiling of root microRNAs reveals novel insights into taproot thickening in radish (Raphanus sativus L.). BMC Plant Biol. 15:4723. doi: 10.1186/s12870-015-0427-3
Zhang, F. L., Wang, M., Liu, X. C., Zhao, X. Y., and Yang, J. P. (2007). Quantitative trait loci analysis for resistance against Turnip mosaic virus based on a doubled-haploid population in Chinese cabbage. Plant Breed. 127, 82–86. doi: 10.1111/j.1439-0523.2007.01431.x
Keywords: radish, Chinese cabbage, distant hybrid, organs, stress, pistil development, RT-qPCR, reference gene
Citation: Duan M, Wang J, Zhang X, Yang H, Wang H, Qiu Y, Song J, Guo Y and Li X (2017) Identification of Optimal Reference Genes for Expression Analysis in Radish (Raphanus sativus L.) and Its Relatives Based on Expression Stability. Front. Plant Sci. 8:1605. doi: 10.3389/fpls.2017.01605
Received: 11 March 2017; Accepted: 31 August 2017;
Published: 15 September 2017.
Edited by:
Roger Deal, Emory University, United StatesReviewed by:
Qiusheng Kong, Huazhong Agricultural University, ChinaHuipeng Pan, South China Agricultural University, China
Copyright © 2017 Duan, Wang, Zhang, Yang, Wang, Qiu, Song, Guo and Li. 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.
*Correspondence: Xixiang Li, lixixiang@caas.cn