Curcuma wenyujin is the source plant of three traditional Chinese medicines, which have been widely used in clinical treatment over 1000 years. The content of terpenes, the major medicinal active ingredients, is relatively low in this plant. Studies have shown that MeJA can promote terpenes biosynthesis in plants. However, the mechanism underlying the effect of MeJA in C. wenyujin remains unclear. In this work, the transcriptome of C. wenyujin leaves with MeJA treatment was analyzed to elucidate the regulation mechanism of MeJA-mediated terpene biosynthesis. Based on the RNA-seq data, 7,246 unigenes were differentially expressed with MeJA treatment. Expression pattern clustering of DEGs revealed that unigenes, related to JA biosynthesis and signal transduction, responded to exogenous MeJA stimulation on the early stage and maintained throughout the process. Subsequently, unigenes related to terpene biosynthesis pathway showed a significant up-regulation with 6 h treatment. The analysis results suggested that MeJA induced the expression of JA biosynthesis genes (such as LOXs, AOSs, AOCs, OPRs, and MFPs) and JA signal transduction core genes (JAZs and MYCs) to activate JA signaling pathway. Meanwhile, downstream JA-responsive genes presented up-regulated expression levels such as AACT, HMGSs, HMGRs, DXSs, DXRs, MCTs, HDSs, and HDRs, thus promoting terpenes biosynthesis. The transcriptional expressions of these genes were validated by qRT-PCR. In addition, six CwTPS genes in response to MeJA were identified. With MeJA treatment, the expression levels of CwTPSs were increased as well as those of the transcription factors MYB, NAC, bZIP, WRKY, AP2/ERF, and HLH. These TFs might potentially regulate terpenes biosynthesis. These results provide insights for regulation mechanism of terpenes biosynthesis.

MYB superfamily is one of the most abundant families in plants, and plays critical role in plant growth, development, metabolism regulation, and stress response. Curcuma wenyujin is the main source plant of three traditional Chinese medicines, which are widely used in clinical treatment due to its diverse pharmacological activities. In present study, 88 CwMYBs were identified and analyzed in C. wenyujin, including 43 MYB-related genes, 42 R2R3-MYB genes, two 3R-MYB genes, and one 4R-MYB gene. Forty-three MYB-related proteins were classified into several types based on conserved domains and specific motifs, including CCA1-like type, R-R type, Myb-CC type, GARP-like type, and TBR-like type. The analysis of motifs in MYB DBD and no-MYB regions revealed the relevance of protein structure and function. Comparative phylogeny analysis divided 42 R2R3-MYB proteins into 19 subgroups and provided a reference for understanding the functions of some CwMYBs based on orthologs of previously characterized MYBs. Expression profile analysis of CwMYB genes revealed the differentially expressed genes responding to various abiotic stresses. Four candidate MYB genes were identified by combining the results of phylogeny analysis and expression analysis. CwMYB10, CwMYB18, CwMYB39, and CwMYB41 were significantly induced by cold, NaCl, and MeJA stress treatments. CwMYB18 and CwMYB41 were proved as regulators with activity of transcriptional activation, whereas CwMYB39 and CwMYB10 were not. They may participate in the response to abiotic stresses through different mechanisms in C. wenyujin. This study was the first step toward understanding the CwMYB family and the response to abiotic stresses in C. wenyujin.