SmMYB4 Is a R2R3-MYB Transcriptional Repressor Regulating the Biosynthesis of Phenolic Acids and Tanshinones in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge is one of the most famous traditional Chinese medicinal plants. The two most important classes of pharmaceutically relevant compounds in S. miltiorrhiza are phenolic acids and tanshinones. The MYB family of transcription factors may efficiently regulate the secondary metabolism in plants. In this study, a subgroup 4 R2R3MYB transcription factor gene, SmMYB4, was isolated from S. miltiorrhiza and functionally characterized using overexpression and a RNAi-mediated silencing. We achieved a total of six overexpressions and eight RNAi transgenic lines from the Agrobacterium leaf disc method. The content of the total phenolics, rosmarinic acid, and salvianolic acid B markedly decreased in the SmMYB4-overexpressing lines but increased in the SmMYB4-RNAi lines. The content of the total tanshinones, cryptotanshinone, and tanshinone IIA decreased in the SmMYB4-overexpressing transgenic lines but increased in the SmMYB4-RNAi lines. A gene expression analysis demonstrated that SmMYB4 negatively regulated the transcription of the critical enzyme genes involved in the phenolic acid and tanshinone biosynthesis. The genetic control of this transcriptional repressor may be used to improve the content of these bioactive compounds in the cultivated S. miltiorrhiza.


Introduction
Salvia miltiorrhiza Bunge, called 'Dan-Shen' in Chinese, is one of the most widely used traditional herbal medicines for various cardiovascular and cerebrovascular diseases. Water-soluble phenolic acids and lipid-soluble tanshinones are the two major classes of pharmaceutical ingredients in S. miltiorrhiza [1]. With an increasing market demand for S. miltiorrhiza, the wild resources are rapidly decreasing and becoming extremely scarce. However, the quality of cultivated S. miltiorrhiza is deteriorating. Therefore, improving the content of the bioactive components in cultivated S. miltiorrhiza using modern plant genetic engineering and metabolic engineering technology is important. The biosynthesis, accumulation, and regulation of these active ingredients have become intensively attractive research areas.
The transcription factors (TFs) can bind the cis-elements of a target gene's promoter region to activate or inhibit its transcription. They offer an effective strategy for enhancing the production of pharmaceutically active metabolites in medicinal plants [2]. Following the publication of the S. miltiorrhiza genome [3], many transcription factor gene families, proportion of vermiculite and nutrient soil and germinated at 25 ± 2 • C, using a cycle of light (16 h at 8000 lx) and dark (8 h). The two-month-old seedlings were frozen in liquid nitrogen for the subsequent total RNA extraction. The seeds were sterilized and cultured on a MS (Murashige and Skoog) medium to obtain sterile plantlets, as described by Yan and Wang [49]. The leaves from these plantlets were used for the genetic transformation experiments.

Subcellular Localization of SmMYB4
The cDNA sequence (with no termination codon) of SmMYB4 was obtained from a public database and then cloned and sequenced from the cDNA library of S. miltiorrhiza with the designed primers. The sequence was inserted into pHBT-GFP-NOS vector at KpnI and BamHI restricted sites to form SmMYB4-GFP. The plasmids were transiently transformed into onion epidermal cells using the PDS-1000 Gene Gun (American Bole BIO-RAD, Santa Clara, CA, USA). The onion epidermis cells were isolated and observed under a Leica DM6000B microscope (Leica, Leica Microsystems Co., Ltd., Wetzlar, Germany) at 475 nm.

Construction of the SmMYB4-OE and SmMYB4-RNAi Expression Vectors
According to the manufacturer's protocol, the total RNA was extracted using E.Z.N.A. OMEGA TM Plant RNA Kit and converted into cDNA using PrimeScript TM RT Master Mix Kit (Takara, Bao Biological Engineering (Dalian) Co., Ltd., Dalian, China).
To construct the overexpression vector, the entire open reading frame (ORF) of SmMYB4 was amplified with the primers OESmMYB4-BglII and OESmMYB4-BstPI (Supplementary Material Table S1). The PCR (polymerase chain reaction) products were ligated into the pMD19-T vector (Takara) and transformed into Escherichia coli DH5α cells. The recombinant plasmids were digested with BglII/BstPI and then ligated into the pCAMBIA1302 vector to generate the overexpression plasmid OESmMYB4-1302. To construct the RNAi plasmid, a specific 246 bp fragment of SmMYB4 was selected as the RNAi target fragment. The forward and reverse RNAi target fragments were amplified with the primers ISmMYB4-HindIII/ISmMYB4-BamHI and ISmMYB4-KpnI/ISmMYB4-XhoI, respectively, and inserted into the pKANNIBAL vector. An interfering box was cloned into pART27 to generate the SmMYB4-RNAi vector named ISmMYB4-pART27. To verify whether the overexpression and RNAi vectors were successfully constructed, OESmMYB4-1302 and ISmMYB4-pART27 were digested with BglII /BstPI and XhoI/NotI, respectively.

Transformation and Detection of S. Miltiorrhiza
The OESmMYB4-1302 and ISmMYB4-pART27 vectors were separately transformed into the Agrobacterium tumefaciens strain GV3101 using a freeze-thaw method [50]. The A. tumefaciens GV3101 harboring the recombinant vectors of OESmMYB4-1302 or ISmMYB4-pART27 were used to generate the transgenic plants according to the protocol established in our laboratory [49].
To confirm whether the target sequences had been integrated into the genome of S. miltiorrhiza, the genomic DNA was extracted (following OMEGA DNA extraction kit instructions) from the fresh leaves of one-month-old transgenic plantlets and nontransformed (WT) plants. The primers of 35S promoter and Kan resistance gene were used for identifying the overexpression and the RNAi transgenic plantlets, respectively (Supplementary Material Table S2). The PCR products were detected using 1% agarose gel electrophoresis.
The total RNA was isolated from the young leaves and converted into cDNA using the above method. The reverse transcripts were utilized as templates for analyzing the expression of SmMYB4 in the transgenic lines and the WT lines, using real-time quantitative PCR (qPCR). The primers Q-SmMYB4-F and Q-SmMYB4-R were synthesized by Beijing Aoke peak Biotechnology Co., Ltd., Beijing, China (Supplementary Material Table S2). SmUbiquitin (JF760206.1) served as an internal control to normalize the expression levels. The qPCR reactions were performed in triplicate under the same conditions, as described previously [44]. The relative gene expression was calculated using the comparative C T method.

Determination of the Total Phenolic Acids and Flavonoids
The transgenic and wild S. miltiorrhiza were sub-cultured for 90 d. Then, the roots were dried in the shade. The sample preparation followed the methods described by Zhang [45]. The total phenolic acid and flavonoid contents were measured using the Folin-Ciocalteu and NaNO 2 -AlCl 3 methods, respectively [51,52]. The gallic acid and epicatechin (EC) were used as a standard for determining the total phenolic and flavonoid content, respectively.

Evaluation of the Salvianolic Acid and Tanshinone Contents
The sample preparation followed the previously described method [45]. The contents of the salvianolic acid and tanshinones were determined using a HPLC (highperformance liquid chromatography) system (Agilent Zorbax, Agilent, Santa Clara, CA, USA). The chromatographic separation was performed using a C18 column (Agilent Eclipse XDB-C18-USP L1, 4.6 mm × 250 mm, 5 µm particle size) at 30 • C and a sample injection volume of 20 µL. The detection at 280 nm, using a flow rate of 1.0 mL/min over a gradient of methanol containing 0.4% acetic acid (buffer A), triple distilled water containing 0.4% acetic acid (buffer B), and acetonitrile (Fisher Scientific) containing 0.4% acetic acid (buffer C). The gradient conditions are shown in Supplementary Material Table S3. The standards of rosmarinic acid (RA), salvianolic acid B (SAB), cryptotanshinone (CT), and tanshinone II A (T-II A) were purchased from the National Institute for the Control of Pharmaceutical and Biological Products (Beijing, China).

Analysis of the Target Genes Regulated by SmMYB4
To elucidate the targets of the transcriptional regulation by SmMYB4, the 2.5 kb upstream sequences of the key enzyme genes were extracted from their corresponding scaffolds. The published genome of S. miltiorrhiza was obtained from the National Data Center of Traditional Chinese Medicine of China (https://ngdc.cncb.ac.cn/search/?dbId= gwh&q=+PRJCA003150/accessed accessed on 18 August 2020). The amino acid sequence of SmMYB4 was used as the input to the online profile inference tool (http://jaspar.genereg. net/ accessed on 18 August 2020) to search for its binding sequences. The matrix with the smallest E-value was used to scan the promoter sequences of the candidate genes and to ascertain whether they were targets of SmMYB4 (the threshold was set to the default value of 80%).

Characterization Analysis of SmMYB4
The isolated full-length SmMYB4 contained an ORF comprising 696 bp, encoding a polypeptide of 231 amino acids with a molecular weight of 26 kDa and a pI of 8.865. SmMYB4 is predicted to localize in the nucleus (ProtComp 9.0 server, http://linux1 .softberry.com/ accessed on 3 May 2021), and is demonstrated by the further subcellular localization experiment (Supplementary Material Figure S1). The phylogenetic analysis (MEGA6.0) indicated that all R2R3-MYB repressors from various plants are divided into two clades ( Figure 1). SmMYB4 is grouped with AtMYB4, AtMYB32, NtMYB2, NtMYB3, and VvMYB4a, belonging to subgroup 4 of the R2R3-MYB transcriptional factors. SmMYB39, which was previously functionally characterized as a repressor of phenolic acid biosynthesis, belongs to the same clade as SmMYB4.
its binding sequences. The matrix with the smallest E-value was used to scan the promoter sequences of the candidate genes and to ascertain whether they were targets of SmMYB4 (the threshold was set to the default value of 80%).

Characterization Analysis of SmMYB4
The isolated full-length SmMYB4 contained an ORF comprising 696 bp, encoding a polypeptide of 231 amino acids with a molecular weight of 26 kDa and a pI of 8.865. SmMYB4 is predicted to localize in the nucleus (ProtComp 9.0 server, http://linux1.softberry.com/ accessed on 3 May 2021), and is demonstrated by the further subcellular localization experiment (Supplementary Material Figure S1). The phylogenetic analysis (MEGA6.0) indicated that all R2R3-MYB repressors from various plants are divided into two clades ( Figure 1). SmMYB4 is grouped with AtMYB4, AtMYB32, NtMYB2, NtMYB3, and VvMYB4a, belonging to subgroup 4 of the R2R3-MYB transcriptional factors. SmMYB39, which was previously functionally characterized as a repressor of phenolic acid biosynthesis, belongs to the same clade as SmMYB4. The protein sequence alignment of SmMYB4 and the numerous R2R3-MYB repressors indicated that SmMYB4 consists of conserved R2 and R3 DNA binding domains at the N-terminus. A bHLH binding motif ([D/E]Lx2[R/K]x3Lx6Lx3R) was found in the R3 DNA binding domain of SmMYB4. Additionally, the amino acid sequence of SmMYB4 has two typically conserved motifs at the C-terminus, C1 (KLIsrGIDPxT/SHRxI/L) and C2 (pdLNLD/ELxiG/S) (Figure 2), which were previously identified in the subgroup 4 R2R3-MYB transcriptional factors [26]. The protein sequence alignment of SmMYB4 and the numerous R2R3-MYB repressors indicated that SmMYB4 consists of conserved R2 and R3 DNA binding domains at the N-terminus. A bHLH binding motif ([D/E]Lx 2 [R/K]x 3 Lx 6 Lx 3 R) was found in the R3 DNA binding domain of SmMYB4. Additionally, the amino acid sequence of SmMYB4 has two typically conserved motifs at the C-terminus, C1 (KLIsrGIDPxT/SHRxI/L) and C2 (pdLNLD/ELxiG/S) (Figure 2), which were previously identified in the subgroup 4 R2R3-MYB transcriptional factors [26].

Molecular Identification of the Transgenic Plantlets
The fragments of about 700 bp were captured from the digestion products of the detected OESmMYB4-1302 vectors, and the expected fragments of approximately 1300 bp and 3400 bp were obtained from the cutting products of the detected ISmMYB4-pART27 vectors, which indicated that the overexpression and the RNAi vectors were successfully constructed (Supplementary Material Figure S2). Then, the constructed vectors were sequenced using their respective primers. The transgenic plants were generated to examine the function of SmMYB4. In this research, fourteen overexpression transgenic lines

Molecular Identification of the Transgenic Plantlets
The fragments of about 700 bp were captured from the digestion products of the detected OESmMYB4-1302 vectors, and the expected fragments of approximately 1300 bp and 3400 bp were obtained from the cutting products of the detected ISmMYB4-pART27 vectors, which indicated that the overexpression and the RNAi vectors were successfully

SmMYB4 Suppresses the Accumulation of Phenolic Acids and Flavonoids
To study the function of SmMYB4 in the secondary metabolites biosynthesis, the total phenolic acid and flavonoid contents of S. miltiorrhiza were detected ( Figure 4A,B). Compared with the WT, the total phenolic acid content was significantly lower in the SmMYB4 overexpression lines and markedly higher in the RNAi lines. The total phenolic acid content was minimal (7.

SmMYB4 Negatively Regulates the Salvianolic Acid and Tanshinone Biosynthesis
RA and SAB are the most abundant water-soluble active constituents in S. miltiorrhiza. It was found that their content had a significant change in the transgenic lines compared with the WT, which had 2.88 mg/g RA and 12.77 mg/g SAB (Figure 4). Compared to the WT, the RA and SAB contents were lower in the overexpression lines (O-4-3, O-4-7, O-4-9). Their respective levels in line O-4-9 notably decreased by 0.566-fold and 0.576-fold. Conversely, their content in the interference lines was significantly higher, being 1.41-and 1.80-fold of their respective contents in the WT ( Figure 4C,D).
The total tanshinone content in the overexpression lines was lower than that in the WT, especially in line O-4-9, which was 0.32-fold lower. In contrast, their content in all interference lines was significantly higher, up to 2.98-fold that of WT (Supplementary Material Figure S4). Compared with the WT, the CT content in the three overexpression lines was lower, the lowest being measured in line O-4-9 (0.04 mg/g), whereas the T-IIA content was not significantly changed. It is noteworthy that the content of CT and T-IIA was markedly higher in three interference lines ( Figure 4E,F), with line I-4-2 having the most CT (10.14-fold of WT) and line I-4-1 having the most T-IIA (17.2-fold of WT).

SmMYB4 Inhibits the Expression of Key Enzyme Genes
To uncover the target genes regulated by SmMYB4, the expression of twelve key enzyme genes involved in the salvianolic acid and tanshinone biosynthetic pathways was analyzed. The qPCR analysis demonstrated that silencing SmMYB4 markedly upregulated SmPAL1, SmC4H, Sm4CL2, and SmTAT ( Figure 5). On the contrary, the overexpression of SmMYB4 downregulated these genes except for Sm4CL2. The highest levels of SmPAL1, SmC4H, Sm4CL2, and SmTAT were 5.46, 2.67, 4.68, and 10.28 times their respective levels observed in the WT. The lowest levels of SmPAL1, SmC4H, and SmTAT were 41.45%, 9.89%, and 17.02% of their respective levels in the WT. In the tanshinone biosynthetic pathway, SmDXS1, SmDXS3, SmGGPPS1, and SmGGPPS3 were significantly upregulated in all interference lines, while only SmDXS1 was downregulated in all overexpressing lines. The transcription levels of other genes increased in the overexpressing lines but were not significantly different from their levels in the WT.  and SmHPPR in the salvianolic biosynthetic pathway, respectively. (F-L) represent the relative expression level ofSmGGPPS1, SmGGPPS3, SmDXS1, SmDXS3, SmHMGR1, SmHMGR2, and SmHMGR3 in the tanshinone biosynthetic pathway, respectively. (* p < 0.5, ** p < 0.01, *** p < 0.001).

SmMYB4 Possibly Binds to the Promoter of Key Enzyme Genes
The matrix {(G/A) (G/T) TG (G/T) T (G/A)} was selected to scan the promoters of SmPAL1, SmC4H, Sm4CL2, SmTAT, SmDXS1, SmDXS3, SmGGPPS1, and SmGGPPS3. Numerous possible SmMYB4 binding sites were found to exist within each promoter. (Supplementary Material Table S5). The Y1H assay demonstrated that SmMYB4 directly bind the promoters of SmTAT1, SmC4H, and SmPAL1 ( Figure 6). Therefore, we deduced that SmMYB4 might carry out its inhibitory function in the salvianolic acid and tanshinone biosynthetic pathways by binding to the promoter of these genes.

SmMYB4 Possibly Binds to the Promoter of Key Enzyme Genes
The matrix {(G/A) (G/T) TG (G/T) T (G/A)} was selected to scan the promoters of SmPAL1, SmC4H, Sm4CL2, SmTAT, SmDXS1, SmDXS3, SmGGPPS1, and SmGGPPS3. Numerous possible SmMYB4 binding sites were found to exist within each promoter. (Supplementary Material Table S5). The Y1H assay demonstrated that SmMYB4 directly bind the promoters of SmTAT1, SmC4H, and SmPAL1 ( Figure 6). Therefore, we deduced that SmMYB4 might carry out its inhibitory function in the salvianolic acid and tanshinone biosynthetic pathways by binding to the promoter of these genes. Figure 6. SmMYB4 binds the promoters of SmTAT1, SmC4H, and SmPAL1. The first and second rows represent the negative and positive controls, respectively. NEG control: pGADT7 + p53HIS2; POS control: pGADT7-p53 + p53HIS2. Figure 6. SmMYB4 binds the promoters of SmTAT1, SmC4H, and SmPAL1. The first and second rows represent the negative and positive controls, respectively. NEG control: pGADT7 + p53HIS2; POS control: pGADT7-p53 + p53HIS2.

Discussion
S. miltiorrhiza, an important medicinal material for treating cardiovascular and cerebrovascular diseases, has a growing market demand. The content of the bioactive components in the roots is one of the most crucial quality evaluation indexes. The exclusive information demonstrated that the secondary metabolites biosynthesis is regulated by diverse transcription factors and microRNAs [54,55]. The R2R3-MYB transcriptional factors are a large gene family in plants and have been reported to play crucial roles in regulating secondary metabolites. They can positively and negatively regulate the generation of secondary metabolites.
The R2R3 MYB proteins are more conserved in the N-terminal DNA binding domains, while the protein sequences at the C-terminus are divergent and believed to be responsible for several kinds of regulatory functions. Based on the specific conserved regions in the C-terminus, R2R3 (MYBs TFs) have been classified into many subgroups, in which subgroup 4 R2R3-MYB transcriptional factors function as negative regulators in the phenylpropanoid and flavonoid biosynthesis [56]. Many subgroup 4 R2R3-MYB transcription factors have been identified as transcriptional repressors. AtMYB4 suppressed the expression of AtC4H and resulted in the downregulated accumulation of the sinapic acid [57]. The ectopic overexpression of PvMYB4 from Panicum virgatum in the transgenic switch grass resulted in a reduced lignin content [58]. BrMYB4 negatively regulated the BrC4H gene involved in phenylpropanoid biosynthesis [59]. NtMYB2 and NtMYB3 repressed the biosynthesis of anthocyanin and flavonols in the Chinese narcissus, respectively [26,60]. SmMYB39 and SmMYB36 negatively regulated the phenolic acid and flavonoid biosynthesis in S. miltiorrhiza [45,46]. In the current study, the overexpression of SmMYB4 significantly decreased the salvianolic acid and tanshinone content, while its silencing enhanced their accumulation, suggesting that SmMYB4 negatively regulates their biosynthesis. Moreover, a C2 motif that is believed to be the key characteristic responsible for the transcriptional inhibition was harbored in the C terminus of SmMYB4.
The variations in the salvianolic acid and tanshinone content are closely related to the expression of genes encoding key enzymes involved in their biosynthesis. The expression levels of SmPAL1, SmC4H, Sm4CL2, and SmTAT were significantly downregulated in the SmMYB4-overexpressing transgenic lines and markedly upregulated in the RNAi transgenic lines. Combined with the promoter analysis, it was concluded that SmMYB4 affects the accumulation of the salvianolic acids by binding to the promoters of the genes mentioned above. Nevertheless, SmPAP1 activated the promoters of two key genes, SmPAL1 and SmC4H [43]. SmMYB98 activated the promoter of SmPAL1 and SmRAS1 in the salvianolic acid biosynthesis [47]. The overexpression of SmMYB36 downregulated most genes in the phenylpropanoid pathway (SmPAL1, SmC4H1, Sm4CL2) and tyrosine pathway (SmTAT1, SmHPPR1). However, the expression level of RAS1 and CYP98A14 did not change significantly. A similar phenomenon occurred in the SmMYB36-overexpressing lines [46]. The expression of SmPAL1, SmPAL2, SmPAL3, SmC4H, Sm4CL2, Sm4CL3, SmRAS1, and CYP98A14 was significantly induced in the SmMYB111-overexpressing lines but markedly reduced in the corresponding silenced lines [44]. The overexpression of SmMYB9b in the hairy roots significantly upregulated the transcription of SmDXS2, SmDXR, SmGGPPS, and SmKSL1, suggesting it enhanced the tanshinone synthesis through the stimulation of the MEP (methylerythritol phosphate) pathway [48]. The expression of SmDXS1, SmDXS2, SmDXR, SmMCT, SmMDS, SmHDS, SmCMK, SmHDR1, SmHMGR2, SmGGPPS1, SmCPS1, SmKSL1, and SmCYP76AHI was enhanced when SmMYB36 was overexpressed [46]. SmMYB98 activated the expression of SmGGPPS1 in the tanshinone biosynthesis [47]. In our study, the expression level of SmDXS1, SmDXS3 SmGGPPS1, and SmGGPPS3 increased significantly in the RNAi transgenic lines but was not changed in the overexpressing lines. We conclude that SmMYB4 affects the accumulation of the tanshinones by negatively regulating the expression of these key enzyme genes. The inability to detect any significant changes in the expression of other genes may be due to the competitive binding of the activator-type MYBs and the repressor-type MYBs to those targets (Figure 7).  The secondary metabolism in plants involves a complex three-dimensional regulatory network. A pathway or a key enzyme gene is usually regulated by many different transcription factors, which can interact and coordinatively regulate the gene expression. For example, PavMYB10.1 interacts with PavbHLH and PavWD40, and binds to the promoter regions of the anthocyanin biosynthesis genes PavANS and PavUFGT in sweet The secondary metabolism in plants involves a complex three-dimensional regulatory network. A pathway or a key enzyme gene is usually regulated by many different transcription factors, which can interact and coordinatively regulate the gene expression. For example, PavMYB10.1 interacts with PavbHLH and PavWD40, and binds to the promoter regions of the anthocyanin biosynthesis genes PavANS and PavUFGT in sweet cherry [61]. The anthocyanin synthesis in the lotus is regulated by a NnMYB5-NnbHLH1-NnTTG1 complex [23]. The MYB activator WHITE PETAL1 (WP1) physically interacts with MtTT8 and MtWD40-1 proteins to form the conserved MYB-basic-helix-loop-helix-WD40 regula-tory module, and regulates the floral carotenoid pigmentation through the transcriptional activation of the carotenoid biosynthetic genes in Medicago truncatula [29]. SmPAP1 regulates the phenolic acid biosynthetic pathway in S. miltiorrhiza by interacting with SmMYC2 and activating SmPAL1 and SmC4H [43]. SmMYB111 positively regulates the phenolic acid biosynthesis in S. miltiorrhiza by interacting with SmTTG1 and SmbHLH51 [44]. However, the specific mechanism of SmMYB4 needs further study. Our work will entail not only the elucidation of the secondary metabolic networks, but also enhance the understanding of the gene regulatory networks integrated with the metabolic networks.

Conclusions
In summary, a subgroup 4 R2R3 MYB transcription factor gene named SmMYB4 was cloned and functionally identified in the present study. Six overexpressions and eight RNA interference transgenic lines were successfully obtained via the Agrobacteriummediated leaf-disk transformed method. Based on the analysis of the content of the phenolic acids and tanshinones in transgenic plants and the expression of genes in the related pathway, we speculated that SmMYB4 negatively regulates the biosynthesis of the phenolic acids by repressing the expression of SmPAL1, SmC4H, Sm4CL2, and SmTAT, and simultaneously negatively regulates the biosynthesis of tanshinones by inhibiting the expression of SmDXS1, SmDXS3 SmGGPPS1, and SmGGPPS3. SmMYB4 may directly bind to the promoter of these genes. Our results will be helpful in better understanding the regulation mechanism of the phenolic acids and tanshinones production in S. miltiorrhiza and provide a train of thought to improve the contents of bioactive compounds in this traditional herbal.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/metabo12100968/s1, Table S1: Primers used in the SmMYB4-OE and SmMYB4-RNAi vector construction; Table S2: Primers for the molecular detection of transgenic plants; Table S3: Gradient conditions of the mobile phase of HPLC; Table S4: Primers for the key enzyme genes in the salvianolic acid and tanshinone biosynthetic pathways; Table S5: Information of the SmMYB4 binding sites in the promoters of the key enzyme genes; Figure S1: Subcellular localization of SmMYB4 in onion epidermal cells; Figure S2: Restriction endonuclease digestion of the overexpression and RNAi expression vectors. A represents the Bgl II/BstP I digestion of the overexpression vectors, in which P represents plasmids, M refers to Marker DL2000, and 1-4 represents the digest results of Bgl II/BstP I. B represents the digestion of the RNAi expression vectors, in which P represents plasmids, M refers to Marker DL2000, 1 refers to the digestion of Xho I, and 2 refers to the digestion of Not I. The target fragment is marked with a red line or arrows; Figure S3: Detection of the transgenic lines at the DNA level Detection. Results of the overexpression and interference transgenic lines were individually shown in Figure S3A