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

Versatile CYP98A enzymes catalyse meta-hydroxylation reveals diversity of salvianolic acids biosynthesis

Salvianolic acids (SA), such as rosmarinic acid (RA), danshensu (DSS), and their derivative salvianolic acid B (SAB), etc. widely existed in Lamiaceae and Boraginaceae families, are of interest due to medicinal properties in the pharmaceutical industries. Hundreds of studies in past decades described that 4-coumaroyl-CoA and 4-hydroxyphenyllactic acid (4-HPL) are common substrates to biosynthesize SA with participation of rosmarinic acid synthase (RAS) and cytochrome P450 98A (CYP98A) subfamily enzymes in different plants. However, in our recent study, several acyl donors and acceptors included DSS as well as their ester-forming products all were determined in SA-rich plants, which indicated that previous recognition to SA biosynthesis is insufficient. Here, we used Salvia miltiorrhiza, a representative important medicinal plant rich in SA, to elucidate the diversity of SA biosynthesis. Various acyl donors as well as acceptors are catalysed by SmRAS to form precursors of RA and two SmCYP98A family members, SmCYP98A14 and SmCYP98A75, are responsible for different positions' meta-hydroxylation of these precursors. SmCYP98A75 preferentially catalyses C-3' hydroxylation, and SmCYP98A14 preferentially catalyses C-3 hydroxylation in RA generation. In addition, relative to C-3' hydroxylation of the acyl acceptor moiety in RA biosynthesis, SmCYP98A75 has been verified as the first enzyme that participates in DSS formation. Furthermore, SmCYP98A enzymes knockout resulted in the decrease and overexpression leaded to dramatic increase of SA accumlation. Our study provides new insights into SA biosynthesis diversity in SA-abundant species and versatility of CYP98A enzymes catalytic preference in meta-hydroxylation reactions. Moreover, CYP98A enzymes are ideal metabolic engineering targets to elevate SA content.

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.

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.

DNA methylation regulates biosynthesis of tanshinones and phenolic acids during growth of Salvia miltiorrhiza

DNA methylation plays a crucial role in the regulation of plant growth and the biosynthesis of secondary metabolites. Danshen (Salvia miltiorrhiza) is a valuable Chinese herbal medicine commonly used to treat cardiovascular diseases; its active ingredients are tanshinones and phenolic acids, which primarily accumulate in roots. Here, we conducted a targeted metabolic analysis of S. miltiorrhiza roots at 3 distinct growth stages: 40 d old (r40), 60 d old (r60), and 90 d old (r90). The contents of tanshinones (cryptotanshinone, tanshinone I, tanshinone IIA, and rosmariquinone) and phenolic acids (rosmarinic acid and salvianolic acid B) gradually increased during plant development. Whole-genome bisulfite sequencing and transcriptome sequencing of roots at the 3 growth stages revealed an increased level of DNA methylation in the CHH context (H represents A, T, or C) context at r90 compared with r40 and r60. Increased DNA methylation levels were associated with elevated expression of various genes linked to epigenetic regulations, including CHROMOMETHYLASE2 (SmCMT2), Decrease in DNA Methylation 1 (SmDDM1), Argonaute 4 (SmAGO4), and DOMAINS REARRANGED METHYLTRANSFERASE 1 (SmDRM1). Moreover, expression levels of many genes involved in tanshinone and salvianolic acid biosynthesis, such as copalyldiphosphate synthase 5 (SmCPS5), cytochrome P450-related enzyme (SmCYP71D464), geranylgeranyl diphosphate synthase (SmGGPPS1), geranyl diphosphate synthase (SmGPPS), hydroxyphenylpyruvate reductase (SmHPPR), and hydroxyphenylpyruvate dioxygenase (SmHPPD), were altered owing to hyper-methylation, indicating that DNA methylation plays an important role in regulating tanshinone and phenolic acid accumulation. Our data shed light on the epigenetic regulation of root growth and the biosynthesis of active ingredients in S. miltiorrhiza, providing crucial clues for further improvement of active compound production via molecular breeding in S. miltiorrhiza.

A haplotype-resolved gap-free genome assembly provides novel insight into monoterpenoid diversification in Mentha suaveolens 'Variegata'

Mentha is a commonly used spice worldwide, which possesses medicinal properties and fragrance. These characteristics are conferred, at least partially, by essential oils such as menthol. In this study, a gap-free assembly with a genome size of 414.3 Mb and 31,251 coding genes was obtained for Mentha suaveolens 'Variegata'. Based on its high heterozygosity (1.5%), two complete haplotypic assemblies were resolved, with genome sizes of 401.9 and 405.7 Mb, respectively. The telomeres and centromeres of each haplotype were almost fully annotated. In addition, we detected a total of 41,135 structural variations. Enrichment analysis demonstrated that genes involved in terpenoid biosynthesis were affected by these structural variations. Analysis of volatile metabolites showed that M. suaveolens mainly produces piperitenone oxide rather than menthol. We identified three genes in the M. suaveolens genome which encode isopiperitenone reductase (ISPR), a key rate-limiting enzyme in menthol biosynthesis. However, the transcription levels of ISPR were low. Given that other terpenoid biosynthesis genes were expressed, M. suaveolens ISPRs may account for the accumulation of piperitenone oxide in this species. The findings of this study may provide a valuable resource for improving the detection rate and accuracy of genetic variants, thereby enhancing our understanding of their impact on gene function and expression. Moreover, our haplotype-resolved gap-free genome assembly offers novel insights into molecular marker-assisted breeding of Mentha.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Strategies on biosynthesis and production of bioactive compounds in medicinal plants

Medicinal plants are a valuable source of essential medicines and herbal products for healthcare and disease therapy. Compared with chemical synthesis and extraction, the biosynthesis of natural products is a very promising alternative for the successful conservation of medicinal plants, and its rapid development will greatly facilitate the conservation and sustainable utilization of medicinal plants. Here, we summarize the advances in strategies and methods concerning the biosynthesis and production of natural products of medicinal plants. The strategies and methods mainly include genetic engineering, plant cell culture engineering, metabolic engineering, and synthetic biology based on multiple "OMICS" technologies, with paradigms for the biosynthesis of terpenoids and alkaloids. We also highlight the biosynthetic approaches and discuss progress in the production of some valuable natural products, exemplifying compounds such as vindoline (alkaloid), artemisinin and paclitaxel (terpenoids), to illustrate the power of biotechnology in medicinal plants.

Reference genome of the nutrition-rich orphan crop chia (Salvia hispanica) and its implications for future breeding

Chia (Salvia hispanica L.) is one of the most popular nutrition-rich foods and pseudocereal crops of the family Lamiaceae. Chia seeds are a rich source of proteins, polyunsaturated fatty acids (PUFAs), dietary fibers, and antioxidants. In this study, we present the assembly of the chia reference genome, which spans 303.6 Mb and encodes 48,090 annotated protein-coding genes. Our analysis revealed that ~42% of the chia genome harbors repetitive content, and identified ~3 million single nucleotide polymorphisms (SNPs) and 15,380 simple sequence repeat (SSR) marker sites. By investigating the chia transcriptome, we discovered that ~44% of the genes undergo alternative splicing with a higher frequency of intron retention events. Additionally, we identified chia genes associated with important nutrient content and quality traits, such as the biosynthesis of PUFAs and seed mucilage fiber (dietary fiber) polysaccharides. Notably, this is the first report of in-silico annotation of a plant genome for protein-derived small bioactive peptides (biopeptides) associated with improving human health. To facilitate further research and translational applications of this valuable orphan crop, we have developed the Salvia genomics database (SalviaGDB), accessible at https://salviagdb.org.

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.

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.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

The First Genome-Wide Mildew Locus O Genes Characterization in the Lamiaceae Plant Family

Powdery mildew (PM) is a widespread plant disease that causes significant economic losses in thousands crops of temperate climates, including Lamiaceae species. Multiple scientific studies describe a peculiar form of PM-resistance associated at the inactivation of specific members of the Mildew Locus O (MLO) gene family, referred to as mlo-resistance. The characterization of Lamiaceae MLO genes, at the genomic level, would be a first step toward their potential use in breeding programs. We carried out a genome-wide characterization of the MLO gene family in 11 Lamiaceae species, providing a manual curated catalog of 324 MLO proteins. Evolutionary history and phylogenetic relationships were studied through maximum likelihood analysis and motif patter reconstruction. Our approach highlighted seven different clades diversified starting from an ancestral MLO domain pattern organized in 18 highly conserved motifs. In addition, 74 Lamiaceae putative PM susceptibility genes, clustering in clade V, were identified. Finally, we performed a codon-based evolutionary analysis, revealing a general high level of purifying selection in the eleven Lamiaceae MLO gene families, and the occurrence of few regions under diversifying selection in candidate susceptibility factors. The results of this work may help to address further biological questions concerning MLOs involved in PM susceptibility. In follow-up studies, it could be investigated whether the silencing or loss-of-function mutations in one or more of these candidate genes may lead to PM resistance.

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.

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.

Genome-Wide Analysis of the MADS-box Gene Family and Expression Analysis during Anther Development in Salvia miltiorrhiza

MADS-box genes constitute a large family of transcription factors that play important roles in plant growth and development. However, our understanding of MADS-box genes involved in anther development and male sterility in Salvia miltiorrhiza is still limited. In this study, 63 MADS-box genes were identified from the genome of the male sterility ecotype Sichuan S. miltiorrhiza (S. miltiorrhiza_SC) unevenly distributed among eight chromosomes. Phylogenetic analysis classified them into two types and 17 subfamilies. They contained 1 to 12 exons and 10 conserved motifs. Evolution analysis showed that segmental duplication was the main force for the expansion of the SmMADS gene family, and duplication gene pairs were under purifying selection. Cis-acting elements analysis demonstrated that the promoter of SmMADS genes contain numerous elements associated with plant growth and development, plant hormones, and stress response. RNA-seq showed that the expression levels of B-class and C-class SmMADS genes were highly expressed during anther development, with SmMADS11 likely playing an important role in regulating anther development and male fertility in S. miltiorrhiza_SC. Overall, this study provides a comprehensive analysis of the MADS-box gene family in S. miltiorrhiza, shedding light on their potential role in anther development and male sterility.

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.

Plant terpene specialized metabolism: complex networks or simple linear pathways?

From the perspectives of pathway evolution, discovery and engineering of plant specialized metabolism, the nature of the biosynthetic routes represents a critical aspect. Classical models depict biosynthesis typically from an end-point angle and as linear, for example, connecting central and specialized metabolism. As the number of functionally elucidated routes increased, the enzymatic foundation of complex plant chemistries became increasingly well understood. The perception of linear pathway models has been severely challenged. With a focus on plant terpenoid specialized metabolism, we review here illustrative examples supporting that plants have evolved complex networks driving chemical diversification. The completion of several diterpene, sesquiterpene and monoterpene routes shows complex formation of scaffolds and their subsequent functionalization. These networks show that branch points, including multiple sub-routes, mean that metabolic grids are the rule rather than the exception. This concept presents significant implications for biotechnological production.

Genome-Wide Characterization and Expression of the Hsf Gene Family in Salvia miltiorrhiza (Danshen) and the Potential Thermotolerance of SmHsf1 and SmHsf7 in Yeast

Salvia miltiorrhiza Bunge (Danshen) is a traditional Chinese herb with significant medicinal value. The yield and quality of Danshen are greatly affected by climatic conditions, in particular high temperatures. Heat shock factors (Hsfs) play important regulatory roles in plant response to heat and other environmental stresses. However, little is currently known about the role played by the Hsf gene family in S. miltiorrhiza. Here, we identified 35 SmHsf genes and classified them into three major groups: SmHsfA (n = 22), SmHsfB (n = 11), and SmHsfC (n = 2) using phylogenetic analysis. The gene structure and protein motifs were relatively conserved within subgroups but diverged among the different groups. The expansion of the SmHsf gene family was mainly driven by whole-genome/segmental and dispersed gene duplications. The expression profile of SmHsfs in four distinct organs revealed its members (23/35) are predominantly expressed in the root. The expression of a large number of SmHsfs was regulated by drought, ultraviolet, heat and exogenous hormones. Notably, the SmHsf1 and SmHsf7 genes in SmHsfB2 were the most responsive to heat and are conserved between dicots and monocots. Finally, heterologous expression analysis showed that SmHsf1 and SmHsf7 enhance thermotolerance in yeast. Our results provide a solid foundation for further functional investigation of SmHsfs in Danshen plants as a response to abiotic stresses.

A near complete genome assembly of chia assists in identification of key fatty acid desaturases in developing seeds

Chia is an annual crop whose seeds have the highest content of α-linolenic acid (ALA) of any plant known to date. We generated a high-quality assembly of the chia genome using circular consensus sequencing (CCS) of PacBio. The assembled six chromosomes are composed of 21 contigs and have a total length of 361.7 Mb. Genome annotation revealed a 53.5% repeat content and 35,850 protein-coding genes. Chia shared a common ancestor with Salvia splendens ~6.1 million years ago. Utilizing the reference genome and two transcriptome datasets, we identified candidate fatty acid desaturases responsible for ALA biosynthesis during chia seed development. Because the seed of S. splendens contains significantly lower proportion of ALA but similar total contents of unsaturated fatty acids, we suggest that strong expression of two ShFAD3 genes are critical for the high ALA content of chia seeds. This genome assembly will serve as a valuable resource for breeding, comparative genomics, and functional genomics studies of chia.

The chromosome-scale assembly of the Salvia rosmarinus genome provides insight into carnosic acid biosynthesis

Rosemary (Salvia rosmarinus) is considered a sacred plant because of its special fragrance and is commonly used in cooking and traditional medicine. Here, we report a high-quality chromosome-level assembly of the S. rosmarinus genome of 1.11 Gb in size; the genome has a scaffold N50 value of 95.5 Mb and contains 40 701 protein-coding genes. In contrast to other diploid Labiataceae, an independent whole-genome duplication event occurred in S. rosmarinus at approximately 15 million years ago. Transcriptomic comparison of two S. rosmarinus cultivars with contrasting carnosic acid (CA) content revealed 842 genes significantly positively associated with CA biosynthesis in S. rosmarinus. Many of these genes have been reported to be involved in CA biosynthesis previously, such as genes involved in the mevalonate/methylerythritol phosphate pathways and CYP71-coding genes. Based on the genomes and these genes, we propose a model of CA biosynthesis in S. rosmarinus. Further, comparative genome analysis of the congeneric species revealed the species-specific evolution of CA biosynthesis genes. The genes encoding diterpene synthase and the cytochrome P450 (CYP450) family of CA synthesis-associated genes form a biosynthetic gene cluster (CPSs-KSLs-CYP76AHs) responsible for the synthesis of leaf and root diterpenoids, which are located on S. rosmarinus chromosomes 1 and 2, respectively. Such clustering is also observed in other sage (Salvia) plants, thus suggesting that genes involved in diterpenoid synthesis are conserved in the Labiataceae family. These findings provide new insights into the synthesis of aromatic terpenoids and their regulation.

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.

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.

Identification, classification, and expression profile analysis of heat shock transcription factor gene family in Salvia miltiorrhiza

In response to abiotic stresses, transcription factors are essential. Heat shock transcription factors (HSFs), which control gene expression, serve as essential regulators of plant growth, development, and stress response. As a model medicinal plant, Salvia miltiorrhiza is a crucial component in the treatment of cardiovascular illnesses. But throughout its growth cycle, S.miltiorrhiza is exposed to a series of abiotic challenges, including heat and drought. In this study, 35 HSF genes were identified based on genome sequencing of Salvia miltiorrhiza utilizing bioinformatics techniques. Additionally, 35 genes were classified into three groups by phylogeny and gene structural analysis, comprising 22 HSFA, 11 HSFB, and two HSFC. The distribution and sequence analysis of motif showed that SmHSFs were relatively conservative. In SmHSF genes, analysis of the promoter region revealed the presence of many cis-acting elements linked to stress, hormones, and growth and development, suggesting that these factors have regulatory roles. The majority of SmHSFs were expressed in response to heat and drought stress, according to combined transcriptome and real-time quantitative PCR (qRT-PCR) analyses. In conclusion, this study looked at the SmHSF gene family using genome-wide identification, evolutionary analysis, sequence characterization, and expression analysis. This research serves as a foundation for further investigations into the role of HSF genes and their molecular mechanisms in plant stress responses.

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.

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

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

Integration of Repeatomic and Cytogenetic Data on Satellite DNA for the Genome Analysis in the Genus Salvia (Lamiaceae)

Within the complicated and controversial taxonomy of cosmopolitan genus Salvia L. (Lamiaceae) are valuable species Salvia officinalis L. and Salvia sclarea L., which are important for the pharmaceutical, ornamental horticulture, food, and perfume industries. Genome organization and chromosome structure of these essential oil species remain insufficiently studied. For the first time, the comparative repeatome analysis of S. officinalis and S. sclarea was performed using the obtained NGS data, RepeatExplorer/TAREAN pipelines and FISH-based chromosome mapping of the revealed satellite DNA families (satDNAs). In repeatomes of these species, LTR retrotransposons made up the majority of their repetitive DNA. Interspecific variations in genome abundance of Class I and Class II transposable elements, ribosomal DNA, and satellite DNA were revealed. Four (S. sclarea) and twelve (S. officinalis) putative satDNAs were identified. Based on patterns of chromosomal distribution of 45S rDNA; 5S rDNA and the revealed satDNAs, karyograms of S. officinalis and S. sclarea were constructed. Promising satDNAs which can be further used as chromosome markers to assess inter- and intraspecific chromosome variability in Salvia karyotypes were determined. The specific localization of homologous satDNA and 45S rDNA on chromosomes of the studied Salvia species confirmed their common origin, which is consistent with previously reported molecular phylogenetic data.

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

ETHNOPHARMACOLOGICAL RELEVANCE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

A chromosome-level genome assembly of chia provides insights into high omega-3 content and coat color variation of its seeds

Chia (Salvia hispanica) is a functional food crop for humans. Although its seeds contain high omega-3 fatty acids, the seed yield of chia is still low. Genomic resources available for this plant are limited. We report the first high-quality chromosome-level genome sequence of chia. The assembled genome size was 347.6 Mb and covered 98.1% of the estimated genome size. A total of 31 069 protein-coding genes were predicted. The absence of recent whole-genome duplication and the relatively low intensity of transposable element expansion in chia compared to its sister species contribute to its small genome size. Transcriptome sequencing and gene duplication analysis reveal that the expansion of the fab2 gene family is likely to be related to the high content of omega-3 in seeds. The white seed coat color is determined by a single locus on chromosome 4. This study provides novel insights into the evolution of Salvia species and high omega-3 content, as well as valuable genomic resources for genetic improvement of important commercial traits of chia and its related species.

Loss of anthocyanidin synthase gene is associated with white flowers of Salvia miltiorrhiza Bge. f. alba, a natural variant of S. miltiorrhiza

SmANS deletion leads to white flower mutation in Salvia miltiorrhiza. SmANS deletion leads to white flower mutation in Salvia miltiorrhiza. Abstract Salvia miltiorrhiza is an essential traditional Chinese medicine (TCM) with purple flowers, and S. miltiorrhiza Bge. f. alba is a unique intraspecific variation with white flowers. The molecular mechanism of flower color formation in S. miltiorrhiza will provide vital information for the variation and evolution. Here, we performed HPLC, transcriptomic, and re-sequencing analyses of purple-flowered S. miltiorrhiza line 'Zihua105' (ZH105) and white-flowered S. miltiorrhiza Bge. f. alba line 'Baihua18' (BH18). Delphinidin was the most anthocyanidin in ZH105, which become the main different between ZH105 vs. BH18 flowers. Transcriptome analysis revealed 299 differentially expressed genes (DEGs). SmANS, the anthocyanidin synthase gene in the down-stream anthocyanin biosynthesis pathway, was significantly expressed in ZH105 corollas, suggesting it might play a key role in white petal formation. Whole-genome re-sequencing revealed that a 6.75 kb segment located on chromosome 5, which contains the complete sequence of the SmANS genes, was lost in BH18 and another S. miltiorrhiza Bge. f. alba line. In contrast, the other five purple-flowered S. miltiorrhiza lines both possessed this segment. Further molecular marker identification also confirmed that wild S. miltiorrhiza Bge. f. alba lines lost regions that contained a complete or important part of SmANS sequences. Subsequently, the research showed that the deletion mutant of SmANS genes resulted in the natural white flower color variant of S. miltiorrhiza.

Identification of Abietane-Type Diterpenoids and Phenolic Acids Biosynthesis Genes in Salvia apiana Jepson Through Full-Length Transcriptomic and Metabolomic Profiling

Salvia apiana (S. apiana) Jepson is a medicinal plant that is frequently used by the Chumash Indians in southern California as a diaphoretic, calmative, diuretic, or antimicrobial agent. Abietane-type diterpenoids (ATDs) and phenolic acids (PAs) are the main bioactive ingredients in S. apiana. However, few studies have looked into the biosynthesis of ATDs and PAs in S. apiana. In this study, using metabolic profiling focused on the ATDs and PAs in the roots and leaves of S. apiana, we found a distinctive metabolic feature with all-around accumulation of ATDs, but absence of salvianolic acid B. To identify the candidate genes involved in these biosynthesis pathways, full-length transcriptome was performed by PacBio single-molecule real-time (SMRT) sequencing. A total of 50 and 40 unigenes were predicted to be involved in ATDs and PAs biosynthesis, respectively. Further transcriptional profile using Illumina HiSeq sequencing showed that the transcriptional variations of these pathways were consistent with the accumulation patterns of corresponding metabolites. A plant kingdom-wide phylogenetic analysis of cytochromes (CYPs) identified two CYP76AK and two CYP76AH subfamily genes that might contribute for the specific ATDs biosynthesis in S. apiana. We also noticed that the clade VII laccase gene family was significantly expanded in Salvia miltiorrhiza compared with that of S. apiana, indicating their involvements in the formation of salvianolic acid B. In conclusion, our results will enable the further understanding of ATDs and PAs biosynthesis in S. apiana and Salvia genus.

SmGDB: genome database of Salvia miltiorrhiza, an important TCM Plant

BACKGROUND

Salvia miltiorrhiza is an important traditional Chinese medicinal (TCM) plant and a model plant in the genetic study of TCM. A series of omics related to Danshen have been published. Integrating, managing, storing, and sharing data has become an urgent problem to be solved in S. miltiorrhiza genetic studies.

OBJECTIVES

The genome database is the link for the exchange, acquisition, and use of different omics data between data producers and users, maximizing value and utilization of data.

METHODS

The genome database included DSS3 genome and five RNA-Seq data. The back-end performs data search and retrieval through the LAMP (Linux, Apache, MySQL, PHP) framework.

RESULTS

Here, we present SmGDB (S. miltiorrhiza genome database; http://8.140.162.85/ ), which houses the latest version of genome sequence and annotation data for S. miltiorrhiza, combining three unpublished RNA-Seq data from our group and two released RNA-Seq data. We also identified a novel gene cluster including seven CYP71D genes involved in the tanshinone synthesis pathway based on genome sequences and expression data. Besides, SmGDB provides user-friendly web interfaces for querying and browsing gene annotation, structure, location, and expression profiles for concerned genes. Popular bioinformatics tools such as 'BLAST', 'Search', 'Heatmap', 'JBrowse', etc., were also provided in SmGDB.

CONCLUSIONS

SmGDB will provide utility for characterizing the structure of the S. miltiorrhiza genome and better understanding gene functions and biological processes underlying complex secondary metabolism in Danshen.

Allelopathic Properties of Lamiaceae Species: Prospects and Challenges to Use in Agriculture

Herbicide resistance due to the increasing reliance on herbicides is a near-term challenge for the world’s agriculture. This has led to a desire to develop new herbicides with a novel mode of action, to address resistance in weed species. Lamiaceae, a large dicotyledonous plant family, is very well known for the multitudinous pharmacological and toxicological properties of its member species. Moreover, many species of this family are significant for their allelopathic activity in natural and laboratory settings. Thus, plants in Lamiaceae have the potential to be sources of alternative herbicides. However, gaps in our knowledge need to be addressed prior to adopting these allelopathic activities in agriculture. Therefore, we review the existing state of knowledge about the Lamiaceae family, the reported allelopathic properties of plant extracts, and their isolated allelochemicals under laboratory, greenhouse, and field conditions. In addition, we offer a perspective on existing challenges and future opportunities for adopting the allelopathic properties of Lamiaceae plant species for green agriculture.

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.

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.

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.

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.