Original Research ARTICLE
Selection of Reference Genes for Quantitative Real-Time RT-PCR Studies in Tomato Fruit of the Genotype MT-Rg1
- 1Unidad de Genómica Avanzada – Laboratorio Nacional de Genómica para la Biodiversidad, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Irapuato, Mexico
- 2Laboratório de Genética Molecular Vegetal, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
Quantitative real-time RT-PCR (qRT-PCR) has become one of the most widely used methods for accurate quantification of gene expression. Since there are no universal reference genes for normalization, the optimal strategy to normalize raw qRT-PCR data is to perform an initial comparison of a set of independent reference genes to assess the most stable ones in each biological model. Normalization of a qRT-PCR experiment helps to ensure that the results are both statistically significant and biologically meaningful. Tomato is the model of choice to study fleshy fruit development. The miniature tomato (Solanum lycopersicum L.) cultivar Micro-Tom (MT) is considered a model system for tomato genetics and functional genomics. A new genotype, containing the Rg1 allele, improves tomato in vitro regeneration. In this work, we evaluated the expression stability of four tomato reference genes, namely CAC, SAND, Expressed, and ACTIN2. We showed that the genes CAC and Exp are the best reference genes of the four we tested during fruit development in the MT-Rg1 genotype. Furthermore, we validated the reference genes by showing that the expression profiles of the transcription factors FRUITFULL1 and APETALA2c during fruit development are comparable to previous reports using other tomato cultivars.
Reverse transcription followed by quantitative PCR (qRT-PCR) assay is an extremely sensitive technique that provides accurate and reproducible quantification of nucleic acids based on the exponential incorporation of fluorescent molecules into genetic material (Gachon et al., 2004; Nolan et al., 2006). Nowadays, qRT-PCR analysis has become the method of choice for gene expression studies and validating transcriptomic data. One of the most crucial points in RT-qPCR data analysis is the choice of a proper normalization method. The purpose of normalization is to correct for variability associated with the experimental procedure, such as the amount of starting material, RNA extraction and enzymatic efficiencies, and differences in overall transcriptional activity between tissues or cells (Expósito-Rodríguez et al., 2008; Gutierrez et al., 2008a,b; Udvardi et al., 2008). To date, the parallel quantification of endogenous reference genes is accepted as the most reliable method for sample normalization. The normalization of relative quantities with reference genes relies on the assumption that the reference genes are stably expressed across all tested samples, and that they are not being significantly altered across treatments or conditions (Vandesompele et al., 2002). However, several reports have shown that the transcript levels of commonly used reference genes, known as housekeeping genes, can vary considerably under different experimental conditions (e.g., Thellin et al., 1999; Suzuki et al., 2000; Czechowski et al., 2005; Jain et al., 2006). Moreover, a reference gene with stable expression in one organism may not be suitable for normalization of gene expression in another organism (e.g., Jain et al., 2006; Jian et al., 2008; Manoli et al., 2012; Lambret-Frotté et al., 2015; Kanakachari et al., 2016).
Since there are no universal reference genes, the optimal strategy to normalize raw qRT-PCR data is to perform an initial comparison of a set of independent reference genes to assess the most stable ones in each particular experimental background or biological model. The use of multiple reference genes does not only produce more reliable data, but permits an evaluation of the stability of the reference genes themselves. In this way, normalization of a qPCR assay helps to ensure that the results are both statistically significant and biologically meaningful (Hellemans et al., 2007; Expósito-Rodríguez et al., 2008; Gutierrez et al., 2008a,b; De Spiegelaere et al., 2015).
Fruits are an important evolutionary acquisition of angiosperms, which convey protection to seeds and ensure their optimal dispersal in the environment. Tomato is the model of choice to study fleshy fruit development (Seymour et al., 2013; Pesaresi et al., 2014). The miniature tomato or micro tomato (Solanum lycopersicum L.) cultivar Micro-Tom (MT) is considered a model system for tomato genetics and functional genomics (Scott and Harbaugh, 1989; Meissner et al., 1997; Emmanuel and Levy, 2002; Dan et al., 2006; Sun et al., 2006; Carvalho et al., 2011). This variety displays convenient traits for plant research, such as short life cycle, miniature size, high density growth, and the suitability for transgenic plant production and large scale mutagenesis (Meissner et al., 1997; Emmanuel and Levy, 2002; Watanabe et al., 2007; Carvalho et al., 2011).
In order to enhance tomato genetic transformation, Pino et al. (2010) developed a MT near-isogenic genotype harboring the Rg1 allele from S. peruvianum. This new genotype, MT-Rg1, allows reduced exposure to exogenous hormone applications during transformation protocols, improving in vitro regeneration (Pino et al., 2010).
Although tomato has become an important model for genetic and molecular studies in fleshy fruits, there are few examples in literature where best reference genes have been identified and validated in tomato. The best example is by Expósito-Rodríguez et al. (2008), with a focus on different tomato tissues. A few other studies report on appropriate reference gene selection for expression analysis in tomato seeds (Dekkers et al., 2012) and pathogen infected tomato plants (Lacerda et al., 2015; Muller et al., 2015). Therefore, there is still a need for more studies on appropriate reference gene selection for expression studies during different conditions and developmental processes, for instance tomato fruit development.
In this work, we evaluated the expression stability of four commonly used tomato reference genes, namely CAC, SAND, Expressed, and ACTIN2, and evaluated the expression of the transcription factors FRUITFULL1 (FUL1) and APETALA2c (AP2c) across eight MT-Rg1 fruit developmental stages. Our results allowed us to select suitable reference genes for qRT-PCR studies during fruit development, and show that the expression profiles of FUL1 and AP2c are similar to those previously reported for other tomato fruit cultivars.
Biological Samples and Candidate Reference Genes
Tomato Micro-Tom Rg1 (MT-Rg1) fruits were sampled at eight different developmental stages (Figure 1), according to Gillaspy et al. (1993). The collected fruit stages are as follows: (1) 1.5 cm immature green fruit, (2) 2 cm immature green fruit, (3) mature green fruit, (4) breaker, (5) turning, (6) orange, (7) red firm, and (8) red ripe fruit. Fruits (pericarp tissue only) were frozen and stored at –80°C until further use. Total RNA was isolated using Trizol, followed by several precipitation steps, a DNase I treatment, and subsequently, cDNA was prepared, as described in the Materials and Methods section.
FIGURE 1. MT-Rg1 fruit developmental stages used for gene expression analysis in this study. Scale bar: 1 cm.
Based on previous studies conducted on S. lycopersicum cv. ciliegia (Expósito-Rodríguez et al., 2008), a total of four candidate reference genes were selected for qRT-PCR normalization. These are CAC, SAND, and Expressed (Expósito-Rodríguez et al., 2008), together with the commonly used house-keeping gene ACTIN2 (Bemer et al., 2012) (Table 1). Primers used are listed in Table 2. The primers for the first three genes span intronic regions (Table 2), i.e., a larger amplicon will be seen when DNA contamination is present. Primer melting curves for all genes showed a unique peak corresponding to the expected amplicon (Supplementary Figure 1). The correct size of the amplicons was verified by gel electrophoreses.
TABLE 1. Candidate reference genes and validation genes used for qPCR expression study in Solanum lycopersicum cv. MT-Rg1 during fruit development stages.
TABLE 2. Details of primers of candidate reference genes, validation genes and parameters derivated from qPCR analysis.
Finally, as described at the end of the Results section, relative expression of the transcription factors FUL1 (Bemer et al., 2012) and AP2c (Karlova et al., 2011) was analyzed and compared to previous reports (Karlova et al., 2011; Bemer et al., 2012) in order to validate the normalization procedure using the most stable reference genes identified during fruit development in MT-Rg1.
Expression Stability of Candidate Reference Genes
Data processing is illustrated in Figure 2. To evaluate the individual reaction kinetics, without the need for a standard curve, we used Real-Time PCR Miner software (Zhao and Fernald, 2005). Raw data of the time (cycle) sequence fluorescence values were imported from the real time PCR machine into Miner to calculate primer efficiency and its associated standard deviation, and the cycle threshold (Ct) value. The resulting data was used as input for qBase.
FIGURE 2. Data processing workflow for quantitative PCR data. Each program generates an input file needed for next step in the workflow.
The qBase software (Hellemans et al., 2007) processes data modules as independent experiments. Thus Ct values for each individual reaction, and each tested reference gene, were imported as independent experiments. The first step in the qBase workflow combines raw data from all individual run files of the same experiment into a single data table, where data points with identical sample and gene names are automatically identified as technical replicates. Then, we executed a data quality control, where replicated reactions that differ between them in more than 0.5 Ct were excluded from further analysis as a potential outlier, totaling 19 excluded cases (6.59% of the complete data set). Next, average efficiency of genes and its associated standard deviation were replaced with the values obtained from the Miner software for each gene. Finally, we performed the quantification of relative quantities, which lie in the conversion of quantification cycle values into relative quantities based on the gene specific amplification efficiency.
Then a normalization factor based on the expression levels of the tested reference genes was calculated using both geNorm (Vandesompele et al., 2002) and NormFinder (Andersen et al., 2004), which are Excel based software packages.
geNorm calculates a gene stability measure M as the average pairwise variation of a particular gene against other reference genes (Vandesompele et al., 2002) (Figure 3A). Genes are ranked according to increasing expression stability, with CAC and Exp displaying the most stable expression profiles and ACT and SAND less stable profiles. The CAC and Exp genes were also reported as one of the most stable reference genes in S. lycopersicum cv. ciliegia (Expósito-Rodríguez et al., 2008). Next, normalization factors were calculated for the most stable reference genes. Additionally, the software calculates the pairwise variation (V) between two sequential normalization factors (NFn and NFn + 1) in order to determine the minimum number of reference genes required for normalization (Figure 3B) (Vandesompele et al., 2002). In other words, pairwise variation estimates how the NF changes if you add another reference gene to the NF calculation. geNorm pairwise comparison shows that adding a third gene has no significant contribution to the newly calculated normalization factor (V2/3), meaning that using two reference genes is sufficient to normalize our qRT-PCR data. Its value is lower than the cutoff of 0.15 proposed by Vandesompele et al. (2002), which suggests the inclusion of additional reference genes are not necessary. From this analysis we concluded the use of the two most stable reference genes CAC and Exp is sufficient for accurate normalization (V2/3 = 0.103) in MT-Rg1 fruits.
FIGURE 3. Expression stability values (M) and pairwise variation analysis of candidate reference genes calculated by geNorm. (A) geNorm average expression stability (M) of the four candidate reference genes. Lower values indicate increased gene stability across samples. In total 16 samples (including both biological replicates) were included in this analysis. (B) The geNorm pairwise variation (V) was analyzed between the normalization factors (NF). Using two reference genes instead of three genes (V2/3) resulted in a value below the cutoff of 0.15, i.e., that the use of two reference gene is sufficient for normalization.
NormFinder (Andersen et al., 2004) is an algorithm to identify the optimal normalization gene among a set of candidates. The validity of this approach is related to the number of samples and candidates analyzed because it considers intra- and intergroup variation to estimate a stability value of the tested reference genes. The intragroup refers to the biological replicates and the intergroup refers to the different fruit developmental stages. We used the input file containing the qBase-calculated relative quantity values for each gene, sample and biological replica and ran the analysis. The results revealed CAC as the most stable reference gene with the best stability value (i.e., the lowest value) of 0.049 (Figure 4). The best combination of genes was CAC and Exp. These results match those proposed by the geNorm software.
FIGURE 4. NormFinder expression stability values (M) and ranking of the candidate reference genes. Lower values indicate more stable expression.
Consequently, we decided to continue the experiment using the best reference genes CAC and Exp to normalize the gene expression data of FUL1 and AP2c and compare it with normalized expression data with the less stable reference gene SAND in all tested samples.
AP2c and FUL1 Expression Profiles
To assess the validity of the MT-Rg1 genotype as a model tomato species, we calculated the relative expression of AP2c and FUL1 at different developmental stages, and compared them with previous reports (Karlova et al., 2011; Bemer et al., 2012). Accurate measurement of gene expression requires a normalization by multiple, carefully selected reference genes. We therefore decided to normalize gene expression using both the most stable reference gene pair identified by GeNorm and NormFinder, CAC and Exp, but also using the less stable reference gene, SAND.
AP2 is a transcription factor which plays important roles in development, ethylene response and pathogen resistance (Karlova et al., 2011). Five distinct tomato cDNAs encoding AP2 putative tomato orthologs were identified by Karlova et al. (2011) in the Moneymaker tomato cultivar. From this report, we chose the AP2c sequence, which shows the highest expression in early stage fruits, but no significant expression beyond the mature green stage. Our results in MT-Rg1 fruits (Figure 5A) were concordant with the previously reported expression profile in the tomato cultivar Moneymaker (Karlova et al., 2011). In the first three fruit stages, we observed a high expression and then expression levels dropped from the breaker stage on.
FIGURE 5. Expression profiles of (A) AP2c and (B) FUL1 in MT-Rg1 fruits. Dataset was normalized using the two best reference genes (CAC and Exp) and compared with normalized data using the poorly ranked reference gene (SAND). Error bars represent the standard error of the mean.
Tomato has two orthologs of FUL, a transcription factor involved in fruit ripening (Bemer et al., 2012). We also calculated the relative expression of FUL1 in order to validate our MT-Rg1 qRT-PCR results. FUL1 expression was found to be very low during early stages of fruit development, but rapidly increased from the breaker stage, reaching its maximum in the orange ripe stage (Figure 5B). A comparable trend of expression was reported in tomato MT fruits (Bemer et al., 2012). Normalizing the qRT-PCR results for AP2c and FUL1 with the less stable reference gene SAND did result in differences in expression levels in the fruit stages with high AP2c and FUL1 expression levels, although the trend of the expression profile did not change substantively.
Several reports have shown the importance of selecting proper reference genes for data normalization, and how the identity of these genes will vary depending on the model of study (Jain et al., 2006; Expósito-Rodríguez et al., 2008; Jian et al., 2008; Manoli et al., 2012; Lambret-Frotté et al., 2015; Kanakachari et al., 2016). Although reference genes have been identified for the S. lycopersicum cv. ciliegia, there are no reports about the most reliable reference genes for the MT-Rg1 tomato genotype.
This work constitutes an effort to validate appropriate reference genes for the quantification of transcript levels by quantitative RT-PCR in the MT-Rg1 tomato genotype. We have tested the expression stabilities of four reference candidate genes (CAC, SAND, Expressed, and ACT2) in a set of eight developmental tomato MT-Rg1 fruit stages.
Correct sample normalization is an absolute prerequisite for reliable and accurate measurement of gene expression. The optimal strategy to normalize raw qPCR data is to perform an initial comparison of a set of independently regulated reference genes to assess the most stable ones in each specific experiment or biological setting (e.g., Czechowski et al., 2005; Expósito-Rodríguez et al., 2008; Gutierrez et al., 2008a; De Spiegelaere et al., 2015).
Because of the increasing attention on a proper normalization of qRT-PCR data, there are increasing number of methods and software packages that have been developed for the validation of the most stable reference genes. Specialized software like geNorm (Vandesompele et al., 2002) and NormFinder (Andersen et al., 2004) are excellent tools to determine best reference genes when a new system or tissue is used. Reference genes should be established for each tested tissue, allowing a better interpretation and biological significance. Unstable reference genes, if used for normalization, can radically change the expression pattern of a given gene under study causing errors in results and, thereby, the interpretation or understanding of gene function (Gutierrez et al., 2008a,b). This emphasizes the importance of preliminary evaluation studies, aimed to identify the most stable reference genes in different organism and also between different tissues of the same species.
The suitability of reference genes identified in this study (CAC and Exp) was validated through an assessment of the expression profiles of two transcription factors. The expression of AP2c and FUL1 during fruit development in the MT-Rg1 genotype showed an activity similar to other tomato cultivars (Karlova et al., 2011; Bemer et al., 2012).
In summary, here we showed that the genes CAC and Exp are appropriate reference genes during fruit development in the MT-Rg1 genotype. Furthermore, we validated the normalization method for the tomato fruit developmental genes AP2c and FUL1.
Materials and Methods
Plant Growth Conditions
Tomato (S. lycopersicum) MT-Rg1 genotype (Pino et al., 2010) seeds were sown in flat trays containing a 3:1:1:1 mixture of leaf soil, soil, sand, and perlite, respectively. Plants were grown in a greenhouse under local conditions (20°43′1″ N; 101°19′56″) in the summer (around 16 h light/8 h dark) at an average mean temperature of 30°C.
Tomato fruits were sampled at eight different developmental stages, based on the fruit stage division proposed by Gillaspy et al. (1993). In total, four fruits of each developmental stage were collected: 1.5 cm immature green fruits, 2 cm immature green fruit, mature green fruit, breaker, turning, orange, red firm, and red ripe (Figure 1). Fruits devoid of seeds and placental tissue (only pericarp tissue) were frozen in liquid nitrogen and stored at -80°C until RNA extraction.
RNA Extraction and cDNA Synthesis
Total RNA extraction was performed using the TRIzol® (Ambion) protocol, followed by precipitation using 0.8 M sodium citrate, 1.2 M NaCl and isoproponol, a second precipitation using 8 M LiCl, and a third precipitation with 3 M NaAc pH 5.2 and 100% ethanol. Subsequently, a DNase I (Life Technologies) treatment was performed according to manufacturer’s specifications, and RNA was recovered using phenol/chloroform extraction followed by a precipitation with 3 M NaAc pH 5.2 and 100% ethanol. RNA was resuspended in 10 μl of DEPC water and quantified by measuring its absorbance at 260 nm. RNA integrity was evaluated by the 260/280 and 260/230 ratios, and confirmed by agarose gel electrophoresis.
cDNA was synthesized from 1.3 μg of total RNA using the SuperScript® III System (Life Technologies). In summary, total RNA was mixed with 1 μl of oligo dT (50 μM), 1 μl of dNTPs (10 μM) and MQ water, giving a total volume of 14 μl, and incubated for 5 min at 65°C and then chilled on ice. Subsequently, 4 μl of First Strand Buffer (5x), 1 μl of DTT (0.1 M) and 1 μl of SuperScript reverse transcriptase III (200 units/ μl) were added, each reaction was incubated for 2 h at 50°C and, finally, inactivated for 5 min at 70°C. cDNA was diluted 1:50 for use in quantitative real-time PCR experiments.
Quantitative Real-Time PCR
Quantitative real-time PCR amplification reactions were performed and run in technical triplicate on 96-wells plates on an ABI 7500 Fast Real Time PCR System using SYBR® Green I (Life Technologies). Each PCR reaction mix consisted of 2 μl of SYBR Green (1:10000), 0.4 μl of forward and reverse oligos (10 mM), 2 μl of PCR buffer (10x), 0.05 μl of dNTPs (10 mM), 1.2 μl of MgCl2 (50 mM) and 0.05 μl of Platinum® Taq DNA Polymerase (Life Technologies; 2 U/rxn) in a total volume of 10 μl. Finally, 10 μl of 1:50 diluted template cDNA was added, resulting in a total volume of 20 μl per PCR reaction. PCR cycling was performed as follows: 5 min at 94°C followed by 40 rounds of 15 s at 94°C, 10 s at 60°C, 15 s at 72°C, and finally 1 round of 35 s at 60°C. Melting curve cycling consisted of: 15 s at 95°C, 1 min at 60°C, 30 s at 95°C, and 15 s at 60°C.
The qRT-PCR experiment was performed using two biological replicates. For the first one, total RNA extracted from one fruit of each stage was used. cDNA for the second biological replicate was prepared from total RNA extracted from a pool of three fruits of each stage. For each analysis, three technical replicates were performed. The Supplementary Figure 2 shows the results of the comparison between the two biological replicates and the standard deviation of the three technical replicates for each fruit stage. In general, the detected Ct values did not differ much, but slightly more stability was observed in those where the cDNA was prepared from the RNA of the three pooled fruits. No difference in the standard deviation was observed between the two strategies (Supplementary Figure 2).
Raw fluorescence values from quantitative PCR experiments were imported into the Real-Time PCR Miner software (Zhao and Fernald, 2005). Ct values, average efficiency and standard deviation calculated by the Miner program were used as input for the qBase software (Hellemans et al., 2007). The qBase program calculates relative quantities of expression and produced an input file for the geNorm software (Vandesompele et al., 2002). The geNorm program determines the minimum number of genes required to calculate a reliable normalization factor. Finally, the qBase data file, with few modifications, was also used as an input file for the NormFinder program (Andersen et al., 2004), which considers the average expression stability between replicates and also between samples (intra- and inter-group variation, respectively) to suggest the best candidates for reference genes.
KG-A and CS performed experimental work and data analysis. KG-A, CS, MA-F, and SdF conceived the project and designed the experiments. KG-A and RACM made the figures. KG-A, RACM, and SdF drafted the manuscript. All authors read and approved the final manuscript version.
The work in the de Folter laboratory was supported by the National Council of Science and Technology of Mexico (CONACyT) grant CB-2012-177739, and SdF and MA-F acknowledge support of the Bilateral Cooperation México-Brazil (CONACyT-CNPq) project 2011-2013 (no. 147675), and CAPES and Carlos Chagas Filho Foundation for Research Support of the State of Rio de Janeiro (FAPERJ).
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.
We thank Lázaro E. P. Peres for the tomato MT-Rg1 seeds.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01386
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
Bemer, M., Karlova, R., Ballester, A. R., Tikunov, Y. M., Bovy, A. G., Wolters-Arts, M., et al. (2012). The tomato FRUITFULL homologs TDR4/FUL1 and MBP7/FUL2 regulate ethylene-independent aspects of fruit ripening. Plant Cell 24, 4437–4451. doi: 10.1105/tpc.112.103283
Carvalho, R. F., Campos, M. L., Pino, L. E., Crestana, S. L., Zsögön, A., Lima, J. E., et al. (2011). Convergence of developmental mutants into a single tomato model system: “Micro-Tom” as an effective toolkit for plant development research. Plant Methods 7:18. doi: 10.1186/1746-4811-7-18
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
Dan, Y., Yan, H., Munyikwa, T., Dong, J., Zhang, Y., and Armstrong, C. L. (2006). MicroTom–a high-throughput model transformation system for functional genomics. Plant Cell Rep. 25, 432–441. doi: 10.1007/s00299-005-0084-3
De Spiegelaere, W., Dern-Wieloch, J., Weigel, R., Schumacher, V., Schorle, H., Nettersheim, D., et al. (2015). Reference gene validation for RT-qPCR, a note on different available software packages. PLoS ONE 10:e0122515. doi: 10.1371/journal.pone.0122515
Dekkers, B. J., Willems, L., Bassel, G. W., Van Bolderen-Veldkamp, R. P., Ligterink, W., Hilhorst, H. W., et al. (2012). Identification of reference genes for RT-qPCR expression analysis in Arabidopsis and tomato seeds. Plant Cell Physiol. 53, 28–37. doi: 10.1093/pcp/pcr113
Expósito-Rodríguez, M., Borges, A. A., Borges-Pé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
Gutierrez, L., Mauriat, M., Guénin, S., Pelloux, J., Lefebvre, J.-F., Louvet, R., et al. (2008a). 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
Gutierrez, L., Mauriat, M., Pelloux, J., Bellini, C., and Van Wuytswinkel, O. (2008b). Towards a systematic validation of references in real-time rt-PCR. Plant Cell 20, 1734–1735. doi: 10.1105/tpc.108.059774
Hellemans, J., Mortier, G., De Paepe, A., Speleman, F., and Vandesompele, J. (2007). qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 8:R19. doi: 10.1186/gb-2007-8-2-r19
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
Jian, B., 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
Kanakachari, M., Solanke, A. U., Prabhakaran, N., Ahmad, I., Dhandapani, G., Jayabalan, N., et al. (2016). Evaluation of suitable reference genes for normalization of qPCR gene expression studies in Brinjal (Solanum melongena L.) during fruit developmental stages. Appl. Biochem. Biotechnol. 178, 433–450. doi: 10.1007/s12010-015-1884-8
Karlova, R., Rosin, F. M., Busscher-Lange, J., Parapunova, V., Do, P. T., Fernie, A. R., et al. (2011). Transcriptome and metabolite profiling show that APETALA2a is a major regulator of tomato fruit ripening. Plant Cell 23, 923–941. doi: 10.1105/tpc.110.081273
Lacerda, A. L., Fonseca, L. N., Blawid, R., Boiteux, L. S., Ribeiro, S. G., and Brasileiro, A. C. (2015). Reference Gene Selection for qPCR Analysis in Tomato-Bipartite Begomovirus Interaction and Validation in Additional Tomato-Virus Pathosystems. PLoS ONE 10:e0136820. doi: 10.1371/journal.pone.0136820
Lambret-Frotté, J., de Almeida, L. C. S., de Moura, S. M., Souza, F. L. F., Linhares, F. S., and Alves-Ferreira, M. (2015). Validating internal control genes for the accurate normalization of qPCR expression analysis of the novel model plant Setaria viridis. PLoS ONE 10:e0135006. doi: 10.1371/journal.pone.0135006
Manoli, A., Sturaro, A., Trevisan, S., Quaggiotti, S., and Nonis, A. (2012). Evaluation of candidate reference genes for qPCR in maize. J. Plant Physiol. 169, 807–815. doi: 10.1016/j.jplph.2012.01.019
Muller, O. A., Grau, J., Thieme, S., Prochaska, H., Adlung, N., Sorgatz, A., et al. (2015). Genome-wide identification and validation of reference genes in infected tomato leaves for quantitative RT-PCR analyses. PLoS ONE 10:e0136499. doi: 10.1371/journal.pone.0136499
Pesaresi, P., Mizzotti, C., Colombo, M., and Masiero, S. (2014). Genetic regulation and structural changes during tomato fruit development and ripening. Front. Plant Sci. 5:124. doi: 10.3389/fpls.2014.00124
Pino, L. E., Lombardi-Crestana, S., Azevedo, M. S., Scotton, D. C., Borgo, L., Quecini, V., et al. (2010). The Rg1 allele as a valuable tool for genetic transformation of the tomato “Micro-Tom” model system. Plant Methods 6:23. doi: 10.1186/1746-4811-6-23
Sun, H.-J., Uchii, S., Watanabe, S., and Ezura, H. (2006). A highly efficient transformation protocol for Micro-Tom, a model cultivar for tomato functional genomics. Plant Cell Physiol. 47, 426–431. doi: 10.1093/pcp/pci251
Thellin, O., Zorzi, W., Lakaye, B., De Borman, B., 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
Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., 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
Watanabe, S., Mizoguchi, T., Aoki, K., Kubo, Y., Mori, H., Imanishi, S., et al. (2007). Ethylmethanesulfonate (EMS) mutagenesis of Solanum lycopersicum cv. Micro-Tom for large-scale mutant screens. Plant Biotechnol. 24, 33–38. doi: 10.5511/plantbiotechnology.24.33
Keywords: qRT-PCR, expression analysis, reference genes, tomato (Solanum lycopersicum), fruit development
Citation: González-Aguilera KL, Saad CF, Chávez Montes RA, Alves-Ferreira M and de Folter S (2016) Selection of Reference Genes for Quantitative Real-Time RT-PCR Studies in Tomato Fruit of the Genotype MT-Rg1. Front. Plant Sci. 7:1386. doi: 10.3389/fpls.2016.01386
Received: 06 July 2016; Accepted: 31 August 2016;
Published: 13 September 2016.
Edited by:Robert G. Franks, North Carolina State University, USA
Copyright © 2016 González-Aguilera, Saad, Chávez Montes, Alves-Ferreira and de Folter. 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: Stefan de Folter, email@example.com