Screening and Verification of Reference Genes for Analysis of Gene Expression in Garlic (Allium sativum L.) under Cold and Drought Stress

The principal objective of this study was to screen and verify reference genes appropriate for gene expression evaluation during plant growth and development under distinct growth conditions. Nine candidate reference genes were screened based on garlic transcriptome sequence data. RT-qPCR was used to detect the expression levels of the aforementioned reference genes in specific tissues under drought and cold stress. Then, geNorm, NormFinder, BestKeeper, and ReFinder were used to consider the consistency of the expression levels of candidate reference genes. Finally, the stress-responsive gene expression of ascorbate peroxidase (APX) was quantitatively evaluated to confirm the chosen reference genes. Our results indicated that there were variations in the abundance and stability of nine reference gene transcripts underneath cold and drought stress, among which ACT and UBC-E2 had the highest transcript abundance, and 18S rRNA and HIS3 had the lowest transcript abundance. UBC and UBC-E2 were the most stably expressed genes throughout all samples; UBC and UBC-E2 were the most stably expressed genes during cold stress, and ACT and UBC were the most stably expressed genes under drought stress. The most stably expressed genes in roots, pseudostems, leaves, and cloves were EF1, ACT, HIS3, UBC, and UBC-E2, respectively, while GAPDH was the most unstable gene during drought and cold stress conditions and in exclusive tissues. Taking the steady reference genes UBC-E2, UBC, and ACT as references during drought and cold stress, the reliability of the expression levels was further demonstrated by detecting the expression of AsAPX. Our work thereby offers a theoretical reference for the evaluation of gene expression in garlic in various tissues and under stress conditions.