Next Article in Journal
Anti-Inflammatory Effect of Protopine through MAPK and NF-κB Signaling Regulation in HepG2 Cell
Previous Article in Journal
Metabolic Profiling and Transcriptional Analysis of Carotenoid Accumulation in a Red-Fleshed Mutant of Pummelo (Citrus grandis)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 27
Issue 14
10.3390/molecules27144599
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analyses Revealed the Key Metabolic Genes and Transcription Factors Involved in Terpenoid Biosynthesis in Sacred Lotus

by
Lili Qin
1,
Fei Du
1,
Ningning Yang
1,
Chen Zhang
1,
Zhiwen Wang
1,
Xingwen Zheng
2,
Jiawei Tang
1,
Liangbo Yang
2 and
Chen Dong
1,*
1
College of Biological Engineering, Henan University of Technology, Zhengzhou 450001, China
2
White Lotus Industrial Development Center of Guangchang County, Fuzhou 344900, China
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(14), 4599; https://doi.org/10.3390/molecules27144599
Submission received: 21 June 2022 / Revised: 18 July 2022 / Accepted: 18 July 2022 / Published: 19 July 2022
Browse Figures
Versions Notes

Abstract

:
As the largest group of structurally diverse metabolites, terpenoids are versatile natural compounds that act as metabolism mediators, plant volatiles, and ecological communicators. However, few terpenoid compounds have been identified in plant parts of sacred lotus (Nelumbo nucifera Gaertn.). To elucidate the molecular genetic basis of the terpene biosynthetic pathway, terpenes from different parts of the plant, including seeds (S), young leaves (YL), mature leaves (ML), white flowers (WF), yellow flowers (YF), and red flowers (RF), were identified by LC-MS/MS and the relative contents of the same terpenes in different parts were compared. The results indicate that all plant parts primarily consist of triterpenes, with only minor quantities of sesquiterpenes and diterpenes, and there were differences in the terpene content detected in different plant parts. To illustrate the biosynthesis of various terpenoids, RNA sequencing was performed to profile the transcriptomes of various plant parts, which generated a total of 126.95 GB clean data and assembled into 29,630 unigenes. Among these unigenes, 105 candidate unigenes are involved in the mevalonate (MVA) pathway, methyl-erythritol phosphate (MEP) pathway, terpenoid backbone biosynthesis pathway, and terpenoid synthases pathway. Moreover, the co-expression network between terpene synthase (TPS) and WRKY transcription factors provides new information for the terpene biosynthesis pathway.

1. Introduction

Sacred lotus belongs to Nelumbonaceae, which contains only two species, sacred lotus (Nelumbo nucifera Gaertn.) and American lotus (Nelumbo lutea Pers.) [1]. Sacred lotus is one of the economically best-known aquatic plants in East Asia, with a history of cultivation for nearly 2000 years in China, and is commonly used for medicine and ornamentation [2,3]. A large number of metabolites, including alkaloids, steroids, flavonoids, triterpenoids, glycosides, and polyphenolshave, were widely detected in sacred lotus and are tightly related to their pharmacological activities, such as anti-ischaemia, antioxidant, anticancer, antiviral, as well as anti-obesity [4,5,6,7,8,9].
As the largest group of structurally diverse metabolites, terpenoids are versatile natural compounds that act as metabolism mediators, plant volatiles, and ecological communicators. Triterpenoids are important plant secondary metabolites with pharmacological effects, such as anticancer, antiviral, and cholesterol lowering. According to the number of rings in the structure, triterpenoids can be classified into monocyclic triterpenes, bicyclic triterpenes, tricyclic triterpenes, tetracyclic triterpenes, and pentacyclic triterpenes, etc. Among them, pentacyclic triterpenoids have received more attention because of more types and functions [10,11,12,13]. Biologically, triterpene saponins are considered to be defensive compounds against external stresses [14,15]. Terpenoids are sequentially biosynthesized from the universal C5 precursors as dimethylallyl diphosphate (DMAPP) and isopentenyl diphosphate (IPP) [16,17,18]. These C5 precursors are sequentially catalyzed by prenyltransferases to form prenyl diphosphates, including geranyl pyrophosphate (GPP), farnesyl pyrophosphate (FPP), and geranylgeranyl pyrophosphate (GGPP). These prenyl diphosphates are further used as the immediate precursors for the biosynthesis of various terpenoids, such as monoterpenes (C10), sesquiterpenes (C15), diterpenes (C20), triterpenes (C30), tetraterpenes (C40), as well as polyterpenes (more than C40) by the action of terpene synthase (TPS). Transcription factors (TFs) regulate the accumulation of terpenoids by activating or repressing the promoters of TPS genes to control their expression and five TFs have been identified as being involved in the regulation of terpene synthesis [19,20].
Despite the scientific and industrial interest in sacred lotus, terpenoid components have not been identified. Moreover, the genes involved in terpenoid biosynthesis were been completely studied in aquatic botany. In recent years, RNA sequencing (RNA-Seq) and metabolite analysis were developed as powerful tools for identifying genes involved in the biosynthesis of various secondary metabolites in higher plants, including Salvia officinalis Linn. [21], Cinnamomum camphora (L.) Presl. [10], and Salvia guaranitica St. Hil. [22].
As one of the best-known medicinal plants, few terpenoid compounds have been identified in plant parts of N. nucifera, and the molecular genetic basis of terpenoid biosynthesis pathways is still unveiled. Considering the specific characteristics of sacred lotus, it was of interest to profile the terpenoid biosynthesis and key genes related with the terpenoid metabolic pathway. In this study, the contents and compositions of terpenoids in seeds, young leaves, mature leaves, and flowers with different colors were analyzed with LC-MS/MS. In addition, the key genes and transcription factors involved in terpenoid biosynthesis were screened by RNA-seq. We believe these results will shed light on the mechanism of terpenoid biosynthesis in sacred lotus.

2. Results

2.1. Identification of Terpenoid Components in Various Plant Parts

For unveiling the distribution of terpenoids in flowers with different colors, the widely planted cultivars with a red flower (Jinlinghuodu), white flower (Baiyinlian), and yellow flower (Jinsenianhua) were selected. In order to compare terpenoids in various plant parts, the seeds, young leaves, and mature leaves from Taikonglian 36 were harvested (Figure 1A). In total, 909 metabolites were detected in six sets of samples (three biological replicates per set) of lotus based on the UPLC-MS/MS detection platform and the self-constructed database. Because a broad-target metabolomic analysis was performed, which had some limitations for detecting each class of triterpenoids, cluster analysis of metabolites revealed that the 909 metabolites could be classified into 12 classes (Figure 1B). In addition, a total of 21 terpenoids were identified in different plant parts, and among these 21 terpenoids, 16 terpenoids could be detected in all plant parts, but their contents appeared different in these plant parts (Table 1).

2.2. RNA-Seq and Transcriptomic Assembly

To identify genes involved in terpenoid biosynthesis in sacred lotus, 18 RNA libraries were prepared and analyzed. In total, 249,600,000 reads were obtained from these 18 libraries (accession number: PRJNA857167). After filtering the original data, clean reads were obtained for subsequent analysis (Table S1). Additionally, clean reads were mapped to the sacred lotus genome; the percentage of sequencing reads generated by each sample successfully aligned to the genome was higher than 80%. These results indicated that the quality of these mapped genes was sufficient to conduct the subsequent analysis (Table S2).

2.3. Gene Annotation and Functional Classification

BLAST alignment was utilized to annotate the 29,630 unigenes of sacred lotus with an E-value threshold of 1e−5 in the public databases: NR (NCBI non-redundant protein sequences), KOG (euKaryotic Ortholog Groups), PFAM (Protein family), GO (Gene Ontology), and KEGG (Kyoto Encyclopedia of Genes and Genomes). In summary, 76,070 unigenes were successfully annotated by at least one database and 18,045 unigenes shared annotation in all databases (Figure 2A).
Based on sequence homology, 22,171 annotated unigenes were categorized into three ontologies with 58 GO terms (Figure 2B). Within the category of cellular component (CC), genes matched to 18 GO terms, the most highly represented of which were ‘cell’ (15,881 unigenes), ‘cell part’ (15,817 unigenes), and ‘cellular process’ (13,402 unigenes). For the molecular function (MF) category, ‘binding’ (12,464 unigenes) and ‘catalytic activity’ (10,985 unigenes) were the two most abundant of 12 GO terms. The largest associated term within the 28 GO terms of the biological process (BP) was ‘cellular process’ (13,402 unigenes).
In total, 25,512 unigenes were matched in 134 metabolic pathways and the biosynthesis of secondary metabolites (3964 unigenes), metabolic pathways (2194 unigenes), and plant–pathogen interaction (534) were the three richest pathways. In addition, 122 unigenes were mapped to the ‘Metabolism of terpenoids’, including ‘Terpenoid backbone biosynthesis’ (ko00900, 64 unigenes), ‘Sesquiterpenoid and triterpenoid biosynthesis’ (ko00909, 25 unigenes), and ‘Diterpenoid biosynthesis’ (ko00904, 33 unigenes) (Figure S1).

2.4. Identification of Differential Expression Genes (DEGs)

To fully explore potential DEGs in various plant parts, 22,336 DEGs were identified by comparing the eleven groups (S vs. ML, S vs. RF, S vs. WF, S vs. YF, S vs. YL, WF vs. RF, WF vs. YF, YL vs. ML, YL vs. WF, and YL vs. RF). The total number of DEGs, up-regulated DEGs, and down-regulated DEGs in each group are presented in Table S3. To further characterize the expression patterns of differential genes, we performed k-means clustering analysis on all DEGs. The results showed that 22,336 DEGs were assigned to 12 clusters with different numbers of DEGs per cluster, ranging from 1012 to 4266 (Figure 3).

2.5. Functional Classification of DEGs

The 50 GO-Terms with the lowest q-value in the enrichment analysis results were selected and a column chart of the enrichment items was drawn (Figure 4A). Among the 33 GO terms of biological process (BP), the largest related term was “cell process”. For the Cellular component (CC) category, the DEGs mainly participated in the chloroplast thylakoid, plastid thylakoid, and other organelles. In the category of molecular function (MF), “structural molecule activity” was the largest related term.
According to the KEGG pathway enrichment analysis, biosynthesis of secondary metabolites (1750) and metabolic pathways (3205) were the two richest pathways. In addition, 99 DEGs were positioned as “terpenoid metabolism”, including “Terpenoid backbone biosynthesis” (ko00900, 55), “sesquiterpene and triterpene biosynthesis” (ko00909, 17), and “diterpene biosynthesis” (ko00904, 27) (Figure 4B).

2.6. DEGs Involved in Terpenoid biosynthesis

To unveil the regulatory mechanism of terpenoid accumulation patterns in different plant parts of sacred lotus, we analyzed the expression profiles of genes involved in terpenoid biosynthesis. In total, 105 DEGs related to terpenoid biosynthesis were identified in sacred lotus and their expression patterns are shown in Figure 5. In total, 55 DEGs were annotated to the terpenoid backbone synthesis pathway (ko00900), of which 13 and 11 DEGs participated in the MVA and MEP pathways, respectively. Furthermore, DEGs encoding key enzymes in the MVA pathways exhibited a higher expression level in YL compared with S and ML. Interestingly, AACT2 and HMGS2 exhibited the highest mRNA level in S, while HMGR3 was the highest expressed in ML (Figure 5 and Table S4). Similar to the MVA pathway, the expression levels of DEGs involved in MEP in leaves was still higher than seeds, and the expression levels of DEGs in ML were significantly higher than YL (Figure 5 and Table S4). Moreover, six unigenes encoding IPTS were identified, including geranyl diphosphate synthase (GPPS), farnesyl diphosphate synthase (FPPS), and geranylgeranyl diphosphate synthase (GGPPS). The expression level of FPPS in YL was higher than that in S and ML and also had higher transcript levels in the three differently colored flowers, which was consistent with the expression pattern of DEGs in the MVA pathway. The expression level of GPPS in S was lower than in other parts, but there was no obvious difference in leaves and flowers. The different expression patterns in the three GGPPS transcripts are worth noting. The GGPPS2 transcript exhibited the highest expression level in mature leaves, with the highest mRNA level of GGPPS3 in WF, RF, and YF. Moreover, the expression level of GGPPS1 was relatively lower than that of GGPPS2 and GGPPS3 (Figure 5 and Table S4).
As key enzymes in the synthesis of terpenoids, 49 DEGs encoding terpenoid synthase (TPS) were found, including 15 genes involved in sesquiterpene and triterpene synthesis (ko00909), 26 genes involved in diterpene synthesis (ko00904), and 8 genes involved in monoterpene synthesis (K07385 and K15095) (Table S5). Correlation network analysis was used to model the relationship of the selected 49 genes encoding TPS with the detected terpenoids (Figure 6). As shown in the figure, 40 TPS genes and 20 terpenoids showed high correlation, for example, (-)-germacrene D synthase (GERD3) was positively correlated with five terpenoids, including Triacetate (lmjp04533), 12,13-Dihydroursolic acid (Zmjp01361), Progesterone (Lmqp010784), 2α,3α,23-trihydroxyolean-12-en-28-oic acid (Smpn009230), and Madasiatic acid (Hmjn003948), with correlation coefficients varying from 0.83 to 0.97. Gibberellin 2 beta-dioxygenase (GA2-5), 1,8-cineole synthase (CINS3), and ent-kaurene oxidase (KO) were positively correlated with the same three terpenoids. 24,30-Dihydroxy-12(13)-enolupinol (pmn001700) was positively correlated with nine TPS genes and negatively correlated with four TPS genes. Beta-amyrin synthase (AMYS1, AMYS2, AMYS3) exhibited a strong positive correlation with 24,30-Dihydroxy-12(13)-enolupinol (pmn001700) and 12,13-Dihydroursolic acid (Zmjp013616).

2.7. Co-Expression Network of TPS and WRKY Transcription Factors

Transcription factors (TFs) usually control the expression of TPS genes by activating or repressing their promoters, thereby regulating the accumulation of terpenoids. In total, 47 DEGs were annotated as WRKY transcription factors, 27 of which were highly correlated with 15 terpenoids (Table S6, Figure S2). A co-expression network of these 27 WRKY transcription factors and the 40 genes encoding TPS was constructed (Figure 7). The results showed that WRKY9 (gene-LOC104588731) was negatively correlated with GERD2, GA2-2, GA2-7, and KAO2-2, which was positively correlated with GERD3, (+)-neomenthol dehydrogenase (NEOD4, NEOD5), AMYS3, and AMYS4. WRKY72 (LOC104596936) was closely related to 15 TPS genes. Among them, WRKY72 was positively correlated with NEOD5, AMYS1, AMYS2, and AMYS3, and the Pearson correlation coefficient exceeded 0.9. By contrast, WRKY72 was negatively correlated with GA2-2, GA2-3, GA2-4, GA2-7, gibberellin 3-beta-dioxygenase (GA3-5), GERD2, and ent-kaurenoic acid hydroxylase (KAO2-2). WRKY17 (LOC104603161) was positively correlated with AMYS1, AMYS3, and AMYS4.

2.8. qRT-PCR Validation of DEGs

To verify the reliability of the RNA-Seq data, the expression of 12 unigenes was detected by qRT-PCR (Figure 8). The expression levels of these genes determined by qRT-PCR were almost consistent with those inferred from the RNA-Seq FPKM data (Figure 8), except for the expression pattern of AACT1 and FPPS1 in white flower.

3. Discussion

The remarkable health and disease alleviation activities of the scared lotus are associated with the high content of bioactive compounds, including polyphenols, flavonoids, phenolic acids, alkaloids, terpenoids, steroids, fatty acids, and glycosides [23]. However, most studies focus on flavonoids, alkaloids, and phenolic acids, while few have been conducted on terpenoids [24,25,26]. In this study, in total, 21 terpenoids were detected (Table 1). Although the genes involved in terpenoid biosynthesis are widely studied in different plants [20,27,28,29,30,31,32,33], genes related to the synthesis of terpenoids in sacred lotus have not been studied so far. The precursors of monoterpenes and diterpenes come from the MEP pathway in plastids, while the precursors of sesquiterpenes and triterpenes come from the MVA pathway in the cytoplasm [34]. Considering the two rate-limiting enzymes, HMGR and DXS played important roles in the overall regulation of the MVA and MEP pathways, respectively [35,36,37]. Among five genes encoding HMGR identified in this study, HMGR1, HMGR2, and HMGR4 exhibited the highest expression in YL, with the lowest expression in ML (Figure 5). In total, three genes encoding DXS were detected; DXS2 and DXS3 had the same expression pattern, with higher expression in leaves than in seeds, and there was no significant difference in expression in the three colors of flowers, which corresponded to the metabolic content results. This result indicated that the distribution of terpenoids in different plant parts was closely related to the expression pattern of genes involved in the MVA and MEP pathways. However, the relative contribution of each pathway to the biosynthesis of various terpenoids was still uncertain.
IPTS are the key enzymes that connect the upstream MVA and MEP pathways with downstream isoprenoid biosynthesis branch points of different structures. In recent years, GPPS, GGPPS, FPPS, and CPPS, as the key enzymes in the biosynthesis of isoprenoid compounds, were identified in higher plants. In addition, most studies showed that the transcription levels of these genes were closely related to the content of corresponding terpenoids [38,39,40,41]. For example, FPPS overexpression increased the production of ganoderic acid [42]. A high number of monoterpenes was produced by overexpression of peppermint (Mentha × piperita) GPPS.SSU in transgenic tobacco plants [43] and the content of abietane diterpenes in Salvia sclarea L. hairy roots was increased by engineering the GGPPS and CPPS [44]. In this research, six DEGs encoding IPTS were identified, including one GPPS, two FPPS, and three GGPPS. We found that the expression of FPPS was higher in young leaves than in seeds and mature leaves, which is consistent with the results detected for the metabolic content of triterpenoids. It has been shown that the GGPPS transcript was mainly concentrated in the above-ground parts, such as leaves, flowers, or fruits, and the expression patterns of different GGPPS homologous genes in the same plant were various [45,46,47,48,49,50,51]; the expression patterns of GGPPS2 and GGPPS3 identified in this study are consistent with this conclusion.
As a key enzyme playing an important role in the synthesis of terpenoids, the TPS genes were identified in A. thaliana, tomato, rice, and tobacco [52,53,54,55,56]. Further, a large number of different TPS and their products were the chief reasons for the variety in terpenes [57,58]. In this study, 40 DEGs encoding TPS exhibited a high correlation with 20 terpenes (Figure 6). In previous studies, multiple TPS in most plants were shown to be multi-product enzymes [58]. Our study indicated that most TPS in lotus also showed a high correlation with multiple products and some terpenoids were also highly correlated with multiple TPS. Therefore, it is speculated that these terpenoids may be regulated by multiple TPS. WRKY transcription factors have been reported to regulate terpenoid production by activating or inhibiting promoter binding of TPS [20,59,60,61]. In our study, our results showed that 26 members of WRKY were closely related to 34 TPS genes (Figure 7). These phenomena suggest that WRKY transcription factors play an important role in the regulation of terpenoid production, but the exact mode of action needs to be further explored and confirmed.

4. Materials and Methods

4.1. Plants Materials

The widely planted cultivars “Taikonglian 36”, “Jinsenianhua”, “Jinlinghuodu”, and “Baiyinlian” were planted in pools at Henan University of Technology, with a photoperiod of light for 16 h and dark for 8 h. Sacred loti were cultivated to the full-bloom stage. The materials including seeds (Taikonglian 36), young leaves (Taikonglian 36), mature leaves (Taikonglian 36), and blooming flowers with red (Jinlinghuodu), white (Baiyinlian), and yellow (Jinsenianhua) were freeze dried in liquid nitrogen and kept at −80 °C.

4.2. LC-MS/MS Analysis

Biological samples were freeze dried using a vacuum freeze dryer (Scientz-100F) and crushed using a mixer mill (MM 400, Retsch) with a zirconia bead for 1.5 min at 30 Hz. Then 100 mg of lyophilized powder was dissolved with 1.2 mL 70% methanol solution, and vortexed 30 s every 30 min 6 times in total. The samples were refrigerated 4 °C overnight and centrifugated at 12,000 rpm for 10 min. The extracts were filtrated (SCAA-104, 0.22 μm pore size; ANPEL, Shanghai, China, http://www.anpel.com.cn/, accessed on 5 September 2021) before UPLC-MS/MS analysis. Chromatography and mass spectrometry were performed with the assistance of Wuhan Metware Biotechnology Co., Ltd. (Wuhan, China) as the method described [62]. Metabolite quantification was accomplished by multiple reaction monitoring (MRM) analysis using triple-quadrupole mass spectrometry. The characteristic ions of each substance were screened by triple quadruple rods, the signal intensity (CPS) of the characteristic ions was obtained in the detector. The sample offline mass spectrometry file was opened with MultiaQuant software. The integration and calibration of the peaks were performed, and the peak area (Area) of each peak represented the relative content of the corresponding substance. The internal standard used was 2-chloro-L-phenylalanine (CAS: 103616-89-3, Bailiwick, 106151-100 mg). The HCA (hierarchical cluster analysis) results of metabolites were presented as heatmaps and normalized signal intensities of metabolites (unit variance scaling) were visualized as a color spectrum.

4.3. Total RNA Extraction and Transcriptome Sequencing

Total RNA was extracted using the Total Plant RNA Extraction Kit (Tsingke, Beijing, China). RNA concentration and purity were determined using NanoPhotometer® spectrophotometer (IMPLEN, München, BY, GER) and RNA concentration was measured using Qubit® RNA Assay Kit in Qubit®2.0 Flurometer (Life Technologies, Waltham, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, Palo Alto, CA, USA). A total amount of 1 µg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, Ipswich, MA, USA) following manufacturer’s recommendations. The cDNA libraries were sequenced on the Illumina sequencing platform by Metware Biotechnology Co., Ltd. (Wuhan, China).

4.4. RNA-Seq Analysis

The original data were filtered to remove reads with adapters. HISAT version 2.1.0 (Daehwan Kim, Baltimore, MA, USA) was used to construct the index and compare clean reads to the reference genome (https://www.ncbi.nlm.nih.gov/genome/?term=AQOG01, accessed on 25 September 2021) [63], and StringTie version 1.3.4d (Mihaela Pertea, Baltimore, MA, USA) was used for new gene prediction [64]. We used featureCounts version 1.6.2 (Yang Liao, Parkville, VIC, AUS) to calculate the gene alignment and then calculated the transcript per million fragments mapped (FPKM) of each gene based on the gene length. DESeq2 version 1.22.1 (Michael I Love, Heidelberg, BW, GER) [65,66] was used to analyze the original count data and to screen for DEGs, and genes satisfying |log2Fold Change| ≥ 1 and False Discovery Rate (FDR) < 0.05 were defined as differentially expressed genes (DEGs). The enrichment analysis was performed based on the hypergeometric test. The gene co-expression network analysis was performed using the Metware Cloud, a free online platform for data analysis (https://cloud.metware.cn, accessed on 3 January 2022), and Pearson correlation coefficient greater than |0.8|.

4.5. Validation of DEGs by qRT-PCR Analysis

The key DEGs in the terpenoid synthesis pathway were selected to validate their expression levels by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Synthesis of cDNA for qRT-PCR was performed with First Strand cDNA Synthesis Kit (ThermoFisher Scientific, Shanghai, China). The primers were synthesized by Tsingke Biotechnology Co., Ltd. (Wuhan, China) (Table S7). The reaction system of qRT-PCR was as follows: cDNA template 0.5 μL at 3000 ng/μL, gene-specific upstream and downstream primers 0.4 μL each at 10 μM, 2 × ChamQ Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China) 10 μL, 8.7 μL of double-distilled water (ddH2O). qRT-PCR was performed by CFX96 Real Time System (Bio-Rad, Hercules, CA, USA). Young leaves were selected as the control and 26s rRNA as the internal reference genes. The relative expression of genes in the qRT-PCR experiment was analyzed by the 2−△△CT method [67]. At least three biological replicates were set for each sample.

5. Conclusions

In this study, 21 terpenoids were detected in sacred lotus and 105 genes involved in terpenoid synthesis were identified. WRKY, an important transcription factor, is highly correlated with the TPS gene family and might play an important regulatory role in the terpenoid synthesis process. This study illustrated the metabolic diversity of terpenoids in sacred lotus and provided a global overview of the gene expression profiles related with terpenoid biosynthesis in sacred lotus.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27144599/s1, Figure S1: KEGG annotation of unigenes. The unigenes were divided into five branches according to the KEGG metabolic pathway: Cellular Processes, Environmental Information Processing, Genetic Information Processing, Metabolism, and Organismal Systems, Figure S2: Correlation analysis of WRKY TFs and released amounts of terpene compounds (Pearson correlation coefficient > 0.8). The circular nodes are WRKY TFs, the triangular nodes are terpenes, the red line represents positive correlation, and the blue line represents negative correlation, Table S1: Sequencing output statistics, Table S2: Comparison efficiency statistics, Table S3: DEGS of groups, Table S4: The fpkm value of genes involved in terpene backbone synthesis, Table S5: The fpkm value of genes involved in terpenoid synthesis, Table S6: The fpkm value of genes annotated as WRKY transcription factors, Table S7: Primers used in qPCR.

Author Contributions

Conceptualization, L.Q. and F.D.; methodology, N.Y.; software, C.Z.; validation, X.Z., J.T. and Z.W.; formal analysis, L.Q.; investigation, C.D.; resources, L.Y.; data curation, C.D.; writing—original draft preparation, L.Q.; writing—review and editing, C.D.; visualization, L.Y.; supervision, C.D.; project administration, C.D.; funding acquisition, C.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the China Postdoctoral Science Foundation (grant number 2020M672308), National Natural Science Foundation of China (grant number no. 31601366), Natural Science Foundation of Henan Province (grant number 202300410114), Cultivation Programme for Young Backbone Teachers in Henan University of Technology, the Innovative Funds Plan of Henan University of Technology (grant number 2020ZKCJ23), and Program for Innovative Research Team (in Science and Technology) in University of Henan Province (grant number 19IRTSTHN008).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw sequencing data were deposited in the NCBI Bioproject database (accession number: PRJNA857167, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA857167, accessed on 17 July 2022).

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples are not available.

References

  1. Ferrer-Gallego, P.P.; Boisset, F.; Jarvis, C. Typification of the sacred lotus nelumbo nucifera (nelumbonaceae). Taxon 2015, 64, 156–159. [Google Scholar] [CrossRef]
  2. Lin, Z.; Damaris, R.N.; Shi, T.; Li, J.; Yang, P. Transcriptomic analysis identifies the key genes involved in stamen petaloid in lotus (Nelumbo nucifera). BMC Genom. 2018, 19, 554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Li, J.; Shi, T.; Huang, L.; He, D.l.; Nyong, A.T.M.; Yang, P. MaSystematic transcriptomic analysis provides insights into lotus (Nelumbo nucifera) seed development. Plant Growth Regul. 2018, 86, 339–350. [Google Scholar] [CrossRef]
  4. Huang, B.; Ban, X.; He, J.; Jing, T.; Tian, J.; Wang, Y. Hepatoprotective and antioxidant activity of ethanolic extracts of edible lotus (Nelumbo nucifera Gaertn.) leaves. Food Chem. 2010, 120, 873–878. [Google Scholar] [CrossRef]
  5. Huang, B.; He, J.; Ban, X.; Zeng, H.; Yao, X.; Wang, Y. Antioxidant activity of bovine and porcine meat treated with extracts from edible lotus (Nelumbo nucifera) rhizome knot and leaf. Meat Sci. 2011, 87, 46–53. [Google Scholar] [CrossRef] [PubMed]
  6. Chen, S.; Wu, B.H.; Fang, J.B.; Liu, Y.L.; Zhang, H.H.; Fang, L.C.; Guan, L.; Li, S.H. Analysis of flavonoids from lotus (Nelumbo nucifera) leaves using high performance liquid chromatography/photodiode array detector tandem electrospray ionization mass spectrometry and an extraction method optimized by orthogonal design. J. Chromatogr. A 2012, 1227, 145–153. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, G.; Zhu, M.; Guo, M. Research advances in traditional and modern use of Nelumbo nucifera: Phytochemicals, health promoting activities and beyond. Crit. Rev. Food Sci. Nutr. 2019, 59, S189–S209. [Google Scholar] [CrossRef]
  8. Mukherjee, P.K.; Mukherjee, D.; Maji, A.K.; Rai, S.; Heinrich, M. The sacred lotus (Nelumbo nucifera)-phytochemical and therapeutic profile. J. Pharm. Pharmacol. 2009, 61, 407–422. [Google Scholar] [CrossRef]
  9. Neill, S.O.; Gould, K.S.; Kilmartin, P.A.; Mitchell, K.A.; Markham, K.R. Antioxidant capacities of green and cyanic leaves in the sun species, Quintinia serrata. Funct. Plant Biol. 2002, 29, 1437–1443. [Google Scholar] [CrossRef]
  10. Chen, C.; Zheng, Y.; Zhong, Y.; Wu, Y.; Li, Z.; Xu, L.A.; Xu, M. Transcriptome analysis and identification of genes related to terpenoid biosynthesis in Cinnamomum camphora. BMC Genom. 2018, 19, 550. [Google Scholar] [CrossRef] [Green Version]
  11. Huang, L.; Li, J.; Ye, H.; Li, C.; Wang, H.; Liu, B.; Zhang, Y. Molecular characterization of the pentacyclic triterpenoid biosynthetic pathway in catharanthus roseus. Planta 2012, 236, 1571–1581. [Google Scholar] [CrossRef] [PubMed]
  12. Zwawiak, J.; Pawelczyk, A.; Olender, D.; Zaprutko, L. Structure and activity of pentacyclic triterpenes codrugs. A review. Mini Rev. Med. Chem. 2021, 21, 1509–1526. [Google Scholar] [CrossRef] [PubMed]
  13. Miranda, R.S.; de Jesus, B.; da Silva Luiz, S.; Viana, C.B.; Adao, M.C.; Figueiredo, F.S.; Carvalho, T.; Silva, M.L.; Londero, V.S.; de Souza Figueiredo, C.T.; et al. Antiinflammatory activity of natural triterpenes—An overview from 2006 to 2021. Phytother. Res. 2022, 36, 1459–1506. [Google Scholar] [CrossRef] [PubMed]
  14. Yao, L.; Lu, J.; Wang, J.; Gao, W.Y. Advances in biosynthesis of triterpenoid saponins in medicinal plants. Chin. J. Nat. Med. 2020, 18, 417–424. [Google Scholar] [CrossRef]
  15. Szakiel, A.; Pączkowski, C.; Henry, M. Influence of environment-al abiotic factors on the content of saponins in plants. Phyto-chem. Rev. 2010, 10, 471–491. [Google Scholar] [CrossRef]
  16. Bouvier, F.; Rahier, A.; Camara, B. Biogenesis, molecular regulation and function of plant isoprenoids. Prog. Lipid Res. 2005, 44, 357–429. [Google Scholar] [CrossRef]
  17. Cazzonelli, C.I.; Pogson, B.J. Source to sink: Regulation of carotenoid biosynthesis in plants. Trends Plant Sci. 2010, 15, 266–274. [Google Scholar] [CrossRef]
  18. Pulido, P.; Perello, C.; Rodriguez-Concepcion, M. New insights into plant isoprenoid metabolism. Mol. Plant 2012, 5, 964–967. [Google Scholar] [CrossRef] [Green Version]
  19. Gantet, P.; Memelink, J. Transcription factors: Tools to engineer the production of pharmacologically active plant metabolites. Trends Pharmacol. Sci. 2002, 23, 563–569. [Google Scholar] [CrossRef]
  20. He, X.Y.; Wang, H.; Yang, J.F.; Deng, K.; Wang, T. RNA sequencing on Amomum villosum Lour. induced by MeJA identifies the genes of WRKY and terpene synthases involved in terpene biosynthesis. Genome 2018, 61, 91–102. [Google Scholar] [CrossRef] [Green Version]
  21. Ali, M.; Li, P.; She, G.; Chen, D.; Wan, X.; Zhao, J. Transcriptome and metabolite analyses reveal the complex metabolic genes involved in volatile terpenoid biosynthesis in garden sage (Salvia officinalis). Sci. Rep. 2017, 7, 16074. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Ali, M.; Hussain, R.M.; Rehman, N.U.; She, G.; Li, P.; Wan, X.; Guo, L.; Zhao, J. De novo transcriptome sequencing and metabolite profiling analyses reveal the complex metabolic genes involved in the terpenoid biosynthesis in Blue Anise Sage (Salvia guaranitica L.). DNA Res. 2018, 25, 597–617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bishayee, A.; Patel, P.A.; Sharma, P.; Thoutireddy, S.; Das, N. Lotus (Nelumbo nucifera Gaertn.) and Its Bioactive Phytocompounds: A Tribute to Cancer Prevention and Intervention. Cancers 2022, 14, 529. [Google Scholar] [CrossRef] [PubMed]
  24. Limwachiranon, J.; Huang, H.; Shi, Z.; Li, L.; Luo, Z. Lotus Flavonoids and Phenolic Acids: Health Promotion and Safe Consumption Dosages. Compr. Rev. Food Sci. Food Saf. 2018, 17, 458–471. [Google Scholar] [CrossRef] [Green Version]
  25. Menéndez-Perdomo, I.M.; Facchini, P.J. Benzylisoquinoline Alkaloids Biosynthesis in Sacred Lotus. Molecules 2018, 23, 2899. [Google Scholar] [CrossRef] [Green Version]
  26. Tungmunnithum, D.; Pinthong, D.; Hano, C. Flavonoids from Nelumbo nucifera Gaertn., a Medicinal Plant: Uses in Traditional Medicine, Phytochemistry and Pharmacological Activities. Medicines 2018, 5, 127. [Google Scholar] [CrossRef] [Green Version]
  27. Damasco, G.; Shivakumar, V.S.; Misciewicz, T.M.; Daly, D.C.; Fine, P. Leaf Transcriptome Assembly of Protium copal (Burseraceae) and Annotation of Terpene Biosynthetic Genes. Genes 2019, 10, 392. [Google Scholar] [CrossRef] [Green Version]
  28. George, J.K.; Shelvy, S.; Fayad, A.M.; Umadevi, P.; Angadi, U.B.; Iquebal, M.A.; Jaiswal, S.; Rai, A.; Kumar, D. De novo transcriptome sequencing assisted identification of terpene synthases from black pepper (Piper nigrum) berry. Physiol. Mol. Biol. Plants 2021, 27, 1153–1161. [Google Scholar] [CrossRef]
  29. Han, X.J.; Wang, Y.D.; Chen, Y.C.; Lin, L.Y.; Wu, Q.K. Transcriptome sequencing and expression analysis of terpenoid biosynthesis genes in Litsea cubeba. PLoS ONE 2013, 8, e76890. [Google Scholar] [CrossRef] [Green Version]
  30. 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. [Google Scholar] [CrossRef]
  31. Qiu, F.; Wang, X.; Zheng, Y.; Wang, H.; Liu, X.; Su, X. Full-Length Transcriptome Sequencing and Different Chemotype Expression Profile Analysis of Genes Related to Monoterpenoid Biosynthesis in Cinnamomum porrectum. Int. J. Mol. Sci. 2019, 20, 6230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Sun, L.; Zhu, B.; Zhang, X.; Zhang, G.; Yan, A.; Wang, H.; Wang, X.; Xu, H. Transcriptome profiles of three Muscat table grape cultivars to dissect the mechanism of terpene biosynthesis. Sci. Data 2019, 6, 89. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Yang, Z.; Xie, C.; Huang, Y.; An, W.; Liu, S.; Huang, S.; Zheng, X. Metabolism and transcriptome profiling provides insight into the genes and transcription factors in volved in monoterpene biosynthesis of borneol chemotype of Cinnamomum camphora induced by mechanical damage. PeerJ 2021, 9, e11465. [Google Scholar] [CrossRef]
  34. Tholl, D. Biosynthesis and biological functions of terpenoids in plants. Adv. Biochem. Eng./Biotechnol. 2015, 148, 63–106. [Google Scholar] [CrossRef]
  35. Cordoba, E.; Salmi, M.; León, P. Unravelling the regulatory mechanisms that modulate the MEP pathway in higher plants. J. Exp. Bot. 2009, 60, 2933–2943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Friesen, J.A.; Rodwell, V.W. The 3-hydroxy-3-methylglutaryl coenzyme-A (HMG-CoA) reductases. Genome Biol. 2004, 5, 248. [Google Scholar] [CrossRef] [Green Version]
  37. Liao, P.; Hemmerlin, A.; Bach, T.J.; Chye, M.L. The potential of the mevalonate pathway for enhanced isoprenoid production. Biotechnol. Adv. 2016, 34, 697–713. [Google Scholar] [CrossRef]
  38. Srivastava, Y.; Tripathi, S.; Mishra, B.; Sangwan, N.S. Cloning and homologous characterization of geranylgeranyl pyrophosphate synthase (GGPPS) from Withania somnifera revealed alterations in metabolic flux towards gibberellic acid biosynthesis. Planta 2022, 256, 4. [Google Scholar] [CrossRef]
  39. Rai, A.; Smita, S.S.; Singh, A.K.; Shanker, K.; Nagegowda, D.A. Heteromeric and homomeric geranyl diphosphate synthases from Catharanthus roseus and their role in monoterpene indole alkaloid biosynthesis. Mol. Plant 2013, 6, 1531–1549. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, G.; Dixon, R.A. Heterodimeric geranyl(geranyl)diphosphate synthase from hop (Humulus lupulus) and the evolution of monoterpene biosynthesis. Proc. Natl. Acad. Sci. USA 2009, 106, 9914–9919. [Google Scholar] [CrossRef] [Green Version]
  41. Zhu, X.F.; Suzuki, K.; Saito, T.; Okada, K.; Tanaka, K.; Nakagawa, T.; Matsuda, H.; Kawamukai, M. Geranylgeranyl pyrophosphate synthase encoded by the newly isolated gene GGPS6 from Arabidopsis thaliana is localized in mitochondria. Plant Mol. Biol. 1997, 35, 331–341. [Google Scholar] [CrossRef] [PubMed]
  42. Fei, Y.; Li, N.; Zhang, D.H.; Xu, J.W. Increased production of ganoderic acids by overexpression of homologous farnesyl diphosphate synthase and kinetic modeling of ganoderic acid production in Ganoderma lucidum. Microb. Cell Factories 2019, 18, 115. [Google Scholar] [CrossRef] [PubMed]
  43. Yin, J.L.; Wong, W.S.; Jang, I.C.; Chua, N.H. Co-expression of peppermint geranyl diphosphate synthase small subunit enhances monoterpene production in transgenic tobacco plants. New Phytol. 2017, 213, 1133–1144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Vaccaro, M.C.; Alfieri, M.; De Tommasi, N.; Moses, T.; Goossens, A.; Leone, A. Boosting the Synthesis of Pharmaceutically Active Abietane Diterpenes in S. sclarea Hairy Roots by Engineering the GGPPS and CPPS Genes. Front. Plant Sci. 2020, 11, 924. [Google Scholar] [CrossRef]
  45. Engprasert, S.; Taura, F.; Kawamukai, M.; Shoyama, Y. Molecular cloning and functional expression of geranylgeranyl pyrophosphate synthase from Coleus forskohlii Briq. BMC Plant Biol. 2004, 4, 18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Kai, G.; Liao, P.; Zhang, T.; Zhou, W.; Wang, J.; Xu, H.; Liu, Y.; Lin, Z. Characterization, expression profiling, and functional identification of a gene encoding geranylgeranyl diphosphate synthase from Salvia miltiorrhiza. Biotechnol. Bioprocess Eng. 2010, 15, 236–245. [Google Scholar] [CrossRef]
  47. Rahimi, S.; Kim, Y.J.; Devi, B.S.R.; Khorolragchaa, A.; Sukweenadhi, J.; Yang, D.C. Isolation and characterization of Panax ginseng geranylgeranyl-diphosphate synthase genes responding to drought stress. Eur. J. Plant Pathol. 2015, 142, 747–758. [Google Scholar] [CrossRef] [Green Version]
  48. Takaya, A.; Zhang, Y.W.; Asawatreratanakul, K.; Wititsuwannakul, D.; Wititsuwannakul, R.; Takahashi, S.; Koyama, T. Cloning, expression and characterization of a functional cDNA clone encoding geranylgeranyl diphosphate synthase of Hevea brasiliensis. Biochim. Biophys. Acta 2003, 1625, 214–220. [Google Scholar] [CrossRef]
  49. Wang, Y.; Miao, Z.; Tang, K. Molecular cloning and functional expression analysis of a new gene encoding geranylgeranyl diphosphate synthase from hazel (Corylus avellana L. Gasaway). Mol. Biol. Rep. 2010, 37, 3439–3444. [Google Scholar] [CrossRef]
  50. Beck, G.; Coman, D.; Herren, E.; Ruiz-Sola, M.A.; Rodríguez-Concepción, M.; Gruissem, W.; Vranová, E. Characterization of the GGPP synthase gene family in Arabidopsis thaliana. Plant Mol. Biol. 2013, 82, 393–416. [Google Scholar] [CrossRef] [Green Version]
  51. Okada, K.; Saito, T.; Nakagawa, T.; Kawamukai, M.; Kamiya, Y. Five geranylgeranyl diphosphate synthases expressed in different organs are localized into three subcellular compartments in Arabidopsis. Plant Physiol. 2000, 122, 1045–1056. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Aubourg, S.; Lecharny, A.; Bohlmann, J. Genomic analysis of the terpenoid synthase (AtTPS) gene family of Arabidopsis thaliana. Mol. Genet. Genom. 2002, 267, 730–745. [Google Scholar] [CrossRef] [PubMed]
  53. Lee, S.; Chappell, J. Biochemical and genomic characterization of terpene synthases in Magnolia grandiflora. Plant Physiol. 2008, 147, 1017–1033. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Yang, H.L.; Liu, Y.J.; Wang, C.L.; Zeng, Q.Y. Molecular evolution of trehalose-6-phosphate synthase (TPS) gene family in Populus, Arabidopsis and rice. PLoS ONE 2012, 7, e42438. [Google Scholar] [CrossRef]
  55. Zang, B.; Li, H.; Li, W.; Deng, X.W.; Wang, X. Analysis of trehalose-6-phosphate synthase (TPS) gene family suggests the formation of TPS complexes in rice. Plant Mol. Biol. 2011, 76, 507–522. [Google Scholar] [CrossRef] [PubMed]
  56. Zhou, F.; Pichersky, E. The complete functional characterisation of the terpene synthase family in tomato. New Phytol. 2020, 226, 1341–1360. [Google Scholar] [CrossRef] [Green Version]
  57. Karunanithi, P.S.; Zerbe, P. Terpene Synthases as Metabolic Gatekeepers in the Evolution of Plant Terpenoid Chemical Diversity. Front. Plant Sci. 2019, 10, 1166. [Google Scholar] [CrossRef] [Green Version]
  58. 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. Plant J. 2011, 66, 212–229. [Google Scholar] [CrossRef]
  59. Suttipanta, N.; Pattanaik, S.; Kulshrestha, M.; Patra, B.; Singh, S.K.; Yuan, L. The transcription factor CrWRKY1 positively regulates the terpenoid indole alkaloid biosynthesis in Catharanthus roseus. Plant Physiol. 2011, 57, 2081–2093. [Google Scholar] [CrossRef] [Green Version]
  60. Wang, M.; Qiu, X.; Pan, X.; Li, C. Transcriptional Factor-Mediated Regulation of Active Component Biosynthesis in Medicinal Plants. Curr. Pharm. Biotechnol. 2021, 22, 848–866. [Google Scholar] [CrossRef]
  61. Patra, B.; Schluttenhofer, C.; Wu, Y.; Pattanaik, S.; Yuan, L. Transcriptional regulation of secondary metabolite biosynthesis in plants. Biochim. Biophys. Acta 2013, 1829, 1236–1247. [Google Scholar] [CrossRef] [PubMed]
  62. Su, Y.; Zhang, J.; Xu, Z.; Li, J.; Wang, J. Integrative analysis of metabolome and transcriptome reveals the mechanism of color formation in white root (Salvia miltiorrhiza). Ind. Crops Prod. 2021, 170, 113784. [Google Scholar] [CrossRef]
  63. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [Green Version]
  65. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Varet, H.; Brillet-Guéguen, L.; Coppée, J.Y.; Dillies, M.A. SARTools: A DESeq2- and EdgeR-Based R Pipeline for Comprehensive Differential Analysis of RNA-Seq Data. PLoS ONE 2016, 11, e0157022. [Google Scholar] [CrossRef] [Green Version]
  67. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Molecules 27 04599 g001 550
Figure 1. Metabolic analysis of various plant parts in sacred lotus. (A) Sample morphology of various plant parts. (B) Sample cluster map. The horizontal is the sample name, the vertical is the metabolite information, and the different colors are the values obtained after the relative content standardization process.
Figure 1. Metabolic analysis of various plant parts in sacred lotus. (A) Sample morphology of various plant parts. (B) Sample cluster map. The horizontal is the sample name, the vertical is the metabolite information, and the different colors are the values obtained after the relative content standardization process.
Molecules 27 04599 g001
Molecules 27 04599 g002 550
Figure 2. Annotation information of assembled unigenes in sacred lotus. (A) Venn diagram of the distribution of annotation information from different public databases. (B) GO classification of unigenes. The abscissa represents the secondary GO entry and the ordinate represents the number of genes in the GO entry. The secondary GO entry includes three parts: molecular function, biological process, and cell composition.
Figure 2. Annotation information of assembled unigenes in sacred lotus. (A) Venn diagram of the distribution of annotation information from different public databases. (B) GO classification of unigenes. The abscissa represents the secondary GO entry and the ordinate represents the number of genes in the GO entry. The secondary GO entry includes three parts: molecular function, biological process, and cell composition.
Molecules 27 04599 g002
Molecules 27 04599 g003 550
Figure 3. Clusters of DEGs. Twelve clusters were classified based on gene expression pattern: Sub Class1–Sub Class12. The horizontal coordinate represents the sample and the vertical coordinate represents the centralized and normalized expressions.
Figure 3. Clusters of DEGs. Twelve clusters were classified based on gene expression pattern: Sub Class1–Sub Class12. The horizontal coordinate represents the sample and the vertical coordinate represents the centralized and normalized expressions.
Molecules 27 04599 g003
Molecules 27 04599 g004 550
Figure 4. Annotation information of DEGs. (A) GO enrichment histogram of DEGs. The abscissa indicates the ratio of genes annotated to the entry to the total number of genes and the ordinate indicates the name of the GO entry. All DEGs were divided into three GO entries: biological process, cellular component, and molecular function. (B) KEGG classification column chart. DEGs were divided into five branches according to the KEGG metabolic pathway: Cellular Processes, Environmental Information Processing, Genetic Information Processing, Metabolism, and Organismal Systems.
Figure 4. Annotation information of DEGs. (A) GO enrichment histogram of DEGs. The abscissa indicates the ratio of genes annotated to the entry to the total number of genes and the ordinate indicates the name of the GO entry. All DEGs were divided into three GO entries: biological process, cellular component, and molecular function. (B) KEGG classification column chart. DEGs were divided into five branches according to the KEGG metabolic pathway: Cellular Processes, Environmental Information Processing, Genetic Information Processing, Metabolism, and Organismal Systems.
Molecules 27 04599 g004
Molecules 27 04599 g005 550
Figure 5. DEGs in terpenoid biosynthesis of sacred lotus. Typical terpenoid biosynthetic pathway with homology heat map of gene transcription level in transcriptome data with substrates and products. The arrows connect the substrate to their corresponding product. The expression pattern of each unigene is displayed in a six-column grid; from left to right are seed, young leaf, mature leaf, red flower, white flower, and yellow flower. Transcript-level data is represented by FPKM. AACT, Acetyl-CoA C-acetyltransferase; HMGS, hydroxymethylglutaryl-CoA synthase; HMGR, hydroxymethylglutaryl-CoA reductase (NADPH); MVK, mevalonate kinase; PMK, phosphomevalonate kinase; MVD, diphosphomevalonate decarboxylase; DXS, 1-deoxy-D-xylulose-5-phosphate synthase; DXR, 1-deoxy-D-xylulose-5-phosphate reductoisomerase; MCT, 2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase; CMK, 4-diphosphocytidyl-2-C-methyl-D-erythritol kinase; MDS, 2-C-methyl-D- erythritol 2,4-cyclodiphosphate synthase; HDS, (E)-4-hydroxy-3-methylbut-2- enyl-diphosphatesynthase; HDR, 4-hydroxy-3-methylbut-2-en-1-yl diphosphate reductase; IDI, isopentenyl-diphosphate Delta-isomerase; DMAPP, dimethylallyl diphosphate; GPPS, geranyl diphosphate synthase; FPPS, farnesyl diphosphate synthase; GGPPS, geranylgeranyl diphosphate synthase; NRS, (3S,6E)-nerolidol synthase; GERD, (-)-germacrene D synthase; FAR-DT, farnesyl-diphosphate farnesyltransferase; SQUS, squalene monooxygenase; AMYS, beta-amyrin synthase; CINS, 1,8-cineole synthase; NEOD, (+)-neomenthol dehydrogenase; CPS, ent-copalyl diphosphate synthase; KS, ent-kaurene synthase; KO, ent-kaurene oxidase; KAO2, ent-kaurenoic acid oxidase 2; GA3, gibberellin 3-beta-dioxygenase; GA2, gibberellin 2-oxidase; GA20, gibberellin 20-oxidase; MAS, momilactone-A synthase.
Figure 5. DEGs in terpenoid biosynthesis of sacred lotus. Typical terpenoid biosynthetic pathway with homology heat map of gene transcription level in transcriptome data with substrates and products. The arrows connect the substrate to their corresponding product. The expression pattern of each unigene is displayed in a six-column grid; from left to right are seed, young leaf, mature leaf, red flower, white flower, and yellow flower. Transcript-level data is represented by FPKM. AACT, Acetyl-CoA C-acetyltransferase; HMGS, hydroxymethylglutaryl-CoA synthase; HMGR, hydroxymethylglutaryl-CoA reductase (NADPH); MVK, mevalonate kinase; PMK, phosphomevalonate kinase; MVD, diphosphomevalonate decarboxylase; DXS, 1-deoxy-D-xylulose-5-phosphate synthase; DXR, 1-deoxy-D-xylulose-5-phosphate reductoisomerase; MCT, 2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase; CMK, 4-diphosphocytidyl-2-C-methyl-D-erythritol kinase; MDS, 2-C-methyl-D- erythritol 2,4-cyclodiphosphate synthase; HDS, (E)-4-hydroxy-3-methylbut-2- enyl-diphosphatesynthase; HDR, 4-hydroxy-3-methylbut-2-en-1-yl diphosphate reductase; IDI, isopentenyl-diphosphate Delta-isomerase; DMAPP, dimethylallyl diphosphate; GPPS, geranyl diphosphate synthase; FPPS, farnesyl diphosphate synthase; GGPPS, geranylgeranyl diphosphate synthase; NRS, (3S,6E)-nerolidol synthase; GERD, (-)-germacrene D synthase; FAR-DT, farnesyl-diphosphate farnesyltransferase; SQUS, squalene monooxygenase; AMYS, beta-amyrin synthase; CINS, 1,8-cineole synthase; NEOD, (+)-neomenthol dehydrogenase; CPS, ent-copalyl diphosphate synthase; KS, ent-kaurene synthase; KO, ent-kaurene oxidase; KAO2, ent-kaurenoic acid oxidase 2; GA3, gibberellin 3-beta-dioxygenase; GA2, gibberellin 2-oxidase; GA20, gibberellin 20-oxidase; MAS, momilactone-A synthase.
Molecules 27 04599 g005
Molecules 27 04599 g006 550
Figure 6. Correlation analysis of TPS genes and released amounts of terpene compounds (Pearson correlation coefficient > 0.8). The circular nodes are TPS genes, the triangular nodes are terpenes, the red line represents positive correlation, and the blue line represents negative correlation.
Figure 6. Correlation analysis of TPS genes and released amounts of terpene compounds (Pearson correlation coefficient > 0.8). The circular nodes are TPS genes, the triangular nodes are terpenes, the red line represents positive correlation, and the blue line represents negative correlation.
Molecules 27 04599 g006
Molecules 27 04599 g007 550
Figure 7. Correlation analysis of TPS genes and WRKY transcription factors (Pearson correlation coefficient > 0.8). The circular nodes are TPS genes, the triangular nodes are WRKY, the red line represents positive correlation, and the blue line represents negative correlation.
Figure 7. Correlation analysis of TPS genes and WRKY transcription factors (Pearson correlation coefficient > 0.8). The circular nodes are TPS genes, the triangular nodes are WRKY, the red line represents positive correlation, and the blue line represents negative correlation.
Molecules 27 04599 g007
Molecules 27 04599 g008 550
Figure 8. qRT-PCR and RNA-Seq analysis of 12 terpenoid-pathway-related candidate unigenes in sacred lotus. FPKM terms of various plants parts were determined by RNA-seq. Relative expression of genes was calculated using young leaves as control. Standard error of the mean for three biological replicates is represented by the error bars. AACT, Acetyl-CoA C-acetyltransferase; HMGS, hydroxymethylglutaryl-CoA synthase; HMGR, hydroxymethylglutaryl-CoA reductase (NADPH); FPPS, farnesyl diphosphate synthase; MVK, mevalonate kinase; MDS, 2-C-methyl-D- erythritol 2,4-cyclodiphosphate synthase; PMK, phosphomevalonate kinase; DXR, 1-deoxy-D-xylulose-5-phosphate reductoisomerase; IDI, isopentenyl-diphosphate Delta-isomerase; DXS, 1-deoxy-D-xylulose-5-phosphate synthase.
Figure 8. qRT-PCR and RNA-Seq analysis of 12 terpenoid-pathway-related candidate unigenes in sacred lotus. FPKM terms of various plants parts were determined by RNA-seq. Relative expression of genes was calculated using young leaves as control. Standard error of the mean for three biological replicates is represented by the error bars. AACT, Acetyl-CoA C-acetyltransferase; HMGS, hydroxymethylglutaryl-CoA synthase; HMGR, hydroxymethylglutaryl-CoA reductase (NADPH); FPPS, farnesyl diphosphate synthase; MVK, mevalonate kinase; MDS, 2-C-methyl-D- erythritol 2,4-cyclodiphosphate synthase; PMK, phosphomevalonate kinase; DXR, 1-deoxy-D-xylulose-5-phosphate reductoisomerase; IDI, isopentenyl-diphosphate Delta-isomerase; DXS, 1-deoxy-D-xylulose-5-phosphate synthase.
Molecules 27 04599 g008
Table
Table 1. Relative content of terpenoids.
Table 1. Relative content of terpenoids.
IndexCompoundsClass IIQ1 (Da)Q3 (Da)LevelSYLMLWFRFYF
Lmyp004084Perillyl alcoholMonoterpenoids1.53 × 1021.07 × 10233.77 × 1047.46 × 1042.75 × 1051.42 × 1041.28 × 1053.54 × 104
Hmlp008487Blumenol CMonoterpenoids2.11 × 1028.11 × 10125.35 × 1043.68 × 1044.29 × 1044.13 × 1044.54 × 1044.00 × 104
MWSslk208Kaurenoic AcidDitepenoids3.01 × 1023.01 × 10221.39 × 1053.99 × 1043.45 × 1041.92 × 1051.25 × 1052.12 × 104
Lmbn014696Pimaric acidDitepenoids3.01 × 1023.01 × 10221.31 × 1054.26 × 1052.61 × 1056.73 × 1055.65 × 1052.27 × 105
Lmqp010784ProgesteroneDitepenoids3.15 × 1021.09 × 10224.26 × 1041.78 × 1061.10 × 1061.65 × 1052.05 × 1062.96 × 105
Lmdn0015606-DeoxyCatalpolSesquiterpenoids3.45 × 1021.65 × 10227.79 × 1042.12 × 1053.96 × 1059.90 × 1041.53 × 1052.15 × 105
Hmmn0039647-Deoxyloganic acidSesquiterpenoids3.59 × 1021.97 × 10231.17 × 1052.09 × 1044.41 × 1046.32 × 1044.45 × 1047.78 × 103
MWSslk197SecoxyloganinSesquiterpenoids4.05 × 1021.65 × 10221.31 × 1059.96 × 1035.50 × 1043.76 × 1042.46 × 1042.01 × 104
Lmjp004533Kisasagenol A TriacetateDitepenoids4.47 × 1022.85 × 10229.00 × 1006.08 × 1041.13 × 1053.12 × 1041.89 × 1046.32 × 104
MWSmce134Betulonic acidTriterpene4.53 × 1024.53 × 10234.16 × 1037.04 × 1033.15 × 1031.13 × 1039.60 × 1029.00 × 100
pmn00170024,30-Dihydroxy-12(13)-enolupinolTriterpene4.55 × 1024.55 × 10211.84 × 1045.13 × 1041.33 × 1043.60 × 1054.79 × 1053.42 × 105
HJN110Betulinic acidTriterpene4.55 × 1024.55 × 10221.30 × 1054.16 × 1051.77 × 1055.79 × 1041.64 × 1051.84 × 105
mws4053Ursolic acidTriterpene4.55 × 1024.55 × 10235.60 × 1054.24 × 1061.08 × 1065.73 × 1052.41 × 1054.01 × 105
MWSmce0523-Epiursolic acidTriterpene4.55 × 1024.55 × 10235.39 × 1054.05 × 1061.14 × 1065.63 × 1052.45 × 1053.94 × 105
Zmjp01361612,13-Dihydroursolic acidTriterpene4.59 × 1024.23 × 10239.00 × 1003.98 × 1041.62 × 1041.49 × 1041.89 × 1051.15 × 105
pmn0017062-Hydroxyoleanolic acidTriterpene4.71 × 1024.71 × 10216.68 × 1051.94 × 1075.08 × 1061.10 × 1061.00 × 1071.41 × 106
Li512114Corosolic Acid Methyl EsterTriterpene4.85 × 1024.53 × 10233.10 × 1035.62 × 1035.89 × 1031.18 × 1032.29 × 1031.32 × 103
MWSmce394Tormentic acidTriterpene4.87 × 1024.69 × 10233.62 × 1033.18 × 1041.52 × 1046.23 × 1026.18 × 1039.00 × 100
Hmjn003948Madasiatic acidTriterpene4.87 × 1024.87 × 10223.49 × 1041.87 × 1065.72 × 1052.35 × 1052.15 × 1066.11 × 105
pmn001426Euscaphic acidTriterpene4.87 × 1024.69 × 10233.18 × 1033.23 × 1041.60 × 1043.93 × 1026.52 × 1039.00 × 100
Smpn0092302α,3α,23-trihydroxyolean-12-en-28-oic acidTriterpene4.87 × 1024.87 × 10223.34 × 1041.88 × 1065.09 × 1052.27 × 1052.01 × 1065.99 × 105
Notes: Index, Metabolite number, Class II, secondary classification of substances. Q1, the molecular weight of the parent ion after the substance was added to the ion by the electrospray ion source, Q3, the characteristic fragment ion, level, substance identification level, “1”, sample substance secondary mass spectra, RT and database substance matching score of 0.7 or more; “2”, sample substance secondary mass spectra, RT and database substance matching score of 0.5–0.7, “3”, sample substance’s five detection parameters Q1, Q3, RT, DP, CE, and database substance check are consistent. The values in the table are the relative content of metabolites, without units, calculated by calculating the peak area formed in the detector by the characteristic ions of each substance. It is not the absolute content of the substance, but the detection conditions are consistent and can be used to compare the differences in the same substance in different samples.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qin, L.; Du, F.; Yang, N.; Zhang, C.; Wang, Z.; Zheng, X.; Tang, J.; Yang, L.; Dong, C. Transcriptome Analyses Revealed the Key Metabolic Genes and Transcription Factors Involved in Terpenoid Biosynthesis in Sacred Lotus. Molecules 2022, 27, 4599. https://doi.org/10.3390/molecules27144599

AMA Style

Qin L, Du F, Yang N, Zhang C, Wang Z, Zheng X, Tang J, Yang L, Dong C. Transcriptome Analyses Revealed the Key Metabolic Genes and Transcription Factors Involved in Terpenoid Biosynthesis in Sacred Lotus. Molecules. 2022; 27(14):4599. https://doi.org/10.3390/molecules27144599

Chicago/Turabian Style

Qin, Lili, Fei Du, Ningning Yang, Chen Zhang, Zhiwen Wang, Xingwen Zheng, Jiawei Tang, Liangbo Yang, and Chen Dong. 2022. "Transcriptome Analyses Revealed the Key Metabolic Genes and Transcription Factors Involved in Terpenoid Biosynthesis in Sacred Lotus" Molecules 27, no. 14: 4599. https://doi.org/10.3390/molecules27144599

APA Style

Qin, L., Du, F., Yang, N., Zhang, C., Wang, Z., Zheng, X., Tang, J., Yang, L., & Dong, C. (2022). Transcriptome Analyses Revealed the Key Metabolic Genes and Transcription Factors Involved in Terpenoid Biosynthesis in Sacred Lotus. Molecules, 27(14), 4599. https://doi.org/10.3390/molecules27144599

Article Metrics

Back to TopTop