TcbHLH14 a Jasmonate Associated MYC2-like Transcription Factor Positively Regulates Pyrethrin Biosynthesis in Tanacetum cinerariifolium

Natural pyrethrins have high application value, and are widely used as a green pesticide in crop pest prevention and control. Pyrethrins are mainly extracted from the flower heads of Tanacetum cinerariifolium; however, the natural content is low. Therefore, it is essential to understand the regulatory mechanisms underlying the synthesis of pyrethrins through identification of key transcription factors. We identified a gene encoding a MYC2-like transcription factor named TcbHLH14 from T. cinerariifolium transcriptome, which is induced by methyl jasmonate. In the present study, we evaluated the regulatory effects and mechanisms of TcbHLH14 using expression analysis, a yeast one-hybrid assay, electrophoretic mobility shift assay, and overexpression/virus-induced gene silencing experiments. We found that TcbHLH14 can directly bind to the cis-elements of the pyrethrins synthesis genes TcAOC and TcGLIP to activate their expression. The transient overexpression of TcbHLH14 enhanced expression of the TcAOC and TcGLIP genes. Conversely, transient silencing of TcbHLH14 downregulated the expression of TcAOC and TcGLIP and reduced the content of pyrethrins. In summary, these results indicate that the potential application of TcbHLH14 in improving the germplasm resources and provide a new insight into the regulatory network of pyrethrins biosynthesis of T. cinerariifolium to further inform the development of engineering strategies for increasing pyrethrins contents.


Introduction
Tanacetum cinerariifolium is a perennial herb of the genus Tanacetum belonging to the Asteraceae family, which is an economically important species and widely used for the extraction of essential oils, mainly comprising pyrethrins that have fast-acting, broadspectrum insecticidal activities with high efficiency. Most importantly, pyrethrins have a short half-life (approximately 2 h) and thus do not leave residue or enrich the environment, offering a safe alternative to many commonly used synthetic pesticides that have potential to cause harm to human health and the ecological environment [1,2]. Thus, pyrethrins have emerged as the most promising botanical insecticides. However, the naturally available amount of pyrethrins in the flower heads is too low with the content of 0.36% to 1.35% in T. cinerariifolium populations [3], thereby limiting the planting and application. Therefore, it is necessary to study the biosynthesis and regulation of pyrethrins to improve the productivity.
The bHLH (basic Helix-Loop-Helix) transcription factor family represents one of the most extensive families of transcription factors in plants. This family comprises two distinct regions with different functions: the basic region and the HLH region. The basic region, which specifically binds to the E-box (CANNTG) or G-box (CACGTG) of the DNA sequence [15]. Conversely, the HLH region has the capability to form heterodimers or homodimers with other proteins [16]. bHLH transcription factors play their roles mainly through the JA signal pathway in plants, and inhibit/activate the expression of JA-responsive genes [17,18], in turn regulating plant growth and development, the stress response, and synthesis of secondary metabolites. CrMYC2 can upregulate the expression of octadecanoid-response Catharanthus AP2 proteins (ORCA) to increase the biosynthesis of terpenoid indole alkaloids (TIAs) in Catharanthus roseus [19]. The heterologous expression of AtMYC2 can enhance the transcription of genes involved in the MEP pathway, thereby accumulating more Abietane diterpenes in Salvia striata [20]. AabHLH112 promotes artemisinin biosynthesis in Artemisia annua [21], and AabHLH1 responds to ABA and activates the expression of artemisinin synthesis-related genes, thereby increasing artemisinin production [22].
We screened TcbHLH14 from T. cinerariifolium as a TF acting in response to the methyl jasmonate (MeJA) signal, and previous studies also found that MeJA can up-regulate the expression of genes related to pyrethrins synthesis to ultimately increase the pyrethrins content [12]. Together, these findings suggest that TcbHLH14 may be involved in the regulation of pyrethrins biosynthesis mediated by MeJA. To further elucidate the mechanism of pyrethrins transcriptional regulation, in the present study, we explored the regulatory effect of the TF TcbHLH14 induced by MeJA on the key genes of pyrethrins synthesis. These results have important theoretical and practical significance for developing an appropriate genetic engineering strategy to enhance the biosynthesis of pyrethrins and cultivate varieties with high pyrethrins contents.

Domain and Homology Analysis of the TcbHLH14 Gene
The bHLH amino acid sequence of A. thaliana was downloaded from the NCBI database, and the phylogenetic tree of TcbHLH14 pyrethrum and Arabidopsis bHLH TFs was constructed ( Figure 1A). The results showed that Arabidopsis AtbHLH14 has the closest genetic relationship with pyrethrum bHLH (GEW34613.1) ( Figure 1B), indicating that T. cinerariifolium GEW34613.1 may have a function similar to AtbHLH14. Therefore, the ORF sequence of GEW34613.1 was cloned and named TcbHLH14. independent evolutionary branches ( Figure 1C): grape, Asteraceae, and other (Cruciferae, Leguminosae, Labiatae, Apocynaceae, Solanaceae) branches, marked with green, blue, and pink background frames, respectively. T. cinerariifolium, Chrysanthemum × morifolium, and A. annua in the Asteraceae branch clustered in the same subfamily, which indicated that they had high homology, revealing that TcbHLH14 protein and its homologs are closely related within the Asteraceae family. (C) Phylogenetic analysis of TcbHLH14 and bHLH proteins of other plants; TcbHLH14 is marked with the red triangle; the green box indicates Vitaceae, the blue box represents Asteraceae, and the pink box includes Brassicaceae Burnett, Fabaceae, Apocynaceae, and Solanaceae. The tree was constructed by the neighbor-joining method with 1000 bootstrap replicates.
The bHLH proteins of other species with high similarity to the pyrethrum TF TcbHLH14 were searched in the NCBI/Blast database (https://www.ncbi.nlm.nih.gov/ Blast.cgi (accessed on 26 September 2022)). The selected homologous proteins were compared with TcbHLH14 proteins using MUSCLE and BoxShade software. The TcbHLH14 protein was compared with A. annua, PWA37438.1, and artichoke (Cynara cardunculus var). The bHLH proteins in Scolymus (Cynara cardunculus, XP_024996993.1), chrysanthemum (Chrysanthemum × morifolium, QVL29796.1), sunflower (Helianthus annuus, XP_022029453.1), and lettuce (Lactuca sativa, XP_023742395.1) were highly similar, reaching 88.54%, 65.12%, 88.78%, 71.43%, and 65.80% similarity, respectively. Compared to TcMYC2 and AtMYC2, TcbHLH14 displays a significant deletion in the bHLH MYC-N To further explore the evolutionary relationship of the T. cinerariifolium TF TcbHLH14, a molecular evolutionary tree was established based on the amino acid sequences of TcbHLH14 and bHLH_MYC proteins of other plant species using MEGAX software. Among the 18 selected plants, bHLH_MYC proteins mainly clustered into three independent evolutionary branches ( Figure 1C): grape, Asteraceae, and other (Cruciferae, Leguminosae, Labiatae, Apocynaceae, Solanaceae) branches, marked with green, blue, and pink background frames, respectively. T. cinerariifolium, Chrysanthemum × morifolium, and A. annua in the Asteraceae branch clustered in the same subfamily, which indicated that they had high homology, revealing that TcbHLH14 protein and its homologs are closely related within the Asteraceae family.

TcbHLH14 Gene Expression in Tissues
The transcriptional level of TcbHLH14 in the flowers, leaves, and floral organs of T. cinerariifolium at different developmental stages was detected by RT-qPCR. The expression level of TcbHLH14 in T. cinerariifolium decreased during the S1-S3 stages of flower development, reached the peak in the S4 stage, decreased gradually due to the gradual withering of flowers, and was barely expressed by the S6 stage ( Figure 2A). In addition, the expression level of TcbHLH14 was the highest in the L1 stage, followed by the L3 stage, and was the lowest in the L2 stage ( Figure 2B). After dissecting the flower head of the S2 stage into five parts (ray floret, disc floret, bract, receptacle, and pedicel) and the S4 stage into two parts (ray floret and disc floret), we found that the expression level of TcbHLH14 was the highest in the S2 bract, followed by the S4 closed-tube flower ( Figure 2C). domain. Additionally, the N-terminal region of TcbHLH14 exhibits numerous deleted protein residues in contrast to AtbHLH14 ( Figure S2).

TcbHLH14 Gene Expression in Tissues
The transcriptional level of TcbHLH14 in the flowers, leaves, and floral organs of T. cinerariifolium at different developmental stages was detected by RT-qPCR. The expression level of TcbHLH14 in T. cinerariifolium decreased during the S1-S3 stages of flower development, reached the peak in the S4 stage, decreased gradually due to the gradual withering of flowers, and was barely expressed by the S6 stage ( Figure 2A). In addition, the expression level of TcbHLH14 was the highest in the L1 stage, followed by the L3 stage, and was the lowest in the L2 stage ( Figure 2B). After dissecting the flower head of the S2 stage into five parts (ray floret, disc floret, bract, receptacle, and pedicel) and the S4 stage into two parts (ray floret and disc floret), we found that the expression level of TcbHLH14 was the highest in the S2 bract, followed by the S4 closed-tube flower ( Figure 2C).
MeJA is an indispensable signaling molecule in plants, which not only regulates secondary plant metabolism but also participates in plant defense responses to biotic and abiotic stresses. We found that the TcbHLH14 gene could be induced by MeJA at different time periods ( Figure 2D).  MeJA is an indispensable signaling molecule in plants, which not only regulates secondary plant metabolism but also participates in plant defense responses to biotic and abiotic stresses. We found that the TcbHLH14 gene could be induced by MeJA at different time periods ( Figure 2D).

Subcellular Localization of TcbHLH14
To identify the subcellular localization of TcbHLH14 protein, the TcbHLH14-GFP fusion protein, driven by the Super promoter ( Figure 3A), was transiently expressed in N. benthamiana and detected by laser confocal scanning microscopy. The green fluorescence signal detected in the TcbHLH14-GFP fusion protein overlapped with the red fluorescence signal in the nuclear localization plasmid under specific excitation light, indicating that TcbHLH14 is a nuclear-localized protein, which corresponds with exerting its function as a TF in regulating gene expression ( Figure 3B).

Subcellular Localization of TcbHLH14
To identify the subcellular localization of TcbHLH14 protein, the TcbHLH14-GFP fusion protein, driven by the Super promoter ( Figure 3A), was transiently expressed in N. benthamiana and detected by laser confocal scanning microscopy. The green fluorescence signal detected in the TcbHLH14-GFP fusion protein overlapped with the red fluorescence signal in the nuclear localization plasmid under specific excitation light, indicating that TcbHLH14 is a nuclear-localized protein, which corresponds with exerting its function as a TF in regulating gene expression ( Figure 3B).

TcbHLH14 Directly Binds to the Promoters of the TcAOC and TcGLIP Genes
The yeast one-hybrid (Y1H) assay, luciferase reporter assay, and EMSA were used to verify whether TcbHLH14 interacts with the promoters of the TcAOC and TcGLIP genes in the pyrethrins biosynthesis pathway. The promoter of each gene was cloned and inserted into the pHis2.1 vector to generate reporter vectors. TcbHLH14 was fused to the GAL4 activation domain to generate the effector construct pGADT7::TcbHLH14 ( Figure  4A). After co-transformation of pGADT7::TcbHLH14/pHis2.1::proTcAOC/proTcGLIP, the Y187 yeast strain showed normal growth on the three-deficiency culture plate SD/-Trp/-Leu/-His, indicating that there was an interaction between TcbHLH14 and the TcAOC and TcGLIP promoters ( Figure 4B). To verify the interaction of TcbHLH14 with the TcAOC and TcGLIP promoters in vivo, reporter vector (TcAOC/TcGLIPpro:LUC) and effector vector (35S: TcbHLH14) were transiently coexpressed in N. benthamiana leaves ( Figure 4C). The firefly luciferase activity was significantly higher than that of the control group after coinfection of pGreenII62-SK::TcbHLH14/pGreenII0800-LUC::pTcAOC/TcGLIP into N. benthamiana leaves, indicating that the TF TcbHLH14 could activate the transcription of TcAOC and TcGLIP gene promoters ( Figure 4D). In the EMSA, TcbHLH14-His was transferred into the E. coli Rosetta strain for induction and purification. A labeled probe containing the G/E-box motif was designed in the TcAOC/TcGLIP promoter regions and the EMSA results showed that TcbHLH14 protein could bind to the probe along with a significant blocking band. When a 50× competitive cold probe was added, the blocking binding band was close to zero, whereas when the 50× mutant cold probe was added, the blocking band was restored, indicating that TcbHLH14 could bind to the G/E-box element in the TcAOC/TcGLIP promoters ( Figure 4E). In summary, these assays collectively

TcbHLH14 Directly Binds to the Promoters of the TcAOC and TcGLIP Genes
The yeast one-hybrid (Y1H) assay, luciferase reporter assay, and EMSA were used to verify whether TcbHLH14 interacts with the promoters of the TcAOC and TcGLIP genes in the pyrethrins biosynthesis pathway. The promoter of each gene was cloned and inserted into the pHis2.1 vector to generate reporter vectors. TcbHLH14 was fused to the GAL4 activation domain to generate the effector construct pGADT7::TcbHLH14 ( Figure 4A). After co-transformation of pGADT7::TcbHLH14/pHis2.1::proTcAOC/proTcGLIP, the Y187 yeast strain showed normal growth on the three-deficiency culture plate SD/-Trp/-Leu/-His, indicating that there was an interaction between TcbHLH14 and the TcAOC and TcGLIP promoters ( Figure 4B). To verify the interaction of TcbHLH14 with the TcAOC and TcGLIP promoters in vivo, reporter vector (TcAOC/TcGLIPpro:LUC) and effector vector (35S: TcbHLH14) were transiently coexpressed in N. benthamiana leaves ( Figure 4C). The firefly luciferase activity was significantly higher than that of the control group after co-infection of pGreenII62-SK::TcbHLH14/pGreenII0800-LUC::pTcAOC/TcGLIP into N. benthamiana leaves, indicating that the TF TcbHLH14 could activate the transcription of TcAOC and TcGLIP gene promoters ( Figure 4D). In the EMSA, TcbHLH14-His was transferred into the E. coli Rosetta strain for induction and purification. A labeled probe containing the G/E-box motif was designed in the TcAOC/TcGLIP promoter regions and the EMSA results showed that TcbHLH14 protein could bind to the probe along with a significant blocking band. When a 50× competitive cold probe was added, the blocking binding band was close to zero, whereas when the 50× mutant cold probe was added, the blocking band was restored, indicating that TcbHLH14 could bind to the G/E-box element in the TcAOC/TcGLIP promoters ( Figure 4E). In summary, these assays collectively demonstrated that TcbHLH14 can directly bind to the TcAOC and TcGLIP promoters to activate these pyrethrin biosynthesis genes.

Transient Overexpression and VIGS of TcbHLH14 in T. cinerariifolium
To further study the effect of TcbHLH14 on pyrethrins biosynthesis, pGreenII62SK::TcbHLH14 was transiently transformed into T. cinerariifolium leaves for overexpression. The expression levels of pyrethrins biosynthesis pathway genes TcGLIP, TcAOC, TcCHS, TcALDH1, and the target gene TcbHLH14 were detected by RT-qPCR. Compared with that of the control, the expression of TcbHLH14 was significantly up-

Transient Overexpression and VIGS of TcbHLH14 in T. cinerariifolium
To further study the effect of TcbHLH14 on pyrethrins biosynthesis, pGreenII62SK::Tcb HLH14 was transiently transformed into T. cinerariifolium leaves for overexpression. The expression levels of pyrethrins biosynthesis pathway genes TcGLIP, TcAOC, TcCHS, TcALDH1, and the target gene TcbHLH14 were detected by RT-qPCR. Compared with that of the control, the expression of TcbHLH14 was significantly up-regulated on the second, fourth, and sixth day after transfection. As a whole, the expression profile showed a changing trend of rising first and then falling. The expression levels of TcGLIP, TcAOC, TcCHS, and TcALDH1 were significantly increased on the fourth day, and the expression of TcAOC, TcCHS, and TcALDH1 was up-regulated by more than 2.96 times. The TcGLIP expression level increased by 1.57 times ( Figure 5A). However, no significant difference (p > 0.05) was observed in the total pyrethrins content between the control and treatment groups ( Figure 5B, Table S1). regulated on the second, fourth, and sixth day after transfection. As a whole, the expression profile showed a changing trend of rising first and then falling. The expression levels of TcGLIP, TcAOC, TcCHS, and TcALDH1 were significantly increased on the fourth day, and the expression of TcAOC, TcCHS, and TcALDH1 was up-regulated by more than 2.96 times. The TcGLIP expression level increased by 1.57 times ( Figure 5A). However, no significant difference (p > 0.05) was observed in the total pyrethrins content between the control and treatment groups ( Figure 5B, Table S1).
To further verify the role of TcbHLH14 in pyrethrins biosynthesis, pTRV2::TcbHLH14 and pTRV1 virus-mediated silencing vectors were transiently transformed into pyrethrum leaves, and the expression levels of the above five genes were determined. The results of RT-qPCR showed that the expression level of TcbHLH14 in the VIGS group was significantly lower than that of the control. The expression of TcGLIP and TcCHS was down-regulated by more than 2 times, while the expression of TcAOC and TcALDH1 was slightly up-regulated after TcbHLH14 silencing, but there was no significant change (Figure 5C). Moreover, the content of pyrethrins in the leaves decreased after TcbHLH14 silencing ( Figure 5D, Table S2).  To further verify the role of TcbHLH14 in pyrethrins biosynthesis, pTRV2::TcbHLH14 and pTRV1 virus-mediated silencing vectors were transiently transformed into pyrethrum leaves, and the expression levels of the above five genes were determined. The results of RT-qPCR showed that the expression level of TcbHLH14 in the VIGS group was significantly lower than that of the control. The expression of TcGLIP and TcCHS was down-regulated by more than 2 times, while the expression of TcAOC and TcALDH1 was slightly up-regulated after TcbHLH14 silencing, but there was no significant change ( Figure 5C). Moreover, the content of pyrethrins in the leaves decreased after TcbHLH14 silencing ( Figure 5D, Table S2).

Discussion
The plant hormone MeJA is a ubiquitous and conservative plant secondary metabolic elicitor [23]. Previous studies have found that MeJA can effectively promote pyrethrins accumulation [24], but its regulatory mechanism is unclear. To explore the effect of JA on pyrethrins synthesis, transcriptome analysis of pyrethrum leaves treated with MeJA was performed, which showed that MeJA could trigger transcriptional re-editing, resulting in extensive metabolic changes in pyrethrum. Moreover, the expression of MYC, the core factor of the JA signal transduction pathway, was significantly up-regulated after MeJA treatment. In the present study, a MYC2-like TF named TcbHLH14 was isolated and cloned from T. cinerariifolium. MYC2 was reported to regulate the synthesis of plant secondary metabolites by activating the expression of its target genes. In Aquilaria sinensis, AsMYC2 regulates the biosynthesis of sesquiterpene by controlling the expression of the ASS1 gene [25]. In response to JA, MdMYC2 could significantly up-regulate the transcription level of the anthocyanin biosynthesis gene and increase the content of anthocyanin in apples (Malus domestica) [26]. In view of the important role of the TF MYC2 in the synthesis of secondary metabolites induced by JA, it is speculated that pyrethrum TcbHLH14 may also induce the transcription of key genes in the pyrethrins metabolic pathway to consequently regulate pyrethrins synthesis.
We found that TcbHLH14 exerts this regulatory role by binding to the E-box or G-box region in the target gene promoter. TFs increase or weaken the target gene transcription mainly by specifically recognizing the binding motif in the promoter. E-box (CANNTG) or G-box (CACGTG) is the key element of transcriptional regulation of bHLH proteins [15]. The bHLH_MYC transcription factor TwMYC2a/b of Tripterygium wilfordii regulates triptolide synthesis by binding to E-box (CACATG) and T/G-box (CACGTT) motifs in the TwTPS27a/b promoter [27]. S. miltiorrhiza SmbHLH3 can specifically recognize and bind the G-box of the promoters of TAT, HPPR, KSL1, and CYP76AH1, which are key genes in the synthesis of phenolic acid and tanshinone [28].
The promoter sequences of genes in the pyrethrins biosynthesis pathway were isolated and cloned in our laboratory. Further analysis revealed that the TcAOC and TcGLIP promoters contained six and eight bHLH TF-binding elements, respectively. Through EMSA and Y1H assays, we further found that TcbHLH14 could bind directly to the G-box and E-box motif in the TcAOC and TcGLIP promoter, respectively. The results of the luciferase reporter assay further showed that TcbHLH14 could activate the expression of the TcAOC and TcGLIP promoters. These results suggest that TcAOC and TcGLIP are the target genes of the TF TcbHLH14. Transient overexpression of TcbHLH14 could effectively up-regulate the expression of the TcAOC and TcGLIP genes; however, the content of pyrethrins was not significantly up-regulated. The reason may be that TcbHLH14 has a fine regulation function similar to that of AtMYC2 [29], which exerts negative regulation of its own expression. Unlike our previous research on TcMYC2, the results of Y1H and Dual-LUC assays indicate that TcbHLH14 does not exhibit binding activity to the TcCHS promoter. Furthermore, while transient overexpression of TcMYC2 in leaves up-regulated TcCHS expression by more than 20-fold [30], transient overexpression of TcbHLH14 results in only a marginal increase in TcCHS expression. This outcome could potentially lead to a reduction in DMAPP flux towards chrysanthemol, resulting in a lack of significant increase in the final accumulation of pyrethrins. In addition, some studies have found that members of the MYC TF family can co-express and regulate the synthesis of secondary metabolites, as is the case for rice OsMYC2 and the OsMYC2-like proteins OsMYL1 and OsMYL2 (both OsMYC2-like proteins), which promote the synthesis of sakuranetin [31]. Thus, lack of cooperation between TcbHLH14 and other MYC proteins may also be the reason for the lack of an effective increase in pyrethrins content under overexpression, although this hypothesis requires further exploration. Transient silencing of TcbHLH14 could significantly reduce the expression of the TcGLIP gene and the content of pyrethrins, whereas the expression level of the TcAOC gene was slightly enhanced, which may be due to the competitive binding of other proteins to TcbHLH14 and activation of TcAOC gene expression. In addition, the expression of the pyrethrins biosynthesis pathway genes TcALDH1 and TcCHS is also regulated by TcbHLH14. In summary, our results demonstrate that TcbHLH14 is a positive regulator to promote pyrethrins biosynthesis.
The synthesis of plant secondary metabolites is closely related to the expression level of a series of key structural enzyme genes, and TFs can accurately regulate structural enzyme genes. In a previous study, Dianthus caryophyllus transgenic plants were established with overexpression of the linalool synthase gene of Clarkia breweri by transgenic technology. Gas chromatography-mass spectrometry analysis of the petal extract showed no significant difference in the linalool content between overexpression and wild-type plants [32]. This suggests that the transfer of a structural gene alone does not necessarily result in the significant accumulation of bioactive products. TFs can often activate the co-expression of multiple key enzyme-encoding genes in a metabolic pathway, which is more efficient than multiple or co-transformation of one or more synthase genes [33]. For example, transferring the TF AaORA into A. annua up-regulated the expression of the structural enzyme genes ADS, CYP71AV1, and DBR2, and significantly increased the contents of artemisinin and dihydroartemisinic acid [34]. S. miltiorrhiza SmMYB52 directly binds to and activates the promoters of genes encoding four enzymes in the salvianolic acid B biosynthesis pathway (SmTAT1, Sm4CL9, SmC4H1, SmHPPR1) to promote salvianolic acid B accumulation [35]. In citrus, CsMADS5 positively regulates carotenoid synthesis by promoting the transcription of the PSY, PDS, and LCYb1 genes [36]. Similarly, we found that TcbHLH14 can specifically bind to the pyrethrins synthesis pathway genes TcGLIP and TcAOC to consequently up-regulate their expression. Therefore, using TFs to regulate the biosynthesis of secondary plant metabolites has great potential for producing valuable bioactive compounds.

Sequence and Phylogenetic Analysis
The bHLH amino acid sequence of Arabidopsis thaliana was downloaded from the National Center for Biotechnology Information (NCBI) database, and the phylogenetic tree of TcbHLH14 pyrethrum and Arabidopsis bHLH TFs was constructed by MEGA X v10.1.8 [37], and visualized by the online software iTOL v5 [38].

Plant Materials
The tissue culture seedlings of T. cinerariifolium 'W99 used in this study are preserved in the plant culture room. The seedlings were cultured under a temperature of 25 ± 2 • C with light for 16 h and dark for 8 h. Late the seedlings were planted in the field of flower base of Huazhong Agricultural University. The flowers at the S1-S6 stages [39] and the leaves at the L1-L3 stages [40] were harvested and stored at −80 • C until analysis. The genotype of 'W99 seedlings grown in culture flask for 1 month were sprayed with 5 mL of 2 mM MeJA (Macklin Biochemical Co., Ltd., Shanghai, China) dissolved in 0.8% ethanol, and the leaves were collected in triplicate at 0 (control), 2, 4, 6, 8, 12, and 24 h.

RNA Extraction and Reverse Transcription-Quantitative Polymerase Chain Reaction (RT-qPCR)
Total RNA was extracted from T. cinerariifolium using TRIzol Reagent (Invitrogen, CA, USA), and reverse-transcribed to synthesize cDNA using TransScript ® One-Step gDNA Removal kit and cDNA Synthesis SuperMix kit (TransGen Biotech, Beijing, China). The TcbHLH14 gene sequence was screened from the MeJA-treated transcriptome determined in our laboratory, and the specific primers for amplifying the open reading frame (ORF) fragment were designed by Primer Premier 5 v1.0 software.
The cDNA was used as the template for qPCR on the Roche LightCycler ® 96 Real-time PCR System (Roche, Basel, Switzerland), using SYBR Premix Ex Taq II kits (Takara, Kusatsu, Japan) with three biological and two technical replicates. The qPCR procedure was similar to that reported previously [41]. TcGADPH was used as the internal reference, and the sequences of the specific primers are shown in Supplementary Table S3 The 2 −∆∆CT method was used to analyze the relative gene expression [42].

Subcellular Localization Analysis
After digesting the pSuper1300-green fluorescent protein (GFP) vector with Swa1 and Kpn1 restriction enzymes, the PCR fragment containing the TcbHLH14 coding region (without the stop codon) was ligated to the linearized pSuper1300-GFP vector by ClonExpress ® II One Step Cloning Kit (Vazyme Biotech, Nanjing, China). The recombinant plasmid was introduced into Agrobacterium tumefaciens strain GV3101 and then co-transformed into tobacco (Nicotiana benthamiana) leaves with a red fluorescent protein (RFP)-tagged nuclear marker. After 72 h of weak light exposure, fluorescence signals were observed with a Leica TCS-SP8 fluorescence microscope (Leica, Wetzlar, Germany).

Yeast One-Hybrid (Y1H) Assay
The TcbHLH14 sequence was integrated in the pGADT7 vector, which was digested with NdeI and XhoI restriction enzymes to obtain the pGADT7:bHLH14 recombinant plasmid. The pH2.1:TcGLIP and pH2.1:TcAOC vectors were obtained as described previously [30]. The effector and reporter plasmids were co-transformed into the yeast strain Y187 using the Super Yeast Transformation kit (Coolaber, Beijing, China). The transformed yeast cells were coated on a deficient medium supplemented with 3-amino-1,2,4-triazol (3-AT) at various concentrations to determine the appropriate concentration of 3-AT, the transformed yeast cells were selected on DDO (SD/-Leu/-Trp) medium, and the Y1H interaction assay was carried out on 30 • C TDO (SD/-Leu/-Trp/-His) medium for 3 days.

Dual-Luciferase (LUC) Transient Expression Assays
To obtain the recombinant vector, the ORF sequence of TcbHLH14 was integrated in the HindIII linearized pGreenII62-SK vector by the ClonExpress II One Step Cloning Kit. Together with the previously constructed pGreenII0800-LUC:AOC and GLIP vectors [30], the reporter vector was transformed into A. tumefaciens GV3101 (containing the pSoup19 vector) and then transformed into N. benthamiana leaves. The negative control vector was transformed into the corresponding part of the same leaf. Leaf materials were collected after 48 h under low light. The luminescence was detected under the NightSHADE LB 985 imaging system (Berthold, Bad Wildbad, Germany) by the Dual-Luciferase reporter assay system kit (Promega, Madison, WI, USA).

Electrophoretic Mobility Shift Assay (EMSA)
The ORF of TcbHLH14 without the stop codon was integrated in BamHI and XhoI -linearized pET28a vector, which was transformed in the Escherichia coli Rosetta (DE3) strain. The recombinant protein was induced by 0.5 mM isopropyl β-d-1-thiogalactopyranoside, incubated at 18 • C for 16 h, and purified on a Ni-NTA column. The cis-elements (E-box motif, G-box motif) in the TcAOC and TcGLIP promoter fragments were labeled at the 5 end with the fluorophore 6-FAM, and 50-fold of unlabeled wild-type and cis-element mutant DNA fragments were used as competitors. EMSA detection was then performed as described previously [43].

Virus-Induced Gene Silencing (VIGS)
After linearizing the pTRV2 vector with BamHI enzyme, The 276-bp specific fragment from TcbHLH14 was designed by siRNA Finder software [44], and integrated into the pTRV2 vector by homologous recombination, and the obtained plasmids pTRV2:TcbHLH14 and pTRV1 were co-transformed into the chemically active A. tumefaciens GV3101 strain.
The VIGS experiment was carried out in T. cinerariifolium as reported previously [45]. Cotransformation of the empty vectors pTRV1 and pTRV2 was used as a control. After 14 days, the samples were used for RT-qPCR and high-performance liquid chromatography (HPLC) analyses.

Transient Overexpression of TcbHLH14
The expression vector pGreenII62-SK was driven by the CaMV35S promoter and the constructed pGreenII62-SK:TcbHLH14 vector was transformed into A. tumefaciens GV3101 (containing the pSoup19 vector). Transient transformation of T. cinerariifolium leaves was performed as described previously [46]. After the RT-qPCR assay, 4-day-old leaves were taken for HPLC analysis.

Pyrethrin Content Measurement by HPLC
The pyrethrins were extracted as described previously [47] and filtered with a 0.22-µm mesh filter, then determined using the Waters HPLC system (Waters Corporation, Milford, MA, USA) equipped with a photodiode array detector using three biological replicates as described previously [48]. pyrethrum extract (Sigma-Aldrich, St. Louis, MO, USA) was used as the standard.

Statistical Analysis
All experiments were performed with at least three biological replicates. One-way ANOVA and Student's t-test (or Tukey post hoc test) were performed using SPSS v.20 software. p-values < 0.05 were considered statistically significant.