Overexpression of AtMYB2 Promotes Tolerance to Salt Stress and Accumulations of Tanshinones and Phenolic Acid in Salvia miltiorrhiza

Salvia miltiorrhiza is a prized traditional Chinese medicinal plant species. Its red storage roots are primarily used for the treatment of cardiovascular and cerebrovascular diseases. In this study, a transcription factor gene AtMYB2 was cloned and introduced into Salvia miltiorrhiza for ectopic expression. Overexpression of AtMYB2 enhanced salt stress resistance in S. miltiorrhiza, leading to a more resilient phenotype in transgenic plants exposed to high-salinity conditions. Physiological experiments have revealed that overexpression of AtMYB2 can decrease the accumulation of reactive oxygen species (ROS) during salt stress, boost the activity of antioxidant enzymes, and mitigate oxidative damage to cell membranes. In addition, overexpression of AtMYB2 promotes the synthesis of tanshinones and phenolic acids by upregulating the expression of biosynthetic pathway genes, resulting in increased levels of these secondary metabolites. In summary, our findings demonstrate that AtMYB2 not only enhances plant tolerance to salt stress, but also increases the accumulation of secondary metabolites in S. miltiorrhiza. Our study lays a solid foundation for uncovering the molecular mechanisms governed by AtMYB2 and holds significant implications for the molecular breeding of high-quality S. miltiorrhiza varieties.

DNA methylation regulates biosynthesis of tanshinones and phenolic acids during growth of Salvia miltiorrhiza

DNA methylation plays a crucial role in the regulation of plant growth and the biosynthesis of secondary metabolites. Danshen (Salvia miltiorrhiza) is a valuable Chinese herbal medicine commonly used to treat cardiovascular diseases; its active ingredients are tanshinones and phenolic acids, which primarily accumulate in roots. Here, we conducted a targeted metabolic analysis of S. miltiorrhiza roots at 3 distinct growth stages: 40 d old (r40), 60 d old (r60), and 90 d old (r90). The contents of tanshinones (cryptotanshinone, tanshinone I, tanshinone IIA, and rosmariquinone) and phenolic acids (rosmarinic acid and salvianolic acid B) gradually increased during plant development. Whole-genome bisulfite sequencing and transcriptome sequencing of roots at the 3 growth stages revealed an increased level of DNA methylation in the CHH context (H represents A, T, or C) context at r90 compared with r40 and r60. Increased DNA methylation levels were associated with elevated expression of various genes linked to epigenetic regulations, including CHROMOMETHYLASE2 (SmCMT2), Decrease in DNA Methylation 1 (SmDDM1), Argonaute 4 (SmAGO4), and DOMAINS REARRANGED METHYLTRANSFERASE 1 (SmDRM1). Moreover, expression levels of many genes involved in tanshinone and salvianolic acid biosynthesis, such as copalyldiphosphate synthase 5 (SmCPS5), cytochrome P450-related enzyme (SmCYP71D464), geranylgeranyl diphosphate synthase (SmGGPPS1), geranyl diphosphate synthase (SmGPPS), hydroxyphenylpyruvate reductase (SmHPPR), and hydroxyphenylpyruvate dioxygenase (SmHPPD), were altered owing to hyper-methylation, indicating that DNA methylation plays an important role in regulating tanshinone and phenolic acid accumulation. Our data shed light on the epigenetic regulation of root growth and the biosynthesis of active ingredients in S. miltiorrhiza, providing crucial clues for further improvement of active compound production via molecular breeding in S. miltiorrhiza.

Endophytic fungus Penicillium steckii DF33 promoted tanshinones biosynthesis in Salvia miltiorrhiza by regulating the expression of CYP450 genes

Salvia miltiorrhiza, a prominent traditional Chinese medicinal resource, has been extensively employed in the management of cardiovascular and cerebrovascular ailments. Ensuring the consistency of S. miltiorrhiza raw materials revolves around the imperative task of maintaining stable tanshinones content and composition. An effective approach in this regard involves the utilization of endophytic fungi as inducers. Within this context, our study spotlights an endophytic fungus, Penicillium steckii DF33, isolated from the roots of S. miltiorrhiza. Remarkably, this fungus has demonstrated a significant capacity to boost the biosynthesis and accumulation of tanshinones. The primary objective of this investigation is to elucidate the underlying regulatory mechanism by which DF33 enhances and regulates the biosynthesis and accumulation of tanshinones. This is achieved through its influence on the differential expression of crucial CYP450 genes within the S. miltiorrhiza hairy roots system. The results revealed that the DF33 elicitor not only promotes the growth of hairy roots but also enhances the accumulation of tanshinones. Notably, the content of cryptotanshinone was reached 1.6452 ± 0.0925 mg g-1, a fourfold increase compared to the control group. Our qRT-PCR results further demonstrate that the DF33 elicitor significantly up-regulates the expression of most key enzyme genes (GGPPS, CPS1, KSL1, CYP76AH1, CYP76AH3, CYP76AK1, CYP71D411) involved in the tanshinone biosynthesis pathway. This effect is particularly pronounced in certain critical CYP450 genes and Tanshinone ⅡA synthase (SmTⅡAS), with their expression levels peaking at 7 days or 14 days, respectively. In summary, endophytic P. steckii DF33 primarily enhances tanshinone biosynthesis by elevating the expression levels of pivotal enzyme genes associated with the modification and transformation stages within the tanshinone biosynthesis pathway. These findings underscore the potential of employing plant probiotics, specifically endophytic and root-associated microbes, to facilitate the biosynthesis and transformation of vital constituents in medicinal plants, and this approach holds promise for enhancing the quality of traditional Chinese medicinal materials.

Comparative transcriptomics of two Salvia subg. Perovskia species contribute towards molecular background of abietane-type diterpenoid biosynthesis

Tanshinones, are a group of diterpenoid red pigments present in Danshen - an important herbal drug of Traditional Chinese Medicine which is a dried root of Salvia miltiorrhiza Bunge. Some of the tanshinones are sought after as pharmacologically active natural products. To date, the biosynthetic pathway of tanshinones has been only partially elucidated. These compounds are also present in some of the other Salvia species, i.a. from subgenus Perovskia, such as S. abrotanoides (Kar.) Sytsma and S. yangii B.T. Drew. Despite of the close genetic relationship between these species, significant qualitative differences in their diterpenoid profile have been discovered. In this work, we have used the Liquid Chromatography-Mass Spectrometry analysis to follow the content of diterpenoids during the vegetation season, which confirmed our previous observations of a diverse diterpenoid profile. As metabolic differences are reflected in different transcript profile of a species or tissues, we used metabolomics-guided transcriptomic approach to select candidate genes, which expression possibly led to observed chemical differences. Using an RNA-sequencing technology we have sequenced and de novo assembled transcriptomes of leaves and roots of S. abrotanoides and S. yangii. As a result, 134,443 transcripts were annotated by UniProt and 56,693 of them were assigned as Viridiplantae. In order to seek for differences, the differential expression analysis was performed, which revealed that 463, 362, 922 and 835 genes indicated changes in expression in four comparisons. GO enrichment analysis and KEGG functional analysis of selected DEGs were performed. The homology and expression of two gene families, associated with downstream steps of tanshinone and carnosic acid biosynthesis were studied, namely: cytochromes P-450 and 2-oxoglutarate-dependend dioxygenases. Additionally, BLAST analysis revealed existence of 39 different transcripts related to abietane diterpenoid biosynthesis in transcriptomes of S. abrotanoides and S. yangii. We have used quantitative real-time RT-PCR analysis of selected candidate genes, to follow their expression levels over the vegetative season. A hypothesis of an existence of a multifunctional CYP76AH89 in transcriptomes of S. abrotanoides and S. yangii is discussed and potential roles of other CYP450 homologs are speculated. By using the comparative transcriptomic approach, we have generated a dataset of candidate genes which provides a valuable resource for further elucidation of tanshinone biosynthesis. In a long run, our investigation may lead to optimization of diterpenoid profile in S. abrotanoides and S. yangii, which may become an alternative source of tanshinones for further research on their bioactivity and pharmacological therapy.

A high-efficient protoplast transient system for screening gene editing elements in Salvia miltiorrhiza

KEY MESSAGE

A high-efficiency protoplast transient system was devised to screen genome editing elements in Salvia miltiorrhiza. Medicinal plants with high-value pharmaceutical ingredients have attracted research attention due to their beneficial effects on human health. Cell wall-free protoplasts of plants can be used to evaluate the efficiency of genome editing mutagenesis. The capabilities of gene editing in medicinal plants remain to be fully explored owing to their complex genetic background and shortfall of suitable transformation. Here, we took the Salvia miltiorrhiza as a representative example for developing a method to screen favorable gene editing elements with high editing efficiency in medical plants by a PEG-mediated protoplast transformation. Results indicated that using the endogenous SmU6.1 of S. miltiorrhiza to drive sgRNA and the plant codon-optimized Cas9 driven by the promoter SlEF1α can enhance the efficiency of editing. In summary, we uncover an efficacious transient method for screening editing elements and shed new light on increasing gene editing efficiency in medicinal plants.

Identification and characterization of a novel tyrosine aminotransferase gene (SmTAT3-2) promotes the biosynthesis of phenolic acids in Salvia miltiorrhiza Bunge

Rosmarinic acid (RA) and salvianolic acid B (SAB) are main phenolic acids in Salvia miltiorrhiza Bunge have been widely used in the treatment of cardiovascular and cerebrovascular diseases due to their excellent pharmacological activity. RA is a precursor of SAB, and tyrosine transaminase (TAT, EC 2.6.1.5) is a crucial rate-limiting enzyme in their metabolism pathway. This study identified a novel TAT gene, SmTAT3-2, and found that it is a new transcript derived from unconventional splicing of SmTAT3. We used different substrates for enzymatic reaction with SmTAT1, SmTAT3 and SmTAT3-2. Subcellular localization of SmTAT1 and SmTAT3-2 was completed based on submicroscopic techniques. In addition, they were overexpressed and CRISPR/Cas9 gene edited in hairy roots of S. miltiorrhiza. Revealed SmTAT3-2 and SmTAT1 showed a stronger affinity for L-tyrosine than SmTAT3, localized in the cytoplasm, and promoted the synthesis of phenolic acid. In overexpressed SmTAT3-2 hairy roots, the content of RA and SAB was significantly increased by 2.53 and 3.38 fold, respectively, which was significantly higher than that of overexpressed SmTAT1 strain compared with EV strain. These findings provide a valuable key enzyme gene for the phenolic acids metabolism pathway and offer a theoretical basis for the clinical application.

Characteristics of CXE family of Salvia miltiorrhiza and identification of interactions between SmGID1s and SmDELLAs

Carboxylesterase (CXE) is a class of hydrolases that contain an α/β folding domain, which plays critical roles in plant growth, development, and stress responses. Based on the genomic and transcriptomic data of Salvia miltiorrhiza, the SmCXE family was systematically analyzed using bioinformatics. The results revealed 34 SmCXE family members in S. miltiorrhiza, and the SmCXE family could be divided into five groups (Group I, Group II, Group III, Group IV, and Group V). Cis-regulatory elements indicated that the SmCXE promoter region contained tissue-specific and development-related, hormone-related, stress-related, and photoresponsive elements. Transcriptome analysis revealed that the expression levels of SmCXE2 were highest in roots and flowers (SmCXE8 was highest in stems and SmCXE19 was highest in leaves). Further, two GA receptors SmCXE1 (SmGID1A) and SmCXE2 (SmGID1B) were isolated from the SmCXE family, which are homologous to other plants. SmGID1A and SmGID1B have conserved HGGSF motifs and active amino acid sites (Ser-Asp-Val/IIe), which are required to maintain their GA-binding activities. SmGID1A and SmGID1B were significantly responsive to gibberellic acid (GA3) and methyl jasmonate (MeJA) treatment. A subcellular assay revealed that SmCXE1 and SmCXE2 resided within the nucleus. SmGID1B can interact with SmDELLAs regardless of whether GA3 exists, whereas SmGID1A can only interact with SmDELLAs in the presence of GA3. A Further assay showed that the GRAS domain mediated the interactions between SmGID1s and SmDELLAs. This study lays a foundation for further elucidating the role of SmCXE in the growth and development of S. miltiorrhiza.

Effect of interfering TOR signaling pathway on the biosynthesis of terpenoids in Salvia miltiorrhiza Bge

The TOR (Target of Rapamycin) signaling pathway, which takes TOR kinase as the core, regulates the absorption, distribution, and recycling of nutrients by integrating metabolic network and other signaling pathways, thus participating in the plant growth-defense trade-off. While terpenoids play an important role in plant growth, development, stress response, and signal transduction. The effect of the TOR signaling pathway on terpenoid biosynthesis in plants has yet to be studied in detail. In this study, the tissue culture seedlings of Salvia miltiorrhiza were treated with the TOR inhibitor AZD8055. The results show that the roots of the control group had begun to grow on the 8th day, while the seedlings treated with AZD8055 had no rooting signs. Combined with the expression changes of genes related to the TOR signaling pathway in the first 8 days, samples on the 3rd, 6th, and 8th days were selected for RNA-Seq analysis. Through RNA-Seq analysis, a total of 50,689 unigenes were obtained from the samples of these three periods, of which 4088 unigenes showed differential expression. The function enrichment and time-series analysis of differentially expressed genes (DEGs) showed that the main influence of the TOR signal pathway on plant growth-related processes was gradually transmitted with treatment time after TOR was inhibited. Pathway enrichment analysis of DEGs showed that the genes in the biosynthesis of terpenoids, such as diterpenoid and carotenoid biosynthetic pathways, could be regulated. Compared with other stages, DEGs related to terpenoid biosynthesis were mainly regulated in the S2 stage. In addition, the genes involved in terpenoid skeleton biosynthesis was also considerably enriched in the S2 stage, according to the results of gene set enrichment analysis (GSEA) of unigenes. Inhibition of the TOR signaling pathway may affect the biosynthesis of terpenoid signaling molecules, inhibit gibberellin's biosynthesis, and promote abscisic acid's biosynthesis. This study has discussed the effect of interfering with the TOR pathway on terpenoid biosynthesis in S. miltiorrhiza from the perspective of omics and provides new insight into the interaction between the terpenoid biosynthesis pathway and the growth-defense trade-off of medicinal plants.

Drug monomers from Salvia miltiorrhiza Bge. promoting tight junction protein expression for therapeutic effects on lung cancer

Salvia miltiorrhiza Bge. is a traditional Chinese medicine (TCM) that has been used for treatment of various diseases, including cancer by activating blood circulation and removing blood stasis. Tanshinone (TanIIA) and cryptotanshinone (CPT) are major lipophilic compounds extracted from the root of Salvia miltiorrhiza Bge., which are considered to be the effective compounds affecting the efficacy of the anti-tumor therapy of Salvia miltiorrhiza Bge. We have explored the mechanism of CPT and TanIIA exerting inhibition in non-small cell lung cancer (NSCLC) to provide experimental data support for guiding the translational development and clinical application of anti-tumor components of TCM. The subcutaneous tumor model and in vitro culture model of A549 cells was constructed to evaluate CPT and TanIIA's tumour-inhibitory effect respectively. RNA sequencing (RNA-seq) and bioinformatics analysis were conducted to identify differentially expressed genes (DEGs) and signalling pathways related to CPT and TanIIA treatment. qRT-PCR and Western blot were used to explore the mechanism of CPT and TanIIA intervention on NSCLC. Both CPT and TanIIA significantly inhibited the proliferation of A549 tumor cells and tumor growth in animal models. After intervention, the migration ability decreased and the level of apoptosis increased. RNA-seq results showed that both CPT and TanIIA could cause gene differential expression, miR-21-5p as one of the most significant gene expression differences between the two groups, and could act on cell connectivity. CPT and TanIIA play a regulatory role in regulating tight junction proteins (Occludin and ZO1), and Occludin mRNA and protein levels were reduced in an in vitro miR-21-5p overexpression A549 cell model. The mechanisms may be related to the reduction of miR-21-5p expression to increase the level of promoted tight junction protein expression for the purpose of inhibiting proliferation and invasion of NSCLC.

Genotype-Controlled Vertical Transmission Exerts Selective Pressure on Community Assembly of Salvia miltiorrhiza

The plant's endophytic fungi play an important role in promoting host development and metabolism. Studies have shown that the factors affecting the assembly of the endophyte community mainly include host genotype, vertical transmission, and soil origin. However, we do not know the role of vertically transmitted endohytic fungi influences on the host-plant's endophytic community assembly. Salvia miltiorrhiza from three production areas were used as research objects; we constructed three production area genotypes of S. miltiorrhiza regenerated seedlings simultaneously. Based on high-throughput sequencing, we analyzed the effects of genotype, soil origin, and vertical transmission on endophytic fungal communities. The results show that the community of soil origins significantly affected the endophytic fungal community in the regenerated seedlings of S. miltiorrhiza. The influence of genotype on community composition occurs through a specific mechanism. Genotype may selectively screen certain communities into the seed, thereby exerting selection pressure on the community composition process of offspring. As the number of offspring increases gradually, the microbiota, controlled by genotype and transmitted vertically, stabilizes, ultimately resulting in a significant effect of genotype on community composition.Furthermore, we observed that the taxa influencing the active ingredients are also selected as the vertically transmitted community. Moreover, the absence of an initial vertically transmitted community in S. miltiorrhiza makes it more vulnerable to infection by pathogenic fungi. Therefore, it is crucial to investigate and comprehend the selection model of the vertically transmitted community under varying genotypes and soil conditions. This research holds significant implications for enhancing the quality and yield of medicinal plants and economic crops.

Different techniques reveal the difference of community structure and function of fungi from root and rhizosphere of Salvia miltiorrhiza Bunge

Fungi have essential functions in plant health and performance. However, the plant-associated functions of many cultured fungi have not been established in detail. Here, the fungal species diversity in Salvia miltiorrhiza roots and rhizosphere was assessed for the first time using culturomics and high-throughput sequencing. We present a comprehensive functional metagenomic analysis of these fungi and verified activity of cellulase and chitinase predicted in the metagenomic analysis. We first collected and cultured fungi from the root and rhizosphere of S. miltiorrhiza. We found 92 species across 37 families and five phyla, with Ascomycota being dominant. Many rDNA internal transcribed spacer sequences could not be assigned to lower taxonomic levels. There were 19 genera of endophytic fungi and 37 genera of rhizosphere fungi. The culturomics approach had lower taxonomic diversity than high-throughput sequencing, but some fungi were only found in cultures. Structural analyses indicated that the dominant species differed in cultured and non-cultured samples at other levels, apart from the phylum level. Functional analysis mapped 223 carbohydrate enzyme families and 393 pathways in the CAZy and KEGG databases, respectively. The most abundant families were glycoside hydrolases and those involved in carbohydrate metabolism. As predicted by metagenomics, we experimentally verified cellulase and chitinase activity for 29 and 74 fungi, respectively. We provide the first evidence of biomass recycling by fungi that are associated with plants. Culturing is essential to reveal the hidden microbial community and critical functions in plant-microbe interactions.

Salicylic acid regulates phenolic acid biosynthesis via SmNPR1-SmTGA2/SmNPR4 modules in Salvia miltiorrhiza

Phenolic acids are the main active ingredients in Salvia miltiorrhiza, which can be used for the treatment of many diseases, particularly cardiovascular diseases. It is known that salicylic acid (SA) can enhance phenolic acid content, but the molecular mechanism of its regulation is still unclear. Nonexpresser of PR genes 1 (NPR1) plays a positive role in the SA signaling pathway. In this study, we identified a SmNPR1 gene that responds to SA induction and systematically investigated its function. We found that SmNPR1 positively affected phenolic acid biosynthesis. Then, we identified a novel TGA transcription factor, SmTGA2, which interacts with SmNPR1. SmTGA2 positively regulates phenolic acid biosynthesis by directly up-regulating SmCYP98A14 expression. After double-gene transgenic analysis and other biochemical assays, it was found that SmNPR1 and SmTGA2 work synergistically to regulate phenolic acid biosynthesis. In addition, SmNPR4 forms a heterodimer with SmNPR1 to inhibit the function of SmNPR1, and SA can alleviate this effect. Collectively, these findings elucidate the molecular mechanism underlying the regulation of phenolic acid biosynthesis by SmNPR1-SmTGA2/SmNPR4 modules and provide novel insights into the SA signaling pathway regulating plant secondary metabolism.

Comparative transcriptome analysis reveals the regulatory effects of exogenous auxin on lateral root development and tanshinone accumulation in Salvia miltiorrhiza

MAIN CONCLUSION

The physiological and transcriptome analysis revealed that auxin was a positive regulator of lateral root development and tanshinone accumulation in Salvia miltiorrhiza. Roots of S. miltiorrhiza are widely used as medicinal materials in China, and the root morphology and content of bioactive compounds [such as phenolic acids and diterpenoid quinones (tanshinones)] are the main factors to determine the quality of this herb. Auxin regulates root development and secondary metabolism in many plant species, but little is known about its function in S. miltiorrhiza. In this study, S. miltiorrhiza seedlings were treated (exogenous application) with the auxin indole-3-acetic acid (IAA) and the polar auxin transport inhibitor N-1-naphthylphthalamic acid (NPA) to investigate the regulatory roles of auxin in S. miltiorrhiza. The results indicated that exogenous IAA promoted both lateral root development and tanshinones biosynthesis in S. miltiorrhiza. The NPA application suppressed the lateral root development but showed no obvious effects on tanshinones accumulation. Based on the RNA-seq analysis, expressions of genes related to auxin biosynthesis and signaling transduction were altered in both treated groups. Coincidental with the enhanced content of tanshinones, transcripts of several key enzyme genes in the tanshinones biosynthetic pathway were stimulated after the exogenous IAA application. The expression profiles of seven common transcription factor domain-containing gene families were analyzed, and the results implied that some AP2/ERF genes were probably responsible for the auxin-induced lateral root development in S. miltiorrhiza. These findings shed new light on the regulatory roles of auxin on root development and bioactive compounds biosynthesis in S. miltiorrhiza, and lay the groundwork for future research into the detailed molecular mechanism underlying these biological functions.

SmJAZ4 interacts with SmMYB111 or SmMYC2 to inhibit the synthesis of phenolic acids in Salvia miltiorrhiza

Jasmonic acid (JA), as an important plant hormone, can induce the synthesis of phenolic acids in Salvia miltiorrhiza Bunge, a model medicinal plant, but the specific mechanism remains to be further elucidated. JA-responsive SmMYB111 positively regulates the biosynthesis of salvianolic acid B (SalB), but the molecular mechanism is unclear. Here, we found that SmMYB111 directly binds to the promoters of SmTAT1 and SmCYP98A14 and activates their transcription. Yeast two hybrid and bimolecular fluorescent complementation assay indicated that SmMYB111 interacts with SmJAZ4. Furthermore, we systematically characterized the function of SmJAZ4, which was highly expressed in flowers and roots and located in the nucleus and cell membrane. The contents of phenolic acids in the SmJAZ4-overexpressed transgenic plantlets and SmJAZ4-overexpressed transgenic hairy roots decreased significantly. SmJAZ4 interacts with SmMYC2 or SmMYB111 to repress their transcriptional activation activity on target enzyme genes of the biosynthesis pathway of phenolic acids. Overall, the molecular mechanism of SmJAZ4-SmMYC2/SmMYB111 module participating in JA signaling regulation of SalB biosynthesis was elucidated, which give a clue for the molecular regulation of phenolic acids biosynthesis in S. miltiorrhiza.

A tracking work on how Sm4CL2 re-directed the biosynthesis of salvianolic acids and tanshinones in Salvia miltiorrhiza hairy roots

Overexpression and antisense expression of Sm4CL2 re-directed the biosynthesis of salvianolic acids and tanshinones in Salvia miltiorrhiza hairy roots. Danshen (Salvia miltiorrhiza Bunge) is a widely used traditional Chinese medicine and its main active ingredients are water-soluble phenolic acids and lipophilic diterpenoids which are produced through the phenylpropanoid pathway and terpenoid pathway, respectively. 4-Coumaric acid: Coenzyme A ligase (4CL) is a key enzyme in the phenylpropanoid metabolism. We had obtained Sm4CL2-overexpressing (Sm4CL2-OE) and antisense Sm4CL2-expressing (anti-Sm4CL2) danshen hairy roots over ten years ago. In the follow-up study, we found that total salvianolic acids in Sm4CL2-OE-4 hairy roots increased to 1.35 times of the control-3, and that in anti-Sm4CL2-1 hairy roots decreased to 37.32% of the control-3, but tanshinones in anti-Sm4CL2-1 was accumulated to 1.77 ± 0.16 mg/g of dry weight, compared to undetectable in Sm4CL2-OE-4 and the control-3 hairy roots. Interestingly, Sm4CL2-OE-4 hairy roots contained more lignin, 1.36 times of the control-3, and enhanced cell wall and xylem lignification. Transcriptomic analysis revealed that overexpression of Sm4CL2 caused the upregulation of other phenylpropanoid pathway genes and antisense Sm4CL2 expression resulted in the downregulation of other phenylpropanoid pathway genes but activated the expression of terpenoid pathway genes like SmCYP76AK5, SmGPPS.SSUII.1 and SmDXS2. Protein-protein interaction analysis suggested that Sm4CL2 might interact with PAL, PAL4, CSE, CCoAOMT and SmCYP84A60, and appeared to play a key role in the interaction network. The tracking work in this study proved that Sm4CL2 could redirect both salvianolic acids and tanshinones biosynthesis possibly through synergistically regulating other pathway genes. It also indicated that genetic modification of plant secondary metabolism with biosynthetic gene might cause other responses through protein-protein interactions.

Genome-wide identification and co-expression network analysis of Aux/IAA gene family in Salvia miltiorrhiza

The auxin/indole-3-acetic acid (Aux/IAA) gene family serves as a principal group of genes responsible for modulating plant growth and development through the auxin signaling pathway. Despite the significance of this gene family, the identification and characterization of members within the well-known Chinese medicinal herb Salvia miltiorrhiza (S. miltiorrhiza) have not been thoroughly investigated. In this study, we employed bioinformatics methods to identify 23 Aux/IAA genes within the genome of S. miltiorrhiza. These genes were classified into typical IAA and atypical IAA based on their domain structure. Our analysis of the promoter regions revealed that the expression of these genes is regulated not only by auxins, but also by other hormones and environmental factors. Furthermore, we found that the expression patterns of these genes varied across various tissues of S. miltiorrhiza. While our initial hypothesis suggested that the primary function of these genes was the interaction between SmIAA and ARF, gene co-expression network analysis revealed that they are also influenced by various other transcription factors, such as WRKY and ERF. The findings establish a sturdy basis for future investigations into the function of the Aux/IAA gene family and exhibit promising prospects for enhancing the genetics of this medicinal flora and its associated species.

Overexpression of SmMYC2 enhances salt resistance in Arabidopsis thaliana and Salvia miltiorrhiza hairy roots

Soil salinity significantly affects both Salvia miltiorrhiza growth and development as well as seed germination throughout field cultivation and production. The basic helix-loop-helix (bHLH) transcription factor (TF) MYC2 contributes significantly to plant stress resistance as a key regulator of the jasmonic acid signaling pathway. In transgenic S. miltiorrhiza hairy roots, SmMYC2 has been shown to promote the accumulation of tanshinone and salvianolic acid, but its role in S. miltiorrhiza of resistance to abiotic stress is unclear. Herein, we found methyl jasmonate (MeJA), NaCl, and PEG treatment all significantly increased SmMYC2 expression. In response to salt stress, SmMYC2 overexpression in yeast increased its rate of growth. Additionally, overexpression of SmMYC2 transgenic Arabidopsis thaliana and S. miltiorrhiza hairy root showed that it might improve salt resistance in transgenic plant. In particular, compared to WT, overexpression of SmMYC2 transgenic Arabidopsis had higher levels of three antioxidant enzymes (superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT)), proline (Pro) content, and ABA-dependent and ABA-independent genes expression. They also had lower levels of malondialdehyde (MDA) and reactive oxygen species (ROS) accumulation. What's more, overexpression of SmMYC2 increases the expression of flavonoid synthesis genes and the accumulation of related components in Arabidopsis. These findings imply that SmMYC2 functions as a positive regulator that regulates plant tolerance to salt through ABA-dependent and independent signaling pathways.

Overexpression of SmLAC25 promotes lignin accumulation and decreases salvianolic acid content in Salvia miltiorrhiza

Laccase (LAC) is a blue multicopper oxidase that contains four copper ions, which is involved in lignin polymerization and flavonoid biosynthesis in plants. Although dozens of LAC genes have been identified in Salvia miltiorrhiza Bunge (a model medicinal plant), most have not been functionally characterized. Here, we explored the expression patterns and the functionality of SmLAC25 in S. miltiorrhiza. SmLAC25 has a higher expression level in roots and responds to methyl jasmonate, auxin, abscisic acid, and gibberellin stimuli. The SmLAC25 protein is localized in the cytoplasm and chloroplasts. Recombinant SmLAC25 protein could oxidize coniferyl alcohol and sinapyl alcohol, two monomers of G-lignin and S-lignin. To investigate its function, we generated SmLAC25-overexpressed S. miltiorrhiza plantlets and hairy roots. The lignin content increased significantly in all SmLAC25-overexpressed plantlets and hairy roots, compared with the controls. However, the concentrations of rosmarinic acid and salvianolic acid B decreased significantly in all the SmLAC25-overexpressed lines. Further studies revealed that the transcription levels of some key enzyme genes in the lignin synthesis pathway (e.g., SmCCR and SmCOMT) were significantly improved in the SmLAC25-overexpressed lines, while the expression levels of multiple enzyme genes in the salvianolic acid biosynthesis pathway were inhibited. We speculated that the overexpression of SmLAC25 promoted the metabolic flux of lignin synthesis, which resulted in a decreased metabolic flux to the salvianolic acid biosynthesis pathway.

Transcriptome analysis of the immunomodulatory effects of Salvia miltiorrhiza polysaccharide on hemocyte immune response in Procambarus clarkii

Natural plant polysaccharide as immune modulator is considered an effective strategy for healthy aquaculture to reduce medicine treatment. Salvia miltiorrhiza polysaccharides (SMP) had applications to regulate immune activity and enhance antioxidant in vertebrates, but the potential function has been rarely reported in crustaceans. In this study, the immunological effects of SMP on hemocytes of Procambarus clarkii were analyzed. Results showed that total superoxide dismutase (T-SOD), phenoloxidase (PO) activity and respiratory burst were up-regulated after SMP treatment. After high-throughput sequencing, 2170 differentially expressed genes (DEGs) including 1294 up-regulated and 876 down-regulated genes were identified. KEGG function enrichment analysis indicated that DEGs are involved in crustaceans cellular immune-related signaling pathways, including lysosome, phagosome and endocytosis. Transcriptome mining and qRT-PCR showed that SMP up-regulated humoral immunity factors gene expression. Diets supplemented with 0.8% SMP significantly up-regulated the total number of hemocytes (THC), T-SOD and PO activity, improved the survival of crayfish after Citrobacter freundii infection. This study suggested that SMP could improve the cellular and humoral immunity of P. clarkii. Furthermore, this finding supplied a molecular foundation for further comprehending the immunopotentiator effects of plant polysaccharides in crustaceans.

Spatiotemporal and Transcriptional Characterization on Tanshinone Initial Synthesis in Salvia miltiorrhiza Roots

Tanshinones are the bioactive constituents of Danshen (Salvia miltiorrhiza Bunge), which is used in Traditional Chinese Medicine to treat cardiovascular and other diseases, and they synthesize and accumulate in the root periderm of S. miltiorrhiza. However, there is no relevant report on the initial stage of tanshinone synthesis, as well as the root structure and gene expression characteristics. The present study aims to provide new insights into how these bioactive principles begin to synthesize by characterizing possible differences in their biosynthesis and accumulation during early root development from both spatial and temporal aspects. The morphological characteristics and the content of tanshinones in roots of S. miltiorrhiza were investigated in detail by monitoring the seedlings within 65 days after germination (DAGs). The ONT transcriptome sequencing was applied to investigate gene expression patterns. The periderm of the S. miltiorrhiza storage taproot initially synthesized tanshinone on about 30 DAGs. Three critical stages of tanshinone synthesis were preliminarily determined: preparation, the initial synthesis, and the continuous rapid synthesis. The difference of taproots in the first two stages was the smallest, and the differentially expressed genes (DEGs) were mainly enriched in terpene synthesis. Most genes involved in tanshinone synthesis were up regulated during the gradual formation of the red taproot. Plant hormone signal transduction and ABC transport pathways were widely involved in S. miltiorrhiza taproot development. Five candidate genes that may participate in or regulate tanshinone synthesis were screened according to the co-expression pattern. Moreover, photosynthetic ferredoxin (FD), cytochrome P450 reductase (CPR), and CCAAT binding transcription factor (CBF) were predicted to interact with the known downstream essential enzyme genes directly. The above results provide a necessary basis for analyzing the initial synthesis and regulation mechanism of Tanshinones.

Identification of Salvia miltiorrhiza Bunge with high and low cadmium accumulation and insight into the mechanisms of cadmium accumulation

Salvia miltiorrhiza Bunge is used as a Chinese herbal medicine (CHM) particularly its roots. No relevant reports at home and abroad have been made on the mechanism of cadmium (Cd) accumulation in S. miltiorrhiza. The Cd accumulation characteristics of S. miltiorrhiza from main cultivation areas in China were evaluated for the first time to obtain high and low Cd accumulation in S. miltiorrhiza roots. Results showed obvious differences in the Cd enrichment capacity of S. miltiorrhiza from different cultivation areas. We took the lead in identifying the germplasm resources of S. miltiorrhiza with high and low Cd accumulation, that is, S. miltiorrhiza roots from Pingyi Shangdong (SDPY) belongs to the resource with high Cd accumulation (SDPY-H) and that from Zhongjiang Sichuan (SCZJ) is the resources with low Cd accumulation (SCZJ-L) based on relevant physiological and biochemical indexes. Although the Cd content of SDPY-H was apparently higher than that from SCZJ-L, its translocation factor from root to aboveground part is significantly lower than that in SCZJ-L. Therefore, planting SCZJ-L is not only an economic and effective way to use Cd in slightly and moderately polluted soil, but also its aboveground part can be used for phytoremediation. Changes in Cd content before and after the use of transpiration inhibitor indicate that SDPY-H enriched Cd through the symplastic pathway, whereas SCZJ-L mainly enriched Cd through the apoplastic pathways. In addition, the role of the symplastic pathway in SCZJ-L is weaker than that in SDPY-H, which were preliminarily revealed by fluorescent quantitative polymerase chain reaction. The significant reduction of the SmNramps transcription expression amount is one of the reasons for the low Cd accumulation of SCZJ-L.

[Gene cloning, induction, and prokaryotic expression of a Sm14-3-3 protein from Salvia miltiorrhiza]

14-3-3 proteins are important proteins in plants, as they regulate plant growth and development and the response to biotic or abiotic stresses. In this study, a 14-3-3 gene(GenBank accession: OM683281) was screened from the cDNA library of the medicinal species Salvia miltiorrhiza by yeast two-hybrid and cloned. The open reading frame(ORF) was 780 bp, encoding 259 amino a cids. Bioinformatics analysis predicted that the protein was a non-transmembrane protein with the molecular formula of C_(1287)H_(2046)N_(346)O_(422)S_9, relative molecular weight of 29.4 kDa, and no signal peptide. Homologous sequence alignment and phylogenetic tree analysis proved that the protein belonged to 14-3-3 family and had close genetic relationship with the 14-3-3 proteins from Arabidopsis thaliana, Oryza sativa, and Nicotiana tabacum. The 14-3-3 gene was ligated to the prokaryotic expression vector pGEX-4 T-1 and then transformed into Escherichia coli BL21 for the expression of recombinant protein. Real-time fluorescent quantitative PCR showed that the expression of this gene was different among roots, stems, leaves, and flowers of S. miltiorrhiza. To be specific, the highest expression was found in leaves, followed by stems, and the lowest expression was detected in flowers. S. miltiorrhiza plants were treated with 15% PEG(simulation of drought), and hormones salicylic acid, methyl jasmonate, and ethephon, respectively, and the expression of 14-3-3 gene peaked at the early stage of induction. Therefore, the gene can quickly respond to abiotic stresses such as drought and plant hormone treatments such as salicylic acid, jasmonic acid, and ethylene. This study lays the foundation for revealing the molecular mechanism of 14-3-3 protein regulating tanshinone biosynthesis and responding to biotic and abiotic stresses.

Integrative transcriptomic and metabolomic analyses unveil tanshinone biosynthesis in Salvia miltiorrhiza root under N starvation stress

Salvia miltiorrhiza is a model plant for Chinese herbal medicine with significant pharmacologic effects due to its tanshinone components. Our previous study indicated that nitrogen starvation stress increased its tanshinone content. However, the molecular mechanism of this low nitrogen-induced tanshinone biosynthesis is still unclear. Thus, this study aimed to elucidate the molecular mechanism of tanshinone biosynthesis in S. miltiorrhiza under different N conditions [N-free (N0), low-N (Nl), and full-N (Nf, as control) conditions] by using transcriptome and metabolome analyses. Our results showed 3,437 and 2,274 differentially expressed unigenes between N0 and Nf as well as Nl and Nf root samples, respectively. N starvation (N0 and Nl) promoted the expression of the genes involved in the MVA and MEP pathway of tanshinone and terpenoid backbone biosynthesis. Gene ontology and KEGG analyses revealed that terpenoid backbone biosynthesis, hormone signal transduction, and phenylpropanoid biosynthesis were promoted under N starvation conditions, whereas starch and sucrose metabolisms, nitrogen and phosphorus metabolisms, as well as membrane development were inhibited. Furthermore, metabolome analysis showed that metabolite compounds and biosynthesis of secondary metabolites were upregulated. This study provided a novel insight into the molecular mechanisms of tanshinone production in S. miltiorrhiza in response to nitrogen stress.

Loss of anthocyanidin synthase gene is associated with white flowers of Salvia miltiorrhiza Bge. f. alba, a natural variant of S. miltiorrhiza

SmANS deletion leads to white flower mutation in Salvia miltiorrhiza. SmANS deletion leads to white flower mutation in Salvia miltiorrhiza. Abstract Salvia miltiorrhiza is an essential traditional Chinese medicine (TCM) with purple flowers, and S. miltiorrhiza Bge. f. alba is a unique intraspecific variation with white flowers. The molecular mechanism of flower color formation in S. miltiorrhiza will provide vital information for the variation and evolution. Here, we performed HPLC, transcriptomic, and re-sequencing analyses of purple-flowered S. miltiorrhiza line 'Zihua105' (ZH105) and white-flowered S. miltiorrhiza Bge. f. alba line 'Baihua18' (BH18). Delphinidin was the most anthocyanidin in ZH105, which become the main different between ZH105 vs. BH18 flowers. Transcriptome analysis revealed 299 differentially expressed genes (DEGs). SmANS, the anthocyanidin synthase gene in the down-stream anthocyanin biosynthesis pathway, was significantly expressed in ZH105 corollas, suggesting it might play a key role in white petal formation. Whole-genome re-sequencing revealed that a 6.75 kb segment located on chromosome 5, which contains the complete sequence of the SmANS genes, was lost in BH18 and another S. miltiorrhiza Bge. f. alba line. In contrast, the other five purple-flowered S. miltiorrhiza lines both possessed this segment. Further molecular marker identification also confirmed that wild S. miltiorrhiza Bge. f. alba lines lost regions that contained a complete or important part of SmANS sequences. Subsequently, the research showed that the deletion mutant of SmANS genes resulted in the natural white flower color variant of S. miltiorrhiza.

Genetic diversity of the cultivated Salvia miltiorrhiza populations revealed by four intergenic spacers

For better understanding the genetic diversity and phylogeny of the cultivated Salvia miltiorrhiza populations, four intergenic spacer sequences, ETS, psbA-trnH, trnL-trnF, and ycf1-rps15 of the 40 populations collected from China were Polymerase Chain Reaction (PCR) amplified, analyzed both individually and in combination. Haplotype diversity analysis showed that the cultivated S. miltiorrhiza populations had a very rich genetic diversity and an excellent capacity to resist environmental pressure. The best-fit nucleotide substitution models for ETS, psbA-trnH, trnL-trnF, ycf1-rps15, and their combined sequences were HKY+I, T92, T92, T92+G, and T92+G, respectively; the nucleotide conversion frequency in the combined sequences was lower than the transversion, and the relatively high nucleotide substitution frequencies suggests its high genetic variability. Neutral tests showed that the spacer sequences of the populations conform with the neutral evolution model, and there has been no current expansion events occurred. Phylogeny analyses based on both the individual and the combined sequences showed that the 40 populations were clustered in two clades with a very similar topological structure. The discrimination rate of the combined sequence marker is significantly increased to 52.5% (21 populations) over the highest 35% (13 populations) by the single marker of ETS, though still inadequate but a big step forward. Further exploration of more DNA markers is needed. This study for the first time revealed the rich genetic diversity and phylogeny of the currently cultivated S. miltiorrhiza populations in China and provides novel alternative molecular markers for the genetic identification and resources evaluation of the cultivated S. miltiorrhiza populations.

A 2-oxoglutarate-dependent dioxygenase converts dihydrofuran to furan in Salvia diterpenoids

Tanshinone ⅡA (TⅡA), a diterpene quinone with a furan ring, is a bioactive compound found in the medicinal herb redroot sage (Salvia miltiorrhiza Bunge), in which both furan and dihydrofuran analogs are present in abundance. Progress has been made recently in elucidating the tanshinone biosynthetic pathway, including heterocyclization of the dihydrofuran D-ring by cytochrome P450s; however, dehydrogenation of dihydrofuran to furan, a key step of furan ring formation, remains uncharacterized. Here, by differential transcriptome mining, we identified six 2-oxoglutarate-dependent dioxygenase (2-ODD) genes whose expressions corresponded to tanshinone biosynthesis. We showed that Sm2-ODD14 acts as a dehydrogenase catalyzing the furan ring aromatization. In vitro Sm2-ODD14 converted cryptotanshinone to TⅡA and thus was designated TⅡA synthase (SmTⅡAS). Furthermore, SmTⅡAS showed a strict substrate specificity, and repression of SmTⅡAS expression in hairy root by RNAi led to increased accumulation of total dihydrofuran-tanshinones and decreased production of furan-tanshinones. We conclude that SmTⅡAS controls the metabolite flux from dihydrofuran- to furan-tanshinones, which influences medicinal properties of S. miltiorrhiza.

De novo Biosynthesis of Salvianolic Acid B in Saccharomyces cerevisiae Engineered with the Rosmarinic Acid Biosynthetic Pathway

Salvianolic acid B (SAB), also named lithospermic acid B, belongs to a class of water-soluble phenolic acids, originating from plants such as Salvia miltiorrhiza. SAB exhibits a variety of biological activities and has been clinically used to treat cardio- and cerebrovascular diseases and also has great potential as a health care product and medicine for other disorders. However, its biosynthetic pathway has not been completely elucidated. Here, we report the de novo biosynthesis of SAB in Saccharomyces cerevisiae engineered with the heterologous rosmarinic acid (RA) biosynthetic pathway. The created pathway contains seven genes divided into three modules on separate plasmids, pRS424-FjTAL-Sm4CL2, pRS425-SmTAT-SmHPPR or pRS425-SmTAT-CbHPPR, and pRS426-SmRAS-CbCYP-CbCPR. These three modules were cotransformed into S. cerevisiae, resulting in the recombinant strains YW-44 and YW-45. Incubation of the recombinant strains in a basic medium without supplementing any substrates yielded 34 and 30 μg/L of SAB. The findings in this study indicate that the created heterologous RA pathway cooperates with the native metabolism of S. cerevisiae to enable the de novo biosynthesis of SAB. This provides a novel insight into a biosynthesis mechanism of SAB and also lays the foundation for the production of SAB using microbial cell factories.

Characterization of NAC family genes in Salvia miltiorrhiza and NAC2 potentially involved in the biosynthesis of tanshinones

The NAC (NAM, ATAF, and CUC) family members are specific transcription factors in plants. The large family is involved in many plant growth and developmental processes, as well as in abiotic/biotic stress responses. It has been well studied in the genomes of various plants, including Arabidopsis thaliana, tomato, and quinoa. However, identification and functional studies of NAC family members in medicinal Salvia miltiorrhiza are limited. Here, we systematically identified 84 NAC genes and named them according to their gene IDs in the recently sequenced genome. The phylogeny of NAC family protein sequences was analyzed using bioinformatics methods, which divided them into nine subfamilies. Then, their chromosomal locations, gene structures and conserved domains were analyzed comprehensively. To further investigate the regulatory functions of NACs in S. miltiorrhiza, we analyzed the response of 10 selected NAC genes to methyl jasmonate and used NAC2 for transgenic experiments. The overexpression of Sm-NAC2 decreased the tanshinone I and IIA contents by 56% and 62%, respectively. However, Sm-NAC2-RNAi promoted the accumulation of four tanshinones, tanshinone I, tanshinone IIA, cryptotanshinone, and dihydrotanshinone I, which increased 3.68-, 4.1-, 3.13- and 5.9- fold, respectively, compared with wild type. In the tanshinone biosynthetic pathways, the overexpression of Sm-NAC2 down-regulated CYP76AH1, and the silencing of Sm-NAC2 up-regulated the expression levels of HMGR1, DXS2, KSL2, and CYP76AH1. This study provides information on the evolution of Sm-NAC genes and their possible functions, and it lays a foundation for further research into the NAC family-associated regulation of tanshinone biosynthesis.

Transcription factor SmSPL7 promotes anthocyanin accumulation and negatively regulates phenolic acid biosynthesis in Salvia miltiorrhiza

Plant-specific SQUAMOSA promoter-binding protein-like (SPL) transcription factors play critical regulatory roles during plant growth and development. However, the functions of SPLs in Salvia miltiorrhiza (SmSPLs; a model medicinal plant) have not been reported. Here, the expression patterns and functions of SmSPL7 were characterized in S. miltiorrhiza. SmSPL7 was expressed in all parts of S. miltiorrhiza, with the highest expression level in the leaves, and could be inhibited by multiple hormones, including methyl jasmonate, auxin, abscisic acid, and gibberellin. SmSPL7 is localized within the nucleus and exhibits robust transcriptional activation activity. Transgenic lines overexpressing SmSPL7 demonstrated pronounced growth inhibition, accompanied by increased anthocyanin accumulation via the genetic activation of the anthocyanin biosynthesis pathway. However, SmSPL7 overexpression significantly decreased salvianolic acid B (SalB) production by inhibiting the transcripts of genes implicated in its biosynthesis pathway. Further analysis indicated that SmSPL7 directly binds to SmTAT1 and Sm4CL9 promoters and blocks their expression to inhibit the biosynthesis of SalB. Taken together, these results indicate that SmSPL7 is a negative regulator of SalB biosynthesis but positively regulates anthocyanin accumulation in S. miltiorrhiza. These findings provide new insights into the functionality of the SPL family while establishing an important foundation for further uncovering the crucial roles of SmSPL7 in the growth of S. miltiorrhiza.

The ERF-VII transcription factor SmERF73 coordinately regulates tanshinone biosynthesis in response to stress elicitors in Salvia miltiorrhiza

Here, we investigate the role of SmERF73, a group VII ETHYLENE RESPONSE FACTOR stress response transcription factor, in the regulation of post-modification of the skeleton precursors of diterpene tanshinones in Salvia miltiorrhiza. Most genes found to be involved in tanshinone biosynthesis are located on chromosome 6, and five of these genes comprise a large gene cluster in S. miltiorrhiza. We found that SmERF73 overexpression in S. miltiorrhiza coordinately up-regulated the transcription of seven tanshinone biosynthetic genes, four of which were located in the tanshinone gene cluster, consequently increasing tanshinone accumulation, while SmERF73 silencing reduced corresponding gene transcription and tanshinone accumulation. SmERF73 recognizes GCC-box promoter elements of four tanshinone-associated genes (DXR1, CPS1, KSL1 and CYP76AH3) and activates their expression. Moreover, SmERF73 and its targets were up-regulated by stress elicitors; SmERF73 appears to be at least partly mediated by the jasmonic acid (JA) signaling pathway via interaction with SmJAZ3. SmERF73 coordinately regulates tanshinone biosynthetic gene expression, suggesting a potential link between tanshinone production and plant stress responses.

MiR408-SmLAC3 Module Participates in Salvianolic Acid B Synthesis in Salvia miltiorrhiza

MicroRNAs (miRNAs) are important regulators of gene expression involved in plant development and abiotic stress responses. Recently, miRNAs have also been reported to be engaged in the regulation of secondary plant metabolism. However, there are few functional studies of miRNAs in medicinal plants. For this study, we obtained Sm-miR408 interference lines to investigate the function of Sm-miR408 in a medicinal model plant (Salvia miltiorrhiza). It was found that inhibiting the expression of Sm-miR408 could increase the content of salvianolic acid B and rosmarinic acid in the roots. The SmLAC3 and Sm-miR408 expression patterns were analyzed by qRT-PCR. A 5’ RLM-RACE assay confirmed that Sm-miR408 targets and negatively regulates SmLAC3. Moreover, the overexpression of SmLAC3 in S. miltiorrhiza promoted the accumulation of salvianolic acids in the roots. Furthermore, the lignin content of the roots in overexpressed SmLAC3 lines was decreased. Taken together, these findings indicated that Sm-miR408 modulates the accumulation of phenolic acids in S. miltiorrhiza by targeting SmLAC3 expression levels.

Late embryogenesis abundant (LEA) gene family in Salvia miltiorrhiza: identification, expression analysis, and response to drought stress

Late embryogenesis abundant (LEA) proteins play important roles in plant defense response to drought stress. However, genome-wide identification of the LEA gene family was not revealed in Salvia miltiorrhiza. In this study, 61 SmLEA genes were identified from S. miltiorrhiza and divided into 7 subfamilies according to their conserved domains and phylogenetic relationships. SmLEA genes contained the LEA conserved motifs and few introns. SmLEA genes of the same subfamilies had similar gene structures and predicted subcellular locations. Our results indicated that the promoters of SmLEA genes contained various cis-acting elements associated with abiotic stress response. In addition, RNA-seq and real-time PCR results suggested that SmLEA genes are specifically expressed in different tissue, and most SmLEA genes can be induced by drought stress. These results provide a valuable foundation for future functional investigations of SmLEA genes and drought stress-resistant breeding of S. miltiorrhiza.

Genome-wide analysis of ATP-binding cassette transporter provides insight to genes related to bioactive metabolite transportation in Salvia miltiorrhiza

BACKGROUND

ATP-binding cassette (ABC) transporters have been found to play important roles in metabolic transport in plant cells, influencing subcellular compartmentalisation and tissue distribution of these metabolic compounds. Salvia miltiorrhiza Bunge, known as Danshen in traditional Chinese medicine, is a highly valued medicinal plant used to treat cardiovascular and cerebrovascular diseases. The dry roots and rhizomes of S. miltiorrhiza contain biologically active secondary metabolites of tanshinone and salvianolic acid. Given an assembled and annotated genome and a set of transcriptome data of S. miltiorrhiza, we analysed and identified the candidate genes that likely involved in the bioactive metabolite transportation of this medicinal plant, starting with the members of the ABC transporter family.

RESULTS

A total of 114 genes encoding ABC transporters were identified in the genome of S. miltiorrhiza. All of these ABC genes were divided into eight subfamilies: 3ABCA, 31ABCB, 14ABCC, 2ABCD, 1ABCE, 7ABCF, 46ABCG, and 10 ABCI. Gene expression analysis revealed tissue-specific expression profiles of these ABC transporters. In particular, we found 18 highly expressed transporters in the roots of S. miltiorrhiza, which might be involved in transporting the bioactive compounds of this medicinal plant. We further investigated the co-expression profiling of these 18 genes with key enzyme genes involved in tanshinone and salvianolic acid biosynthetic pathways using quantitative reverse transcription polymerase chain reaction (RT-qPCR). From this RT-qPCR validation, we found that three ABC genes (SmABCG46, SmABCG40, and SmABCG4) and another gene (SmABCC1) co-expressed with the key biosynthetic enzymes of these two compounds, respectively, and thus might be involved in tanshinone and salvianolic acid transport in root cells. In addition, we predicted the biological functions of S. miltiorrhiza ABC transporters using phylogenetic relationships and analysis of the transcriptome to find biological functions.

CONCLUSIONS

Here, we present the first systematic analysis of ABC transporters in S. miltiorrhiza and predict candidate transporters involved in bioactive compound transportation in this important medicinal plant. Using genome-wide identification, transcriptome profile analysis, and phylogenetic relationships, this research provides a new perspective on the critical functions of ABC transporters in S. miltiorrhiza.

Genome-wide analysis of ATP-binding cassette transporter provides insight to genes related to bioactive metabolite transportation in Salvia miltiorrhiza

BACKGROUND

ATP-binding cassette (ABC) transporters have been found to play important roles in metabolic transport in plant cells, influencing subcellular compartmentalisation and tissue distribution of these metabolic compounds. Salvia miltiorrhiza Bunge, known as Danshen in traditional Chinese medicine, is a highly valued medicinal plant used to treat cardiovascular and cerebrovascular diseases. The dry roots and rhizomes of S. miltiorrhiza contain biologically active secondary metabolites of tanshinone and salvianolic acid. Given an assembled and annotated genome and a set of transcriptome data of S. miltiorrhiza, we analysed and identified the candidate genes that likely involved in the bioactive metabolite transportation of this medicinal plant, starting with the members of the ABC transporter family.

RESULTS

A total of 114 genes encoding ABC transporters were identified in the genome of S. miltiorrhiza. All of these ABC genes were divided into eight subfamilies: 3ABCA, 31ABCB, 14ABCC, 2ABCD, 1ABCE, 7ABCF, 46ABCG, and 10 ABCI. Gene expression analysis revealed tissue-specific expression profiles of these ABC transporters. In particular, we found 18 highly expressed transporters in the roots of S. miltiorrhiza, which might be involved in transporting the bioactive compounds of this medicinal plant. We further investigated the co-expression profiling of these 18 genes with key enzyme genes involved in tanshinone and salvianolic acid biosynthetic pathways using quantitative reverse transcription polymerase chain reaction (RT-qPCR). From this RT-qPCR validation, we found that three ABC genes (SmABCG46, SmABCG40, and SmABCG4) and another gene (SmABCC1) co-expressed with the key biosynthetic enzymes of these two compounds, respectively, and thus might be involved in tanshinone and salvianolic acid transport in root cells. In addition, we predicted the biological functions of S. miltiorrhiza ABC transporters using phylogenetic relationships and analysis of the transcriptome to find biological functions.

CONCLUSIONS

Here, we present the first systematic analysis of ABC transporters in S. miltiorrhiza and predict candidate transporters involved in bioactive compound transportation in this important medicinal plant. Using genome-wide identification, transcriptome profile analysis, and phylogenetic relationships, this research provides a new perspective on the critical functions of ABC transporters in S. miltiorrhiza.

JA-Responsive Transcription Factor SmMYB97 Promotes Phenolic Acid and Tanshinone Accumulation in Salvia miltiorrhiza

Phenolic acids and tanshinones are active principles in Salvia miltiorrhiza Bunge administered for cardiovascular and cerebrovascular diseases. Jasmonic acid (JA) promotes secondary metabolite accumulation, but the regulatory mechanism is unknown in S. miltiorrhiza. We identified and characterized the JA-responsive gene SmMYB97. Multiple sequence alignment and phylogenetic tree analyses showed that SmMYB97 was clustered with AtMYB11, AtMYB12, and ZmP1 in the subgroup S7 regulating flavonol biosynthesis. SmMYB97 was highly expressed in S. miltiorrhiza leaves and induced by methyl jasmonate (MeJA). SmMYB97 was localized in the nucleus and had strong transcriptional activation activity. SmMYB97 overexpression increased phenolic acid and tanshinone biosynthesis and upregulated the genes implicated in these processes. Yeast one-hybrid and transient transcriptional activity assays disclosed that SmMYB97 binds the PAL1, TAT1, CPS1, and KSL1 promoter regions. SmJAZ8 interacts with SmMYB97 and downregulates the genes that it controls. This study partially clarified the regulatory network of MeJA-mediated secondary metabolite biosynthesis in S. miltiorrhiza.

A high-quality reference genome sequence of Salvia miltiorrhiza provides insights into tanshinone synthesis in its red rhizomes

Salvia miltiorrhiza Bunge, also known as red sage or Danshen, is an important traditional Chinese medicine (TCM) that has been used for thousands of years to treat cardiovascular and other diseases. It is also considered an important model TCM plant. Here, a high-quality reference genome of S. miltiorrhiza was generated by combining PacBio long-read sequencing and chromatin interaction mapping (Hi-C) technologies, resulting in the chromosome-scale assembly of a 594.75-Mb genome sequence with a contig N50 of 2.70 Mb. This assembly shows the highest level of continuity for a Danshen genome generated thus far. The S. miltiorrhiza genome contained 32,483 protein-coding genes, with a repetitive DNA content of approximately 64.84%. The high percentage of young LTRs suggests that multiple TE transposition bursts occurred recently in S. miltiorrhiza. Genes unique to secondary metabolism pathways were expanded in the S. miltiorrhiza genome. A new CYP450 gene cluster was identified in the phloem of red roots where active components were synthesized. This reference genome sequence will facilitate future studies aimed at the elucidation of the secondary metabolism synthesis pathway and the genetic improvement of S. miltiorrhiza.

Integrative omic and transgenic analyses reveal the positive effect of ultraviolet-B irradiation on salvianolic acid biosynthesis through upregulation of SmNAC1

Salvianolic acids (SalAs), a group of secondary metabolites in Salvia miltiorrhiza, are widely used for treating cerebrovascular diseases. Their biosynthesis is modulated by a variety of abiotic factors, including ultraviolet-B (UV-B) irradiation; however, the underlying mechanisms remain largely unknown. Here, an integrated metabolomic, proteomic, and transcriptomic approach coupled with transgenic analyses was employed to dissect the mechanisms underlying UV-B irradiation-induced SalA biosynthesis. Results of metabolomics showed that 28 metabolites, including 12 SalAs, were elevated in leaves of UV-B-treated S. miltiorrhiza. Meanwhile, the contents of several phytohormones, including jasmonic acid and salicylic acid, which positively modulate the biosynthesis of SalAs, also increased in UV-B-treated S. miltiorrhiza. Consistently, 20 core biosynthetic enzymes and numerous transcription factors that are involved in SalA biosynthesis were elevated in treated samples as indicated by a comprehensive proteomic analysis. Correlation and gene expression analyses demonstrated that the NAC1 gene, encoding a NAC transcriptional factor, was positively involved in UV-B-induced SalA biosynthesis. Accordingly, overexpression and RNA interference of NAC1 increased and decreased SalA contents, respectively, through regulation of key biosynthetic enzymes. Furthermore, ChIP-qPCR and Dual-LUC assays showed that NAC1 could directly bind to the CATGTG and CATGTC motifs present in the promoters of the SalA biosynthesis-related genes PAL3 and TAT3, respectively, and activate their expression. Our results collectively demonstrate that NAC1 plays a crucial role in UV-B irradiation-induced SalA biosynthesis. Taken together, our findings provide mechanistic insights into the UV-B-induced SalA biosynthesis in S. miltiorrhiza, and shed light on a great potential for the development of SalA-abundant varieties through genetic engineering.

smi-miR396b targeted SmGRFs, SmHDT1, and SmMYB37/4 synergistically regulates cell growth and active ingredient accumulation in Salvia miltiorrhiza hairy roots

MIR396b had been cloned and overexpressed in Salvia miltiorrhiza hairy roots. MiR396b targets SmGRFs, SmHDT1, and SmMYB37/4 to regulate cell growth and secondary metabolism in S. miltiorrhiza hairy roots. Danshen (Salvia miltiorrhiza Bunge) is a valuable medicinal herb with two kinds of clinically used natural products, salvianolic acids and tanshinones. miR396 is a conserved microRNA and plays extensive roles in plants. However, it is still unclear how miR396 works in S. miltiorrhiza. In this study, an smi-MIR396b has been cloned from S. miltiorrhiza. Overexpression of miR396b in danshen hairy roots inhibited hairy root growth, reduced salvianolic acid concentration, but enhanced tanshinone accumulation, resulting in the biomass and total salvianolic acids respectively reduced to 55.5 and 72.1% of the control and total tanshinones increased up to 1.91-fold of the control. Applied degradome sequencing, 5'RLM-RACE, and qRT-PCR, 13 targets for miR396b were identified including seven conserved SmGRF1-7 and six novel ones. Comparative transcriptomics and microRNomics analysis together with qRT-PCR results confirmed that miR396b targets SmGRFs, SmHDT1, and SmMYB37/4 to mediate the phytohormone, especially gibberellin signaling pathways and consequentially resulted in the phenotype variation of miR396b-OE hairy roots. Furthermore, miR396b could be activated by methyl jasmonate, abscisic acid, gibberellin, salt, and drought stresses. The findings in this study indicated that smi-miR396b acts as an upstream and central regulator in cell growth and the biosynthesis of tanshinones and salvianolic acids, shedding light on the coordinated regulation of plant growth and biosynthesis of active ingredients in S. miltiorrhiza.

Systematic analysis of SmWD40s, and responding of SmWD40-170 to drought stress by regulation of ABA- and H2O2-induced stomal movement in Salvia miltiorrhiza bunge

WD40 proteins play crucial roles in response to abiotic stress. By screening the genome sequences of Salvia miltiorrhiza Bunge, 225 SmWD40 genes were identified and divided into 9 subfamilies (I-IX). Physiological, biochemical, gene structure, conserved protein motif and GO annotation analyses were performed on SmWD40 family members. The SmWD40-170 was found in 110 SmWD40 genes that contain drought response elements, SmWD40-170 was one of these genes whose response in terms of expression under drought was significant. The expression of SmWD40-170 was also up-regulated by ABA and H₂O₂. Through observed the stomatal phenotype of SmWD40-170 transgenic lines, the stomatal closure was abolished under dehydration, ABA and H₂O₂ treatment in SmWD40-170 knockdown lines. Abscisic acid (ABA), as the key phytohormone, elevates reactive oxygen species (ROS) levels under drought stress. The ABA-ROS interaction mediated the generation of H₂O₂ and the activation of anion channel in guard cells. The osmolality alteration of guard cells further accelerated the stomatal closure. As a second messenger, nitric oxide (NO) regulated ABA signaling, the NO stimulated protein kinase activity inhibited the K⁺ influx which result in stomatal closure. These NO-relevant events were essential for ABA-induced stomatal closure. The reduction of NO production was also observed in the guard cells of SmWD40-170 knockdown lines. The abolished of stomatal closure attributed to the SmWD40-170 deficiency induced the reduction of NO content. In general, the SmWD40-170 is a critical drought response gene in SmWD40 gene family and regulates ABA- and H₂O₂-induced stomatal movement by affecting the synthesis of NO.

Integrated metabolomic and transcriptomic analysis of the anthocyanin regulatory networks in Salvia miltiorrhiza Bge. flowers

BACKGROUND

The objectives of this study were to reveal the anthocyanin biosynthesis metabolic pathway in white and purple flowers of Salvia miltiorrhiza using metabolomics and transcriptomics, to identify different anthocyanin metabolites, and to analyze the differentially expressed genes involved in anthocyanin biosynthesis.

RESULTS

We analyzed the metabolomics and transcriptomics data of S. miltiorrhiza flowers. A total of 1994 differentially expressed genes and 84 flavonoid metabolites were identified between the white and purple flowers of S. miltiorrhiza. Integrated analysis of transcriptomics and metabolomics showed that cyanidin 3,5-O-diglucoside, malvidin 3,5-diglucoside, and cyanidin 3-O-galactoside were mainly responsible for the purple flower color of S. miltiorrhiza. A total of 100 unigenes encoding 10 enzymes were identified as candidate genes involved in anthocyanin biosynthesis in S. miltiorrhiza flowers. Low expression of the ANS gene decreased the anthocyanin content but enhanced the accumulation of flavonoids in S. miltiorrhiza flowers.

CONCLUSIONS

Our results provide valuable information on the anthocyanin metabolites and the candidate genes involved in the anthocyanin biosynthesis pathways in S. miltiorrhiza.

Systematic Analysis of Kelch Repeat F-box (KFB) Protein Gene Family and Identification of Phenolic Acid Regulation Members in Salvia miltiorrhiza Bunge

S. miltiorrhiza is a well-known Chinese herb for the clinical treatment of cardiovascular and cerebrovascular diseases. Tanshinones and phenolic acids are the major secondary metabolites and significant pharmacological constituents of this plant. Kelch repeat F-box (KFB) proteins play important roles in plant secondary metabolism, but their regulation mechanism in S. miltiorrhiza has not been characterized. In this study, we systematically characterized the S. miltiorrhiza KFB gene family. In total, 31 SmKFB genes were isolated from S. miltiorrhiza. Phylogenetic analysis of those SmKFBs indicated that 31 SmKFBs can be divided into four groups. Thereinto, five SmKFBs (SmKFB1, 2, 3, 5, and 28) shared high homology with other plant KFBs which have been described to be regulators of secondary metabolism. The expression profile of SmKFBs under methyl jasmonate (MeJA) treatment deciphered that six SmKFBs (SmKFB1, 2, 5, 6, 11, and 15) were significantly downregulated, and two SmKFBs (SmKFB22 and 31) were significantly upregulated. Tissue-specific expression analysis found that four SmKFBs (SmKFB4, 11, 16, and 17) were expressed preferentially in aerial tissues, while two SmKFBs (SmKFB5, 25) were predominantly expressed in roots. Through a systematic analysis, we speculated that SmKFB1, 2, and 5 are potentially involved in phenolic acids biosynthesis.

The cytosolic protein GRP1 facilitates abscisic acid- and darkness-induced stomatal closure in Salvia miltiorrhiza

By screening an expressed sequence tag (EST) library of Salvia miltiorrhiza, we detected an acidic protein, SmGRP1, with no significant similarities to the other sequences in public databases. SmGRP1 encodes a peptide of 151 amino acids, 33.77 % of which are glutamic acid residues, and the peptide was positive according to "stains-all" staining. Expression analysis revealed that SmGRP1 was expressed in all examined tissues of S. miltiorrhiza but was most highly expressed in the leaves and stems. Without a signal peptide, SmGRP1 localized to the cytoplasm in protoplasts in subcellular localization experiments. SmGRP1 expression was prominently enhanced by ABA and darkness treatments; the protein could also be induced by high temperature, NaCl, and dehydration treatments, while low temperature suppressed its expression. Furthermore, although there were no obvious phenotypic differences in SmGRP1 overexpression and SmGRP1 knockdown mutants compared with control plants under normal culture conditions, the stomata of the knockdown lines remained open when treated with ABA, darkness, NO, and H₂O₂. In addition, the water loss rate of the knockdown mutants was faster than that of the control lines and overexpression mutants when exposed to air. These observations indicate that SmGRP1 is a novel acidic protein with potential calcium-binding capability and is involved in stomatal movement and stress resistance.

Salvia miltiorrhiza-derived miRNAs suppress vascular remodeling through regulating OTUD7B/KLF4/NMHC IIA axis

Objective: Abnormal proliferation and migration of vascular smooth muscle cells (VSMCs) are essential for vascular remodeling. Natural compounds with diterpene chinone or phenolic acid structure from Salvia miltiorrhiza, an eminent medicinal herb widely used to treat cardiovascular diseases in China, can effectively attenuate vascular remodeling induced by vascular injury. However, it remains unknown whether Salvia miltiorrhiza-derived miRNAs can protect VSMCs from injury by environmental stimuli. Here, we explored the role and underlying mechanisms of Salvia miltiorrhiza-derived Sal-miR-1 and 3 in the regulation of VSMC migration and monocyte adhesion to VSMCs induced by thrombin. Methods: A mouse model for intimal hyperplasia was established by the ligation of carotid artery and the injured carotid arteries were in situ-transfected with Sal-miR-1 and 3 using F-127 pluronic gel. The vascular protective effects of Sal-miR-1 and 3 were assessed via analysis of intimal hyperplasia with pathological morphology. VSMC migration and adhesion were analyzed by the wound healing, transwell membrane assays, and time-lapse imaging experiment. Using loss- and gain-of-function approaches, Sal-miR-1 and 3 regulation of OTUD7B/KLF4/NMHC IIA axis was investigated by using luciferase assay, co-immunoprecipitation, chromatin immunoprecipitation, western blotting, etc. Results:Salvia miltiorrhiza-derived Sal-miR-1 and 3 can enter the mouse body after intragastric administration, and significantly suppress intimal hyperplasia induced by carotid artery ligation. In cultured VSMCs, these two miRNAs inhibit thrombin-induced the migration of VSMCs and monocyte adhesion to VSMCs. Mechanistically, Sal-miR-1 and 3 abrogate OTUD7B upregulation by thrombin via binding to the different sites of the OTUD7B 3'UTR. Most importantly, OTUD7B downregulation by Sal-miR-1 and 3 attenuates KLF4 protein levels via decreasing its deubiquitylation, whereas decreased KLF4 relieves its repression of transcription of NMHC IIA gene and thus increases NMHC IIA expression levels. Further, increased NMHC IIA represses VSMC migration and monocyte adhesion to VSMCs via maintaining the contractile phenotype of VSMCs. Conclusions: Our studies not only found the novel bioactive components from Salvia miltiorrhiza but also clarified the molecular mechanism underlying Sal-miR-1 and 3 inhibition of VSMC migration and monocyte adhesion to VSMCs. These results add important knowledge to the pharmacological actions and bioactive components of Salvia miltiorrhiza. Sal-miR-1 and 3-regulated OTUD7B/KLF4/NMHC IIA axis may represent a therapeutic target for vascular remodeling.

The ethylene response factor SmERF8 regulates the expression of SmKSL1 and is involved in tanshinone biosynthesis in Saliva miltiorrhiza hairy roots

Saliva miltiorrhiza ethylene response factor (SmERF), predicted to be expressed genome-wide, is the potential regulator of tanshinone biosynthesis. However, few studies have investigated its transcriptional regulation pathways in tanshinone biosynthesis. Here, we report an ethylene response factor (SmERF8), which was screened by the SmKSL1 (a key gene in tanshinone biosynthesis) promoter from the S. miltiorrhiza cDNA library. The SmERF8, highly expressed in S. miltiorrhiza root head, is sensitive to Eth stress, and its protein was enriched in the nucleus. The SmERF8 recognizes the GCC-box in the SmKSL1 promoter. Overexpression and RNAi of SmERF8 in S. miltiorrhiza transgenic hairy roots showed that the tanshinone contents were significantly increased in the overexpression transgenic lines and decreased in RNAi lines. These results suggest that the SmERF8 may be a central activator that regulates the expression of SmKSL1 by binding the GCC-box and then promoting tanshinone biosynthesis. Thus, the SmERF8 may functionally accelerate tanshinone biosynthesis by the transcriptional regulation of its key gene.

SmPPT, a 4-hydroxybenzoate polyprenyl diphosphate transferase gene involved in ubiquinone biosynthesis, confers salt tolerance in Salvia miltiorrhiza

SmPPT, which encodes 4-hydroxybenzoate polyprenyl diphosphate transferase involved in ubiquinone biosynthesis, confers salt tolerance to S. miltiorrhiza through enhancing the activities of POD and CAT to scavenge ROS. Ubiquinone (UQ), also known as coenzyme Q (CoQ), is a key electron transporter in the mitochondrial respiratory system. UQ is composed of a benzene quinone ring and a polyisoprenoid side chain. Attachment of polyisoprenoid side chain to the benzene quinone ring is a rate-limiting step catalyzed by 4-hydroxybenzoate polyprenyl diphosphate transferase (PPT). So far, only a few plant PPT-encoding genes have been functionally analyzed. Through genome-wide analysis and subsequent molecular cloning, a PPT-encoding gene, termed SmPPT, was identified from an economically and academically important medicinal model plant, Salvia miltiorrhiza. SmPPT contained many putative cis-elements associated with abiotic stresses in the promoter region and were responsive to PEG-6000 and methyl jasmonate treatments. The deduced SmPPT protein contains the PT_UbiA conserved domain of polyprenyl diphosphate transferase and an N-terminal mitochondria transit peptide. Transient expression assay of SmPPT-GFP fusion protein showed that SmPPT was mainly localized in the mitochondria. SmPPT could functionally complement coq2 mutation and catalyzed UQ6 production in yeast cells. Overexpression of SmPPT increased UQ production and enhanced salt tolerance in S. miltiorrhiza. Under salinity stress conditions, transgenic plants accumulated less H₂O₂ and malondialdehyde and exhibited higher peroxidase (POD) and catalase (CAT) activities compared with wild-type plants. It indicates that SmPPT confers salt tolerance to S. miltiorrhiza at least partially through enhancing the activities of POD and CAT to scavenge ROS.

Transcriptome sequencing of Salvia miltiorrhiza after infection by its endophytic fungi and identification of genes related to tanshinone biosynthesis

Context: Salvia miltiorrhiza Bunge (Labiatae) is a traditional Chinese herb. Endophytic fungi, which are biotic elicitors, can induce accumulation of secondary metabolites in their host plants.Objective: To analyze the interaction mechanism between S. miltiorrhiza and endophytic fungi.Materials and methods: Endophytic fungi U104 producing tanshinone IIA were isolated from the healthy disease-free tissue of root of S. miltiorrhiza by conventional methods. The endophytic fungus U104 of S. miltiorrhiza was co-cultured with the sterile seedlings of S. miltiorrhiza for 20 d (temp:day/night = 26 °C/18 °C, photoperiod:12/12 h, illuminance:2000 Lx). Transcriptome sequencing of S. miltiorrhiza seedlings after 20 d of co-cultivation was performed using the Illumina platform.Results: A total of 3713 differentially expressed genes (DEGs) were obtained. These different expression genes, such as STPII, LTP2, MYB transcription factors, CNGC, CDPK, Rboh, CaM, MAP2K1/MEK1, WRKY33, SGT1/SGT and Hsp90/htpG, showed that host S. miltiorrhiza had biological defence response in the initial stage of interaction. Under the induction of endophytic fungi, 14 key enzyme genes were up-regulated in the tanshinone biosynthesis pathway: DXS, DXS2, DXR, HMGR3, AACT, MK, PMK, GGPPS2, GPPS, KSL, IDI, IPII, FDPS and CPS.Discussion and conclusions: A total of 14 key genes were obtained from the tanshinone component synthesis and metabolic pathways, providing a reasonable explanation for the accumulation of tanshinone components, an accumulation induced by endophytic fungi, in the host plants. The large amounts of data generated in this study provide a strong and powerful platform for future functional and molecular studies of interactions between host plants and their endophytic fungi.

[Cloning,subcellular localization and expression pattern analysis of transcription factor SmGRAS3 from Salvia miltiorrhiza]

GRAS transcription factors play important roles in the regulation of plant root growth and GA signaling. In this study,SmGRAS3 gene was cloned,which open reading frame was 2 247 bp,and encoding 748 amino acids. The physicochemical properties and structure of SmGRAS3 and its encoded protein were analyzed by bioinformatics software. This gene belongs to the SCL9 subfamily of the GRAS family,and its promoter sequence mainly contains the light response,stress response,and hormone response elements. It may interact with the GA signal pathway and anti-stress related proteins. The subcellular localization showed that SmGRAS3 protein was mainly located in the nucleus. The expression pattern analysis showed that the expression of Sm GRAS3 was the highest in the root and the lowest in the stem,and both light and low temperature could induce the high expression level of SmGRAS3. This study provides a foundation for further study on the roles of SmGRAS3 gene in the root growth and stress tolerance of Salvia miltiorrhiza.

The AP2/ERF transcription factor SmERF128 positively regulates diterpenoid biosynthesis in Salvia miltiorrhiza

The novel AP2/ERF transcription factor SmERF128 positively regulates diterpenoid tanshinone biosynthesis by activating the expression of SmCPS1, SmKSL1, and SmCYP76AH1 in Salvia miltiorrhiza. Certain members of the APETALA2/ethylene-responsive factor (AP2/ERF) family regulate plant secondary metabolism. Although it is clearly documented that AP2/ERF transcription factors (TFs) are involved in sesquiterpenoid biosynthesis, the regulation of diterpenoid biosynthesis by AP2/ERF TFs remains elusive. Here, we report that the novel AP2/ERF TF SmERF128 positively regulates diterpenoid tanshinone biosynthesis in Salvia miltiorrhiza. Overexpression of SmERF128 increased the expression levels of copalyl diphosphate synthase 1 (SmCPS1), kaurene synthase-like 1 (SmKSL1) and cytochrome P450 monooxygenase 76AH1 (SmCYP76AH1), whereas their expression levels were decreased when SmERF128 was silenced. Accordingly, the content of tanshinone was reduced in SmERF128 RNA interference (RNAi) hairy roots and dramatically increased in SmERF128 overexpression hairy roots, as demonstrated through Ultra Performance Liquid Chromatography (UPLC) and Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) analysis. Furthermore, SmERF128 activated the expression of SmCPS1, SmKSL1, and SmCYP76AH1 by binding to the GCC box, and to the CRTDREHVCBF2 (CBF2) and RAV1AAT (RAA) motifs within their promoters during in vivo and in vitro assays. Our findings not only reveal the molecular basis of how the AP2/ERF transcription factor SmERF128 regulates diterpenoid biosynthesis, but also provide useful information for improving tanshinone production through genetic engineering.

Heterologous Expression of Salvia miltiorrhiza MicroRNA408 Enhances Tolerance to Salt Stress in Nicotiana benthamiana

MicroRNAs (miRNAs) are a class of endogenous small RNAs that regulate the expression of target genes post-transcriptionally; they are known to play major roles in development and responses to abiotic stress. MicroRNA408 (miR408) is a conserved small RNA in plants; it was reported that miR408 genes were involved in abiotic stress in Arabidopsis. However, miR408 in Salvia miltiorrhiza has been rarely investigated. In this study, we cloned Sm-MIR408, the miR408 precursor sequence, and its promoter sequence from S. miltiorrhiza and the role in tolerance to salt stress is described. The effects of salt stress on miR408 expression were studied by using β-glucuronidase (GUS) staining. Our data indicated that transgenic tobacco overexpressing Sm-MIR408 promoted seed germination and reduced the accumulation of reactive oxygen species under salt stress. Transcript levels of antioxidative genes, i.e., NbSOD, NbPOD, and NbCAT, and their enzyme activities increased in salinity-stressed transgenic tobacco plants, suggesting a better antioxidant system to cope the oxidative damage caused by salinity stress. Taken together, these findings indicated that miR408 functions in positive responses to salt tolerance in tobacco.

Overexpression of SmbHLH148 induced biosynthesis of tanshinones as well as phenolic acids in Salvia miltiorrhiza hairy roots

SmbHLH148 activated the whole biosynthetic pathways of phenolic acids and tanshinones, thus upregulated the production of both the two groups of pharmaceutical ingredients in Salvia miltiorrhiza. Phenolic acids and tanshinones are the two important groups of pharmaceutical ingredients presented in Salvia miltiorrhiza Bunge. The bHLH transcription factors could regulate secondary metabolism efficiently in plants. However, there are only some MYCs have been studied on regulation of either phenolic acids or tanshinones biosynthesis. In this study, a bHLH TF named SmbHLH148, which is homologous to AtbHLH148, AtbHLH147 and CubHLH1, was isolated and functionally characterized from S. miltiorrhiza. Transcription of SmbHLH148 could be intensely induced by ABA and also be moderately induced by MeJA and GA. SmbHLH148 is present in all the six tissues and mostly expressed in fibrous root and flowers. Subcellular localization analysis found that SmbHLH148 was localized in the nucleus. Overexpression of SmbHLH148 significantly increased not only three phenolic acids components accumulation but also three tanshinones content. Content of caffeic acid, rosmarinic acid and salvianolic acid B were reached to 2.87-, 4.00- and 5.99-fold of the control in the ObHLH148-3, respectively. Content of dihydrotanshinone I, cryptotanshinone, and tanshinone I were also present highest in ObHLH148-3, reached 2.5-, 5.04- and 3.97-fold of the control, respectively. Expression analysis of pathway genes of phenolic acids and tanshinones in transgenic lines showed that most of them were obviously upregulated. Moreover, transcription of AREB and JAZs were also induced in SmbHLH148 overexpression lines. These results suggested that SmbHLH148 might be taken part in ABA and MeJA signaling and activated almost the whole biosynthetic pathways of phenolic acids and tanshinones, thus the production of phenolic acids and tanshinones were upregulated.