Abstract
Origanum vulgare, belonging to the Lamiaceae family, is a principal culinary herb used worldwide which possesses great antioxidant and antibacterial properties corresponding to various volatile organic components (VOCs). However, the metabolite profiles and underlying biosynthesis mechanisms of elaborate tissues (stems, leaves, bracts, sepals, petals) of Origanum vulgare have seldom been reported. Here, solid-phase microextraction-gas chromatography/mass spectrometry results showed that Origanum vulgare 'Hot and Spicy' (O. vulgare 'HS') was extremely rich in carvacrol and had the tissue dependence characteristic. Moreover, a full-length transcriptome analysis revealed carvacrol biosynthesis and its tissue-specific expression patterns of 'upstream' MVA/MEP pathway genes and 'downstream' modifier genes of TPSs, CYPs, and SDRs. Furthermore, the systems biology method of modular organization analysis was applied to cluster 16,341 differently expressed genes into nine modules and to identify significant carvacrol- and peltate glandular trichome-correlated modules. In terms of these positive and negative modules, weighted gene co-expression network analysis results showed that carvacrol biosynthetic pathway genes are highly co-expressed with TF genes, such as ZIPs and bHLHs, indicating their involvement in regulating the biosynthesis of carvacrol. Our findings shed light on the tissue specificity of VOC accumulation in O. vulgare 'HS' and identified key candidate genes for carvacrol biosynthesis, which would allow metabolic engineering and breeding of Origanum cultivars.