Transcriptomics integrated with metabolomics provides a new strategy for mining key genes in response to low temperature stress in Helictotrichon virescens

H. virescens is a perennial herbaceous plant with highly tolerant to cold weather, but the key genes that respond to low temperature stress still remain unclear. Hence, RNA-seq was performed using leaves of H. virescens treated at 0 °C and 25 °C for 12 h, 36 h, and 60 h, respectively, and a total of 9416 DEGs were significantly enriched into seven KEGG pathways. The LC-QTRAP platform was performed using leaves of H. virescens leaves at 0 °C and 25 °C for 12 h, 36 h, and 60 h, respectively, and a total of 1075 metabolites were detected, which were divided into 10 categories. Additionally, 18 major metabolites, two key pathways, and six key genes were mined using a multi-omics analytical strategy. The RT-PCR results showed that with the extension of treatment time, the expression levels of key genes in the treatment group gradually increased, and the difference between the treatment group and the control group was extremely significant. Notably, the functional verification results showed that the key genes positively regulated cold tolerance of H. virescens. These results can lay a foundation for the in-depth analysis of the mechanism of response of perennial herbs to low temperature stress.