Expansion within the CYP71D subfamily drives the heterocyclization of tanshinones synthesis 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.

High-quality chromosome-level genome assembly and multi-omics analysis of rosemary (Salvia rosmarinus) reveals new insights into the environmental and genome adaptation

High-quality genome of rosemary (Salvia rosmarinus) represents a valuable resource and tool for understanding genome evolution and environmental adaptation as well as its genetic improvement. However, the existing rosemary genome did not provide insights into the relationship between antioxidant components and environmental adaptability. In this study, by employing Nanopore sequencing and Hi-C technologies, a total of 1.17 Gb (97.96%) genome sequences were mapped to 12 chromosomes with 46 121 protein-coding genes and 1265 non-coding RNA genes. Comparative genome analysis reveals that rosemary had a closely genetic relationship with Salvia splendens and Salvia miltiorrhiza, and it diverged from them approximately 33.7 million years ago (MYA), and one whole-genome duplication occurred around 28.3 MYA in rosemary genome. Among all identified rosemary genes, 1918 gene families were expanded, 35 of which are involved in the biosynthesis of antioxidant components. These expanded gene families enhance the ability of rosemary adaptation to adverse environments. Multi-omics (integrated transcriptome and metabolome) analysis showed the tissue-specific distribution of antioxidant components related to environmental adaptation. During the drought, heat and salt stress treatments, 36 genes in the biosynthesis pathways of carnosic acid, rosmarinic acid and flavonoids were up-regulated, illustrating the important role of these antioxidant components in responding to abiotic stresses by adjusting ROS homeostasis. Moreover, cooperating with the photosynthesis, substance and energy metabolism, protein and ion balance, the collaborative system maintained cell stability and improved the ability of rosemary against harsh environment. This study provides a genomic data platform for gene discovery and precision breeding in rosemary. Our results also provide new insights into the adaptive evolution of rosemary and the contribution of antioxidant components in resistance to harsh environments.

Metabolomic and transcriptomic analyses highlight metabolic regulatory networks of Salvia miltiorrhiza in response to replant disease

BACKGROUND

Salvia miltiorrhiza, a well-known traditional Chinese medicine, frequently suffers from replant diseases that adversely affect its quality and yield. To elucidate S. miltiorrhiza's metabolic adaptations to replant disease, we analyzed its metabolome and transcriptome, comparing normal and replant diseased plants for the first time.

RESULTS

We identified 1,269 metabolites, 257 of which were differentially accumulated metabolites, and identified 217 differentially expressed genes. Integrated transcriptomic and metabolomic analyses revealed a significant up-regulation and co-expression of metabolites and genes associated with plant hormone signal transduction and flavonoid biosynthesis pathways in replant diseases. Within plant hormone signal transduction pathway, plants afflicted with replant disease markedly accumulated indole-3-acetic acid and abscisic acid, correlating with high expression of their biosynthesis-related genes (SmAmidase, SmALDH, SmNCED, and SmAAOX3). Simultaneously, changes in hormone concentrations activated plant hormone signal transduction pathways. Moreover, under replant disease, metabolites in the local flavonoid metabolite biosynthetic pathway were significantly accumulated, consistent with the up-regulated gene (SmHTC1 and SmHTC2). The qRT-PCR analysis largely aligned with the transcriptomic results, confirming the trends in gene expression. Moreover, we identified 10 transcription factors co-expressed with differentially accumulated metabolites.

CONCLUSIONS

Overall, we revealed the key genes and metabolites of S. miltiorrhiza under replant disease, establishing a robust foundation for future inquiries into the molecular responses to combat replant stress.

The genomes of seven economic Caesalpinioideae trees provide insights into polyploidization history and secondary metabolite biosynthesis

The Caesalpinioideae subfamily contains many well-known trees that are important for economic sustainability and human health, but a lack of genomic resources has hindered their breeding and utilization. Here, we present chromosome-level reference genomes for the two food and industrial trees Gleditsia sinensis (921 Mb) and Biancaea sappan (872 Mb), the three shade and ornamental trees Albizia julibrissin (705 Mb), Delonix regia (580 Mb), and Acacia confusa (566 Mb), and the two pioneer and hedgerow trees Leucaena leucocephala (1338 Mb) and Mimosa bimucronata (641 Mb). Phylogenetic inference shows that the mimosoid clade has a much higher evolutionary rate than the other clades of Caesalpinioideae. Macrosynteny comparison suggests that the fusion and breakage of an unstable chromosome are responsible for the difference in basic chromosome number (13 or 14) for Caesalpinioideae. After an ancient whole-genome duplication (WGD) shared by all Caesalpinioideae species (CWGD, ∼72.0 million years ago [MYA]), there were two recent successive WGD events, LWGD-1 (16.2-19.5 MYA) and LWGD-2 (7.1-9.5 MYA), in L. leucocephala. Thereafter, ∼40% gene loss and genome-size contraction have occurred during the diploidization process in L. leucocephala. To investigate secondary metabolites, we identified all gene copies involved in mimosine metabolism in these species and found that the abundance of mimosine biosynthesis genes in L. leucocephala largely explains its high mimosine production. We also identified the set of all potential genes involved in triterpenoid saponin biosynthesis in G. sinensis, which is more complete than that based on previous transcriptome-derived unigenes. Our results and genomic resources will facilitate biological studies of Caesalpinioideae and promote the utilization of valuable secondary metabolites.

Genome-Wide Identification and Functional Analysis of Salvia miltiorrhiza MicroRNAs Reveal the Negative Regulatory Role of Smi-miR159a in Phenolic Acid Biosynthesis

MicroRNAs (miRNAs) are a group of endogenous small non-coding RNAs in plants. They play critical functions in various biological processes during plant growth and development. Salvia miltiorrhiza is a well-known traditional Chinese medicinal plant with significant medicinal, economic, and academic values. In order to elucidate the role of miRNAs in S. miltiorrhiza, six small RNA libraries from mature roots, young roots, stems, mature leaves, young leaves and flowers of S. miltiorrhiza and one degradome library from mixed tissues were constructed. A total of 184 miRNA precursors, generating 137 known and 49 novel miRNAs, were genome-widely identified. The identified miRNAs were predicted to play diversified regulatory roles in plants through regulating 891 genes. qRT-PCR and 5' RLM-RACE assays validated the negative regulatory role of smi-miR159a in SmMYB62, SmMYB78, and SmMYB80. To elucidate the function of smi-miR159a in bioactive compound biosynthesis, smi-miR159a transgenic hairy roots were generated and analyzed. The results showed that overexpression of smi-miR159a caused a significant decrease in rosmarinic acid and salvianolic acid B contents. qRT-PCR analysis showed that the targets of smi-miR159a, including SmMYB62, SmMYB78, and SmMYB80, were significantly down-regulated, accompanied by the down-regulation of SmPAL1, SmC4H1, Sm4CL1, SmTAT1, SmTAT3, SmHPPR1, SmRAS, and SmCYP98A14 genes involved in phenolic acid biosynthesis. It suggests that smi-miR159a is a significant negative regulator of phenolic acid biosynthesis in S. miltiorrhiza.

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.

Genome-Wide Identification and Expression Analysis of Sucrose Nonfermenting 1-Related Protein Kinase (SnRK) Genes in Salvia miltiorrhiza in Response to Hormone

The SnRK gene family is the chief component of plant stress resistance and metabolism through activating the phosphorylation of downstream proteins. S. miltiorrhiza is widely used for the treatment of cardiovascular diseases in Asian countries. However, information about the SnRK gene family of S. miltiorrhiza is not clear. The aim of this study is to comprehensively analyze the SnRK gene family of S. miltiorrhiza and its response to phytohormone. Here, 33 SmSnRK genes were identified and divided into three subfamilies (SmSnRK1, SmSnRK2 and SmSnRK3) according to phylogenetic analysis and domain. SmSnRK genes within same subgroup shared similar protein motif composition and were unevenly distributed on eight chromosomes of S. miltiorrhiza. Cis-acting element analysis showed that the promoter of SmSnRK genes was enriched with ABRE motifs. Expression pattern analysis revealed that SmSnRK genes were preferentially expressed in leaves and roots. Most SmSnRK genes were induced by ABA and MeJA treatment. Correlation analysis showed that SmSnRK3.15 and SmSnRK3.18 might positively regulate tanshinone biosynthesis; SmSnRK3.10 and SmSnRK3.12 might positively regulate salvianolic acid biosynthesis. RNAi-based silencing of SmSnRK2.6 down-regulated the biosynthesis of tanshinones and biosynthetic genes expression. An in vitro phosphorylation assay verified that SmSnRK2.2 interacted with and phosphorylated SmAREB1. These findings will provide a valuable basis for the functional characterization of SmSnRK genes and quality improvement of S. miltiorrhiza.

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

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

Characterization of CYP82 genes involved in the biosynthesis of structurally diverse benzylisoquinoline alkaloids in Corydalis yanhusuo

Benzylisoquinoline alkaloids (BIAs) represent a significant class of secondary metabolites with crucial roles in plant physiology and substantial potential for clinical applications. CYP82 genes are involved in the formation and modification of various BIA skeletons, contributing to the structural diversity of compounds. In this study, Corydalis yanhusuo, a traditional Chinese medicine rich in BIAs, was investigated to identify the catalytic function of CYP82s during BIA formation. Specifically, 20 CyCYP82-encoding genes were cloned, and their functions were identified in vitro. Ten of these CyCYP82s were observed to catalyze hydroxylation, leading to the formation of protopine and benzophenanthridine scaffolds. Furthermore, the correlation between BIA accumulation and the expression of CyCYP82s in different tissues of C. yanhusuo was assessed their. The identification and characterization of CyCYP82s provide novel genetic elements that can advance the synthetic biology of BIA compounds such as protopine and benzophenanthridine, and offer insights into the biosynthesis of BIAs with diverse structures in C. yanhusuo.

Molecular regulation of the key specialized metabolism pathways in medicinal plants

The basis of modern pharmacology is the human ability to exploit the production of specialized metabolites from medical plants, for example, terpenoids, alkaloids, and phenolic acids. However, in most cases, the availability of these valuable compounds is limited by cellular or organelle barriers or spatio-temporal accumulation patterns within different plant tissues. Transcription factors (TFs) regulate biosynthesis of these specialized metabolites by tightly controlling the expression of biosynthetic genes. Cutting-edge technologies and/or combining multiple strategies and approaches have been applied to elucidate the role of TFs. In this review, we focus on recent progress in the transcription regulation mechanism of representative high-value products and describe the transcriptional regulatory network, and future perspectives are discussed, which will help develop high-yield plant resources.

Genome-Wide Identification and Characterization of the Salvia miltiorrhiza Histone Deacetylase (HDAC) Family in Response to Multiple Abiotic Stresses

Salvia miltiorrhiza is a plant commonly used in traditional Chinese medicine. Its material bases for treating diseases are tanshinones and phenolic acids, including salvianolic acids. Histone deacetylase proteins (HDACs) are a class of specific functional enzymes that interact with acetylation groups on the N-terminal lysine of histone proteins further regulate gene transcription through structural changes at the chromatin level. HDACs involved in the growth and development of various plants, and induced by plant hormones to regulate the internal environment of plants to resist stress, at the same time affect the accumulation of some secondary metabolites. However, the role of SmHDACs on the accumulation of salvianolic acid in S. miltiorrhiza remains unclear. In this study, 16 SmHDACs genes were identified from the high-quality S. miltiorrhiza genome, their physicochemical properties were predicted. In phylogenetic trees co-constructed with HDACs proteins from other plants, SmHDACs was divided into three subfamilies, each with similar motif and conserved domain composition. The distribution of the three subfamilies is similar to that of dicotyledonous plants. Chromosome localization analysis showed that SmHDACs genes were randomly located. Cis-acting element analysis predicted that SmHDACs gene expression may be related to and induced by various phytohormones, such as MeJA and ABA. By combining the expression pattern and co-expression network induced by phytohormones, we speculate that SmHDACs may further influence the synthesis of salvianolic acid, and identified SmHDA5, a potential functional gene, then speculate its downstream target based on the co-expression network. In summary, we analyzed the SmHDACs gene family of S. miltiorrhiza and screened out the potential functional gene SmHDA5. From the perspective of epigenetics, we proposed the molecular mechanism of plant hormone promoting salvianolic acid synthesis, which filled the gap in the subdivision of histone deacetylase in S. miltiorrhiza research, provided a theoretical basis for the culture and transformation of S. miltiorrhiza germplasm resources.

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.

IMP: bridging the gap for medicinal plant genomics

Medicinal plants have garnered significant attention in ethnomedicine and traditional medicine due to their potential antitumor, anti-inflammatory and antioxidant properties. Recent advancements in genome sequencing and synthetic biology have revitalized interest in natural products. Despite the availability of sequenced genomes and transcriptomes of these plants, the absence of publicly accessible gene annotations and tabular formatted gene expression data has hindered their effective utilization. To address this pressing issue, we have developed IMP (Integrated Medicinal Plantomics), a freely accessible platform at https://www.bic.ac.cn/IMP. IMP curated a total of 8 565 672 genes for 84 high-quality genome assemblies, and 2156 transcriptome sequencing samples encompassing various organs, tissues, developmental stages and stimulations. With the integrated 10 analysis modules, users could simply examine gene annotations, sequences, functions, distributions and expressions in IMP in a one-stop mode. We firmly believe that IMP will play a vital role in enhancing the understanding of molecular metabolic pathways in medicinal plants or plants with medicinal benefits, thereby driving advancements in synthetic biology, and facilitating the exploration of natural sources for valuable chemical constituents like drug discovery and drug production.

Multiomics analyses of two Leonurus species illuminate leonurine biosynthesis and its evolution

The Lamiaceae family is renowned for its terpenoid-based medicinal components, but Leonurus, which has traditional medicinal uses, stands out for its alkaloid-rich composition. Leonurine, the principal active compound found in Leonurus, has demonstrated promising effects in reducing blood lipids and treating strokes. However, the biosynthetic pathway of leonurine remains largely unexplored. Here, we present the chromosome-level genome sequence assemblies of Leonurus japonicus, known for its high leonurine production, and Leonurus sibiricus, characterized by very limited leonurine production. By integrating genomics, RNA sequencing, metabolomics, and enzyme activity assay data, we constructed the leonurine biosynthesis pathway and identified the arginine decarboxylase (ADC), uridine diphosphate glucosyltransferase (UGT), and serine carboxypeptidase-like (SCPL) acyltransferase enzymes that catalyze key reactions in this pathway. Further analyses revealed that the UGT-SCPL gene cluster evolved by gene duplication in the ancestor of Leonurus and neofunctionalization of SCPL in L. japonicus, which contributed to the accumulation of leonurine specifically in L. japonicus. Collectively, our comprehensive study illuminates leonurine biosynthesis and its evolution in Leonurus.

Cytotoxic diterpenoids from Salvia glutinosa and comparison with the tanshinone profile of danshen (Salvia miltiorrhiza)

The roots of Salvia miltiorrhiza are the source of the traditional Chinese medicine danshen and the class of tanshinones, particular quinoid nor-diterpenoids of the abietane type. Of these compounds, cryptotanshinone, dihydrotanshinone I, tanshinone I, and tanshinone IIA, have been extensively studied for their anticancer potential, not only but as well because of their high abundance in S. miltiorrhiza and their thus easy availability. However, also additional Salvia species are known to contain tanshinones, mainly such of the subgenus Glutinaria, of which S. glutinosa is the only species widely occurring in Europe. Using UHPLC-DAD-MS, the tanshinone profile of S. glutinosa roots collected from two different locations was compared to the profile in S. miltiorrhiza roots. In addition, tanshinone IIA and another six diterpenoids from S. glutinosa were investigated for their antiproliferative and cytotoxic potential against MDA-MB-231 and HL-60 cells. Apart from dihydrotanshinone I, which has been previously characterized due to its anticancer properties, we determined danshenol A as a highly antiproliferative and cytotoxic agent, significantly surpassing the effects of dihydrotanshinone I. With regard to the diterpenoid profile, S. miltiorrhiza showed a higher concentration for most of the tanshinones, except for (+)-danshexinkun A, which was present in comparable amounts in both species. Danshenol A, in contrast, was only present in S. glutinosa as were dehydroabietic acid and (+)-pisiferic acid. The results of our study underlines the long traditional use of danshen due to its high amount on tanshinones, but also demonstrates the potential value of investigating closely related species for the discovery of new biologically active lead compounds.

Functional characterization of CYP81C16 involved in the tanshinone biosynthetic pathway in Salvia miltiorrhiza

Danshen, the dried roots and rhizomes of Salvia miltiorrhiza Bunge (S. miltiorrhiza), is widely used in the treatment of cardiovascular and cerebrovascular diseases. Tanshinones, the bioactive compounds from Danshen, exhibit a wide spectrum of pharmacological properties, suggesting their potential for future therapeutic applications. Tanshinone biosynthesis is a complex process involving at least six P450 enzymes that have been identified and characterized, most of which belong to the CYP76 and CYP71 families. In this study, CYP81C16, a member of the CYP71 clan, was identified in S. miltiorrhiza. An in vitro assay revealed that it could catalyze the hydroxylation of four para-quinone-type tanshinones, namely neocryptotanshinone, deoxyneocryptotanshinone, and danshenxinkuns A and B. SmCYP81C16 emerged as a potential broad-spectrum oxidase targeting the C-18 position of para-quinone-type tanshinones with an impressive relative conversion rate exceeding 90%. Kinetic evaluations andin vivo assays underscored its highest affinity towards neocryptotanshinone among the tested substrates. The overexpression of SmCYP81C16 promoted the accumulation of (iso)tanshinone in hairy root lines. The characterization of SmCYP81C16 in this study accentuates its potential as a pivotal tool in the biotechnological production of tanshinones, either through microbial or plant metabolic engineering.

Functions of Representative Terpenoids and Their Biosynthesis Mechanisms in Medicinal Plants

Terpenoids are the broadest and richest group of chemicals obtained from plants. These plant-derived terpenoids have been extensively utilized in various industries, including food and pharmaceuticals. Several specific terpenoids have been identified and isolated from medicinal plants, emphasizing the diversity of biosynthesis and specific functionality of terpenoids. With advances in the technology of sequencing, the genomes of certain important medicinal plants have been assembled. This has improved our knowledge of the biosynthesis and regulatory molecular functions of terpenoids with medicinal functions. In this review, we introduce several notable medicinal plants that produce distinct terpenoids (e.g., Cannabis sativa, Artemisia annua, Salvia miltiorrhiza, Ginkgo biloba, and Taxus media). We summarize the specialized roles of these terpenoids in plant-environment interactions as well as their significance in the pharmaceutical and food industries. Additionally, we highlight recent findings in the fields of molecular regulation mechanisms involved in these distinct terpenoids biosynthesis, and propose future opportunities in terpenoid research, including biology seeding, and genetic engineering in medicinal plants.

Advanced application of nanotechnology in active constituents of Traditional Chinese Medicines

Traditional Chinese Medicines (TCMs) have been used for centuries for the treatment and management of various diseases. However, their effective delivery to targeted sites may be a major challenge due to their poor water solubility, low bioavailability, and potential toxicity. Nanocarriers, such as liposomes, polymeric nanoparticles, inorganic nanoparticles and organic/inorganic nanohybrids based on active constituents from TCMs have been extensively studied as a promising strategy to improve the delivery of active constituents from TCMs to achieve a higher therapeutic effect with fewer side effects compared to conventional formulations. This review summarizes the recent advances in nanocarrier-based delivery systems for various types of active constituents of TCMs, including terpenoids, polyphenols, alkaloids, flavonoids, and quinones, from different natural sources. This review covers the design and preparation of nanocarriers, their characterization, and in vitro/vivo evaluations. Additionally, this review highlights the challenges and opportunities in the field and suggests future directions for research. Nanocarrier-based delivery systems have shown great potential in improving the therapeutic efficacy of TCMs, and this review may serve as a comprehensive resource to researchers in this field.

(±)-Salvicatone A: A Pair of C27-Meroterpenoid Enantiomers with Skeletons from the Roots and Rhizomes of Salvia castanea Diels f. tomentosa Stib

(±)-Salvicatone A (1), a C27-meroterpenoid featuring a unique 6/6/6/6/6-pentacyclic carbon skeleton with a 7,8,8a,9,10,10a-hexahydropyren-1 (6H)-one motif, was isolated from the roots and rhizomes of Salvia castanea Diels f. tomentosa Stib. Its structure was characterized by comprehensive spectroscopic analyses along with computer-assisted structure elucidation, including ACD/structure elucidator and quantum chemical calculations with 1H/13C NMR and electronic circular dichroism. Biogenetically, compound 1 was constructed from decarboxylation following [4 + 2] Diels-Alder cycloaddition reaction between caffeic acid and miltirone analogue. Bioassays showed that (-)-1 and (+)-1 inhibited nitric oxide production in lipopolysaccharide-induced RAW264.7 macrophage cells with an IC50 value of 6.48 ± 1.25 and 15.76 ± 5.55 μM, respectively. The structure-based virtual screening based on the pharmacophores in ePharmaLib, as well as the molecular docking and molecular dynamics simulations study, implied that (-)-1 and (+)-1 may potentially bind to retinoic acid receptor-related orphan receptor C to exert anti-inflammatory activities.

CYP72D19 from Tripterygium wilfordii catalyzes C-2 hydroxylation of abietane-type diterpenoids

KEY MESSAGE

CYP72D19, the first functional gene of the CYP72D subfamily, catalyzes the C-2 hydroxylation of abietane-type diterpenoids. The abietane-type diterpenoids, e.g., triptolide, tripdiolide, and 2-epitripdiolide, are the main natural products for the anti-tumor, anti-inflammatory, and immunosuppressive activities of Tripterygium wilfordii, while their biosynthetic pathways are not resolved. Here, we cloned and characterized the CYP72D19-catalyzed C-2 hydroxylation of dehydroabietic acid, a compound that has been proven to be a biosynthetic intermediate in triptolide biosynthesis. Through molecular docking and site-directed mutagenesis, L386, L387, and I493 near the active pocket were found to have an important effect on the enzyme activity, which also indicates that steric hindrance of residues plays an important role in function. In addition, CYP72D19 also catalyzed a variety of abietane-type diterpenoids with benzene ring, presumably because the benzene ring of the substrate molecule stabilized the C-ring, allowing the protein and the substrate to form a relatively stable spatial structure. This is the first demonstration of CYP72D subfamily gene function. Our research provides important genetic elements for the structural modification of active ingredients and the heterologous production of other 2-hydroxyl abietane-type natural products.

Identification of novel flavin-dependent monooxygenase from Strobilanthes Cusia reveals molecular basis of indoles' biosynthetic logic

BACKGROUND

Strobilanthes cusia (Nees) Kuntze is a traditional medical plant distributed widely in south China. The indole compounds that originated from the plant are responsible for its pharmacological activities. However, the reason why indole ingredients are accumulated in this herb and how it is biosynthesized has remained largely unknown.

RESULTS

In this study, metabolic and transcriptional profiling measurement experiments of different S. cusia organs were carried out to understand the underlying molecular basis of indoles' biosynthetic logic. A metabolic investigation demonstrated that the indoles are primarily accumulated mainly in aerial parts, particularly in leaves. RNA-seq was employed to reveal the organ specific accumulation of indoles in different S. cusia organs. Meanwhile, a flavin-dependent monooxygenase gene (ScFMO1) was found in S. cusia, and it has capacity to produce indoxyl from indole by the fermentation assay. Finally, we assessed the outcomes of transient expression experiment in tobacco and confirmed that ScFMO1 localizes in cytoplasm.

CONCLUSIONS

Our results suggest that ScFMO1 plays a key role in biosynthesis of indoles (Indigo, indirubin, indican, etc.), it will be useful for illuminating the molecular basis of the medicinal indoles' biosynthesis and developing strategies for improving their yields.

Phenotyping of Salvia miltiorrhiza Roots Reveals Associations between Root Traits and Bioactive Components

Plant phenomics aims to perform high-throughput, rapid, and accurate measurement of plant traits, facilitating the identification of desirable traits and optimal genotypes for crop breeding. Salvia miltiorrhiza (Danshen) roots possess remarkable therapeutic effect on cardiovascular diseases, with huge market demands. Although great advances have been made in metabolic studies of the bioactive metabolites, investigation for S. miltiorrhiza roots on other physiological aspects is poor. Here, we developed a framework that utilizes image feature extraction software for in-depth phenotyping of S. miltiorrhiza roots. By employing multiple software programs, S. miltiorrhiza roots were described from 3 aspects: agronomic traits, anatomy traits, and root system architecture. Through K-means clustering based on the diameter ranges of each root branch, all roots were categorized into 3 groups, with primary root-associated key traits. As a proof of concept, we examined the phenotypic components in a series of randomly collected S. miltiorrhiza roots, demonstrating that the total surface of root was the best parameter for the biomass prediction with high linear regression correlation (R2 = 0.8312), which was sufficient for subsequently estimating the production of bioactive metabolites without content determination. This study provides an important approach for further grading of medicinal materials and breeding practices.

Comparing genomes of Fructus Amomi-producing species reveals genetic basis of volatile terpenoid divergence

Wurfbainia longiligularis and Wurfbainia villosa are both rich in volatile terpenoids and are 2 primary plant sources of Fructus Amomi used for curing gastrointestinal diseases. Metabolomic profiling has demonstrated that bornyl diphosphate (BPP)-related terpenoids are more abundant in the W. villosa seeds and have a wider tissue distribution in W. longiligularis. To explore the genetic mechanisms underlying the volatile terpenoid divergence, a high-quality chromosome-level genome of W. longiligularis (2.29 Gb, contig N50 of 80.39 Mb) was assembled. Functional characterization of 17 terpene synthases (WlTPSs) revealed that WlBPPS, along with WlTPS 24/26/28 with bornyl diphosphate synthase (BPPS) activity, contributes to the wider tissue distribution of BPP-related terpenoids in W. longiligularis compared to W. villosa. Furthermore, transgenic Nicotiana tabacum showed that the GCN4-motif element positively regulates seed expression of WvBPPS and thus promotes the enrichment of BPP-related terpenoids in W. villosa seeds. Systematic identification and analysis of candidate TPS in 29 monocot plants from 16 families indicated that substantial expansion of TPS-a and TPS-b subfamily genes in Zingiberaceae may have driven increased diversity and production of volatile terpenoids. Evolutionary analysis and functional identification of BPPS genes showed that BPP-related terpenoids may be distributed only in the Zingiberaceae of monocot plants. This research provides valuable genomic resources for breeding and improving Fructus Amomi with medicinal and edible value and sheds light on the evolution of terpenoid biosynthesis in Zingiberaceae.

Prokaryotic expression and characterization of artificial self-sufficient CYP120A monooxygenases

Cytochrome P450 monooxygenases CYP120As are the unique non-membrane P450s, which are extensively involved in retinoid biodegradation. As the O-functionalized 1,3,3-trimethylcyclohex-1-ene moiety exists in many bioactive compounds which could only be catalyzed by Class II P450s, exploration of the catalytic repertoire of CYP120As is therefore highly attractive. However, up to date, only one bacteriogenic candidate (CYP120A1) was demonstrated for the hydroxylation of C16 and C17 of retinoic acid, by utilizing the integral membrane protein cytochrome P450 reductase redox partner for the electron transfer. Herein, we provided an efficient prokaryotic functional expression system of CYP120As in E. coli by expression of the CYP120A1 coupled with several reductase partners. Fusion redox partners to the C-terminal of the heme-domain are also working on other CYP120A members. Among them, the fusion protein of CYP120A29 and FAD/FMN reductase from Bacillus megaterium P450BM3 (CYP101A2) showed the highest expression level. Based on the available translational fusion systems, the regioselectivity and the substrate scope of the CYP120As have also been explored. This work represents a good starting point for further expanding the catalytic potential of CYP120 family. KEY POINTS: • Characterization of CYP120As in E. coli is firstly achieved by constructing fusion proteins. • The feasibility of three P450 reductase domains to CYP120As was evaluated. • Hydroxylated products of retinoic acid by six CYP120As were sorted and analyzed.

Functional divergence of CYP76AKs shapes the chemodiversity of abietane-type diterpenoids in genus Salvia

The genus Salvia L. (Lamiaceae) comprises myriad distinct medicinal herbs, with terpenoids as one of their major active chemical groups. Abietane-type diterpenoids (ATDs), such as tanshinones and carnosic acids, are specific to Salvia and exhibit taxonomic chemical diversity among lineages. To elucidate how ATD chemical diversity evolved, we carried out large-scale metabolic and phylogenetic analyses of 71 Salvia species, combined with enzyme function, ancestral sequence and chemical trait reconstruction, and comparative genomics experiments. This integrated approach showed that the lineage-wide ATD diversities in Salvia were induced by differences in the oxidation of the terpenoid skeleton at C-20, which was caused by the functional divergence of the cytochrome P450 subfamily CYP76AK. These findings present a unique pattern of chemical diversity in plants that was shaped by the loss of enzyme activity and associated catalytic pathways.

A bacterial-like Pictet-Spenglerase drives the evolution of fungi to produce β-carboline glycosides together with separate genes

Diverse β-carboline (βC) alkaloids are produced by microbes, plants, and animals with myriad bioactivities and drug potentials. However, the biosynthetic mechanism of βCs remains largely elusive, especially regarding the hydroxyl and glucosyl modifications of βCs. Here, we report the presence of the bacterial-like Pictet-Spenglerase gene Fcs1 in the entomopathogenic Beauveria fungi that can catalyze the biosynthesis of the βC skeleton. The overexpression of Fcs1 in Beauveria bassiana led to the identification of six βC methyl glycosides, termed bassicarbosides (BCSs) A-F. We verified that the cytochrome P450 (CYP) genes adjacent to Fcs1 cannot oxidize βCs. Alternatively, the separated CYP684B2 family gene Fcs2 was identified to catalyze βC hydroxylation together with its cofactor gene Fcs3. The functional homologue of Fcs2 is only present in the Fcs1-containing fungi and highly similar to the Fcs1-connected yet nonfunctional CYP. Both evolved quicker than those from fungi without Fcs1 homologues. Finally, the paired methyl/glucosyl transferase genes were verified to mediate the production of BCSs from hydroxy-βCs. All these functionally verified genes are located on different chromosomes of Beauveria, which is in contrast to the typical content-clustered feature of fungal biosynthetic gene clusters (BGCs). We also found that the production of BCSs selectively contributed to fungal infection of different insect species. Our findings shed light on the biosynthetic mechanism of βC glycosides, including the identification of a βC hydroxylase. The results of this study also propose an evolving process of fungal BGC formation following the horizontal transfer of a bacterial gene to fungi.

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.

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

MAIN CONCLUSION

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

Functional Study and Efficient Catalytic Element Mining of CYP76AHs in Salvia Plants

Salvia is a large genus with hundreds of species used in traditional Chinese medicine. Tanshinones are a highly representative class of exclusive compounds found in the Salvia genus that exhibit significant biological activity. Tanshinone components have been identified in 16 Salvia species. The CYP76AH subfamily (P450) is crucial for the synthesis of tanshinone due to its catalytic generation of polyhydroxy structures. In this study, a total of 420 CYP76AH genes were obtained, and phylogenetic analysis showed their clear clustering relationships. Fifteen CYP76AH genes from 10 Salvia species were cloned and studied from the perspectives of evolution and catalytic efficiency. Three CYP76AHs with significantly improved catalytic efficiency compared to SmCYP76AH3 were identified, providing efficient catalytic elements for the synthetic biological production of tanshinones. A structure-function relationship study revealed several conserved residues that might be related to the function of CYP76AHs and provided a new mutation direction for the study of the directed evolution of plant P450.

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.

A developmental gradient reveals biosynthetic pathways to eukaryotic toxins in monocot geophytes

Numerous eukaryotic toxins that accumulate in geophytic plants are valuable in the clinic, yet their biosynthetic pathways have remained elusive. A lead example is the >150 Amaryllidaceae alkaloids (AmAs) including galantamine, an FDA-approved treatment for Alzheimer's disease. We show that while AmAs accumulate to high levels in many tissues in daffodils, biosynthesis is localized to nascent, growing tissue at the base of leaves. A similar trend is found for the production of steroidal alkaloids (e.g. cyclopamine) in corn lily. This model of active biosynthesis enabled elucidation of a complete set of biosynthetic genes for the production of AmAs. Taken together, our work sheds light on the developmental and enzymatic logic of diverse alkaloid biosynthesis in daffodil. More broadly, it suggests a paradigm for biosynthesis regulation in monocot geophytes where plants are protected from herbivory through active charging of newly formed cells with eukaryotic toxins that persist as aboveground tissue develops.

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.

Identification and functional characterization of three cytochrome P450 genes for the abietane diterpenoid biosynthesis in Isodon lophanthoides

We identify two ferruginol synthases and a 11-hydroxyferruginol synthase from a traditional Chinese medicinal herb Isodon lophanthoides and propose their involvement in two independent abietane diterpenoids biosynthetic pathways. Isodon lophanthoides is a traditional Chinese medicinal herb rich in highly oxidized abietane-type diterpenoids. These compounds exhibit a wide range of pharmaceutical activities, yet the biosynthesis is barely known. Here, we describe the screening and functional characterization of P450s that oxidize the abietane skeleton abietatriene. We mainly focused on CYP76 family and identified 12 CYP76AHs by mining the RNA-seq data of I. lophanthoides. Among the 12 CYP76AHs, 6 exhibited similar transcriptional expression features as upstream diterpene synthases, including root or leaf-preferential expression pattern and highly MeJA inducibility. These six P450s were considered as first-tier candidates and functionally characterized in yeast and plant cells. In yeast assays showed that both CYP76AH42 and CYP76AH43 were ferruginol synthases hydroxylating the C12 position of abietatriene, whereas CYP76AH46 was characterized as a 11-hydroxyferruginol synthase which catalyzes two successive oxidations at C12 and C11 of abietatriene. Heterologous expression of three CYP76AHs in Nicotiana benthamiana resulted in the formation of ferruginol. qPCR analysis showed CYP76AH42 and CYP76AH43 were mainly expressed in the root, which was consistent with the distribution of ferruginol in the root periderms. CYP76AH46 was primarily expressed in the leaves where barely ferruginol or 11-hydroxyferruginol was detected. In addition to distinct organ-specific expression pattern, three CYP76AHs exhibited different genomic structures (w or w/o introns), low protein sequence identities (51-63%) and were placed in separate subclades in the phylogenetic tree. These results suggest that the identified CYP76AHs may be involved in at least two independent abietane biosynthetic pathways in the aerial and underground parts of I. lophanthoides.

A chromosome-level genome assembly reveals that a bipartite gene cluster formed via an inverted duplication controls monoterpenoid biosynthesis in Schizonepeta tenuifolia

Biosynthetic gene clusters (BGCs) are regions of a genome where genes involved in a biosynthetic pathway are in proximity. The origin and evolution of plant BGCs as well as their role in specialized metabolism remain largely unclear. In this study, we have assembled a chromosome-scale genome of Japanese catnip (Schizonepeta tenuifolia) and discovered a BGC that contains multiple copies of genes involved in four adjacent steps in the biosynthesis of p-menthane monoterpenoids. This BGC has an unprecedented bipartite structure, with mirrored biosynthetic regions separated by 260 kilobases. This bipartite BGC includes identical copies of a gene encoding an old yellow enzyme, a type of flavin-dependent reductase. In vitro assays and virus-induced gene silencing revealed that this gene encodes the missing isopiperitenone reductase. This enzyme evolved from a completely different enzyme family to isopiperitenone reductase from closely related Mentha spp., indicating convergent evolution of this pathway step. Phylogenomic analysis revealed that this bipartite BGC has emerged uniquely in the S. tenuifolia lineage and through insertion of pathway genes into a region rich in monoterpene synthases. The cluster gained its bipartite structure via an inverted duplication. The discovered bipartite BGC for p-menthane biosynthesis in S. tenuifolia has similarities to the recently described duplicated p-menthane biosynthesis gene pairs in the Mentha longifolia genome, providing an example of the convergent evolution of gene order. This work expands our understanding of plant BGCs with respect to both form and evolution, and highlights the power of BGCs for gene discovery in plant biosynthetic pathways.

A chromosome-level genome assembly reveals that tandem-duplicated CYP706V oxidase genes control oridonin biosynthesis in the shoot apex of Isodon rubescens

The ent-kaurenoids (e.g., oridonin and enmein) from the Isodon genus (Lamiaceae) are one class of diterpenoids with rich structural diversity and intriguing pharmaceutical activity. In contrast to the well-established gibberellin pathway, oxidative modifications diversifying the ent-kaurene skeleton in Isodon have remained undetermined for half a century. Here we report a chromosome-level genome assembly of I. rubescens, a well-recognized oridonin producer long favored by Asian people as a traditional herb with antitumor effects. The shoot apex was confirmed to be the actual region actively producing ent-kaurene diterpenoids. Through comparative genomics and phylogenetic analyses, we discovered a cluster of tandem-duplicated CYP706V oxygenase-encoding genes located on an ancient genomic block widely distributed in eudicots, whereas almost exclusively emerged in Isodon plants. In the shoot apex, IrCYP706V2 and IrCYP706V7 oxidized the ent-kaurene core in the initial stage of oridonin biosynthesis. Loss of CYP706Vs in other Lamiaceae plants offered an explanation for the specific kaurenoid production in Isodon plants. Moreover, we found that the Isodon genomes encode multiple diterpenoid synthases that are potentially involved in generating diterpenoid diversity. These findings provided new insights into the evolution of the lineage-specific diterpenoid pathway and laid a foundation for improving production of bioactive ent-kaurene-type diterpenoids by molecular breeding and synthetic biology approaches.

Plant (di)terpenoid evolution: from pigments to hormones and beyond

Covering: up to 2014-2022.Diterpenoid biosynthesis in plants builds on the necessary production of (E,E,E)-geranylgeranyl diphosphate (GGPP) for photosynthetic pigment production, with diterpenoid biosynthesis arising very early in land plant evolution, enabling stockpiling of the extensive arsenal of (di)terpenoid natural products currently observed in this kingdom. This review will build upon that previously published in the Annual Review of Plant Biology, with a stronger focus on enzyme structure-function relationships, as well as additional insights into the evolution of (di)terpenoid metabolism since generated.

Promiscuity, a Driver of Plant Cytochrome P450 Evolution?

Plant cytochrome P450 monooxygenases were long considered to be highly substrate-specific, regioselective and stereoselective enzymes, in this respect differing from their animal counterparts. The functional data that have recently accumulated clearly counter this initial dogma. Highly promiscuous P450 enzymes have now been reported, mainly in terpenoid pathways with functions in plant adaptation, but also some very versatile xenobiotic/herbicide metabolizers. An overlap and predictable interference between endogenous and herbicide metabolism are starting to emerge. Both substrate preference and permissiveness vary between plant P450 families, with high promiscuity seemingly favoring retention of gene duplicates and evolutionary blooms. Yet significant promiscuity can also be observed in the families under high negative selection and with essential functions, usually enhanced after gene duplication. The strategies so far implemented, to systematically explore P450 catalytic capacity, are described and discussed.

Tandemly duplicated CYP82Ds catalyze 14-hydroxylation in triptolide biosynthesis and precursor production in Saccharomyces cerevisiae

Triptolide is a valuable multipotent antitumor diterpenoid in Tripterygium wilfordii, and its C-14 hydroxyl group is often selected for modification to enhance both the bioavailability and antitumor efficacy. However, the mechanism for 14-hydroxylation formation remains unknown. Here, we discover 133 kb of tandem duplicated CYP82Ds encoding 11 genes on chromosome 12 and characterize CYP82D274 and CYP82D263 as 14-hydroxylases that catalyze the metabolic grid in triptolide biosynthesis. The two CYP82Ds catalyze the aromatization of miltiradiene, which has been repeatedly reported to be a spontaneous process. In vivo assays and evaluations of the kinetic parameters of CYP82Ds indicate the most significant affinity to dehydroabietic acid among multiple intermediates. The precursor 14-hydroxy-dehydroabietic acid is successfully produced by engineered Saccharomyces cerevisiae. Our study provides genetic elements for further elucidation of the downstream biosynthetic pathways and heterologous production of triptolide and of the currently intractable biosynthesis of other 14-hydroxyl labdane-type secondary metabolites.

Structure-function analysis of CYP719As involved in methylenedioxy bridge-formation in the biosynthesis of benzylisoquinoline alkaloids and its de novo production

Benzylisoquinoline alkaloids (BIAs) are a type of secondary metabolite with clinical application value. (S)-stylopine is a special BIA which contains methylenedioxy bridge structures. CYP719As could catalyze the methylenedioxy bridge-formation on the A or D rings of protoberberine alkaloids, while displaying significant substrate regiospecificity. To explore the substrate preference of CYP719As, we cloned and identified five CyCYP719A candidates from Corydalis yanhusuo. Two CyCYP719As (CyCYP719A39 and CyCYP719A42) with high catalytic efficiency for the methylenedioxy bridge-formation on the D or A rings were characterized, respectively. The residues (Leu 294 for CyCYP719A42 and Asp 289 for CyCYP719A39) were identified as the key to controlling the regioselectivity of CYP719As affecting the methylenedioxy bridge-formation on the A or D rings by homology modeling and mutation analysis. Furthermore, for de novo production of BIAs, CyCYP719A39, CyCYP719A42, and their mutants were introduced into the (S)-scoulerine-producing yeast to produce 32 mg/L (S)-stylopine. These results lay a foundation for understanding the structure-function relationship of CYP719A-mediated methylenedioxy bridge-formation and provide yeast strains for the BIAs production by synthetic biology.

Leveraging yeast to characterize plant biosynthetic gene clusters

Plant biosynthetic gene clusters (BGCs) contain multiple physically clustered non-homologous genes that encode enzymes catalyzing diverse reactions in one plant natural product biosynthetic pathway. A growing number of plant BGCs have emerged as an underlying resource for understanding plant specialized metabolism and evolution, but the characterization remains challenging. Recent studies have demonstrated that baker's yeast can serve as a versatile platform for the characterization of plant BGCs, from single-gene characterization to multiple genes and hitherto unknown putative BGC validation and elucidation. In this review, we will summarize the strategies and examples of the applications of yeast in plant BGC characterization and share our perspective on the development of a systematic pipeline to fully leverage yeast to advance the understanding of plant BGCs and plant natural product biomanufacturing.

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.

Metabolomics-centered mining of plant metabolic diversity and function: Past decade and future perspectives

Plants are natural experts in organic synthesis, being able to generate large numbers of specific metabolites with widely varying structures that help them adapt to variable survival challenges. Metabolomics is a research discipline that integrates the capabilities of several types of research including analytical chemistry, statistics, and biochemistry. Its ongoing development provides strategies for gaining a systematic understanding of quantitative changes in the levels of metabolites. Metabolomics is usually performed by targeting either a specific cell, a specific tissue, or the entire organism. Considerable advances in science and technology over the last three decades have propelled us into the era of multi-omics, in which metabolomics, despite at an earlier developmental stage than genomics, transcriptomics, and proteomics, offers the distinct advantage of studying the cellular entities that have the greatest influence on end phenotype. Here, we summarize the state of the art of metabolite detection and identification, and illustrate these techniques with four case study applications: (i) comparing metabolite composition within and between species, (ii) assessing spatio-temporal metabolic changes during plant development, (iii) mining characteristic metabolites of plants in different ecological environments and upon exposure to various stresses, and (iv) assessing the performance of metabolomics as a means of functional gene identification , metabolic pathway elucidation, and metabolomics-assisted breeding through analyzing plant populations with diverse genetic variations. In addition, we highlight the prominent contributions of joint analyses of plant metabolomics and other omics datasets, including those from genomics, transcriptomics, proteomics, epigenomics, phenomics, microbiomes, and ion-omics studies. Finally, we discuss future directions and challenges exploiting metabolomics-centered approaches in understanding plant metabolic diversity.

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.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

De Novo Hybrid Assembly of the Salvia miltiorrhiza Mitochondrial Genome Provides the First Evidence of the Multi-Chromosomal Mitochondrial DNA Structure of Salvia Species

Salvia miltiorrhiza has been an economically important medicinal plant. Previously, an S. miltiorrhiza mitochondrial genome (mitogenome) assembled from Illumina short reads, appearing to be a single circular molecule, has been published. Based on the recent reports on the plant mitogenome structure, we suspected that this conformation does not accurately represent the complexity of the S. miltiorrhiza mitogenome. In the current study, we assembled the mitogenome of S. miltiorrhiza using the PacBio and Illumina sequencing technologies. The primary structure of the mitogenome contained two mitochondrial chromosomes (MC1 and MC2), which corresponded to two major conformations, namely, Mac1 and Mac2, respectively. Using two approaches, including (1) long reads mapping and (2) polymerase chain reaction amplification followed by Sanger sequencing, we observed nine repeats that can mediate recombination. We predicted 55 genes, including 33 mitochondrial protein-coding genes (PCGs), 3 rRNA genes, and 19 tRNA genes. Repeat analysis identified 112 microsatellite repeats and 3 long-tandem repeats. Phylogenetic analysis using the 26 shared PCGs resulted in a tree that was congruent with the phylogeny of Lamiales species in the APG IV system. The analysis of mitochondrial plastid DNA (MTPT) identified 16 MTPTs in the mitogenome. Moreover, the analysis of nucleotide substitution rates in Lamiales showed that the genes atp4, ccmB, ccmFc, and mttB might have been positively selected. The results lay the foundation for future studies on the evolution of the Salvia mitogenome and the molecular breeding of 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.

Integration of genome and transcriptome reveal molecular regulation mechanism of early flowering trait in Prunus genus (Prunus mume and Prunus persica)

Flowering time is crucial for the survival and reproduction. Prunus genus belongs to the Rosaceae family and includes several hundred species of flowering trees and shrubs with important ornamental and economic values. However, the molecular mechanism underlying early flowering in Prunus genus is unclear. Here, we utilized the genome and transcriptome of P. mume and P. persica to explore the transcriptional regulation mechanism of early flowering. Comparative genomics found that genes accounting for 92.4% of the total P. mume genome and 91.2% of the total P. persica genome belonged to orthogroups. A total of 19,169 orthogroups were found between P. mume and P. persica, including 20,431 corresponding orthologues and 20,080 collinearity gene pairs. A total of 305 differentially expressed genes (DEGs) associated with early flowering were found, among which FT, TLI65, and NAP57 were identified as hub genes in the early flowering regulation pathway. Moreover, we identified twenty-five transcription factors (TFs) from nine protein families, including MADS-box, AP2/ERF, and MYB. Our results provide insights into the underlying molecular model of flowering time regulation in Prunus genus and highlight the utility of multi-omics in deciphering the properties of the inter-genus plants.