Genome-Wide Identification and Characterization of Salvia miltiorrhiza Laccases Reveal Potential Targets for Salvianolic Acid B Biosynthesis

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.

Enhancement of specialized metabolites using CRISPR/Cas gene editing technology in medicinal plants

Plants are the richest source of specialized metabolites. The specialized metabolites offer a variety of physiological benefits and many adaptive evolutionary advantages and frequently linked to plant defense mechanisms. Medicinal plants are a vital source of nutrition and active pharmaceutical agents. The production of valuable specialized metabolites and bioactive compounds has increased with the improvement of transgenic techniques like gene silencing and gene overexpression. These techniques are beneficial for decreasing production costs and increasing nutritional value. Utilizing biotechnological applications to enhance specialized metabolites in medicinal plants needs characterization and identification of genes within an elucidated pathway. The breakthrough and advancement of CRISPR/Cas-based gene editing in improving the production of specific metabolites in medicinal plants have gained significant importance in contemporary times. This article imparts a comprehensive recapitulation of the latest advancements made in the implementation of CRISPR-gene editing techniques for the purpose of augmenting specific metabolites in medicinal plants. We also provide further insights and perspectives for improving metabolic engineering scenarios in medicinal plants.

Transcriptome Analysis of Salvia miltiorrhiza under Drought Stress

Phenolic acids are one of the major secondary metabolites accumulated in Salvia miltiorrhiza with various pharmacological activities. Moderate drought stress can promote the accumulation of phenolic acids in S. miltiorrhiza, while the mechanism remains unclear. Therefore, we performed transcriptome sequencing of S. miltiorrhiza under drought treatment. A total of 47,169 unigenes were successfully annotated in at least one of the six major databases. Key enzyme genes involved in the phenolic acid biosynthetic pathway, including SmPAL, SmC4H, Sm4CL, SmTAT, SmHPPR, SmRAS and SmCYP98A14, were induced. Unigenes annotated as laccase correlated with SmRAS and SmCYP98A14 were analyzed, and seven candidates that may be involved in the key step of SalB biosynthesis by RA were obtained. A total of 15 transcription factors significantly up-regulated at 2 h and 4 h potentially regulating phenolic acid biosynthesis were screened out. TRINITY_DN14213_c0_g1 (AP2/ERF) significantly transactivated the expression of SmC4H and SmRAS, suggesting its role in the regulation of phenolic acid biosynthesis. GO and KEGG enrichment analysis of differential expression genes showed that phenylpropanoid biosynthesis and plant hormone signal transduction were significantly higher. The ABA-dependent pathway is essential for resistance to drought and phenolic acid accumulation. Expression patterns in drought and ABA databases showed that four PYLs respond to both drought and ABA, and three potential SnRK2 family members were annotated and analyzed. The present study presented a comprehensive transcriptome analysis of S. miltiorrhiza affected by drought, which provides a rich source for understanding the molecular mechanism facing abiotic stress in S. miltiorrhiza.

miR397-LACs mediated cadmium stress tolerance in Arabidopsis thaliana

Cadmium (Cd) is a non-essential heavy metal, assimilated in plant tissue with other nutrients, disturbing the ions' homeostasis in plants. The plant develops different mechanisms to tolerate the hazardous environmental effects of Cd. Recently studies found different miRNAs that are involved in Cd stress. In the current study, miR397 mutant lines were constructed to explore the molecular mechanisms of miR397 underlying Cd tolerance. Compared with the genetically modified line of overexpressed miR397 (artificial miR397, amiR397), the lines of downregulated miR397 (Short Tandem Target Mimic miR397, STTM miR397) showed more substantial Cd tolerance with higher chlorophyll a & b, carotenoid and lignin content. ICP-OES revealed higher cell wall Cd and low total Cd levels in STTM miR397 than in the wild-type and amiR397 plants.Further, the STTM plants produced fewer reactive oxygen species (ROS) and lower activity of antioxidants enzymes (e.g., catalase [CAT], malondialdehyde [MDA]) compared with amiR397 and wild-type plants after stress, indicating that silencing the expression of miR397 can reduce oxidative damage. In addition, the different family transporters' gene expression was much higher in the amiR397 plants than in the wild type and STTM miRNA397. Our results suggest that miR397 plays a role in Cd tolerance in Arabidopsis thaliana. Overexpression of miR397 could decrease Cd tolerance in plants by regulating the expression of LAC 2/4/17, changing the lignin content, which may play an important role in inducing different stress-tolerant mechanisms and protecting the cell from a hazardous condition. This study provides a basis to elucidate the functions of miR397 and the Cd stress tolerance mechanism in Arabidopsis thaliana.

Molecular identification of a laccase that catalyzes the oxidative coupling of a hydroxycinnamic acid amide for hordatine biosynthesis in barley

Plants produce dimerized phenolic compounds as secondary metabolites. Hordatine A (HA), a dehydrodimer of p-coumaroylagmatine (pCA), is an antifungal compound accumulated at high levels in young barley (Hordeum vulgare) seedlings. The enzyme responsible for the oxidative dimerization of pCA, which is the final step of the hordatine biosynthetic pathway, has not been identified. In this study, we first verified the presence of this enzyme activity in the crude extract of barley seedlings. Because the enzyme activity was not dependent on H2 O2 , the responsible enzyme was not peroxidase, which was previously implicated in HA biosynthesis. The analysis of the dissection lines of wheat (Triticum aestivum) carrying aberrant barley 2H chromosomes detected HA in the wheat lines carrying the distal part of the 2H short arm. This chromosomal region contains two laccase genes (HvLAC1 and HvLAC2) that are highly expressed at the seedling stage and may encode enzymes that oxidize pCA during the formation of HA. Changes in the HvLAC transcript levels coincided with the changes in the HA biosynthesis-related enzyme activities in the crude extract and the HA content in barley seedlings. Moreover, HvLAC genes were heterologously expressed in Nicotiana benthamiana leaves and in bamboo (Phyllostachys nigra) suspension cells and HA biosynthetic activities were detected in the crude extract of transformed N. benthamiana leaves and bamboo suspension cells. The HA formed by the enzymatic reaction had the same stereo-configuration as the naturally occurring HA. These results demonstrate that HvLAC enzymes mediate the oxidative coupling of pCA during HA biosynthesis.

Methyl jasmonate elicitation effect on the metabolic profile of cambial meristematic cells culture derived from sweet basil (Ocimum basilicum L.) in relation to antioxidant activity: Untargeted metabolomics study in a time-based approach

The undifferentiated cambial meristematic cell (CMC) has been recognized as a value-added production platform for plant natural products in comparison to the dedifferentiated plant cell line (DDC). In a time-based approach at 0, 24, 48, and 72 h, the present study aimed at investigating the phytochemical metabolome of methyl jasmonate (MeJA)-elicited CMC cultures derived from sweet basil (Ocimum basilicum L.), including primary and secondary metabolites analyzed using GC/TOF-MS post-silylation and RP-UPLC-C₁₈-FT-MS/MS, respectively, as well as the analysis of aroma composition using headspace SPME-GC-MS. The results revealed a stress response in primary metabolism manifested by an increase in amino and organic acids reaching their maximum levels after 48 (1.3-fold) and 72 (1.7-fold) h, respectively. In addition, phenolic acids (e.g., sagerinic acid, rosmarinic acid, and 3-O-methylrosmarinic acid) followed by flavonoid aglycones (e.g., salvigenin and 5,6,4'-trihydroxy-7,3'-dimethoxyflavone) were the most abundant with prominent increases at 48 (1.2-fold) and 72 (2.1-fold) h, respectively. The aroma was intensified by the elicitation along the time, especially after 48 and 72 h. Furthermore, multivariate data analyses, including principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) confirmed elicitation effect, especially post 48 and 72 h. The study further assessed the effect of MeJA elicitation on the antioxidant and polyphenolic content. The cultures at 48 h demonstrated a significant (p < 0.05) antioxidant activity concurrently with correlation with total polyphenolic content using Pearson's correlation. Our study provides new insights to the elicitation impact on primary and secondary metabolism, in addition to aroma profile, to orchestrate the stress response and in relation to antioxidant effect.

Overexpression of SmLAC25 promotes lignin accumulation and decreases salvianolic acid content in Salvia miltiorrhiza

Laccase (LAC) is a blue multicopper oxidase that contains four copper ions, which is involved in lignin polymerization and flavonoid biosynthesis in plants. Although dozens of LAC genes have been identified in Salvia miltiorrhiza Bunge (a model medicinal plant), most have not been functionally characterized. Here, we explored the expression patterns and the functionality of SmLAC25 in S. miltiorrhiza. SmLAC25 has a higher expression level in roots and responds to methyl jasmonate, auxin, abscisic acid, and gibberellin stimuli. The SmLAC25 protein is localized in the cytoplasm and chloroplasts. Recombinant SmLAC25 protein could oxidize coniferyl alcohol and sinapyl alcohol, two monomers of G-lignin and S-lignin. To investigate its function, we generated SmLAC25-overexpressed S. miltiorrhiza plantlets and hairy roots. The lignin content increased significantly in all SmLAC25-overexpressed plantlets and hairy roots, compared with the controls. However, the concentrations of rosmarinic acid and salvianolic acid B decreased significantly in all the SmLAC25-overexpressed lines. Further studies revealed that the transcription levels of some key enzyme genes in the lignin synthesis pathway (e.g., SmCCR and SmCOMT) were significantly improved in the SmLAC25-overexpressed lines, while the expression levels of multiple enzyme genes in the salvianolic acid biosynthesis pathway were inhibited. We speculated that the overexpression of SmLAC25 promoted the metabolic flux of lignin synthesis, which resulted in a decreased metabolic flux to the salvianolic acid biosynthesis pathway.

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.

De novo Biosynthesis of Salvianolic Acid B in Saccharomyces cerevisiae Engineered with the Rosmarinic Acid Biosynthetic Pathway

Salvianolic acid B (SAB), also named lithospermic acid B, belongs to a class of water-soluble phenolic acids, originating from plants such as Salvia miltiorrhiza. SAB exhibits a variety of biological activities and has been clinically used to treat cardio- and cerebrovascular diseases and also has great potential as a health care product and medicine for other disorders. However, its biosynthetic pathway has not been completely elucidated. Here, we report the de novo biosynthesis of SAB in Saccharomyces cerevisiae engineered with the heterologous rosmarinic acid (RA) biosynthetic pathway. The created pathway contains seven genes divided into three modules on separate plasmids, pRS424-FjTAL-Sm4CL2, pRS425-SmTAT-SmHPPR or pRS425-SmTAT-CbHPPR, and pRS426-SmRAS-CbCYP-CbCPR. These three modules were cotransformed into S. cerevisiae, resulting in the recombinant strains YW-44 and YW-45. Incubation of the recombinant strains in a basic medium without supplementing any substrates yielded 34 and 30 μg/L of SAB. The findings in this study indicate that the created heterologous RA pathway cooperates with the native metabolism of S. cerevisiae to enable the de novo biosynthesis of SAB. This provides a novel insight into a biosynthesis mechanism of SAB and also lays the foundation for the production of SAB using microbial cell factories.

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

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

HPTLC and ATR/FTIR Characterization of Antioxidants in Different Rosemary Extracts

The effect of spontaneous fermentation by lactic acid bacteria on the extraction yield of bioactive compounds and antioxidant activity from rosemary leaf extracts was investigated using high-performance thin-layer chromatography (HPTLC). Brining and spontaneous fermentation with lactic acid bacteria more than doubled extraction of polyphenolics and antioxidants from the rosemary leaves. The results show that lactic acid fermentation enhances antioxidant activity in extracts by increasing the total phenolic content but does not increase extraction of phytosterols. Increased extraction of phenolic oxidants during fermentation assisted extraction, results from the in situ generated natural eutectic solvent from the plant sample. ATR-FTIR spectra from the bioactive bands suggests that this increased antioxidant activity is associated with increased extraction of rosmarinic acid, depolymerised lignin, abietane diterpenoids and 15-hydroxy-7-oxodehydroabietic acid.

MiR408-SmLAC3 Module Participates in Salvianolic Acid B Synthesis in Salvia miltiorrhiza

MicroRNAs (miRNAs) are important regulators of gene expression involved in plant development and abiotic stress responses. Recently, miRNAs have also been reported to be engaged in the regulation of secondary plant metabolism. However, there are few functional studies of miRNAs in medicinal plants. For this study, we obtained Sm-miR408 interference lines to investigate the function of Sm-miR408 in a medicinal model plant (Salvia miltiorrhiza). It was found that inhibiting the expression of Sm-miR408 could increase the content of salvianolic acid B and rosmarinic acid in the roots. The SmLAC3 and Sm-miR408 expression patterns were analyzed by qRT-PCR. A 5’ RLM-RACE assay confirmed that Sm-miR408 targets and negatively regulates SmLAC3. Moreover, the overexpression of SmLAC3 in S. miltiorrhiza promoted the accumulation of salvianolic acids in the roots. Furthermore, the lignin content of the roots in overexpressed SmLAC3 lines was decreased. Taken together, these findings indicated that Sm-miR408 modulates the accumulation of phenolic acids in S. miltiorrhiza by targeting SmLAC3 expression levels.

Multiplexed CRISPR/Cas9-Mediated Knockout of Laccase Genes in Salvia miltiorrhiza Revealed Their Roles in Growth, Development, and Metabolism

Laccases are multicopper-containing glycoproteins related to monolignol oxidation and polymerization. These properties indicate that laccases may be involved in the formation of important medicinal phenolic acid compounds in Salvia miltiorrhiza such as salvianolic acid B (SAB), which is used for cardiovascular disease treatment. To date, 29 laccases have been found in S. miltiorrhiza (SmLACs), and some of which (SmLAC7 and SmLAC20) have been reported to influence the synthesis of phenolic acids. Because of the functional redundancy of laccase genes, their roles in S. miltiorrhiza are poorly understood. In this study, the CRISPR/Cas9 system was used for targeting conserved domains to knockout multiple genes of laccase family in S. miltiorrhiza. The expressions of target laccase genes as well as the phenolic acid biosynthesis key genes decrease dramatically in editing lines. Additionally, the growth and development of hairy roots was significantly retarded in the gene-edited lines. The cross-sections examination of laccase mutant hairy roots showed that the root development was abnormal and the xylem cells in the edited lines became larger and looser than those in the wild type. Additionally, the accumulation of RA as well as SAB was decreased, and the lignin content was nearly undetectable. It suggested that SmLACs play key roles in development and lignin formation in the root of S. miltiorrhiza and they are necessary for phenolic acids biosynthesis.

Genome-wide identification and expression analysis of the superoxide dismutase (SOD) gene family in Salvia miltiorrhiza

Adverse environmental conditions, such as salinity, cold, drought, heavy metals, and pathogens affect the yield and quality of Salvia miltiorrhiza, a well-known medicinal plant used for the treatment of cardiovascular and cerebrovascular diseases. Superoxide dismutase (SOD), a key enzyme of antioxidant system in plants, plays a vital role in protecting plants against various biotic and abiotic stresses via scavenging the reactive oxygen species produced by organisms. However, little is known about the SOD gene family in S. miltiorrhiza. In this study, eight SOD genes, including three Cu/Zn-SODs, two Fe-SODs and three Mn-SODs, were identified in the S. miltiorrhiza genome. Their gene structures, promoters, protein features, phylogenetic relationships, and expression profiles were comprehensively investigated. Gene structure analysis implied that most SmSODs have different introns/exons distrbution patterns. Many cis-elements related to different stress responses or plant hormones were found in the promoter of each SmSOD. Expression profile analysis indicated that SmSODs exhibited diverse responses to cold, salt, drought, heavy metal, and plant hormones. Additionally, 31 types of TFs regulating SmSODs were predicted and analyzed. These findings provided valuable information for further researches on the functions and applications of SmSODs in S. miltiorrhiza growth and adaptation to stress.

Genome-wide characterization of the laccase gene family in Setaria viridis reveals members potentially involved in lignification

Five laccase genes are potentially involved in developmental lignification in the model C4 grass Setaria viridis and their different tissue specificities suggest subfunctionalization events. Plant laccases are copper-containing glycoproteins involved in monolignol oxidation and, therefore, their activity is essential for lignin polymerization. Although these enzymes belong to large multigene families with highly redundant members, not all of them are thought to be involved in lignin metabolism. Here, we report on the genome-wide characterization of the laccase gene family in the model C4 grass Setaria viridis and further identification of the members potentially involved in monolignol oxidation. A total of 52 genes encoding laccases (SvLAC1 to SvLAC52) were found in the genome of S. viridis, and phylogenetic analyses showed that these genes were heterogeneously distributed among the characteristic six subclades of the family and are under relaxed selective constraints. The observed expansion in the total number of genes in this species was mainly caused by tandem duplications within subclade V, which accounts for 68% of the whole family. Comparative phylogenetic analyses showed that the expansion of subclade V is specifically observed for the Paniceae tribe within the Panicoideae subfamily in grasses. Five SvLAC genes (SvLAC9, SvLAC13, SvLAC15, SvLAC50, and SvLAC52) fulfilled the criteria established to identify lignin-related candidates: (1) phylogenetic proximity to previously characterized lignin-related laccases from other species, (2) similar expression pattern to that observed for lignin biosynthetic genes in the S. viridis elongating internode, and (3) high expression in S. viridis tissues undergoing active lignification. In addition, in situ hybridization experiments not only confirmed that these selected SvLAC genes were expressed in lignifying cells, but also that their expression showed different tissue specificities, suggesting subfunctionalization events within the family. These five laccase genes are strong candidates to be involved in lignin polymerization in S. viridis and might be good targets for lignin bioengineering strategies.