SmMYB36, a Novel R2R3-MYB Transcription Factor, Enhances Tanshinone Accumulation and Decreases Phenolic Acid Content in Salvia miltiorrhiza Hairy Roots

The MpNAC72/MpERF105-MpMYB10b module regulates anthocyanin biosynthesis in Malus 'Profusion' leaves infected with Gymnosporangium yamadae

Apple rust is a serious fungal disease affecting Malus plants worldwide. Infection with the rust pathogen Gymnosporangium yamadae induces the accumulation of anthocyanins in Malus to resist rust disease. However, the mechanism of anthocyanin biosynthesis regulation in Malus against apple rust is still unclear. Here, we show that MpERF105 and MpNAC72 are key regulators of anthocyanin biosynthesis via the ethylene-dependent pathway in M. 'Profusion' leaves under rust disease stress. Exogenous ethephon treatment promoted high expression of MpERF105 and MpNAC72 and anthocyanin accumulation in G. yamadae-infected M. 'Profusion' leaves. Overexpression of MpERF105 increased the total anthocyanin content of Malus plant material and acted by positively regulating its target gene, MpMYB10b. MpNAC72 physically interacted with MpERF105 in vitro and in planta, and the two form a protein complex. Coexpression of the two leads to higher transcript levels of MpMYB10b and higher anthocyanin accumulation. In addition, overexpression of MpERF105 or MpNAC72 enhanced the resistance of M. 'Profusion' leaves to apple rust. In conclusion, our results elucidate the mechanism by which MpERF105 and MpNAC72 are induced by ethylene in G. yamadae-infected M. 'Profusion' leaves and promote anthocyanin accumulation by mediating the positive regulation of MpMYB10b expression.

The Smi-miR858a-SmMYB module regulates tanshinone and phenolic acid biosynthesis in Salvia miltiorrhiza

Tanshinones and phenolic acids are two major classes of bioactive compounds in Salvia miltiorrhiza. Revealing the regulatory mechanism of their biosynthesis is crucial for quality improvement of S. miltiorrhiza medicinal materials. Here we demonstrated that Smi-miR858a-Smi-miR858c, a miRNA family previously known to regulate flavonoid biosynthesis, also played critical regulatory roles in tanshinone and phenolic acid biosynthesis in S. miltiorrhiza. Overexpression of Smi-miR858a in S. miltiorrhiza plants caused significant growth retardation and tanshinone and phenolic acid reduction. Computational prediction and degradome and RNA-seq analyses revealed that Smi-miR858a could directly cleave the transcripts of SmMYB6, SmMYB97, SmMYB111, and SmMYB112. Yeast one-hybrid and transient transcriptional activity assays showed that Smi-miR858a-regulated SmMYBs, such as SmMYB6 and SmMYB112, could activate the expression of SmPAL1 and SmTAT1 involved in phenolic acid biosynthesis and SmCPS1 and SmKSL1 associated with tanshinone biosynthesis. In addition to directly activating the genes involved in bioactive compound biosynthesis pathways, SmMYB6, SmMYB97, and SmMYB112 could also activate SmAOC2, SmAOS4, and SmJMT2 involved in the biosynthesis of methyl jasmonate, a significant elicitor of plant secondary metabolism. The results suggest the existence of dual signaling pathways for the regulation of Smi-miR858a in bioactive compound biosynthesis in S. miltiorrhiza.

Transcriptome and volatile compounds analyses of floral development provide insight into floral scent formation in Paeonia lactiflora 'Wu Hua Long Yu'

Herbaceous peony (Paeonia lactiflora) is a well-known ornamental plant in China, celebrated for its beautiful flowers that can emit fragrances. However, exact molecular mechanisms governing synthesis of floral volatiles within herbaceous peony remain unclear. To address this gap in knowledge, our study focused on analyzing the transcriptome and the levels of floral volatile compounds in P. lactiflora 'Wu Hua Long Yu' at different stages of flower development. Using gas chromatography-mass spectrometry (GC-MS), we obtained eighteen major volatile compounds, with monoterpenes being the dominant components among them. Our transcriptome analysis, based on pooled sequencing data, revealed the most differentially expressed genes (DEGs) existed between stages S1 and S3 of flower development. Among these DEGs, we identified 89 functional genes associated with the synthesis of volatile monoterpenes, with 28 of these genes showing a positive correlation with the release of monoterpenes. Specifically, key regulators of monoterpene synthesis in herbaceous peony appear to be 1-deoxy-D-xylulose 5-phosphate synthase (DXS), geranyl pyrophosphate synthase (GPPS), and terpene synthase (TPS). Additionally, our study identified some transcription factors (TFs) that may be involved in the biosynthesis of monoterpenes. These discoveries offer invaluable illumination into the intricate molecular underpinnings orchestrating the generation of floral fragrances in herbaceous peonies, and they offer a foundation for further research to identify and utilize candidate gene resources for this purpose.

Transcriptomic Analyses Reveal Insights into the Shared Regulatory Network of Phenolic Compounds and Steviol Glycosides in Stevia rebaudiana

Stevia rebaudiana (Bertoni) is a highly valuable crop for the steviol glycoside content in its leaves, which are no-calorie sweeteners hundreds of times more potent than sucrose. The presence of health-promoting phenolic compounds, particularly flavonoids, in the leaf of S. rebaudiana adds further nutritional value to this crop. Although all these secondary metabolites are highly desirable in S. rebaudiana leaves, the genes regulating the biosynthesis of phenolic compounds and the shared gene network between the regulation of biosynthesis of steviol glycosides and phenolic compounds still need to be investigated in this species. To identify putative candidate genes involved in the synergistic regulation of steviol glycosides and phenolic compounds, four genotypes with different contents of these compounds were selected for a pairwise comparison RNA-seq analysis, yielding 1136 differentially expressed genes. Genes that highly correlate with both steviol glycosides and phenolic compound accumulation in the four genotypes of S. rebaudiana were identified using the weighted gene co-expression network analysis. The presence of UDP-glycosyltransferases 76G1, 76H1, 85C1, and 91A1, and several genes associated with the phenylpropanoid pathway, including peroxidase, caffeoyl-CoA O-methyltransferase, and malonyl-coenzyme A:anthocyanin 3-O-glucoside-6″-O-malonyltransferase, along with 21 transcription factors like SCL3, WRK11, and MYB111, implied an extensive and synergistic regulatory network involved in enhancing the production of such compounds in S. rebaudiana leaves. In conclusion, this work identified a variety of putative candidate genes involved in the biosynthesis and regulation of particular steviol glycosides and phenolic compounds that will be useful in gene editing strategies for increasing and steering the production of such compounds in S. rebaudiana as well as in other species.

Comparative transcriptome analysis and identification of candidate R2R3-MYB genes involved in anthraquinone biosynthesis in Rheum palmatum L

BACKGROUND

Rheum palmatum L. has important medicinal value because it contains biologically active anthraquinones. However, the key genes and TFs involved in anthraquinone biosynthesis and regulation in R. palmatum remain unclear.

METHODS

Based on full length transcriptome data, in this study, we screened the differentially expressed genes in the anthraquinone biosynthesis pathway. The R2R3-MYB family genes of R. palmatum were systematically identified based on full-length transcriptome sequencing followed by bioinformatics analyses. The correlation analysis was carried out by using co-expression analysis, protein interaction analysis, and real-time fluorescence quantitative analysis after MeJA treatment. The RpMYB81 and RpMYB98 genes were amplified by RT-PCR, and their subcellular localization and self-activation characteristics were analyzed.

RESULTS

Comparative transcriptome analysis results revealed a total of 3525 upregulated and 6043 downregulated DEGs in the CK versus MeJA group; 28 DEGs were involved in the anthraquinone pathway. Eleven CHS genes that belonged to the PKS family were differentially expressed and involved in anthraquinone biosynthesis. Twelve differentially expressed MYBs genes were found to be co-expressed and interact with CHS genes. Furthermore, 52 MYB genes were identified as positive regulators of anthraquinone biosynthesis and were further characterized. Three MYB genes including RpMYB81, RpMYB98, and RpMYB100 responded to MeJA treatment in R. palmatum, and the levels of these genes were verified by qRT-PCR. RpMYB81 was related to anthraquinone biosynthesis. RpMYB98 had an interaction with genes in the anthraquinone biosynthesis pathway. RpMYB81 and RpMYB98 were mainly localized in the nucleus. RpMYB81 had self-activation activity, while RpMYB98 had no self-activation activity.

CONCLUSION

RpMYB81, RpMYB98, and RpMYB100 were significantly induced by MeJA treatment. RpMYB81 and RpMYB98 are located in the nucleus, and RpMYB81 has transcriptional activity, suggesting that it might be involved in the transcriptional regulation of anthraquinone biosynthesis in R. palmatum.

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

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

Establishment and validation of a callus tissue transformation system for German chamomile (Matricaria chamomilla L.)

BACKGROUND

German chamomile (Matricaria chamomilla L.) is an important medicinal plant, and the essential oils in the flowers have various biological activities. Genetic transformation systems are important for plant quality improvement and molecular research. To the best of our knowledge, a genetic transformation system has not yet been reported for German chamomile.

RESULTS

In this study, we developed Agrobacterium-mediated transformation protocols for German chamomile callus tissues. This involved optimizing key parameters, such as hygromycin and cefotaxime concentrations, bacterial density, and infection and co-culture durations. We also performed gas chromatography-mass spectrometry analysis to identify volatile compounds in non-transgenic and transgenic callus and hairy root tissues. Furthermore, to compare and verify the callus transformation system of German chamomile, we transferred McFPS to the hairy roots of German chamomile. The results showed that the optimal conditions for Agrobacterium-mediated callus tissue transformation were as follows: explant, petiole; cefotaxime concentration, 300 mg/L; hygromycin concentration, 10 mg/L; and bacterial solution concentration, OD600 = 0.6; callus transformation efficiency was the highest when the co-culture time was 3 days.

CONCLUSIONS

Establishment of a high-efficiency callus transformation system will lay the foundation for gene function identification in German chamomile.

LcMYB106 suppresses monoterpene biosynthesis by negatively regulating LcTPS32 expression in Litsea cubeba

Litsea cubeba, the core species of the Lauraceae family, is valuable for the production of essential oils due to its high concentration of monoterpenes (90%). The key monoterpene synthase and metabolic regulatory network of monoterpene biosynthesis have provided new insights for improving essential oil content. However, there are few studies on the regulation mechanism of monoterpenes in L. cubeba. In this study, we investigated LcTPS32, a member of the TPS-b subfamily, and identified its function as an enzyme for the synthesis of monoterpenes, including geraniol, α-pinene, β-pinene, β-myrcene, linalool and eucalyptol. The quantitative real-time PCR analysis showed that LcTPS32 was highly expressed in the fruits of L. cubeba and contributed to the characteristic flavor of its essential oil. Overexpression of LcTPS32 resulted in a significant increase in the production of monoterpenes in L. cubeba by activating both the MVA and MEP pathways. Additionally, the study revealed that LcMYB106 played a negative regulatory role in monoterpenes biosynthesis by directly binding to the promoter of LcTPS32. Our study indicates that LcMYB106 could serve as a crucial target for metabolic engineering endeavors, aiming at enhancing the monoterpene biosynthesis in L. cubeba.

Contemporary understanding of transcription factor regulation of terpenoid biosynthesis in plants

KEY MESSAGE

This review summarized how TFs function independently or in response to environmental factors to regulate terpenoid biosynthesis via fine-tuning the expression of rate-limiting enzymes. Terpenoids are derived from various species and sources. They are essential for interacting with the environment and defense mechanisms, such as antimicrobial, antifungal, antiviral, and antiparasitic properties. Almost all terpenoids have high medicinal value and economic performance. Recently, the control of enzyme genes on terpenoid biosynthesis has received a great deal of attention, but transcriptional factors regulatory network on terpenoid biosynthesis and accumulation has yet to get a thorough review. Transcription factors function as activators or suppressors independently or in response to environmental stimuli, fine-tuning terpenoid accumulation through regulating rate-limiting enzyme expression. This study investigates the advancements in transcription factors related to terpenoid biosynthesis and systematically summarizes previous works on the specific mechanisms of transcription factors that regulate terpenoid biosynthesis via hormone signal-transcription regulatory networks in plants. This will help us to better comprehend the regulatory network of terpenoid biosynthesis and build the groundwork for terpenoid development and effective utilization.

Analysis of WAK Genes in Nine Cruciferous Species with a Focus on Brassica napus L

The wall-associated kinase family contributes to plant cell elongation and pathogen recognition. Nine Cruciferous species were studied for identification and molecular evolution of the WAK gene family. Firstly, 178 WAK genes were identified. A phylogenetic tree was constructed of the Cruciferous WAK proteins into four categories, of which the Brassica rapa, Brassica oleracea and Brassica napus genes in the U's triangle were more closely related. The WAK gene family was unevenly distributed in B. napus chromosomal imaging, with the largest number of BnWAK genes located on chromosome C08. In the expression analysis, the expression patterns of the WAK gene family varied under different stress treatments, and some members of BnWAKs were significantly different under stress treatments. This study lays a foundation for further revealing the functional mechanisms of the WAK gene family in Brassica napus.

Designing plant flavonoids: harnessing transcriptional regulation and enzyme variation to enhance yield and diversity

Plant synthetic biology has emerged as a powerful and promising approach to enhance the production of value-added metabolites in plants. Flavonoids, a class of plant secondary metabolites, offer numerous health benefits and have attracted attention for their potential use in plant-based products. However, achieving high yields of specific flavonoids remains challenging due to the complex and diverse metabolic pathways involved in their biosynthesis. In recent years, synthetic biology approaches leveraging transcription factors and enzyme diversity have demonstrated promise in enhancing flavonoid yields and expanding their production repertoire. This review delves into the latest research progress in flavonoid metabolic engineering, encompassing the identification and manipulation of transcription factors and enzymes involved in flavonoid biosynthesis, as well as the deployment of synthetic biology tools for designing metabolic pathways. This review underscores the importance of employing carefully-selected transcription factors to boost plant flavonoid production and harnessing enzyme promiscuity to broaden flavonoid diversity or streamline the biosynthetic steps required for effective metabolic engineering. By harnessing the power of synthetic biology and a deeper understanding of flavonoid biosynthesis, future researchers can potentially transform the landscape of plant-based product development across the food and beverage, pharmaceutical, and cosmetic industries, ultimately benefiting consumers worldwide.

Advances in the Biosynthesis of Terpenoids and Their Ecological Functions in Plant Resistance

Secondary metabolism plays an important role in the adaptation of plants to their environments, particularly by mediating bio-interactions and protecting plants from herbivores, insects, and pathogens. Terpenoids form the largest group of plant secondary metabolites, and their biosynthesis and regulation are extremely complicated. Terpenoids are key players in the interactions and defense reactions between plants, microorganisms, and animals. Terpene compounds are of great significance both to plants themselves and the ecological environment. On the one hand, while protecting plants themselves, they can also have an impact on the environment, thereby affecting the evolution of plant communities and even ecosystems. On the other hand, their economic value is gradually becoming clear in various aspects of human life; their potential is enormous, and they have broad application prospects. Therefore, research on terpenoids is crucial for plants, especially crops. This review paper is mainly focused on the following six aspects: plant terpenes (especially terpene volatiles and plant defense); their ecological functions; their biosynthesis and transport; related synthesis genes and their regulation; terpene homologues; and research and application prospects. We will provide readers with a systematic introduction to terpenoids covering the above aspects.

Global dissection of R2R3-MYB in Pogostemon cablin uncovers a species-specific R2R3-MYB clade

MYB family is one of the largest transcription factor families in plants and plays a crucial role in regulating plant biochemical and physiological processes. However, R2R3-MYBs in patchouli have not been systematically investigated. Here, based on the gene annotation of patchouli genome sequence, 484 R2R3-MYB transcripts were detected. Further in-depth analysis of the gene structure and expression of R2R3-MYBs supported the tetraploid hybrid origin of patchouli. When combined with R2R3-MYBs from Arabidopsis, a phylogenetic tree of patchouli R2R3-MYBs was constructed and divided into 31 clades. Interestingly, a patchouli-specific R2R3-MYB clade was found and confirmed by homologous from other Lamiaceae species. The syntenic analysis demonstrated that tandem duplication contributed to its evolution. This study systematically analysed the R2R3-MYB family in patchouli, providing information on its gene characterization, functional prediction, and species evolution.

Lavandula Species, Their Bioactive Phytochemicals, and Their Biosynthetic Regulation

Lavandula species are one of the most useful aromatic and medicinal plants and have great economic potential. The phytopharmaceutical contribution of the secondary metabolites of the species is unquestionable. Most recent studies have been focusing on the elucidation of the genetic background of secondary metabolite production in lavender species. Therefore, knowledge of not only genetic but especially epigenetic mechanisms for the regulation of secondary metabolites is necessary for the modification of those biosynthesis processes and the understanding of genotypic differences in the content and compositional variability of these products. The review discusses the genetic diversity of Lavandula species in relation to the geographic area, occurrence, and morphogenetic factors. The role of microRNAs in secondary-metabolites biosynthesis is described.

The Blinin Accumulation Promoted by CbMYB32 Involved in Conyza blinii Resistance to Nocturnal Low Temperature

Blinin, a unique terpenoid from Conyza blinii (C. blinii), benefits our health even though this is not its primary function. Physiological and ecological studies have found that the great secondary metabolites participate in important biological processes and relate to species evolution, environmental adaptation, and so on. Moreover, our previous studies have shown that the metabolism and accumulation of blinin has a close correspondence with nocturnal low temperature (NLT). To find out the transcriptional regulation linker in the crosstalk between blinin and NLT, RNA-seq, comparative analysis, and co-expression network were performed. The results indicated that CbMYB32 is located in a nucleus without independent transcriptional activation activity and is probably involved in the metabolism of blinin. Furthermore, we compared the silence and overexpression of CbMYB32 with wild C. blinii. Compared with the overexpression and the wildtype, the CbMYB32 silence line lost more than half of the blinin and detected more peroxide under NLT. Finally, as a characteristic secret of C. blinii, it is reasonable to infer that blinin participates in the NLT adaptation mechanism and has contributed to the systematic evolution of C. blinii.

Transcriptional regulatory network of high-value active ingredients in medicinal plants

High-value active ingredients in medicinal plants have attracted research attention because of their benefits for human health, such as the antimalarial artemisinin, anticardiovascular disease tanshinones, and anticancer Taxol and vinblastine. Here, we review how hormones and environmental factors promote the accumulation of active ingredients, thereby providing a strategy to produce high-value drugs at a low cost. Focusing on major hormone signaling events and environmental factors, we review the transcriptional regulatory network mediating biosynthesis of representative active ingredients. In this network, many transcription factors (TFs) simultaneously control multiple synthase genes; thus, understanding the molecular mechanisms affecting transcriptional regulation of active ingredients will be crucial to developing new breeding possibilities.

Jasmonic acid regulates the biosynthesis of medicinal metabolites via the JAZ9-MYB76 complex in Salvia miltiorrhiza

Jasmonic acid (JA) signaling pathway plays an important role in tanshinone and phenolic acid biosynthesis in Salvia miltiorrhiza. However, the specific regulatory mechanism remains largely unclear. Previous work showed that a JASMONATE ZIM-domain (JAZ) protein, SmJAZ9, acted as a repressor of tanshinone production in S. miltiorrhiza. In this study, we revealed that SmJAZ9 reduced both phenolic acid accumulation and related biosynthetic gene expression, confirming that SmJAZ9 also negatively affected phenolic acid biosynthesis. Then, we identified a novel MYB transcription factor, SmMYB76, which interacted with SmJAZ9. SmMYB76 repressed phenolic acid biosynthesis by directly downregulating SmPAL1, Sm4CL2, and SmRAS1. Further investigation demonstrated that JA mediated phenolic acids biosynthesis via SmJAZ9-SmMYB76 complex. Taken together, these findings state the molecular mechanism that SmJAZ9-SmMYB76 regulated phenolic acid biosynthesis at the transcriptional and protein levels, which provided new insights into JA signaling pathway regulating plant metabolism.

SbMYB3 transcription factor promotes root-specific flavone biosynthesis in Scutellaria baicalensis

Scutellaria baicalensis Georgi produces abundant root-specific flavones (RSFs), which provide various benefits to human health. We have elucidated the complete biosynthetic pathways of baicalein and wogonin. However, the transcriptional regulation of flavone biosynthesis in S. baicalensis remains unclear. We show that the SbMYB3 transcription factor functions as a transcriptional activator involved in the biosynthesis of RSFs in S. baicalensis. Yeast one-hybrid and transcriptional activation assays showed that SbMYB3 binds to the promoter of flavone synthase II-2 (SbFNSII-2) and enhances its transcription. In S. baicalensis hairy roots, RNAi of SbMYB3 reduced the accumulation of baicalin and wogonoside, and SbMYB3 knockout decreased the biosynthesis of baicalein, baicalin, wogonin, and wogonoside, whereas SbMYB3 overexpression enhanced the contents of baicalein, baicalin, wogonin, and wogonoside. Transcript profiling by qRT-PCR demonstrated that SbMYB3 activates SbFNSII-2 expression directly, thus leading to more abundant accumulation of RSFs. This study provides a potential target for metabolic engineering of RSFs.

A Novel R2R3-MYB Transcription Factor SbMYB12 Positively Regulates Baicalin Biosynthesis in Scutellaria baicalensis Georgi

Scutellaria baicalensis Georgi is an annual herb from the Scutellaria genus that has been extensively used as a traditional medicine for over 2000 years in China. Baicalin and other flavonoids have been identified as the principal bioactive ingredients. The biosynthetic pathway of baicalin in S. baicalensis has been elucidated; however, the specific functions of R2R3-MYB TF, which regulates baicalin synthesis, has not been well characterized in S. baicalensis to date. Here, a S20 R2R3-MYB TF (SbMYB12), which encodes 263 amino acids with a length of 792 bp, was expressed in all tested tissues (mainly in leaves) and responded to exogenous hormone methyl jasmonate (MeJA) treatment. The overexpression of SbMYB12 significantly promoted the accumulation of flavonoids such as baicalin and wogonoside in S. baicalensis hairy roots. Furthermore, biochemical experiments revealed that SbMYB12 is a nuclear-localized transcription activator that binds to the SbCCL7-4, SbCHI-2, and SbF6H-1 promoters to activate their expression. These results illustrate that SbMYB12 positively regulates the generation of baicalin and wogonoside. In summary, this work revealed a novel S20 R2R3-MYB regulator and enhances our understanding of the transcriptional and regulatory mechanisms of baicalin biosynthesis, as well as sheds new light on metabolic engineering in S. baicalensis.

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.

Transcription Factor SmSPL2 Inhibits the Accumulation of Salvianolic Acid B and Influences Root Architecture

The SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE (SPL) transcription factor play vital roles in plant growth and development. Although 15 SPL family genes have been recognized in the model medical plant Salvia miltiorrhiza Bunge, most of them have not been functionally characterized to date. Here, we performed a careful characterization of SmSPL2, which was expressed in almost all tissues of S. miltiorrhiza and had the highest transcriptional level in the calyx. Meanwhile, SmSPL2 has strong transcriptional activation activity and resides in the nucleus. We obtained overexpression lines of SmSPL2 and rSmSPL2 (miR156-resistant SmSPL2). Morphological changes in roots, including longer length, fewer adventitious roots, decreased lateral root density, and increased fresh weight, were observed in all of these transgenic lines. Two rSmSPL2-overexpressed lines were subjected to transcriptome analysis. Overexpression of rSmSPL2 changed root architectures by inhibiting biosynthesis and signal transduction of auxin, while triggering that of cytokinin. The salvianolic acid B (SalB) concentration was significantly decreased in rSmSPL2-overexpressed lines. Further analysis revealed that SmSPL2 binds directly to the promoters of Sm4CL9, SmTAT1, and SmPAL1 and inhibits their expression. In conclusion, SmSPL2 is a potential gene that efficiently manipulate both root architecture and SalB concentration in S. miltiorrhiza.

The SmMYB36-SmERF6/SmERF115 module regulates the biosynthesis of tanshinones and phenolic acids in salvia miltiorrhiza hairy roots

Tanshinone and phenolic acids are the most important active substances of Salvia miltiorrhiza, and the insight into their transcriptional regulatory mechanisms is an essential process to increase their content in vivo. SmMYB36 has been found to have important regulatory functions in the synthesis of tanshinone and phenolic acid; paradoxically, its mechanism of action in S. miltiorrhiza is not clear. Here, we demonstrated that SmMYB36 functions as a promoter of tanshinones accumulation and a suppressor of phenolic acids through the generation of SmMYB36 overexpressed and chimeric SmMYB36-SRDX (EAR repressive domain) repressor hairy roots in combination with transcriptomic-metabolomic analysis. SmMYB36 directly down-regulate the key enzyme gene of primary metabolism, SmGAPC, up-regulate the tanshinones biosynthesis branch genes SmDXS2, SmGGPPS1, SmCPS1 and down-regulate the phenolic acids biosynthesis branch enzyme gene, SmRAS. Meanwhile, SmERF6, a positive regulator of tanshinone synthesis activating SmCPS1, was up-regulated and SmERF115, a positive regulator of phenolic acid biosynthesis activating SmRAS, was down-regulated. Furthermore, the seven acidic amino acids at the C-terminus of SmMYB36 are required for both self-activating domain and activation of target gene expression. As a consequence, this study contributes to reveal the potential relevance of transcription factors synergistically regulating the biosynthesis of tanshinone and phenolic acid.

Mining of the CULLIN E3 ubiquitin ligase genes in the whole genome of Salvia miltiorrhiza

CULLIN (CUL) proteins are E3 ubiquitin ligases that are involved in a wide variety of biological processes as well as in response to stress in plants. In Salvia miltiorrhiza, CUL genes have not been characterized and its role in plant development, stress response and secondary metabolite synthesis have not been studied. In this study, genome-wide analyses were performed to identify and to predict the structure and function of CUL of S. miltiorrhiza. Eight CUL genes were identified from the genome of S. miltiorrhiza. The CUL genes were clustered into four subgroups according to phylogenetic relationships. The CUL domain was highly conserved across the family of CUL genes. Analysis of cis-acting elements suggested that CUL genes might play important roles in a variety of biological processes, including abscission reaction acid (ABA) processing. To investigate this hypothesis, we treated hairy roots of S. miltiorrhiza with ABA. The expression of CUL genes varied obviously after ABA treatment. Co-expression network results indicated that three CUL genes might be involved in the biosynthesis of phenolic acid or tanshinone. In summary, the mining of the CUL genes in the whole genome of S. miltiorrhiza contribute novel information to the understanding of the CUL genes and its functional roles in plant secondary metabolites, growth and development.

The diurnal emission of floral scent in Oncidium hybrid orchid is controlled by CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) through the direct regulation on terpene synthase

BACKGROUND

To adapt the periodic fluctuation of environmental factors, plants are subtle to monitor the natural variation for the growth and development. The daily activities and physiological functions in coordination with the natural variation are regulated by circadian clock genes. The circadian emission of floral scents is one of the rhythmic physiological activities controlled by circadian clock genes. Here, we study the molecular mechanism of circadian emission pattern of ocimene and linalool compounds in Oncidium Sharry Baby (Onc. SB) orchid.

RESULTS

GC-Mass analysis revealed that Onc. SB periodically emitted ocimene and linalool during 6 to 14 o'clock daily. Terpene synthase, one of the key gene in the terpenoid biosynthetic pathway is expressed in coordination with scent emission. The promoter structure of terpene synthase revealed a circadian binding sequence (CBS), 5'-AGATTTTT-3' for CIRCADIAN CLOCK ASSOCIATED1 (CCA1) transcription factor. EMSA data confirms the binding affinity of CCA1. Transactivation assay further verified that TPS expression is regulated by CCA1. It suggests that the emission of floral scents is controlled by CCA1.

CONCLUSIONS

The work validates that the mechanism of circadian emission of floral scents in Onc. Sharry Baby is controlled by the oscillator gene, CCA1(CIRCADIAN CLOCK ASSOCIATED 1) under light condition. CCA1 transcription factor up-regulates terpene synthase (TPS) by binding on CBS motif, 5'-AGATTTTT-3' of promoter region to affect the circadian emission of floral scents in Onc. SB.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

The sage genome provides insight into the evolutionary dynamics of diterpene biosynthesis gene cluster in plants

The widely cultivated medicinal and ornamental plant sage (Salvia officinalis L.) is an evergreen shrub of the Lamiaceae family, native to the Mediterranean. We assembled a high-quality sage genome of 480 Mb on seven chromosomes, and identified a biosynthetic gene cluster (BGC) encoding two pairs of diterpene synthases (diTPSs) that, together with the cytochromes P450 (CYPs) genes located inside and outside the cluster, form two expression cascades responsible for the shoot and root diterpenoids, respectively, thus extending BGC functionality from co-regulation to orchestrating metabolite production in different organs. Phylogenomic analysis indicates that the Salvia clades diverged in the early Miocene. In East Asia, most Salvia species are herbaceous and accumulate diterpenoids in storage roots. Notably, in Chinese sage S. miltiorrhiza, the diterpene BGC has contracted and the shoot cascade has been lost. Our data provide genomic insights of micro-evolution of growth type-associated patterning of specialized metabolite production in plants.

Comparative proteomics reveals biochemical changes in Salvia miltiorrhiza Bunge during sweating processing

ETHNOPHARMACOLOGICAL RELEVANCE

Salvia miltiorrhiza Bunge is a bulk medicinal material used in traditional Chinese medicine, that can cure cardiovascular diseases, neurasthenia, and other conditions. Sweating is a frequently used method of processing S. miltiorrhiza for medical applications. We previously demonstrated changes to the metabolic profile of linoleic acid, glyoxylate, and dicarboxylate after Sweating. However, this alteration has not been explained at the molecular level.

MATERIALS AND METHODS

Fresh roots of Salvia miltiorrhiza Bunge were treated by the Sweating processing, and then the tandem mass tag technique was used to compare the proteome difference between Sweating S. miltiorrhiza and non-Sweating S. miltiorrhiza.

RESULTS

We identified a total of 850 differentially expressed proteins after Sweating treatment in S. miltiorrhiza, including 529 upregulated proteins and 321 downregulated proteins. GO enrichment analysis indicated that these differentially expressed proteins are involved in external encapsulating structure, cell wall, oxidoreductase activity, ligase activity, and others. Further analysis showed that CYP450, the pathogenesis-related protein Bet v 1 family, and the peroxidase domain were the major protein domains. KEGG enrichment identified 18 pathways, of which phenylpropanoid biosynthesis is the most important one related to the metabolite profile and is the principal chemical component of S. miltiorrhiza.

CONCLUSION

This study addressed potential molecular mechanisms in S. miltiorrhiza after Sweating, and the findings provide reasons for the changes in biochemical properties and metabolites changes which might cause pharmacological variation at the proteome level.

Genome-wide analysis of MYB family genes in Tripterygium wilfordii and their potential roles in terpenoid biosynthesis

Terpenoids are a class of significant bioactive components in the woody vine of Tripterygium wilfordii. Previous studies have shown that MYB transcription factors play important roles in plant secondary metabolism, growth, and developmental processes. However, the MYB involved in terpenoid biosynthesis in Tripterygium wilfordii are unknown. To identify Tripterygium wilfordii MYB (TwMYB) genes that are involved in terpenoid biosynthesis, we conducted the genome-wide analysis of the TwMYB gene family. A total of 207 TwMYBs were identified including 84 1R-TwMYB, 117 R2R3-TwMYB, four 3R-TwMYB, and two 4R-TwMYB genes. The most abundant R2R3-TwMYBs together with their Arabidopsis homologs were categorized into 26 subgroups. Intraspecific collinearity analysis found that the 74.9% of the TwMYBs may be generated by segmental duplication events, and 36.7% of duplicated gene pairs were derived from the specific whole genome duplication (WGD) event in Tripterygium wilfordii. In addition, interspecies collinearity analysis found that 16 TwMYB genes formed homologous gene pairs with MYB genes in seven representative species, which indicated they may have a key role in evolution. Notably, we found that the TwMYB genes were differentially expressed in various tissues by expression pattern analysis. In order to further select the candidate genes related to terpenoid biosynthesis, the assay of Methyl jasmonate (MeJA) induction and analysis of phylogenetic tree was conducted. It was speculated that six candidate TwMYB genes (TwMYB33, TwMYB34, TwMYB45, TwMYB67, TwMYB102, and TwMYB103) are involved in regulating terpenoid biosynthesis. This study is the first systematic analysis of the TwMYB gene family and will lay a foundation for the functional characterization of TwMYB genes in the regulation of terpenoid biosynthesis.

5-Azacytidine increases tanshinone production in Salvia miltiorrhiza hairy roots through epigenetic modulation

Recent studies have indicated strong connections between epigenetic modulation and secondary metabolites in plants. It is vital to understand the roles of epigenetics in the production of secondary metabolites. In this study, the inhibitor of DNA methylation 5-azacytidine (5-Az) was used on the hairy roots of the medicinal plant Salvia miltiorrhiza to investigate its effect on secondary metabolite production, gene expression, methylation levels in genomic DNA and promoter regions. Our results showed that the contents of tanshinones in S. miltiorrhiza hairy roots increased by 1.5–5 times, and some genes in the biosynthesis pathway showed an upward trend. According to our NGS analysis, the methylation pattern in the promotor of the gene encoding copalyl diphosphate synthase (CPS) was altered, and 51 out of 145 cytosines were demethylated during 5-Az treatment. A total of 36 putative transcription factors (TFs) binding cites were identified in these demethylation sites. Among these TFs binding cites, cis-regulatory elements for the binding of NF-Y and MYB were frequently found in our results. This is the first report to demonstrate a possible mechanism of DNA methylation participating in tanshinone biosynthesis in S. miltiorrhiza hairy roots by modulating the CPS promoter and TFs binding sites.

Increasing Expression of PnGAP and PnEXPA4 Provides Insights Into the Enlargement of Panax notoginseng Root Size From Qing Dynasty to Cultivation Era

Root size is a key trait in plant cultivation and can be influenced by the cultivation environment. However, physical evidence of root size change in a secular context is scarce due to the difficulty in preserving ancient root samples, and how they were modified during the domestication and cultivation stays unclear. About 100 ancient root samples of Panax notoginseng, preserved as tribute in the Palace Museum (A.D. 1636 to 1912, Qing dynasty), provided an opportunity to investigate the root size changes during the last 100 years of cultivation. The dry weight of ancient root samples (~120 tou samples, tou represents number of roots per 500 g dry weight) is 0.22-fold of the modern samples with the biggest size (20 tou samples). Transcriptome analysis revealed that PnGAP and PnEXPA4 were highly expressed in 20 tou samples, compared with the 120 tou samples, which might contribute to the thicker cell wall and a higher content of lignin, cellulose, and callose in 20 tou samples. A relatively lower content of dencichine and higher content of ginsenoside Rb1 in 20 tou samples are also consistent with higher expression of ginsenoside biosynthesis-related genes. PnPHL8 was filtrated through transcriptome analysis, which could specifically bind the promoters of PnGAP, PnCYP716A47, and PnGGPPS3, respectively. The results in this study represent the first physical evidence of root size changes in P. notoginseng in the last 100 years of cultivation and contribute to a comprehensive understanding of how the cultivation environment affected root size, chemical composition, and clinical application.

CRISPR/Cas9-Mediated Targeted Mutagenesis of FtMYB45 Promotes Flavonoid Biosynthesis in Tartary Buckwheat (Fagopyrum tataricum)

The clustered regularly interspaced short palindromic repeat/CRISPR-associated protein 9 (CRISPR/Cas9) technology is an efficient genome editing tool used in multiple plant species. However, it has not been applied to Tartary buckwheat (Fagopyrum tataricum), which is an important edible and medicinal crop rich in rutin and other flavonoids. FtMYB45 is an R2R3-type MYB transcription factor that negatively regulates flavonoid biosynthesis in Tartary buckwheat. Here, the CRISPR/Cas9 system polycistronic tRNA-sgRNA (PTG)/Cas9 was employed to knock out the FtMYB45 gene in Tartary buckwheat. Two single-guide RNAs (sgRNAs) were designed to target the second exon of the FtMYB45 gene. Twelve transgenic hairy roots were obtained using Agrobacterium rhizogenes-mediated transformation. Sequencing data revealed that six lines containing six types of mutations at the predicted double-stranded break site were generated using sgRNA1. The mutation frequency reached 50%. A liquid chromatography coupled with triple quadrupole mass spectrometry (LC-QqQ-MS) based metabolomic analysis revealed that the content of rutin, catechin, and other flavonoids was increased in hairy root mutants compared with that of lines transformed with the empty vector. Thus, CRISPR/Cas9-mediated targeted mutagenesis of FtMYB45 effectively increased the flavonoids content of Tartary buckwheat. This finding demonstrated that the CRISPR/Cas9 system is an efficient tool for precise genome editing in Tartary buckwheat and lays the foundation for gene function research and quality improvement in Tartary buckwheat.

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.

Comparative transcriptome analyses of different Salvia miltiorrhiza varieties during the accumulation of tanshinones

Salvia miltiorrhiza (Labiatae) is an important medicinal plant in traditional Chinese medicine. Tanshinones are one of the main active components of S. miltiorrhiza. It has been found that the intraspecific variation of S. miltiorrhiza is relatively large and the content of tanshinones in its roots of different varieties is also relatively different. To investigate the molecular mechanisms that responsible for the differences among these varieties, the tanshinones content was determined and comparative transcriptomics analysis was carried out during the tanshinones accumulation stage. A total of 52,216 unigenes were obtained from the transcriptome by RNA sequencing among which 23,369 genes were differentially expressed among different varieties, and 2,016 genes including 18 diterpenoid biosynthesis-related genes were differentially expressed during the tanshinones accumulation stage. Functional categorization of the differentially expressed genes (DEGs) among these varieties revealed that the pathway related to photosynthesis, oxidative phosphorylation, secondary metabolite biosynthesis, diterpenoid biosynthesis, terpenoid backbone biosynthesis, sesquiterpenoid and triterpenoid biosynthesis are the most differentially regulated processes in these varieties. The six tanshinone components in these varieties showed different dynamic changes in tanshinone accumulation stage. In addition, combined with the analysis of the dynamic changes, 277 DEGs (including one dehydrogenase, three CYP450 and 24 transcription factors belonging to 12 transcription factor families) related to the accumulation of tanshinones components were obtained. Furthermore, the KEGG pathway enrichment analysis of these 277 DEGs suggested that there might be an interconnection between the primary metabolic processes, signaling processes and the accumulation of tanshinones components. This study expands the vision of intraspecific variation and gene regulation mechanism of secondary metabolite biosynthesis pathways in medicinal plants from the “omics” perspective.

Integrative Omics Analyses Reveal the Effects of Copper Ions on Salvianolic Acid Biosynthesis

Salvianolic acids, a group of secondary metabolites produced by Salvia miltiorrhiza, are widely used for treating cerebrovascular diseases. Copper is recognized as a necessary microelement and plays an essential role in plant growth. At present, the effect of copper on the biosynthesis of SalAs is unknown. Here, an integrated metabolomic and transcriptomic approach, coupled with biochemical analyses, was employed to dissect the mechanisms by which copper ions induced the biosynthesis of SalAs. In this study, we identified that a low concentration (5 μM) of copper ions could promote growth of S. miltiorrhiza and the biosynthesis of SalAs. Results of the metabolomics analysis showed that 160 metabolites (90 increased and 70 decreased) were significantly changed in S. miltiorrhiza treated with low concentration of copper ions. The differential metabolites were mainly involved in amino acid metabolism, the pentose phosphate pathway, and carbon fixation in photosynthetic organisms. The contents of chlorophyll a, chlorophyll b, and total chlorophyll were significantly increased in leaves of low concentration of copper-treated S. miltiorrhiza plants. Importantly, core SalA biosynthetic genes (laccases and rosmarinic acid synthase), SalA biosynthesis-related transcription factors (MYBs and zinc finger CCCH domain-containing protein 33), and chloroplast proteins-encoding genes (blue copper protein and chlorophyll-binding protein) were upregulated in the treated samples as indicated by a comprehensive transcriptomic analysis. Bioinformatics and enzyme activity analyses showed that laccase 20 contained copper-binding motifs, and its activity in low concentration of copper ions-treated S. miltiorrhiza was much higher than that in the control. Our results demonstrate that enhancement of copper ions of the accumulation of SalAs might be through regulating laccase 20, MYBs, and zinc finger transcription factors, and photosynthetic genes.

Transcriptome analysis of two Pogostemon cablin chemotypes reveals genes related to patchouli alcohol biosynthesis

Pogostemon cablin, a medicinally and economically important perennial herb, is cultivated around the world due to its medicinal and aromatic properties. Different P. cablin cultivars exhibit different morphological traits and patchouli oil components and contents (especially patchouli alcohol (PA) and pogostone (PO)). According to the signature constituent of the leaf, P. cablin was classified into two different chemotypes, including PA-type and PO-type. To better understand the molecular mechanisms of PA biosynthesis, the transcriptomes of Chinese-cultivated P. cablin cv. PA-type “Nanxiang” (NX) and PO-type “Paixiang” (PX) were analyzed and compared with ribonucleic acid sequencing (RNA-Seq) technology. We obtained a total of 36.83 G clean bases from the two chemotypes, compared them with seven databases and revealed 45,394 annotated unigenes. Thirty-six candidate unigenes participating in the biosynthesis of PA were found in the P. cablin transcriptomes. Overall, 8,390 differentially expressed unigenes were identified between the chemotypes, including 2,467 upregulated and 5,923 downregulated unigenes. Furthermore, six and nine differentially expressed genes (DEGs) were mapped to the terpenoid backbone biosynthetic and sesquiterpenoid and triterpenoid biosynthetic pathways, respectively. One key sesquiterpene synthase gene involved in the sesquiterpenoid and triterpenoid biosynthetic pathways, encoding patchoulol synthase variant 1, was significantly upregulated in NX. Additionally, GC-MS analysis of the two chemotypes in this study showed that the content of PA in NX was significantly higher than that of PX, while the content of PO showed the opposite phenotype. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that the DEG expression tendency was consistent with the transcriptome sequencing results. Overall, 23 AP2/ERF, 13 bHLH, 11 MYB, 11 NAC, three Trihelix, 10 WRKY and three bZIP genes that were differentially expressed may act as regulators of terpenoid biosynthesis. Altogether, 8,314 SSRs were recognized within 6,825 unigenes, with a distribution frequency of 18.32%, among which 1,202 unigenes contained more than one SSR. The transcriptomic characteristics of the two P. cablin chemotypes are comprehensively reported in this study, and these results will contribute to a better understanding of the molecular mechanism of PA biosynthesis. Our transcriptome data also provide a valuable genetic resource for further studies on P. cablin.

Recent Advanced Metabolic and Genetic Engineering of Phenylpropanoid Biosynthetic Pathways

The MYB transcription factors (TFs) are evolving as critical role in the regulation of the phenylpropanoid and tanshinones biosynthetic pathway. MYB TFs relate to a very important gene family, which are involved in the regulation of primary and secondary metabolisms, terpenoids, bioactive compounds, plant defense against various stresses and cell morphology. R2R3 MYB TFs contained a conserved N-terminal domain, but the domain at C-terminal sorts them different regarding their structures and functions. MYB TFs suppressors generally possess particular repressive motifs, such as pdLNLD/ELxiG/S and TLLLFR, which contribute to their suppression role through a diversity of complex regulatory mechanisms. A novel flower specific "NF/YWSV/MEDF/LW" conserved motif has a great potential to understand the mechanisms of flower development. In the current review, we summarize recent advanced progress of MYB TFs on transcription regulation, posttranscriptional, microRNA, conserved motif and propose directions to future prospective research. We further suggest there should be more focus on the investigation for the role of MYB TFs in microalgae, which has great potential for heterologous protein expression system for future perspectives.

R2R3-MYB Transcription Factor SmMYB52 Positively Regulates Biosynthesis of Salvianolic Acid B and Inhibits Root Growth in Salvia miltiorrhiza

The dried root of Salvia miltiorrhiza is a renowned traditional Chinese medicine that was used for over 1000 years in China. Salvianolic acid B (SalB) is the main natural bioactive product of S. miltiorrhiza. Although many publications described the regulation mechanism of SalB biosynthesis, few reports simultaneously focused on S. miltiorrhiza root development. For this study, an R2R3-MYB transcription factor gene (SmMYB52) was overexpressed and silenced, respectively, in S. miltiorrhiza sterile seedlings. We found that SmMYB52 significantly inhibited root growth and indole-3-acetic acid (IAA) accumulation, whereas it activated phenolic acid biosynthesis and the jasmonate acid (JA) signaling pathway. Quantitative real-time polymerase chain reaction (qRT-PCR) analyses revealed that SmMYB52 suppressed the transcription levels of key enzyme-encoding genes involved in the IAA biosynthetic pathway and activated key enzyme-encoding genes involved in the JA and phenolic acid biosynthesis pathways. In addition, yeast one-hybrid (Y1H) and dual-luciferase assay showed that SmMYB52 directly binds to and activates the promoters of several key enzyme genes for SalB biosynthesis, including SmTAT1, Sm4CL9, SmC4H1, and SmHPPR1, to promote the accumulation of SalB. This is the first report of a regulator that simultaneously affects root growth and the production of phenolic acids in S. miltiorrhiza.

Comparative transcriptome analysis provides insights into steviol glycoside synthesis in stevia (Stevia rebaudiana Bertoni) leaves under nitrogen deficiency

Transcriptome analysis revealed the potential mechanism of nitrogen regulating steviol glycosides synthesis via shifting of leaf carbon metabolic flux or inducing certain transcription factors. Nitrogen (N) plays key regulatory roles in both stevia (Stevia rebaudiana) growth and the synthesis of its functional metabolite steviol glycosides (SGs), but the mechanism by which this nutrient regulates SGs synthesis remains to be elucidated. To address this question, a pot experiment was performed in a greenhouse where stevia plants fertilized with N (the control as CK plants) and compared with plants without the supply of N. Physiological and biochemical analyses were conducted to test the growth and metabolic responses of plants to N regimes. Our results showed that N deficiency significantly inhibited plant growth and leaf photosynthesis, while increased leaf SGs contents in stevia (49.97, 46.64 and 84.80% respectively for rebaudioside A, stevioside, and rebaudioside C), which may be partly due to "concentration effect". Then, transcriptome analysis was conducted to understand the underlying mechanisms. A total of 535 differentially expressed genes were identified, and carbon metabolism-related events were highlighted by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. Many of these genes were significantly upregulated by N-deficiency, including those involved in "phenylpropanoid biosynthesis", "flavonoid biosynthesis" and "starch and sucrose metabolism". Our study also analyzed the expression patterns of SGs synthesis-related genes under two N regimes and the potential transcription factors linking N nutrition and SG metabolism. N-deficiency may promote SGs synthesis by changing the carbon metabolism flux or inducing certain transcription factors. Our results provide deeper insight into the relationship between N nutrition and SGs synthesis in stevia plants.

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.

Transcriptome Analysis of Light-Regulated Monoterpenes Biosynthesis in Leaves of Mentha canadensis L

Light is a key environmental aspect that regulates secondary metabolic synthesis. The essential oil produced in mint (Mentha canadensis L.) leaves is used widely in the aromatics industry and in medicine. Under low-light treatment, significant reductions in peltate glandular trichome densities were observed. GC-MS analysis showed dramatically reduced essential oil and menthol contents. Light affected the peltate glandular trichomes’ development and essential oil yield production. However, the underlying mechanisms of this regulation were elusive. To identify the critical genes during light-regulated changes in oil content, following a 24 h darkness treatment and a 24 h recovery light treatment, leaves were collected for transcriptome analysis. A total of 95,579 unigenes were obtained, with an average length of 754 bp. About 56.58% of the unigenes were annotated using four public protein databases: 10,977 differentially expressed genes (DEGs) were found to be involved in the light signaling pathway and monoterpene synthesis pathway. Most of the TPs showed a similar expression pattern: downregulation after darkness treatment and upregulation after the return of light. In addition, the genes involved in the light signal transduction pathway were analyzed. A series of responsive transcription factors (TFs) were identified and could be used in metabolic engineering as an effective strategy for increasing essential oil yields.

Plasma Membrane H+-ATPase SmPHA4 Negatively Regulates the Biosynthesis of Tanshinones in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge has been widely used in the treatment of cardiovascular and cerebrovascular diseases, due to the pharmacological action of its active components such as the tanshinones. Plasma membrane (PM) H⁺-ATPase plays key roles in numerous physiological processes in plants. However, little is known about the PM H⁺-ATPase gene family in S. miltiorrhiza (Sm). Here, nine PM H⁺-ATPase isoforms were identified and named SmPHA1-SmPHA9. Phylogenetic tree analysis showed that the genetic distance of SmPHAs was relatively far in the S. miltiorrhiza PM H⁺-ATPase family. Moreover, the transmembrane structures were rich in SmPHA protein. In addition, SmPHA4 was found to be highly expressed in roots and flowers. HPLC revealed that accumulation of dihydrotanshinone (DT), cryptotanshinone (CT), and tanshinone I (TI) was significantly reduced in the SmPHA4-OE lines but was increased in the SmPHA4-RNAi lines, ranging from 2.54 to 3.52, 3.77 to 6.33, and 0.35 to 0.74 mg/g, respectively, suggesting that SmPHA4 is a candidate regulator of tanshinone metabolites. Moreover, qRT-PCR confirmed that the expression of tanshinone biosynthetic-related key enzymes was also upregulated in the SmPHA4-RNAi lines. In summary, this study highlighted PM H⁺-ATPase function and provided new insights into regulatory candidate genes for modulating secondary metabolism biosynthesis in S. miltiorrhiza.

Transcription Factor: A Powerful Tool to Regulate Biosynthesis of Active Ingredients in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge is a common Chinese herbal medicine, and its major active ingredients are phenolic acids and tanshinones, which are widely used to treat vascular diseases. However, the wild form of S. miltiorrhiza possess low levels of these important pharmaceutical agents; thus, improving their levels is an active area of research. Transcription factors, which promote or inhibit the expressions of multiple genes involved in one or more biosynthetic pathways, are powerful tools for controlling gene expression in biosynthesis. Several families of transcription factors have been reported to participate in regulating phenolic acid and tanshinone biosynthesis and influence their accumulation. This review summarizes the current status in this field, with focus on the transcription factors which have been identified in recent years and their functions in the biosynthetic regulation of phenolic acids and tanshinones. Otherwise, the new insight for further research is provided. Finally, the application of the biosynthetic regulation of active ingredients by the transcription factors in S. miltiorrhiza are discussed, and new insights for future research are explored.

The methyl jasmonate-responsive transcription factor SmMYB1 promotes phenolic acid biosynthesis in Salvia miltiorrhiza

Water-soluble phenolic acids are major bioactive compounds in the medicinal plant species Salvia miltiorrhiza. Phenolic acid biosynthesis is induced by methyl jasmonate (MeJA) in this important Chinese herb. Here, we investigated the mechanism underlying this induction by analyzing a transcriptome library of S. miltiorrhiza in response to MeJA. Global transcriptome analysis identified the MeJA-responsive R2R3-MYB transcription factor-encoding gene SmMYB1. Overexpressing SmMYB1 significantly promoted phenolic acid accumulation and upregulated the expression of genes encoding key enzymes in the phenolic acid biosynthesis pathway, including cytochrome P450-dependent monooxygenase (CYP98A14). Dual-luciferase (dual-LUC) assays and/or an electrophoretic mobility shift assays (EMSAs) indicated that SmMYB1 activated the expression of CYP98A14, as well as the expression of genes encoding anthocyanin biosynthesis pathway enzymes, including chalcone isomerase (CHI) and anthocyanidin synthase (ANS). In addition, SmMYB1 was shown to interact with SmMYC2 to additively promote CYP98A14 expression compared to the action of SmMYB1 alone. Taken together, these results demonstrate that SmMYB1 is an activator that improves the accumulation of phenolic acids and anthocyanins in S. miltiorrhiza. These findings lay the foundation for in-depth studies of the molecular mechanism underlying MeJA-mediated phenolic acid biosynthesis and for the metabolic engineering of bioactive ingredients in S. miltiorrhiza.

The methyl jasmonate-responsive transcription factor SmMYB1 promotes phenolic acid biosynthesis in Salvia miltiorrhiza

Water-soluble phenolic acids are major bioactive compounds in the medicinal plant species Salvia miltiorrhiza. Phenolic acid biosynthesis is induced by methyl jasmonate (MeJA) in this important Chinese herb. Here, we investigated the mechanism underlying this induction by analyzing a transcriptome library of S. miltiorrhiza in response to MeJA. Global transcriptome analysis identified the MeJA-responsive R2R3-MYB transcription factor-encoding gene SmMYB1. Overexpressing SmMYB1 significantly promoted phenolic acid accumulation and upregulated the expression of genes encoding key enzymes in the phenolic acid biosynthesis pathway, including cytochrome P450-dependent monooxygenase (CYP98A14). Dual-luciferase (dual-LUC) assays and/or an electrophoretic mobility shift assays (EMSAs) indicated that SmMYB1 activated the expression of CYP98A14, as well as the expression of genes encoding anthocyanin biosynthesis pathway enzymes, including chalcone isomerase (CHI) and anthocyanidin synthase (ANS). In addition, SmMYB1 was shown to interact with SmMYC2 to additively promote CYP98A14 expression compared to the action of SmMYB1 alone. Taken together, these results demonstrate that SmMYB1 is an activator that improves the accumulation of phenolic acids and anthocyanins in S. miltiorrhiza. These findings lay the foundation for in-depth studies of the molecular mechanism underlying MeJA-mediated phenolic acid biosynthesis and for the metabolic engineering of bioactive ingredients in S. miltiorrhiza.

Lipopolysaccharide Enhances Tanshinone Biosynthesis via a Ca2+-Dependent Manner in Salvia miltiorrhiza Hairy Roots

Tanshinones, the major bioactive components in Salvia miltiorrhiza Bunge (Danshen), are synthesized via the mevalonic acid (MVA) pathway or the 2-C-methyl-D-erythritol-4-phosphate (MEP) pathway and the downstream biosynthesis pathway. In this study, the bacterial component lipopolysaccharide (LPS) was utilized as a novel elicitor to induce the wild type hairy roots of S. miltiorrhiza. HPLC analysis revealed that LPS treatment resulted in a significant accumulation of cryptotanshinone (CT) and dihydrotanshinone I (DTI). qRT-PCR analysis confirmed that biosynthesis genes such as SmAACT and SmHMGS from the MVA pathway, SmDXS and SmHDR from the MEP pathway, and SmCPS, SmKSL and SmCYP76AH1 from the downstream pathway were markedly upregulated by LPS in a time-dependent manner. Furthermore, transcription factors SmWRKY1 and SmWRKY2, which can activate the expression of SmDXR, SmDXS and SmCPS, were also increased by LPS. Since Ca²⁺ signaling is essential for the LPS-triggered immune response, Ca²⁺ channel blocker LaCl₃ and CaM antagonist W-7 were used to investigate the role of Ca²⁺ signaling in tanshinone biosynthesis. HPLC analysis demonstrated that both LaCl₃ and W-7 diminished LPS-induced tanshinone accumulation. The downstream biosynthesis genes including SmCPS and SmCYP76AH1 were especially regulated by Ca²⁺ signaling. To summarize, LPS enhances tanshinone biosynthesis through SmWRKY1- and SmWRKY2-regulated pathways relying on Ca²⁺ signaling. Ca²⁺ signal transduction plays a key role in regulating tanshinone biosynthesis 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.

SmMYB2 promotes salvianolic acid biosynthesis in the medicinal herb Salvia miltiorrhiza

MYB transcription factors play vital roles in plant growth and metabolism. The phytohormone methyl jasmonate (MeJA) promotes phenolic acid accumulation in the medicinal herb Salvia miltiorrhiza, but the regulatory mechanism is poorly understood. Here, we identified the MeJA-responsive R2R3-MYB transcription factor gene SmMYB2 from a transcriptome library produced from MeJA-treated S. miltiorrhiza hairy roots. SmMYB2 expression was tightly correlated with the expression of key salvianolic acid biosynthetic genes including CYP98A14. SmMYB2 was highly expressed in the periderm of S. miltiorrhiza and SmMYB2 localized to the nucleus. Overexpressing SmMYB2 in S. miltiorrhiza hairy roots significantly increased the levels of salvianolic acids (including rosmarinic acid and salvianolic acid B) by upregulating salvianolic acid biosynthetic genes such as CYP98A14. SmMYB2 binds to the MYB-binding motifs in the promoter of CYP98A14, as confirmed by a dual-luciferase assay and electrophoretic mobility shift assays. Anthocyanin contents were significantly higher in SmMYB2-overexpressing hairy root lines than the control, primarily due to the increased expression of CHI, DFR, and ANS. These findings reveal the novel regulatory role of SmMYB2 in MeJA-mediated phenolic acid biosynthesis, providing a useful target gene for metabolic engineering and shedding light on the salvianolic acid regulatory network.

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.