Integrated Analysis of Transcriptome and Metabolome and Evaluation of Antioxidant Activities in Lavandula pubescens

Lavandula pubescens, belonging to the Labiatae family, is a newly discovered strongly aromatic species of lavender that is potentially beneficial for human health. Given the economic importance of lavender species, we sought in this study to characterize the terpenoid biosynthesis of L. pubescens by obtaining transcriptomic and metabolic datasets. Transcriptome analysis of L. pubescens grown aseptically in tissue culture medium yielded 124,233 unigenes with an average length of 470 bp and N50 value of 522 bp from 9,476,122,928 raw reads. In order to provide relevant biological information, the unigenes were annotated using the following public databases: National Center for Biotechnology Information (NCBI) nucleotide (NT) and non-redundant protein (NR), Brassica (BRAD), Arabidopsis Information Resource (TAIR), Clusters of Orthologous Groups (COG), and Gene Ontology (GO). NR annotation results revealed that L. pubescens is genetically closely related to Sesamum indicum. On the basis of the transcriptome data, a total of 14 cDNA clones encoding the terpene biosynthetic genes LpDXS, LpMCT, LpMCS, LpHDR, LpIDI, LpAACT, LpHMGS, LpHMGR, LpMVK, LpPMK, LpMVD, LpGPPS, LpSQS, and LpGGPPS were identified in L. pubescens. These were quantified in the roots, stems, and leaves of L. pubescens using quantitative real-time polymerase chain reaction (qRT-PCR), which revealed that the gene expression levels were higher in the leaves and stems than in the roots, which was found to be consistent with the levels of ursolic and oleanolic acids in the different organs using high-performance liquid chromatography (HPLC). A total of 48 hydrophilic metabolites were identified and quantified in the organs using gas chromatography time-of-flight mass spectrometry (GC-TOFMS). Furthermore, the antioxidant activity of an ethyl acetate extract of L. pubescens leaves was examined using different methods to determine the potential therapeutic properties. A reducing power assay revealed that the absorbance values increased in a concentration-dependent manner, whereas a 2,2-diphenyl-1-picrylhydrazyl radical scavenging assay indicated the strong activity (60.4 ± 0.9%) of the ethyl acetate extract at a concentration of 100 µg/mL, which also showed strong hydrogen peroxide (57.4 ± 2.7%), superoxide radical (62.1 ± 0.7%), and hydroxyl radical (58.6 ± 0.4%) scavenging activities.