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

Transcriptome Analysis of Stephania yunnanensis and Functional Validation of CYP80s Involved in Benzylisoquinoline Alkaloid Biosynthesis

The medicinal plant Stephania yunnanensis is rich in aporphine alkaloids, a type of benzylisoquinoline alkaloid (BIA), with aporphine being the representative and most abundant compound, but our understanding of the biosynthesis of BIAs in this plant has been relatively limited. Previous research reported the genome of S. yunnanensis and preliminarily identified the norcoclaurine synthase (NCS), which is involved in the early stages of the BIA biosynthetic pathways. However, the key genes promoting the formation of the aporphine skeleton have not yet been reported. In this study, based on the differences in the content of crebanine and several other BIAs in different tissues, we conducted transcriptome sequencing of roots, stems, and leaves. We then identified candidate genes through functional annotation and sequence alignment and further analyzed them in combination with the genome. Based on this analysis, we identified three CYP80 enzymes (SyCYP80Q5-1, SyCYP80Q5-3, and SyCYP80G6), which exhibited different activities toward (S)- and (R)-configured substrates in S. yunnanensis and demonstrated strict stereoselectivity enroute to aporphine. This study provides metabolomic and transcriptomic information on the biosynthesis of BIAs in S. yunnanensis, offers valuable insights into the elucidation of BIA biosynthesis, and lays the foundation for the complete analysis of pathways for more aporphine alkaloids.

Transcriptomics and metabolomics revealed the molecular basis of the color formation in the roots of Panax notoginseng

Panax notoginseng is a traditional Chinese medicine rich in many pharmacological components. The root of the 'Miaoxiang Sanqi No. 2' is yellow or greenish yellow, while a novel cultivar-'Wenyuan Ziqi No. 1' shows purple root and is thought to have high medicinal value. Little information is available about the anthocyanin biosynthesis in P. notoginseng root. In this study, we compared the 'Miaoxiang Sanqi No. 2' and 'Wenyuan Ziqi No. 1' in morphological, transcriptional and metabolic levels. The results showed that purple rich in the periderm, rhizome and phloem around cambium of the 'Wenyuan Ziqi No. 1' root and cyanidin 3-O-galactoside was the main anthocyanin causing purple. Moreover, 'Wenyuan Ziqi No. 1' highly accumulated in 155 metabolites, including flavones, phenylpropanoids and lipids. Transcriptome data showed that phenylpropanoid biosynthesis pathway genes are highly expressed in 'Wenyuan Ziqi No. 1'. Conjoint analysis showed that anthocyanin biosynthesis pathway substances were highly accumulated in 'Wenyuan Ziqi No. 1', and the expression level of structural genes involved in anthocyanin biosynthesis pathway was higher in 'Wenyuan Ziqi No. 1'. Meanwhile, eight R2R3-MYB genes that might be involved in anthocyanin biosynthesis were identified. The comprehensive analysis of two cultivars provides new insights into the understanding of root coloration in P. notoginseng.

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.

Characterization of the physiological parameters, effective components, and transcriptional profiles of Polygonum multiflorum Thunb. Under pH stress

Polygonum multiflorum Thunb. is a traditional Chinese medicine with extensive distribution and robust adaptability, but comprehensive research on its acid and alkali resistance is presently lacking. This study aimed to analyze the effects of 5 months of continuous pH stress on the physiological and photosynthetic parameters of P. multiflorum, and the content of effective components. Results revealed that pH stress significantly influenced the normal growth, physiological functions, and photosynthetic indicators of P. multiflorum. At soil pH 4.5, the tubers of P. multiflorum exhibited the highest levels of 2,3,5,4'-tetrahydroxy stilbene-2-O-β-d-glucoside (THSG) and total anthraquinones at 5.41% and 0.38%, respectively. However, increased soil pH significantly reduced the content of THSG and total anthraquinones. Reference-free transcriptome analysis was further conducted on P. multiflorum treated at pH 4.5 and 9.5, generating a total of 47,305 unigenes with an N50 of 2118 bp, of which 31,058 (65.65%) were annotated. Additionally, 2472 differentially expressed genes (DEGs) were identified. Among them, 17 DEGs associated with the biosynthesis of THSG and anthraquinones were screened. A comprehensive analysis of differential gene expression and effective component content demonstrated a significant positive correlation between the content of effective components and the 14 DEGs' expression but a negative correlation with soil pH. This study highlighted the influence of varying soil pH values on the effective component content of P. multiflorum. Specific acidic conditions proved beneficial for the synthesis and accumulation of THSG and total anthraquinones in P. multiflorum, thereby enhancing the quality of the medicinal material.

Quality variation and biosynthesis of anti-inflammatory compounds for Capparis spinosa based on the metabolome and transcriptome analysis

Introduction

Capparis spinosa L. fruits as edible and medicinal plant, has anti-inflammatory activities. The different morphological characteristics of C. spinosa fruits from Ili, Turpan, and Karamay may affect their anti-inflammatory components and functions.

Methods

The anti-inflammatory activity of C. spinosa fruit was assessed using an LPS-induced inflammatory cell model. Furthermore, the differences in anti-inflammatory compounds were analyzed by metabolome and RNA-seq. Additionally, the anti-inflammatory mechanism was elucidated using network pharmacology.

Results

In the study, we found that the 95% ethanol extracts (CSE) obtained from the three kinds of fruits showed remarkable anti-inflammatory effects both in vivo and in vitro. However, the CSE derived from Ili fruits significantly reduced CD86 levels on DCs. As a result of metabolomic analysis, the metabolic profiles of Ili fruits differed significantly from those of the other two habitats, which were consistent with transcriptome analysis. A total of 15 compounds exhibiting anti-inflammatory activity were subjected to screening, revealing a greater accumulation of flavonoids in the Turpan and Karamay districts. Notably, phenolic compounds were identified as the principal anti-inflammatory components in C. spinosa.

Conclusion

There were significant differences in the morphology, metabolites, transcriptional levels, and anti-inflammatory activity of C. spinosa from the three districts.

Chromosome-level genome assembly of Salvia miltiorrhiza with orange roots uncovers the role of Sm2OGD3 in catalyzing 15,16-dehydrogenation of tanshinones

Salvia miltiorrhiza is well known for its clinical practice in treating heart and cardiovascular diseases. Its roots, used for traditional Chinese medicine materials, are usually brick-red due to accumulation of red pigments, such as tanshinone IIA and tanshinone I. Here we report a S. miltiorrhiza line (shh) with orange roots. Compared with the red roots of normal S. miltiorrhiza plants, the contents of tanshinones with a single bond at C-15,16 were increased, whereas those with a double bond at C-15,16 were significantly decreased in shh. We assembled a high-quality chromosome-level genome of shh. Phylogenomic analysis showed that the relationship between two S. miltiorrhiza lines with red roots was closer than the relationship with shh. It indicates that shh could not be the mutant of an extant S. miltiorrhiza line with red roots. Comparative genomic and transcriptomic analyses showed that a 1.0 kb DNA fragment was deleted in shh Sm2OGD3m. Complementation assay showed that overexpression of intact Sm2OGD3 in shh hairy roots recovered furan D-ring tanshinone accumulation. Consistently, in vitro protein assay showed that Sm2OGD3 catalyzed the conversion of cyptotanshinone, 15,16-dihydrotanshinone I and 1,2,15,16-tetrahydrotanshinone I into tanshinone IIA, tanshinone I and 1,2-dihydrotanshinone I, respectively. Thus, Sm2OGD3 functions as tanshinone 15,16-dehydrogenase and is a key enzyme in tanshinone biosynthesis. The results provide novel insights into the metabolic network of medicinally important tanshinone compounds.

Spatiotemporal and Transcriptional Characterization on Tanshinone Initial Synthesis in Salvia miltiorrhiza Roots

Tanshinones are the bioactive constituents of Danshen (Salvia miltiorrhiza Bunge), which is used in Traditional Chinese Medicine to treat cardiovascular and other diseases, and they synthesize and accumulate in the root periderm of S. miltiorrhiza. However, there is no relevant report on the initial stage of tanshinone synthesis, as well as the root structure and gene expression characteristics. The present study aims to provide new insights into how these bioactive principles begin to synthesize by characterizing possible differences in their biosynthesis and accumulation during early root development from both spatial and temporal aspects. The morphological characteristics and the content of tanshinones in roots of S. miltiorrhiza were investigated in detail by monitoring the seedlings within 65 days after germination (DAGs). The ONT transcriptome sequencing was applied to investigate gene expression patterns. The periderm of the S. miltiorrhiza storage taproot initially synthesized tanshinone on about 30 DAGs. Three critical stages of tanshinone synthesis were preliminarily determined: preparation, the initial synthesis, and the continuous rapid synthesis. The difference of taproots in the first two stages was the smallest, and the differentially expressed genes (DEGs) were mainly enriched in terpene synthesis. Most genes involved in tanshinone synthesis were up regulated during the gradual formation of the red taproot. Plant hormone signal transduction and ABC transport pathways were widely involved in S. miltiorrhiza taproot development. Five candidate genes that may participate in or regulate tanshinone synthesis were screened according to the co-expression pattern. Moreover, photosynthetic ferredoxin (FD), cytochrome P450 reductase (CPR), and CCAAT binding transcription factor (CBF) were predicted to interact with the known downstream essential enzyme genes directly. The above results provide a necessary basis for analyzing the initial synthesis and regulation mechanism of Tanshinones.

Plant metabolomics: a new strategy and tool for quality evaluation of Chinese medicinal materials

The present quality control method of Chinese medicinal materials (CMM) has obvious deficiency, which cannot be compatible with the multi-target and multi-component characteristics and production process of CMM. Plant metabolomics with a huge impetus to comprehensively characterize the metabolites and clarify the complexity and integrity of CMM, has been widely used in the research of CMM. This article comprehensively reviewed the application of plant metabolomics in the quality control of CMM. It introduced the concept, technique, and application examples, discussed the prospects, limitations, improvements of plant metabolomics. MS and NMR, as important techniques for plant metabolomics, are mainly highlighted in the case references. The purpose of this article is to clarify the advantage of plants metabolomics for promoting the optimization of the CMM quality control system and proposing a system approach to realize the overall quality control of CMM based on plant metabolomics combined with multidisciplinary method.

Pharmacokinetics, distribution and excretion of inulin-type fructan CPA after oral or intravenous administration to mice

The aim of this work has been to establish and validate a simple and efficient method to detect the concentration of inulin-type fructan CPA from the roots of Codonopsis pilosula (Franch.) Nannf. in biosamples, and then apply it to evaluate the pharmacokinetics behavior, distribution character in tissue and excretion in mice. In this work, fluorescein isothiocyanate (FITC) was used to label CPA. Then FCPA was intravenously and orally administered to mice at different doses. In both i.v and p.o administration, FCPA concentration slowly declined in the circulatory system with a much longer T_1/2 and MRT. After p.o administration, the area under the time curve (AUC0-∞) was dose-dependently increased. Taken together, FCPA showed poor absorption and wide tissue distribution. These pharmacokinetic results yield helpful insights into the pharmacological actions of FCPA.

Integration of Transcriptome and Metabolome Provides New Insights to Flavonoids Biosynthesis in Dendrobium huoshanense

Dendrobium huoshanense is both a traditional herbal medicine and a plant of high ornamental and medicinal value. We used transcriptomics and metabolomics to investigate the effects of growth year on the secondary metabolites of D. huoshanense stems obtained from four different years of cultivation. In this study, a total of 428 differentially accumulated metabolites (DAMs) and 1802 differentially expressed genes (DEGs) were identified. The KEGG enrichment analysis of DEGs and DAMs revealed significant differences in "Flavonoid biosynthesis", "Phenylpropanoid biosynthesis" and "Flavone and flavonol biosynthesis". We summarize the biosynthesis pathway of flavonoids in D. huoshanense, providing new insights into the biosynthesis and regulation mechanisms of flavonoids in D. huoshanense. Additionally, we identified two candidate genes, FLS (LOC110107557) and F3'H (LOC110095936), which are highly involved in flavonoid biosynthesis pathway, by WGCNA analysis. The aim of this study is to investigate the effects of growth year on secondarily metabolites in the plant and provide a theoretical basis for determining a reasonable harvesting period for D. huoshanense.

Comparative transcriptome analyses of different Salvia miltiorrhiza varieties during the accumulation of tanshinones

Salvia miltiorrhiza (Labiatae) is an important medicinal plant in traditional Chinese medicine. Tanshinones are one of the main active components of S. miltiorrhiza. It has been found that the intraspecific variation of S. miltiorrhiza is relatively large and the content of tanshinones in its roots of different varieties is also relatively different. To investigate the molecular mechanisms that responsible for the differences among these varieties, the tanshinones content was determined and comparative transcriptomics analysis was carried out during the tanshinones accumulation stage. A total of 52,216 unigenes were obtained from the transcriptome by RNA sequencing among which 23,369 genes were differentially expressed among different varieties, and 2,016 genes including 18 diterpenoid biosynthesis-related genes were differentially expressed during the tanshinones accumulation stage. Functional categorization of the differentially expressed genes (DEGs) among these varieties revealed that the pathway related to photosynthesis, oxidative phosphorylation, secondary metabolite biosynthesis, diterpenoid biosynthesis, terpenoid backbone biosynthesis, sesquiterpenoid and triterpenoid biosynthesis are the most differentially regulated processes in these varieties. The six tanshinone components in these varieties showed different dynamic changes in tanshinone accumulation stage. In addition, combined with the analysis of the dynamic changes, 277 DEGs (including one dehydrogenase, three CYP450 and 24 transcription factors belonging to 12 transcription factor families) related to the accumulation of tanshinones components were obtained. Furthermore, the KEGG pathway enrichment analysis of these 277 DEGs suggested that there might be an interconnection between the primary metabolic processes, signaling processes and the accumulation of tanshinones components. This study expands the vision of intraspecific variation and gene regulation mechanism of secondary metabolite biosynthesis pathways in medicinal plants from the “omics” perspective.

Permissive action of H2O2 mediated ClUGT75 expression for auxin glycosylation and Al3+- tolerance in watermelon

Although Al³⁺-toxicity is one of the limiting factors for crop production in acidic soils, little is known about the Al³⁺-tolerance mechanism in watermelon, a fairly acid-tolerant crop. This work aimed to identify the interaction between the H₂O₂ scavenging pathway and auxin glycosylation relevant to watermelon Al³⁺-tolerance. By analyzing expressions of hormone-related ClUGTs and antioxidant enzyme genes in Al³⁺-tolerant (ZJ) and Al³⁺-sensitive (NBT) cultivars, we identified ClUGT75s (B1, B2, and D1) and ClSOD1-2-ClCAT as crucial components associated with Al³⁺-tolerance. Al³⁺-stress significantly increased H₂O₂ content by 92.7% in NBT and 42.3% in ZJ, accompanied by less Al³⁺-, auxin (IAA and IBA), and MDA contents in ZJ than NBT. These findings coincided with significant ClSOD1-2 expression and stable dismutation activity in NBT than ZJ. Hence, higher H₂O₂ content in the root apex of NBT than ZJ correlated with a significant increase in auxin content and ClSOD1-2 up-regulation. Moreover, Al³⁺-activated ClUGT75D1 and ClUGT75B2 in ZJ coincided with no considerable change in IBA content, suggesting that glycosylation-mediated changes in IBA content might be relevant to Al³⁺-tolerance in watermelon. Furthermore, exogenous H₂O₂ and IBA indicated ClUGT75D1 modulating IBA is likely dependent on H₂O₂ background. We hypothesize that a higher H₂O₂ level in NBT represses ClUGT75, resulting in increased auxin than those in ZJ roots. Thus, excess in both H₂O₂ and auxin aggravated the inhibition of root elongation under Al³⁺-stress. Our findings provide insights on the permissive action of H₂O₂ in the mediation of auxin glycosylation by ClUGT75 in root apex for Al³⁺-tolerance in watermelon.

Comparative transcriptome analysis reveals key insights into male sterility in Salvia miltiorrhiza Bunge

Background

Large-scale heterosis breeding depends upon stable, inherited male sterility lines. We accidentally discovered a male sterility line (SW-S) in the F₁progeny of a Salvia miltiorrhiza Bunge from Shandong, China (purple flowers) crossed with a S. miltiorrhiza f. alba from Sichuan, China (white flowers). We sought to provide insights into the pollen development for male sterility in S. miltiorrhiza.

Methods

The phenotypic and cytological features of the SW-S and fertile control SW-F were observed using scanning electron microscopy and paraffin sections to identify the key stage of male sterility. Transcriptome profiles were recorded for anthers at the tetrad stage of SW-S and SW-F using Illumina RNA-Seq.

Results

The paraffin sections showed that sterility mainly occurred at the tetrad stage of microspore development, during which the tapetum cells in the anther compartment completely fell off and gradually degraded in the sterile line. There was little-to-no callose deposited around the microspore cells. The tetrad microspore was shriveled and had abnormal morphology. Therefore, anthers at the tetrad stage of SW-S and fertile control SW-F were selected for comparative transcriptome analysis. In total, 266,722,270 clean reads were obtained from SW-S and SW-F, which contained 36,534 genes. There were 2,571 differentially expressed genes (DEGs) in SW-S and SW-F, of which 63.5% were downregulated. Gene Ontology (GO) enrichment analysis indicated that the differentially expressed genes were enriched in 56 functional groups (GO terms); of these, all DEGs involved in microgametogenesis and developmental maturation were downregulated in SW-S. These results were confirmed by quantitative RT-PCR. The two GO terms contained 18 DEGs, among which eight DEGs (namely: GPAT3, RHF1A, phosphatidylinositol, PFAS, MYB96, MYB78, Cals5, and LAT52) were related to gamete development. There were 10 DEGs related to development and maturation, among which three genes were directly related to pollen development (namely: ACT3, RPK2, and DRP1C). Therefore, we believe that these genes are directly or indirectly involved in the pollen abortion of SW-S. Our study provides insight into key genes related to sterility traits in S. miltiorrhiza, and the results can be further exploited in functional and mechanism studies.

Strepimidazoles A-G from the Plant Endophytic Streptomyces sp. PKU-EA00015 with Inhibitory Activities against a Plant Pathogenic Fungus

Seven new 4-acyl-2-aminoimidazoles, designated strepimidazoles A-G (1-7), were discovered from the endophytic Streptomyces sp. PKU-EA00015 isolated from Salvia miltiorrhiza Bunge, whose dry root "Danshen" is one of the most widely used traditional Chinese medicines. The resonance signals of the 2-aminoimidazole moiety in 1-7 were absent in the NMR spectra due to tautomerization, and the structures of 1-7 were identified after preparation of their acetylation products 1a-7a, respectively. Compounds 1-7 represent a new family of 2-aminoimidazole-containing natural products, enriching the structural diversity of natural products from endophytic origin. Compounds 1-7 showed different degrees of inhibitory activities against the plant pathogenic fungus Verticillium dahliae V991, revealing structure-activity relationships on the acyl moieties. The plant pathogenic fungus V. dahliae has been confirmed to cause serious chlorosis of cultivated S. miltiorrhiza Bunge in China. This study opens the door for further investigation of mutualistic relationships between S. miltiorrhiza Bunge and their endophytic actinomycetes and for possible antifungal agent development for biological control of V. dahliae in the future.

LC-MS-based multivariate statistical analysis for the screening of potential thrombin/factor Xa inhibitors from Radix Salvia Miltiorrhiza

BACKGROUND

The dry root and rhizome of Salvia miltiorrhiza Bunge, or Danshen, is a well-known traditional Chinese medicine with anticoagulant activity. Taking into account that thrombin (THR) and factor Xa (FXa) play crucial roles in the coagulation cascade, it is reasonable and meaningful to screening THR and/or FXa inhibitors from Danshen.

METHODS

Four extracts [butanol (BA), ethyl acetate (EA) and remained extract (RE) from 75% ethanol extract, and water extract (WE)] of Danshen were prepared, and their THR/FXa inhibitory activities were assessed in vitro. Then, the active EA extract was further separated by silica-gel column chromatography (SC), and its fractions (SC1-SC5) were analyzed by LC-MS. The principal component analysis (PCA) and orthogonal partial least squares discriminate analysis (OPLS-DA) were employed for predicting the specific marker compounds. The chemical structures of targeted compounds were identified by LC-MS/MS and their interactions with THR/FXa were analyzed by the molecular docking analysis.

RESULTS

Danshen EA extract showed strong activity against THR and FXa, and its fractions (SC1-SC5) exhibited obvious difference in inhibitory activity against these two enzymes. Furthermore, four marker compounds with potential THR/FXa inhibitory activity were screened by PCA and OPLS-DA, and were identified as cryptotanshinone, tanshinone I, dihydrotanshinone I and tanshinone IIA. The molecular docking study showed that all these four tanshinones can interact with some key amino acid residues of the THR/FXa active cavities, such as HIS57 and SER195, which were considered to be promising candidates targeting THR and/or FXa with low binding energy (< - 7 kcal mol-1).

CONCLUSIONS

LC-MS combined with multivariate statistical analysis can effectively screen potential THR/FXa inhibitory components in Danshen.

Transcriptomic and Metabolomic Changes Triggered by Fusarium solani in Common Bean (Phaseolus vulgaris L.)

Common bean (Phaseolus vulgaris L.) is a major legume and is frequently attacked by fungal pathogens, including Fusarium solani f. sp. phaseoli (FSP), which cause Fusarium root rot. FSP substantially reduces common bean yields across the world, including China, but little is known about how common bean plants defend themselves against this fungal pathogen. In the current study, we combined next-generation RNA sequencing and metabolomics techniques to investigate the changes in gene expression and metabolomic processes in common bean infected with FSP. There were 29,722 differentially regulated genes and 300 differentially regulated metabolites between control and infected plants. The combined omics approach revealed that FSP is perceived by PAMP-triggered immunity and effector-triggered immunity. Infected seedlings showed that common bean responded by cell wall modification, ROS generation, and a synergistic hormone-driven defense response. Further analysis showed that FSP induced energy metabolism, nitrogen mobilization, accumulation of sugars, and arginine and proline metabolism. Importantly, metabolic pathways were most significantly enriched, which resulted in increased levels of metabolites that were involved in the plant defense response. A correspondence between the transcript pattern and metabolite profile was observed in the discussed pathways. The combined omics approach enhances our understanding of the less explored pathosystem and will provide clues for the development of common bean cultivars' resistant to FSP.