Selection of Novel Reference Genes by RNA-Seq and Their Evaluation for Normalising Real-Time qPCR Expression Data of Anthocyanin-Related Genes in Lettuce and Wild Relatives

Lettuce is a popular vegetable source of bioactive compounds, like anthocyanins, powerful antioxidants present in red and semi-red varieties. Selection of reliable reference genes (RGs) for the normalization of real-time quantitative PCR (qPCR) data is crucial to obtain accurate gene expression results. Among the genes with totally unrelated biological functions, six candidate RGs (ADF2, CYB5, iPGAM, SCL13, TRXL3-3, and VHA-H) with low variation in expression according to RNA-seq analyses, were selected for future expression studies of anthocyanin-related genes in three different experiments: leaf colour comparison (green vs. red) in commercial varieties; tissue comparison (leaf vs. stem) in a wild relative; and drought stress experiment in commercial and traditional varieties, and a wild relative. Expression profiles of the candidate RGs were obtained by qPCR and their stability was assessed by four different analytical tools, geNorm, NormFinder, BestKeeper, and Delta Ct method, all integrated in RefFinder. All results considered, we recommend CYB5 to be used as RG for the leaf colour experiment and TRXL3-3 for the tissue and drought stress ones, as they were the most stable genes in each case. RNA-seq is useful to preselect novel RGs although validation by qPCR is still advisable. These results provide helpful information for gene expression studies in Lactuca spp. under the described conditions.