Expression Patterns and Functional Analysis of Three SmTAT Genes Encoding Tyrosine Aminotransferases in Salvia miltiorrhiza

Tyrosine aminotransferase (TAT, E.C. 2.6.1.5) is a pyridoxal phosphate-dependent aminotransferase that is widely found in living organisms. It catalyzes the transfer of the amino group on tyrosine to α-ketoglutarate to produce 4-hydroxyphenylpyruvic acid (4-HPP) and is the first enzyme for tyrosine degradation. Three SmTATs have been identified in the genome of Salvia miltiorrhiza (a model medicinal plant), but their information is very limited. Here, the expression profiles of the three SmTAT genes (SmTAT1, SmTAT2, and SmTAT3) were studied. All three genes expressed in different tissues and responded to methyl jasmonate stimuli. SmTAT proteins are localized in the cytoplasm. The recombinant SmTATs were subjected to in vitro biochemical properties. All three recombinant enzymes had TAT activities and SmTAT1 had the highest catalytic activity for tyrosine, followed by SmTAT3. Also, SmTAT1 preferred the direction of tyrosine deamination to 4-HPP, while SmTAT2 preferred transamination of 4-HPP to tyrosine. In parallel, transient overexpression of SmTATs in tobacco leaves revealed that all three SmTAT proteins catalyzed tyrosine to 4-HPP in vivo, with SmTAT1 exhibiting the highest enzymatic activity. Overall, our results lay a foundation for the production of tyrosine-derived secondary metabolites via metabolic engineering or synthetic biology in the future.


Introduction
Salvia miltiorrhiza Bunge is a well-known medicinal plants in the family Lamiaceae, and its dried root has used for the treatment of chest pains and cardiovascular diseases [1].Its main active ingredients include fat-soluble tanshinones and water-soluble phenolic acids [2].Phenolic acids, especially rosmarinic acid (RA) and salvianolic acid B (SalB), have strong scavenging ability for free radicals and hypolipidemic, anti-inflammatory, and analgesic effects [3].Identifying genes involved in the biosynthesis and regulation of phenolic acids in S. miltiorrhiza has attracted increasing attention recently [4].As far as we know, 29 genes involved in the biosynthesis of SalB have been identified in the S. miltiorrhiza genome, including three SmTATs [5].
In most microbes, 4-HPP is the intermediate of the tyrosine biosynthetic pathway, and TATs are usually responsible for the final step of tyrosine biosynthesis from 4-HPP [11].However, TAT activity has been proposed to be involved in tyrosine catabolism by deamination rather than tyrosine synthesis in most plants [12].The broad substrate specificity of TAT allows it to use glutamate and phenylalanine as amino group donors and p-hydroxyphenylpyruvate, phenylpyruvate, and alpha-ketocaproic acid as amino group acceptors [13].The catalytic condition of the enzyme is mild and the catalytic efficiency is high without regeneration of the coenzyme cycle [14].
The biochemical and structural properties of TATs from microbes and mammals have been understood in considerable detail [13].However, the characterization and function of TAT in plants has not been well studied.Genes encoding TATs belong to a small multigene family.There are eight putative TAT genes in the Arabidopsis genome [12], four TAT genes in Malus domestica [11], and three TATs in S. miltiorrhiza [5].In Arabidopsis, AtTAT1 had a higher activity towards tyrosine compared with AtTAT2.AtTAT1 and AtTAT2 have both distinct and shared functions in tyrosine metabolism and AtTAT1 plays a major role in tocopherol biosynthesis compared to AtTAT2 [15].
TAT plays an important role in plant growth and development.TAT gene is responds to multiple abiotic stresses and is induced by methyl jasmonate (MeJA) [16].Drought and low nitrogen stress increase TAT activity in poplar roots, thereby affecting RA content [17].Most recent studies show that AtTAT1 is involved in tyrosine catabolism and contributes to plant survival in the dark [18].MdTAT2 overexpression increases resistance to drought and osmotic stress in Malus domestica [11].
In this study, three SmTAT genes in S. miltiorrhiza were cloned and their expression patterns were analyzed.The recombinant proteins of the three SmTATs were subjected to in vitro catalytic reactions and their substrate specificity was compared.Furthermore, their catalytic activities toward tyrosine were confirmed in vivo by transient over-expression system in the leaves of Nicotiana benthamiana.To our knowledge, this is the first study about biochemical characterization of TAT in S. miltiorrhiza.The results provide valuable information for the production of tyrosine-derived secondary metabolites, such as SalB and RA, via metabolic engineering in the future.

Bioinformatics Analysis of SmTATs in S. miltiorrhiza
Gene structure analysis of the three SmTATs was shown in Figure 1A.There are 5 introns in the SmTAT1, and 6 in SmTAT2 and SmTAT3.The ORF of SmTAT1, SmTAT2, and SmTAT3 was 1236 bp, 1308 bp, and 1260 bp, respectively.Amino acid sequences of the three SmTATs and TATs from Arabidopsis, N. benthamiana, Oryza sativa, and Malus domestica were used to create a phylogenetic tree.As shown in Figure 1B, all the TATs were clustered into three groups.SmTAT1, SmTAT3 and AtTAT1 (At5g53970) were clustered in group II, while SmTAT2 and AtTAT2 (At5g36160) were clustered together.The promoter sequences of SmTATs were analyzed using the PlantCare database.Cis-elements such as ABRE (abscisic acid response element) and some light response elements were found in the promoter regions of all the three SmTATs (Supplementary Table S1), suggesting that SmTATs may respond to phytohormone and environmental stress during the growth and development of S. miltiorrhiza.

Expression Patterns of SmTATs
The qRT-PCR results showed that SmTAT1, SmTAT2, and SmTAT3 were expressed in different tissues, including roots, stems, leaves, and flowers.And all the three members had higher transcription levels in the leaves (Figure 2A).In addition, the expression levels of all the three SmTATs were significantly increased under the treatment of MeJA (Figure 2B).

Expression Patterns of SmTATs
The qRT-PCR results showed that SmTAT1, SmTAT2, and SmTAT3 were expressed in different tissues, including roots, stems, leaves, and flowers.And all the three members had higher transcription levels in the leaves (Figure 2A).In addition, the expression levels of all the three SmTATs were significantly increased under the treatment of MeJA (Figure 2B).
Furthermore, we obtained transgenic Arabidopsis expressing proSmTATs::GUS to investigate the spatiotemporal expression patterns of the three SmTATs.GUS staining indicated that the spatiotemporal expression patterns of the three SmTATs were basically consistent, and GUS signals could be observed at different developmental stages (Figures 3-5).Strong GUS signals were obvious in all the young seedlings of Arabidopsis expressing GUS driven by proSmTAT1, proSmTAT2, or proSmTAT3.At the flowering stage, the GUS expression profiles in the root was consistent for Arabidopsis expressing proSmTATs::GUS, with obvious signal in the lateral root, while no GUS signal was observed in the taproot; The GUS signal in the stem leaf was stronger than that in the rosette leaf for Arabidopsis expressing proSmTAT1::GUS; the GUS signal in the flower driven by proSmTAT2 was comparatively weaker, while the signal in the stem driven by proSmTAT3 was comparatively stronger.
The qRT-PCR results showed that SmTAT1, SmTAT2, and SmTAT3 were expressed in different tissues, including roots, stems, leaves, and flowers.And all the three members had higher transcription levels in the leaves (Figure 2A).In addition, the expression levels of all the three SmTATs were significantly increased under the treatment of MeJA (Figure 2B).Tukey's comparisons) tested for significant differences between means (indicated by different letters at p < 0.05).* and ** represent a significant difference at p < 0.05 and p < 0.01 levels, respectively, compared with the control.
Furthermore, we obtained transgenic Arabidopsis expressing proSmTATs::GUS to investigate the spatiotemporal expression patterns of the three SmTATs.GUS staining indicated that the spatiotemporal expression patterns of the three SmTATs were basically consistent, and GUS signals could be observed at different developmental stages (Figures 3-5).Strong GUS signals were obvious in all the young seedlings of Arabidopsis expressing GUS driven by proSmTAT1, proSmTAT2, or proSmTAT3.At the flowering stage, the GUS expression profiles in the root was consistent for Arabidopsis expressing proSmTATs::GUS, with obvious signal in the lateral root, while no GUS signal was observed in the taproot; The GUS signal in the stem leaf was stronger than that in the rosette leaf for Arabidopsis expressing proSmTAT1::GUS; the GUS signal in the flower driven by proSmTAT2 was comparatively weaker, while the signal in the stem driven by proSmTAT3 was comparatively stronger.

SmTATs Are Located in the Cytoplasm
To determine the subcellular localization of the three SmTAT proteins, we created the fusion expression vector 35Spro::SmTATs-eGFP.The fusion expression vector was transferred into Arabidopsis protoplasts, with 35Spro::eGFP serving as a positive control.As shown in Figure 6, the green fluorescence of GFP was uniformly distributed throughout the cells, whereas the green fluorescence of the SmTAT1/2/3-GFP fusion pro-

SmTATs Are Located in the Cytoplasm
To determine the subcellular localization of the three SmTAT proteins, we created the fusion expression vector 35Spro::SmTATs-eGFP.The fusion expression vector was transferred into Arabidopsis protoplasts, with 35Spro::eGFP serving as a positive control.As shown in Figure 6, the green fluorescence of GFP was uniformly distributed throughout the cells, whereas the green fluorescence of the SmTAT1/2/3-GFP fusion pro-

SmTATs Are Located in the Cytoplasm
To determine the subcellular localization of the three SmTAT proteins, we created the fusion expression vector 35Spro::SmTATs-eGFP.The fusion expression vector was transferred into Arabidopsis protoplasts, with 35Spro::eGFP serving as a positive control.As shown in Figure 6, the green fluorescence of GFP was uniformly distributed throughout the cells, whereas the green fluorescence of the SmTAT1/2/3-GFP fusion proteins was specifically localized in the cytoplasm.Our results indicated that SmTAT1/2/3 are cytoplasmic proteins.
teins was specifically localized in the cytoplasm.Our results indicated that SmTAT are cytoplasmic proteins.

Prokaryotic Expression and Enzyme Assays of SmTATs
To analyze the enzyme activities of SmTATs, the recombinant proteins we pressed in E. coli and purified by affinity chromatography column.The recombina sion proteins GST-SmTAT1, His-SmTAT2, and His-SmTAT3 were analyzed by SDS-(Figure 7).

Prokaryotic Expression and Enzyme Assays of SmTATs
To analyze the enzyme activities of SmTATs, the recombinant proteins were expressed in E. coli and purified by affinity chromatography column.The recombinant fusion proteins GST-SmTAT1, His-SmTAT2, and His-SmTAT3 were analyzed by SDS-PAGE (Figure 7).
teins was specifically localized in the cytoplasm.Our results indicated that SmTA are cytoplasmic proteins.

Prokaryotic Expression and Enzyme Assays of SmTATs
To analyze the enzyme activities of SmTATs, the recombinant proteins we pressed in E. coli and purified by affinity chromatography column.The recombina sion proteins GST-SmTAT1, His-SmTAT2, and His-SmTAT3 were analyzed by SDS-(Figure 7).When α-ketoglutarate was used as the amino acceptor, the affinity of SmTAT1 for tyrosine (K m of 0.31 mM) was higher than that of SmTAT2 and SmTAT3 (Table 1, Figure 8A).The catalytic efficiency of SmTAT1 for tyrosine was 1240 and 744 times higher than that of SmTAT2 and SmTAT3, respectively.When phenylalanine was used as the amino donor, the affinity of SmTAT1 for phenylalanine (K m of 3.08 mM) was higher than that of SmTAT2 and SmTAT3 (K m of 6.93 and 5.47 mM) (Table 1, Figure 8B), and the catalytic efficiency of SmTAT1 for phenylalanine were 6.88 and 39 times higher than that of SmTAT2 and SmTAT3, respectively.In the reverse reaction, when glutamate was used as the amino donor, the affinity of SmTAT2 for 4-HPP was highest (K m of 0.48 mM) (Table 1, Figure 9A), and the catalytic efficiency of SmTAT2 for 4-HPP was 0.92 and 3.33 times higher than that of SmTAT1 and SmTAT3, respectively.The affinity of SmTAT1 for phenylpyruvate was higher (K m of 0.71 mM) than that of SmTAT2 and SmTAT3 (Table 1, Figure 9B).The catalytic efficiency of SmTAT2 for phenylpyruvate was 326 and 203 times higher than that of SmTAT2 and SmTAT3, respectively.The catalytic activities of recombinant SmTATs to transfer the amino group on tyrosine to α-ketoglutarate in vitro were confirmed by HPLC.As shown in chromatogram (Figure 10), all the three SmTATs had tyrosine aminotransferase activity and SmTAT1 had the highest catalytic efficiency for tyrosine transamination to 4-HPP.UPLC-QTOF-MS was used to identify the reaction product 4-HPP.Overall, the catalytic efficiency of SmTAT1 was higher than SmTAT2 and SmTAT3.SmTAT1 showed the highest substrate affinity for tyrosine and preferred to transfer the amino group on tyrosine to α-ketoglutarate, while SmTAT2 preferred the reverse reaction to convert 4-HPP to tyrosine.
The catalytic activities of recombinant SmTATs to transfer the amino group on tyrosine to α-ketoglutarate in vitro were confirmed by HPLC.As shown in chromatogram (Figure 10), all the three SmTATs had tyrosine aminotransferase activity and SmTAT1 had the highest catalytic efficiency for tyrosine transamination to 4-HPP.UPLC-QTOF-MS was used to identify the reaction product 4-HPP.

Transient Overexpression of SmTATs Increased the Content of 4-HPP in the Leaves of N. benthamiana
To further verify the enzyme activity of SmTATs in vivo, we transiently overexpressed SmTATs in the leaves of N. benthamiana.Our qRT-PCR showed that a high transcription level of SmTAT1/2/3 was detected in the transgenic tobacco leaves (Figure 11A).Meanwhile, we determined the content of 4-HPP by UPLC-QTOF-MS.The results showed that the content of pHPL was significantly higher in SmTAT1/2/3-overexpressed leaves than the control (Figure 11B).The content of 4-HPP in the control group was 3.5 μg/g, whereas the content of 4-HPP in the leaves overexpressing SmTAT1, SmTAT2, and SmTAT3 were 18.24 μg/g, 10.11 μg/g, and 15.19 μg/g, respectively.The tyrosine aminotransferase activity of the three SmTATs was confirmed in N. benthamiana.

Transient Overexpression of SmTATs Increased the Content of 4-HPP in the Leaves of N. benthamiana
To further verify the enzyme activity of SmTATs in vivo, we transiently overexpressed SmTATs in the leaves of N. benthamiana.Our qRT-PCR showed that a high transcription level of SmTAT1/2/3 was detected in the transgenic tobacco leaves (Figure 11A).Meanwhile, we determined the content of 4-HPP by UPLC-QTOF-MS.The results showed that the content of pHPL was significantly higher in SmTAT1/2/3-overexpressed leaves than the control (Figure 11B).The content of 4-HPP in the control group was 3.5 µg/g, whereas the content of 4-HPP in the leaves overexpressing SmTAT1, SmTAT2, and SmTAT3 were 18.24 µg/g, 10.11 µg/g, and 15.19 µg/g, respectively.The tyrosine aminotransferase activity of the three SmTATs was confirmed in N. benthamiana.

Discussion
Since the publication of the genome sequence of S. miltiorrhiza, a well-known medicinal plant, there has been an increasing number of studies on the function of genes involved in the biosynthesis and regulation of active ingredients [19].Three genes encoding TAT have been identified in S. miltiorrhiza [5], but their information is very limited.In the present study, we systematically compared the expression patterns of the three genes and their biochemical characteristics were also investigated.The results lay a foundation for the production of tyrosine-derived nature compounds by metabolic engineering or synthetic biology.
Our qRT-PCR results showed that all the three SmTATs were highly expressed at different tissues of the two-year-old S. miltiorrhiza, with a higher level in the leaves.The TAT-mediated tyrosine metabolic pathway generates plastoquinone and ubiquinone, which function as electron and proton carriers in the electron transport chain of photosynthesis and respiration [20].We speculate that SmTAT1, SmTAT2, and SmTAT3 may play a role in photosynthesis and respiration.SmTAT1 showed the highest transcription level in the stems of two-month-old S. miltiorrhiza plantlet [21].We found that the expression level was lowest in the stem of two-year-old plant (Figure 2A).The promoter sequence determines the expression pattern of a gene [22].In this study, we obtained proSmTATs::GUS transgenic Arabidopsis and performed GUS staining.The GUS signals were basically consistent with the qRT-PCR results.
Jasmonates are signal molecules of defense/stress pathways and MeJA is an effective elicitor to enhance the production of phenolic acids in S. miltiorrhiza [23].Previous study indicated that the transcription levels of all the three SmTATs in the hairy roots of S. miltiorrhiza were significantly up-regulated under the treatment of MeJA [5], which are consistent with our results (Figure 2B).TAT functions in abiotic stress environment.For instance, three MdTATs in apple were induced by drought stress and confer tolerance to drought and osmotic stresses in plants [11].Whether SmTATs are responsive to drought stress in S. miltiorrhiza deserves further investigation.
Prokaryotic expression system has been widely used to obtain large amounts of recombinant proteins [24].However, the recombinant proteins expressed by prokaryotic

Discussion
Since the publication of the genome sequence of S. miltiorrhiza, a well-known medicinal plant, there has been an increasing number of studies on the function of genes involved in the biosynthesis and regulation of active ingredients [19].Three genes encoding TAT have been identified in S. miltiorrhiza [5], but their information is very limited.In the present study, we systematically compared the expression patterns of the three genes and their biochemical characteristics were also investigated.The results lay a foundation for the production of tyrosine-derived nature compounds by metabolic engineering or synthetic biology.
Our qRT-PCR results showed that all the three SmTATs were highly expressed at different tissues of the two-year-old S. miltiorrhiza, with a higher level in the leaves.The TAT-mediated tyrosine metabolic pathway generates plastoquinone and ubiquinone, which function as electron and proton carriers in the electron transport chain of photosynthesis and respiration [20].We speculate that SmTAT1, SmTAT2, and SmTAT3 may play a role in photosynthesis and respiration.SmTAT1 showed the highest transcription level in the stems of two-month-old S. miltiorrhiza plantlet [21].We found that the expression level was lowest in the stem of two-year-old plant (Figure 2A).The promoter sequence determines the expression pattern of a gene [22].In this study, we obtained proSmTATs::GUS transgenic Arabidopsis and performed GUS staining.The GUS signals were basically consistent with the qRT-PCR results.
Jasmonates are signal molecules of defense/stress pathways and MeJA is an effective elicitor to enhance the production of phenolic acids in S. miltiorrhiza [23].Previous study indicated that the transcription levels of all the three SmTATs in the hairy roots of S. miltiorrhiza were significantly up-regulated under the treatment of MeJA [5], which are consistent with our results (Figure 2B).TAT functions in abiotic stress environment.For instance, three MdTATs in apple were induced by drought stress and confer tolerance to drought and osmotic stresses in plants [11].Whether SmTATs are responsive to drought stress in S. miltiorrhiza deserves further investigation.
Prokaryotic expression system has been widely used to obtain large amounts of recombinant proteins [24].However, the recombinant proteins expressed by prokaryotic cells are expressed too fast, resulting in aggregation of the expressed protein molecules into inactive inclusion bodies [25].To obtain recombinant SmTATs, three E. coli expression systems (pGEX-4T-1, pET-28a, and pET-28a-MBP) were used in this study.Soluble SmTAT1 was successfully obtained by the prokaryotic expression vector pGEX-4T-1.The recombinant SmTAT2 and SmTAT3 tended to be expressed in the form of inclusion bodies.We finally obtain bioactive SmTAT2 and SmTAT3 from inclusion bodies expressed in pET-28a-MBP via direct dilution.
TAT catalyzes the reversible transamination between tyrosine and 4-HPP [7].TATs are responsible for the tyrosine biosynthesis from 4-HPP in most microbes [11], while TAT is involved in tyrosine catabolism by deamination rather than tyrosine synthesis in most plants [12].TAT can also use glutamate and phenylalanine as amino group donors and p-hydroxy-phenylpyruvate, phenylpyruvate, and alpha-ketocaproic acid as amino group acceptors [13].PvTAT from Prunella vulgaris had substrates preference for tyrosine and favored the conversion of tyrosine to 4-HPP [26].In Arabidopsis, AtTAT1 and AtTAT2 have a distinct substrate specificity and preferred reaction direction.AtTAT1 had a higher activity towards tyrosine than AtTAT2.Also, AtTAT1 favored the direction of tyrosine deamination to 4-HPP, whereas AtTAT2 preferred the reverse reaction [12].Here, we found that the catalytic efficiency of SmTAT1 is higher than SmTAT2 and SmTAT3.SmTAT1 showed the highest substrate affinity for tyrosine and preferred the transfer of amino group on tyrosine to α-ketoglutarate, while SmTAT2 preferred the reverse reaction to convert 4-HPP to tyrosine (Table 1 and Figure 8).Interesting, as shown in Figure 1, AtTAT1 and SmTAT1 (two members preferred the conversion of tyrosine to 4-HPP) are clustered in the same group, while AtTAT2 and SmTAT2 (preferred the reaction of 4-HPP to tyrosine) are clustered in another group (Figure 1).The biochemical characterization of the TATs support the phylogenetic evolutionary tree.It should be mentioned that the low concentrations of substrate are missing in kinetic analysis because of shortcomings of the experimental design.In future studies, we will overcome the shortcoming to obtain more accurate kinetic parameters.
TAT is a key enzyme in the biosynthesis of tyrosine-derived secondary metabolites in plants, such as phenolic acids, uronic acids, tocopherols, and alkaloids [27].The overexpression of PvTAT in Prunella vulgaris hairy roots increased the content of RA [26].Virus-induced silencing of TAT led to the reduction in isoquinoline alkaloid content in opium poppy [8].In the hairy roots of S. miltiorrhiza, overexpressing SmTAT1 alone did not significantly increase RA content, but the co-expression of SmTAT1 and SmHPPR1 resulted in a 4.3 times increase in RA content [28].For the two homologous TAT proteins in Arabidopsis, AtTAT1 plays a major role in tocopherol biosynthesis compared to AtTAT2 [15].In the present study, in vitro and in vivo experiments supported that SmTAT1 exhibited higher activity in converting tyrosine to 4-HPP.Many MeJA responsive transcription factors, such as SmMYC2 [29], SmMYB97 [30], and SmMYB111 [31], regulate the production of phenolic acids in S. miltiorrhiza by binding to the promoter of SmTAT1 and activate its expression.We speculated that SmTAT1 plays a dominant role in the biosynthesis of phenolic acids among the three homologous TAT proteins in S. miltiorrhiza.
In conclusion, three SmTAT genes in S. miltiorrhiza were cloned and their spatiotemporal expression profiles were comparatively analyzed.All the three genes significantly responded to MeJA stimuli and the SmTAT proteins are all located in the cytoplasm.In vitro catalytic reactions of the recombinant SmTATs indicated that all the three SmTATs had tyrosine aminotransferase activity and SmTAT1 had the highest catalytic efficiency for tyrosine transamination to 4-HPP.Transient overexpression of SmTATs in tobacco leaves indicated that all three SmTAT proteins have catalytic activities toward tyrosine.To further confirm their role to secondary metabolites, we will overexpress and knockout SmTAT genes in S. miltiorrhiza to analyze whether they affect the biosynthesis of rosmarinic acid in the future.

Plant Materials and Growth Conditions
Two-year-old Salvia miltiorrhiza in the experimental field of Shaanxi Normal University was collected at the flowering stage and the roots, stems, leaves, and flowers were separately frozen in liquid nitrogen for RNA extraction.Two-month-old S. miltiorrhiza plantlets were used for the treatment of methyl jasmonate (MeJA) as we described previously [32].
Arabidopsis thaliana ecotype Columbia-0 and Nicotiana benthamiana seedlings were cultured in a culture room at 23 • C and a humidity of 55% under 16 h of light and 8 h of darkness.While Arabidopsis used for subcellular localization experiment was cultured under short light condition (8 h light, 16 h dark).

Isolation of SmTAT1/2/3 and Sequence Analysis
Total RNA was extracted from two-month-old S. miltiorrhiza using the Plant Total RNA Extraction Kit (Vazyme, Nanjing, China) and then reverse transcribed into cDNA using Hiscript II reverse transcriptase (Vazyme, Nanjing, China).Genomic DNA was extracted from two-month-old S. miltiorrhiza using the DN15-Plant DNA Mini Kit (Aidlab, Beijing, China).The open reading frames (ORFs) of SmTAT1, SmTAT2, and SmTAT3 were amplified with SmTATs-F/R primers (Supplementary Table S2).The promoter regions of SmTAT1, SmTAT2, and SmTAT3, which is 1350 bp, 1125 bp, and 1495 bp, respectively, were amplified with ProSmTATs-F/R primers (Supplementary Table S2).The amplified fragments were inserted into the vector pMD19-T (TaKaRa, Dalian, China) and confirmed by sequencing.
A phylogenetic evolutionary tree of TATs from different species was created using MEGA X software (1000 bootstraps) by the neighbor-joining method, and then embellished by the online EvolView tool (http://www.evolgenius.info/evolview,accessed on 3 June 2023).The cis-elements in the promoter regions were predicted using the PlantCare online website (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/,accessed on 5 December 2022).
All primers used in this study are listed in Supplementary Table S2.

Expression Analysis via qRT-PCR
Total RNA was extracted from different tissues of S. miltiorrhiza and seedlings treated with MeJA and then reverse transcribed to cDNA.Real-time quantitative PCR (qRT-PCR) was performed on a Roche LightCycler 96 system (Roche Diagnostics GmbH, Basel, Switzerland) using SYBR Green qPCR mixture (Vazyme, Nanjing, China).SmUbiquitin was used as internal reference gene.Expression levels of SmTAT1/2/3 were analyzed via 2 −∆∆CT analysis [33].

Expression and Purification of SmTATs from Escherichia coli
The ORF of SmTAT1 was inserted into the EcoR I and Not I sites of the pGEX-4T-1 vector.The obtained vector pGEX-4T-1-SmTAT1 was transformed into Arctic Express (DE3) strain for protein expression.The transgenic Arctic Express strain was cultured, harvested, and lysed by ultrasound.The recombinant SmTAT1 was isolated and purified according to the previous method [34].
The ORFs of SmTAT2 and SmTAT3 were, respectively, inserted into the Sal I\Xho I sites and the BamH I\Sal I sites of the modified pET28a which carries a MBP tag (pET28a-MBP).The recombinant pET28a-MBP-SmTAT2 and pET28a-MBP-SmTAT3 were transformed into BL21 (DE3) strain for protein expression, respectively.The transgenic BL21 strain was cultured in 100 mL LB medium with 100 mg/L kanamycin at 37 • C until OD 600 = 0.3-0.4,and then isopropyl-β-d-thiogalactoside was added with the final concentration of 1mM to induce the expression of SmTAT2 and SmTAT3 at 16 • C until OD 600 = 0.8.The cell cultures were harvested by centrifugation and lysed with 30 mL buffer A (500 mM Tric-HCL, 0.15 mM NaCl, 5 mM EDTA, 1 mg/mL lysozyme, pH = 8).The lysed cells were subjected to centrifugation at 12,000× g for 15 min at 4 • C. The resulting pellets were washed with 20 mL buffer B (500 mM Tric-HCl, 0.15 mM NaCl, 2.5 M urea, 5 mM EDTA, pH = 8), and purified using the NI-NTA column (QualitYard, Beijing, China).The purified proteins were washed from the column with 2 mL 8 M urea, then denatured overnight at 4 • C with 40 mL buffer D (500 mM Tric-HCl, 0.15 mM NaCl, 5 mM EDTA, pH = 8).
The recombinant proteins were identified via SDS-PAGE.
Agroinfltration of N. benthamiana leaves with Rhizobium radiobacter strain EHA105, carrying either pEarleyGate202-SmTATs or pEarleyGate202 empty vector, was performed according to the previous protocols [36].The infiltrated leaves were cultured in the dark for 48 h and then collected for qRT-PCR to detect the expression levels of SmTATs.NtActin was used as the reference gene.
The conditions for MS analysis were as follows: The MSE scanning range was m/z 80−800, and the scanning time was 0.2 s.Ion source was set at the temperature of 120 • C, capillary voltage of 1.0 kV, cone voltage of 40 V, desolvent temperature of 400 • C, and collision voltage of 20 eV (4-HPP, t R = 0.562 min, m/z 181.14 [M-H] + ).

Heterologous Expression of proSmHPPRs::GUS and GUS Staining
The promoter regions of SmTAT1, SmTAT2, and SmTAT3 were, respectively, inserted into the pCAMBIA1391Z via double digestion to drive the expression of GUS.The Arabidopsis expressing ProSmHPPRs::GUS was acquired through the inflorescence dip method [38].The T 2 generation transgenic lines at different development stages were used for GUS staining as previously described [39].

Figure 2 .
Figure 2. Expression profiles of SmTATs in S. miltiorrhiza.(A) Transcription levels of SmTATs in different tissues.(B) Expression changes in SmTATs in response to MeJA treatments.All data are the means of three biological replicates, with error bars indicating SD.One-way ANOVA (followed by

Figure 2 .
Figure 2. Expression profiles of SmTATs in S. miltiorrhiza.(A) Transcription levels of SmTATs in different tissues.(B) Expression changes in SmTATs in response to MeJA treatments.All data are the means of three biological replicates, with error bars indicating SD.One-way ANOVA (followed by

Figure 2 .
Figure 2. Expression profiles of SmTATs in S. miltiorrhiza.(A) Transcription levels of SmTATs in different tissues.(B) Expression changes in SmTATs in response to MeJA treatments.All data are the means of three biological replicates, with error bars indicating SD.One-way ANOVA (followed by Tukey's comparisons) tested for significant differences between means (indicated by different letters at p < 0.05).* and ** represent a significant difference at p < 0.05 and p < 0.01 levels, respectively, compared with the control.

Figure 11 .
Figure 11.Transient over-expression of SmTAT1/2/3 in N. benthamiana leaves.(A) Transcription levels of SmTATs by qRT-PCR.(B) Determination of 4-HPP by UPLC-QTOF-MS.All data are the means of three biological replicates, with error bars indicating SD.One-way ANOVA (followed by Tukey's comparisons) tested for significant differences between means.Different letters indicate significant difference between the two samples at the level of p < 0.05.

Figure 11 .
Figure 11.Transient over-expression of SmTAT1/2/3 in N. benthamiana leaves.(A) Transcription levels of SmTATs by qRT-PCR.(B) Determination of 4-HPP by UPLC-QTOF-MS.All data are the means of three biological replicates, with error bars indicating SD.One-way ANOVA (followed by Tukey's comparisons) tested for significant differences between means.Different letters indicate significant difference between the two samples at the level of p < 0.05.