Genome-wide identification and characterization of novel genes involved in terpenoid biosynthesis in Salvia miltiorrhiza

Biosynthesis and signal transduction of plant growth regulators and their effects on bioactive compound production in Salvia miltiorrhiza (Danshen)

Plant growth regulators (PGRs) are involved in multiple aspects of plant life, including plant growth, development, and response to environmental stimuli. They are also vital for the formation of secondary metabolites in various plants. Salvia miltiorrhiza is a famous herbal medicine and has been used commonly for > 2000 years in China, as well as widely used in many other countries. S. miltiorrhiza is extensively used to treat cardiovascular and cerebrovascular diseases in clinical practices and has specific merit against various diseases. Owing to its outstanding medicinal and commercial potential, S. miltiorrhiza has been extensively investigated as an ideal model system for medicinal plant biology. Tanshinones and phenolic acids are primary pharmacological constituents of S. miltiorrhiza. As the growing market for S. miltiorrhiza, the enhancement of its bioactive compounds has become a research hotspot. S. miltiorrhiza exhibits a significant response to various PGRs in the production of phenolic acids and tanshinones. Here, we briefly review the biosynthesis and signal transduction of PGRs in plants. The effects and mechanisms of PGRs on bioactive compound production in S. miltiorrhiza are systematically summarized and future research is discussed. This article provides a scientific basis for further research, cultivation, and metabolic engineering in S. miltiorrhiza.

The methyl jasmonate-responsive transcription factor SmERF106 promotes tanshinone accumulation in Salvia miltiorrhiza

Understanding the regulatory biosynthesis mechanisms of active compounds in herbs is vital for the preservation and sustainable use of natural medicine resources. Diterpenoids, which play a key role in plant growth and resistance, also serve as practical products for humans. Tanshinone, a class of abietane-type diterpenes unique to the Salvia genus, such as Salvia miltiorrhiza, is an excellent model for studying diterpenoids. In this study, we discovered that a transcription factor, SmERF106, responds to MeJA induction and is located in the nucleus. It exhibits a positive correlation with the expression of SmKSL1 and SmIDI1, which are associated with tanshinone biosynthesis. We performed DNA affinity purification sequencing (DAP-seq) to predict genes that may be transcriptionally regulated by SmERF106. Our cis-elements analysis suggested that SmERF106 might bind to GCC-boxes in the promoters of SmKSL1 and SmIDI1. This indicates that SmKSL1 and SmIDI1 could be potential target genes regulated by SmERF106 in the tanshinone biosynthesis pathway. Their interaction was then demonstrated through a series of in vitro and in vivo binding experiments, including Y1H, EMSA, and Dual-LUC. Overexpression of SmERF106 in the hairy root of S. miltiorrhiza led to a significant increase in tanshinone content and the transcriptional levels of SmKSL1 and SmIDI1. In summary, we found that SmERF106 can activate the transcription of SmKSL1 and SmIDI1 in response to MeJA induction, thereby promoting tanshinone biosynthesis. This discovery provides new insights into the regulatory mechanisms of tanshinones in response to JA and offers a potential gene tool for tanshinone metabolic engineering strategy.

Metabolic engineering of terpene metabolism in lavender

Background

Several members of the Lamiaceae family of plants produce large amounts of essential oil [EO] that find extensive applications in the food, cosmetics, personal hygiene, and alternative medicine industries. There is interest in enhancing EO metabolism in these plants.

Main body

Lavender produces a valuable EO that is highly enriched in monoterpenes, the C10 class of the isoprenoids or terpenoids. In recent years, substantial effort has been made by researchers to study terpene metabolism and enhance lavender EO through plant biotechnology. This paper reviews recent advances related to the cloning of lavender monoterpene biosynthetic genes and metabolic engineering attempts aimed at improving the production of lavender monoterpenes in plants and microbes.

Conclusion

Metabolic engineering has led to the improvement of EO quality and yield in several plants, including lavender. Furthermore, several biologically active EO constituents have been produced in microorganisms.

Plant terpenoid biosynthetic network and its multiple layers of regulation

Terpenoids constitute one of the largest and most chemically diverse classes of primary and secondary metabolites in nature with an exceptional breadth of functional roles in plants. Biosynthesis of all terpenoids begins with the universal five‑carbon building blocks, isopentenyl diphosphate (IPP) and its allylic isomer dimethylallyl diphosphate (DMAPP), which in plants are derived from two compartmentally separated but metabolically crosstalking routes, the mevalonic acid (MVA) and methylerythritol phosphate (MEP) pathways. Here, we review the current knowledge on the terpenoid precursor pathways and highlight the critical hidden constraints as well as multiple regulatory mechanisms that coordinate and homeostatically govern carbon flux through the terpenoid biosynthetic network in plants.

Changes in secondary metabolites contents and stress responses in Salvia miltiorrhiza via ScWRKY35 overexpression: Insights from a wild relative Salvia castanea

Salvia castanea Diels, a close wild relative to the medicinal plant, Salvia miltiorrhiza Bunge, primarily grows in high-altitude regions. While the two species share similar active compounds, their content varies significantly. WRKY transcription factors are key proteins, which regulate plant growth, stress response, and secondary metabolism. We identified 46 ScWRKY genes in S. castanea and found that ScWRKY35 was a highly expressed gene associated with secondary metabolites accumulation. This study aimed to explore the role of ScWRKY35 gene in regulating the accumulation of secondary metabolites and its response to UV and cadmium (Cd) exposure in S. miltiorrhiza. It was found that transgenic S. miltiorrhiza hairy roots overexpressing ScWRKY35 displayed upregulated expression of genes related to phenolic acid synthesis, resulting in increased salvianolic acid B (SAB) and rosmarinic acid (RA) contents. Conversely, tanshinone pathway gene expression decreased, leading to lower tanshinone levels. Further, overexpression of ScWRKY35 upregulated Cd transport protein HMA3 in root tissues inducing Cd sequestration. In contrast, the Cd uptake gene NRAMP1 was downregulated, reducing Cd absorption. In response to UV radiation, ScWRKY35 overexpression led to an increase in the accumulation of phenolic acid and tanshinone contents, including upregulation of genes associated with salicylic acid (SA) and jasmonic acid (JA) synthesis. Altogether, these findings highlight the role of ScWRKY35 in enhancing secondary metabolites accumulation, as well as in Cd and UV stress modulation in S. miltiorrhiza, which offers a novel insight into its phytochemistry and provides a new option for the genetic improvement of the plants.

Physiological and molecular mechanisms of ZnO quantum dots mitigating cadmium stress in Salvia miltiorrhiza

This study delved into the physiological and molecular mechanisms underlying the mitigation of cadmium (Cd) stress in the model medicinal plant Salvia miltiorrhiza through the application of ZnO quantum dots (ZnO QDs, 3.84 nm). A pot experiment was conducted, wherein S. miltiorrhiza was subjected to Cd stress for six weeks with foliar application of 100 mg/L ZnO QDs. Physiological analyses demonstrated that compared to Cd stress alone, ZnO QDs improved biomass, reduced Cd accumulation, increased the content of photosynthetic pigments (chlorophyll and carotenoids), and enhanced the levels of essential nutrient elements (Ca, Mn, and Cu) under Cd stress. Furthermore, ZnO QDs significantly lowered Cd-induced reactive oxygen species (ROS) content, including H2O2, O2-, and MDA, while enhancing the activity of antioxidant enzymes (SOD, POD, APX, and GSH-PX). Additionally, ZnO QDs promoted the biosynthesis of primary and secondary metabolites, such as total protein, soluble sugars, terpenoids, and phenols, thereby mitigating Cd stress in S. miltiorrhiza. At the molecular level, ZnO QDs were found to activate the expression of stress signal transduction-related genes, subsequently regulating the expression of downstream target genes associated with metal transport, cell wall synthesis, and secondary metabolite synthesis via transcription factors. This activation mechanism contributed to enhancing Cd tolerance in S. miltiorrhiza. In summary, these findings shed light on the mechanisms underlying the mitigation of Cd stress by ZnO QDs, offering a potential nanomaterial-based strategy for enhancing Cd tolerance in medicinal plants.

Genome-Wide Identification and Characterization of the PHT1 Gene Family and Its Response to Mycorrhizal Symbiosis in Salvia miltiorrhiza under Phosphate Stress

Phosphorus (P) is a vital nutrient element that is essential for plant growth and development, and arbuscular mycorrhizal fungi (AMF) can significantly enhance P absorption. The phosphate transporter protein 1 (PHT1) family mediates the uptake of P in plants. However, the PHT1 gene has not yet been characterized in Salvia miltiorrhiza. In this study, to gain insight into the functional divergence of PHT1 genes, nine SmPHT1 genes were identified in the S. miltiorrhiza genome database via bioinformatics tools. Phylogenetic analysis revealed that the PHT1 proteins of S. miltiorrhiza, Arabidopsis thaliana, and Oryza sativa could be divided into three groups. PHT1 in the same clade has a similar gene structure and motif, suggesting that the features of each clade are relatively conserved. Further tissue expression analysis revealed that SmPHT1 was expressed mainly in the roots and stems. In addition, phenotypic changes, P content, and PHT1 gene expression were analyzed in S. miltiorrhiza plants inoculated with AMF under different P conditions (0 mM, 0.1 mM, and 10 mM). P stress and AMF significantly affected the growth and P accumulation of S. miltiorrhiza. SmPHT1;6 was strongly expressed in the roots colonized by AMF, implying that SmPHT1;6 was a specific AMF-inducible PHT1. Taken together, these results provide new insights into the functional divergence and genetic redundancy of the PHT1 genes in response to P stress and AMF symbiosis in S. miltiorrhiza.

Chromosome-level genome assembly of Prunella vulgaris L. provides insights into pentacyclic triterpenoid biosynthesis

Prunella vulgaris is one of the bestselling and widely used medicinal herbs. It is recorded as an ace medicine for cleansing and protecting the liver in Chinese Pharmacopoeia and has been used as the main constitutions of many herbal tea formulas in China for centuries. It is also a traditional folk medicine in Europe and other countries of Asia. Pentacyclic triterpenoids are a major class of bioactive compounds produced in P. vulgaris. However, their biosynthetic mechanism remains to be elucidated. Here, we report a chromosome-level reference genome of P. vulgaris using an approach combining Illumina, ONT, and Hi-C technologies. It is 671.95 Mb in size with a scaffold N50 of 49.10 Mb and a complete BUSCO of 98.45%. About 98.31% of the sequence was anchored into 14 pseudochromosomes. Comparative genome analysis revealed a recent WGD in P. vulgaris. Genome-wide analysis identified 35 932 protein-coding genes (PCGs), of which 59 encode enzymes involved in 2,3-oxidosqualene biosynthesis. In addition, 10 PvOSC, 358 PvCYP, and 177 PvUGT genes were identified, of which five PvOSCs, 25 PvCYPs, and 9 PvUGTs were predicted to be involved in the biosynthesis of pentacyclic triterpenoids. Biochemical activity assay of PvOSC2, PvOSC4, and PvOSC6 recombinant proteins showed that they were mixed amyrin synthase (MAS), lupeol synthase (LUS), and β-amyrin synthase (BAS), respectively. The results provide a solid foundation for further elucidating the biosynthetic mechanism of pentacyclic triterpenoids in P. vulgaris.

ZnO quantum dots alleviate salt stress in Salvia miltiorrhiza by enhancing growth, scavenging reactive oxygen species, and modulating stress-responsive genes

Experiments were conducted to investigate the alleviating effects of ZnO quantum dots (ZnO QDs) on salt stress in Salvia miltiorrhiza by comparing them with conventional ZnO nanoparticles (ZnO NPs). The results demonstrated that compared with salt stress alone, foliar application of ZnO QDs significantly improved the biomass as well as the total chlorophyll and carotenoids contents under salt stress. ZnO QDs reduced H2O2 and MDA levels, decreased non-enzymatic antioxidant (ASA and GSH) content, and improved antioxidant enzyme (POD, SOD, CAT, PAL, and PPO) activity under salt stress. Metal elemental analysis further demonstrated that the ZnO QDs markedly increased Zn and K contents while decreasing Na content, resulting in a lower Na/K ratio compared to salt stress alone. Finally, RNA sequencing results indicated that ZnO QDs primarily regulated genes associated with stress-responsive pathways, including plant hormone signal transduction, the MAPK signaling pathway, and metabolic-related pathways, thereby alleviating the adverse effects of salt stress. In comparison, ZnO NPs did not exhibit similar effects in terms of improving plant growth, enhancing the antioxidant system, or regulating stress-responsive genes under salt stress. These findings highlight the distinct advantages of ZnO QDs and suggest their potential as a valuable tool for mitigating salt stress in plants.

Integrating Transcriptome and Chemical Analyses to Provide Insights into Biosynthesis of Terpenoids and Flavonoids in the Medicinal Industrial Crop Andrographis paniculate and Its Antiviral Medicinal Parts

Andrographis paniculata is a medicinal plant traditionally used to produce diterpene lactones and flavonoids, which possess various biological activities. Widely distributed in China, India, and other Southeast Asia countries, A. paniculata has become an important economic crop, significantly treating SARS-CoV-2, and is being cultivated on a large scale in southern China. The biosynthesis of active ingredients in A. paniculata are regulated and controlled by genes, but their specific roles are still not fully understood. To further explore the growth regulation factors and utilization of its medicinal parts of this industrial crop, chemical and transcriptome analyses were conducted on the roots, stems, and leaves of A. paniculata to identify the biosynthesis pathways and related candidate genes of the active ingredients. The chemical analysis revealed that the main components of A. paniculata were diterpene lactones and flavonoids, which displayed potential ability to treat SARS-CoV-2 through molecular docking. Moreover, the transcriptome sequencing annotated a total of 40,850 unigenes, including 7962 differentially expressed genes. Among these, 120 genes were involved in diterpene lactone biosynthesis and 60 genes were involved in flavonoid biosynthesis. The expression of diterpene lactone-related genes was the highest in leaves and the lowest in roots, consistent with our content determination results. It is speculated that these highly expressed genes in leaves may be involved in the biosynthesis pathway of diterpenes. Furthermore, two class Ⅰ terpene synthases in A. paniculata transcriptome were also annotated, providing reference for the downstream pathway of the diterpene lactone biosynthesis. With their excellent market value, our experiments will promote the study of the biosynthetic genes for active ingredients in A. paniculata and provide insights for subsequent in vitro biosynthesis.

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.

Antioxidants of Non-Enzymatic Nature: Their Function in Higher Plant Cells and the Ways of Boosting Their Biosynthesis

Plants are exposed to a variety of abiotic and biotic stresses leading to increased formation of reactive oxygen species (ROS) in plant cells. ROS are capable of oxidizing proteins, pigments, lipids, nucleic acids, and other cell molecules, disrupting their functional activity. During the process of evolution, numerous antioxidant systems were formed in plants, including antioxidant enzymes and low molecular weight non-enzymatic antioxidants. Antioxidant systems perform neutralization of ROS and therefore prevent oxidative damage of cell components. In the present review, we focus on the biosynthesis of non-enzymatic antioxidants in higher plants cells such as ascorbic acid (vitamin C), glutathione, flavonoids, isoprenoids, carotenoids, tocopherol (vitamin E), ubiquinone, and plastoquinone. Their functioning and their reactivity with respect to individual ROS will be described. This review is also devoted to the modern genetic engineering methods, which are widely used to change the quantitative and qualitative content of the non-enzymatic antioxidants in cultivated plants. These methods allow various plant lines with given properties to be obtained in a rather short time. The most successful approaches for plant transgenesis and plant genome editing for the enhancement of biosynthesis and the content of these antioxidants are discussed.

Tanshinone IIA targeting cell signaling pathways: a plausible paradigm for cancer therapy

Natural compounds originating from plants offer a wide range of pharmacological potential and have traditionally been used to treat a wide range of diseases including cancer. Tanshinone IIA (Tan IIA), a bioactive molecule found in the roots of the Traditional Chinese Medicine (TCM) herb Salvia miltiorrhiza, has been shown to have remarkable anticancer properties through several mechanisms, such as inhibition of tumor cell growth and proliferation, metastasis, invasion, and angiogenesis, as well as induction of apoptosis and autophagy. It has demonstrated excellent anticancer efficacy against cell lines from breast, cervical, colorectal, gastric, lung, and prostate cancer by modulating multiple signaling pathways including PI3K/Akt, JAK/STAT, IGF-1R, and Bcl-2-Caspase pathways. This review focuses on the role of Tan IIA in the treatment of various cancers, as well as the underlying molecular mechanisms.

Characteristics of Salvia miltiorrhiza methylome and the regulatory mechanism of DNA methylation in tanshinone biosynthesis

Salvia miltiorrhiza is a model medicinal plant with significant economic and medicinal value. Its roots produce a group of diterpenoid lipophilic bioactive components, termed tanshinones. Biosynthesis and regulation of tanshinones has attracted widespread interest. However, the methylome of S. miltiorrhiza has not been analysed and the regulatory mechanism of DNA methylation in tanshinone production is largely unknown. Here we report single-base resolution DNA methylomes from roots and leaves. Comparative analysis revealed differential methylation patterns for CG, CHG, and CHH contexts and the association between DNA methylation and the expression of genes and small RNAs. Lowly methylated genes always had higher expression levels and 24-nucleotide sRNAs could be key players in the RdDM pathway in S. miltiorrhiza. DNA methylation variation analysis showed that CHH methylation contributed mostly to the difference. Go enrichment analysis showed that diterpenoid biosynthetic process was significantly enriched for genes with downstream overlapping with hypoCHHDMR in July_root when comparing with those in March_root. Tanshinone biosynthesis-related enzyme genes, such as DXS2, CMK, IDI1, HMGR2, DXR, MDS, CYP76AH1, 2OGD25, and CYP71D373, were less CHH methylated in gene promoters or downstream regions in roots collected in July than those collected in March. Consistently, gene expression was up-regulated in S. miltiorrhiza roots collected in July compared with March and the treatment of DNA methylation inhibitor 5-azacytidine significantly promoted tanshinone production. It suggests that DNA methylation plays a significant regulatory role in tanshinone biosynthesis in S. miltiorrhiza through changing the levels of CHH methylation in promoters or downstreams of key enzyme genes.

Chromosome-level genome assembly of Salvia miltiorrhiza with orange roots uncovers the role of Sm2OGD3 in catalyzing 15,16-dehydrogenation of tanshinones

Salvia miltiorrhiza is well known for its clinical practice in treating heart and cardiovascular diseases. Its roots, used for traditional Chinese medicine materials, are usually brick-red due to accumulation of red pigments, such as tanshinone IIA and tanshinone I. Here we report a S. miltiorrhiza line (shh) with orange roots. Compared with the red roots of normal S. miltiorrhiza plants, the contents of tanshinones with a single bond at C-15,16 were increased, whereas those with a double bond at C-15,16 were significantly decreased in shh. We assembled a high-quality chromosome-level genome of shh. Phylogenomic analysis showed that the relationship between two S. miltiorrhiza lines with red roots was closer than the relationship with shh. It indicates that shh could not be the mutant of an extant S. miltiorrhiza line with red roots. Comparative genomic and transcriptomic analyses showed that a 1.0 kb DNA fragment was deleted in shh Sm2OGD3m. Complementation assay showed that overexpression of intact Sm2OGD3 in shh hairy roots recovered furan D-ring tanshinone accumulation. Consistently, in vitro protein assay showed that Sm2OGD3 catalyzed the conversion of cyptotanshinone, 15,16-dihydrotanshinone I and 1,2,15,16-tetrahydrotanshinone I into tanshinone IIA, tanshinone I and 1,2-dihydrotanshinone I, respectively. Thus, Sm2OGD3 functions as tanshinone 15,16-dehydrogenase and is a key enzyme in tanshinone biosynthesis. The results provide novel insights into the metabolic network of medicinally important tanshinone compounds.

Phosphomevalonate kinase regulates the MVA/MEP pathway in mango during ripening

Mango is a popular tropical fruit with a great diversity in taste and aroma, contributed primarily by terpenoids. Phosphomevalonate kinase (PMK) is a key enzyme for isoprenoid biosynthesis in the mevalonic acid (MVA) pathway responsible for terpenoids. In this study, two cultivars of mango, "Dashehari" and "Banganpalli", showing opposite spatio-temporal patterns of ripening polarity, were investigated for studying the role of MiPMK in aroma production. MiPMK transcription and enzyme activity increased during ripening in both varieties. Expression in the early-ripening inner zones preceded that in the later-ripening outer zones of "Dashehari" while it was higher in the early ripening outer zones in "Banganpalli". Polypeptide sequences of the two enzymes showed differences in a few amino acids that were also reflected in kinetic properties such as specific activity and pH optima. Silencing of MiPMK in "Dashehari" fruit by VIGS suppressed the kinase activity and led to changes in relative contributions of the mevalonic acid (MVA) and methylerythritol 4-phosphate (MEP) pathways. This also altered the fruit metabolite profile with a reduction/disappearance of sesquiterpenes such as geranyl geraniol, trans-farnesol, β-caryophyllene, β-pinene, bisabolene and guaiane but the appearance of menthol and d-limonene in silenced fruit. The study shows that MiPMK levels may control downstream metabolite flux of the MVA pathway in mango.

Genome-Wide Identification and Abiotic Stress Response Analysis of PP2C Gene Family in Woodland and Pineapple Strawberries

Protein phosphatase 2C (PP2C) is a negative regulator of serine/threonine residue protein phosphatase and plays an important role in abscisic acid (ABA) and abiotic-stress-mediated signaling pathways in plants. The genome complexity of woodland strawberry and pineapple strawberry is different due to the difference in chromosome ploidy. This study conducted a genome-wide investigation of the FvPP2C (Fragaria vesca) and FaPP2C (Fragaria ananassa) gene family. Fifty-six FvPP2C genes and 228 FaPP2C genes were identified from the woodland strawberry and pineapple strawberry genomes, respectively. FvPP2Cs were distributed on seven chromosomes, and FaPP2Cs were distributed on 28 chromosomes. The size of the FaPP2C gene family was significantly different from that of the FvPP2C gene family, but both FaPP2Cs and FvPP2Cs were localized in the nucleus, cytoplasm, and chloroplast. Phylogenetic analysis revealed that 56 FvPP2Cs and 228 FaPP2Cs could be divided into 11 subfamilies. Collinearity analysis showed that both FvPP2Cs and FaPP2Cs had fragment duplication, and the whole genome duplication was the main cause of PP2C gene abundance in pineapple strawberry. FvPP2Cs mainly underwent purification selection, and there were both purification selection and positive selection effects in the evolution of FaPP2Cs. Cis-acting element analysis found that the PP2C family genes of woodland and pineapple strawberries mainly contained light responsive elements, hormone responsive elements, defense and stress responsive elements, and growth and development-related elements. The results of quantitative real-time PCR (qRT-PCR) showed that the FvPP2C genes showed different expression patterns under ABA, salt, and drought treatment. The expression level of FvPP2C18 was upregulated after stress treatment, which may play a positive regulatory role in ABA signaling and abiotic stress response mechanisms. This study lays a foundation for further investigation on the function of the PP2C gene family.

Uncovering a miltiradiene biosynthetic gene cluster in the Lamiaceae reveals a dynamic evolutionary trajectory

The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris, Plectranthus barbatus, and Leonotis leonurus. We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites.

Integrated analysis of lncRNAs, mRNAs, and TFs to identify network modules underlying diterpenoid biosynthesis in Salvia miltiorrhiza

Long non-coding RNAs (lncRNAs) are transcripts of more than 200 nucleotides (nt) in length, with minimal or no protein-coding capacity. Increasing evidence indicates that lncRNAs play important roles in the regulation of gene expression including in the biosynthesis of secondary metabolites. Salvia miltiorrhiza Bunge is an important medicinal plant in China. Diterpenoid tanshinones are one of the main active components of S. miltiorrhiza. To better understand the role of lncRNAs in regulating diterpenoid biosynthesis in S. miltiorrhiza, we integrated analysis of lncRNAs, mRNAs, and transcription factors (TFs) to identify network modules underlying diterpenoid biosynthesis based on transcriptomic data. In transcriptomic data, we obtained 6,651 candidate lncRNAs, 46 diterpenoid biosynthetic pathway genes, and 11 TFs involved in diterpenoid biosynthesis. Combining the co-expression and genomic location analysis, we obtained 23 candidate lncRNA-mRNA/TF pairs that were both co-expressed and co-located. To further observe the expression patterns of these 23 candidate gene pairs, we analyzed the time-series expression of S. miltiorrhiza induced by methyl jasmonate (MeJA). The results showed that 19 genes were differentially expressed at least a time-point, and four lncRNAs, two mRNAs, and two TFs formed three lncRNA-mRNA and/or TF network modules. This study revealed the relationship among lncRNAs, mRNAs, and TFs and provided new insight into the regulation of the biosynthetic pathway of S. miltiorrhiza diterpenoids.

Functional Characterization of Three Diterpene Synthases Responsible for Tetracyclic Diterpene Biosynthesis in Scoparia dulcis

Scoparia dulcis produces unique biologically active diterpenoids such as scopadulcic acid B (SDB). They are biosynthesized from geranylgeranyl diphosphate (GGPP) via syn-copalyl diphosphate (syn-CPP) and scopadulanol as an important key intermediate. In this paper, we functionally characterized three diterpene synthases, SdCPS2, SdKSL1 and SdKSL2, from S. dulcis. The SdCPS2 catalyzed a cyclization reaction from GGPP to syn-CPP, and SdKSL1 did from syn-CPP to scopadulan-13α-ol. On the other hand, SdKSL2 was found to incorporate a non-sense mutation at 682. Therefore, we mutated the nucleotide residue from A to G in SdKSL2 to produce SdKSL2mut, and it was able to recover the catalytic function from syn-CPP to syn-aphidicol-16-ene, the precursor to scopadulin. From our results, SdCPS2 and SdKSL1 might be important key players for SDB biosynthesis in S. dulcis.

SmDXS5, acting as a molecular valve, plays a key regulatory role in the primary and secondary metabolism of tanshinones in Salvia miltiorrhiza

Red sage, the dry root and rhizome of the herbaceous plant Salvia miltiorrhiza Bunge, is widely used for treating various diseases. The low content of tanshinones (terpenoids) has always restricted development of the S. miltiorrhiza industry. Here, we found that SmDXS5, a rate-limiting enzyme-coding gene located at the intersection of primary and secondary metabolism, can effectively change the transcription level and secondary metabolome profile of hairy roots of S. miltiorrhiza, and significantly increase the content of tanshinones. Agrobacterium rhizogenes was used to infuse S. miltiorrhiza explants, and hairy roots of S. miltiorrhiza expressing the SmDXS5 gene were obtained successfully. We identified 39 differentially accumulated metabolites (DAMs) by metabolomics based on ultra-high performance liquid chromatography quadrupole exactive mass spectrometry and multivariate statistics. These DAMs might be key metabolites of SmDXS5 gene regulation. RNA sequencing was used to compare gene expression between the hairy roots of the SmDXS5 overexpressing group and the blank control (BC) group. Compared with the BC group, 18,646 differentially expressed genes were obtained: 8994 were upregulated and 9,652 downregulated. The combined transcriptome and metabolome analyses revealed that the mevalonate and methylerythritol phosphate pathways and synthase gene expression levels in the SmDXS5 overexpressing group were upregulated significantly, and the accumulation of tanshinone components was increased significantly, which promoted the process of glycolysis and promoted the transformation of carbohydrates to secondary metabolism. Moreover, the expression of SmPAL, the first rate-limiting enzyme gene of the phenylpropane pathway, decreased, reducing the accumulation of phenolic acid, another secondary metabolite. Therefore, SmDXS5 can be defined as a 'valve' gene, mainly responsible for regulating the distribution of primary and secondary metabolic flow of tanshinones in S. miltiorrhiza, and for other secondary metabolic pathways. The discovery of SmDXS5 and its molecular valve function in regulating primary and secondary metabolism will provide a basis for the industrial production of tanshinone components, and cultivation of high quality S. miltiorrhiza.

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.

Targeting autophagy in atherosclerosis: Advances and therapeutic potential of natural bioactive compounds from herbal medicines and natural products

Atherosclerosis (AS) is the most common causes of cardiovascular disease characterized by the formation of atherosclerotic plaques in the arterial wall, and it has become a dominant public health problem that seriously threaten people worldwide. Autophagy is a cellular self-catabolism process, which is critical to protect cellular homeostasis against harmful conditions. Emerging evidence suggest that dysregulated autophagy is involved in the development of AS. Therefore, pharmacological interventions have been developed to inhibit the AS via autophagy induction. Among various AS treating methods, herbal medicines and natural products have been applied as effective complementary and alternative medicines to ameliorate AS and its associated cardiovascular disease. Recently, mounting evidence revealed that natural bioactive compounds from herbs and natural products could induce autophagy to suppress the occurrence and development of AS, by promoting cholesterol efflux, reducing plaque inflammation, and inhibiting apoptosis or senescence. In the present review, we highlight recent findings regarding possible effects and molecular mechanism of natural compounds in autophagy-targeted mitigation of atherosclerosis, aiming to provide new potential therapeutic strategies for the atherosclerosis treatment preclinically and clinically.

Comparative transcriptomics of Pinus massoniana organs provides insights on terpene biosynthesis regulation

Terpenoids are the most important natural products collected from conifer species. However, the molecular mechanisms and core factors underlying terpenoid biosynthesis in Pinus massoniana remain unclear. To clarify these mechanisms, this study aimed to identify potential genes that might participate in the terpenoid biosynthesis of P. massoniana. In this study, single molecule real-time (SMRT) sequencing and expression analysis were used to confirm the expression patterns of genes involved in the cones, immature needles, mature needles, immature branches, and mature branches of P. massoniana. A total of 31,331 lncRNAs and 71,240 mRNAs were identified from these organs, and the greatest number of differentially expressed genes (DEGs) was discovered between needles and branches. Weighted gene coexpression network analysis (WGCNA) classified all expressed genes into nine typical modules with 11 kinds of transcription factors (TFs), namely, AP2-ERF, ARF, AUX-IAA, C2H2, Dof, F-box, SBP, WRKY, bHLH, bZIP, and GRAS, and seven kinds of functional genes, namely, ABC transporter, cellulose synthase (CesA), leucine-rich repeats (LRR), cytochrome P450 (CYT P450), pathogenesis-related protein (PR), terpene synthase (TPS), and chlorophyllase enzyme. A molecular network was constructed for hub genes, TFs, and functional genes in three modules. The potential function of eight candidate genes, including PmbHLH2, PmERF1, PmRGA, PmGAI, PmbZIP1, PmLOB1, PmMADS1, and PmMYB1, was validated through correlation analysis between terpenoid contents and expression levels, subcellular localization, and transcriptional activation activity, which provides us with probable regulators of terpenoid biosynthesis in conifers.

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.

Identification and Characterization of Jasmonic Acid Biosynthetic Genes in Salvia miltiorrhiza Bunge

Jasmonic acid (JA) is a vital plant hormone that performs a variety of critical functions for plants. Salvia miltiorrhiza Bunge (S. miltiorrhiza), also known as Danshen, is a renowned traditional Chinese medicinal herb. However, no thorough and systematic analysis of JA biosynthesis genes in S. miltiorrhiza exists. Through genome-wide prediction and molecular cloning, 23 candidate genes related to JA biosynthesis were identified in S. miltiorrhiza. These genes belong to four families that encode lipoxygenase (LOX), allene oxide synthase (AOS), allene oxide cyclase (AOC), and 12-OPDA reductase3 (OPR3). It was discovered that the candidate genes for JA synthesis of S. miltiorrhiza were distinct and conserved, in contrast to related genes in other plants, by evaluating their genetic structures, protein characteristics, and phylogenetic trees. These genes displayed tissue-specific expression patterns concerning to methyl jasmonate (MeJA) and wound tests. Overall, the results of this study provide valuable information for elucidating the JA biosynthesis pathway in S. miltiorrhiza by comprehensive and methodical examination.

Isolation and Comprehensive in Silico Characterisation of a New 3-Hydroxy-3-Methylglutaryl-Coenzyme A Reductase 4 (HMGR4) Gene Promoter from Salvia miltiorrhiza: Comparative Analyses of Plant HMGR Promoters

Salvia miltiorrhiza synthesises tanshinones with multidirectional therapeutic effects. These compounds have a complex biosynthetic pathway, whose first rate limiting enzyme is 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR). In the present study, a new 1646 bp fragment of the S. miltiorrhiza HMGR4 gene consisting of a promoter, 5′ untranslated region and part of a coding sequence was isolated and characterised in silico using bioinformatics tools. The results indicate the presence of a TATA box, tandem repeat and pyrimidine-rich sequence, and the absence of CpG islands. The sequence was rich in motifs recognised by specific transcription factors sensitive mainly to light, salicylic acid, bacterial infection and auxins; it also demonstrated many binding sites for microRNAs. Moreover, our results suggest that HMGR4 expression is possibly regulated during flowering, embryogenesis, organogenesis and the circadian rhythm. The obtained data were verified by comparison with microarray co-expression results obtained for Arabidopsis thaliana. Alignment of the isolated HMGR4 sequence with other plant HMGRs indicated the presence of many common binding sites for transcription factors, including conserved ones. Our findings provide valuable information for understanding the mechanisms that direct transcription of the S. miltiorrhiza HMGR4 gene.

The Expression Profiles of the Salvia miltiorrhiza 3-Hydroxy-3-methylglutaryl-coenzyme A Reductase 4 Gene and Its Influence on the Biosynthesis of Tanshinones

Salvia miltiorrhiza is a medicinal plant that synthesises biologically-active tanshinones with numerous therapeutic properties. An important rate-limiting enzyme in the biosynthesis of their precursors is 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR). This study presents the organ-specific expression profile of the S. miltiorrhiza HMGR4 gene and its sensitivity to potential regulators, viz. gibberellic acid (GA₃), indole-3-acetic acid (IAA) and salicylic acid (SA). In addition, it demonstrates the importance of the HMGR4 gene, the hormone used, the plant organ, and the culture environment for the biosynthesis of tanshinones. HMGR4 overexpression was found to significantly boost the accumulation of dihydrotanshinone I (DHTI), cryptotanshinone (CT), tanshinone I (TI) and tanshinone IIA (TIIA) in roots by 0.44 to 5.39 mg/g dry weight (DW), as well as TIIA in stems and leaves. S. miltiorrhiza roots cultivated in soil demonstrated higher concentrations of the examined metabolites than those grown in vitro. GA₃ caused a considerable increase in the quantity of CT (by 794.2 µg/g DW) and TIIA (by 88.1 µg/g DW) in roots. In turn, IAA significantly inhibited the biosynthesis of the studied tanshinones in root material.

Carotenoid Biosynthesis: Genome-Wide Profiling, Pathway Identification in Rhodotorula glutinis X-20, and High-Level Production

Rhodotorula glutinis, as a member of the family Sporidiobolaceae, is of great value in the field of biotechnology. However, the evolutionary relationship of R. glutinis X-20 with Rhodosporidiobolus, Sporobolomyces, and Rhodotorula are not well understood, and its metabolic pathways such as carotenoid biosynthesis are not well resolved. Here, genome sequencing and comparative genome techniques were employed to improve the understanding of R. glutinis X-20. Phytoene desaturase (crtI) and 15-cis-phytoene synthase/lycopene beta-cyclase (crtYB), key enzymes in carotenoid pathway from R. glutinis X-20 were more efficiently expressed in S. cerevisiae INVSc1 than in S. cerevisiae CEN.PK2-1C. High yielding engineered strains were obtained by using synthetic biology technology constructing carotenoid pathway in S. cerevisiae and optimizing the precursor supply after fed-batch fermentation with palmitic acid supplementation. Genome sequencing analysis and metabolite identification has enhanced the understanding of evolutionary relationships and metabolic pathways in R. glutinis X-20, while heterologous construction of carotenoid pathway has facilitated its industrial application.

Transcriptome and metabolome reveal the accumulation of secondary metabolites in different varieties of Cinnamomum longepaniculatum

BACKGROUND

Cinnamomum longepaniculatum (Gamble) N. Chao ex H. W. Li, whose leaves produce essential oils, is a traditional Chinese medicine and economically important tree species. In our study, two C. longepaniculatum varieties that have significantly different essential oil contents and leaf phenotypes were selected as the materials to investigate secondary metabolism.

RESULT

The essential oil content and leaf phenotypes were different between the two varieties. When the results of both transcriptome and metabolomic analyses were combined, it was found that the differences were related to phenylalanine metabolic pathways, particularly the metabolism of flavonoids and terpenoids. The transcriptome results based on KEGG pathway enrichment analysis showed that pathways involving phenylpropanoids, tryptophan biosynthesis and terpenoids significantly differed between the two varieties; 11 DEGs (2 upregulated and 9 downregulated) were associated with the biosynthesis of other secondary metabolites, and 12 DEGs (2 upregulated and 10 downregulated) were related to the metabolism of terpenoids and polyketides. Through further analysis of the leaves, we detected 196 metabolites in C. longepaniculatum. The abundance of 49 (26 downregulated and 23 upregulated) metabolites differed between the two varieties, which is likely related to the differences in the accumulation of these metabolites. We identified 12 flavonoids, 8 terpenoids and 8 alkaloids and identified 4 kinds of PMFs from the leaves of C. longepaniculatum.

CONCLUSIONS

The combined results of transcriptome and metabolomic analyses revealed a strong correlation between metabolite contents and gene expression. We speculate that light leads to differences in the secondary metabolism and phenotypes of leaves of different varieties of C. longepaniculatum. This research provides data for secondary metabolite studies and lays a solid foundation for breeding ideal C. longepaniculatum plants.

Integrated metabolite profiling and transcriptome analysis reveals tissue-specific regulation of terpenoid biosynthesis in Artemisia argyi

Artemisia argyi L. is a widely distributed medicinal plant in China. The major bioactive substances of essential oils extracted from leaves are terpenoids. Although many researches have studied the pharmacological effects of the essential oils, the tissue-specific accumulation of terpenoid biosynthesis and the regulatory networks in A. argyi are poorly understood. This study conducted an integrated metabolomic and transcriptomic analysis of roots, stems, and leaves to investigate the tissue-specific regulatory network of terpenoid biosynthesis in A. argyi. We identified 77 unigenes putatively involved in terpenoid backbone biosynthesis. Three rate-determining enzyme genes (DXS, DXR, and HDR) of the methylerythritol phosphate pathway were predominantly expressed in leaves, and strongly co-expressed with eight transcription factors (2 MYBs, 4 WRKYs, and 2 AP2s). An metabolite-transcript correlation analysis revealed 26 putative cytochrome P450s related to terpenoid metabolism in leaves. These results provide a foundation for the future metabolic engineering of useful terpenoids in A. argyi.

Tanshinones: Leading the way into Lamiaceae labdane-related diterpenoid biosynthesis

Tanshinones are the bioactive diterpenoid constituents of the traditional Chinese medicinal herb Danshen (Salvia miltiorrhiza), and are examples of the phenolic abietanes widely found within the Lamiaceae plant family. Due to the significant interest in these labdane-related diterpenoid natural products, their biosynthesis has been intensively investigated. In addition to providing the basis for metabolic engineering efforts, this work further yielded pioneering insights into labdane-related diterpenoid biosynthesis in the Lamiaceae more broadly. This includes stereochemical foreshadowing of aromatization, with novel protein domain loss in the relevant diterpene synthase, as well as broader phylogenetic conservation of the relevant enzymes. Beyond such summary of more widespread metabolism, formation of the furan ring that characterizes the tanshinones also has been recently elucidated. Nevertheless, the biocatalysts for the pair of demethylations remain unknown, and the intriguing potential connection of these reactions to the further aromatization observed in the tanshinones are speculated upon here.

Methyl Jasmonate Activates the 2C Methyl-D-erithrytol 2,4-cyclodiphosphate Synthase Gene and Stimulates Tanshinone Accumulation in Salvia miltiorrhiza Solid Callus Cultures

The present study characterizes the 5′ regulatory region of the SmMEC gene. The isolated fragment is 1559 bp long and consists of a promoter, 5′UTR and 31 nucleotide 5′ fragments of the CDS region. In silico bioinformatic analysis found that the promoter region contains repetitions of many potential cis-active elements. Cis-active elements associated with the response to methyl jasmonate (MeJa) were identified in the SmMEC gene promoter. Co-expression studies combined with earlier transcriptomic research suggest the significant role of MeJa in SmMEC gene regulation. These findings were in line with the results of the RT-PCR test showing SmMEC gene expression induction after 72 h of MeJa treatment. Biphasic total tanshinone accumulation was observed following treatment of S. miltiorrhiza solid callus cultures with 50–500 μM methyl jasmonate, with peaks observed after 10–20 and 50–60 days. An early peak of total tanshinone concentration (0.08%) occurred after 20 days of 100 μM MeJa induction, and a second, much lower one, was observed after 50 days of 50 μM MeJa stimulation (0.04%). The dominant tanshinones were cryptotanshinone (CT) and dihydrotanshinone (DHT). To better understand the inducing effect of MeJa treatment on tanshinone biosynthesis, a search was performed for methyl jasmonate-responsive cis-active motifs in the available sequences of gene proximal promoters associated with terpenoid precursor biosynthesis. The results indicate that MeJa has the potential to induce a significant proportion of the presented genes, which is in line with available transcriptomic and RT-PCR data.

Integrated Transcriptomics and Proteomics to Reveal Regulation Mechanism and Evolution of SmWRKY61 on Tanshinone Biosynthesis in Salvia miltiorrhiza and Salvia castanea

Tanshinones found in Salvia species are the main active compounds for the treatment of cardiovascular and cerebrovascular diseases, but their contents are hugely different in different species. For example, tanshinone IIA content in Salvia castanea Diels f. tomentosa Stib. is about 49 times higher than that in Salvia miltiorrhiza Bunge. The molecular mechanism responsible for this phenomenon remains largely unknown. To address this, we performed comparative transcriptomic and proteomic analyses of S. miltiorrhiza and S. castanea. A total of 296 genes in S. castanea and 125 genes in S. miltiorrhiza were highly expressed at both the transcriptional and proteome levels, including hormone signal regulation, fungus response genes, transcription factors, and CYP450. Among these differentially expressed genes, the expression of SmWRKY61 was particularly high in S. castanea. Overexpression of SmWRKY61 in S. miltiorrhiza could significantly increase the content of tanshinone I and tanshinone IIA, which were 11.09 and 33.37 times of the control, respectively. Moreover, SmWRKY61 had a strong regulatory effect, elevating the expression levels of tanshinone pathway genes such as DXS2, CMK, HMGS2, 1, KSL1, KSL2, CYP76AH1, and CYP76AK3. For the WRKY family, 79 SmWRKYs were originally obtained and classified into three main groups. Collinearity analysis indicated a more specific extension of WRKY gene family in Salvia genus. In 55 Salvia species, only 37 species contained the WRKY61 sequence, and high SmWRKY61 expression in some Salvia L. species was often accompanied by high tanshinone accumulation. The above results suggest that SmWRKY61 is a highly effective regulator of tanshinone accumulation and may be a key factor resulting in high tanshinone accumulation in S. castanea.

Endophytic fungus Cladosporium tenuissimum DF11, an efficient inducer of tanshinone biosynthesis in Salvia miltiorrhiza roots

Salvia miltiorrhiza is a traditional medicinal plant mainly used for cardiovascular and cerebrovascular disease treatment. Tanshinones are the main bioactive constituents of S. miltiorrhiza, which mainly accumulate around its root periderm tissue. Endophytic fungi are important bioelicitors or probiotics that can promote the accumulation of secondary metabolites and sustainable cultivation of medicinal plants. Among them, endophytic Cladosporium spp., possessing a variety of biotransformation and metabolic abilities, is an ideal elicitor source. Here, we used a gnotobiotic system to investigate the effects of the endophytic fungus Cladosporium tenuissimum DF11 on tanshinone biosynthesis in S. miltiorrhiza roots. The results showed that C. tenuissimum DF11 mainly colonizes the intercellular space of the root tissues and promotes tanshinone biosynthesis and accumulation in S. miltiorrhiza roots by upregulating the expression of the genes encoding for key enzymes HMGR, DXS, DXR, GGPPS, CPS, KSL and CYP76AH1 of the tanshinone biosynthesis pathway. The expression levels of almost all genes encoding for key enzymes reached the response peak in the first or second week after DF11 colonization. Taken together, the endophytic fungus C. tenuissimum DF11 could promote secondary metabolite accumulation in S. miltiorrhiza roots. These results indicate that DF11 will be a potential biofertilizer fungus to regulate and stabilize the quality of cultivated S. miltiorrhiza medicinal materials.

Transcriptome and HPLC Analysis Reveal the Regulatory Mechanisms of Aurantio-Obtusin in Space Environment-Induced Senna obtusifolia Lines

Senna obtusifolia is a famous medicinal plant that is widely used in Asian countries. Its seed plays an important role in the treatment of many diseases because it contains various anthraquinones and flavonoids. Our previous studies have indicated that three space environment-induced S. obtusifolia lines (SP-lines) i.e., QC10, QC29, and QC46, have higher seed yield and aurantio-obtusin (AO) content. However, the underlying mechanism of higher AO content in SP-lines is still unknown. Herein, transcriptome sequencing and HPLC were employed to analyze the differences between SP-lines and ground control (GC3) and elucidate the regulatory mechanisms of AO accumulation in SP-lines. The results show that 4002 differentially expressed genes (DEGs) were identified in SP-lines versus (vs.) GC3. DEGs in the QC10 vs. GC3, QC29 vs. GC3, and QC46 vs. GC3 comparisons were classified into 28, 36, and 81 GO terms and involved in 63, 74, and 107 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. KEGG pathway and gene expression analysis revealed that DEGs involved in anthraquinone pathways were significantly elevated in QC10 and QC46. Integrating the results of GO annotation, KEGG enrichment, and gene expression analysis, we propose that the elevated genes such as DAHPS, DHQS, and MenB enhance the metabolic flux in the anthraquinone pathway and promote AO content in QC10 and QC46. Taken together, this study elucidated the mechanism of AO content in SP-lines and provides valuable genetic information for S. obtusifolia. In addition, to the best of our knowledge, this study presents the first transcriptome analysis of environment-induced medicinal plants and paves the way to select elite S. obtusifolia varieties in the future.

Isoprenoid biosynthesis regulation in poplars by methylerythritol phosphate and mevalonic acid pathways

It is critical to develop plant isoprenoid production when dealing with human-demanded industries such as flavoring, aroma, pigment, pharmaceuticals, and biomass used for biofuels. The methylerythritol phosphate (MEP) and mevalonic acid (MVA) plant pathways contribute to the dynamic production of isoprenoid compounds. Still, the cross-talk between MVA and MEP in isoprenoid biosynthesis is not quite recognized. Regarding the rate-limiting steps in the MEP pathway through catalyzing 1-deoxy-D-xylulose5-phosphate synthase and 1-deoxy-D-xylulose5-phosphate reductoisomerase (DXR) and also the rate-limiting step in the MVA pathway through catalyzing 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), the characterization and function of HMGR from Populus trichocarpa (PtHMGR) were analyzed. The results indicated that PtHMGR overexpressors (OEs) displayed various MEP and MVA-related gene expressions compared to NT poplars. The overexpression of PtDXR upregulated MEP-related genes and downregulated MVA-related genes. The overexpression of PtDXR and PtHMGR affected the isoprenoid production involved in both MVA and MEP pathways. Here, results illustrated that the PtHMGR and PtDXR play significant roles in regulating MEP and MVA-related genes and derived isoprenoids. This study clarifies cross-talk between MVA and MEP pathways. It demonstrates the key functions of HMGR and DXR in this cross-talk, which significantly contribute to regulate isoprenoid biosynthesis in poplars.

Overexpression of SmSCR1 Promotes Tanshinone Accumulation and Hairy Root Growth in Salvia miltiorrhiza

Lipid-soluble tanshinone is one of the main bioactive substances in the medicinal plant Salvia miltiorrhiza, and its medicinal demand is growing rapidly. Yeast extract (YE) modulates the tanshinone biosynthesis, but the underlying regulatory network remains obscure. In this study, a YE-responsive transcriptional factor Scarecrow1 (SCR1) was identified in S. miltiorrhiza from the YE-induced transcriptome dataset. SmSCR1 is located in the nucleus. Overexpression of SmSCR1 in S. miltiorrhiza roots resulted in a significantly higher accumulation of tanshinone than the control, with the highest 1.49-fold increase. We also detected upregulation of tanshinone biosynthetic genes, SmSCR1 and SmHMGR1, and distinct alteration of growth and development of the hairy roots in the overexpression lines compared to the control. An inverse phenotype was observed in SmSCR1-SRDX suppression expression lines. We found that SmSCR1 can bind to the promoter of SmCPS1 to induce its expression. This study provides new insight into the regulatory mechanism on the growth and development of hairy roots, tanshinone accumulation, and the metabolic engineering of bioactive compounds in S. miltiorrhiza.

Chitosan elicitation of saponin accumulation in Psammosilene tunicoides hairy roots by modulating antioxidant activity, nitric oxide production and differential gene expression

Elicitation is one of the most effective strategies for enhancing plant bioactive compounds, such as triterpenoid saponins. Chitosan gained worldwide attention for biostimulant activity, but little is known about its roles in the elicitation of triterpenoid saponin in medicinal plants. Here, we explored the regulatory network of chitosan on saponin accumulation in hairy root cultures of Psammosilene tunicoides, a valuable medicinal herb known for its pain-relieving properties endemic to China. Compared with control, the highest total saponin accumulation exhibited a 4.55-fold enhancement in hairy roots elicited by 200 mg L^-1 chitosan for nine days. High-performance liquid chromatography (HPLC) revealed the yields of quillaic acid, gypsogenin and gypsogenin-3-O-β-D-glucuronopyranoside were significantly increased after chitosan treatments. Moreover, exogenous chitosan application dramatically triggered the reactive oxygen species (ROS) scavenging enzyme activities and nitric oxide (NO) content in hairy roots. Comparative transcriptome analysis from chitosan-treated (1 and 9 d) or control groups revealed that differentially expressed genes (DEGs) were greatly enriched in plant-pathogen interaction and metabolic processes. The transcriptions of candidate DEGs involved in chitosan-elicited saponin metabolism were increased, especially genes encoding antioxidant enzymes (SOD, POD and GR), stress-responsive transcription factors (WRKYs and NACs) and terpenoid biosynthetic enzymes (DXS, GPPS and SE). Taken together, these results indicate that chitosan elicitor promotes triterpenoid saponin biosynthesis by enhancing antioxidant activities, NO production and differential gene expression in P. tunicoides hairy roots.

Cloning and expression analysis of mevalonate kinase and phosphomevalonate kinase genes associated with the MVA pathway in Santalum album

Sandalwood (Santalum album L.) is highly valued for its fragrant heartwood and extracted oil. Santalols, which are the main components of that oil, are terpenoids, and these are biosynthesized via the mevalonic acid (MVA) pathway. Mevalonate kinase (MK) and phosphomevalonate kinase (PMK) are key enzymes in the MVA pathway. Little is known about the genes that encode MK and PMK in S. album or the mechanism that regulates their expression. To isolate and identify the functional genes involved in santalol biosynthesis in S. album, an MK gene designated as SaMK, and a PMK gene designated as SaPMK, were cloned from S. album. The sequences of these genes were analyzed. A bioinformatics analysis was conducted to assess the homology of SaMK and SaPMK with MK and PMK genes from other plants. The subcellular localization of SaMK and SaPMK proteins was also investigated, as was the functional complementation of SaMK and SaPMK in yeast. Our results show that the full-length cDNA sequences of SaMK and SaPMK were 1409 bp and 1679 bp long, respectively. SaMK contained a 1381 bp open reading frame (ORF) encoding a polypeptide of 460 amino acids and SaPMK contained a 1527 bp ORF encoding a polypeptide of 508 amino acids. SaMK and SaPMK showed high homology with MK and PMK genes of other plant species. Functional complementation of SaMK in a MK-deficient mutant yeast strain YMR208W and SaPMK in a PMK-deficient mutant yeast strain YMR220W confirmed that cloned SaMK and SaPMK cDNA encode a functional MK and PMK, respectively, mediating MVA biosynthesis in yeast. An analysis of tissue expression patterns revealed that SaMK and SaPMK were constitutively expressed in all the tested tissues. SaMK was highly expressed in young leaves but weakly expressed in sapwood. SaPMK was highly expressed in roots and mature leaves, but weakly expressed in young leaves. Induction experiments with several elicitors showed that SaMK and SaPMK expression was upregulated by methyl jasmonate. These results will help to further study the role of MK and PMK genes during santalol biosynthesis in S. album.

Chia (Salvia hispanica) Gene Expression Atlas Elucidates Dynamic Spatio-Temporal Changes Associated With Plant Growth and Development

Chia (Salvia hispanica L.), now a popular superfood and a pseudocereal, is one of the richest sources of dietary nutrients such as protein, fiber, and polyunsaturated fatty acids (PUFAs). At present, the genomic and genetic information available in the public domain for this crop are scanty, which hinders an understanding of its growth and development and genetic improvement. We report an RNA-sequencing (RNA-Seq)-based comprehensive transcriptome atlas of Chia sampled from 13 tissue types covering vegetative and reproductive growth stages. We used ~355 million high-quality reads of total ~394 million raw reads from transcriptome sequencing to generate de novo reference transcriptome assembly and the tissue-specific transcript assemblies. After the quality assessment of the merged assemblies and implementing redundancy reduction methods, 82,663 reference transcripts were identified. About 65,587 of 82,663 transcripts were translated into 99,307 peptides, and we were successful in assigning InterPro annotations to 45,209 peptides and gene ontology (GO) terms to 32,638 peptides. The assembled transcriptome is estimated to have the complete sequence information for ~86% of the genes found in the Chia genome. Furthermore, the analysis of 53,200 differentially expressed transcripts (DETs) revealed their distinct expression patterns in Chia's vegetative and reproductive tissues; tissue-specific networks and developmental stage-specific networks of transcription factors (TFs); and the regulation of the expression of enzyme-coding genes associated with important metabolic pathways. In addition, we identified 2,411 simple sequence repeats (SSRs) as potential genetic markers from the transcripts. Overall, this study provides a comprehensive transcriptome atlas, and SSRs, contributing to building essential genomic resources to support basic research, genome annotation, functional genomics, and molecular breeding of Chia.

Consequences of LED Lights on Root Morphological Traits and Compounds Accumulation in Sarcandra glabra Seedlings

This study evaluated the effects of different light spectra (white light; WL, blue light; BL and red light; RL) on the root morphological traits and metabolites accumulation and biosynthesis in Sarcandra glabra. We performed transcriptomic and metabolomic profiling by RNA-seq and ultra-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UPLC-ESI-MS/MS), respectively. When morphological features were compared to WL, BL substantially increased under-ground fresh weight, root length, root surface area, and root volume, while RL inhibited these indices. A total of 433 metabolites were identified, of which 40, 18, and 68 compounds differentially accumulated in roots under WL (WG) vs. roots under BL (BG), WG vs. roots under RL (RG), and RG vs. BG, respectively. In addition, the contents of sinapyl alcohol, sinapic acid, fraxetin, and 6-methylcoumarin decreased significantly in BG and RG. In contrast, chlorogenic acid, rosmarinyl glucoside, quercitrin and quercetin were increased considerably in BG. Furthermore, the contents of eight terpenoids compounds significantly reduced in BG. Following transcriptomic profiling, several key genes related to biosynthesis of phenylpropanoid-derived and terpenoids metabolites were differentially expressed, such as caffeic acid 3-O-methyltransferase) (COMT), hydroxycinnamoyl-CoA shikimate hydroxycinnamoyl transferase (HCT), O-methyltransferase (OMT), and 1-deoxy-D-xylulose-5-phosphate synthetase (DXS). In summary, our findings showed that BL was suitable for growth and accumulation of bioactive metabolites in root tissue of S. glabra. Exposure to a higher ratio of BL might have the potential to improve the production and quality of S. glabra seedlings, but this needs to be confirmed further.

SmKFB5 protein regulates phenolic acid biosynthesis by controlling the degradation of phenylalanine ammonia-lyase in Salvia miltiorrhiza

Phenolic acids are the major secondary metabolites and significant bioactive constituents of the medicinal plant Salvia miltiorrhiza. Many enzyme-encoding genes and transcription factors involved in the biosynthesis of phenolic acids have been identified, but the underlying post-translational regulatory mechanisms are poorly understood. Here, we demonstrate that the S. miltiorrhiza Kelch repeat F-box protein SmKFB5 physically interacts with three phenylalanine ammonia-lyase (PAL) isozymes and mediates their proteolytic turnover via the ubiquitin-26S proteasome pathway. Disturbing the expression of SmKFB5 reciprocally affected the abundance of SmPAL protein and the accumulation of phenolic acids, suggesting that SmKFB5 is a post-translational regulator responsible for the turnover of PAL and negatively controlling phenolic acids. Furthermore, we discovered that treatment of the hairy root of S. miltiorrhiza with methyl jasmonate suppressed the expression of SmKFB5 while inducing the transcription of SmPAL1 and SmPAL3. These data suggested that methyl jasmonate consolidated both transcriptional and post-translational regulation mechanisms to enhance phenolic acid biosynthesis. Taken together, our results provide insights into the molecular mechanisms by which SmKFB5 mediates the regulation of phenolic acid biosynthesis by jasmonic acid, and suggest valuable targets for plant breeders in tailoring new cultivars.

Expressional diversity of grapevine 3-Hydroxy-3-methylglutaryl-CoA reductase (VvHMGR) in different grapes genotypes

BACKGROUND

3-Hydroxy-3-methylglutaryl-CoA reductase (HMGR) is a key enzyme in the mevalonate (MVA) pathway, which regulates the metabolism of terpenoids in the cytoplasm and determines the type and content of downstream terpenoid metabolites.

RESULTS

Results showed that grapevine HMGR family has three members, such as VvHMGR1, VvHMGR2, and VvHMGR3. The expression of VvHMGRs in 'Kyoho' has tissue specificity, for example, VvHMGR1 keeps a higher expression, VvHMGR2 is the lowest, and VvHMGR3 gradually decreases as the fruit development. VvHMGR3 is closely related to CsHMGR1 and GmHMGR9 and has collinearity with CsHMGR2 and GmHMGR4. By the prediction of interaction protein, it can interact with HMG-CoA synthase, MVA kinase, FPP/GGPP synthase, diphosphate mevalonate decarboxylase, and participates in the synthesis and metabolism of terpenoids. VvHMGR3 have similar trends in expression with some of the genes of carotenoid biosynthesis and MEP pathways. VvHMGR3 responds to various environmental and phytohormone stimuli, especially salt stress and ultraviolet (UV) treatment. The expression level of VvHMGRs is diverse in grapes of different colors and aroma. VvHMGRs are significantly higher in yellow varieties than that in red varieties, whereas rose-scented varieties showed significantly higher expression than that of strawberry aroma. The expression level is highest in yellow rose-scented varieties, and the lowest in red strawberry scent varieties, especially 'Summer Black' and 'Fujiminori'.

CONCLUSION

This study confirms the important role of VvHMGR3 in the process of grape fruit coloring and aroma formation, and provided a new idea to explain the loss of grape aroma and poor coloring during production. There may be an additive effect between color and aroma in the HMGR expression aspect.

Biochemistry of Terpenes and Recent Advances in Plant Protection

Biodiversity is adversely affected by the growing levels of synthetic chemicals released into the environment due to agricultural activities. This has been the driving force for embracing sustainable agriculture. Plant secondary metabolites offer promising alternatives for protecting plants against microbes, feeding herbivores, and weeds. Terpenes are the largest among PSMs and have been extensively studied for their potential as antimicrobial, insecticidal, and weed control agents. They also attract natural enemies of pests and beneficial insects, such as pollinators and dispersers. However, most of these research findings are shelved and fail to pass beyond the laboratory and greenhouse stages. This review provides an overview of terpenes, types, biosynthesis, and their roles in protecting plants against microbial pathogens, insect pests, and weeds to rekindle the debate on using terpenes for the development of environmentally friendly biopesticides and herbicides.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Crosstalk of Multi-Omics Platforms with Plants of Therapeutic Importance

From time immemorial, humans have exploited plants as a source of food and medicines. The World Health Organization (WHO) has recorded 21,000 plants with medicinal value out of 300,000 species available worldwide. The promising modern "multi-omics" platforms and tools have been proven as functional platforms able to endow us with comprehensive knowledge of the proteome, genome, transcriptome, and metabolome of medicinal plant systems so as to reveal the novel connected genetic (gene) pathways, proteins, regulator sequences and secondary metabolite (molecule) biosynthetic pathways of various drug and protein molecules from a variety of plants with therapeutic significance. This review paper endeavors to abridge the contemporary advancements in research areas of multi-omics and the information involved in decoding its prospective relevance to the utilization of plants with medicinal value in the present global scenario. The crosstalk of medicinal plants with genomics, transcriptomics, proteomics, and metabolomics approaches will be discussed.

Transcriptomic Analysis, Cloning, Characterization, and Expression Analysis of Triterpene Biosynthetic Genes and Triterpene Accumulation in the Hairy Roots of Platycodon grandiflorum Exposed to Methyl Jasmonate

Platycodon grandiflorum is a perennial plant that has been used for medicinal purposes. Specifically, it is widely used in Northern China and Korea for the treatment of various diseases. Terpenoids belong to a group called secondary metabolites and have attracted a wide range of interest. Here, we determined the expressed sequence tag (EST) library of the methyl jasmonate (MeJA)-treated hairy root of P. grandiflorum. In total, 5760 ESTs were obtained, but after deleting the vector sequences and removing poor-quality sequences, a total of 2536 ESTs were attained. Of these, 811 contigs and 1725 singletons were annotated. The data were further analyzed with a focus on the gene families of the terpenoid biosynthetic pathway (TBP). We identified and characterized four TBP genes; among these were three full-length cDNAs encoding PgHMGS, PgMK, and PgMVD, whereas PgHMGR had a partial sequence based on the EST library database. Phylogenetic analysis and a pairwise identity matrix showed that these identified genes were closely related to those of other higher plants. Moreover, the tertiary structure and multiple alignment analysis showed that several distinct conserved motifs were present. Quantitative reverse transcription-polymerase chain reaction results revealed that TBP genes were constitutively expressed in all organs of P. grandiflorum, while the expression of transcript levels of these genes varied distinctly among the organs. Additionally, the total amount of platycosides was highly detected in the root, accumulating in concentrations 7.8 and 2.6 times higher than in the hairy root and stem, respectively, and 1.4 times higher than in the leaf and flower. The highest amount of total phytosterols was found to accumulate in the leaves at 9.3, 9.1, 1.8, and 1.6 times higher than that of the stem, root, hairy root, and flower, respectively. After the hairy root was exposed to the MeJA treatment, the transcript levels of PgHMGS, PgHMGR, PgMK, and PgMVD had significantly increased. The highest level of transcript accumulation occurred at 3 h after initial exposure for most of the genes. Meanwhile, triterpene saponin accumulation increased with the increase in the time of exposure, and at 48 h after initial exposure, the total saponin content was the highest recorded.

Comparative transcriptome analysis reveals key insights into male sterility in Salvia miltiorrhiza Bunge

Background

Large-scale heterosis breeding depends upon stable, inherited male sterility lines. We accidentally discovered a male sterility line (SW-S) in the F₁progeny of a Salvia miltiorrhiza Bunge from Shandong, China (purple flowers) crossed with a S. miltiorrhiza f. alba from Sichuan, China (white flowers). We sought to provide insights into the pollen development for male sterility in S. miltiorrhiza.

Methods

The phenotypic and cytological features of the SW-S and fertile control SW-F were observed using scanning electron microscopy and paraffin sections to identify the key stage of male sterility. Transcriptome profiles were recorded for anthers at the tetrad stage of SW-S and SW-F using Illumina RNA-Seq.

Results

The paraffin sections showed that sterility mainly occurred at the tetrad stage of microspore development, during which the tapetum cells in the anther compartment completely fell off and gradually degraded in the sterile line. There was little-to-no callose deposited around the microspore cells. The tetrad microspore was shriveled and had abnormal morphology. Therefore, anthers at the tetrad stage of SW-S and fertile control SW-F were selected for comparative transcriptome analysis. In total, 266,722,270 clean reads were obtained from SW-S and SW-F, which contained 36,534 genes. There were 2,571 differentially expressed genes (DEGs) in SW-S and SW-F, of which 63.5% were downregulated. Gene Ontology (GO) enrichment analysis indicated that the differentially expressed genes were enriched in 56 functional groups (GO terms); of these, all DEGs involved in microgametogenesis and developmental maturation were downregulated in SW-S. These results were confirmed by quantitative RT-PCR. The two GO terms contained 18 DEGs, among which eight DEGs (namely: GPAT3, RHF1A, phosphatidylinositol, PFAS, MYB96, MYB78, Cals5, and LAT52) were related to gamete development. There were 10 DEGs related to development and maturation, among which three genes were directly related to pollen development (namely: ACT3, RPK2, and DRP1C). Therefore, we believe that these genes are directly or indirectly involved in the pollen abortion of SW-S. Our study provides insight into key genes related to sterility traits in S. miltiorrhiza, and the results can be further exploited in functional and mechanism studies.

Proteomic reveals the influences of smoke-water and karrikinolide on the biosynthesis of salvianolic acids and lignins in Salvia miltiorrhiza hairy roots

The proteins related to the biosynthesis of salvianolic acids and lignins were regulated by smoke-water and karrikinolide in Salvia miltiorrhiza hairy roots. The effects of smoke-water (SW) and karrikinolide (KAR₁) on the biosynthesis of salvianolic acids and lignins in Salvia miltiorrhiza hairy roots have been studied using proteomic technology. The results showed that a total of 1290 and 1678 differentially expressed proteins were respectively obtained in SW and KAR₁ comparing to the control. Bioinformatics analysis indicated the differentially expressed proteins responding to SW and KAR₁ treatments mainly involved in macromolecule metabolic process, cell part, binding, etc., and most of the proteins were located at the cytoplasm and cell membrane, followed by nuclear. In addition, the proteins involved in salvianolic acids biosynthesis were up-regulated, including 4-coumarate-CoA ligase (EC 6.2.1.12) and shikimate O-hydroxycinnamoyl-transferase (EC 2.3.1.133). Enzymes involved in lignins biosynthesis were also identified, e.g. cinnamyl-alcohol dehydrogenase (EC 1.1.1.195) and peroxidase (EC 1.11.1.7). The results indicated that proteins related to the biosynthesis of salvianolic acids and lignins were regulated by SW and KAR₁ in S. miltiorrhiza hairy roots. This study will enhance our understanding of the mechanism by which SW and KAR₁ on the biosynthesis of salvianolic acids and lignins in S. miltiorrhiza hairy roots.