Analysis of the Genome Sequence of the Medicinal Plant Salvia miltiorrhiza

Genome-Wide Characterization and Expression Analysis of GATA Transcription Factors in Response to Methyl Jasmonate in Salvia miltiorrhiza

Salvia miltiorrhiza is an important medicinal plant, which is mainly used for treatment of cardiovascular and cerebrovascular diseases. GATA transcription factors are evolutionarily conser-ved proteins that play essential roles in biological process of plants. In this study, we systematically characterized the GATA transcription factors in S. miltiorrhiza. A total 28 SmGATA genes were identified and divided into four subfamilies based on phylogenetic analysis and domain. SmGATA genes being clustered into a subfamily have similar conserved motifs and exon-intron patterns, and unevenly distribute on eight chromosomes of S. miltiorrhiza. Tissue-specific expression analysis based on transcriptome datasets showed that the majority of SmGATA genes were preferentially expressed in roots. Under methyl jasmonate (MeJA) treatment, the quantitative real-time PCR (qRT-PCR) analysis indicated that several SmGATA genes in roots showed distinct upregulation post-MeJA treatment, especially SmGATA08, which was highly responsive to MeJA, and might be involved in the jasmonate signal, thereby affecting root growth, development, tolerance to various stresses, or secondary metabolites biosynthesis. The study found that several SmGATAs, like SmGATA08, are highly responsive to MeJA, indicating that these SmGATAs might be vital in the biosynthesis of tanshinones and phenolic acids by regulating the response to MeJA in S. miltiorrhiza. Our results laid the foundation for understanding their biological roles and quality improvement in S. miltiorrhiza.

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

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

CRISPR-Cas gene editing technology and its application prospect in medicinal plants

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas gene editing technology has opened a new era of genome interrogation and genome engineering because of its ease operation and high efficiency. An increasing number of plant species have been subjected to site-directed gene editing through this technology. However, the application of CRISPR-Cas technology to medicinal plants is still in the early stages. Here, we review the research history, structural characteristics, working mechanism and the latest derivatives of CRISPR-Cas technology, and discussed their application in medicinal plants for the first time. Furthermore, we creatively put forward the development direction of CRISPR technology applied to medicinal plant gene editing. The aim is to provide a reference for the application of this technology to genome functional studies, synthetic biology, genetic improvement, and germplasm innovation of medicinal plants. CRISPR-Cas is expected to revolutionize medicinal plant biotechnology in the near future.

Chromosome-Scale Assembly of the Dendrobium nobile Genome Provides Insights Into the Molecular Mechanism of the Biosynthesis of the Medicinal Active Ingredient of Dendrobium

Orchids constitute approximately 10% of flowering plant species. However, only about 10 orchid genomes have been published. Metabolites are the main way through which orchids respond to their environment. Dendrobium nobile, belonging to Dendrobium, the second largest genus in Orchidaceae, has high ornamental, medicinal, and ecological value. D. nobile is the source of many popular horticultural varieties. Among the Dendrobium species, D. nobile has the highest amount of dendrobine, which is regarded as one of the criteria for evaluating medicinal quality. Due to lack of data and analysis at the genomic level, the biosynthesis pathways of dendrobine and other related medicinal ingredients in D. nobile are unknown. In this paper, we report a chromosome-scale reference genome of D. nobile to facilitate the investigation of its genomic characteristics for comparison with other Dendrobium species. The assembled genome size of D. nobile was 1.19 Gb. Of the sequences, 99.45% were anchored to 19 chromosomes. Furthermore, we identified differences in gene number and gene expression patterns compared with two other Dendrobium species by integrating whole-genome sequencing and transcriptomic analysis [e.g., genes in the polysaccharide biosynthesis pathway and upstream of the alkaloid (dendrobine) biosynthesis pathway]. Differences in the TPS and CYP450 gene families were also found among orchid species. All the above differences might contribute to the species-specific medicinal ingredient biosynthesis pathways. The metabolic pathway-related analysis will provide further insight into orchid responses to the environment. Additionally, the reference genome will provide important insights for further molecular elucidation of the medicinal active ingredients of Dendrobium and enhance the understanding of orchid evolution.

Comparative genomics reveal the convergent evolution of CYP82D and CYP706X members related to flavone biosynthesis in Lamiaceae and Asteraceae

Distant species producing the same secondary metabolites is an interesting and common phenomenon in nature. A classic example of this is scutellarein whose derivatives have been used clinically for more than 30 years. Scutellarein occurs in significant amounts in species of two different orders, Scutellaria baicalensis and Erigeron breviscapus, which diverged more than 100 million years ago. Here, according to the genome-wide selection and functional identification of 39 CYP450 genes from various angiosperms, we confirmed that only seven Scutellaria-specific CYP82D genes and one Erigeron CYP706X gene could perform the catalytic activity of flavone 6-hydroxylase (F6H), suggesting that the convergent evolution of scutellarein production in these two distant species was caused by two independently evolved CYP450 families. We also identified seven Scutellaria-specific CYP82D genes encoding flavone 8-hydroxylase (F8H). The evolutionary patterns of CYP82 and CYP706 families via kingdom-wide comparative genomics highlighted the evolutionary diversity of CYP82D and the specificity of CYP706X in angiosperms. Multi-collinearity and phylogenetic analysis of CYP82D in Scutellaria confirmed that the function of F6H evolved from F8H. Furthermore, the SbaiCYP82D1^A319D , EbreCYP706X^R130A , EbreCYP706X^F312D and EbreCYP706X^A318D mutants can significantly decrease the catalytic activity of F6H, revealing the contribution of crucial F6H amino acids to the scutellarein biosynthesis of distant species. This study provides important insights into the multi-origin evolution of the same secondary metabolite biosynthesis in the plant kingdom.

Genome-Wide Analysis of U-box E3 Ubiquitin Ligase Family in Response to ABA Treatment in Salvia miltiorrhiza

Plant U-box (PUB) proteins are ubiquitin ligases (E3) involved in multiple biological processes and in response to plant stress. However, the various aspects of the genome and the differences in functions between the U-box E3 (UBE3) ubiquitin ligases remain quite obscure in Salvia miltiorrhiza. The 60 UBE3 genes in the S. miltiorrhiza genome were recognized in the present study. The phylogenetic analysis, gene structure, motifs, promoters, and physical and chemical properties of the genes were also examined. Based on the phylogenetic relationship, the 60 UBE3 genes were categorized under six different groups. The U-box domain was highly conserved across the family of UBE3 genes. Analysis of the cis-acting element revealed that the UBE3 genes might play an important role in a variety of biological processes, including a reaction to the abscisic acid (ABA) treatment. To investigate this hypothesis, an ABA treatment was developed for the hairy roots of S. miltiorrhiza. Thirteen out of the UBE3 genes significantly increased after the ABA treatment. The co-expression network revealed that nine UBE3 genes might be associated with phenolic acids or tanshinone biosynthesis. The findings of the present study brought fresh and new understanding to the participation of the UBE3 gene family in plants, specifically in their biological responses mediated by the ABA. In S. miltiorrhiza, this gene family may be crucial during the ABA treatment. Significantly, the results of this study contribute novel information to the understanding of the ubiquitin ligase gene and its role in plant growth.

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

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

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

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

Genome survey sequencing and characterization of simple sequence repeat (SSR) markers in Platostoma palustre (Blume) A.J.Paton (Chinese mesona)

Platostoma palustre (Blume) A.J.Paton is an annual herbaceous persistent plant of the Labiatae family. However, there is a lack of genomic data for this plant, which severely restricts its genetic improvement. In this study, we performed genome survey sequencing of P. palustre and developed simple sequence repeat (SSR) markers based on the resulting sequence. K-mer analysis revealed that the assembled genome size was approximately 1.21 Gb. A total of 15,498 SSR motifs were identified and characterized in this study; among them, dinucleotide, and hexanucleotide repeats had the highest and lowest, respectively. Among the dinucleotide repeat motifs, AT/TA repeat motifs were the most abundant, and GC/CG repeat motifs were rather rare, accounting for 44.28% and 0.63%, respectively. Genetic similarity coefficient analysis by the UPMGA methods clustered 12 clones, of P. palustre and related species into two subgroups. These results provide helpful information for further research on P. palustre resources and variety improvements.

P2Y12 antagonists: Approved drugs, potential naturally isolated and synthesised compounds, and related in-silico studies

P2Y₁₂ is a platelet surface protein which is responsible for the amplification of P2Y₁ response. It plays a crucial role in platelet aggregation and thrombus formation through an ADP-induced platelet activation mechanism. Despite that P2Y₁₂ platelets' receptor is an excellent target for developing antiplatelet agents, only five approved medications are currently in clinical use which are classified into thienopyridines and nucleoside-nucleotide derivatives. In the past years, many attempts for developing new candidates as P2Y₁₂ inhibitors have been made. This review highlights the importance and the role of P2Y₁₂ receptor as part of the coagulation cascade, its reported congenital defects, and the type of assays which are used to verify and measure its activity. Furthermore, an overview is given of the clinically approved medications, the potential naturally isolated inhibitors, and the synthesised candidates which were tested either in-vitro, in-vivo and/or clinically. Finally, we outline the in-silico attempts which were carried out using virtual screening, molecular docking and dynamics simulations in efforts of designing novel P2Y₁₂ antagonists. Various phytochemical classes might be considered as a corner stone for the discovery of novel P2Y₁₂ inhibitors, whereas a wide range of ring systems can be deliberated as leading scaffolds in that area synthetically and theoretically.

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

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

Genome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza

Circular RNAs (circRNAs) play various roles in cellular functions. However, no studies have been reported on the potential involvement of circRNAs in the biosynthesis of secondary metabolites in plants. Here, we performed a genome-wide discovery of circRNAs from root, stem and leaf samples of Salvia miltiorrhiza using RNA-Seq. We predicted a total of 2,476 circRNAs with at least two junction reads using circRNA_finder and CIRI, of which 2,096, 151 and 229 were exonic, intronic and intergenic circRNAs, respectively. Sequence similarity analysis showed that 294 out of 2,476 circRNAs were conserved amongst multiple plants. Of the 55 predicted circRNAs, 31 (56%) were validated successfully by PCR and Sanger sequencing using convergent and divergent primer pairs. Alternative circularisation analysis showed that most parental genes produced two circRNAs. Functional enrichment analyses of the parental genes showed that the primary metabolism pathways were significantly enriched, particularly the carbon metabolism. Differential expression analysis showed that the expression profiles of circRNAs were tissue-specific. Co-expression analysis showed 275 circRNAs, and their parental genes had significantly positive correlations. However, 14 had significantly negative correlations. Weighted gene co-expression network analysis showed that nine circRNAs were co-expressed with four modules of protein-coding genes. Next, we found 416 exonic circRNAs with miRNA-binding sites, suggesting possible interactions between circRNAs and miRNAs. Lastly, we found six validated circRNAs, namely, SMscf2473-46693-46978, SMscf3091-29256-29724, SMscf16-111773-112193, SMscf432-13232-13866, SMscf7007-10563-10888 and SMscf1730-1749-2013, which were originated from the genes involved in the biosynthesis of secondary metabolites. Their parental genes were acetyl-CoA C-acetyltransferase 1 (SmAACT1), 1-deoxy-d-xylulose-5-phosphate synthase 2 (SmDXS2), 4-hydroxy-3-methylbut-2-enyl diphosphate reductase 1 (SmHDR1), kaurene synthase-like 2 (SmKSL2), DWF4 and CYP88A3, respectively. In particular, the correlation coefficient of SMscf2473-46693-46978 and SmDXS2 gene was 0.86 (p = 0.003), indicating a potential interaction between this pair of circRNA and its parent gene. Our results provided the first comprehensive catalogue of circRNAs in S. miltiorrhiza and identified one circRNA that might play important roles in the biosynthesis of secondary metabolites.

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.

Sclareol and linalyl acetate are produced by glandular trichomes through the MEP pathway

Sclareol, an antifungal specialized metabolite produced by clary sage, Salvia sclarea, is the starting plant natural molecule used for the hemisynthesis of the perfume ingredient ambroxide. Sclareol is mainly produced in clary sage flower calyces; however, the cellular localization of the sclareol biosynthesis remains unknown. To elucidate the site of sclareol biosynthesis, we analyzed its spatial distribution in the clary sage calyx epidermis using laser desorption/ionization mass spectrometry imaging (LDI-FTICR-MSI) and investigated the expression profile of sclareol biosynthesis genes in isolated glandular trichomes (GTs). We showed that sclareol specifically accumulates in GTs' gland cells in which sclareol biosynthesis genes are strongly expressed. We next isolated a glabrous beardless mutant and demonstrate that more than 90% of the sclareol is produced by the large capitate GTs. Feeding experiments, using 1-¹³C-glucose, and specific enzyme inhibitors further revealed that the methylerythritol-phosphate (MEP) biosynthetic pathway is the main source of isopentenyl diphosphate (IPP) precursor used for the biosynthesis of sclareol. Our findings demonstrate that sclareol is an MEP-derived diterpene produced by large capitate GTs in clary sage emphasing the role of GTs as biofactories dedicated to the production of specialized metabolites.

Metabolomics integrated with transcriptomics reveals the distribution of iridoid and crocin metabolic flux in Gardenia jasminoides Ellis

Gardenia jasminoides Ellis (G. jasminoides) fruits are used as a resource for obtaining natural colorants and in traditional Chinese herbal medicine. However, G. jasminoides presents a relatively long flowering period and different ripening periods, so there are significant differences in the accumulation of metabolites in fruits of different colors. In addition, the complete metabolic pathways of iridoidsand crocins, which are used as medicinal composition of G. jasminoides, are poorly understood at present. In this research, we comprehensively compared the transcriptome and metabolites profiles of the developmental stages and locations of iridoid and crocin biosynthesis. A large number of differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were detected in four groups of samples, and clear variation in the pattern of metabolite abundance and gene expression were observed among different fruit colors and parts. Geniposide and gardenoside mainly accumulated in the sarcocarp of green fruit (GFS) and the sarcocarp of red fruit (FS), respectively. Crocin mainly accumulated in the peel and sarcocarp of red fruits. In the iridoid pathway, we hypothesized that there was a transport mechanism from the sarcocarp to the peel of G. jasminoides because of the inconsistent expression of G8O, 10-HGO and IS associated with differences in fruit ripening. UGTs play an important role in the biosynthesis of the active components of G. jasminoides. Combined transcriptome and metabonomics analysis showed a negative correlation between the biosynthesis of geniposide and crocin. The redirection of the metabolic flux and the regulation of key enzymes may be the main reasons for the changes in the biosynthesis of iridoid and crocin in G. jasminoides fruit. Our study expended valuable information for functional genomic library and provided new insights for metabolic engineering of secondary metabolite in G. Jasminoides.

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

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

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

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

Chromosome-scale assembly and evolution of the tetraploid Salvia splendens (Lamiaceae) genome

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species.

Plant Metabolic Gene Clusters: Evolution, Organization, and Their Applications in Synthetic Biology

Plants are a remarkable source of high-value specialized metabolites having significant physiological and ecological functions. Genes responsible for synthesizing specialized metabolites are often clustered together for a coordinated expression, which is commonly observed in bacteria and filamentous fungi. Similar to prokaryotic gene clustering, plants do have gene clusters encoding enzymes involved in the biosynthesis of specialized metabolites. More than 20 gene clusters involved in the biosynthesis of diverse metabolites have been identified across the plant kingdom. Recent studies demonstrate that gene clusters are evolved through gene duplications and neofunctionalization of primary metabolic pathway genes. Often, these clusters are tightly regulated at nucleosome level. The prevalence of gene clusters related to specialized metabolites offers an attractive possibility of an untapped source of highly useful biomolecules. Accordingly, the identification and functional characterization of novel biosynthetic pathways in plants need to be worked out. In this review, we summarize insights into the evolution of gene clusters and discuss the organization and importance of specific gene clusters in the biosynthesis of specialized metabolites. Regulatory mechanisms which operate in some of the important gene clusters have also been briefly described. Finally, we highlight the importance of gene clusters to develop future metabolic engineering or synthetic biology strategies for the heterologous production of novel metabolites.

The first draft genome of Picrorhiza kurrooa, an endangered medicinal herb from Himalayas

Picrorhiza kurrooa is an endangered medicinal herb which is distributed across the Himalayan region at an altitude between 3000–5000 m above mean sea level. The medicinal properties of P. kurrooa are attributed to monoterpenoid picrosides present in leaf, rhizome and root of the plant. However, no genomic information is currently available for P. kurrooa, which limits our understanding about its molecular systems and associated responses. The present study brings the first assembled draft genome of P. kurrooa by using 227 Gb of raw data generated by Illumina and PacBio RS II sequencing platforms. The assembled genome has a size of n = ~ 1.7 Gb with 12,924 scaffolds. Four pronged assembly quality validations studies, including experimentally reported ESTs mapping and directed sequencing of the assembled contigs, confirmed high reliability of the assembly. About 76% of the genome is covered by complex repeats alone. Annotation revealed 24,798 protein coding and 9789 non-coding genes. Using the assembled genome, a total of 710 miRNAs were discovered, many of which were found responsible for molecular response against temperature changes. The miRNAs and targets were validated experimentally. The availability of draft genome sequence will aid in genetic improvement and conservation of P. kurrooa. Also, this study provided an efficient approach for assembling complex genomes while dealing with repeats when regular assemblers failed to progress due to repeats.

Insights into salvianolic acid B biosynthesis from chromosome-scale assembly of the Salvia bowleyana genome

Salvia bowleyana is a traditional Chinese medicinal plant that is a source of nutritional supplements rich in salvianolic acid B and a potential experimental system for the exploration of salvianolic acid B biosynthesis in the Labiatae. Here, we report a high-quality chromosome-scale genome assembly of S. bowleyana covering 462.44 Mb, with a scaffold N50 value of 57.96 Mb and 44,044 annotated protein-coding genes. Evolutionary analysis revealed an estimated divergence time between S. bowleyana and its close relative S. miltiorrhiza of ~3.94 million years. We also observed evidence of a whole-genome duplication in the S. bowleyana genome. Transcriptome analysis showed that SbPAL1 (PHENYLALANINE AMMONIA-LYASE1) is highly expressed in roots relative to stem and leaves, paralleling the location of salvianolic acid B accumulation. The laccase gene family in S. bowleyana outnumbered their counterparts in both S. miltiorrhiza and Arabidopsis thaliana, suggesting that the gene family has undergone expansion in S. bowleyana. Several laccase genes were also highly expressed in roots, where their encoded proteins may catalyze the oxidative reaction from rosmarinic acid to salvianolic acid B. These findings provide an invaluable genomic resource for understanding salvianolic acid B biosynthesis and its regulation, and will be useful for exploring the evolution of the Labiatae.

Recent Advances in Molecular Marker-Assisted Breeding for Quality Improvement of Traditional Chinese Medicine

BACKGROUND

The quality of Traditional Chinese Medicine (TCM), reflected by its bioactive compounds and associated contents, is directly linked to its clinical efficacy. Therefore, it is of great importance to improve the quality of TCM by increasing the bioactive compound content.

METHODS

Mapping the active component content-associated QTLs in TCM and further markerassisted breeding has enabled us to rapidly and effectively cultivate new varieties with high bioactive compound contents, which has opened the door for genetic breeding studies on medicinal plants.

RESULTS

In this paper, a strategy and technical molecular breeding method for TCM are discussed. The development of four methods and progress in functional marker development, as well as the applications of such markers in TCM, are reviewed.

CONCLUSION

The progress in, challenges of, and future of marker-assisted breeding for quality improvement of TCM are discussed, which provide valuable scientific references for future molecular breeding.

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.

Genome-wide identification of the Tubby-Like Protein (TLPs) family in medicinal model plant Salvia miltiorrhiza

Tubby-Like Proteins (TLPs) are important transcription factors with many functions and are found in both animals and plants. In plants, TLPs are thought to be involved in the abiotic stress response. To reveal the potential function of TLPs in the medicinal model plant Salvia miltiorrhiza, we identified 12 S. miltiorrhiza TLPs (SmTLPs) and conducted a comprehensive analysis. We examined SmTLP gene structure, protein structure, phylogenetics, and expression analysis. Our results show that all SmTLPs, except SmTLP11, have a complete typical Tub domain. Promoter analysis revealed that most SmTLPs are involved in hormone and abiotic stress responses. Expression analysis revealed that the 12 SmTLPs could be divided into three categories: those specifically expressed in roots, those specifically expressed in stems, and those specifically expressed in leaves. Additional studies have shown that SmTLP10 may play an important role in the plant cold resistance, while SmTLP12 may be involved in the S. miltiorrhiza ABA metabolic pathway. Our study represents the first comprehensive investigation of TLPs in S. miltiorrhiza. These data may provide useful clues for future studies and may support the hypotheses regarding the role of TLPs in plant abiotic stress process. All in all, we may provide a reference for improving S. miltiorrhiza quality using genetic engineering technology.

Crosstalk of Multi-Omics Platforms with Plants of Therapeutic Importance

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

Conservation Genomics of Wild Red Sage (Salvia miltiorrhiza) and Its Endangered Relatives in China: Population Structure and Interspecific Relationships Revealed From 2b-RAD Data

Red sage (Salvia miltiorrhiza) is a widely used medicinal plant for treatment of cardiovascular and cerebrovascular diseases. Because of excessive excavation by huge market demand and habitat loss by human activities, the wild population resources of S. miltiorrhiza have reduced drastically in recent years. Meanwhile, population status of two closely related species S. bowleyana and S. paramiltiorrhiza were in a trend of decreasing due to their potential replacement of S. miltiorrhiza. Particularly, S. paramiltiorrhiza was threatened and endemic to a small region in eastern China. However, to date there has been no conservation genetic research reported for wild S. miltiorrhiza population and its endangered relatives. Assess the wild germplasm diversity for S. miltiorrhiza and its related species would provide fundamental genetic background for cultivation and molecular breeding of this medicinally important species. In the present study, we investigated the genetic diversity, population structure, and intra/inter-specific differentiation of S. miltiorrhiza and above two relatives using 2b-RAD genome-wide genotyping method. By investigating 81 individuals of S. miltiorrhiza, 55 individuals of S. bowleyana and 15 individuals of S. paramiltiorrhiza from 23 locations in China, we obtained 23,928 SNPs in total. A comparatively high genetic diversity was observed in S. miltiorrhiza (π = 0.0788, H e = 0.0783 ± 0.0007). The observed and expected heterozygosity in populations of these three species ranged from 0.0297 to 0.1481 and 0.0251 to 0.831, respectively. Two major lineage groups were detected in the examined S. miltiorrhiza populations. The results indicated that Dabie Mountain as a genetic diversity center of S. miltiorrhiza and possible complex inter-specific genetic exchange/hybridization occurred between S. miltiorrhiza and the two relatives. We suggest that strategic conservation and germplasm preservation should be considered not only for wild populations of S. miltiorrhiza, but also for its related S. bowleyana and S. paramiltiorrhiza.

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

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

Recent progress and new perspectives for diterpenoid biosynthesis in medicinal plants

Diterpenoids, including more than 18,000 compounds, represent an important class of metabolites that encompass both phytohormones and some industrially relevant compounds. These molecules with complex, diverse structures and physiological activities, have high value in the pharmaceutical industry. Most medicinal diterpenoids are extracted from plants. Major advances in understanding the biosynthetic pathways of these active compounds are providing unprecedented opportunities for the industrial production of diterpenoids by metabolic engineering and synthetic biology. Here, we summarize recent developments in the field of diterpenoid biosynthesis from medicinal herbs. An overview of the pathways and known biosynthetic enzymes is presented. In particular, we look at the main findings from the past decade and review recent progress in the biosynthesis of different groups of ringed compounds. We also discuss diterpenoid production using synthetic biology and metabolic engineering strategies, and draw on new technologies and discoveries to bring together many components into a useful framework for diterpenoid production.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

PCPD: Plant cytochrome P450 database and web-based tools for structural construction and ligand docking

Plant cytochrome P450s play key roles in the diversification and functional modification of plant natural products. Although over 200,000 plant P450 gene sequences have been recorded, only seven crystalized P450 genes severely hampered the functional characterization, gene mining and engineering of important P450s. Here, we combined Rosetta homologous modeling and MD-based refinement to construct a high-resolution P450 structure prediction process (PCPCM), which was applied to 181 plant P450s with identified functions. Furthermore, we constructed a ligand docking process (PCPLD) that can be applied for plant P450s virtual screening. 10 examples of virtual screening indicated the process can reduce about 80% screening space for next experimental verification. Finally, we constructed a plant P450 database (PCPD: http://p450.biodesign.ac.cn/), which includes the sequences, structures and functions of the 181 plant P450s, and a web service based on PCPCM and PCPLD. Our study not only developed methods for the P450-specific structure analysis, but also introduced a universal approach that can assist the mining and functional analysis of P450 enzymes.

Genome-Wide Analysis of Terpene Synthase Gene Family in Mentha longifolia and Catalytic Activity Analysis of a Single Terpene Synthase

Terpenoids are a wide variety of natural products and terpene synthase (TPS) plays a key role in the biosynthesis of terpenoids. Mentha plants are rich in essential oils, whose main components are terpenoids, and their biosynthetic pathways have been basically elucidated. However, there is a lack of systematic identification and study of TPS in Mentha plants. In this work, we genome-widely identified and analyzed the TPS gene family in Mentha longifolia, a model plant for functional genomic research in the genus Mentha. A total of 63 TPS genes were identified in the M. longifolia genome sequence assembly, which could be divided into six subfamilies. The TPS-b subfamily had the largest number of genes, which might be related to the abundant monoterpenoids in Mentha plants. The TPS-e subfamily had 18 members and showed a significant species-specific expansion compared with other sequenced Lamiaceae plant species. The 63 TPS genes could be mapped to nine scaffolds of the M. longifolia genome sequence assembly and the distribution of these genes is uneven. Tandem duplicates and fragment duplicates contributed greatly to the increase in the number of TPS genes in M. longifolia. The conserved motifs (RR(X)8W, NSE/DTE, RXR, and DDXXD) were analyzed in M. longifolia TPSs, and significant differentiation was found between different subfamilies. Adaptive evolution analysis showed that M. longifolia TPSs were subjected to purifying selection after the species-specific expansion, and some amino acid residues under positive selection were identified. Furthermore, we also cloned and analyzed the catalytic activity of a single terpene synthase, MlongTPS29, which belongs to the TPS-b subfamily. MlongTPS29 could encode a limonene synthase and catalyze the biosynthesis of limonene, an important precursor of essential oils from the genus Mentha. This study provides useful information for the biosynthesis of terpenoids in the genus Mentha.

Introduction of mutations in plants with prime editing

The emergence of a clustered regularly interspaced short palindromic repeats/CRISPR-associated proteins (CRISPR/Cas) system has had a revolutionary impact on plant biology. However, this system and further developed base editing are limited by their inherent imperfection. Prime editing, a just arrival technology based on CRISPR/Cas, can directly and precisely edit a specified DNA site without double strand breaks and donor DNA by integrating an engineered reverse transcriptase (RT) with a catalytically impaired Cas9 endonuclease and introducing genetic information into prime editing guide RNA (pegRNA). In addition, it has a wider range of editing types than base editing and can install all types of editing theoretically. Prime editing was originally developed in mammalian cells and has recently been applied to plants. Here, we describe the origin of prime editing and compare it with traditional CRISPR/Cas9 and base editing; then, we exemplify it in plants, including strategies and methods. Accordingly, we generate the overall procedures of prime editing to provide instructions for its application. Furthermore, we summarize its improvements in the approach, such as optimizing the length of a primer binding site and RT template, as well as pursuing an optimal nicking site in the unedited sequence. Finally, we discuss the potential impact on domestication and improvement of agricultural crops, sustainable utilization of medicinal plants, cultivation of varieties of horticultural plants, and revelation of the genetic code, in order to offer a reference for the further study and development of prime editing.

The chromosome-based lavender genome provides new insights into Lamiaceae evolution and terpenoid biosynthesis

The aromatic shrub Lavandula angustifolia produces various volatile terpenoids that serve as resources for essential oils and function in plant-insect communication. To better understand the genetic basis of the terpenoid diversity in lavender, we present a high-quality reference genome for the Chinese lavender cultivar "Jingxun 2" using PacBio and Hi-C technologies to anchor the 894.50 Mb genome assembly into 27 pseudochromosomes. In addition to the γ triplication event, lavender underwent two rounds of whole-genome duplication (WGD) during the Eocene-Oligocene (29.6 MYA) and Miocene-Pliocene (6.9 MYA) transitions. As a result of tandem duplications and lineage-specific WGDs, gene families related to terpenoid biosynthesis in lavender are substantially expanded compared to those of five other species in Lamiaceae. Many terpenoid biosynthesis transcripts are abundant in glandular trichomes. We further integrated the contents of ecologically functional terpenoids and coexpressed terpenoid biosynthetic genes to construct terpenoid-gene networks. Typical gene clusters, including TPS-TPS, TPS-CYP450, and TPS-BAHD, linked with compounds that primarily function as attractants or repellents, were identified by their similar patterns of change during flower development or in response to methyl jasmonate. Comprehensive analysis of the genetic basis of the production of volatiles in lavender could serve as a foundation for future research into lavender evolution, phytochemistry, and ecology.

Genome-wide identification and functional characterization of natural antisense transcripts in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most widely used traditional medicines. Natural antisense transcripts (NATs) are a class of long noncoding RNAs that can regulate gene expression. Here, we identified 812 NATs, including 168 cis-NATs and 644 trans-NATs from twelve root, flower, and leaf samples of S. miltiorrhiza using RNA-seq. The expression profiles for 41 of 50 NATs and their sense transcripts (STs) obtained from RNA-Seq were validated using qRT-PCR. The expression profiles of 17 NATs positively correlated with their STs. GO and KEGG pathway analyses mapped the STs for cis-NATs to pathways for biosynthesis of secondary metabolites. We characterized four NATs in detail, including NAT0001, NAT0002, NAT0004, and NAT00023. Their STs are kaurene synthase-like 1 and the homologs of UDP-glucose flavonoid 3-O-glucosyltransferase 6, UDP-glycosyltransferase 90A1, and beta-glucosidase 40, respectively. The first gene is involved in the biosynthesis of bioactive tanshinones, the next two are involved in anthocyanin biosynthesis, whereas the last is involved in phenylpropanoid biosynthesis. Besides, we found seven STs that are potential targets of miRNAs. And we found two miRNAs including miR156a and miR7208, might originate from NATs, NAT0112 and NAT0086. The results suggest that S. miltiorrhiza NATs might interact with STs, produce miRNAs, and be regulated by miRNAs. They potentially play significant regulatory roles in the biosynthesis of bioactive compounds.

Genome-Wide Identification of the TIFY Family in Salvia miltiorrhiza Reveals That SmJAZ3 Interacts With SmWD40-170, a Relevant Protein That Modulates Secondary Metabolism and Development

Salvia miltiorrhiza Bunge (S. miltiorrhiza), a traditional Chinese medicinal herb, contains numerous bioactive components with broad range of pharmacological properties. By increasing the levels of endogenous jasmonate (JA) in plants or treating them with methyl jasmonate (MeJA), the level of tanshinones and salvianolic acids can be greatly enhanced. The jasmonate ZIM (JAZ) proteins belong to the TIFY family, and act as repressors, releasing targeted transcriptional factors in the JA signaling pathway. Herein, we identified and characterized 15 TIFY proteins present in S. miltiorrhiza. Quantitative reverse transcription PCR analysis indicated that the JAZ genes were all constitutively expressed in different tissues and were induced by MeJA treatments. SmJAZ3, which negatively regulates the tanshinones biosynthesis pathway in S. miltiorrhiza and the detailed molecular mechanism is poorly understood. SmJAZ3 acts as a bait protein to capture and identify a WD-repeat containing the protein SmWD40-170. Further molecular and genetic analysis revealed that SmWD40-170 is a positive regulator, promoting the accumulation of secondary metabolites in S. miltiorrhiza. Our study systematically analyzed the TIFY family and speculated a module of the JAZ-WD40 complex provides new insights into the mechanisms regulating the biosynthesis of secondary metabolites in S. miltiorrhiza.

Expansion within the CYP71D subfamily drives the heterocyclization of tanshinones synthesis in Salvia miltiorrhiza

Tanshinones are the bioactive nor-diterpenoid constituents of the Chinese medicinal herb Danshen (Salvia miltiorrhiza). These groups of chemicals have the characteristic furan D-ring, which differentiates them from the phenolic abietane-type diterpenoids frequently found in the Lamiaceae family. However, how the 14,16-epoxy is formed has not been elucidated. Here, we report an improved genome assembly of Danshen using a highly homozygous genotype. We identify a cytochrome P450 (CYP71D) tandem gene array through gene expansion analysis. We show that CYP71D373 and CYP71D375 catalyze hydroxylation at carbon-16 (C16) and 14,16-ether (hetero)cyclization to form the D-ring, whereas CYP71D411 catalyzes upstream hydroxylation at C20. In addition, we discover a large biosynthetic gene cluster associated with tanshinone production. Collinearity analysis indicates a more specific origin of tanshinones in Salvia genus. It illustrates the evolutionary origin of abietane-type diterpenoids and those with a furan D-ring in Lamiaceae.

Genome-wide identification and analysis of DNA methyltransferase and demethylase gene families in Dendrobium officinale reveal their potential functions in polysaccharide accumulation

BACKGROUND

DNA methylation is a conserved and important epigenetic modification involved in the regulation of numerous biological processes, including plant development, secondary metabolism, and response to stresses. However, no information is available regarding the identification of cytosine-5 DNA methyltransferase (C5-MTase) and DNA demethylase (dMTase) genes in the orchid Dendrobium officinale.

RESULTS

In this study, we performed a genome-wide analysis of DoC5-MTase and DodMTase gene families in D. officinale. Integrated analysis of conserved motifs, gene structures and phylogenetic analysis showed that eight DoC5-MTases were divided into four subfamilies (DoCMT, DoDNMT, DoDRM, DoMET) while three DodMTases were divided into two subfamilies (DoDML3, DoROS1). Multiple cis-acting elements, especially stress-responsive and hormone-responsive ones, were found in the promoter region of DoC5-MTase and DodMTase genes. Furthermore, we investigated the expression profiles of DoC5-MTase and DodMTase in 10 different tissues, as well as their transcript abundance under abiotic stresses (cold and drought) and at the seedling stage, in protocorm-like bodies, shoots, and plantlets. Interestingly, most DoC5-MTases were downregulated whereas DodMTases were upregulated by cold stress. At the seedling stage, DoC5-MTase expression decreased as growth proceeded, but DodMTase expression increased.

CONCLUSIONS

These results provide a basis for elucidating the role of DoC5-MTase and DodMTase in secondary metabolite production and responses to abiotic stresses in D. officinale.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.

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.

Pinellia ternata (Thunb.) Breit: A review of its germplasm resources, genetic diversity and active components

ETHNOPHARMACOLOGICAL RELEVANCE

The medicinal plant Pinellia ternata has been widely used in China, Korea, and Japan and has been demonstrated to be highly effective for treating cough, vomiting, infection, and inflammatory diseases. Modern pharmacological investigations have demonstrated its multiple activities, such as antitussive, expectorant, antiemetic, antitumor, antibacterial, and sedative-hypnotic activities.

AIM OF THE REVIEW

This review aims to summarize the information about the biological traits, genetic diversity, active components, and continuous cropping obstacle of P. ternata in order to improve its use.

MATERIALS AND METHODS

In this review, the relevant literature was gathered by using Pinellia ternata, genetic diversity, active components, and continuous cropping obstacle as the keywords from Google Scholar, PubMed, Springer Link, the Wiley online library, SciFinder, SCOPUS, Baidu Scholar, China national knowledge infrastructure (CNKI), and WANFANF DATA (up to April 2020).

RESULTS

P. ternata is the most widely used herb in the Pinellia genus to treat several diseases. The genetic diversity of P. ternata has been extensively studied, and its high genetic diversity level in China has been demonstrated. Modern pharmacological research has indicated that amino acids, alkaloids, and polysaccharides are the main active components supporting P. ternata's medicinal effects. However, an efficient method for determining its active components is still unavailable. The method used to evaluate Pinelliae Rhizoma (PR) quality standards should be further optimized. The continuous cropping obstacle has a significant effect on the quantity and quality of P. ternata. The underlying mechanism of the continuous cropping obstacle needs to be further explored.

CONCLUSIONS

P. ternata has emerged as a valuable source of traditional medicine. Some uses of P. ternata in medicine have been validated by pharmacological investigations. However, a more efficient analytical method should be established to evaluate the quality of PR based on multiple quality markers. Furthermore, high-performance liquid chromatography (HPLC) and DNA barcoding should be introduced to identify the authenticity of PR. In addition, the genes involved in the metabolic synthesis pathways of the main active components, population genetic relationships, the quality control of processed PR, and the continuous cropping obstacle need to be further elucidated. We hope this review will allow for better utilization of this valuable herb.

Comparative Genome Analysis of Scutellaria baicalensis and Scutellaria barbata Reveals the Evolution of Active Flavonoid Biosynthesis

Scutellaria baicalensis (S. baicalensis) and Scutellaria barbata (S. barbata) are common medicinal plants of the Lamiaceae family. Both produce specific flavonoid compounds, including baicalein, scutellarein, norwogonin, and wogonin, as well as their glycosides, which exhibit antioxidant and antitumor activities. Here, we report chromosome-level genome assemblies of S. baicalensis and S. barbata with quantitative chromosomal variation (2n = 18 and 2n = 26, respectively). The divergence of S. baicalensis and S. barbata occurred far earlier than previously reported, and a whole-genome duplication (WGD) event was identified. The insertion of long terminal repeat elements after speciation might be responsible for the observed chromosomal expansion and rearrangement. Comparative genome analysis of the congeneric species revealed the species-specific evolution of chrysin and apigenin biosynthetic genes, such as the S. baicalensis-specific tandem duplication of genes encoding phenylalanine ammonia lyase and chalcone synthase, and the S. barbata-specific duplication of genes encoding 4-CoA ligase. In addition, the paralogous duplication, colinearity, and expression diversity of CYP82D subfamily members revealed the functional divergence of genes encoding flavone hydroxylase between S. baicalensis and S. barbata. Analyzing these Scutellaria genomes reveals the common and species-specific evolution of flavone biosynthetic genes. Thus, these findings would facilitate the development of molecular breeding and studies of biosynthesis and regulation of bioactive compounds.

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.

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

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

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

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

Characterization of Chloroplast Genomes From Two Salvia Medicinal Plants and Gene Transfer Among Their Mitochondrial and Chloroplast Genomes

Salvia species have been widely used as medicinal plants and have played an important role in the treatment and recovery of individuals with COVID-19. In this study, we reported two newly identified whole chloroplast genome sequences of Salvia medicinal plants (Salvia yangii and Salvia miltiorrhiza f. alba) and compared them with those of seven other reported Salvia chloroplast genomes. These were proven to be highly similar in terms of overall size, genome structure, gene content, and gene order. We identified 10 mutation hot spots (trnK-rps16, atpH-atpI, psaA-ycf3, ndhC-trnV, ndhF, rpl32-trnL, ndhG-ndhI, rps15-ycf1, ycf1a, and ycf1b) as candidate DNA barcodes for Salvia. Additionally, we observed the transfer of nine large-sized chloroplast genome fragments, with a total size of 49,895 bp (accounting for 32.97% of the chloroplast genome), into the mitochondrial genome as they shared >97% sequence similarity. Phylogenetic analyses of the whole chloroplast genome provided a high resolution of Salvia. This study will pave the way for the identification and breeding of Salvia medicinal plants and further phylogenetic evolutionary research on them as well.

ABA-responsive transcription factor bZIP1 is involved in modulating biosynthesis of phenolic acids and tanshinones in Salvia miltiorrhiza

Phenolic acids and tanshinones are major bioactive ingredients in Salvia miltiorrhiza, which possess pharmacological activities with great market demand. However, transcriptional regulation of phenolic acid and tanshinone biosynthesis remains poorly understood. Here, a basic leucine zipper transcription factor (TF) named SmbZIP1 was screened from the abscisic acid (ABA)-induced transcriptome library. Overexpression of SmbZIP1 positively promoted phenolic acid biosynthesis by enhancing expression of biosynthetic genes such as cinnamate-4-hydroxylase (C4H1). Furthermore, biochemical experiments revealed that SmbZIP1 bound the G-Box-like1 element in the promoter of the C4H1 gene. Meanwhile, SmbZIP1 inhibited accumulation of tanshinones mainly by suppressing the expression of biosynthetic genes including geranylgeranyl diphosphate synthase (GGPPS) which was confirmed as a target gene by in vitro and in vivo experiments. In contrast, the phenolic acid content was reduced and tanshinone was enhanced in CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9]-mediated knockout lines. In addition, the previously reported positive regulator of tanshinone biosynthesis, SmERF1L1, was found to be inhibited in SmbZIP1 overexpression lines indicated by RNA sequencing, and was proven to be the target of SmbZIP1. In summary, this work uncovers a novel regulator and deepens our understanding of the transcriptional and regulatory mechanisms of phenolic acid and tanshinone biosynthesis, and also sheds new light on metabolic engineering in S. miltiorrhiza.

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.

Genome-wide characterization and expression profiling of MAPK cascade genes in Salvia miltiorrhiza reveals the function of SmMAPK3 and SmMAPK1 in secondary metabolism

BACKGROUND

The contribution of mitogen-activated protein kinase (MAPK) cascades to plant growth and development has been widely studied, but this knowledge has not yet been extended to the medicinal plant Salvia miltiorrhiza, which produces a number of pharmacologically active secondary metabolites.

RESULTS

In this study, we performed a genome-wide survey and identified six MAPKKK kinases (MAPKKKKs), 83 MAPKK kinases (MAPKKKs), nine MAPK kinases (MAPKKs) and 18 MAPKs in the S. miltiorrhiza genome. Within each class of genes, a small number of subfamilies were recognized. A transcriptional analysis revealed differences in the genes' behaviour with respect to both their site of transcription and their inducibility by elicitors and phytohormones. Two genes were identified as strong candidates for playing roles in phytohormone signalling. A gene-to-metabolite network was constructed based on correlation analysis, highlighting the likely involvement of two of the cascades in the synthesis of two key groups of pharmacologically active secondary metabolites: phenolic acids and tanshinones.

CONCLUSION

The data provide insight into the functional diversification and conservation of MAPK cascades in S. miltiorrhiza.