Abstract
The popular leafy vegetable lettuce (Lactuca sativa L.) is susceptible to cold stress during the growing season, which slows growth rate, causes leaf yellowing and necrosis, and reduced yield and quality. In this study, transcriptomic and metabolomic analyses of two cold-resistant lettuce cultivars (GWAS-W42 and F11) and two cold-sensitive lettuce cultivars (S13K079 and S15K058) were performed to identify the mechanisms involved in the cold response of lettuce. Overall, transcriptome analysis identified 605 differentially expressed genes (DEGs), including significant enrichment of genes involved in the flavonoid and flavonol (CHS, CHI, F3H, FLS, CYP75B1, HCT, etc.) biosynthetic pathways related to oxidation-reduction and catalytic activity. Untargeted metabolomic analysis identified fifteen flavonoid metabolites and 28 other metabolites potentially involved in the response to cold stress; genistein, quercitrin, quercetin derivatives, kaempferol derivatives, luteolin derivatives, apigenin and their derivatives accumulate at higher levels in cold-resistant cultivars. Moreover, MYBs, bHLHs, WRKYs and Dofs also play positive role in the low temperature response, which affected the expression of structural genes contributing to the variation of metabolites between the resistant and sensitive. These results provide valuable evidence that the metabolites and genes involved in the flavonoid biosynthetic pathway play important roles in the response of lettuce to cold stress.