De novo transcriptome sequencing in Salvia miltiorrhiza to identify genes involved in the biosynthesis of active ingredients

CDCS-TCM: A framework based on complex network theory to analyze the causality and dynamic correlation of substances in the metabolic process of traditional Chinese medicine

ETHNOPHARMACOLOGICAL RELEVANCE

Traditional Chinese medicine, with the feature of synergistic effects of multi-component, multi-pathway and multi-target, plays an important role in the treatment of cancer, cardiovascular and cerebrovascular diseases, etc. However, chemical components in traditional Chinese medicine are complex and most of the pharmacological mechanisms remain unclear, especially the relationships of chemical components change during the metabolic process.

AIM OF STUDY

Our aim is to provide a method based on complex network theory to analyze the causality and dynamic correlation of substances in the metabolic process of traditional Chinese medicine.

MATERIALS AND METHODS

We proposed a framework named CDCS-TCM to analyze the causality and dynamic correlation between substances in the metabolic process of traditional Chinese medicine. Our method mainly consists two parts. The first part is to discover the local and global causality by the causality network. The second part is to investigate the dynamic correlations and identify the essential substance by dynamic substance correlation network.

RESULTS

We developed a CDCS-TCM method to analyze the causality and dynamic correlation of substances. Using the XiangDan Injection for ischemic stroke as an example, we have identified the important substances in the metabolic process including substance pairs with strong causality and the dynamic changes of the core effector substance clusters.

CONCLUSION

The proposed framework will be useful for exploring the correlations of active ingredients in traditional Chinese medicine more effectively and will provide a new perspective for the elucidation of drug action mechanisms and the new drug discovery.

Transcriptome sequencing and metabolome analysis reveal the molecular mechanism of Salvia miltiorrhiza in response to drought stress

Salvia miltiorrhiza is commonly used as a Chinese herbal medicine to treat different cardiovascular and cerebrovascular illnesses due to its active ingredients. Environmental conditions, especially drought stress, can affect the yield and quality of S. miltiorrhiza. However, moderate drought stress could improve the quality of S. miltiorrhiza without significantly reducing the yield, and the mechanism of this initial drought resistance is still unclear. In our study, transcriptome and metabolome analyses of S. miltiorrhiza under different drought treatment groups (CK, A, B, and C groups) were conducted to reveal the basis for its drought tolerance. We discovered that the leaves of S. miltiorrhiza under different drought treatment groups had no obvious shrinkage, and the malondialdehyde (MDA) contents as well as superoxide dismutase (SOD) and peroxidase (POD) activities dramatically increased, indicating that our drought treatment methods were moderate, and the leaves of S. miltiorrhiza began to initiate drought resistance. The morphology of root tissue had no significant change under different drought treatment groups, and the contents of four tanshinones significantly enhanced. In all, 5213, 6611, and 5241 differentially expressed genes (DEGs) were shared in the A, B, and C groups compared with the CK group, respectively. The results of KEGG and co-expression analysis showed that the DEGs involved in plant-pathogen interactions, the MAPK signaling pathway, phenylpropanoid biosynthesis, flavonoid biosynthesis, and plant hormone signal transduction responded to drought stress and were strongly correlated with tanshinone biosynthesis. Furthermore, the results of metabolism analysis indicated that 67, 72, and 92 differentially accumulated metabolites (DAMs), including fumarate, ferulic acid, xanthohumol, and phytocassanes, which were primarily involved in phenylpropanoid biosynthesis, flavonoid biosynthesis, and diterpenoid biosynthesis pathways, were detected in these groups. These discoveries provide valuable information on the molecular mechanisms by which S. miltiorrhiza responds to drought stress and will facilitate the development of drought-resistant and high-quality S. miltiorrhiza production.

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.

Genome-wide identification and expression analysis of the GASA gene family in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

The Gibberellic Acid-Stimulated Arabidopsis (GASA) gene family is widely involved in the regulation of plant growth, development, and stress response. However, information on the GASA gene family has not been reported in Chinese cabbage (Brassica rapa L. ssp. pekinensis).

RESULTS

Here, we conducted genome-wide identification and analysis of the GASA genes in Chinese cabbage. In total, 15 GASA genes were identified in the Chinese cabbage genome, and the physicochemical property, subcellular location, and tertiary structure of the corresponding GASA proteins were elucidated. Phylogenetic analysis, conserved motif, and gene structure showed that the GASA proteins were divided into three well-conserved subfamilies. Synteny analysis proposed that the expansion of the GASA genes was influenced mainly by whole-genome duplication (WGD) and transposed duplication (TRD) and that duplication gene pairs were under negative selection. Cis-acting elements of the GASA promoters were involved in plant development, hormonal and stress responses. Expression profile analysis showed that the GASA genes were widely expressed in different tissues of Chinese cabbage, but their expression patterns appeared to diverse. The qRT-PCR analysis of nine GASA genes confirmed that they responded to salt stress, heat stress, and hormonal triggers.

CONCLUSIONS

Overall, this study provides a theoretical basis for further exploring the important role of the GASA gene family in the functional genome of Chinese cabbage.

Transcriptomic and proteomic analysis reveals (E)-2-hexenal modulates tomato resistance against Botrytis cinerea by regulating plant defense mechanism

In a previous study, we observed that (E)-2-hexenal stimulated systemic resistance against B. cinerea in tomato plants. However, the molecular mechanisms underlying (E)-2-hexenal-mediated regulation of systemic immunity against B. cinerea remained unclear. In the current study, the global mechanism underlying (E)-2-hexenal-meidated regulation of biotic stress tolerance in tomato was investigated using RNA-seq- and LC-MS/MS- integrated transcriptomic and proteomic analyses. Compared to control plants, (E)-2-hexenal-treated plants exhibited reduced susceptibility to B. cinerea, with a 50.51% decrease in lesion diameters. Meanwhile, (E)-2-hexenal vapor fumigation significantly increased total phenolic content and activities of various antioxidant enzymes peroxidase (POD), phenylalanine ammonia lyase (PAL), and lipoxygenase (LOX). A total of 233 differentially expressed genes (DEGs) and 400 differentially expressed proteins (DEPs), respectively, were identified. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that (E)-2-hexenal treatment markedly affected the expression of genes involved in multiple metabolic pathways, especially glutathione metabolism, phenylpropanoid biosynthesis, plant hormone signal transduction, and MAPK signaling pathway. Notably, proteomic analysis revealed modulation of the activities of several defense response proteins, such as pathogenesis-related (PR) proteins (Solyc02g031950.3.1, Solyc02g031920.4.1, and Solyc04g064870.3.1), peroxidases (Solyc06g050440.3.1, Solyc01g105070.3.1, Solyc01g015080.3.1, Solyc03g025380.3.1 and Solyc06g076630.3.1). Our results provide a comprehensive analysis of the effects of (E)-2-hexenal treatment on the transcriptome and proteome of tomato plants, which might be used as a reference in further studies on plant defense responses against pathogens.

Integrative Analysis of Transcriptome, Proteome, and Phosphoproteome Reveals Potential Roles of Photosynthesis Antenna Proteins in Response to Brassinosteroids Signaling in Maize

Brassinosteroids are a recently discovered group of substances that promote plant growth and productivity. Photosynthesis, which is vital for plant growth and high productivity, is strongly influenced by brassinosteroid signaling. However, the molecular mechanism underlying the photosynthetic response to brassinosteroid signaling in maize remains obscure. Here, we performed integrated transcriptome, proteome, and phosphoproteomic analyses to identify the key photosynthesis pathway that responds to brassinosteroid signaling. Transcriptome analysis suggested that photosynthesis antenna proteins and carotenoid biosynthesis, plant hormone signal transduction, and MAPK signaling in CK VS EBR and CK VS Brz were significantly enriched in the list of differentially expressed genes upon brassinosteroids treatment. Consistently, proteome and phosphoproteomic analyses indicated that photosynthesis antenna and photosynthesis proteins were significantly enriched in the list of differentially expressed proteins. Thus, transcriptome, proteome, and phosphoproteome analyses showed that major genes and proteins related to photosynthesis antenna proteins were upregulated by brassinosteroids treatment in a dose-dependent manner. Meanwhile, 42 and 186 transcription factor (TF) responses to brassinosteroid signals in maize leaves were identified in the CK VS EBR and CK VS Brz groups, respectively. Our study provides valuable information for a better understanding of the molecular mechanism underlying the photosynthetic response to brassinosteroid signaling in maize.

Transcriptome sequencing and DEG analysis in different developmental stages of floral buds induced by potassium chlorate in Dimocarpus longan

Potassium chlorate can promote off-season flowering in longan, but the molecular mechanisms are poorly understood. In this study, four-year-old 'Shixia' longan trees were injected in the trunk with potassium chlorate, and terminal buds were sampled and analyzed using transcriptomics and bioinformatics tools. To generate a reference longan transcriptome, we obtained 207,734 paired-end reads covering a total of 58,514,149 bp, which we assembled into 114,445 unigenes. Using this resource, we identified 3,265 differentially expressed genes (DEGs) that were regulated in longan terminal buds in response to potassium chlorate treatment for 2, 6 or 30 days, including 179 transcription factor genes. By reference to the Arabidopsis literature, we then defined 38 longan genes involved in flowering, from which we constructed the longan flowering pathway. According to RNA-seq data, at least 24 of these genes, which participate in multiple signaling pathways, are involved in potassium chlorate-stimulated floral induction, and the differential regulation in terminal buds of ten floral pathway genes (GI, CO, GID1, GA4, GA5, FLC, AP1, LFY, FT and SOC1) was confirmed by qRT-PCR. These data will contribute to an improved understanding of the functions of key genes involved in longan floral induction by potassium chlorate.

Integrative transcriptomic and metabolomic analyses unveil tanshinone biosynthesis in Salvia miltiorrhiza root under N starvation stress

Salvia miltiorrhiza is a model plant for Chinese herbal medicine with significant pharmacologic effects due to its tanshinone components. Our previous study indicated that nitrogen starvation stress increased its tanshinone content. However, the molecular mechanism of this low nitrogen-induced tanshinone biosynthesis is still unclear. Thus, this study aimed to elucidate the molecular mechanism of tanshinone biosynthesis in S. miltiorrhiza under different N conditions [N-free (N0), low-N (Nl), and full-N (Nf, as control) conditions] by using transcriptome and metabolome analyses. Our results showed 3,437 and 2,274 differentially expressed unigenes between N0 and Nf as well as Nl and Nf root samples, respectively. N starvation (N0 and Nl) promoted the expression of the genes involved in the MVA and MEP pathway of tanshinone and terpenoid backbone biosynthesis. Gene ontology and KEGG analyses revealed that terpenoid backbone biosynthesis, hormone signal transduction, and phenylpropanoid biosynthesis were promoted under N starvation conditions, whereas starch and sucrose metabolisms, nitrogen and phosphorus metabolisms, as well as membrane development were inhibited. Furthermore, metabolome analysis showed that metabolite compounds and biosynthesis of secondary metabolites were upregulated. This study provided a novel insight into the molecular mechanisms of tanshinone production in S. miltiorrhiza in response to nitrogen stress.

Molecular characterization of a Novel NAD+-dependent farnesol dehydrogenase SoFLDH gene involved in sesquiterpenoid synthases from Salvia officinalis

Salvia officinalis is one of the most important medicinal and aromatic plants in terms of nutritional and medicinal value because it contains a variety of vital active ingredients. Terpenoid compounds, particularly monoterpenes (C10) and sesquiterpenes, are the most important and abundant among these active substances (C15). Terpenes play a variety of roles and have beneficial biological properties in plants. With these considerations, the current study sought to clone theNAD+-dependent farnesol dehydrogenase (SoFLDH, EC: 1.1.1.354) gene from S. officinalis. Functional analysis revealed that, SoFLDH has an open reading frame of 2,580 base pairs that encodes 860 amino acids.SoFLDH has two conserved domains and four types of highly conserved motifs: YxxxK, RXR, RR (X8) W, TGxxGhaG. However, SoFLDH was cloned from Salvia officinalis leaves and functionally overexpressed in Arabidopsis thaliana to investigate its role in sesquiterpenoid synthases. In comparison to the transgenic plants, the wild-type plants showed a slight delay in growth and flowering formation. To this end, a gas chromatography-mass spectrometry analysis revealed that SoFLDH transgenic plants were responsible for numerous forms of terpene synthesis, particularly sesquiterpene. These results provide a base for further investigation on SoFLDH gene role and elucidating the regulatory mechanisms for sesquiterpene synthesis in S. offcinalis. And our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. offcinalis.

Identification of differentially expressed genes at different post-natal development stages of longissimus dorsi muscle in Tianzhu white yak

The regulatory mechanisms controlling post-natal muscle development in the yak (Bos grunniens) are still largely unknown, yet the growth and development of muscle is a complex process that plays a crucial role in determining the yield and quality of an animal's meat. In this study, we performed a transcriptome analysis based on the RNA sequencing (RNA-Seq) of yak longissimus dorsi muscle tissue obtained from calves (6 months of age; 6 M), young adults (30 months of age; 30 M) and adult (54 months of age; 54 M) to identify which genes are differentially expressed and to investigate their temporal expression profiles. In total, 1788 differentially expressed genes (DEGs) (|log2FC| ≥ 1, P-adjusted < 0.05) were detected by pairwise comparisons between the different age groups. The expression levels of 10 of the DEGs were confirmed using reverse transcription-quantitative PCR (RT-qPCR), and the results were consistent with the transcriptome profile. A time-series expression profile analysis clustered the DEGs into four groups that could be divided into two classes (P < 0.05): class 1 profiles, which had up-regulated patterns of gene expression and class 2 profiles, which featured down-regulated patterns. Based on that cluster analysis, GO enrichment analysis revealed 1073, 127, and 184 terms as significantly enriched in biological process (BP), cellular component (CC), and molecular function (MF) categories in the class 1 profiles, while 714, 66, and 206 terms were significantly enriched in BP, CC, and MF in the class 2 profiles. A KEGG pathway analysis revealed that DEGs from the class 1 profiles were enriched in 62 pathways, with the most enriched being the phosphoinositide 3-kinase (PI3K) - protein kinase B (Akt)-signaling pathway. The DEGs from the class 2 profiles were enriched in 16 pathways, of which forkhead box protein O (FoxO) - signaling was the most enriched. Taken together, these results provide insight into the mechanisms of skeletal muscle development, as well suggesting some potential genes of importance for yak meat production.

Profiling and Functional Analysis of mRNAs during Skeletal Muscle Differentiation in Goats

Skeletal myogenesis is a complicated biological event that involves a succession of tightly controlled gene expressions. In order to identify novel regulators of this process, we performed mRNA-Seq studies of goat skeletal muscle satellite cells (MuSCs) cultured under proliferation (GM) and differentiation (DM1/DM5) conditions. A total of 19,871 goat genes were expressed during these stages, 198 of which represented novel transcripts. Notably, in pairwise comparisons at the different stages, 2551 differentially expressed genes (DEGs) were identified (p < 0.05), including 1560 in GM vs. DM1, 1597 in GM vs. DM5, and 959 in DM1 vs. DM5 DEGs. The time-series expression profile analysis clustered the DEGs into eight gene groups, three of which had significantly upregulated and downregulated patterns (p < 0.05). Functional enrichment analysis showed that DEGs were enriched for essential biological processes such as muscle structure development, muscle contraction, muscle cell development, striated muscle cell differentiation, and myofibril assembly, and were involved in pathways such as the MAPK, Wnt and PPAR signaling pathways. Moreover, the expression of eight DEGs (MYL2, DES, MYOG, FAP, PLK2, ADAM, WWC1, and PRDX1) was validated. These findings offer novel insights into the transcriptional regulation of skeletal myogenesis in goats.

De Novo Transcriptome Analysis Reveals Putative Genes Involved in Anthraquinone Biosynthesis in Rubia yunnanensis

Rubia yunnanensis Diels (R. yunnanensis), a Chinese perennial plant, is well-known for its medicinal values such as rheumatism, contusion, and anemia. It is rich in bioactive anthraquinones, but the biosynthetic pathways of anthraquinones in R. yunnanensis remain unknown. To investigate genes involved in anthraquinone biosynthesis in R. yunnanensis, we generated a de novo transcriptome of R. yunnanensis using the Illumina HiSeq 2500 sequencing platform. A total of 636,198 transcripts were obtained, in which 140,078 transcripts were successfully annotated. A differential gene expression analysis identified 15 putative genes involved in anthraquinone biosynthesis. Additionally, the hairy roots of R. yunnanensis were treated with 200 µM Methyl Jasmonate (MeJA). The contents of six bioactive anthraquinones and gene expression levels of 15 putative genes were measured using ultra performance liquid chromatography coupled with mass spectrometry (UPLC-MS/MS) and real-time quantitative polymerase chain reaction (RT-qPCR), respectively. The results showed that the expressions levels for 11 of the 15 genes and the contents of two of six anthraquinones significantly increased by MeJA treatment. Pearson’s correlation analyses indicated that the expressions of 4 of the 15 putative genes were positively correlated with the contents of rubiquinone (Q3) and rubiquinone-3-O-β-d-xylopranosyl-(1→6)-β-d-glucopyranoside (Q20). This study reported the first de novo transcriptome of R. yunnanensis and shed light on the anthraquinone biosynthesis and genetic information for R. yunnanensis.

SmbHLH60 and SmMYC2 antagonistically regulate phenolic acids and anthocyanins biosynthesis in Salvia miltiorrhiza

INTRODUCTION

Salvia miltiorrhiza is a renowned traditional Chinese medicinal plant with extremely high medicinal value, especially for cardiovascular and cerebrovascular diseases. The jasmonic acid (JA) signaling pathway plays an important role in the improved biosynthesis of secondary metabolites, which is mediated by a major transcriptional regulator, MYC2. However, the JA regulatory mechanism of secondary metabolites biosynthesis in S. miltiorrhiza is still largely unknown.

OBJECTIVES

Our work focuses on the dissection of the molecular mechanism of transcriptional regulation in MeJA-mediated biosynthesis of medicinal components of S. miltiorrhiza. We examined the role of MeJA-responsive bHLH transcription factors (TFs) in improving bioactive secondary metabolites accumulation in S. miltiorrhiza.

METHODS

Hairy root transformation based on CRISPR/Cas9 technique was used to decipher gene function(s). Changes in the content of phenolic acids were evaluated by HPLC. Y1H, EMSA and dual-LUC assays were employed to analyze the molecular mechanism of SmbHLH60 in the regulation on the biosynthesis of phenolic acids and anthocyanins. Y2H, BiFC and pull-down affinity assays were used to corroborate the interaction between SmbHLH60 and SmMYC2.

RESULTS

Being one of the most significantly negatively regulated bHLH genes by MeJA, a new transcription factor SmbHLH60 was discovered and characterized. Over-expression of SmbHLH60 resulted in significant inhibition of phenolic acid and anthocyanin biosynthesis in S. miltiorrhiza by transcriptionally repressing of target genes such as SmTAT1 and SmDFR, whereas CRISPR/Cas9-generated knockout of SmbHLH60 resulted in the opposite effect. In addition, SmbHLH60 and SmMYC2 formed a heterodimer to antagonistically regulate phenolic acid and anthocyanin biosynthesis.

CONCLUSION

Our results clarified that SmbHLH60 is a negativeregulator on the biosynthesis of phenolic acids and anthocyanins. SmbHLH60 competed with SmMYC2 in an antagonistic manner, providing new insights for the molecular mechanism of MeJA-mediated regulation on the biosynthesis of secondary metabolites in S. miltiorrhiza.

Insights into the plateau adaptation of Salvia castanea by comparative genomic and WGCNA analyses

INTRODUCTION

Salvia castanea, a wild plant species is adapted to extreme Qinghai-Tibetan plateau (QTP) environments. It is also used for medicinal purposes due to high ingredient of tanshinone IIA (T-IIA). Despite its importance to Chinese medicinal industry, the mechanisms associated with secondary metabolites accumulation (i.e. T-IIA and rosmarinic acid (RA)) in this species have not been characterized. Also, the role of special underground tissues in QTP adaptation of S. castanea is still unknown.

OBJECTIVES

We explored the phenomenon of periderm-like structure in underground stem center of S. castanea with an aim to unravel the molecular evolutionary mechanisms of QTP adaptation in this species.

METHODS

Morphologic observation and full-length transcriptome of S. castanea plants were conducted. Comparative genomic analyses of S. castanea with other 14 representative species were used to reveal its phylogenetic position and molecular evolutionary mechanisms. RNA-seq and WGCNA analyses were applied to understand the mechanisms of high accumulations of T-IIA and RA in S. castanea tissues.

RESULTS

Based on anatomical observations, we proposed a "trunk-branches" developmental model to explain periderm-like structure in the center of underground stem of S. castanea. Our study suggested that S. castanea branched off from cultivated Danshen around 16 million years ago. During the evolutionary process, significantly expanded orthologous gene groups, 24 species-specific and 64 positively selected genes contributed to morphogenesis and QTP adaptation in S. castanea. RNA-seq and WGCNA analyses unraveled underlying mechanisms of high accumulations of T-IIA and RA in S. castanea and identified NAC29 and TGA22 as key transcription factors.

CONCLUSION

We proposed a "trunk-branches" developmental model for the underground stem in S. castanea. Adaptations to extreme QTP environment in S. castanea are associated with accumulations of high secondary metabolites in this species.

Tissue-specific transcriptome for Rheum tanguticum reveals candidate genes related to the anthraquinones biosynthesis

Rheum tanguticum (Maxim. ex Regel) Maxim. ex Balf. is a herbaceous perennial plant indigenous to China, and its root and rhizomes were usually used as an important traditional Chinese medicine. However, the genomic resources are still scarce for R. tanguticum and even for Rheum genus. Transcriptome datasets from different tissues of R. tanguticum were obtained to screen the genes related to anthraquinones biosynthesis, and five free anthraquinones were also determined. Nine cDNA libraries of roots, stems and leaves were generated, and a total of 272 million high-quality reads were assembled into 257,942 unigenes. Based on the functional annotation, A total of 227 candidate enzyme genes involved in the MVA, MEP, shikimate and polyketide pathways were identified, and several differentially expressed genes found functionally associated with anthraquinones biosynthesis showed distinct tissue-specific expression patterns. Especially, we found that the expression levels of PKS III genes might result in the content differences of free anthraquinones in different tissues of R. tanguticum. Besides, 137,400 SSR loci were identified, and 64,081 SSR primer pairs were successfully designed based on these loci. Our results not only provide cues for the genetic mechanism of anthraquinone content differences in different tissues of R. tanguticum, but also lay genomic foundation for the subsequent genetic engineering and breeding for Rheum species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01099-8.

Transcriptome analysis reveals important candidate gene families related to oligosaccharides biosynthesis in Morindaofficinalis

Morinda officinalis How (MO) is one of the best-known traditional herbs and is widely cultivated in subtropical and tropical areas for many years, especially in southern China. Oligosaccharides are the major constituents in the roots of MO, which is well known for its therapeutic effects with anti-depression, anti-osteoporosis, memory-enhancing, ect. To date, the main gene families that regulate the biosynthetic pathway of MO oligosaccharides metabolism yet have been published. In our study, six cDNA libraries generated from six plants of MO were sequenced utilizing an Illumina HiSeq 4000 platform. Corresponding totals of more than 132.60 million clean reads were obtained from the six libraries and assembled into 25,812 unigenes with an average length of 1288 bp. Moreover, 6036 unigenes were found to be allocated to 26 pathways maps using several public databases, and 2538 differential expression genes (DEGs) were screened. Among them, 25 genes from three families were selected as the mainly candidate genes related to MO oligosaccharides biosynthesis. Then, the expression patterns of six DEGs closely related to MO oligosaccharides biosynthesis were verified by quantitative real-time PCR (qRT-PCR). Besides, the MO was clustered more closely to Coffea arabica of Rubiaceae. In summary, the transcriptomic analysis was used to investigate the differences in expression genes of oligosaccharides biosynthesis, with the notable outcome that several key gene families were closely linked to oligosaccharides biosynthesis.

Transcriptome and proteome of the corm, leaf and flower of Hypoxis hemerocallidea (African potato)

The corm of Hypoxis hemerocallidea, commonly known as the African potato, is used in traditional medicine to treat several medical conditions such as urinary infections, benign prostate hyperplasia, inflammatory conditions and testicular tumours. The metabolites contributing to the medicinal properties of H. hemerocallidea have been identified in several studies and, more recently, the active terpenoids of the plant were profiled. However, the biosynthetic pathways and the enzymes involved in the production of the terpene metabolites in H. hemerocallidea have not been characterised at a transcriptomic or proteomic level. In this study, total RNA extracted from the corm, leaf and flower tissues of H. hemerocallidea was sequenced on the Illumina HiSeq 2500 platform. A total of 143,549 transcripts were assembled de novo using Trinity and 107,131 transcripts were functionally annotated using the nr, GO, COG, KEGG and SWISS-PROT databases. Additionally, the proteome of the three tissues were sequenced using LC-MS/MS, revealing aspects of secondary metabolism and serving as data validation for the transcriptome. Functional annotation led to the identification of numerous terpene synthases such as nerolidol synthase, germacrene D synthase, and cycloartenol synthase amongst others. Annotations also revealed a transcript encoding the terpene synthase phytoalexin momilactone A synthase. Differential expression analysis using edgeR identified 946 transcripts differentially expressed between the three tissues and revealed that the leaf upregulates linalool synthase compared to the corm and the flower tissues. The transcriptome as well as the proteome of Hypoxis hemerocallidea presented here provide a foundation for future research.

Chia (Salvia hispanica) Gene Expression Atlas Elucidates Dynamic Spatio-Temporal Changes Associated With Plant Growth and Development

Chia (Salvia hispanica L.), now a popular superfood and a pseudocereal, is one of the richest sources of dietary nutrients such as protein, fiber, and polyunsaturated fatty acids (PUFAs). At present, the genomic and genetic information available in the public domain for this crop are scanty, which hinders an understanding of its growth and development and genetic improvement. We report an RNA-sequencing (RNA-Seq)-based comprehensive transcriptome atlas of Chia sampled from 13 tissue types covering vegetative and reproductive growth stages. We used ~355 million high-quality reads of total ~394 million raw reads from transcriptome sequencing to generate de novo reference transcriptome assembly and the tissue-specific transcript assemblies. After the quality assessment of the merged assemblies and implementing redundancy reduction methods, 82,663 reference transcripts were identified. About 65,587 of 82,663 transcripts were translated into 99,307 peptides, and we were successful in assigning InterPro annotations to 45,209 peptides and gene ontology (GO) terms to 32,638 peptides. The assembled transcriptome is estimated to have the complete sequence information for ~86% of the genes found in the Chia genome. Furthermore, the analysis of 53,200 differentially expressed transcripts (DETs) revealed their distinct expression patterns in Chia's vegetative and reproductive tissues; tissue-specific networks and developmental stage-specific networks of transcription factors (TFs); and the regulation of the expression of enzyme-coding genes associated with important metabolic pathways. In addition, we identified 2,411 simple sequence repeats (SSRs) as potential genetic markers from the transcripts. Overall, this study provides a comprehensive transcriptome atlas, and SSRs, contributing to building essential genomic resources to support basic research, genome annotation, functional genomics, and molecular breeding of Chia.

Crosstalk of Multi-Omics Platforms with Plants of Therapeutic Importance

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

Seasonal Variations of Rosmarinic Acid and Its Glucoside and Expression of Genes Related to Their Biosynthesis in Two Medicinal and Aromatic Species of Salvia subg. Perovskia

Simple Summary

Here, we studied two closely related medicinal and aromatic plants from Asia, called Russian sage or from their previously used Latin name–Perovskia. These plants contain various specialized metabolites called phenylpropanoids that contribute to their medicinal uses. In our experiments, several different specialized phytochemicals were traced down in the roots and leaves with the major metabolite called rosmarinic acid, known for health beneficial properties. In order to check if the composition of these plants is regulated by specific genes encoding proteins that assemble these phytochemicals, we analyzed their expression during the growth season (spring, summer and fall). Despite being the closest kin, the two species of Russian sage displayed different seasonal changes in the composition of bioactive metabolites and the activity of genes responsible for their production. The genes’ activity was correlated with rosmarinic acid content in the roots but not in the green parts of the plants. Two genes pointed out were linked to the regulation of rosmarinic acid biosynthesis, called RAS (for Rosmarinic Acid-Synthase) and a newly reported version of an oxidizing enzyme called Cyp98A14. These discoveries broaden our understanding of relationships between the genes’ activity and production of bioactive constituents in herbs such as the two studied species of Russian sages.

Abstract

Salvia abrotanoides Kar. and Salvia yangii B.T. Drew are medicinal and aromatic plants belonging to the subgenus Perovskia and used as herbal medicines in Asia. Derivatives of caffeic acid, mainly rosmarinic acid (RA), are the major phenolic compounds identified in these plants. Understanding the factors and molecular mechanisms regulating the accumulation of pharmacologically and ecologically relevant phenolic metabolites is essential for future biotechnological and medical applications. Up to date, no studies of phenylpropanoid biosynthetic pathway at the transcriptional level has been performed in the Perovskia subgenus. Using a combined qRT-PCR transcriptional activity analysis with LC-MS based metabolic profiling of roots and leaves at the beginning, in the middle and at the end of vegetation season, we have identified the following gene candidates with properties correlating to phenolic acid biosynthesis in S. abrotanoides and S. yangii: PAL, C4H, 4CL, TAT, HPPR, RAS1, RAS2 and Cyp98A14. A comparison of phenolic acid profiles with gene transcript levels revealed the transcriptional regulation of RA biosynthesis in the roots but not the leaves of the studied species. Additionally, RAS1 and Cyp98A14 were identified as rate-limiting steps regulating phenylpropanoid biosynthesis on a transcription level. In the future, this will facilitate the gene-based metabolic enhancement of phenolic compounds production in these promising medicinal herbs.

Transcription Factor: A Powerful Tool to Regulate Biosynthesis of Active Ingredients in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge is a common Chinese herbal medicine, and its major active ingredients are phenolic acids and tanshinones, which are widely used to treat vascular diseases. However, the wild form of S. miltiorrhiza possess low levels of these important pharmaceutical agents; thus, improving their levels is an active area of research. Transcription factors, which promote or inhibit the expressions of multiple genes involved in one or more biosynthetic pathways, are powerful tools for controlling gene expression in biosynthesis. Several families of transcription factors have been reported to participate in regulating phenolic acid and tanshinone biosynthesis and influence their accumulation. This review summarizes the current status in this field, with focus on the transcription factors which have been identified in recent years and their functions in the biosynthetic regulation of phenolic acids and tanshinones. Otherwise, the new insight for further research is provided. Finally, the application of the biosynthetic regulation of active ingredients by the transcription factors in S. miltiorrhiza are discussed, and new insights for future research are explored.

Reprogramming plant specialized metabolism by manipulating protein kinases

Being sessile, plants have evolved sophisticated mechanisms to balance between growth and defense to survive in the harsh environment. The transition from growth to defense is commonly achieved by factors, such as protein kinases (PKs) and transcription factors, that initiate signal transduction and regulate specialized metabolism. Plants produce an array of lineage-specific specialized metabolites for chemical defense and stress tolerance. Some of these molecules are also used by humans as drugs. However, many of these defense-responsive metabolites are toxic to plant cells and inhibitory to growth and development. Plants have, thus, evolved complex regulatory networks to balance the accumulation of the toxic metabolites. Perception of external stimuli is a vital part of the regulatory network. Protein kinase-mediated signaling activates a series of defense responses by phosphorylating the target proteins and translating the stimulus into downstream cellular signaling. As biosynthesis of specialized metabolites is triggered when plants perceive stimuli, a possible connection between PKs and specialized metabolism is well recognized. However, the roles of PKs in plant specialized metabolism have not received much attention until recently. Here, we summarize the recent advances in understanding PKs in plant specialized metabolism. We aim to highlight how the stimulatory signals are transduced, leading to the biosynthesis of corresponding metabolites. We discuss the post-translational regulation of specialized metabolism and provide insights into the mechanisms by which plants respond to the external signals. In addition, we propose possible strategies to increase the production of plant specialized metabolites in biotechnological applications using PKs.

Transcriptome analysis of Curcuma wenyujin from Haikou and Wenzhou, and a comparison of the main constituents and related genes of Rhizoma Curcumae

For more than a thousand years, Rhizoma Curcumae (known as E zhu), a Chinese herbal medicine, has been used to eradicate blood stasis and relieve aches. The plant Curcuma wenyujin, which is grown primarily in Wenzhou, China, is considered the best source of Rhizoma Curcumae. In this study, we sought to ascertain differences in transcript profiles of C. wenyujin grown in traditional (Wenzhou) and recently established (Haikou) production areas based on Illumina and RNA (RNA-seq) sequencing. We also examined differences in the main components of the volatile oil terpene; curcumin, polysaccharide, and starch constituents and related genes in the corresponding pathways, in C. wenyujin cultivated in the two production areas. We accordingly found that the essential oil (2.05%), curcumin (1.46%), and polysaccharide (8.90%) content in Wenzhou rhizomes was higher than that in the rhizomes of plants from Haikou (1.60%, 0.91%, and 6.15%, respectively). In contrast, the starch content of Wenzhou rhizomes (17.0%) was lower than that of Haikou rhizomes (23.8%). Furthermore, we detected significant differences in the oil components of Haikou and Wenzhou rhizomes, with curzerene (32.34%), curdione (21.35%), and germacrene B (9.39%) being the primary components of the essential oil derived from Wenzhou rhizomes, and curzerene (20.13%), curdione (14.73%), and cineole (9.76%) being the main constituents in Haikou rhizomes. Transcriptome and qPCR analyses revealed considerable differences in gene expression between Wenzhou and Haikou rhizomes. The expression of terpene, curcumin, and polysaccharide pathway-related genes in Wenzhou rhizomes was significantly up-regulated, whereas the expression of starch-associated genes was significantly down-regulated, compared with those in Haikou rhizomes. Difference in the content of terpene, curcumin, polysaccharides, and starch in rhizomes from the two production areas could be explained in terms of differences in expression of the related genes.

Cross-Species Annotation of Expressed Genes and Detection of Different Functional Gene Modules Between 10 Cold- and 10 Hot-Propertied Chinese Herbal Medicines

According to the traditional Chinese medicine (TCM) system, Chinese herbal medicines (HMs) can be divided into four categories: hot, warm, cold, and cool. A cool nature usually is categorized as a cold nature, and a warm nature is classified as a hot nature. However, the detectable characteristics of the gene expression profile associated with the cold and hot properties have not been studied. To address this question, a strategy for the cross-species annotation of conserved genes was established in the present study by using transcriptome data of 20 HMs with cold and hot properties. Functional enrichment analysis was performed on group-specific expressed genes inferred from the functional genome of the reference species (i.e., Arabidopsis). Results showed that metabolic pathways relevant to chrysoeriol, luteolin, paniculatin, and wogonin were enriched for cold-specific genes, and pathways of inositol, heptadecane, lauric acid, octanoic acid, hexadecanoic acid, and pentadecanoic acid were enriched for hot-specific genes. Six functional modules were identified in the HMs with the cold property: nucleotide biosynthetic process, peptidy-L-cysteine S-palmitoylation, lipid modification, base-excision repair, dipeptide transport, and response to endoplasmic reticulum stress. For the hot HMs, another six functional modules were identified: embryonic meristem development, embryonic pattern specification, axis specification, regulation of RNA polymerase II transcriptional preinitiation complex assembly, mitochondrial RNA modification, and cell redox homeostasis. The research provided a new insight into HMs’ cold and hot properties from the perspective of the gene expression profile of plants.

Polyphenolic Composition and Anti-Melanoma Activity of White Forsythia (Abeliophyllum distichum Nakai) Organ Extracts

Abeliophyllum distichum Nakai, commonly called white forsythia, is a monotypic genus endemic to Korea. Although A. distichum is mainly used as an ornamental plant because of its horticultural value, recent studies have demonstrated its bioactivities, including antioxidant and anti-inflammatory activities, prompting us to investigate the potential anticancer effect of A. distichum organ extracts (leaves, fruit, and branches) against human melanoma SK-MEL-2 cells. The methanol extract of A. distichum leaves (AL) exhibited dose- and time-dependent cytotoxicities against SK-MEL-2 cells but not against HDFa human dermal fibroblasts. Based on high-performance liquid chromatography analysis, we identified 18 polyphenolic compounds from A. distichum organ extracts and suggest that differences in anticancer activity between organ extracts should be caused by different compositions of polyphenolic compounds. Additionally, the Annexin V/propidium iodide staining assay and analysis of caspase activity and expression indicated that AL induced cell death, including early and late apoptosis, as well as necrosis, by inducing the extrinsic pathway. Furthermore, we analyzed the differentially expressed genes between mock- and AL-treated cells using RNA-seq technology, suggesting that the anti-melanoma action of AL is mediated by down-regulation of the phosphoinositide 3-kinase (PI3K)/Akt signaling pathway. Taken together, these results shed light on the potential use of A. distichum as a green resource with potent anti-melanoma activity.

Full-length transcriptome analysis of Coptis deltoidea and identification of putative genes involved in benzylisoquinoline alkaloids biosynthesis based on combined sequencing platforms

The study carry out comprehensive transcriptome analysis of C. deltoidea and exploration of BIAs biosynthesis and accumulation based on UHPLC-MS/MS and combined sequencing platforms. Coptis deltoidea is an important medicinal plant with a long history of medicinal use, which is rich in benzylisoquinoline alkaloids (BIAs). In this study, Ultra performance liquid chromatography-electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) and combined sequencing platforms were performed for exploration of BIAs biosynthesis, accumulation and comprehensive transcriptome analysis of C. deltoidea. By metabolism profiling, the accumulation of ten BIAs was analyzed using UHPLC-MS/MS and different contents were observed in different organs. From transcriptome sequencing result, we applied single-molecule real-time (SMRT) sequencing to C. deltoidea and generated a total of 75,438 full-length transcripts. We proposed the candidate biosynthetic pathway of tyrosine, precursor of BIAs, and identified 64 full length-transcripts encoding enzymes putatively involved in BIAs biosynthesis. RNA-Seq data indicated that the majority of genes exhibited relatively high expression level in roots. Transport of BIAs was also important for their accumulation. Here, 9 ABC transporters and 2 MATE transporters highly homologous to known alkaloid transporters related with BIAs transport in roots and rhizomes were identified. These findings based on the combined sequencing platforms provide valuable genetic information for C. deltoidea and the results of transcriptome combined with metabolome analysis can help us better understand BIAs biosynthesis and transport in this medicinal plant. The information will be critical for further characterization of C. deltoidea transcriptome and molecular-assisted breeding for this medicinal plant with scarce resources.

The ethylene response factor SmERF8 regulates the expression of SmKSL1 and is involved in tanshinone biosynthesis in Saliva miltiorrhiza hairy roots

Saliva miltiorrhiza ethylene response factor (SmERF), predicted to be expressed genome-wide, is the potential regulator of tanshinone biosynthesis. However, few studies have investigated its transcriptional regulation pathways in tanshinone biosynthesis. Here, we report an ethylene response factor (SmERF8), which was screened by the SmKSL1 (a key gene in tanshinone biosynthesis) promoter from the S. miltiorrhiza cDNA library. The SmERF8, highly expressed in S. miltiorrhiza root head, is sensitive to Eth stress, and its protein was enriched in the nucleus. The SmERF8 recognizes the GCC-box in the SmKSL1 promoter. Overexpression and RNAi of SmERF8 in S. miltiorrhiza transgenic hairy roots showed that the tanshinone contents were significantly increased in the overexpression transgenic lines and decreased in RNAi lines. These results suggest that the SmERF8 may be a central activator that regulates the expression of SmKSL1 by binding the GCC-box and then promoting tanshinone biosynthesis. Thus, the SmERF8 may functionally accelerate tanshinone biosynthesis by the transcriptional regulation of its key gene.

Assembly and Annotation of a Draft Genome of the Medicinal Plant Polygonum cuspidatum

Polygonum cuspidatum (Japanese knotweed, also known as Huzhang in Chinese), a plant that produces bioactive components such as stilbenes and quinones, has long been recognized as important in traditional Chinese herbal medicine. To better understand the biological features of this plant and to gain genetic insight into the biosynthesis of its natural products, we assembled a draft genome of P. cuspidatum using Illumina sequencing technology. The draft genome is ca. 2.56 Gb long, with 71.54% of the genome annotated as transposable elements. Integrated gene prediction suggested that the P. cuspidatum genome encodes 55,075 functional genes, including 6,776 gene families that are conserved in the five eudicot species examined and 2,386 that are unique to P. cuspidatum. Among the functional genes identified, 4,753 are predicted to encode transcription factors. We traced the gene duplication history of P. cuspidatum and determined that it has undergone two whole-genome duplication events about 65 and 6.6 million years ago. Roots are considered the primary medicinal tissue, and transcriptome analysis identified 2,173 genes that were expressed at higher levels in roots compared to aboveground tissues. Detailed phylogenetic analysis demonstrated expansion of the gene family encoding stilbene synthase and chalcone synthase enzymes in the phenylpropanoid metabolic pathway, which is associated with the biosynthesis of resveratrol, a pharmacologically important stilbene. Analysis of the draft genome identified 7 abscisic acid and water deficit stress-induced protein-coding genes and 14 cysteine-rich transmembrane module genes predicted to be involved in stress responses. The draft de novo genome assembly produced in this study represents a valuable resource for the molecular characterization of medicinal compounds in P. cuspidatum, the improvement of this important medicinal plant, and the exploration of its abiotic stress resistance.

Metabolome and transcriptome analyses reveal quality change in the orange-rooted Salvia miltiorrhiza (Danshen) from cultivated field

Background

The dry root and rhizome of Salvia miltiorrhiza Bunge, or Danshen, is a well-known, traditional Chinese medicine. Tanshinones are active compounds that accumulate in the periderm, resulting in red-colored roots. However, lines with orange roots have been observed in cultivated fields. Here, we performed metabolome and transcriptome analyses to investigate the changes of orange-rooted Danshen.

Methods

Metabolome analysis was performed by ultra-high-performance liquid chromatography-quadrupole time-of-flight mass spectrometry (UPLC/Q-Tof–MS) to investigate the metabolites variation between orange Danshen and normal Danshen. RNA sequencing and KEGG enrichment analysis were performed to analyzing the differentially expressed genes between orange-rooted and normal Danshen.

Results

In total, 40 lipophilic components were detected in metabolome analysis, and seven compounds were significantly decreased in the orange Danshen, including the most abundant active compounds, tanshinone IIA and tanshinone I in normal Danshen. Systematic analysis of transcriptome profiles revealed that the down-regulated genes related to catalytic dehydrogenation was not detected. However, two genes related to stress resistance, and four genes related to endoplasmic reticulum (ER)-associated degradation of proteins were up-regulated in orange Danshen.

Conclusions

Decreases in the content of dehydrogenated furan ring tanshinones such as tanshinone IIA resulted in phenotypic changes and quality degradation of Danshen. Transcriptome analysis indicated that incorrect folding and ER-associated degradation of corresponding enzymes, which could catalyze C₁₅-C₁₆ dehydrogenase, might be contributed to the decrease in dehydrogenated furan ring tanshinones, rather than lower expression of the relative genes. This limited dehydrogenation of cryptotanshinone and dihydrotanshinone I into tanshinones IIA and I products, respectively, led to a reduced quality of Danshen in cultivated fields.

Polysaccharide Fraction Extracted from Endophytic Fungus Trichoderma atroviride D16 Has an Influence on the Proteomics Profile of the Salvia miltiorrhiza Hairy Roots

Trichoderma atroviride develops a symbiont relationship with Salvia miltiorrhiza and this association involves a number of signaling pathways and proteomic responses between both partners. In our previous study, we have reported that polysaccharide fraction (PSF) of T. atroviride could promote tanshinones accumulation in S.miltiorrhiza hairy roots. Consequently, the present data elucidates the broad proteomics changes under treatment of PSF. Furthermore, we reported several previously undescribed and unexpected responses, containing gene expression patterns consistent with biochemical stresses and metabolic patterns inside the host. In summary, the PSF-induced tanshinones accumulation in S.miltiorrhiza hairy roots may be closely related to Ca2+ triggering, peroxide reaction, protein phosphorylation, and jasmonic acid (JA) signal transduction, leading to an increase in leucine-rich repeat (LRR) protein synthesis. This results in the changes in basic metabolic flux of sugars, amino acids, and protein synthesis, along with signal defense reactions. The results reported here increase our understanding of the interaction between T.atroviride and S.miltiorrhiza and specifically confirm the proteomic responses underlying the activities of PSF.

Computer-Aided Saturation Mutagenesis of Arabidopsis thaliana Ent-Copalyl Diphosphate Synthase

Ent-copalyl diphosphate synthase controls the biosynthesis of gibberellin plant hormones, which in turn coordinate the expression of numerous enzymes. Some gibberellin-dependent genes encode enzymes coordinating the biosynthesis of tanshinones: diterpene derivatives with broad medical applications. New biotechnological approaches, such as metabolic engineering using naturally occurring or mutated enzymes, have been proposed to meet the growing demand for tanshinones which is currently met by the Chinese medicinal plant Salvia miltiorrhiza Bunge. These mutants may be prepared by directed evolution, saturation mutagenesis or rational enzyme design. In the presented paper, 15,257 non-synonymous variants of Arabidopsis thaliana ent-copalyl diphosphate synthase were obtained using the SNAP2 tool. The obtained forms were screened to isolate variants with potentially improved biological functions. A group of 455 mutants with potentially improved stability was isolated and subjected to further screening on the basis of ligand–substrate affinity, and both secondary structure and active site structure stability. Finally, a group of six single mutants was obtained, which were used to construct double mutants with potentially improved stability and ligand affinity. The potential influence of single mutations on protein stability and ligand affinity was evaluated by double mutant cycle analysis. Finally, the procedure was validated by in silico assessment of the experimentally verified enzyme mutants with reduced enzymatic activity.

Transcriptome Analysis Shows That IFN-I Treatment and Concurrent SAV3 Infection Enriches MHC-I Antigen Processing and Presentation Pathways in Atlantic Salmon-Derived Macrophage/Dendritic Cells

Type I interferons (IFNs) have been shown to play an important role in shaping adaptive immune responses in addition to their antiviral properties in immune cells. To gain insight into the impact of IFN-I-induced pathways involved in early adaptive immune responses, i.e., antigen-presenting pathways, in an Atlantic salmon-derived (Salmo salar L.) macrophage cell line (TO-cells), we used a comparative de novo transcriptome analysis where cells were treated with IFN-I or kept untreated and concurrently infected with salmonid alphavirus subtype 3 (SAV3). We found that concurrent treatment of TO-cells with IFN-I and SAV3 infection (SAV3/IFN⁺) significantly enriched the major histocompatibility complex class I (MHC-I) pathway unlike the non-IFN-I treated TO-cells (SAV3/IFN⁻) that had lower expression levels of MHC-I pathway-related genes. Genes such as the proteasomal activator (PA28) and β-2 microglobulin (β2M) were only differentially expressed in the SAV3/IFN⁺ cells and not in the SAV3/IFN⁻ cells. MHC-I pathway genes like heat shock protein 90 (Hsp90), transporter of antigen associated proteins (TAPs) and tapasin had higher expression levels in the SAV3/IFN⁺ cells than in the SAV3/IFN⁻ cells. There were no MHC-II pathway-related genes upregulated in SAV3/IFN⁺-treated cells, and cathepsin S linked to the degradation of endosomal antigens in the MHC-II pathway was downregulated in the SAV3/IFN⁻ cells. Overall, our findings show that concurrent IFN-I treatment of TO-cells and SAV3 infection enriched gene expression linked to the MHC-I antigen presentation pathway. Data presented indicate a role of type I IFNs in strengthening antigen processing and presentation that may facilitate activation particularly of CD8+ T-cell responses following SAV3 infection, while SAV3 infection alone downplayed MHC-II pathways.

De Novo Assembled Transcriptome Analysis and Identification of Genic SSR Markers in Red-Flowered Strawberry

Red-flowered strawberry is a new ornamental flower derived from intergeneric hybridization (Fragaria × Potentilla). To date, few molecular markers have been reported for this plant. RNA sequencing provides a relatively fast and low-cost approach for large-scale detection of simple sequence repeats (SSRs). In the present study, we profiled the transcriptome of red-flowered strawberry by Illumina HiSeq 2500 to identify SSRs related to petal color. Based on 2 million clean reads of red and white flowers from red-flowered strawberry hybrids, we assembled 91,835 unigenes with an average length of 717 bp. After functional annotation and prediction, there were 47,058 unigenes; of these, 26,861 had a gene ontology annotation, with 14,264 SSR loci. Mononucleotide SSRs were the predominant repeat type (47.20%, n = 6724), followed by di- (32.50%, n = 4641), tri- (19.10%, n = 2729), tetra- (0.90%, n = 132), hexa- (0.2%, n = 21), and penta- (0.10%, n = 16) nucleotide repeats. The most frequent di-, tri-, and tetra-nucleotide repeats were AG/CT, AAG/CTT, and AAAG/CTTT, respectively. PCR amplification with 105 SSR primer pairs yielded four bands specific to red flowers, namely UgRFsr57622, UgRFsr94149, UgRFsr40142, and UgRFsr54608; corresponding 4 trait-specific markers were found to co-segregate with white and red flower color in hybrid population, demonstrating that the genic SSR marker is useful to discriminate between white and red flowers in strawberry. Markers to discriminate flower color in red-flowered strawberry will be useful for early selection of progeny and for breeding management.

Comparative Transcriptome Analysis Identifies Genes Involved in Diosgenin Biosynthesis in Trigonella foenum-graecum L

Trigonella foenum-graecum L. (fenugreek) is a valuable resource of producing diosgenin which serves as a substrate for synthesizing more than two hundred kinds of steroidal drugs. Phytochemical analysis indicated that methyl jasmonate (MeJA) efficiently induced diosgenin biosynthesis in fenugreek seedlings. Though early steps up to cholesterol have recently been elucidated in plants, cytochrome P450 (CYP)- and glycosyltransferase (GT)-encoding genes involved in the late steps from cholesterol to diosgenin remain unknown. This study established comparative fenugreek transcriptome datasets from the MeJA-treated seedlings and the corresponding control lines. Differential gene expression analysis identified a number of MeJA-induced CYP and GT candidate genes. Further gene expression pattern analysis across a different MeJA-treating time points, together with a phylogenetic analysis, suggested specific family members of CYPs and GTs that may participate in the late steps during diosgenin biosynthesis. MeJA-induced transcription factors (TFs) that may play regulatory roles in diosgenin biosynthesis were also discussed. This study provided a valuable genetic resource to functionally characterize the genes involved in diosgenin biosynthesis, which will push forward the production of diosgenin in microbial organisms using a promising synthetic biology strategy.

The biosynthesis of phenolic acids is positively regulated by the JA-responsive transcription factor ERF115 in Salvia miltiorrhiza

Phenolic acids are important secondary metabolites produced in the Chinese medicinal plant Salvia miltiorrhiza, but little is known about the transcription factors involved in the regulation of tanshinone and phenolic acid biosynthesis. Here, a novel AP2/ERF transcription factor SmERF115 was isolated and functionally characterized. SmERF115 was most responsive to methyl jasmonate (MeJA) treatment and was localized in the nucleus. The phenolic acid production was increased in SmERF115-overexpressing hairy roots, but with a decrease in tanshinone content. In contrast, silencing of SmERF115 reduced the phenolic acid level, but increased tanshinone content. The expression of the key biosynthetic gene SmRAS1 was up-regulated in SmERF115 overexpression lines but was down-regulated in SmERF115-RNAi lines. Yeast one-hybrid (Y1H) assay and EMSA showed that SmERF115 directly binds to the promoter of SmRAS1, while dual-luciferase assays showed that SmERF115 could activate expression of SmRAS1 in vivo. Furthermore, global transcriptomic analysis by RNA sequencing revealed that expression of other genes such as PAL3, 4CL5, TAT3, and RAS4 was also increased in the overexpression line, implying that they were potentially involved in the SmERF115-mediated pathway. Our data show that SmERF115 is a positive regulator of phenolic acid biosynthesis, and may be a potential target for further metabolic engineering of phenolic acid biosynthesis in S. miltiorrhiza.

De novo transcriptome sequencing and metabolite profiling analyses reveal the complex metabolic genes involved in the terpenoid biosynthesis in Blue Anise Sage (Salvia guaranitica L.)

Many terpenoid compounds have been extracted from different tissues of Salvia guaranitica. However, the molecular genetic basis of terpene biosynthesis pathways is virtually unknown. In this study, approximately 4 Gb of raw data were generated from the transcriptome of S. guaranitica leaves using Illumina HiSeq 2000 sequencing. After filtering and removing the adapter sequences from the raw data, the number of reads reached 32 million, comprising 186 million of high-quality nucleotide bases. A total of 61,400 unigenes were assembled de novo and annotated for establishing a valid database for studying terpenoid biosynthesis. We identified 267 unigenes that are putatively involved in terpenoid metabolism (including, 198 mevalonate and methyl-erythritol phosphate (MEP) pathways, terpenoid backbone biosynthesis genes and 69 terpene synthases genes). Moreover, three terpene synthase genes were studied for their functions in terpenoid biosynthesis by using transgenic Arabidopsis; most transgenic Arabidopsis plants expressing these terpene synthetic genes produced increased amounts of terpenoids compared with wild-type control. The combined data analyses from the transcriptome and metabolome provide new insights into our understanding of the complex metabolic genes in terpenoid-rich blue anise sage, and our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. guaranitica.

Molecular Characterization and Overexpression of SmJMT Increases the Production of Phenolic Acids in Salvia miltiorrhiza

Jasmonic acid (JA) carboxyl methyltransferase (JMT), a key enzyme in jasmonate-regulated plant responses, may be involved in plant defense and development by methylating JA to MeJA, thus influencing the concentrations of MeJA in plant. In this study, we isolated the JMT gene from Salvia miltiorrhiza, an important medicinal plant widely used to treat cardiovascular disease. We present a genetic manipulation strategy to enhance the production of phenolic acids by overexpresion SmJMT in S. miltiorrhiza. Global transcriptomic analysis using RNA sequencing showed that the expression levels of genes involved in the biosynthesis pathway of phenolic acids and MeJA were upregulated in the overexpression lines. In addition, the levels of endogenous MeJA, and the accumulation of rosmarinic acid (RA) and salvianolic acid (Sal B), as well as the concentrations of total phenolics and total flavonoids in transgenic lines, were significantly elevated compared with the untransformed control. Our results demonstrate that overexpression of SmJMT promotes the production of phenolic acids through simultaneously activating genes encoding key enzymes involved in the biosynthesis pathway of phenolic acids and enhancing the endogenous MeJA levels in S. miltiorrhiza.

AtPAP1 Interacts With and Activates SmbHLH51, a Positive Regulator to Phenolic Acids Biosynthesis in Salvia miltiorrhiza

Phenolic acids from Salvia miltiorrhiza have drawn considerable attention in recent years because of their remarkable pharmacological activities. We previously reported that Arabidopsis thaliana transcription factor production of anthocyanin pigment 1 (AtPAP1) has strong capability to promote the production of phenolic acids in S. miltiorrhiza. However, the responsible molecular mechanism is unclear. Here, we analyzed the transcriptome of transgenic S. miltiorrhiza that over-expressed AtPAP1. Transcriptome analysis revealed 4,152 genes that were differentially expressed due to ectopic AtPAP1 overexpression. SmbHLH51, a novel bHLH gene significantly up-regulated by constitutive expression of AtPAP1, was isolated from S. miltiorrhiza for detailed functional characterization. SmbHLH51 localizes in the nuclei and interacts with AtPAP1, indicating that they probably comprise a regulatory transcription complex. Enhanced or reduced expression of SmbHLH51 was achieved in S. miltiorrhiza by gain- or loss-of-function assays, respectively, revealing that SmbHLH51 is a positive transcriptional regulator of the pathway for phenolic acid biosynthesis. We propose that applying this functional genomics approach through the transcriptomic analyses is an efficient means for identifying novel genes involved in plant secondary metabolism.

Identification of Glutathione Peroxidase (GPX) Gene Family in Rhodiola crenulata and Gene Expression Analysis under Stress Conditions

Glutathione peroxidases (GPXs) are important enzymes in the glutathione-ascorbate cycle for catalyzing the reduction of H₂O₂ or organic hydroperoxides to water. GPXs play an essential role in plant growth and development by participating in photosynthesis, respiration, and stress tolerance. Rhodiola crenulata is a popular traditional Chinese medicinal plant which displays an extreme energy of tolerance to harsh alpine climate. The GPXs gene family might provide R. crenulata for extensively tolerance to environment stimulus. In this study, five GPX genes were isolated from R. crenulata. The protein amino acid sequences were analyzed by bioinformation softwares with the results that RcGPXs gene sequences contained three conserve cysteine residues, and the subcellular location predication were in the chloroplast, endoplasmic reticulum, or cytoplasm. Five RcGPXs members presented spatial and temporal specific expression with higher levels in young and green organs. And the expression patterns of RcGPXs in response to stresses or plant hormones were investigated by quantitative real-time PCR. In addition, the putative interaction proteins of RcGPXs were obtained by yeast two-hybrid with the results that RcGPXs could physically interact with specific proteins of multiple pathways like transcription factor, calmodulin, thioredoxin, and abscisic acid signal pathway. These results showed the regulation mechanism of RcGPXs were complicated and they were necessary for R. crenulata to adapt to the treacherous weather in highland.

Tanshinone IIA Attenuates Insulin Like Growth Factor 1 -Induced Cell Proliferation in PC12 Cells through the PI3K/Akt and MEK/ERK Pathways

The insulin like growth factor 1 (IGF-1) and its receptor (IGF-1R) facilitate tumor proliferation and progression. Tanshinone IIA (TSN) is an active diterpene quinone isolated from the roots of the herbal plant Salvia miltiorrhiza. TSN inhibits the proliferation of various types of cancer cells but its role in the IGF-1R-induced proliferation of pheochromocytoma (PC12) cells and the potential mechanisms are largely unknown. This study aims to investigate the anti-proliferative effect of TSN in PC12 cells and its role on IGF-1R signaling transduction. PC12 cells were treated with IGF-1 with or without TSN, methyl thiazolytetrazolium (MTT) assay, and cell counting kit-8 and flow cytometry were used to evaluate the proliferation of PC12 cells. The role of TSN on the apoptosis of PC12 cells were detected by flow cytometry as well. The effects of TSN and IGF-1 on the phosphorylation of IGF-1R, protein kinase B (Akt), extracellular-signal related kinase 1/2 (ERK1/2) and other downstream targets were analyzed by Western blotting analysis. Our results showed that IGF-1 promoted the growth of PC12 cells in a dose-dependent manner and increased the phosphorylation of IGF-1R, whereas TSN attenuated the effect of IGF-1. Interestingly, TSN did not induce cell apoptosis in PC12 cells. Moreover, TSN attenuated the phosphorylation of Akt and ERK1/2 induced by IGF-1, and the phosphorylation of glycogen synthase kinase-3β, forkhead box O3a (FOXO3a) and c-Raf were also inhibited by TSN. Furthermore, TSN inhibited cell growth induced by IGF-1 and blocked the activation of IGF-1R in SH-SY5Y cells. Taken together, TSN has an inhibitory effect on the proliferation of PC12 cells via down-regulation of the phosphorylated IGF-1R and its downstream signaling.

Anticancer potential of rosmarinic acid and its improved production through biotechnological interventions and functional genomics

Rosmarinic acid (RA) is a highly valued natural phenolic compound that is very commonly found in plants of the families Lamiaceae and Boraginaceae, including Coleus blumei, Heliotropium foertherianum, Rosmarinus officinalis, Perilla frutescens, and Salvia officinalis. RA is also found in other members of higher plant families and in some fern and horned liverwort species. The biosynthesis of RA is catalyzed by the enzymes phenylalanine ammonia lyase and cytochrome P450-dependent hydroxylase using the amino acids tyrosine and phenylalanine. Chemically, RA can be produced via methods involving the esterification of 3,4-dihydroxyphenyllactic acid and caffeic acid. Some of the derivatives of RA include melitric acid, salvianolic acid, lithospermic acid, and yunnaneic acid. In plants, RA is known to have growth-promoting and defensive roles. Studies have elucidated the varied pharmacological potential of RA and its derived molecules, including anticancer, antiangiogenic, anti-inflammatory, antioxidant, and antimicrobial activities. The demand for RA is therefore, very high in the pharmaceutical industry, but this demand cannot be met by plants alone because RA content in plant organs is very low. Further, many plants that synthesize RA are under threat and near extinction owing to biodiversity loss caused by unscientific harvesting, over-collection, environmental changes, and other inherent features. Moreover, the chemical synthesis of RA is complicated and expensive. Alternative approaches using biotechnological methodologies could overcome these problems. This review provides the state of the art information on the chemistry, sources, and biosynthetic pathways of RA, as well as its anticancer properties against different cancer types. Biotechnological methods are also discussed for producing RA using plant cell, tissue, and organ cultures and hairy-root cultures using flasks and bioreactors. The recent developments and applications of the functional genomics approach and heterologous production of RA in microbes are also highlighted. This chapter will be of benefit to readers aiming to design studies on RA and its applicability as an anticancer agent.

Transcriptome Analyses from Mutant Salvia miltiorrhiza Reveals Important Roles for SmGASA4 during Plant Development

Salvia miltiorrhiza (S. miltiorrhiza) is an important Chinese herb that is derived from the perennial plant of Lamiaceae, which has been used to treat neurasthenic insomnia and cardiovascular disease. We produced a mutant S. miltiorrhiza (MT), from breeding experiments, that possessed a large taproot, reduced lateral roots, and defective flowering. We performed transcriptome profiling of wild type (WT) and MT S. miltiorrhiza using second-generation Illumina sequencing to identify differentially expressed genes (DEGs) that could account for these phenotypical differences. Of the DEGs identified, we investigated the role of SmGASA4, the expression of which was down-regulated in MT plants. SmGASA4 was introduced into Arobidopsis and S. militiorrhiza under the control of a CaMV35S promoter to verify its influence on abiotic stress and S. miltiorrhiza secondary metabolism biosynthesis. SmGASA4 was found to promote flower and root development in Arobidopsis. SmGASA4 was also found to be positively regulated by Gibberellin (GA) and significantly enhanced plant resistance to salt, drought, and paclobutrazol (PBZ) stress. SmGASA4 also led to the up-regulation of the genes involved in salvianolic acid biosynthesis, but inhibited the expression of the genes involved in tanshinone biosynthesis. Taken together, our results reveal SmGASA4 as a promising candidate gene to promote S. miltiorrhiza development.

Roles of Plant-Associated Microbiota in Traditional Herbal Medicine

The microbiome of medicinal plants may directly impact the metabolome of the host, and thus could influence the efficacy of herbal medicine. We advocate a herb for traditional Chinese medicine (TCM), Salvia miltiorrhiza, as a prime model system to study how microbes may interact with medicinal plants to modify phytochemical production.

Identification of genes associated with stress tolerance in moth bean [Vigna aconitifolia (Jacq.) Marechal], a stress hardy crop

Moth bean is the most drought and heat tolerant cultigens among Asian Vigna. We performed comparative transcriptome analysis of moth bean cultivar "Marumoth" under control and stress condition. De novo transcriptome assembly was carried out by using Velvet followed by Oases softwares. Differential expression analyses, SSR identification and validation and mapping of pathways and transcription factors were conducted. A total of 179,979 and 201,888 reads were generated on Roche 454 platform and 48,617,205 and 45,449,053 reads were generated on ABI Solid platform for the control and stressed samples. Combined assembly from Roche and ABI Solid platforms generated 16,090 and 15,096 transcripts for control and stressed samples. We found 1287 SSRs and 5606 transcripts involved in 179 pathways. The 55 transcription factor families represented 19.42% of total mothbean transcripts. In expression profiling, ten transcripts were found to be up-regulated and 41 down-regulated while 490 showed no major change under moisture stress condition. Stress inducible genes like Catalase, Cyt P450 monooxygenase, heat shock proteins (HSP 90 and HSP 70), oxidoreductase, protein kinases, dehydration responsive protein (DRP), universal stress protein and ferridoxin NADH oxidoreductase genes were up-regulated in stressed sample. Genes which might be involved in moisture stress tolerance in moth bean were identified and these might be useful for stress tolerance breeding in moth bean and other related crops.

The ethylene response factor SmERF6 co-regulates the transcription of SmCPS1 and SmKSL1 and is involved in tanshinone biosynthesis in Salvia miltiorrhiza hairy roots

The SmERF6, which recognizes the GCC-box of SmCPS1 and SmKSL1 promoter in nucleus, regulates the tanshinone biosynthesis in Salvia miltiorrhiza hairy roots. Tanshinone, an important medicinal ingredient in Salvia miltiorrhiza, is best known for its use in medicine. However, the transcription factor regulation of tanshinone biosynthesis is unclear. Here, we isolated and identified a transcription factor in the ERF family of S. miltiorrhiza, SmERF6, which was screened from an S. miltiorrhiza cDNA library by the promoters of two key tanshinone synthesis genes (SmKSL1 and SmCPS1); this factor regulated tanshinone biosynthesis. The gene was highly expressed in the root and responded to ethylene treatment. SmERF6 modulated tanshinone biosynthesis by directly binding to an ethylene-responsive element (GCC-box) of the SmKSL1 and SmCPS1 promoters and activating their transcription. Overexpression of SmERF6 in the hairy roots increased their tanshinone accumulation, and SmERF6 silencing by RNAi led to a lower tanshinone content. Furthermore, tanshinone accumulation maintained homeostasis with the total phenolic acid and flavonoid contents in S. miltiorrhiza. These findings elucidated how SmERF6 directly co-regulates the transcription of SmCPS1 and SmKSL1 and modulates tanshinone synthesis to accelerate the metabolic flux of tanshinone accumulation in S. miltiorrhiza.

Genome survey sequencing of Dioscorea zingiberensis

Dioscorea zingiberensis (Dioscoreceae) is the main plant source of diosgenin (steroidal sapogenins), the precursor for the production of steroid hormones in the pharmaceutical industry. Despite its large economic value, genomic information of the genus Dioscorea is currently unavailable. Here, we present an initial survey of the D. zingiberensis genome performed by next-generation sequencing technology together with a genome size investigation inferred by flow cytometry. The whole genome survey of D. zingiberensis generated 31.48 Gb of sequence data with approximately 78.70× coverage. The estimated genome size is 800 Mb, with a high level of heterozygosity based on K-mer analysis. These reads were assembled into 334 288 contigs with a N50 length of 1079 bp, which were further assembled into 92 163 scaffolds with a total length of 173.46 Mb. A total of 4935 genes, 81 tRNAs, 69 rRNAs, and 661 miRNAs were predicted by the genome analysis, and 263 484 repeated sequences were obtained with 419 372 simple sequence repeats (SSRs). Among these SSRs, the mononucleotide repeat type was the most abundant (up to 54.60% of the total SSRs), followed by the dinucleotide (29.60%), trinucleotide (11.37%), tetranucleotide (3.53%), pentanucleotide (0.65%), and hexanucleotide (0.25%) repeat types. The 1C-value of D. zingiberensis was calibrated against Salvia miltiorrhiza and calculated as 0.87 pg (851 Mb) by flow cytometry, which was very close to the result of the genome survey. This is the first report of genome-wide characterization within this taxon.

Transcriptome and metabolite analyses reveal the complex metabolic genes involved in volatile terpenoid biosynthesis in garden sage (Salvia officinalis)

A large number of terpenoid compounds have been extracted from different tissues of S. officinalis. However, the molecular genetic basis of terpene biosynthesis pathways is virtually unknown. In this study, approximately 6.6 Gb of raw data were generated from the transcriptome of S. officinalis leaves using Illumina HiSeq 2000 sequencing. After filtering and removing the adapter sequences from the raw data, the number of reads reached 21 million, comprising 98 million of high-quality nucleotide bases. 48,671 unigenes were assembled de novo and annotated for establishing a valid database for studying terpenoid biosynthesis. We identified 135 unigenes that are putatively involved in terpenoid metabolism, including 70 mevalonate and methyl-erythritol phosphate pathways, terpenoid backbone biosynthesis genes, and 65 terpene synthase genes. Moreover, five terpene synthase genes were studied for their functions in terpenoid biosynthesis by using transgenic tobacco; most transgenic tobacco plants expressing these terpene synthetic genes produced increased amounts of terpenoids compared with wild-type control. The combined data analyses from the transcriptome and metabolome provide new insights into our understanding of the complex metabolic genes in terpenoid-rich sage, and our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. officinalis.

Dynamic transcriptomic analysis in hircine longissimus dorsi muscle from fetal to neonatal development stages

Muscle growth and development from fetal to neonatal stages consist of a series of delicately regulated and orchestrated changes in expression of genes. In this study, we performed whole transcriptome profiling based on RNA-Seq of caprine longissimus dorsi muscle tissue obtained from prenatal stages (days 45, 60, and 105 of gestation) and neonatal stage (the 3-day-old newborn) to identify genes that are differentially expressed and investigate their temporal expression profiles. A total of 3276 differentially expressed genes (DEGs) were identified (Q value < 0.01). Time-series expression profile clustering analysis indicated that DEGs were significantly clustered into eight clusters which can be divided into two classes (Q value < 0.05), class I profiles with downregulated patterns and class II profiles with upregulated patterns. Based on cluster analysis, GO enrichment analysis found that 75, 25, and 8 terms to be significantly enriched in biological process (BP), cellular component (CC), and molecular function (MF) categories in class I profiles, while 35, 21, and 8 terms to be significantly enriched in BP, CC, and MF in class II profiles. KEGG pathway analysis revealed that DEGs from class I profiles were significantly enriched in 22 pathways and the most enriched pathway was Rap1 signaling pathway. DEGs from class II profiles were significantly enriched in 17 pathways and the mainly enriched pathway was AMPK signaling pathway. Finally, six selected DEGs from our sequencing results were confirmed by qPCR. Our study provides a comprehensive understanding of the molecular mechanisms during goat skeletal muscle development from fetal to neonatal stages and valuable information for future studies of muscle development in goats.

Transcriptome analysis of Paecilomyces hepiali at different growth stages and culture additives to reveal putative genes in cordycepin biosynthesis

Paecilomyces hepialid (PH) is an endoparasitic fungus of Cordyceps sinensis (CS) and has become a substitute for CS due to their similar pharmacological activities. Because the market demand for CS is difficult to satisfy, and cordycepin, the effective compound of CS, is difficult to industrially produce, we produced 5 samples of PH by culturing for different durations and adding different additives to the culture broth, and detected their cordycepin content with UPLC ESI MS/MS. Then we grouped these cultures into five transcriptome comparisons containing 3 time variable groups and 2 additive variable groups. We used next-generation (NG) sequencing methods to acquire transcriptomic information and investigated the response of gene expression to the additives and the influence of different growth stages. This work will contribute to a better understanding of purine metabolism in PH, and possibly in other Cordyceps species. It will provide a useful resource to further advance transcriptomics studies in Cordyceps species.