- Research
- Open access
- Published:
Identification of key genes controlling monoterpene biosynthesis of Citral-type Cinnamomum bodinieri Levl. Based on transcriptome and metabolite profiling
BMC Genomics volume 25, Article number: 540 (2024)
Abstract
The citral-type is the most common chemotype in Cinnamomum bodinieri Levl (C. bodinieri), which has been widely used in the daily necessities, cosmetics, biomedicine, and aromatic areas due to their high citral content. Despite of this economic prospect, the possible gene-regulatory roles of citral biosynthesis in the same geographic environment remains unknown. In this study, the essential oils (EOs) of three citral type (B1, B2, B3) and one non-citral type (B0) varieties of C. bodinieri were identified by GC-MS after hydrodistillation extraction in July. 43 components more than 0.10% were identified in the EOs, mainly composed of monoterpenes (75.8–91.84%), and high content citral (80.63–86.33%) were identified in citral-type. Combined transcriptome and metabolite profiling analysis, plant-pathogen interaction(ko04626), MAPK signaling pathway-plant(ko04016), starch and sucrose metabolism(ko00500), plant hormone signal transduction(ko04075), terpenoid backbone biosynthesis (ko00900) and monoterpenoid biosynthesis (ko00902) pathways were enriched significantly. The gene expression of differential genes were linked to the monoterpene content, and the geraniol synthase (CbGES), alcohol dehydrogenase (CbADH), geraniol 8-hydroxylase-like (CbCYP76B6-like) and 8-hydroxygeraniol dehydrogenase (Cb10HGO) were upregulated in the citral-type, indicating that they were associated with high content of geraniol and citral. The activities of CbGES and CbADH in citral type were higher than in non-citral type, which was corroborated by enzyme-linked immunosorbent assay (ELISA). This study on the accumulation mechanism of citral provides a theoretical basis for the development of essential oil of C. bodinieri.
Introduction
The citral is very beloved worldwide for the sweet lemon fragrance, which is a food additive permitted and a raw material for the synthesis of violetone, widely used in the flavour, fragrance and perfume industry [1, 2]. Moreover, citral has anti-inflammatory, bactericidal, antioxidant, insecticidal activities and cardio protective potential properties [3,4,5,6,7,8,9]. Woody plants have large biomass generally, and it is significant if we can find tree species with high citral content. Natural citral was used to be extracted from the fruits of Litsea cubeba in China, but its fruit was particularly troublesome to pick and the fruit ripening season was also short [10]. The branches and leaves of C. bodinieri were found to be rich in citral, which will likely be a new source of natural citral. Unfortunately, C. bodinieri was divided into 13 types based on the principal components of the EOs, such as borneol-rich, camphor-rich, cadinol-rich, cymol-rich, nerolidol-rich, methyleugenol-rich, safrole-rich, citral-rich, and so on. Among these chemotypes, the camphor-rich and citral-rich are the two most common chemotype in C. bodinieri [11]. Why does same morphology C. bodinieri in the same geographical environment generate different essential oil compositions? A putative inherent genetic effect and differentiation may be existed; therefore, a comprehensive research of the mechanism of citral formation is highly desirable.
The research of C. bodinieri have focused on stress resistance [12], photosynthetic and fluorescence parameters [13], insect control [14], tissue culture and nursery technology [15]. Nevertheless, the formation mechanisms of the precious chemotype C. bodinieri EOs have rarely been studied. Citral is a monoterpenoid, and its synthesis pathways contain three steps. Firstly the terpenoid backbone were biosynthesized via mevalonate pathway and 2-C-Methyl-D-erythritol 4-phosphate pathway to generate isopentenyl-PP and dimethylallyl-PP. Secondly isopentenyl-PP and dimethylallyl-PP were catalyzed to geranyl-PP by geranyl diphosphate synthase [16]. Finally geranyl-PP were modified by different monoterpene synthases to produce various monoterpene compounds [17]. Geranial and neral were two isomers of citral, where does geranial and neral in C. bodinieri come from?
RNA-Seq technology was widely used in mining functional genes, probing active ingredient metabolic pathways and identifying key enzyme genes in Cinnamomum camphora (C. camphora) [18,19,20,21], C. burmannii [19]. Metabolites are the phenotypes in the central dogma of molecular biology, which can be regard as a fingerprint of the joint action of genes and the environment [22, 23]. In our study, essential oil of the leaves were extracted by hydrodistillation and indentified the main components by GC-MS. Transcriptomics and volatile metabonomics in the citral-type and non-citral type varieties of C. bodinieri were performed. In addition, the genes related to the biosynthesis of citral were corroborated by qRT-PCR and ELISA. The elucidation of the formation mechanism of citral would lay the theoretical foundation for improving the yield and quality of essential oil of citral-type C. bodinieri.
Materials and methods
Plant material and reagent
Three citral-type (B1, B2, B3) and one non-citral type (B0) varieties of 15-year-old C. bodinieri, were sourced from forestry in the Tianzhu County, Qiandong nan Miao and Dong Autonomous Prefecture, Guizhou province, China (Latitude: 26° 97′ N, Longitude: 109° 10′ E) in July 2022. The plants were authenticated by professor Zhinong Jin, and deposited in the gene bank of Nanchang institute of technology, and the voucher numbers were C. bodinieri-GZ/TZ/028 (B1); C. bodinieri-GZ/TZ/029 (B2); C. bodinieri-GZ/TZ/031 (B3); C. bodinieri-GZ/TZ/041 (B0). The citral-type C. bodinieri used in our study were validated as Jiangxi provincial improved tree varieties after specificity, consistency and stability tests (http://ly.jiangxi.gov.cn/art/2023/10/16/art_68735_4629779.html).
Essential oils isolation and analysis
Essential oils were extracted from the leaves of 15-year old cuttings of the same height of three citral-type (B1, B2, B3) and one non-citral type (B0) varieties in July 2022, three sample plants were randomly selected as replicates for each variety. 400 g leaves and 2 L distilled water were placed in a hydro-distilled modified Clevenger apparatus [24], then distillated for 1 h. Calculate the essential oils yield based on the following formula:
where: Y is the leaf oil yield (%); M1 is the weight of essential oils; M2 is the weight of leaves.
Analyses of EOs compounds were performed by GC-MS (Agilent 7890B-5975 C, USA), referring to our previous study [7]. The HP-5MS (30 m×250 μm×0.25 μm) column was used. The electron impact of the mass spectrometer was 70 eV, the scanning mass range was set at 50 ~ 650 m/z, the scanning rate was set at 0.5 times/s, the temperature of the conductor was 250℃, the temperature of the ion source was 230℃, the temperature of the quadrupole was 150℃, and the doubling voltage was 1,200 V. Helium was used as the carrier gas (the flow rate was 1.2 mL/min), and 0.1 µL was injected (the separation ratio was 20 : 1). The oven temperature program conditions were as follows: an initial temperature of 80℃ for 5 min with a solvent delay of 3 min, followed by a gradual increase to 120℃ at a rate of 2.5 ℃/min, where it was held for 1 min, and then to 240℃ at a rate of 20 ℃/min for 5 min. The total run time was 60 min. Essential oils were diluted with methanol (1%), filtered, and auto-sampled.
The literature retention index (RI) was determined from the NIST library (https://webbook.nist.gov/chemistry/). The experimental retention index (RI) was determined by analyzing the samples with n-alkane C7-C40 standards under the same chromatographic conditions by injecting the samples, calculating and comparing RI, identifying as the same compound when RI in the experiment and in the literature were within 30:
where: X is the retention time of the unknown compound; Xn and Xn + 1 are the retention times of the corresponding before and after n-alkanes (Xn < Xi < Xn + 1).
RNA sequencing and data analysis
Twelve RNA samples (B0, B1, B2, B3 groups with three replicates) were extracted by the RNeasy Plant Mini Kit (Qiagen, Germany). Sequencing libraries were generated using NEBNext® Ultra™ RNA Library Prep Kit and library quality was inspected by Agilent Bioanalyzer 2100, then sequenced using the Illumina platform, finally, sequence data was filtered to get clean data. Over 83% of filtered reads per sample were mapped to the C. camphora reference genome (GWHBGBX00000000) [25] (Table S1). Final unigenes were aligned to the NCBI non-redundant (Nr), KEGG pathway database, Gene ontology (GO), KOG database, and Swiss-Prot protein database. Feature Counts v1.6.2 was used to calculate Fragments Per Kilobase of Transcript Per Million Mapped Reads (FPKM). The thresholds for significantly differential expression were set at p ≤ 0.01 and |log2(FoldChange)| ≥ 1. The raw sequence data has been submitted to the NCBI (PRJNA970043). PCA (principal component analysis), Heatmap (Visualization map), Inter-sample correlation plot and OPLS-DA were performed by statistics function within R (www.r-project.org).
Candidate genes confirmation
To confirm candidate genes related to formation of citral, the DEGs related to terpenoid backbone biosynthesis and monoterpenoid biosynthesis pathways were selected for quantitative real-time PCR (qRT‑PCR) [26], and three replications were done for each validation. Amplification consisted in an initial denaturation step at 95℃ for 5 min, followed by 40 cycles of denaturating at 95℃ for 10 s, 60℃ annealing temperature for 30 s. The specificity of primer was verified through a melting curve analysis at the end of the qRT‑PCR (0.5℃ ramping for 10 s, from 60℃ to 95℃). Normalization of qRT-PCR expression analysis of target genes was achieved through internal reference gene ACTIN and the primers were designed using the NCBI Primer-BLAST designing tool (Table S2). Relative gene expression was calculated using the 2−△△ct method and the graphic was statistically analyzed by Origin(version 2023b, USA).
Enzyme viability assay
The geraniol synthase (CbGES) and alcohol dehydrogenase (CbADH) enzyme activities of twelve samples (B0, B1, B2, B3 groups with three replicates) were determined by an enzyme-linked immunosorbent assay (ELISA) kit specific for plant(Mlbio, Shanghai, China) according to the manufacturer’s instructions. Double antibody sandwich method was used in the kit assay plant CbGES and CbADH level in the sample. At first, the capture antibody was encapsulated on a solid phase carrier, and antigen was added to bind, and the re-added detection antibody binded to the antigen, turned into antibody-antigen-enzyme-antibody complex, added TMB substrate and became blue color at HRP enzyme-catalyzed, then the color transformed to yellow by the action of acid. Finally, the absorbance (OD) was measured at 450 nm with an enzyme standard meter, and the enzyme activity by comparing the OD of the samples to the standard curve was calculated.
Result
Essential oils yield of citral-type and non-citral type C. Bodinieri leaves
The citral-type and non-citral type leaves were the same in morphology. The citral-type leaves EOs were transparent, whereas the non-citral type leaves EOs appeared camphor crystals precipitated and yellowish (Fig. 1). The oil yield ranged from 0.35 to 0.89% (w/w FW) and from 0.82 to 1.92% (w/w DW), respectively. The EOs yields of citral-type leaves were higher than the non-citral type (p ≤ 0.01) (Table 1).
Essential oils composition of citral-type and non-citral type C. Bodinieri leaves
The GC-MS experiment identified B2 variety attained 94.17% (14 constituents) and B1 variety 93.87% (14 constituents), followed by B0 variety 92.30% (34 constituents), B3 variety 90.56% (15 constituents). Two geometric isomers of citral, geranial and neral were dominated in the EOs composition in leaves of the citral-type C. bodinieri (B1, B2 and B3). In addition, unsaturated citral Z-isocitral and E-isocitral were existed frequently. The four dominant components ranged from 81.65 to 88.23% (Table 2). In the non-citral type C. bodinieri (B0), camphor (46.12%) was the most abundant, and borneol (14.67%) followed. Oxygenated monoterpenes dominated in the chemical composition of B0-B3 C. bodinieri EOs with proportions of 66.90%, 91.64%, 91.80% and 87.31%, respectively.
Transcriptome analysis of citral-type and non-citral type C. Bodinieri leaves
To identify the genes involved in citral-type C. bodinieri monoterpene biosynthesis, RNA-Seq was performed on three citral-type (B1, B2, B3) and one non-citral type (B0) varieties. A total of 12 cDNA libraries were constructed, comprising 558.83 million raw reads. After quality control, 538.02 million high-quality clean reads, including 80.7 Gb of nucleotide sequences.Furthermore, the percentage of Q30 base in each sample was not less than 92.27%. Over 83.26% of filtered reads per sample were mapped to the C. camphora reference genome (Table S1). 5,560 unigenes were annotated in NR, 4,697 in Pfam, 4,357 in Swiss-Prot, 418 in TF, 4,285 in KEGG, 4,842 in GO, 5,555 in TrEMBL and 2,868 in KOG databases (Table S3).The DEGs based on Kyoto Encylopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) were analyzed (Figures S1 and S2).
Three citral-type (B1, B2, B3) and one non-citral type (B0) varieties were clearly separated in principal component analysis (PCA) score plots (Fig. 2a). In B0 Vs B1, B0 Vs B2 and B0 Vs B3 volcano plots, 8337 DEGs (1952 upregulated and 2731 downregulated), 3710 DEGs (1675 upregulated and 2662 downregulated), and 6202 DEGs (1718 upregulated and 2353 downregulated) were displayed respectively (Fig. 2c). Venn diagram showed that 2985 DEGs were identified between citral-type and non-citral type, which may be associated with lemon and camphor fragrance differences (Fig. 2b).
Screening of genes related to synthesis of monoterpenoids
KEGG enrichment scatter plot of DEGs showed that plant-pathogen interaction (ko04626), MAPK signaling pathway-plant (ko04016), starch and sucrose metabolism (ko00500) and plant hormone signal transduction (ko04075) were enriched significantly (Fig. 3).
The transcriptional expression of ko04626, ko04016 and ko04075 pathway genes were investigated by preparing heat maps (Figs. 4, 5 and 6).As expected, only several putative nucleoside-diphosphate kinase (ndk), catalase (CAT), serine/threonine-protein kinase SRK2 (SNRK2), EIN3-binding F-box protein (EBF1_2), transmembrane protein 222 (TMEM222), glycerol kinase(glpk), BRI1 kinase inhibitor 1 (BKI1), two-component response regulator ARR-A family (ARR-A), ABA responsive element binding factor (ABF) and phenylpyruvate tautomerase (MIF) genes were up-regulated in B1, B2 and B3 relative to B0, respectively. Most other genes were signifcantly down-regulated in B1, B2 and B3 relative to B0, respectively (Figs. 4–6, Table S4-S6). It is indicated that most genes in plant-pathogen interaction (ko04626), MAPK signaling pathway-plant (ko04016) and plant hormone signal transduction (ko04075) pathways were actively expressed in camphor-type C. bodinieri, which might be related to the stronger anthelmintic and antibacterial properties of the camphor.
The transcriptional expression of ko00500 pathway genes was also investigated by preparing heat maps (Fig. 7). The starch synthase (glgA), UTP–glucose-1-phosphate uridylyltransferase (galU) and ectonucleotide pyrophosphatase/phosphodiesterase family member 1/3 (ENPP1_3) were signifcantly up-regulated in B1, B2 and B3 relative to B0, respectively. Five trehalose 6-phosphate synthase/phosphatase (TPS), glucan endo-1,3-beta-glucosidase 1/2/3 (GN1_2_3), two glucose-1-phosphate adenylyltransferase (glgC) and seven alpha-amylase (amyA) were signifcantly down-regulated in B1, B2 and B3 relative to B0, respectively (Fig. 7, Table S7). We hypothesized that the active expression of these genes involved in starch and sugar metabolism were related to the need for carbohydrates for secondary metabolism in C. bodinieri.
However, monoterpene was the largest proportion of C. bodinieri essential oils composition (Table 2), so we focused on genes related to monoterpene synthesis when screening for related genes. Four monoterpenoid synthesis pathways were further screened, terpenoid backbone biosynthesis (ko00900), monoterpenoid biosynthesis (ko00902), pinene, camphor and geraniol degradation (ko00907) and sesquiterpenoid and triterpenoid biosynthesis (ko00909).
Geraniol and nerol were catalyzed to geranial and neral by geraniol dehydrogenase (NAD+) or alcohol dehydrogenase(NADP+), respectively. Next, geranial was catalyzed to geranic acid by geranial dehydrogenase (Fig. 8a). Additionally, geranyl diphosphate (GPP) was cycled and hydrolyzed to borneol, and then borneol was oxidized to camphor by borneol dehydrogenase (BDH) (Fig. 8b).
All monoterpene production was dependent on terpenoid backbone biosynthesis, therefore, genes encoding acetyl-CoA C-acetyltransferase (CbAACT), 1-deoxy-D-xylulose-5-phosphate synthase (CbDXS), 1-deoxy-D-xylulose-5-phosphate reductoisomerase (CbDXR), iphosphomevalonate decarboxylase (CbMVD), isopentenyl phosphate kinase (CbIPK), isopentenyl-diphosphate Delta-isomerase (CbIDI), geranyl diphosphate synthase (CbGPS), geranylgeranyl diphosphate synthase, type II (CbGPPS) and farnesyl iphosphate synthase (CbFDPS) were identified. These key genes in MEP and MVA pathways were expressed actively, with high expression in both citral-type and non-citral type of C. bodinieri, favouring the synthesis of monoterpene synthetic precursors (Fig. 8c). In the monoterpenoid biosynthesis pathway, alpha-terpineol synthase (CbATS), (3 S)-linalool synthase (CbTPS14), neomenthol dehydrogenase (CbNDS), geraniol 8-hydroxylase-like (CbCYP76B6-like), 8-hydroxygeraniol dehydrogenase(Cb10HGO) and borneol diphosphate synthase (CbBDH) were identified. Genes encoding of the geraniol metabolic pathway was identified including geraniol synthase (CbGES) and alcohol dehydrogenase (CbADH) (Fig. 8c).
Among these genes, the expression of CbAACT, CbDXS, CbDXR, CbMVD, CbIPK, CbIDI, CbIPPI, CbGPPS, CbFDPS, CbGES, CbADH, CbCYP76B6-like, Cb10HGO and CbATS were significantly higher in the citral-type. The neomenthol dehydrogenase (CbNDS) and borneol diphosphate synthase (CbBDH) showed significantly higher expression in non-citral type in comparison to citral-type leaf expression.
Confirmation of candidate genes related to formation of citral via qRT‑PCR
We further selected 8 genes in terpenoid backbone biosynthesis pathway (CbAACT, CbDXS, CbDXR, CbMVD, CbIPK, CbFDPS, CbGPPS and CbGPS), 7 genes in monoterpenoid biosynthesis pathway (CbATS, CbCYP76B6-like, Cb10HGO, CbGES, CbADH, CbNDS and CbBDH), to validate the transcriptomic results. The expression of these unigenes observed by qRT-PCR were strongly supported the RNA-Seq value (Fig. 9). The expression of CbAACT, CbDXS, CbDXR, CbMVD, CbIPK, CbFDPS, and CbGPPS in citral-type was significantly higher than that in non-citral type. This result was corroborated with the high content of monoterpenes (hydrocarbon and oxygenated) in the citral-type (Table 2). The high expression of these genes provided abundant precursors for the synthesis of monoterpenes. In addition, the expression of CbGPS was lower in citral-type, whereas the CbGPS was a key gene of monoterpenes synthesis universal precursor geranyldiphosphate (GPP). Moreover, the expression of CbATS, CbCYP76B6-like, Cb10HGO, CbGES, CbADH in citral-type were significantly higher than that in non-citral type. Notably, the relative expression of CbNDS and CbBDH in non-citral type were similar to their FPKM (R = 0.998 and 0.988) (Fig. 9).
Enzyme activity of CbGES and CbADH in citral-type and non-citral type
Transcriptome analyses of C. bodinieri and qRT-PCR assays showed that citral content was related to the gene expression of CbGES and CbADH. An enzyme-linked immunosorbent assay (ELISA) can be used for quantitative measurement of protein. The geraniol synthase (CbGES) enzyme activities of three citral-type (B1, B2, B3) and one non-citral type (B0) C. bodinieri varieties were 163.44 ± 0.84, 169.44 ± 0.09, 174.28 ± 2.31 and 104.26 ± 2.84, respectively. The alcohol dehydrogenase (CbADH) enzyme activities of B1, B2, B3 and B0 were 169.21 ± 1.76, 210.21 ± 2.02, 188.09 ± 0.50 and 136.38 ± 2.44, respectively. The enzyme activity of the geraniol synthase (CbGES) and alcohol dehydrogenase (CbADH) in the citral-type were higher than in non-citral type significantly by ELISA (p < 0.01) (Fig. 10).
Discussion
The C. bodinieri is suitable for large-scale cultivation and convenient for harvesting in short rotation period, and rich in EOs, especially contain more than 80% citral in the citral-type, therefore, this plant has the potential to be an unrivalled source of natural citral [7]. The transcriptome is an effective tool for the study of terpene metabolite synthesis mechanisms, and the publication of C. camphora genome has provided good conditions for the study of the molecular mechanisms of Lauraceae plants [25, 27,28,29,30,31], and the reference transcriptome of C. bodinieri were first systemically investigated. The possible genes of citral formation was studied base on reference transcriptomics and metabolomics.
Essential oils yield and composition of citral-type C. Bodinieri
While the leaves oil yield of citral-type C. bodinieri were 0.80–0.89%, the non-citral type were 0.35 ± 0.18%, indicating different chemotypes had significant differences (p ≤ 0.01) and the same chemotype with no significant difference. Similar field test was done in C. bodinieri, presented a different result that the non-citral type (1.18%) was higher than citral-type (0.26%) significantly [11]. In terms of the extract and the main components, there may be a significant difference among individuals, which were related to factors such as environmental, genetic, geographical, seasonal and extraction methods [10, 32, 33]. Therefore, we collected leaves from camphor trees grown in the same environment at the same time to ensure that the differences in results were only due to genetic factors. Metabolites were the basis of an organism’s phenotype and can contribute to understand biological processes and their mechanisms intuitively. We identified an interesting finding, monoterpene volatiles play important roles in aroma generation of C. bodinieri, forty-three constituents were identified comprising more than 90% of the composition of oils. The different chemical types of C. bodinieri do not differ morphologically, but differ in their main components. High content neral (32.57–37.13%) and geranial (48.06–49.15%) were identified in the essential oil of three citral-type (B1, B2, B3), which were higher than the rest of citral-rich plants in the world such as Litsea cubeba [10], Cymbopogon citratus [34], Ocimum gratissimum [35], Backhousia citriodra [36]. Those rare citral-type C. bodinieri varieties were identified through the scenting method from approximately 40,000 seedlings from Guizhou province, China. The citral-types were also high EOs yielding varieties, contained 80.63-86.33% citral and suitable for large-scale cultivation, therefore they will be expected to become new sources of natural citral. The study of the citral and camphor chemotype cannot only consider citral and camphor, but also the metabolic relationship of substances on the synthetic pathway of citral and camphor(Figure.3b, d) [37, 38]. Neral, geranial and geraniol were higher in citral-type, whereas geranic acid, (+)-borneol and (+)-camphor were higher in non-citral type of C. bodinieri (Table 2). These results illustrated the formation mechanism of citral and camphor of C. bodinieri at the level of substance metabolism.
Terpene biosynthesis genes related to the formation of citral
In our research, the monoterpenes in the citral-type (87.70–92.28%) were higher than the non-citral type (75.8%)(Table 2), therefore, the high expression values (log(FPKM + 1)) of key genes CbAACT, CbDXS, CbDXR, CbMVD, CbIPK, CbIDI, CbIPPI, CbGPPS and CbFDPS of terpenoid backbone biosynthesis in citral-type were reasonable. The higher expression of these genes provided more precursors for downstream biosynthesis of monoterpenes, which were consistent with previous transcriptome results of C. camphora [18, 20], Curcuma wenyujin [39], Litsea cubeba [40], tree peony [41]. In many plants, geraniol synthase (GES) could catalyze geranyl pyrophosphate to generate geraniol [37, 42]. Besides, the biosynthesis of geraniol in rose is RhNUDX1 gene, a new member of the Nudix hydrolase family [43, 44]. Geraniol and nerol are synthesized by the action of geraniol dehydrogenase or ethanol dehydrogenase [37, 45]. In our study, Geraniol dehydrogenase (GDH) was not identiffed from the transcriptome database, possibly due to incomplete annotation of structural genes in the reference genome, however the putative alcohol dehydrogenase (CbADH) gene expression in citral type were higher than the non-citral type. The enzyme activities of CbGES and CbADH in the citral-type were significantly higher than in the non-citral type by ELISA. Our research group have cloned the cDNA of the geraniol synthase from C. camphora (CcGES) and performed functional characterization of the protein encoded by CcGES. The terpene genes, such as cytochrome P450s (CbCYP76B6 and Cb10HGO), terpene synthase (CbATS) and geranylgeranyl pyrophosphate synthase (CbGPPS), were significantly selected in the citral-type group but not in the non-citral type group. The cytochrome P450s is a class of enzymes necessary for the functionalization of monoterpenes, among which CYP76B6 catalyze the hydroxylation of geraniol at their C-10 position to generate 8-oxogeraniol [46]. The 10HGO was a key gene of iridoid glycoside synthesis, oxidizing 10-oxogeraniol to 10-oxogeranial, which was a NAPDH-dependent cytochrome P450 monoterpene oxidase [47, 48]. In our study, enhanced CbCYP76B6 and Cb10HGO activity due to increased substrate as a result of increased geraniol content. Borneol, the precursor of camphor, was higher in non-citral type than in citral-type. Then it generated camphor by the catalytic of borneol dehydrogenase (BDH). Borneol diphosphate synthase (CbBDH) gene expression in the non-citral type were significantly higher than in the citral-types, which was associated with the synthesis of camphor [31, 38, 49, 50].
Although further analyses are needed to verificated by functional characterization of the protein encoded and genetic transformation experiments, this study nevertheless provides the molecular evidence for the CbGES and CbADH strong correlation with citral formation, as well as the CbBDH related to camphor formation in C. bodinieri.
Conclusion
In recent years, high throughput sequencing and metabolome have provided several breakthroughs in aromatic plants. This work analyzed transcriptome and metabolite profiling of three citral-types and one non-citral type C. bodinieri. Monoterpenes were identified as main components in the leaves essential oils, neral and geranial were the most abundant components in the citral-type, while (+)-borneol and (+)-camphor were the major component in the non-citral type. Fifteen candidate unigenes about terpenoid biosynthesis were identified and validated. Among these genes, the CbGES and CbADH were correlated with citral formation, and the CbBDH was related to camphor formation in C. bodinieri, however, further functional studies are needed. This first reference transcriptome sequence of C. bodinieri database can contribute to understand the biosynthesis of citral in Lauraceae plants.
Data availability
Data supporting the findings of this work were available within the paper and its supplementary information files. Raw reads yielded from Illumina sequencing have been uploaded to the NCBI Sequence Read Archive (https://submit.ncbi.nlm.nih.gov/subs/sra/) and accession numbers for the twelve samples are as below: B1-1 (SAMN34992896); B1-2 (SAMN34992897); B1-3 (SAMN34992898); B2-1 (SAMN34992899); B2-2 (SAMN34992900); B2-3 (SAMN34992901); B3-1 (SAMN34992902); B3-2 (SAMN34992903); B3-3 (SAMN34992904); B0-1 (SAMN34992905); B0-2 (SAMN34992906); B0-3 (SAMN34992907).
References
Hirai M, Ota Y, Ito M. Diversity in principal constituents of plants with a lemony scent and the predominance of citral. J Nat Med. 2022;76:254–8.
Zhou F, Liu H, Wen Z, Zhang B, Chen G. Toward the efficient synthesis of pseudoionone from Citral in a continuous-Flow Microreactor. Ind Eng Chem Res. 2018;57:11288–98.
Chaouki W, Leger DY, Liagre B, Beneytout J, Hmamouchi M. Citral inhibits cell proliferation and induces apoptosis and cell cycle arrest in MCF-7 cells. Fundam Clin Pharmacol. 2009;23:549–56.
Da Silva Júnior AQ, Da Silva DS, Figueiredo PLB, Sarrazin SLF, Bouillet LEM, De Oliveira RB, et al. Seasonal and circadian evaluation of a citral-chemotype from Lippia alba essential oil displaying antibacterial activity. Biochem Syst Ecol. 2019;85:35–42.
Leite MCA, Bezerra APDB, Sousa JPD, Guerra FQS, Lima EDO. Evaluation of antifungal activity and mechanism of action of Citral against Candida albicans. Evidence-Based Complement Altern Med. 2014;2014:1–9.
Pei C, Tsung S, Ju C, Jie C, Shu C, Yueh H, et al. Anti-inflammatory activity of neral and geranial isolated from fruits of Litsea cubeba Lour. J Funct Foods. 2015;19:248–58.
Ling Q, Zhang B, Wang Y, Xiao Z, Hou J, Xiao C, et al. Chemical composition and antioxidant activity of the essential oils of Citral-Rich Chemotype Cinnamomum camphora and Cinnamomum bodinieri. Molecules. 2022;27:7356.
Sharma S, Habib S, Sahu D, Gupta J. Chemical Properties and therapeutic potential of Citral, a Monoterpene isolated from Lemongrass. MC. 2020;17:2–12.
Silva GDSE, Marques JNDJ, Linhares EPM, Bonora CM, Costa ÉT, Saraiva MF. Review of anticancer activity of monoterpenoids: Geraniol, nerol, geranial and neral. Chemico-Biol Interact. 2022;362:109994.
Fan G, Ning X, Chen S, Zhong L, Guo C, Yang Y, et al. Differences in fruit yields and essential oil contents and composition among natural provenances of Litsea cubeba in China and their relationships with main habitat factors. Ind Crops Prod. 2023;194:116285.
Fu C, Liu X, Liu Q, Qiu F, Yan J, Zhang Y, et al. Variations in essential oils from the leaves of Cinnamomum bodinieri in China. Molecules. 2023;28:3659.
Han H, Zhang L, Li S, Zhao R, Wang F, Dong R et al. Transcriptome and Metabolome Integrated Analysis reveals the mechanism of Cinnamomum bodinieri root response to alkali stress. 2023. https://doi.org/10.21203/rs.3.rs-2487448/v1.
Liu C, Liu Y, Guo K, Li G, Zheng Y, Yu L, et al. Comparative ecophysiological responses to drought of two shrub and four tree species from karst habitats of southwestern China. Trees. 2011;25:537–49.
Zhou J, Zhang J, Tregenza T, Pan Y, Wang Q, Shi H, et al. Larval host preference and suitability for the Sawfly Mesoneura rufonota among five Cinnamomun Tree species. Insects. 2020;11:76.
Xiao Z. Effects of IBA on rooting ability of Cinnamomum Bodinieri citral type micro-shoots from transcriptomics analysis. Plant Biotechnol Rep. 2020;14:467–77.
Croteau Rodney. Biosynthesis and catabolism of monoterpenoids. Chem Rev. 1987;87:929–54.
Chen F, Tholl D, Bohlmann J, Pichersky E. The family of terpene synthases in plants: a mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom: terpene synthase family. Plant J. 2011;66:212–29.
Chen C, Zheng Y, Zhong Y, Wu Y, Li Z, Xu L, et al. Transcriptome analysis and identification of genes related to terpenoid biosynthesis in Cinnamomum camphora. BMC Genomics. 2018;19:550.
Hou C, Zhang Q, Xie P, Lian H, Wang Y, Liang D, et al. Full-length transcriptome sequencing reveals the molecular mechanism of monoterpene and sesquiterpene biosynthesis in Cinnamomum burmannii. Front Genet. 2023;13:1087495.
Hou J, Zhang J, Zhang B, Jin X, Zhang H, Jin Z. Transcriptional Analysis of Metabolic Pathways and Regulatory Mechanisms of Essential Oil Biosynthesis in the leaves of Cinnamomum camphora (L.) Presl. Front Genet. 2020;11:598714.
Jiang X, Wu W, Xiao F, Xiong Z-Y, Xu H. Transcriptome analysis for leaves of five chemical types in Cinnamomum camphora. Hereditas. 2014;36:58–68.
Zhang H, Chen M, Wen H, Wang Z, Chen J. Transcriptomic and metabolomic analyses provide insight into the volatile compounds of citrus leaves and flowers. BMC Plant Biol. 2020;20.
Zhong Y, Chen C, Gong xue, Luan X, Wu zhaoxiang, Li huihu, et al. Transcriptome and metabolome analyses reveal a key role of the anthocyanin biosynthetic pathway cascade in the pigmentation of a Cinnamomum camphora red bark mutant (‘Gantong 1’). Industrial Crops Prod. 2022;175:114236.
Zhang B, Xiao Z, Zhang H, Cao M, Liu Y, Jin Z. Chemical constituents and chemotypes of Fresh Leaf essential oil of wild species belonging to Sect. Camphor (Trew.) Meissn. In Southeastern China. J Essent Oil Bearing Plants. 2019;22:1115–22.
Shen T, Qi H, Luan X, Xu W, Yu F, Zhong Y, et al. The chromosome-level genome sequence of the camphor tree provides insights into Lauraceae evolution and terpene biosynthesis. Plant Biotechnol J. 2022;20:244–6.
Wang M, Jiao Y, Zhao Y, Gao ming, Wang Y. Phytohormone and transcriptome of pericarp reveals jasmonate and LcMYC2 are involved in neral and geranial biosynthesis in Litsea cubeba. Ind Crops Prod. 2022;177:114423.
Chaw S, Liu Y, Wu Y, Wang H, Lin CI, Wu C, et al. Stout Camphor tree genome fills gaps in understanding of flowering plant genome evolution. Nat Plants. 2019;5:63–73.
Chen J, Mao Y, Wang LA, Chromosome-Level. Genome of the Camphor Tree and the underlying genetic and climatic factors for its top-geoherbalism. Front Plant Sci. 2022;13:17.
Jiang R, Chen X, Liao X, Peng D, Han X, Zhu C, et al. A chromosome-level genome of the Camphor Tree and the underlying genetic and climatic factors for its top-geoherbalism. Front Plant Sci. 2022;13:827890.
Sun W, Xiang S, Zhang Q, Xiao L, Zhang D, Zhang P, et al. The camphor tree genome enhances the understanding of magnoliid evolution. J Genet Genomics. 2021;49:249–53.
Wang X, Xu C, Zheng Y, Wu Y, Zhang Y, Zhang T et al. Chromosome-level genome assembly and resequencing of camphor tree (Cinnamomum camphora) provides insight into phylogeny and diversification of terpenoid and triglyceride biosynthesis of Cinnamomum. Hortic Res. 2022;:216–74.
Zhang T, Yang H, Wen S, Qiu F, Liu X. Effects of Harvest season and Storage Time on the essential oil of the Linalool Chemotype of Cinnamomum camphora. J Essent Oil Bearing Plants. 2019;22:1379–85.
Zuo Z, Wang B, Ying B, Zhou L, Zhang R. Monoterpene emissions contribute to thermotolerance in Cinnamomum camphora. Trees. 2017;31:1759–71.
Kaur G, Arya SK, Singh B, Singh S, Sushmita, Saxena G, et al. Comparative transcriptional analysis of metabolic pathways and mechanisms regulating essential oil biosynthesis in four elite Cymbopogon spp. Int J Biol Macromol. 2023;229:943–51.
Kumar A, Lal RK, Chanotiya CS, Dwivedi A. The pre-eminence of agro-parameters and chemical constituents in the influence of harvest interval by traits × environment interaction over the years in lemon-scented basil (Ocimum Africanum Lour). Ind Crops Prod. 2021;172:113989.
Southwell I. Backhousia Citriodora F. Muell. (Lemon Myrtle), an unrivalled source of Citral. Foods. 2021;10:1596–609.
Sangwan NS, Farooqi AHA, Sangwan RS. Effect of drought stress on growth and essential oil metabolism in lemongrasses. New Phytol. 1994;128:173–9.
Sarker LS, Galata M, Demissie ZA, Mahmoud SS. Molecular cloning and functional characterization of borneol dehydrogenase from the glandular trichomes of Lavandula x intermedia. Arch Biochem Biophys. 2012;528:163–70.
Jiang C, Fei X, Pan X, Huang H, Qi Y, Wang X, et al. Tissue-specific transcriptome and metabolome analyses reveal a gene module regulating the terpenoid biosynthesis in Curcuma Wenyujin. Ind Crops Prod. 2021;170:113758.
Zhao Y, Chen Y, Gao M, Yin H, Wu L, Wang Y. Overexpression of geranyl diphosphate synthase small subunit 1 (LcGPPS.SSU1) enhances the monoterpene content and biomass. Ind Crops Prod. 2020;143:111926.
Li R, Li Z, Leng P, Hu Z, Wu J, Dou D. Transcriptome sequencing reveals terpene biosynthesis pathway genes accounting for volatile terpene of tree peony. Planta. 2021;254:67.
Wang M, Gao M, Zhao Y, Chen Y, Wu L, Yin H, et al. LcERF19, an AP2/ERF transcription factor from Litsea cubeba, positively regulates geranial and neral biosynthesis. Hortic Res. 2022;9:uhac093.
Conart C, Bomzan DP, Huang X, Bassard J, Paramita SN, Saint Marcoux D, et al. A cytosolic bifunctional geranyl/farnesyl diphosphate synthase provides MVA-derived GPP for geraniol biosynthesis in rose flowers. Proc Natl Acad Sci USA. 2023;120:e2221440120.
Magnard J, Roccia A, Caissard J, Vergne P, Sun P, Hecquet R, et al. Biosynthesis of monoterpene scent compounds in roses. Science. 2015;349:81–3.
Iijima Y, Koeduka T, Suzuki H, Kubota K. Biosynthesis of geranial, a potent aroma compound in ginger rhizome (Zingiber officinale): molecular cloning and characterization of geraniol dehydrogenase. Plant Biotechnol. 2014;31:525–34.
Höfer R, Dong L, André F, Ginglinger JF, Lugan R, Gavira C, et al. Geraniol hydroxylase and hydroxygeraniol oxidase activities of the CYP76 family of cytochrome P450 enzymes and potential for engineering the early steps of the (seco)iridoid pathway. Metab Eng. 2013;20:221–32.
Awadasseid A, Li W, Liu Z, Qiao C, Pang J, Zhang G, et al. Characterization of Camptotheca acuminata 10-hydroxygeraniol oxidoreductase and iridoid synthase and their application in biological preparation of nepetalactol in Escherichia coli featuring NADP+ - NADPH cofactors recycling. Int J Biol Macromol. 2020;162:1076–85.
Dong T, Song S, Wang Y, Yang R, Chen P, Su J, et al. Effects of 5-azaC on Iridoid Glycoside Accumulation and DNA methylation in Rehmannia Glutinosa. Front Plant Sci. 2022;13:913717.
Ma R, Su P, Jin B, Guo J, Tian M, Mao L, et al. Molecular cloning and functional identification of a high-efficiency (+)-borneol dehydrogenase from Cinnamomum camphora (L.) Presl. Plant Physiol Biochem. 2021;158:363–71.
Tian N, Tang Y, Xiong S, Tian D, Chen Y, Wu D, et al. Molecular cloning and functional identification of a novel borneol dehydrogenase from Artemisia annua L. Ind Crops Prod. 2015;77:190–5.
Funding
This work was supported by Jiangxi Science and Technology Department (grant numbers 20232BAB205055, 20203ABC28W016), Jiangxi provincial department of forestry (grant numbers [2020]07,[2019]21), and Jiangxi Education Department (grant number GJJ2201515).
Author information
Authors and Affiliations
Contributions
Q.L., B.Z. and Y.W. wrote the main manuscript text. Z.J. and Y.L. took responsibility for the main manuscript conceptualization and methodology. Z.X., J.H. and Q.L. prepared data curation. J.Z. and C.X. performed formal analysis and draw figures. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
All authors have read and agreed to the published version of the manuscript. The authors declare no conflict of interest.
Material licences
The citral-type C. bodinieri used in our study were validated as Jiangxi provincial improved tree varieties after specificity, consistency and stability tests. The improved tree varieties belonged to Nanchang Institute of Technology and the breeders were Dr Qingyan Ling and Beihong Zhang. Our collection of C. bodinieri was approved by Nanchang Institute of Technology.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Ling, Q., Zhang, B., Wang, Y. et al. Identification of key genes controlling monoterpene biosynthesis of Citral-type Cinnamomum bodinieri Levl. Based on transcriptome and metabolite profiling. BMC Genomics 25, 540 (2024). https://doi.org/10.1186/s12864-024-10419-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-024-10419-7