Integrated Volatile Metabolome and Transcriptome Analyses Provide Insights into the Formation of Benzenoid–Phenylpropanoid Aroma Substance Eugenol in the Rosa hybrida ‘Lanxing’ Flowering

Abstract

Scent is the key character of the horticultural ornamental plant rose, and benzenoid–phenylpropanoid compounds are the main source of scent. However, the underlying biosynthesis mechanism of these benzenoid–phenylpropanoid scent metabolites during Rosa flowering is poorly understood. In this study, the volatile metabolome and transcriptome conjoint analysis was conducted on the six stages petals of the variety ‘Lanxing’ to investigate the synthesis of benzenoid–phenylpropanoid metabolites. A total of 25 benzenoid–phenylpropanoid volatile compounds were identified, of which eugenol possessed the highest content. Meanwhile, transcriptome analysis produced 87.9 million clean reads and 22,004 differentially expressed genes (DEGs). Group pairwise comparison of gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis found DEGs were enriched into phenylpropanoid compound synthesis related pathway. Weighted gene co-expression network analysis (WGCNA) found a MEgreenyellow gene module (650 DEGs) correlated with phenylpropanoid compounds. Based on the eugenol content variation and gene spatio-temporal expression, a key candidate gene RcEGS32 related to the synthesis of eugenol was identified. Co-expression network analysis found that five transcription factors, RcMYB1, RcBES1, RcERF2, RcbHLH1, and RcTUB, may act as regulators in the eugenol synthesis process by directly binding to RcEGS32 or forming a complex unit. This study provided key insights into the formation of the scent substance eugenol during flowering, offering a valuable volatile metabolome and transcriptome resource for the future target trait-related gene discovery of roses.