Selection and Validation of Reference Genes for RT-qPCR Analysis in Spinacia oleracea under Abiotic Stress

Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is an accurate and convenient method for mRNA quantification. Selection of optimal reference gene(s) is an important step in RT-qPCR experiments. However, the stability of housekeeping genes in spinach (Spinacia oleracea) under various abiotic stresses is unclear. Evaluating the stability of candidate genes and determining the optimal gene(s) for normalization of gene expression in spinach are necessary to investigate the gene expression patterns during development and stress response. In this study, ten housekeeping genes, 18S ribosomal RNA (18S rRNA), actin, ADP ribosylation factor (ARF), cytochrome c oxidase subunit 5C (COX), cyclophilin (CYP), elongation factor 1-alpha (EF1 alpha), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), histone H3 (H3), 50S ribosomal protein L2 (RPL2), and tubulin alpha chain (TUB alpha) from spinach, were selected as candidates in roots, stems, leaves, flowers, and seedlings in response to high temperature, CdCl2, NaCl, NaHCO3, and Na2CO3 stresses. The expression of these genes was quantified by RT-qPCR and evaluated by NormFinder, BestKeeper, and geNorm. 18S rRNA, actin, ARF, COX, CYP, EF1 alpha, GAPDH, H3, and RPL2 were detected as optimal reference genes for gene expression analysis of different organs and stress responses. The results were further confirmed by the expression pattern normalized with different reference genes of two heat-responsive genes. Here, we optimized the detection method of the gene expression pattern in spinach. Our results provide the optimal candidate reference genes which were crucial for RT-qPCR analysis.