Integrative physiological, transcriptome, and metabolome analyses reveal the associated genes and metabolites involved in cold stress response in common vetch (Vicia sativa L.)

Abstract As an economically important legume, Vicia sativa is commonly used as fodder for livestock and in health‐promoting foods for human consumption; however, cold stress is one of the most important limitations to its utilization at high altitudes. The shoot tissues of three cold‐resistant (Lan1, Lan2, and Lan3) and three cold‐sensitive (368, 521, and 538) accessions were collected after cold (4°C) treatment for 0, 6, and 48 h to assess their molecular regulatory mechanisms during cold stress via integrative analyses of their transcriptomes and metabolomes. Our Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that the biosynthesis of phenylpropanoids was the most enriched pathway in the cold stress response of V. sativa. A total of 18 differential metabolites were obtained in the biosynthesis of the phenylpropanoid pathway, and 722 differentially expressed genes (DEGs) showed a significant positive correlation with naringenin and p‐coumaric acid following a weighted gene co‐expression network analysis, while 84 DEGs showed significant negative associations. Moreover, four DEGs were overexpressed in yeast, and it was indicated that their biological function is to confer cold tolerance. These results provide a valuable resource for understanding the underlying molecular mechanism of V. sativa's response to cold stress and furthering the development of breeding research.