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.