Skip to main content

Identification of key genes controlling monoterpene biosynthesis of Citral-type Cinnamomum bodinieri Levl. Based on transcriptome and metabolite profiling

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.

Peer Review reports

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:

$${\rm{Y}}\,{\rm{ = }}\,{\rm{(}}{{\rm{M}}_{\rm{1}}}\,{\rm{/}}\,{{\rm{M}}_{\rm{2}}}{\rm{)}}\,{\rm{ \times }}\,{\rm{100}}\,{\rm{\% }}$$

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:

$${\rm{RI}}\,{\rm{ = }}\,{\rm{100}}\left[ {{{{\rm{lo}}{{\rm{g}}_{{\rm{10}}}}{X_{\rm{i}}}\, - \,{\rm{lo}}{{\rm{g}}_{{\rm{10}}}}{X_{\rm{n}}}} \over {{\rm{lo}}{{\rm{g}}_{{\rm{10}}}}{X_{{\rm{n + 1}}}}\, - \,{\rm{lo}}{{\rm{g}}_{{\rm{10}}}}{X_{\rm{n}}}}}\,{\rm{ + }}\,n} \right]$$

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).

Fig. 1
figure 1

Citral-type and non-citral type samples of C. bodinieri. (a) Two chemotypes of C. bodinieri leaves. (b) Two chemotypes of C. bodinieri leaves EOs

Table 1 Essential oils content of the C. bodinieri leaves (%)

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.

Table 2 Essential oils composition in leaves of the C. bodinieri

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).

Fig. 2
figure 2

Transcriptomic analysis. (a) Principal component analysis (PCA) score plots. (b) Venn diagram of DEGs. (c) Volcano plots of differentially expressed genes (DEGs)

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).

Fig. 3
figure 3

KEGG enrichent scatter plot of DEGs. (a) B0 vs. B1. (b) B0 vs. B2. (c) B0 vsB3

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.

Fig. 4
figure 4

Differential expression of genes involved in plant-pathogen interaction pathway in citral-type and noncitral-type of C. bodinieri, designated as B1, B2, B3 and B0, respectively. Expression values were presented as (FPPM + 1) normalized log2-transformed counts

Fig. 5
figure 5

Differential expression of genes involved in MAPK signaling pathway-plant pathway in citral-type and noncitral-type of C. bodinieri

Fig. 6
figure 6

Differential expression of genes involved in plant hormone signal transduction pathway in citral-type and noncitral-type of C. bodinieri

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.

Fig. 7
figure 7

Differential expression of genes involved in starch and sucrose metabolismpathway in citral-type and noncitral-type of 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).

Fig. 8
figure 8

Exploration of the monoterpene biosynthesis pathway in C. bodinieri. (a) Citral synthesis pathway diagram. (b) Camphor synthesis pathway diagram. (c) Heat map of expression of key genes for the formation of citral and camphor

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).

Fig. 9
figure 9

The expression of 15 genes related to citral and camphor formation was determined by qRT-PCR. Note: R mean pearson correlation coefficients between the relative expression and FPKM

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).

Fig. 10
figure 10

Enzyme activity of the geraniol synthase (CbGES) and alcohol dehydrogenase (CbADH) in the citral-type (B1, B2, B3) and non-citral type (B0) by ELISA (p < 0.01)

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

  1. 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.

    Article  CAS  PubMed  Google Scholar 

  2. 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.

    Article  CAS  Google Scholar 

  3. 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.

    Article  CAS  PubMed  Google Scholar 

  4. 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.

    Article  Google Scholar 

  5. 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.

    Article  Google Scholar 

  6. 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.

    Article  Google Scholar 

  7. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. 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.

    Article  Google Scholar 

  9. 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.

    Article  CAS  Google Scholar 

  10. 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.

    Article  CAS  Google Scholar 

  11. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 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.

  13. 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.

    Article  Google Scholar 

  14. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 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.

    Article  Google Scholar 

  16. Croteau Rodney. Biosynthesis and catabolism of monoterpenoids. Chem Rev. 1987;87:929–54.

    Article  Google Scholar 

  17. 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.

    Article  CAS  PubMed  Google Scholar 

  18. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. 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.

    CAS  PubMed  Google Scholar 

  22. 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.

  23. 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.

    Article  CAS  Google Scholar 

  24. 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.

    Article  CAS  Google Scholar 

  25. 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.

    Article  CAS  PubMed  Google Scholar 

  26. 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.

    Article  CAS  Google Scholar 

  27. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. 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.

    Google Scholar 

  29. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 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.

    Article  PubMed  Google Scholar 

  31. 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.

  32. 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.

    Article  CAS  Google Scholar 

  33. Zuo Z, Wang B, Ying B, Zhou L, Zhang R. Monoterpene emissions contribute to thermotolerance in Cinnamomum camphora. Trees. 2017;31:1759–71.

    Article  CAS  Google Scholar 

  34. 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.

    Article  CAS  PubMed  Google Scholar 

  35. 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.

    Article  CAS  Google Scholar 

  36. Southwell I. Backhousia Citriodora F. Muell. (Lemon Myrtle), an unrivalled source of Citral. Foods. 2021;10:1596–609.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Sangwan NS, Farooqi AHA, Sangwan RS. Effect of drought stress on growth and essential oil metabolism in lemongrasses. New Phytol. 1994;128:173–9.

    Article  Google Scholar 

  38. 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.

    Article  CAS  PubMed  Google Scholar 

  39. 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.

    Article  CAS  Google Scholar 

  40. 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.

    Article  CAS  Google Scholar 

  41. 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.

    Article  CAS  PubMed  Google Scholar 

  42. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  43. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. 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.

    Article  CAS  PubMed  Google Scholar 

  45. 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.

    Article  CAS  Google Scholar 

  46. 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.

    Article  PubMed  Google Scholar 

  47. 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.

    Article  CAS  PubMed  Google Scholar 

  48. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  49. 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.

    Article  CAS  PubMed  Google Scholar 

  50. 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.

    Article  CAS  Google Scholar 

Download references

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

Authors

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

Correspondence to Zhinong Jin or Yuanqiu Liu.

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

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10419-7

Keywords