## 2 Results
### Identification of Unigenes Related to the Biosynthetic Pathways of Terpenoids and Polysaccharides  
In this study, 189 candidate unigenes related to the biosynthetic pathways of terpenoids were screened and identified from the _D. nobile_ transcriptome. For genes related to the biosynthesis of terpenoids, 52, 39, 46, and 52 unigenes were involved in terpenoid backbone (ko00900), sesquiterpenoid and triterpen (ko00909), monoterepenoid (ko00902), and diterpenoid (ko00904) biosynthesis, respectively. Among them, seven unigenes were significantly upregulated in DF14; these genes were the (\(-\))-\(\alpha\)-terpineol synthase gene (_Loc110100458_, _Loc114579671_, _Loc110092374_, _Loc110092382_, and _Loc110097122_), the geraniol 8-hydroxylase gene (_LOC110098600_), and the squalene monooxygenase gene (_LOC110095241_). The correlation analysis (\(\mid\) r \(\mid\)\(\geq\) 0.8) showed that a total of 238 transcription factors (TFs) were identified that may participate in the regulation of sesquiterpenoid and/or triterpenoid biosynthesis. These TFs included twenty-four MYB and MYB-related  
Figure 5: Functional GO and KEGG pathway classification of differentially expressed genes. **(A)** The horizontal axis represents the enriched GO term, and the vertical axis represents the number of differentially expressed genes. **(B)** The three main categories of the directed acyclic graphs are displayed in the thumbnail view. The node colors are based on the q-values. The GO terms appear at the horizontal node locations; the darker the color, the higher the confidence level. TFs, sixteen zinc finger proteins, one ZAT11, one YABBY, eight WRKY, two TINY, four TCP, one SIARF2A, one SBH1, one RINGLET, two REVEILLE, two PTI, two PRR, two protein indeterminate domains, two HK3, one PEP, one PCF, six NF-YA/B/C, three SNAT2, thirteen NAC, one MYC2, one multiprotein-bridging factor 1c, two MTERF, one MADS, one lysine-specific demethylylase JMJ30, three AUX/IAA, one LOB, one LATERAL ROOT PRIMORDIUM 1, three KNAP/KAN, four JAZ, twelve HD-ZIP, six Hsf, one GT-2, one GRF6, four GRAS, one GLK1, one GHD7, five GATA, one GAG2, two FAR1, eight ERF, seven EREBP, one E2 FB, two DREBA, ten DOF, one DIVARICATA, one CPRF2, one BZR1, eight bZIP, six BTB/POZ, eight bHLH, four BEL1, five B-box, six B3, one AS2, one YABBY 2, two TCP, one RAD5A/SMARCA3, two PCF, three ORR, one NPR5, one ATXR7, and eleven other TFs (Supplementary Table S2). Of these, after stimulation with MeJA, the TFs that were increased were mainly the WRKY, NAC, MYB, HD-ZIP, ERF/EREBP, bZIP, bHLH, BTB/POZ, and DOF TF families. In this paper, 99 candidate unigenes related to the biosynthetic pathways of polysaccharides were screened and identified from the _D. nobile_ transcriptome. For genes related to the biosynthesis of polysaccharides, 56, 24, 22 and 20 unigenes participated in starch and sucrose metabolism (ko00500), amino sugar and nucleotide sugar metabolism (ko00520), fructose and mannose metabolism (ko00051), and galactose metabolism (ko00052), respectively. The correlation analysis showed that a total of 313 TFs were identified that may participate in the regulation of polysaccharide biosynthesis.