Regulation Mechanism of MYC Family Transcription Factors in Jasmonic Acid Signalling Pathway on Taxol Biosynthesis

Paclitaxel is an important anticancer drug. The phytohormone jasmonic acid can significantly induce the biosynthesis of paclitaxel in Taxus, but the molecular mechanism has not yet been resolved. To establish the jasmonic acid signalling pathway of Taxus media, based on the gene of the jasmonic acid signalling pathway of Arabidopsis thaliana, sequence analysis was performed to isolate the jasmonic acid signal from the transcriptome, a transcriptional cluster of pathway gene homologs and the full length of 22 genes were obtained by RACE PCR at 5′ and 3′: two EI ubiquitin ligase genes, COI1-1 and COI1-2;7 MYC bHLH type transcription factor (MYC2, MYC3, MYC4, JAM1, JAM2, EGL3, TT8); 12 JAZ genes containing the ZIM domain; and MED25, one of the components of the transcriptional complex. The protein interaction between each were confirmed by yeast two hybridization and bimolecular fluorescence complementation based on similar genes interaction in Arabidopsis. A similar jasmonate signaling pathway was illustrated in T. media. All known paclitaxel biosynthesis genes promoters were isolated by genome walker PCR. To investigate the jasmonate signaling effect on these genes’ expression, the transcription activity of MYC2, MYC3 and MYC4 on these promoters were examined. There are 12, 10 and 11 paclitaxel biosynthesis genes promoters that could be activated by MYC2, MYC3 and MYC4.


Introduction
Paclitaxel is an important anticancer drug, which is mainly produced by chemical semi-synthesis by using paclitaxel, which is biosynthesized in Pacific yew. Due to its complex molecular structure, its raw materials have been dependent on Taxus biomass, and its biosynthesis and regulation have been studied [1]. The phytohormone jasmonic acid can significantly induce the biosynthesis of paclitaxel in Taxus, but the molecular mechanism has not yet been resolved [2][3][4].
The biosynthetic pathway of paclitaxel is very complex, involving at least 19 steps of enzyme-catalysed reactions, and the related enzymes and synthetic pathways have not been fully resolved. The synthetic pathway can be divided into the following four steps: taxane ring synthesis; side hydroxylation, acetylation and phenylacetylation modification; epoxidation modification at C5 position and oxidation at C9; and side chain at C13 position addition [5].
The biosynthesis step of paclitaxel is very complicated, and it is very difficult to directly increase the expression of the biosynthetic pathway gene to increase the paclitaxel content. As the analysis of all synthetic pathways is not complete, the metabolic pathways cannot be precisely regulated to control the metabolic flux. In this case, for such a complex metabolic pathway, a better choice for regulating genes in the synthetic pathway is to regulate upstream transcription factors, thereby comprehensively regulating the expression of genes in the synthetic pathway [9]. It is important to analyse the transcription factors that regulate the expression of paclitaxel synthesis pathway genes.
Transcription factors regulating Taxol biosynthesis in Taxus are being closely studied [2][3][4]. Previous studies have used two methods to isolate related transcription factors: one is to first isolate the promoter of the key gene of the Taxol synthesis pathway, and then to study the transcription factor that can be combined with the promoter; the other is to study jasmonic acid-induced transcription factors in order to analyse whether it can regulate gene expression of paclitaxel.
In the transcription factor-related studies in Taxus, no relevant research results of jasmonic acid signaling pathways in model plants have been used. There is also no systematic and in-depth discussion of the signaling pathway of jasmonic acid in Taxus. Li et al. used high-throughput sequencing to study the transcriptome of jasmonic acid treatment for 16 h, and obtained an AtJAZ-like gene induced by jasmonic acid, and a similar AtMYC2 gene, but no in-depth study [10]. In the transcriptome sequencing study published by the Sun et al., the related gene expression after eight days of jasmonic acid treatment was involved, but no jasmine similar to that in the model plants induced in the early stage of jasmonic acid treatment was involved in the ketoacid signalling pathway transcription factor [11]. Lenka et al. published article on AtMYC2-like genes [8]; except for the TcMYC2 gene, which is known to have a full-length gene, other transcription factors were not found to be similar to AtMYC2 in model plants.
The AtJAZ-binding region also does not contain active regions that activate transcription. Maybe due to the shortcomings of the adopted technical methods, the conclusions obtained were also different from those of similar genes in model organisms [12].
In this study, high-throughput sequencing was used to obtain transcription factors after an early induction with jasmonic acid. In our previous study, the sequence of several genes related to the paclitaxel synthesis pathway and jasmonic acid signaling pathway were obtained by sequencing analysis of the transcriptome and expression profile of Taxus. Based on this, the full-length sequence of jasmonic acid signalling pathway-related transcription factors in T. media was obtained by sequence comparison analysis and RACE PCR. It includes sequences similar to the important genes of the jasmonic acid signalling pathways such as AtCOI1, AtMYC2, AtJAZ and AtMED25 (At, Arabidopsis thaliana, if the source of the gene is not clearly prefixed, it is from the T. media). The interaction between JAZ and MYC transcription factors and other genes was studied, and the jasmonic acid signalling pathway in Taxus was confirmed. The above work further studied the regulation of the promoter of the paclitaxel synthesis pathway gene by the MYC transcription factor, and elucidated the regulation mechanism of the transcription factor of the jasmonic acid signaling pathway on the biosynthesis pathway of paclitaxel.
We also found that some of the genes involved in the paclitaxel synthesis pathway were highly expressed by jasmonic acid at 0.5 h, similar to the transcription factors of the jasmonic acid signalling pathway and the expression of self-synthesis-related genes, but faster. Expression of anthocyanin synthesis-related genes regulated by the jasmonic acid signalling pathway was noted at three hours [13]. This indicated that the jasmonic acid signalling pathway plays an important role in regulating the gene expression of the paclitaxel synthesis pathway.
To analyse the mechanism of jasmonic acid signalling pathway specifically regulating the paclitaxel synthesis pathway gene, 14 promoters of paclitaxel synthesis pathway gene were isolated, and five MYC transcription factors were verified by fluorescence activity assay in Arabidopsis protoplasts. Regulation of the promoter of the paclitaxel biosynthetic pathway gene and the transcription factor mutants with application value were designed, which lays the foundation for the future construction of a high-yield paclitaxel-based Taxus cell line.

JAZ
At least eight JAZ-containing ZIM domains were obtained in the differentially expressed genes by the transcriptome sequencing. For comprehensive detection of JAZ that may be contained in the transcriptome, BLAST analysis was performed using the full length of the Arabidopsis AtJAZ1 gene and the eight JAZs with reference to the Tm_ZK transcriptome, and the threshold of the E value was set to 0.01. Finally, 16 transcription cluster sequences which may be JAZ were obtained, as shown in Table 1. For the 16 possible JAZs, specific primers were designed and the full length of the gene was obtained by RACE PCR at the 3 and 5 ends. Three of the transcription clusters comp60304_c0, comp64513_c0, and comp68543_c0 failed to obtain a full-length sequence. The other 13 obtained the full length sequence, and wherein comp70749_c0 was a non-coding sequence. The remaining 12 transcription clusters obtained the encoded amino acid sequence-see Table 2. Homology analysis and domain analysis in the amino acid sequence were performed (see Figures 2a and A1), and the 12 genes were named JAZ1 to JAZ12.
The primers were designed to amplify the 12 full-length genes, and the cDNA was used as a template for PCR. The bands obtained from the cloning were sequenced to obtain 18 sequences. Among the 18, JAZ1, JAZ2, JAZ9, JAZ11 and JAZ12 all contained varying cleavages. There were three in JAZ2, and the other had two coding sequences (see the Appendices A and B for the specific sequences). Similar to the variable cleavage of the JAZ gene in model plants, it is important for the diversity of the regulatory functions of the JAZ gene. comp64513_c0, and comp68543_c0 failed to obtain a full-length sequence. The other 13 obtained the full length sequence, and wherein comp70749_c0 was a non-coding sequence. The remaining 12 transcription clusters obtained the encoded amino acid sequence-see Table 2.
Homology analysis and domain analysis in the amino acid sequence were performed (see Figure 2a and A1), and the 12 genes were named JAZ1 to JAZ12.

COI1
BLAST alignment analysis of the Tm_ZK transcriptome revealed two AtCOI1-like genes, named COI1. To detect the AtCOI1-like gene that may be contained in the transcriptome, BLAST analysis was performed using the full length of Arabidopsis AtCOI1 and the above two COI1 with reference to transcriptome, and the threshold of the E value was set to 0.00001. Only the two sequences were obtained, and comp78658_c0 and comp69536_c0 may be the transcription cluster sequence of COI1.
For the two possible COI1, specific primers were designed and the full length was obtained by RACE PCR at the 3 end and 5 end. It was found that the two COI1 were actually the full length of the single COI1 gene.
Specific primers to amplify COI1 were produced, and the two highly similar COI1, named COI1-1 and COI1-2, were cloned. Homology analysis was performed with AtCOI1 in Arabidopsis, and the results are shown in Figure A2. There were only 10 amino acid differences between COI1-1 and COI1-2, and most of them were not conservative amino acids. Only the amino acid at position 363 was a conservative amino acid variation.

MYC Transcription Factor
The MYC transcription factor belongs to the bHLH type. Ten bHLH-type transcription factors were obtained in the differentially expressed genes by transcriptome sequencing analysis. For comprehensive detection of MYC-like transcription factors that may be contained in the transcriptome, BLAST analysis was performed using the transcription factors of Arabidopsis AtMYC2 full length and the 10 bHLH types, and the threshold of the E value was set to 0.01. Finally, 19 possible transcription clusters were found, and homology analysis was performed to obtain full length of RACE for seven transcription clusters of comp78752_c0, comp74670_c0, comp77888_c0, comp77131_c0, comp75263_c0, comp77266_c0, and comp78399_c0. It was found that comp75263_c0 and comp78752_c0 were the sequences at both ends of the same gene ( Table 3).
The primers were then designed and cloned using a high-fidelity enzyme to confirm the final gene sequence. The comp77888_c0 was amplified in both cDNA and genomic DNA. In addition to the predicted length bands, a shorter band was obtained. The truncated comp77888_c0 we named comp77888s_c0, so we finally isolated and identified 7 AtMYC2 similar transcription factors in T.media: comp74670_c0, comp78752_c0, comp77888_c0, comp77888s_c0, comp77131_c0, comp78399_c0 and comp77266_c0 (see the Appendices A and B for the sequence). Among them, MYC2, MYC3, MYC4 and JAM1 had higher similarity with AtMYC2 in Arabidopsis (Figure 2b), which may be the transcription factor of the core jasmonic acid signalling pathway. The gene expression levels of MYC2, MYC3, and MYC4 were also induced by jasmonic acid, especially MYC3 and MYC4, which were similar to those in the model plants. However, JAM1, which is highly homologous to MYC4, lacked a domain that binds to JAZ and may be a transcription factor that inhibits jasmonic acid signalling.
JAM2 is similar to Arabidopsis thaliana's inhibitory bHLH-type transcription factor AtJAM1, which is also upregulated by the induction of methyl jasmonate treatment, which may be another transcription factor that inhibits the jasmonic acid signalling pathway.
The similarity of transcription factors between TT8 and EGL3 and the jasmonic acid signalling pathway core in model plants was low, and the amino acid sequence of the relevant domain was also different, which is similar to the transcription factor sequence in the regulation of anthocyanin synthesis in model plants. Like the model plants and not induced by the jasmonic acid signal, it is likely to be a transcription factor involved in the regulation of anthocyanin metabolism in Taxus.

MED25
In model plants, a transcription factor such as AtMYC2 binds to AtMED25 to recruit RNA polymerase and initiate transcription of a downstream gene of the jasmonic acid signal [14]. At MED25 is an important gene for transcription initiation of jasmonic acid downstream pathway genes and plays an important role in the jasmonic acid signalling pathway. The unique and full length MED25 was obtained in the Taxus transcriptome by BLAST analysis. Primers were designed and amplified with high fidelity enzyme to obtain MED25. See the Appendices A and B for the sequence.

Interaction between Jasmonic Acid Signalling Pathway Genes in Taxus
The jasmonic acid signal is transmitted from the receptor to the protein through the interaction between the proteins. It is necessary to further investigate and detect the jasmonic acid signal by exploring whether the genes of the jasmonic acid signalling pathway are interacting, and to study JAZ and COI1, JAZ and MYC, JAZ and JAZ, MYC and MYC, MYC and MED25 protein interactions, determine whether it is the same as in the model plants, and then construct the jasmonic acid signal in the Taxus.
Since the interaction of the genes has been verified by two methods in model plants, in this study, to detect the interaction between related proteins, the yeast two-hybrid method was used to obtain preliminary interactions between proteins.

JAZ and COI1
Jasmonyl isoleucine can mediate the binding of JAZ to COI1, thereby ubiquitin-labelling JAZ and degrading JAZ-released transcription factors. The combination of AtJAZ and AtCOI1 is an important basis for jasmonic acid signal reception in Arabidopsis. This study analysed 12 JAZ and 2 COI1 from the Taxus, and the analog of coronatine, jasmonyl isoleucine. The results are shown in Table 4 ("-"represents no signal, "+" represents a signal, the number of "+"represents the signal's strength). Table 4. Yeast two-hybrid analysis. The results of the interaction between JAZ and COI1.
There is also a similar combination of JAZ protein and COI1 protein in the Taxus. However, it is different from the results of the model plant Arabidopsis.

JAZ and MYC
In model plants, AtJAZ binds to AtMYCs, thereby inhibiting the transcriptional activity of AtMYCs. In this study, the binding of all JAZ to MYC in Taxus was analysed by the yeast two-hybrid method, and the interaction of some JAZ and MYC in tobacco was verified by using bimolecular fluorescence complementation (see Figure 3). Jasmonyl isoleucine can mediate the binding of JAZ to COI1, thereby ubiquitin-labelling JAZ and degrading JAZ-released transcription factors. The combination of AtJAZ and AtCOI1 is an important basis for jasmonic acid signal reception in Arabidopsis. This study analysed 12 JAZ and 2 COI1 from the Taxus, and the analog of coronatine, jasmonyl isoleucine. The results are shown in Table 4("-"represents no signal, "+" represents a signal, the number of "+"represents the signal's strength).
There is also a similar combination of JAZ protein and COI1 protein in the Taxus. However, it is different from the results of the model plant Arabidopsis.

JAZ and MYC
In model plants, AtJAZ binds to AtMYCs, thereby inhibiting the transcriptional activity of AtMYCs. In this study, the binding of all JAZ to MYC in Taxus was analysed by the yeast two-hybrid method, and the interaction of some JAZ and MYC in tobacco was verified by using bimolecular fluorescence complementation (see Figure 3) (a)   Yeast two-hybrid results showed that transcription factors MYC2, MYC3 and MYC4 interact with most of the JAZ proteins. The exceptions with no interactions are present between MYC2 and JAZ1, MYC3 and JAZ1 and JAZ8, MYC4 and JAZ8. Like the predicted results, JAM1 lacking the JID domain has no interaction with 11 JAZs, but, unexpectedly, it can interact with JAZ8. Another inhibitory transcription factor, JAM2, interacts with JAZ4~JAZ12.
Bimolecular fluorescence complementation verified the interaction between JAZ6 and MYC2, MYC3, MYC4 and JAM2. JAZ also interacts with MYC type and JAM type transcription factors.

JAZ and JAZ
The interaction between AtJAZ and the others is an important condition for the transcription factor inhibition complex, and enriches the means of plant regulation of downstream pathway gene expression through the interaction between differently cleaved JAZ proteins.
The yeast two-hybrid assay was used in this study to analyse the interaction between JAZ in Taxus ( Figure 4a). domain has no interaction with 11 JAZs, but, unexpectedly, it can interact with JAZ8. Another inhibitory transcription factor, JAM2, interacts with JAZ4~JAZ12. Bimolecular fluorescence complementation verified the interaction between JAZ6 and MYC2, MYC3, MYC4 and JAM2. JAZ also interacts with MYC type and JAM type transcription factors.

JAZ and JAZ
The interaction between AtJAZ and the others is an important condition for the transcription factor inhibition complex, and enriches the means of plant regulation of downstream pathway gene expression through the interaction between differently cleaved JAZ proteins.
The yeast two-hybrid assay was used in this study to analyse the interaction between JAZ in Taxus (Figure 4a).  Summarize the interaction between JAZ proteins as shown in Table 5("-"represents no signal, "+" represents a signal, the number of "+"represents the signal's strength):  Summarize the interaction between JAZ proteins as shown in Table 5 ("-"represents no signal, "+" represents a signal, the number of "+" represents the signal's strength):  -

MYC and MYC
The MYC transcription factor contains the domain of bHLH, in which the domain of HLH is responsible for the binding of MYC transcription factors to each other, forming a homologous or heterologous bimolecular structure in the plant to regulate the transcription of downstream genes. There is no interaction between activated MYC-like transcription factors and inhibitory MYC-like transcription factors in model plants. In this study, the interaction between MYC2, MYC3, MYC4 and JAM1, JAM2 of Taxus was analysed. The results are shown in Figure 4b.
Yeast two-hybrid results showed an interaction between MYC2 and MYC3, JAM2 interacted with itself, and there was no interaction between other MYC-type transcription factors.

MYC and MED25
The MYC-like transcription factor binds to MED25, thereby recruiting the RNA polymerase complex and activating transcription. Yeast two-hybrid study of the interaction between MYC2 and other transcription factors in Taxus and MED25; the results are shown in Figure 5a. The MYC transcription factor contains the domain of bHLH, in which the domain of HLH is responsible for the binding of MYC transcription factors to each other, forming a homologous or heterologous bimolecular structure in the plant to regulate the transcription of downstream genes. There is no interaction between activated MYC-like transcription factors and inhibitory MYC-like transcription factors in model plants. In this study, the interaction between MYC2, MYC3, MYC4 and JAM1, JAM2 of Taxus was analysed. The results are shown in Figure 4b.
Yeast two-hybrid results showed an interaction between MYC2 and MYC3, JAM2 interacted with itself, and there was no interaction between other MYC-type transcription factors.

MYC and MED25
The MYC-like transcription factor binds to MED25, thereby recruiting the RNA polymerase complex and activating transcription. Yeast two-hybrid study of the interaction between MYC2 and other transcription factors in Taxus and MED25; the results are shown in Figure 5a. The activated transcription factors MYC2, MYC3 and MYC4 in Taxus can bind to MED25, while JAM1 has no interaction with MED25. Surprisingly, JAM2 binds to MED25, indicating that it may have some transcriptional activation activity. EGL3, which may be involved in anthocyanin synthesis, may also bind to MED25, indicating that it may also have transcriptional activation activity.

Detection and Analysis of MYC Transcription Factor Activity
To further confirm the activity of the MYC-type transcription factor, the transcriptional activity of the relevant transcription factors was verified in Arabidopsis protoplasts (see Figure A4). GAL4BD, The activated transcription factors MYC2, MYC3 and MYC4 in Taxus can bind to MED25, while JAM1 has no interaction with MED25. Surprisingly, JAM2 binds to MED25, indicating that it may have some transcriptional activation activity. EGL3, which may be involved in anthocyanin synthesis, may also bind to MED25, indicating that it may also have transcriptional activation activity.

Detection and Analysis of MYC Transcription Factor Activity
To further confirm the activity of the MYC-type transcription factor, the transcriptional activity of the relevant transcription factors was verified in Arabidopsis protoplasts (see Figure A4). GAL4BD, which binds to the GAL4 promoter, was fused to a transcription factor and the transcriptional activities of MYC2, MYC3, MYC4, JAM1 and JAM2 in protoplasts were examined.
Among them, the positive control AtMYB21 had strong transcriptional activation activity, and the negative control AtbHLH17 had no transcriptional activation activity. The transcription factor MYC2 in Taxus had strong transcriptional activation, and the transcriptional activation activities of MYC3 and MYC4 were low. While JAM1 and JAM2 had almost no activation activity, the results of yeast two-hybridization are consistent with the fact that JAM2 is more active than JAM1.

Synthetic Pathway Gene Promoter Separation and Analysis
To analyse the regulation of transcription factor on the expression of paclitaxel synthesis pathway gene, it is important to isolate the promoter of the synthesis pathway gene. In the existing literature, only a few promoters of the paclitaxel synthesis pathway have been isolated but they were not long enough for the analysis of the promoter properties [2,10].
Therefore, a Genome walker PCR was used to isolate the promoters of all 14 known genes in the Taxol synthesis pathway in Taxus. The sequence is shown in the Appendices A and B. To investigate the possible regulation of these genes by MYC transcription factors, we analysed the number of E-box (CANNTG) that can bind to MYC transcription factors, especially the G-box which has the strongest binding ability to MYC transcription factors. CACGTG sequences and variant G-boxes (AACGTG, CACGAG and special E-box, CATGTG CATGTG) are defined as High Affinity (HA) E-box; AACGTG and CACGAG are defined as modified G-box, and other E-boxes are defined as low affinity (LA) E-box [15]. The number of promoters in the gene was analysed, and the promoter information of the obtained gene is shown in Table 6: Table 6. List of cis-element analysis information in the paclitaxel biosynthesis pathway gene promoter sequence. The 14 promoters of th synthetic pathway genes all contain E-box, and all of them are induced by jasmonic acid, and the promoters of six genes contain G-box. The TBT promoter contains the most E-box and the most cis-elements with high affinity for MYC-like transcription factors. The 14βOH promoter contains the least amount of E-box and is low-affinity. The specific location of the cis component is shown in Figure A5.

Activated MYC Transcription Factor Induces a Paclitaxel Biosynthesis Pathway Gene Promoter
Gel electrophoresis migration experiments demonstrated that full-length MYC2 binds to G-box and high-affinity E-box [2]. This indicates that MYC transcription factors have the potential to activate paclitaxel synthesis pathway gene expression. To study the regulation of MYC transcription factors on metabolic pathway genes in this study, plasmids with fluorescent protein and related gene promoter that continuously express transcription factors were transfected into Arabidopsis protoplasts.
The regulation of MYC2, MYC3 and MYC4 on the promoter of paclitaxel synthesis pathway gene was examined. The principle is shown in Figure A6.
For each promoter, the fluorescence intensity data ratio was normalized to the 62sk plasmid as a reference, and the results were as Figure 6: activate paclitaxel synthesis pathway gene expression. To study the regulation of MYC transcription factors on metabolic pathway genes in this study, plasmids with fluorescent protein and related gene promoter that continuously express transcription factors were transfected into Arabidopsis protoplasts. The regulation of MYC2, MYC3 and MYC4 on the promoter of paclitaxel synthesis pathway gene was examined. The principle is shown in Figure A6. For each promoter, the fluorescence intensity data ratio was normalized to the 62sk plasmid as a reference, and the results were as Figure 6: (a) The activation of the same transcription factor for the promoters of different genes is shown in Figure 7a,b,c: (a) In comparison, among the three activating transcription factors of the jasmonic acid signalling pathway, MYC2, MYC3, and MYC4, although they have different effects on each promoter, MYC2 can activate the expression of all promoter-induced genes. MYC3 and MYC4 had no induction effect on the promoter of BAPT and 14βOH, and both had weaker activation effects on the promoters of TAT and TBT.
The promoters of the genes for the paclitaxel synthesis pathway are substantially similar in relative strength induced by the above transcription factors. The classification is as follows: the promoter of DBTNBT, 2αOH, 7βOH, TS is the strongest; followed by 5αOH, 13αOH, DBAT; followed by PAM, 10βOH, GGPPS; while the promoters of TAT and TBT and BAPT, 14βOH are the weakest. For the promoters of the latter two genes, two of them have an effect of inhibiting transcription.

Inhibitory MYC Transcription Factor Inhibits Activation MYC Transcription Factor
The model plant Arabidopsis has four inhibitory MYC transcription factors that bind to the promoter by competitively binding to the activated MYC transcription factor, thereby inhibiting the action of the activated MYC transcription factor [16,17]. For JAM1 and JAM2 isolated from Taxus, it was found that there was no corresponding transcriptional activation in protoplasts. The result is shown in Figure 7d,e. In the future, JAM1 and JAM2 overexpression plasmid and reporter genes regulated by all gene promoters should be transferred into protoplasts, and further investigated for it can activate the transcription of the promoter of paclitaxel synthesis pathway gene.
JAM1 and JAM2 had no significant transcriptional activation of all of the tested promoters. Yeast two-hybrid results showed that the inhibitory MYC transcription factor could not bind to the activated MYC transcription factor, similar to that in the model plant. Whether or not it is similar to the model plant, it inhibits the action of the activated transcription factor by competitively binding the promoter sequence to the activated MYC-like transcription factor. Subsequent studies will co-transfect plasmids expressing inhibitory and activating transcription factors in protoplasts, investigate the strength of transcriptional activity, and analyse whether they have similar inhibitory effects and principles.

Mutant MYC Transcription Factor Activity is not Affected by JAZ
JAZ is an important protein in the signal reception process of jasmonic acid, and it is also a protein that inhibits the activity of transcription factors, and it is also induced by jasmonic acid signal, which plays a feedback inhibition effect in the jasmonic acid signalling pathway. A mutant MYC-like transcription factor in a model plant that does not interact with JAZ but activates transcription of a downstream pathway gene, providing an idea for using biotechnology to increase the signalling pathway of jasmonic acid.
In Taxus, JAZ also inhibits the expression of downstream pathway genes by binding to activated transcription factors. This feedback inhibition is an important means by attenuating the effects of JAZ to attenuate the jasmonic acid signaling pathway. However, the number of JAZ in Taxus is high, and the way of mutating JAZ to regulate downstream signalling pathways is complicated, and mutations and their interacting MYC transcription factors are feasible methods.
By analysing the complex structure of AtJAZ and AtMYC3 in model plants, the MYC2 and MYC3 variants were designed by MYC2YL, MYC2D, MYC2WDY and MYC3YL, which are the conserved amino acids that play an important role in the binding of JAZ-like transcription factors to JAZ. MYC3D and MYC3WDY (for detailed mutation sites, refer to the sequence of each transcription factor variant in the Appendices A and B).
First, the binding of the above six MYC transcription factor variants to MED25 was analysed by yeast two-hybrid. The result is shown in Figure 5b. From comparison of results from the unmutated MYC experimental group, it can be seen that, apart from MYC2YL which does not bind to MED25, other mutations do not affect the binding of MYC transcription factors to MED25, indicating that the transcriptional activities of other mutant MYC transcription factors are likely to be unaffected.
The interaction between 12 JAZ and mutant MYC was further studied by yeast two-hybrid assay to analyse whether the mutant MYC could bind to JAZ. The result is shown in Figure 8. Taxus is high, and the way of mutating JAZ to regulate downstream signalling pathways is complicated, and mutations and their interacting MYC transcription factors are feasible methods. By analysing the complex structure of AtJAZ and AtMYC3 in model plants, the MYC2 and MYC3 variants were designed by MYC2YL, MYC2D, MYC2WDY and MYC3YL, which are the conserved amino acids that play an important role in the binding of JAZ-like transcription factors to JAZ. MYC3D and MYC3WDY (for detailed mutation sites, refer to the sequence of each transcription factor variant in the appendix).
First, the binding of the above six MYC transcription factor variants to MED25 was analysed by yeast two-hybrid. The result is shown in Figure 5b. From comparison of results from the unmutated MYC experimental group, it can be seen that, apart from MYC2YL which does not bind to MED25, other mutations do not affect the binding of MYC transcription factors to MED25, indicating that the transcriptional activities of other mutant MYC transcription factors are likely to be unaffected.
The interaction between 12 JAZ and mutant MYC was further studied by yeast two-hybrid assay to analyse whether the mutant MYC could bind to JAZ. The result is shown in Figure 8. Three forms of mutations all affect the protein interaction between JAZ and MYC transcription factors. Mutation of a single D amino acid site affects the binding of MYC2 and MYC3 to most JAZ, while the mutation of three amino acids of WDY almost prevents MYC2 and MYC3 from binding to all JAZ proteins. The mutation of two amino acids of YL has a stronger effect in MYC2, which makes MYC2 unable to bind to all JAZ, but also cannot bind to MED25; the YL mutation of MYC3 makes it have no interaction with any JAZ except JAZ7.
In summary, the WDY mutation of MYC2 and the YL mutation of MYC3 may be novel transcription factors that are not inhibited by JAZ but have strong transcriptional activity and can be used to improve gene expression of paclitaxel biosynthesis pathway.
The regulation of the MYC-like transcription factor of the mutant on the paclitaxel biosynthetic pathway promoter will be further verified in protoplasts. Three forms of mutations all affect the protein interaction between JAZ and MYC transcription factors. Mutation of a single D amino acid site affects the binding of MYC2 and MYC3 to most JAZ, while the mutation of three amino acids of WDY almost prevents MYC2 and MYC3 from binding to all JAZ proteins. The mutation of two amino acids of YL has a stronger effect in MYC2, which makes MYC2 unable to bind to all JAZ, but also cannot bind to MED25; the YL mutation of MYC3 makes it have no interaction with any JAZ except JAZ7.

Taxus has a Similar Mechanism of Signal Accepting and Conduction of Jasmonic Acid
In summary, the WDY mutation of MYC2 and the YL mutation of MYC3 may be novel transcription factors that are not inhibited by JAZ but have strong transcriptional activity and can be used to improve gene expression of paclitaxel biosynthesis pathway.
The regulation of the MYC-like transcription factor of the mutant on the paclitaxel biosynthetic pathway promoter will be further verified in protoplasts.

Taxus Has a Similar Mechanism of Signal Accepting and Conduction of Jasmonic Acid
As an important phytohormone, jasmonic acid is an important component of plant stress response. It has been proven to be ubiquitous in plants, not only in the model plant Arabidopsis thaliana [14,[18][19][20][21]. Similar acceptance and transmission pathways for jasmonic acid signals should also be present in the Taxus. Prior to this study, a small number of genes for the jasmonic acid signalling pathway have been isolated by high-throughput sequencing [10], but have not conducted in-depth studies.
The jasmonic acid signalling pathway in Taxus was confirmed by high-throughput sequencing and identification of related genes and their function. Twenty-one isolated genes of the jasmonic acid signalling pathways were similar to those found in the model plants, considering characteristic protein domains and protein-protein interactions. It is shown that a similar jasmonic acid signalling pathway may exist in Taxus and play an important role.
Firstly, 12 JAZ interacted with two COI1 in the presence of coronatine, and the interaction intensity of some JAZ and COI1 increased with the increase of coronin concentration, indicating that these are probably active molecules of jasmonic acid signalling in Taxus. Jasmonyl isoleucine also mediates the interaction between the JAZ and COI1. However, JAZ8 and JAZ9 also have strong interaction with COI1 in the absence of coronatine, which is different from that in model plants. This indicates that the signalling pathways in Taxus may be different. The relevant biological significance requires in-depth research.
There is also an interaction between the protein of MYC transcription factor and the JAZ protein in Taxus, which is an important way for JAZ to inhibit the transcriptional activities of MYC2, MYC3 and MYC4. JAM1, lacking the JID domain, does not bind to most JAZ, indicating that the mutual binding is based on the JID domain in the MYC-like transcription factor. JAZ8, which contains two EAR domains and has a long length and JAS domain, is not combined with MYC3 and MYC4, but can be combined with MYC2 and the JAM1 without the JID domain. It may play a special role in the jasmonic acid signalling pathway of Taxus. Subsequent structural analysis of the complexes formed by AtMYC3 and AtJAZ1 of A. thaliana resulted in mutations in MYC2 and MYC3, and mutations in the JID domain resulted in their inability to bind to JAZ. This also suggests that the jasmonic acid signalling pathway in Taxus is similar to that in Arabidopsis [16,[22][23][24]. Interactions exist between JAZ proteins, especially between any combinations of any two of the six JAZs from JAZ7 to JAZ12, but there is less interaction between any combinations of any two of the six JAZs from JAZ1 to JAZ6, indicating that JAZ7-JAZ12 plays an important role in the formation of inhibitory protein complexes [22].
There are also interactions between MYC transcription factors, MYC2 and MYC3, and JAM2 interacts with itself. This differs from the broad interaction between inhibitory MYC transcription factors between activated MYC-like transcription factors in model plants.
The final MYC-type transcription factor recruits RNA polymerase by binding to MED25 to activate transcription of downstream genes. MYC2, MYC3 and MYC4 in Taxus have the same function. JAM1 without activation activity could not bind to MED25. However, JAM2, which has no apparent transcriptional activity, binds to MED25, indicating that it may bind to MED25 competitively, thereby inhibiting jasmonic acid signalling.

Jasmonic Acid Signalling Pathway Transcription Factor Regulates Paclitaxel Biosynthesis Pathway Gene Expression
The transcription factors MYC2, MYC3 and MYC4 of the jasmonic acid signalling pathway play an important role in the paclitaxel biosynthesis pathway. In the promoter region of most paclitaxel synthesis pathway genes, there are G-boxes and similar cis-elements that bind to MYC-type transcription factors, including the promoters of TAT, TBT, DBTNBT, 2αOH, 7βOH, 5αOH, 13αOH, DBAT, 10βOH and PAM, TS. The promoter of BAPT, GGPPS has fewer cis-elements, while the promoter of 14βOH has no corresponding cis-elements.
The activity of the three transcription factors also showed transcriptional activation activity against the promoters of DBTNBT, 2αOH, 7βOH, 5αOH, 13αOH, DBAT, 10βOH, PAM, GGPPS and TS, which contain cis-elements. However, there was no strong transcriptional activation activity for BAPT and 14βOH. It is shown that the G-box cis-element is important for the activity of the above three transcription factors. However, the promoters of TAT and TBT, which contain a large number of cis-like elements similar to G-box, suggest that MYC transcription factors should bind to the transcription factors of the two genes. However, in the Arabidopsis protoplasts, they are not strongly activated by MYC transcription factors; in particular, MYC3 and MYC4 have a weak regulation effect on them. Whether the excessive cis-elements in the promoters of TAT and TBT lead to inhibition requires further analysis.
JAM2 has a significant inhibitory effect on all paclitaxel synthesis pathway genes, and its expression is induced by jasmonic acid signalling. It is likely to be a feedback inhibition signal of the jasmonic acid signalling pathway in Taxus. The inhibitory effect of JAM1 is weaker than JAM2. JAM2 can be combined with MED25. JAM1 and MYC4 almost have the same sequence. Besides the above functions, whether they are different from model plants, these two transcription factors have other unknown functions and need further research.
In summary, the jasmonic acid signalling pathway in Taxus is summarized in Figure 9. BAPT and 14βOH. It is shown that the G-box cis-element is important for the activity of the above three transcription factors. However, the promoters of TAT and TBT, which contain a large number of cis-like elements similar to G-box, suggest that MYC transcription factors should bind to the transcription factors of the two genes. However, in the Arabidopsis protoplasts, they are not strongly activated by MYC transcription factors; in particular, MYC3 and MYC4 have a weak regulation effect on them. Whether the excessive cis-elements in the promoters of TAT and TBT lead to inhibition requires further analysis. JAM2 has a significant inhibitory effect on all paclitaxel synthesis pathway genes, and its expression is induced by jasmonic acid signalling. It is likely to be a feedback inhibition signal of the jasmonic acid signalling pathway in Taxus. The inhibitory effect of JAM1 is weaker than JAM2. JAM2 can be combined with MED25. JAM1 and MYC4 almost have the same sequence. Besides the above functions, whether they are different from model plants, these two transcription factors have other unknown functions and need further research.
In summary, the jasmonic acid signalling pathway in Taxus is summarized in Figure 9.  To confirm the function of the above jasmonic acid signalling pathway-related genes, the construction of related variants is the most important evidence, and the construction of the To confirm the function of the above jasmonic acid signalling pathway-related genes, the construction of related variants is the most important evidence, and the construction of the transgenic method of Taxus cell line is a worldwide problem, which is also related to this study. The study will validate the function of related genes and realize the construction of high-yield Taxus cell lines in the future.

Cell Culture, Jasmonate Treatment, and RNA Isolation
The T. media cell line Zike used in this study was induced from sapling leaf purchased in Beijing. The cells were subcultured using modified SSS medium. For jasmonate treatment, 25 µL 100 mM methyl jasmonate (MeJA) dissolved in ethanol was added to 2 g cells that were suspension cultured in 50 mL of medium 14 days after subculture. Samples were collected after 0, 0.5, 3, and 24 h MeJA treatment using a Büchner funnel. Three biologically independent experiments were performed and a total of 12 samples were obtained. Total RNA was isolated from each of the samples as described before [25]. Generally, the cells were frozen in liquid nitrogen and ground to powder. Every 2 g of cell powder was lysed in 10 mL GTC buffer (62.5 mM Tris-HCl, 12.5 mM EDTA, and 5 M guanidinium isothiocyanate, pH 8.0), then incubated on ice for 15 min and centrifuged at 15,000× g for 15 min. The supernatant was added to an equal volume of isopropanol to precipitate the total RNA. The total RNA was purified using an RNeasy Mini kit (Qiagen, City, Germany) according to the manufacturer's instructions.

RACE PCR
The full length sequence of the 5 and 3 ends of the gene was obtained using the GeneRacer Kit with SuperScript III RT kit based on Invitrogen's RACE PCR technology (CA, USA) using high quality RNA extraction. First, the extracted RNA is dephosphorylated, and then the enzyme is used to remove the 5 end cap of the RNA, and then the known RNA linker sequence is ligated. Using this as a template, the first strand of the cDNA is synthesized by the oligo(T) used as a locking primer to bind the the 3 end ploy(A) tail of mRNA. Design specific primers for nested PCR at the 5 and 3 ends, respectively, clone the obtained DNA fragments, analyse the sequence of the coding region after sequencing, and after splicing, both stop codons and reach the longest coding region to determine the full length sequence at both ends.

Yeast Two-Hybrid
Based on Clontech's MATCHMAKER LexA Two-Hybrid System, JAZ and other related genes are ligated into a pLexA vector containing a LexA promoter binding domain (BD), and another gene for detecting interaction is linked to transcriptional activation domain(AD) of the pB42AD vector. The above two plasmids were then transferred to the lacZ-expressing active EGY48 yeast strain containing LexA control using the Frozen-EZ Yeast Transformation Kit of Zymo Research. After overnight culture in SD/-His/-Trp/-Ura liquid medium, transfer 5-20 µL of yeast culture solution to SD/Gal/Raf/-His/-Trp/-Ura/ for activity detection. In X-gal + BU salts solid medium, the appearance of blue spots was observed after four days of culture to evaluate the interaction between lacZ activity and protein.

Bimolecular Fluorescence Complementation
The sequences of the two genes detecting the interaction were ligated into p35SnYFP and p35ScYFP, respectively, and the fusion protein was formed with the C-terminus and the N-terminus of the yellow fluorescent protein YFP, respectively. The assembled plasmids were then transferred to GV3101 Agrobaterium tumefaciens by electroporation. After induction with acetosyringone for two days, the two plasmids were simultaneously transferred into the tobacco leaves of the 4-week culture by injection, and then the fluorescence signal of the injection site on the leaves was observed using a laser confocal microscope.

Genome Walker PCR Separation Promoter Sequence
To obtain the promoter of the paclitaxel biosynthetic pathway gene, the genome of T. media was extracted using the Plant Genomic DNA Purification Kit (Tiangen, Beijing, China). Based on Clontech's Universal Genome Walker Kit (Clontech, Takara, Japan), genomic DNA was digested with four different blunt-end restriction enzymes DraI, EcoRV, PvuII, and StuI, followed by ligation to the linker. Primers were designed for two rounds of nested PCR by known coding region sequences and linker sequences, and the obtained DNA bands were cloned, sequenced, spliced and analysed to obtain the promoter sequences of the corresponding genes.

Detection of Transcription Factor Activity in Protoplasts
To analyse the transcriptional activity, the transcription factor was transferred to a plasmid 35SGAL4DB containing a GAL4DB domain in a plant cell to form a plasmid expressing a fusion of the GAL4DB domain and a transcription factor. Based on the mature Arabidopsis protoplast transformation method [26], the Arabidopsis protoplasts were transferred by a PEG transformation method, and a plasmid containing a GAL4 promoter-controlled transcription factor was simultaneously transformed. The plasmid also contained the reporter gene for the LUC fluorescent protein in the plasmid and the plasmid for the continuous expression of the REN fluorescent protein internal reference. The fluorescence intensity changes of the two fluorescent proteins were detected using the Promega Dual Luciferase Reporter Assay system kit (Cat. E1910) (Madison, WI, USA), and the activity of the transcription factors in the plant cells was analysed by comparing the changes in the fluorescence intensity ratio.

Detection of Transcription Factor Regulation in Protoplasts
The promoter of the paclitaxel biosynthetic pathway gene was inserted into the coding region of the LUC fluorescent protein in the 0800 plasmid to construct a plasmid, which contained the reporter gene with promoter-controlled LUC fluorescent protein and the continuous expression of REN fluorescent protein as internal reference. At the same time, the transcription factor to be detected is inserted into another plasmid 62sk, and the transcription factor can be continuously expressed in the plant by 35S regulation. The two plasmids were co-transformed into protoplasts by PEG transformation. The fluorescence intensity of the two fluorescent proteins was simultaneously detected using the Promega Dual Luciferase Reporter Assay system (Cat. E1910) kit, and the regulation of the transcription factor on the paclitaxel biosynthesis pathway promoter was analyzed.

Conclusions
In this study, based on the jasmonic acid signalling pathway of the Arabidopsis model, the jasmonic acid signalling pathway of T. media was isolated and identified. The interaction between related proteins was analysed, and the molecular mechanism of the jasmonic acid signalling pathway in Taxus was constructed. Further analysis of the activation or inhibition of the promoter of paclitaxel biosynthesis pathway gene by transcription factors such as MYC2 confirmed the regulation of MYC2 and other transcription factors on the biosynthesis pathway of paclitaxel. Based on this design, the transcriptional activity of MYC2 and MYC3 was not affected by the JAZ inhibition function, which laid a foundation for the construction of a high-yield paclitaxel variant Taxus cell line.

Figure 5.
Location map of cis-elements in the paclitaxel biosynthetic pathway gene promoter sequence. The white rectangles show the promoter region of the corresponding gene in proportion, the dark red arrow indicates the start of the coding frame of the gene; the red triangle indicates the G-box, and the yellow triangle indicates the high affinity E-box or the variant G-Box. The green triangle indicates low affinity or no effect E-box. Figure 6. Detection of protoplasts MYC-like transcription factors regulate the activity of the paclitaxel biosynthetic pathway gene promoter. The plasmid overexpressing the MYC transcription factor was used as the effector, the fluorescent protein LUC regulated by the paclitaxel biosynthesis pathway gene promoter was used as the reporter gene, and the overexpressed red fluorescent protein REN was used as the internal reference to be transferred into the Arabidopsis protoplasts. The ratio between the fluorescent signals. Empty plasmid was used as a negative control.  The plasmid overexpressing the MYC transcription factor was used as the effector, the fluorescent protein LUC regulated by the paclitaxel biosynthesis pathway gene promoter was used as the reporter gene, and the overexpressed red fluorescent protein REN was used as the internal reference to be transferred into the Arabidopsis protoplasts. The ratio between the fluorescent signals. Empty plasmid was used as a negative control. Figure A6. Detection of protoplasts MYC-like transcription factors regulate the activity of the paclitaxel biosynthetic pathway gene promoter. The plasmid overexpressing the MYC transcription factor was used as the effector, the fluorescent protein LUC regulated by the paclitaxel biosynthesis pathway gene promoter was used as the reporter gene, and the overexpressed red fluorescent protein REN was used as the internal reference to be transferred into the Arabidopsis protoplasts. The ratio between the fluorescent signals. Empty plasmid was used as a negative control.