Evaluation of New Reference Genes in Papaya for Accurate Transcript Normalization under Different Experimental Conditions

<div><p>Real-time reverse transcription PCR (RT-qPCR) is a preferred method for rapid and accurate quantification of gene expression studies. Appropriate application of RT-qPCR requires accurate normalization though the use of reference genes. As no single reference gene is universally suitable for all experiments, thus reference gene(s) validation under different experimental conditions is crucial for RT-qPCR analysis. To date, only a few studies on reference genes have been done in other plants but none in papaya. In the present work, we selected 21 candidate reference genes, and evaluated their expression stability in 246 papaya fruit samples using three algorithms, geNorm, NormFinder and RefFinder. The samples consisted of 13 sets collected under different experimental conditions, including various tissues, different storage temperatures, different cultivars, developmental stages, postharvest ripening, modified atmosphere packaging, 1-methylcyclopropene (1-MCP) treatment, hot water treatment, biotic stress and hormone treatment. Our results demonstrated that expression stability varied greatly between reference genes and that different suitable reference gene(s) or combination of reference genes for normalization should be validated according to the experimental conditions. In general, the internal reference genes <em>EIF</em> (Eukaryotic initiation factor 4A), <em>TBP1</em> (TATA binding protein 1) and <em>TBP2</em> (TATA binding protein 2) genes had a good performance under most experimental conditions, whereas the most widely present used reference genes, <em>ACTIN</em> (Actin 2), <em>18S rRNA</em> (18S ribosomal RNA) and <em>GAPDH</em> (Glyceraldehyde-3-phosphate dehydrogenase) were not suitable in many experimental conditions. In addition, two commonly used programs, geNorm and Normfinder, were proved sufficient for the validation. This work provides the first systematic analysis for the selection of superior reference genes for accurate transcript normalization in papaya under different experimental conditions.</p> </div

Summary of the experimental conditions and samples used in present study.

<p>#indicated that the sample dates including two types: two treatments were 6 and one was 7.</p

Average expression stability values (M) of the candidate reference genes.

<p>Average expression stability values (M) of the reference genes were measured during stepwise exclusion of the least stable reference genes. A lower M value indicated more stable expression, as analyzed by the geNorm software in papaya sample sets under different experimental conditions, including different storage temperatures (a), hot water treatment (b), modified atmosphere packaging (c), hormone treatment (d), 1-MCP fumigation treatment (e), different developmental stages (f), different tissues (g), biotic stress (h), postharvest ripening: cultivar of ‘<i>Shuiyou 2′</i> (i), cultivar of ‘<i>Hongri 1′</i> (j), cultivar of ‘<i>Hongri 3′</i> (k). Different cultivars samples (l) and all papaya samples (m) were also given.</p

Determination of the optimal number of reference genes.

<p>Pair-wise variation (V) calculated by geNorm to determine the minimum number of reference genes for accurate normalization in different experiment conditions. Arrow indicates the optimal number of genes for normalization in each sample sets.</p

RT-qPCR CT values for the candidate reference genes.

<p>Expression data displayed as CT values for each reference gene in all papaya samples. A line across the box is depicted as the median. The box indicates the 25th and 75th percentiles. Whiskers represent the maximum and minimum values.</p

Consensus of stability ranking of the reference gene estimated by geNorm and NormFinder.

<p>Consensus of stability ranking of the reference gene estimated by geNorm and NormFinder.</p

Specificity of primer pairs for RT-qPCR amplification.

<p>Equal amounts of cDNAs from all tested samples were mixed as the template. 2.5% non-denaturing agarose gel electrophoresis showed amplification of a specific product of the expected size for each reference gene. M represented DNA size marker.</p

Relative quantification of <i>CpaEXY1</i> expression using validated reference genes for normalization under different experimental conditions.

<p>(a) The validated reference gene(s) used as normalization factors were one (<i>UBCE</i>) or two (<i>UBCE</i>+<i>CYP</i>) most stable reference genes, and one most unstable gene (<i>GAPDH</i>) in postharvest ripening of ‘<i>Shuiyou 2′</i> sample sets. (b) The validated reference gene(s) used as normalization factors were one (<i>UBQ</i>) or two (<i>UBQ</i>+<i>CYP</i>) most stable reference genes, and the most unstable one (<i>CHY</i>) in postharvest ripening of ‘<i>Hongri 3′</i> sample sets. Reference genes validated by geNorm or NormFinder. Each value represented the means of three replicates, and vertical bars indicate the standard deviations (SD). In <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0044405#pone-0044405-g005" target="_blank">Figure 5b</a>, <i>UBQ</i> +<i>CYP</i> and <i>UBQ</i> normalized curves belonged to left y axis, and <i>CHY</i> normalized curve belonged to right y axis.</p