Combining metabolomics and transcriptomics to characterize tanshinone biosynthesis in Salvia miltiorrhiza

Genome-wide analysis, molecular cloning and expression profiling reveal tissue-specifically expressed, feedback-regulated, stress-responsive and alternatively spliced novel genes involved in gibberellin metabolism in Salvia miltiorrhiza

BACKGROUND

Gibberellin (GA), a classical phytohormone, plays significant roles in plant growth and development. It shares the important intermediate diphosphate precursor, GGPP, with the main lipophilic bioactive components, diterpenoid tanshinones in Salvia miltiorrhiza Bunge, one of the most important Traditional Chinese Medicine materials and an emerging model medicinal plant. Analysis of GA metabolism and regulation may help to demonstrate the biological functions of GAs and the crosstalk between GA metabolism and tanshinone biosynthesis in S. miltiorrhiza. However, genes involved in the conversion of ent-kaurene to GAs have not been systematically studied.

RESULTS

Through genome-wide prediction and molecular cloning, twenty two candidate gibberellin metabolism pathway genes were systematically identified for the first time. It includes a SmKO, two SmKAOs, six SmGA20oxs, two SmGA3oxs and eleven SmGA2oxs, of which twenty genes are novel. The deduced proteins showed sequence conservation and divergence. Gibberellin metabolism pathway genes exhibited tissue-specific expression patterns and responded differentially to exogenous GA3 treatment, indicating differential regulation of gibberellin metabolism in different tissue types in S. miltiorrhiza. SmKAO1, SmKAO2, SmGA2ox2, and SmGA2ox4-SmGA2ox7 were significantly up-regulated; SmGA20ox2, SmGA3ox1, SmGA2ox1, SmGA2ox8, SmGA2ox10 and SmGA2ox11 were significantly down-regulated; while the responses of many other genes varied among different tissue-types and time-points of GA3 treatment, suggesting the complexity of feedback regulation. Tanshinone biosynthesis-related genes, such as SmCPS1 and SmKSL1, were up-regulated in response to GA3 treatment. Among the 22 identified genes, nine responded to yeast extract and Ag(+)-treatment in S. miltiorrhiza hairy roots. Moreover, tissue-specifically expressed splice variants were identified for SmKO, SmGA20ox3, SmGA2ox3 and SmGA2ox11, of which SmKOv1, SmGA20ox3v and SmGA2ox11v1 were GA3-responsive, suggesting the importance of alternative splicing in regulating GA metabolism.

CONCLUSIONS

The results show tissue-specifically expressed, feedback-regulated, stress-responsive and alternatively spliced novel genes and reveal multiple layer regulation of GA metabolism and crosstalk between gibberellin metabolism and tanshinone biosynthesis in S. miltiorrhiza.

Enhancing diterpenoid concentration in Salvia miltiorrhiza hairy roots through pathway engineering with maize C1 transcription factor

Tanshinones are valuable natural diterpenoids from danshen (Salvia miltiorrhiza Bunge). Here, it was demonstrated that maize transcription factor C1 improved the accumulation of tanshinones by comprehensively upregulating the pathway genes, especially SmMDC and SmPMK in danshen hairy roots, yielding total tanshinones up to 3.59mg g(-1) of dry weight in line C1-6, a 3.4-fold increase compared with the control. Investigation of 2024bp of the SmMDC promoter fragment revealed that C1-mediated upregulation of terpenoid genes was possibly due to the direct interaction of C1 with its recognition sequences. The increase of tanshinones was accompanied by a decrease of salvianolic acid production, the other bioactive ingredient in danshen, by up to 37% compared with the control. This was the result of the downregulation of SmTAT, the entry-point gene of the tyrosine pathway, which promoted metabolic flow to anthocyanins rather than to salvianolic acids. Based on the findings of the present study, it was concluded that cis-acting elements shared by terpenoid and phenylpropanoid biosynthetic genes are partially responsible for the C1-stimulated variation of tanshinone and salvianolic acid concentrations.

Analysis of Polygala tenuifolia Transcriptome and Description of Secondary Metabolite Biosynthetic Pathways by Illumina Sequencing

Radix polygalae, the dried roots of Polygala tenuifolia and P. sibirica, is one of the most well-known traditional Chinese medicinal plants. Radix polygalae contains various saponins, xanthones, and oligosaccharide esters and these compounds are responsible for several pharmacological properties. To provide basic breeding information, enhance molecular biological analysis, and determine secondary metabolite biosynthetic pathways of P. tenuifolia, we applied Illumina sequencing technology and de novo assembly. We also applied this technique to gain an overview of P. tenuifolia transcriptome from samples with different years. Using Illumina sequencing, approximately 67.2% of unique sequences were annotated by basic local alignment search tool similarity searches against public sequence databases. We classified the annotated unigenes by using Nr, Nt, GO, COG, and KEGG databases compared with NCBI. We also obtained many candidates CYP450s and UGTs by the analysis of genes in the secondary metabolite biosynthetic pathways, including putative terpenoid backbone and phenylpropanoid biosynthesis pathway. With this transcriptome sequencing, future genetic and genomics studies related to the molecular mechanisms associated with the chemical composition of P. tenuifolia may be improved. Genes involved in the enrichment of secondary metabolite biosynthesis-related pathways could enhance the potential applications of P. tenuifolia in pharmaceutical industries.

Combination of transcriptomic and metabolomic analyses reveals a JAZ repressor in the jasmonate signaling pathway of Salvia miltiorrhiza

Jasmonates (JAs) are plant-specific key signaling molecules that respond to various stimuli and are involved in the synthesis of secondary metabolites. However, little is known about the JA signal pathway, especially in economically significant medicinal plants. To determine the functions of novel genes that participate in the JA-mediated accumulation of secondary metabolites, we examined the metabolomic and transcriptomic signatures from Salvia miltiorrhiza. For the metabolome, 35 representative metabolites showing significant changes in rates of accumulation were extracted and identified. We also screened out 2131 differentially expressed unigenes, of which 30 were involeved in the phenolic secondary metabolic pathway, while 25 were in the JA biosynthesis and signal pathways. Among several MeJA-induced novel genes, SmJAZ8 was selected for detailed functional analysis. Transgenic plants over-expressing SmJAZ8 exhibited a JA-insensitive phenotype, suggesting that the gene is a transcriptional regulator in the JA signal pathway of S. miltiorrhiza. Furthermore, this transgenic tool revealed that JAZ genes have novel function in the constitutive accumulation of secondary metabolites. Based on these findings, we propose that the combined strategy of transcriptomic and metabolomic analyses is valuable for efficient discovery of novel genes in plants.

Enzymes for synthetic biology of ambroxide-related diterpenoid fragrance compounds

Ambrox and related ambroxides are highly priced in the fragrance industry, and valued for their delicate odor and fixative properties. Historically, ambrox was obtained from ambergris, a waxy excretion produced by sperm whales, now an endangered species. Synthetic ambroxides have replaced ambergris in perfume manufacture. Plant labdane diterpenoids can serve as starting material for ambroxide synthesis. Among these, the diterpene alcohol sclareol is the major industrial precursor obtained from cultivated clary sage (Salvia sclarea). In plants, a large family of diterpene synthase (diTPS) enzymes controls key reactions in diterpenoid biosynthesis. Advanced metabolite profiling and high-throughput sequencing of fragrant and medicinal plants have accelerated discovery of novel diTPS functions, providing a resource for combinatorial synthetic biology and metabolic engineering approaches. This chapter highlights recent progress on the discovery, characterization, and engineering of plant diTPSs with potential uses in ambroxide production. It features biosynthesis of sclareol, cis-abienol, and diterpene resin acids, as sources of genes and enzymes for diterpenoid bioproducts.

The Biosynthetic Pathways of Tanshinones and Phenolic Acids in Salvia miltiorrhiza

Secondary metabolites from plants play key roles in human medicine and chemical industries. Due to limited accumulation of secondary metabolites in plants and their important roles, characterization of key enzymes involved in biosynthetic pathway will enable metabolic engineering or synthetic biology to improve or produce the compounds in plants or microorganisms, which provides an alternative for production of these valuable compounds. Salvia miltiorrhiza, containing tanshinones and phenolic acids as its active compounds, has been widely used for the treatment of cardiovascular and cerebrovascular diseases. The biosynthetic analysis of secondary metabolites in S. miltiorrhiza has made great progress due to the successful genetic transformation system, simplified hairy roots system, and high-throughput sequencing. The cloned genes in S. miltiorrhiza had provided references for functional characterization of the post-modification steps involved in biosynthesis of tanshinones and phenolic acids, and further utilization of these steps in metabolic engineering. The strategies used in these studies could provide solid foundation for elucidation of biosynthetic pathways of diterpenoids and phenolic acids in other species. The present review systematically summarizes recent advances in biosynthetic pathway analysis of tanshinones and phenolic acids as well as synthetic biology and metabolic engineering applications of the rate-limiting genes involved in the secondary metabolism in S. miltiorrhiza.

Genome-wide characterisation and analysis of bHLH transcription factors related to tanshinone biosynthesis in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge (Labiatae) is an emerging model plant for traditional medicine, and tanshinones are among the pharmacologically active constituents of this plant. Although extensive chemical and pharmaceutical studies of these compounds have been performed, studies on the basic helix-loop-helix (bHLH) transcription factors that regulate tanshinone biosynthesis are limited. In our study, 127 bHLH transcription factor genes were identified in the genome of S. miltiorrhiza, and phylogenetic analysis indicated that these SmbHLHs could be classified into 25 subfamilies. A total of 19 sequencing libraries were constructed for expression pattern analyses using RNA-Seq. Based on gene-specific expression patterns and up-regulated expression patterns in response to MeJA treatment, 7 bHLH genes were revealed as potentially involved in the regulation of tanshinone biosynthesis. Among them, the gene expression of SmbHLH37, SmbHLH74 and SmbHLH92 perfectly matches the accumulation pattern of tanshinone biosynthesis in S. miltiorrhiza. Our results provide a foundation for understanding the molecular basis and regulatory mechanisms of bHLH transcription factors in S. miltiorrhiza.

Identification and characterization of mRNA-like noncoding RNAs in Salvia miltiorrhiza

Identification and characterization of 5,446 mlncRNAs from Salvia miltiorrhiza showed that the majority of identified mlncRNAs were stress responsive, providing a framework for elucidating mlncRNA functions in S. miltiorrhiza. mRNA-like noncoding RNAs (mlncRNAs) are transcribed by RNA polymerase II and are polyadenylated, capped and spliced. They play important roles in plant development and defense responses. However, there is no information available for mlncRNAs in Salvia miltiorrhiza Bunge, the first Chinese medicinal material entering the international market. To perform a transcriptome-wide identification of S. miltiorrhiza mlncRNAs, we assembled over 8 million RNA-seq reads from GenBank database and 5,624 ESTs from PlantGDB into 44422 unigenes. Using a computational identification pipeline, we identified 5446 S. miltiorrhiza mlncRNA candidates from the assembled unigenes. Of the 5446 mlncRNAs, 2 are primary transcripts of conserved miRNAs, and 2030 can be grouped into 470 families with at least two members in a family. Quantitative real-time PCR analysis of mlncRNAs with at least 900 nt showed that the majority were differentially expressed in roots, stems, leaves and flowers and responsive to methyl jasmonate (MeJA) treatment in S. miltiorrhiza. Analysis of published RNA-seq data showed that a total of 3,044 mlncRNAs were expressed in hairy roots of S. miltiorrhiza and the expression of 1,904 of the 3,044 mlncRNAs was altered by yeast extract and Ag(+) treatment. The results indicate that the majority of mlncRNAs are involved in plant response to stress. This study provides a framework for understanding the roles of mlncRNAs in S. miltiorrhiza.

Metabolomics for phytochemical discovery: development of statistical approaches using a cranberry model system

Metabolomics is the qualitative and quantitative analysis of all of the small molecules in a biological sample at a specific time and influence. Technologies for metabolomics analysis have developed rapidly as new analytical tools for chemical separations, mass spectrometry, and NMR spectroscopy have emerged. Plants have one of the largest metabolomes, and it is estimated that the average plant leaf can contain upward of 30 000 phytochemicals. In the past decade, over 1200 papers on plant metabolomics have been published. A standard metabolomics data set contains vast amounts of information and can either investigate or generate hypotheses. The key factors in using plant metabolomics data most effectively are the experimental design, authentic standard availability, extract standardization, and statistical analysis. Using cranberry (Vaccinium macrocarpon) as a model system, this review will discuss and demonstrate strategies and tools for analysis and interpretation of metabolomics data sets including eliminating false discoveries and determining significance, metabolite clustering, and logical algorithms for discovery of new metabolites and pathways. Together these metabolomics tools represent an entirely new pipeline for phytochemical discovery.

Molecular cloning and expression analysis of WRKY transcription factor genes in Salvia miltiorrhiza

Background

WRKY proteins comprise a large family of transcription factors and play important regulatory roles in plant development and defense response. The WRKY gene family in Salvia miltiorrhiza has not been characterized.

Results

A total of 61 SmWRKYs were cloned from S. miltiorrhiza. Multiple sequence alignment showed that SmWRKYs could be classified into 3 groups and 8 subgroups. Sequence features, the WRKY domain and other motifs of SmWRKYs are largely conserved with Arabidopsis AtWRKYs. Each group of WRKY domains contains characteristic conserved sequences, and group-specific motifs might attribute to functional divergence of WRKYs. A total of 17 pairs of orthologous SmWRKY and AtWRKY genes and 21 pairs of paralogous SmWRKY genes were identified. Maximum likelihood analysis showed that SmWRKYs had undergone strong selective pressure for adaptive evolution. Functional divergence analysis suggested that the SmWRKY subgroup genes and many paralogous SmWRKY gene pairs were divergent in functions. Various critical amino acids contributed to functional divergence among subgroups were detected. Of the 61 SmWRKYs, 22, 13, 4 and 1 were predominantly expressed in roots, stems, leaves, and flowers, respectively. The other 21 were mainly expressed in at least two tissues analyzed. In S. miltiorrhiza roots treated with MeJA, significant changes of gene expression were observed for 49 SmWRKYs, of which 26 were up-regulated, 18 were down-regulated, while the other 5 were either up-regulated or down-regulated at different time-points of treatment. Analysis of published RNA-seq data showed that 42 of the 61 identified SmWRKYs were yeast extract and Ag⁺-responsive. Through a systematic analysis, SmWRKYs potentially involved in tanshinone biosynthesis were predicted.

Conclusion

These results provide insights into functional conservation and diversification of SmWRKYs and are useful information for further elucidating SmWRKY functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1411-x) contains supplementary material, which is available to authorized users.

Computational Identification and Systematic Classification of Novel Cytochrome P450 Genes in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most economically important medicinal plants. Cytochrome P450 (CYP450) genes have been implicated in the biosynthesis of its active components. However, only a dozen full-length CYP450 genes have been described, and there is no systematic classification of CYP450 genes in S. miltiorrhiza. We obtained 77,549 unigenes from three tissue types of S. miltiorrhiza using RNA-Seq technology. Combining our data with previously identified CYP450 sequences and scanning with the CYP450 model from Pfam resulted in the identification of 116 full-length and 135 partial-length CYP450 genes. The 116 genes were classified into 9 clans and 38 families using standard criteria. The RNA-Seq results showed that 35 CYP450 genes were co-expressed with CYP76AH1, a marker gene for tanshinone biosynthesis, using r≥0.9 as a cutoff. The expression profiles for 16 of 19 randomly selected CYP450 obtained from RNA-Seq were validated by qRT-PCR. Comparing against the KEGG database, 10 CYP450 genes were found to be associated with diterpenoid biosynthesis. Considering all the evidence, 3 CYP450 genes were identified to be potentially involved in terpenoid biosynthesis. Moreover, we found that 15 CYP450 genes were possibly regulated by antisense transcripts (r≥0.9 or r≤–0.9). Lastly, a web resource (SMCYP450, http://www.herbalgenomics.org/samicyp450) was set up, which allows users to browse, search, retrieve and compare CYP450 genes and can serve as a centralized resource.

Biochemical characterization of the castor bean ent-kaurene synthase(-like) family supports quantum chemical view of diterpene cyclization

It has become apparent that plants have extensively diversified their arsenal of labdane-related diterpenoids (LRDs), in part via gene duplication and neo-functionalization of the ancestral ent-kaurene synthase (KS) required for gibberellin metabolism. For example, castor bean (Ricinus communis) was previously shown to produce an interesting set of biosynthetically related diterpenes, specifically ent-sandracopimaradiene, ent-beyerene, and ent-trachylobane, in addition to ent-kaurene, using four separate diterpene synthases, albeit these remain unidentified. Notably, despite mechanistic similarity of the underlying reaction to that catalyzed by KSs, ent-beyerene and ent-trachylobane synthases have not yet been identified. Given our interest in LRD biosynthesis, and the recent availability of the castor bean genome sequence, a synthetic biology approach was applied to biochemically characterize the four KS(-like) enzymes [KS(L)s] found in Ricinus communis [i.e., the RcKS(L)s]. In particular, using bacteria engineered to produce the relevant ent-copalyl diphosphate precursor and synthetic genes based on the predicted RcKS(L)s, although this ultimately required correction of a "splicing" error in one of the predicted genes, highlighting the dependence of such a synthetic biology approach on accurate gene sequences. Nevertheless, it is possible to assign each of the four RcKS(L)s to one of the previously observed diterpene synthase activities, providing access to functionally enzymes. Intriguingly, the product distribution of the RcKS(L)s seems to support the distinct diterpene synthase reaction mechanism proposed by quantum chemical calculations, rather than the classically proposed pathway.