De novo leaf and root transcriptome analysis identified novel genes involved in steroidal sapogenin biosynthesis in Asparagus racemosus

Arbuscular mycorrhiza and plant growth promoting endophytes facilitates accumulation of saponin under moderate drought stress

Objective

Paris polyphylla var. yunnanensis, one of the important medicinal plant resources in Yunnan, China, usually takes 6-8 years to be harvested. Therefore, it is urgent to find a method that can not only shorten its growth years, but also improve its quality. In this study, we examined the effects of a combination treatment of arbuscular mycorrhizal fungi (AMF) and plant growth-promoting endophytes (PGPE) and drought stress on the accumulation of saponins in it.

Methods

P. polyphylla var. yunnanensis was infected with a mixture of AMF and PGPE under drought stress. The content of saponins, as well as morphological, physiological, and biochemical indicators, were all measured. The UGTs gene related to saponin synthesis was obtained from transcriptome data by homologous comparison, which were used for RT-PCR and phylogenetic analysis.

Results

Regardless of water, AMF treatment could infect the roots of P. polyphylla var. yunnanensis, however double inoculation with AMF and PGPE (AMF + PGPE) would reduce the infection rate of AMF. Plant height, aboveground and underground fresh weight did not differ significantly between the single inoculation AMF and the double inoculation treatment under different water conditions, but the inoculation treatment significantly increased the plant height of P. polyphylla var. yunnanensis compared to the non-inoculation treatment. Single inoculation with AMF considerably increased the net photosynthetic rate, stomatal conductance, and transpiration rate of P. polyphylla var. yunnanensis leaves under various water conditions, but double inoculation with AMF + PGPE greatly increased the intercellular CO2 concentration and chlorophyll fluorescence parameter (Fv/Fm). Under diverse water treatments, single inoculation AMF had the highest proline content, whereas double inoculation AMF + PGPE may greatly improve the amount of abscisic acid (ABA) and indoleacetic acid (IAA) compared to normal water under moderate drought. Double inoculation AMF + PGPE treatment improved the proportion of N, P, and K in the rhizome of P. polyphylla var. yunnanensis under various water conditions. Under moderate drought stress, AMF + PGPE significantly enhanced the contents of P. polyphylla var. yunnanensis saponins I, II, VII, and total saponins as compared to normal water circumstances. Farnesyl diphosphate synthase (FPPS), Geranyl pyrophosphate synthase (GPPS), Cycloartenol synthase (CAS), and Squalene epoxidase (SE1) were the genes that were significantly up-regulated at the same time. The amount of saponins was favorably linked with the expression of CAS, GPPS, and SE1. Saponin VI content and glycosyl transferase (UGT) 010922 gene expression were found to be substantially associated, as was saponin II content and UGT010935 gene expression.

Conclusion

Under moderate drought, AMF + PGPE was more conducive to the increase of hormone content, nutrient absorption, and total saponin content in P. polyphylla var. yunnanensis, and AMF + PGPE could up regulate the expression of key genes and UGTs genes in one or more steroidal saponin synthesis pathways to varying degrees, thereby stimulating the synthesis and accumulation of steroidal saponins in the rhizome of P. polyphylla var. yunnanensis. The combination of AMF and PGPE inoculation, as well as adequate soil drought, reduced the buildup of saponins in P. polyphylla var. yunnanensis and increased its quality.

Transcriptome analysis of three medicinal plants of the genus Polygonatum: identification of genes involved in polysaccharide and steroidal saponins biosynthesis

Polysaccharides and saponins are the main active components of Polygonati Rhizoma. Studying the molecular mechanism of their synthesis pathway is helpful in improving the content of active components at the molecular level. At present, transcriptome analysis of three Polygonatum species (Polygonatum sibiricum Red., Polygonatum cyrtonema Hua, Polygonatum kingianum Coll. et Hemsl.) has been reported, but no comparative study has been found on the transcriptome data of the three species. Transcriptome sequencing was performed on the rhizomes of three Polygonatum species based on high-throughput sequencing technology, and all transcripts were assembled. A total of 168,108 unigenes were generated after the removal of redundancy, of which 121,642 were annotated in seven databases. Through differential analysis and expression analysis of key enzyme genes in the synthesis pathway of three Polygonatum polysaccharides and steroidal saponins, 135 differentially expressed genes encoding 18 enzymes and 128 differentially expressed genes encoding 28 enzymes were identified, respectively. Numerous transcription factors are involved in the carbohydrate synthesis pathway. Quantitative real-time PCR was used to further verify the gene expression level. In this paper, we present a public transcriptome dataset of three medicinal plants of the genus Polygonatum, and analyze the key enzyme genes of polysaccharide and steroidal saponins synthesis pathway, which lays a foundation for improving the active component content of Polygonati Rhizoma by molecular means.

Integrated evolutionary pattern analyses reveal multiple origins of steroidal saponins in plants

Steroidal saponins are a class of specialized metabolites essential for plant's response to biotic and abiotic stresses. They are also important raw materials for the industrial production of steroid drugs. Steroidal saponins are present in some monocots, such as Dioscorea and Paris, but their distribution, origin, and evolution in plants remain poorly understood. By reconstructing the evolutionary history of the steroidal saponin-associated module (SSAM) in plants, we reveal that the steroidal saponin pathway has its origin in Asparagus and Dioscorea. Through evaluating the distribution and evolutionary pattern of steroidal saponins in angiosperms, we further show that steroidal saponins originated multiple times in angiosperms, and exist in early diverged lineages of certain monocot lineages including Asparagales, Dioscoreales, and Liliales. In these lineages, steroidal saponins are synthesized through the high copy and/or high expression mechanisms of key genes in SSAM. Together with shifts in gene evolutionary rates and amino acid usage, these molecular mechanisms shape the current distribution and diversity of steroidal saponins in plants. Consequently, our results provide new insights into the distribution, diversity and evolutionary history of steroidal saponins in plants, and enhance our understanding of plants' resistance to abiotic and biotic stresses. Additionally, fundamental understanding of the steroidal saponin biosynthesis will facilitate their industrial production and pharmacological applications.

De novo sequencing, assembly, and characterization of Asparagus racemosus transcriptome and analysis of expression profile of genes involved in the flavonoid biosynthesis pathway

Asparagus racemosus is known for its diverse content of secondary metabolites, i.e., saponins, alkaloids, and a wide range of flavonoids. Flavonoids, including phenols and polyphenols, have a significant role in plant physiology and are synthesized in several tissues. Despite the diverse role of flavonoids, genetic information is limited for flavonoid biosynthesis pathways in A. racemosus. The current study explores full-scale functional genomics information of A. racemosus by de novo transcriptome sequencing using Illumina paired-end sequencing technology to elucidate the genes involved in flavonoid biosynthesis pathways. The de novo assembly of high-quality paired-end reads resulted in ∼2.3 million high-quality reads with a pooled transcript of 45,647 comprising ∼76 Mb transcriptome with a mean length (bp) of 1,674 and N50 of 1,868bp. Furthermore, the coding sequence (CDS) prediction analysis from 45,647 pooled transcripts resulted in 45,444 CDS with a total length and mean length of 76,398,686 and 1,674, respectively. The Gene Ontology (GO) analysis resulted in a high number of CDSs assigned to 25,342 GO terms, which grouped the predicted CDS into three main domains, i.e., Biological Process (19,550), Molecular Function (19,873), and Cellular Component (14,577). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database was used to categorize 6,353 CDS into 25 distinct biological pathway categories, in which the majority of mapped CDS were shown to be related to translation (645), followed by signal transduction (532), carbohydrate metabolism (524), folding, sorting, and degradation (522). Among these, only ∼64 and 14 CDSs were found to be involved in the phenylpropanoid and flavonoid biosynthesis pathways, respectively. Quantitative Real-time PCR was used to check the expression profile of fourteen potential flavonoid biosynthesis pathway genes. The qRT-PCR analysis result matches the transcriptome sequence data validating the Illumina sequence results. Moreover, a large number of genes associated with the flavonoids biosynthesis pathway were found to be upregulated under the induction of methyl jasmonate. The present-day study on transcriptome sequence data of A. racemosus can be utilized for characterizing genes involved in flavonoid biosynthesis pathways and for functional genomics analysis in A. racemosus using the reverse genetics approach (CRISPR/Cas9 technology).

Identification of key genes involved in secondary metabolite biosynthesis in Digitalis purpurea

The medicinal plant Digitalis purpurea produces cardiac glycosides that are useful in the pharmaceutical industry. These bioactive compounds are in high demand due to ethnobotany's application to therapeutic procedures. Recent studies have investigated the role of integrative analysis of multi-omics data in understanding cellular metabolic status through systems metabolic engineering approach, as well as its application to genetically engineering metabolic pathways. In spite of numerous omics experiments, most molecular mechanisms involved in metabolic pathways biosynthesis in D. purpurea remain unclear. Using R Package Weighted Gene Co-expression Network Analysis, co-expression analysis was performed on the transcriptome and metabolome data. As a result of our study, we identified transcription factors, transcriptional regulators, protein kinases, transporters, non-coding RNAs, and hub genes that are involved in the production of secondary metabolites. Since jasmonates are involved in the biosynthesis of cardiac glycosides, the candidate genes for Scarecrow-Like Protein 14 (SCL14), Delta24-sterol reductase (DWF1), HYDRA1 (HYD1), and Jasmonate-ZIM domain3 (JAZ3) were validated under methyl jasmonate treatment (MeJA, 100 μM). Despite early induction of JAZ3, which affected downstream genes, it was dramatically suppressed after 48 hours. SCL14, which targets DWF1, and HYD1, which induces cholesterol and cardiac glycoside biosynthesis, were both promoted. The correlation between key genes and main metabolites and validation of expression patterns provide a unique insight into the biosynthesis mechanisms of cardiac glycosides in D. purpurea.

Next-generation sequencing in the biodiversity conservation of endangered medicinal plants

Medicinal plants have been used as traditional herbal medicines in the treatment of various types of diseases. However, the increased demand for these plants highlights the importance of conservation specifically for endangered species. Significant advancements in next-generation sequencing (NGS) technologies have accelerated medicinal plant research while reducing costs and time demands. NGS systems enable high-throughput whole genome sequencing as well as direct RNA sequencing and transcriptome analysis. The sequence data sets created can be used in a variety of areas of study, including biodiversity conservation, comparative genomics, transcriptomic analysis, single cell mining, metagenomics, epigenetics, molecular marker discovery, multi genome sequencing, and so on. Commercial sequencing service providers are constantly working to improve technologies to address bioinformatics problems in NGS data analysis. Several genome sequencing projects on medicinal plants have been completed recently and a few more are in the works. In some medicinal plants, massive NGS-based data has been developed. In the present review, we have attempted to briefly discuss advancements in NGS technology on medicinally essential plants in India. The review will also provide ideas for applying NGS technologies for exploring genomes of various endangered medicinal plants whose genome sequences are not normally available and thus provides valuable insights for the conservation of these vulnerable species.

Analysis of the mechanism of saponin biotransformation in Dioscoreae nipponicae rhizoma by the endophytic fungus Fusarium sp. C39 using whole-genome sequencing

Fusarium sp. C39 is an endophytic fungus of Dioscorea nipponica Makino. Symbiosis of Fusarium sp. C39 with Rhizoma Dioscoreae Nipponicae (RDN) can significantly increase the content of saponin, which provides a new approach for saponin production and reduces the pressure on natural sources of saponins. However, the underlying mechanism is not clear, limiting its application. In this study, the genome of Fusarium sp. C39 was sequenced, the gene functions were predicted via gene annotation, and the genome was compared to the genomes of four related species. Fusarium sp. C39 is predicted to encode many key enzyme genes involved in saponin synthesis, which could transform the mevalonate, isopentenyl pyrophosphate, and various intermediate compounds present in the RDN extract into saponins. The Fusarium sp. C39 genome contains specific genes that are conducive to its endophytic lifestyle and can provide abundant raw materials for saponin synthesis. Based on the genomic analysis, we proposed the mechanism by which Fusarium sp. C39 generates saponins and provides a theoretical basis for rapid, efficient, low-cost production of saponins.

Melatonin confers fenugreek tolerance to salinity stress by stimulating the biosynthesis processes of enzymatic, non-enzymatic antioxidants, and diosgenin content

Salinity-induced stress is widely considered a main plant-growth-limiting factor. The positive effects of melatonin in modulating abiotic stresses have led this hormone to be referred to as a growth regulator in plants. This study aims to show how melatonin protects fenugreek against the negative effects of salt stress. Different amounts of melatonin (30, 60, and 90 ppm), salinity stress (150 mM and 300 mM), and the use of both salinity and melatonin were used as treatments. The results showed that applying different melatonin levels to salinity-treated fenugreek plants effectively prevented the degradation of chlorophyll a, chlorophyll b, total chlorophyll, and carotenoid contents compared with salinity treatment without melatonin application. Besides, melatonin increases the biosynthesis of enzymatic and non-enzymatic antioxidants, thereby adjusting the content of reactive oxygen species, free radicals, electrolyte leakage, and malondialdehyde content. It was observed that applying melatonin increased the activity of potassium-carrying channels leading to the maintenance of ionic homeostasis and increased intracellular water content under salinity stress. The results revealed that melatonin activates the defense signaling pathways in fenugreek through the nitric oxide, auxin, and abscisic acid-dependent pathways. Melatonin, in a similar vein, increased the expression of genes involved in the biosynthesis pathway of diosgenin, a highly important steroidal sapogenin in medical and food industries, and hence the diosgenin content. When 150 mM salinity stress and 60 ppm melatonin were coupled, the diosgenin concentration rose by more than 5.5 times compared to the control condition. In conclusion, our findings demonstrate the potential of melatonin to enhance the plant tolerance to salinity stress by stimulating biochemical and physiological changes.

Altering Sterol Composition Implied That Cholesterol Is Not Physiologically Associated With Diosgenin Biosynthesis in Trigonella foenum-graecum

Diosgenin serves as an important precursor of most steroidal drugs in market. Cholesterol was previously deemed as a sterol origin leading to diosgenin biosynthesis. This study reports that cholesterol is not in parallel with diosgenin biosynthesis in Trigonella foenum-graecum. We first perturbed its sterol composition using inhibitors specific for the upstream isoprenoid pathway enzymes, HMGR (3-hydroxy-3-methylgutaryl-CoA reductase) on the mevalonate (MVA) and DXR (1-deoxy-D-xylulose-5-phosphate reductoisomerase) on the 2-C-methyl-D-erythritol-4-phophate (MEP) pathways, and have revealed that diosgenin and cholesterol reversely or differently accumulated in either the MVA or the MEP pathway-suppressed plants, challenging the previously proposed role of cholesterol in diosgenin biosynthesis. To further investigate this, we altered the sterol composition by suppressing and overexpressing the 24-sterol methyltransferase type 1 (SMT1) gene in T. foenum-graecum, as SMT1 acts in the first committed step of diverting the carbon flux of cholesterol toward biosynthesis of 24-alkyl sterols. Knockdown of TfSMT1 expression led to increased cholesterol level but caused a large reduction of diosgenin. Diosgenin was increased upon the TfSMT1-overexpressing, which, however, did not significantly affect cholesterol biosynthesis. These data consistently supported that diosgenin biosynthesis in T. foenum-graecum is not associated with cholesterol. Rather, campesterol, a 24-alkyl sterol, was indicative of being correlative to diosgenin biosynthesis in T. foenum-graecum.

Advances in steroidal saponins biosynthesis

This work reviews recent advances in the pathways and key enzymes of steroidal saponins biosynthesis and sets the foundation for the biotechnological production of these useful compounds through transformation of microorganisms. Steroidal saponins, due to their specific chemical structures and active effects, have long been important natural products and that are irreplaceable in hormone production and other pharmaceutical industries. This article comprehensively reviewed the previous and current research progress and summarized the biosynthesis pathways and key biosynthetic enzymes of steroidal saponins that have been discovered in plants and microoganisms. On the basis of the general biosynthetic pathway in plants, it was found that the starting components, intermediates and catalysing enzymes were diverse between plants and microorganisms; however, the functions of their related enzymes tended to be similar. The biosynthesis pathways of steroidal saponins in microorganisms and marine organisms have not been revealed as clearly as those in plants and need further investigation. The elucidation of biosynthetic pathways and key enzymes is essential for understanding the synthetic mechanisms of these compounds and provides researchers with important information to further develop and implement the massive production of steroidal saponins by biotechnological approaches and methodologies.

The Genetic Regulation of Secondary Metabolic Pathways in Response to Salinity and Drought as Abiotic Stresses

Global development has generated a plethora of unfavorable and adverse environmental factors for the living organisms in the ecosystem. Plants are sessile organisms, and they are crucial to sustain life on earth. Since plants are sessile, they face a great number of environmental challenges related to abiotic stresses, such as temperature fluctuation, drought, salinity, flood and metal contamination. Salinity and drought are considered major abiotic stresses that negatively affect the plants’ growth and production of useful content. However, plants have evolved various molecular mechanisms to increase their tolerance to these environmental stresses. There is a whole complex system of communication (cross-talk) through massive signaling cascades that are activated and modulated in response to salinity and drought. Secondary metabolites are believed to play significant roles in the plant’s response and resistance to salinity and drought stress. Until recently, attempts to unravel the biosynthetic pathways were limited mainly due to the inadequate plant genomics resources. However, recent advancements in generating high-throughput “omics” datasets, computational tools and functional genomics approach integration have aided in the elucidation of biosynthetic pathways of many plant bioactive metabolites. This review gathers comprehensive knowledge of plants’ complex system that is involved in the response and resistance to salinity and water deficit stresses as abiotic stress. Additionally, it offers clues in determining the genes involved in this complex and measures its activity. It covers basic information regarding the signaling molecules involved in salinity and drought resistance and how plant hormones regulate the cross-talking mechanism with emphasis on transcriptional activity. Moreover, it discusses many studies that illustrate the relationship between salinity and drought and secondary metabolite production. Furthermore, several transcriptome analysis research papers of medicinal plants are illustrated. The aim of this review is to be a key for any researcher that is aspiring to study the relationship between salinity and drought stresses and secondary metabolite production at the transcriptome and transcription level.

Assessing and integrating the transcriptome analysis with plant development, trichomes, and secondary metabolites yield potential in Mentha arvensisL

Mentha arvensisL. (corn mint) is well known for the production of menthol, a used commodity in flavouring industries, and provides natural fragrances. Glandular trichomes are responsible for producing specific secondary metabolites in vascular plants having species chemistry. Ten cultivars/varieties of M. arvensis, namely, Saksham, Kosi, Himalaya, Gomti, Sambhav, Kalka, Damroo, Kushal, and Shivalik, were used to study the developmental regulation of trichomes, essential oil yield, chemical constituents of essential oil and morphological parameters were estimated with gene expression using a randomized block design. Simultaneously, RNA sequence-based transcriptome analysis was done to reveal the transcription factors and differential gene analysis, which are responsible for the biosynthesis of essential oil as well as trichome development. Plant growth showed the maximum transition between 35 and 50 days stage, while essential oil and its metabolite bioconversion was observed in between 70 and 100 days stage. Glandular trichomes were maximally increased between 50, 70 days, and 100 days stage in var. Kosi followed by var. Saryu which has rapid growth in oil content. Menthol reductase activity was found to be a regulatory element during development, as it follows the inverse trend of menthol content and leads to menthol accumulation in subcuticular spaces. Transcriptional factors, cog, and nonredundant novel genes were identified. The composition of mintessential oils is regulated at multiple levels, including transcript abundance, catalytic properties of enzyme catalysts, and cell type-specific epigenetic processes.

Tissue-specific signatures of metabolites and proteins in asparagus roots and exudates

Comprehensive untargeted and targeted analysis of root exudate composition has advanced our understanding of rhizosphere processes. However, little is known about exudate spatial distribution and regulation. We studied the specific metabolite signatures of asparagus root exudates, root outer (epidermis and exodermis), and root inner tissues (cortex and vasculature). The greatest differences were found between exudates and root tissues. In total, 263 non-redundant metabolites were identified as significantly differentially abundant between the three root fractions, with the majority being enriched in the root exudate and/or outer tissue and annotated as 'lipids and lipid-like molecules' or 'phenylpropanoids and polyketides'. Spatial distribution was verified for three selected compounds using MALDI-TOF mass spectrometry imaging. Tissue-specific proteome analysis related root tissue-specific metabolite distributions and rhizodeposition with underlying biosynthetic pathways and transport mechanisms. The proteomes of root outer and inner tissues were spatially very distinct, in agreement with the fundamental differences between their functions and structures. According to KEGG pathway analysis, the outer tissue proteome was characterized by a high abundance of proteins related to 'lipid metabolism', 'biosynthesis of other secondary metabolites' and 'transport and catabolism', reflecting its main functions of providing a hydrophobic barrier, secreting secondary metabolites, and mediating water and nutrient uptake. Proteins more abundant in the inner tissue related to 'transcription', 'translation' and 'folding, sorting and degradation', in accord with the high activity of cortical and vasculature cell layers in growth- and development-related processes. In summary, asparagus root fractions accumulate specific metabolites. This expands our knowledge of tissue-specific plant cell function.

The Plausible Role of Indian Traditional Medicine in Combating Corona Virus (SARS-CoV 2): A Mini-Review

SARS-CoV 2 is a novel virus strain of Coronavirus, reported in China in late December 2019. Its highly contagious nature in humans has prompted WHO to designate the ongoing pandemic as a Public Health Emergency of International Concern. At this moment, there is no specific treatment and the therapeutic strategies to deal with the infection are only supportive, with prevention aimed at reducing community transmission. A permanent solution for the pandemic, which has brought the world economy to the edge of collapse, is the need of the hour. This situation has brought intense research in traditional systems of medicine. Indian Traditional System, Ayurveda, has a clear concept of the cause and treatment of pandemics. Through this review, information on the potential antiviral traditional medicines along with their immunomodulatory pathways are discussed. We have covered the seven most important Indian traditional plants with antiviral properties: Withania somnifera (L.) Dunal (family: Solanaceae), Tinospora cordifolia (Thunb.) Miers (family: Menispermaceae), Phyllanthus emblica L. (family: Euphorbiaceae), Asparagus racemosus L. (family: Liliaceae), Glycyrrhiza glabra L. (family: Fabaceae), Ocimum sanctum L. (family: Lamiaceae) and Azadirachta indica A. Juss (family: Meliaceae) in this review. An attempt is also made to bring into limelight the importance of dietary polyphenol, Quercetin, which is a potential drug candidate in the making against the SARS-CoV2 virus.

Transcriptome analysis of Paris polyphylla var. yunnanensis illuminates the biosynthesis and accumulation of steroidal saponins in rhizomes and leaves

Paris polyphylla var. yunnanensis can synthesize Paris saponins with multiple effective therapies, and its rhizome has become an indispensable ingredient in many patented drugs. However, how Paris saponin content changes in tissues at different stages and the molecular mechanisms underlying the production and accumulation of the bioactive compounds are unclear. This study aimed to uncover the mechanisms underlying the biosynthesis and accumulation by integrating transcriptome sequencing and phytochemical investigation of the leaves and rhizomes at different growth stages. Paris saponin content in leaves was lower during the fruiting stage than the vegetative stage, whereas the content in rhizomes increased during the fruiting stage. The candidate genes related to Paris saponin biosynthesis were determined by transcriptome analyses. Most biosynthetic genes were found to be abundantly expressed in the leaves during the vegetative stage in the light of expression profiles and functional enrichment results. The expression patterns of the differentially expressed genes related to the biosynthesis were positively correlated with the accumulation of saponins in tissues. These findings suggest that both leaves and rhizomes are capable of biosynthesizing Paris saponins, and that aerial plant parts can be used to extract them. The different patterns of biosynthesis and accumulation in the leaves and rhizomes were also determined here. This study will help improve our understanding of the mechanisms underlying the biosynthesis and accumulation of Paris saponins, and aid in the comprehensive development and utilization of this medicinal plant.

Comparative transcriptome analysis of root, stem, and leaf tissues of Entada phaseoloides reveals potential genes involved in triterpenoid saponin biosynthesis

Background

Entada phaseoloides (L.) Merr. is an important traditional medicinal plant. The stem of Entada phaseoloides is popularly used as traditional medicine because of its significance in dispelling wind and dampness and remarkable anti-inflammatory activities. Triterpenoid saponins are the major bioactive compounds of Entada phaseoloides. However, genomic or transcriptomic technologies have not been used to study the triterpenoid saponin biosynthetic pathway in this plant.

Results

We performed comparative transcriptome analysis of the root, stem, and leaf tissues of Entada phaseoloides with three independent biological replicates and obtained a total of 53.26 Gb clean data and 116,910 unigenes, with an average N50 length of 1218 bp. Putative functions could be annotated to 42,191 unigenes (36.1%) based on BLASTx searches against the Non-redundant, Uniprot, KEGG, Pfam, GO, KEGG and COG databases. Most of the unigenes related to triterpenoid saponin backbone biosynthesis were specifically upregulated in the stem. A total of 26 cytochrome P450 and 17 uridine diphosphate glycosyltransferase candidate genes related to triterpenoid saponin biosynthesis were identified. The differential expressions of selected genes were further verified by qPT-PCR.

Conclusions

The dataset reported here will facilitate the research about the functional genomics of triterpenoid saponin biosynthesis and genetic engineering of Entada phaseoloides.

De novo transcriptome analysis unravels tissue-specific expression of candidate genes involved in major secondary metabolite biosynthetic pathways of Plumbago zeylanica: implication for pharmacological potential

KEY MESSAGE

The present study provides comparative transcriptome analysis, besides identifying functional secondary metabolite genes of Plumbago zeylanica with pharmacological potential for future functional genomics, and metabolomic engineering of secondary metabolites from this plant towards diversified biomedical applications.

ABSTRACT

Plumbago zeylanica is a widely used medicinal plant of the traditional Indian system of medicine with wide pharmacological potential to treat several disorders. The present study aimed to carry out comparative transcriptome analysis in leaf and root tissue of P. zeylanica using Illumina paired end sequencing to identify tissue-specific functional genes involved in the biosynthesis of secondary metabolites, contributing to its therapeutic efficacy. De novo sequencing assembly resulted in the identification of 62,321 "Unigenes" transcripts with an average size of 1325 bp. Functional annotation using BLAST2GO resulted in the identification of 50,301 annotated transcripts (80.71%) and GO assigned to 18,814 transcripts. KEGG pathway annotation of the "Unigenes" revealed that 2465 transcripts could be assigned to 242 KEGG pathway maps wherein the number of transcripts involved in secondary metabolism was distinct in root and leaf transcriptome. Among the secondary metabolite biosynthesis pathways, the cluster of "Unigenes" encoding enzymes of 'Phenylpropanoid biosynthesis pathway' represents the largest group (84 transcripts) followed by 'Terpenoid Backbone biosynthesis' (48 transcripts). The transcript levels of the candidate unigenes encoding key enzymes of phenylpropanoid (PAL, TAL) and flavanoid biosynthesis (CHS, ANS, FLS) pathways were up-regulated in root, while the expression levels of candidate "Unigenes" transcript for monoterpenoid (DXS, ISPF), diterpenoid biosynthesis (SPS, SDS) and indole alkaloid pathways (STR) were significantly higher in leaf of P. zeylanica. Interestingly, validation of differential gene expression profile by qRT-PCR also confirmed that candidate "Unigenes" enzymes of phenylpropanoid and flavonoid biosynthesis were highly expressed in the root, while the key regulatory enzymes of terpenoid and indole alkaloid compounds were up-regulated in the leaf, suggesting that (differences in) the levels of these functional genes could be attributed to the (differential) pharmacological activity (between root and leaf) in tissues of P. zeylanica.

Asparagus racemosus bZIP transcription factor-regulated squalene epoxidase (ArSQE) promotes germination and abiotic stress tolerance in transgenic tobacco

A. racemosus is a rich source of pharmacologically active steroidal saponins. Most of the studies are related to its chemistry and pharmacology, but the pathway involved in the biosynthesis of steroidal saponin is not much emphasized. Squalene epoxidase acts as a rate-limiting enzyme in this biosynthesis. In this study, we have selected root specific squalene epoxidase ArSQE from A. racemosus for its characterization. ArSQE was able to complement ergosterol auxotrophy in erg1 yeast mutants. Mutants were sensitive to the antifungal drug terbinafine, whereas ArSQE complementation made them tolerant to the same drug. ArSQE plays a significant role in early germination in transgenic tobacco. The transgenic tobacco seedlings overexpressing ArSQE were tolerant to terbinafine and abiotic stress. Expression analysis of transcripts in ArSQE transgenic lines suggests that it mostly affects ABA, GA, stress, and sterol related functions in transgenic tobacco. Further, root specific MeJA responsive A. racemosus bZIP transcription factors (TFs), ArTGA1 and ArTGA2, were identified that bind to MeJA responsive cis-element present in the promoter region of ArSQE. Characterization of ArSQE of A. racemosus provides new information about its regulation through MeJA responsive bZIP TF along with its role in the development and abiotic stress response in transgenic tobacco.

Transcriptome analysis of Asparagus officinalis reveals genes involved in the biosynthesis of rutin and protodioscin

Garden asparagus (Asparagus officinalis L.) is a popular vegetable cultivated worldwide. The secondary metabolites in its shoot are helpful for human health. We analyzed A. officinalis transcriptomes and identified differentially expressed genes (DEGs) involved in the biosynthesis of rutin and protodioscin, which are health-promoting functional compounds, and determined their association with stem color. We sequenced the complete mRNA transcriptome using the Illumina high-throughput sequencing platform in one white, three green, and one purple asparagus cultivars. A gene set was generated by de novo assembly of the transcriptome sequences and annotated using a BLASTx search. To investigate the relationship between the contents of rutin and protodioscin and their gene expression levels, rutin and protodioscin were analyzed using high-performance liquid chromatography. A secondary metabolite analysis using high-performance liquid chromatography showed that the rutin content was higher in green asparagus, while the protodioscin content was higher in white asparagus. We studied the genes associated with the biosynthesis of the rutin and protodioscin. The transcriptomes of the five cultivars generated 336 599 498 high-quality clean reads, which were assembled into 239 873 contigs with an average length of 694 bp, using the Trinity v2.4.0 program. The green and white asparagus cultivars showed 58 932 DEGs. A comparison of rutin and protodioscin biosynthesis genes revealed that 12 of the 57 genes associated with rutin and two of the 50 genes associated with protodioscin showed more than four-fold differences in expression. These DEGs might have caused a variation in the contents of these two metabolites between green and white asparagus. The present study is possibly the first to report transcriptomic gene sets in asparagus. The DEGs putatively involved in rutin and protodioscin biosynthesis might be useful for molecular engineering in asparagus.

Identification of candidate genes involved in steroidal alkaloids biosynthesis in organ-specific transcriptomes of Veratrum nigrum L

Veratrum nigrum is protected plant of Melanthiaceae family, able to synthetize unique steroidal alkaloids important for pharmacy. Transcriptomes from leaves, stems and rhizomes of in vitro maintained V. nigrum plants were sequenced and annotated for genes and markers discovery. Sequencing of samples derived from the different organs resulted in a total of 108,511 contigs with a mean length of 596 bp. Transcripts derived from leaf and stalk were annotated at 28%, and 38% in Nr nucleotide database, respectively. The sequencing revealed 949 unigenes related with lipid metabolism, including 73 transcripts involved in steroids and genus-specific steroid alkaloids biosynthesis. Additionally, 3203 candidate SSRs markers we identified in unigenes with average density of one SSR locus every 6.2 kb sequence. Unraveling of biochemical machinery of the pathway responsible for steroidal alkaloids will open possibility to design and optimize biotechnological process. The transcriptomic data provide valuable resources for biochemical, molecular genetics, comparative transcriptomics, functional genomics, ecological and evolutionary studies of V. nigrum.

The transcriptome enables the identification of candidate genes behind medicinal value of Drumstick tree (Moringa oleifera)

Moringa oleifera is a plant well-known for its nutrition value, drought resistance and medicinal properties. cDNA libraries from five different tissues (leaf, root, stem, seed and flower) of M. oleifera cultivar Bhagya were generated and sequenced. We developed a bioinformatics pipeline to assemble transcriptome, along with the previously published M. oleifera genome, to predict 17,148 gene models. Few candidate genes related to biosynthesis of secondary metabolites, vitamins and ion transporters were identified. Expressions were further confirmed by real-time quantitative PCR experiments for few promising leads. Quantitative estimation of metabolites, as well as elemental analysis, was also carried out to support our observations. Enzymes in the biosynthesis of vitamins and metabolites like quercetin and kaempferol are highly expressed in leaves, flowers and seeds. The expression of iron transporters and calcium storage proteins were observed in root and leaves. In general, leaves retain the highest amount of small molecules of interest.

Transcriptome analyses of Paris polyphylla var. chinensis, Ypsilandra thibetica, and Polygonatum kingianum characterize their steroidal saponin biosynthesis pathway

Steroidal saponins, one of the most diverse groups of plant-derived natural products, elicit biological and pharmacological activities; however, the genes involved in their biosynthesis and the corresponding biosynthetic pathway in monocotyledon plants remain unclear. This study aimed to identify genes involved in the biosynthesis of steroidal saponins by performing a comparative analysis among transcriptomes of Paris polyphylla var. chinensis (PPC), Ypsilandra thibetica (YT), and Polygonatum kingianum (PK). De novo transcriptome assemblies generated 57,537, 140,420, and 151,773 unigenes from PPC, YT, and PK, respectively, of which 56.54, 47.81, and 44.30% were successfully annotated, respectively. Among the transcriptomes for PPC, YT, and PK, we identified 194, 169, and 131; 17, 14, and 26; and, 80, 122, and 113 unigenes corresponding to terpenoid backbone biosynthesis; sesquiterpenoid and triterpenoid biosynthesis; and, steroid biosynthesis pathways, respectively. These genes are putatively involved in the biosynthesis of cholesterol that is the primary precursor of steroidal saponins. Phylogenetic analyses indicated that lanosterol synthase may be exclusive to dicotyledon plant species, and the cytochrome P450 unigenes were closely related to clusters CYP90B1 and CYP734A1, which are UDP-glycosyltransferases unigenes homologous with the UGT73 family. Thus, unigenes of β-glucosidase may be candidate genes for catalysis of later period modifications of the steroidal saponin skeleton. Our data provide evidence to support the hypothesis that monocotyledons biosynthesize steroidal saponins from cholesterol via the cycloartenol pathway.

Comparative analysis of proteomic and metabolomic profiles of different species of Paris

An extract prepared from species of Paris is the most widely consumed herbal product in China. The genus Paris includes a variety of genotypes with different medicinal component contents but only two are defined as official sources. Closely related species have different medicinal properties because of differential expression of proteins and metabolites. To better understand the molecular basis of these differences, we examined proteomic and metabolomic changes in rhizomes of P. polyphylla var. chinensis, P. polyphylla var. yunnanensis, and P. fargesii var. fargesii using a technique known as sequential window acquisition of all theoretical mass spectra as well as gas chromatography-time-of-flight mass spectrometry. In total, 419 proteins showed significant abundance changes, and 33 metabolites could be used to discriminate Paris species. A complex analysis of proteomic and metabolomic data revealed a higher efficiency of sucrose utilization and an elevated protein abundance in the sugar metabolic pathway of P. polyphylla var. chinensis. The pyruvate content and efficiency of acetyl-CoA-utilization in saponin biosynthesis were also higher in P. polyphylla var. chinensis than in the other two species. The results expand our understanding of the proteome and metabolome of Paris and offer new insights into the species-specific traits of these herbaceous plants. SIGNIFICANCE: The traditional Chinese medicine Paris is the most widely consumed herbal product for the treatment of joint pain, rheumatoid arthritis and antineoplastic. All Paris species have roughly the same morphological characteristics; however, different members have different medicinal compound contents. Efficient exploitation of genetic diversity is a key factor in the development of rare medicinal plants with improved agronomic traits and malleability to challenging environmental conditions. Nevertheless, only a partial understanding of physiological and molecular mechanisms of different plants of Paris can be achieved without proteomics. To better understand the molecular basis of these differences and facilitate the use of other Paris species, we examine proteomic metabolomic changes in rhizomes of Paris using the technique known as SWATH-MS and GC/TOF-MS. Our research has provided information that can be used in other studies to compare metabolic traits in different Paris species. Our findings can also serve as a theoretical basis for the selection and cultivation of other Paris species with a higher medicinal value.

Tissue-specific transcriptional biomarkers in medicinal plants: Application of large-scale meta-analysis and computational systems biology

Biosynthesis of secondary metabolites in plant is a complex process, regulated by many genes and influenced by several factors. In recent years, the next-generation sequencing (NGS) technology and advanced statistical analysis such as meta-analysis and computational systems biology have provided novel opportunities to overcome biological complexity. Here, we performed a meta-analysis on publicly available transcriptome datasets of twelve economically significant medicinal plants to identify differentially expressed genes (DEGs) between shoot and root tissues and to find the key molecular features which may be effective in the biosynthesis of secondary metabolites. Meta-analysis identified a total of 880 genes with differential expression between two tissues. Functional enrichment and KEGG pathway analysis indicated that the functions of those DEGs are highly associated with the developmental process, starch metabolic process, response to stimulus, porphyrin and chlorophyll metabolism, biosynthesis of secondary metabolites and phenylalanine metabolism. In addition, systems biology analysis of the DEGs was applied to find protein-protein interaction network and discovery of significant modules. The detected modules were associated with hormone signal transduction, transcription repressor activity, response to light stimulus and epigenetic processes. Finally, analysis was extended to search for putative miRNAs that are associated with DEGs. A total of 31 miRNAs were detected which belonged to 16 conserved families. The present study provides a comprehensive view to better understand the tissue-specific expression of genes and mechanisms involved in secondary metabolites synthesis and may provide candidate genes for future researches to improve yield of secondary metabolites.

Recent advances in steroidal saponins biosynthesis and in vitro production

Steroidal saponins exhibited numerous pharmacological activities due to the modification of their backbone by different cytochrome P450s (P450) and UDP glycosyltransferases (UGTs). Plant-derived steroidal saponins are not sufficient for utilizing them for commercial purpose so in vitro production of saponin by tissue culture, root culture, embryo culture, etc, is necessary for its large-scale production. Saponin glycosides are the important class of plant secondary metabolites, which consists of either steroidal or terpenoidal backbone. Due to the existence of a wide range of medicinal properties, saponin glycosides are pharmacologically very important. This review is focused on important medicinal properties of steroidal saponin, its occurrence, and biosynthesis. In addition to this, some recently identified plants containing steroidal saponins in different parts were summarized. The high throughput transcriptome sequencing approach elaborates our understanding related to the secondary metabolic pathway and its regulation even in the absence of adequate genomic information of non-model plants. The aim of this review is to encapsulate the information related to applications of steroidal saponin and its biosynthetic enzymes specially P450s and UGTs that are involved at later stage modifications of saponin backbone. Lastly, we discussed the in vitro production of steroidal saponin as the plant-based production of saponin is time-consuming and yield a limited amount of saponins. A large amount of plant material has been used to increase the production of steroidal saponin by employing in vitro culture technique, which has received a lot of attention in past two decades and provides a way to conserve medicinal plants as well as to escape them for being endangered.

Comprehensive metabolic and transcriptomic profiling of various tissues provide insights for saponin biosynthesis in the medicinally important Asparagus racemosus

Asparagus racemosus (Shatavari), belongs to the family Asparagaceae and is known as a “curer of hundred diseases” since ancient time. This plant has been exploited as a food supplement to enhance immune system and regarded as a highly valued medicinal plant in Ayurvedic medicine system for the treatment of various ailments such as gastric ulcers, dyspepsia, cardiovascular diseases, neurodegenerative diseases, cancer, as a galactogogue and against several other diseases. In depth metabolic fingerprinting of various parts of the plant led to the identification of 13 monoterpenoids exclusively present in roots. LC-MS profiling led to the identification of a significant number of steroidal saponins (33). However, we have also identified 16 triterpene saponins for the first time in A. racemosus. In order to understand the molecular basis of biosynthesis of major components, transcriptome sequencing from three different tissues (root, leaf and fruit) was carried out. Functional annotation of A. racemosus transcriptome resulted in the identification of 153 transcripts involved in steroidal saponin biosynthesis, 45 transcripts in triterpene saponin biosynthesis, 44 transcripts in monoterpenoid biosynthesis and 79 transcripts in flavonoid biosynthesis. These findings will pave the way for better understanding of the molecular basis of steroidal saponin, triterpene saponin, monoterpenoids and flavonoid biosynthesis in A. racemosus.

De novo sequencing, assembly and characterisation of Aloe vera transcriptome and analysis of expression profiles of genes related to saponin and anthraquinone metabolism

BACKGROUND

Aloe vera is a perennial, succulent, drought-resistant plant that exhibits many pharmacological characteristics such as wound healing ability against skin burns, anti-ulcer, anti-inflammatory, anti-tumor, anti-viral, anti-hypercholesterolemic, anti-hyperglycemic, anti-asthmatic and much more. Despite great medicinal worth, little genomic information is available on Aloe vera. This study is an initiative to explore the full-scale functional genomics of Aloe vera by generating whole transcriptome sequence database, using Illumina HiSeq technology and its progressive annotation specifically with respect to the metabolic specificity of the plant.

RESULTS

Transcriptome sequencing of root and leaf tissue of Aloe vera was performed using Illumina paired-end sequencing technology. De novo assembly of high quality paired-end reads, resulted into 1,61,733 and 2,21,792 transcripts with mean length of 709 and 714 nucleotides for root and leaf respectively. The non-redundant transcripts were clustered using CD-HIT-EST, yielding a total of 1,13,063 and 1,41,310 unigenes for root and leaf respectively. A total of 6114 and 6527 CDS for root and leaf tissue were enriched into 24 different biological pathway categories using KEGG pathway database. DGE profile prepared by calculating FPKM values was analyzed for differential expression of specific gene encoding enzymes involved in secondary metabolite biosynthesis. Sixteen putative genes related to saponin, lignin, anthraquinone, and carotenoid biosynthesis were selected for quantitative expression by real-time PCR. DGE as well as qRT PCR expression analysis represented up-regulation of secondary metabolic genes in root as compared to leaf. Furthermore maximum number of genes was found to be up-regulated after the induction of methyl jasmonate, which stipulates the association of secondary metabolite synthesis with the plant's defense mechanism during stress. Various transcription factors including bHLH, NAC, MYB were identified by searching predicted CDS against PlantTFdb.

CONCLUSIONS

This is the first transcriptome database of Aloe vera and can be potentially utilized to characterize the genes involved in the biosynthesis of important secondary metabolites, metabolic regulation, signal transduction mechanism, understanding function of a particular gene in the biology and physiology of plant of this species as well as other species of Aloe genus.

Identification, characterization and expression analysis of genes involved in steroidal saponin biosynthesis in Dracaena cambodiana

Dracaena cambodiana is a traditional medicinal plant used for producing dragon's blood. The plants and dragon's blood of D. cambodiana contain a rich variety of steroidal saponins. However, little is known about steroidal saponin biosynthesis and its regulation in D. cambodiana. Here, 122 genes encoding enzymes involved in steroidal saponin biosynthesis were identified based on transcriptome data, with 29 of them containing complete open reading frames (ORF). Transcript expression analysis revealed that several genes related to steroidal saponin biosynthesis showed distinct tissue-specific expression patterns; the expression levels of genes encoding the key enzymes involved in the biosynthesis and early modification of steroidal saponins were significantly down-regulated in the stems in response to the inducer of dragon's blood, exhibiting positive correlations with the content of steroidal saponins. These results provide insights on the steroidal saponins biosynthetic pathway and mechanisms underlying induced formation of dragon's blood in D. cambodiana.

Changes in fenugreek transcriptome induced by methyl jasmonate and steroid precursors revealed by RNA-Seq

Changes in fenugreek transcriptome related to enhanced production of steroids were induced by methyl jasmonate, cholesterol and squalene, and recorded using RNA-seq. A total of 112,850 unigenes were obtained after de novo assembling of next generation sequencing data, and used for functional annotations. In steroidal saponins pathway, transcripts involved in mevalonate, terpenoid backbone and plant sterol synthesis were annotated. Overexpression of several transcripts from phytosterol biosynthesis pathway was confirmed by quantitative RT-PCR. In diosgenin biosynthesis pathway, fatty acid ω-hydroxylase (CYP86A2) and steroid 22-alpha-hydroxylase (CYP90B1) genes were annotated in all induced transcriptomes. Moreover, direct sequencing confirmed increased levels of CYP90B1, unspecific monooxygenase and 26-hydroxylase genes in plants with elevated level of diosgenin. New unigenes corresponding to enzymes involved in biosynthesis of diosgenin from cycloartenol via cholesterol were obtained and the role of CYP72A family in steroidal saponin biosynthesis was proposed. Additional support for biosynthetic pathway from cycloartenol to diosgenin was provided.

RNA-sequencing-based transcriptome and biochemical analyses of steroidal saponin pathway in a complete set of Allium fistulosum-A. cepa monosomic addition lines

The genus Allium is a rich source of steroidal saponins, and its medicinal properties have been attributed to these bioactive compounds. The saponin compounds with diverse structures play a pivotal role in Allium’s defense mechanism. Despite numerous studies on the occurrence and chemical structure of steroidal saponins, their biosynthetic pathway in Allium species is poorly understood. The monosomic addition lines (MALs) of the Japanese bunching onion (A. fistulosum, FF) with an extra chromosome from the shallot (A. cepa Aggregatum group, AA) are powerful genetic resources that enable us to understand many physiological traits of Allium. In the present study, we were able to isolate and identify Alliospiroside A saponin compound in A. fistulosum with extra chromosome 2A from shallot (FF2A) and its role in the defense mechanism against Fusarium pathogens. Furthermore, to gain molecular insight into the Allium saponin biosynthesis pathway, high-throughput RNA-Seq of the root, bulb, and leaf of AA, MALs, and FF was carried out using Illumina's HiSeq 2500 platform. An open access Allium Transcript Database (Allium TDB, http://alliumtdb.kazusa.or.jp) was generated based on RNA-Seq data. The resulting assembled transcripts were functionally annotated, revealing 50 unigenes involved in saponin biosynthesis. Differential gene expression (DGE) analyses of AA and MALs as compared with FF (as a control) revealed a strong up-regulation of the saponin downstream pathway, including cytochrome P450, glycosyltransferase, and beta-glucosidase in chromosome 2A. An understanding of the saponin compounds and biosynthesis-related genes would facilitate the development of plants with unique saponin content and, subsequently, improved disease resistance.

De novo leaf and root transcriptome analysis to identify putative genes involved in triterpenoid saponins biosynthesis in Hedera helix L

Hedera helix L. is an important traditional medicinal plant in Europe. The main active components are triterpenoid saponins, but none of the potential enzymes involved in triterpenoid saponins biosynthesis have been discovered and annotated. Here is reported the first study of global transcriptome analyses using the Illumina HiSeq™ 2500 platform for H. helix. In total, over 24 million clean reads were produced and 96,333 unigenes were assembled, with an average length of 1385 nt; more than 79,085 unigenes had at least one significant match to an existing gene model. Differentially Expressed Gene analysis identified 6,222 and 7,012 unigenes which were expressed either higher or lower in leaf samples when compared with roots. After functional annotation and classification, two pathways and 410 unigenes related to triterpenoid saponins biosynthesis were discovered. The accuracy of these de novo sequences was validated by RT-qPCR analysis and a RACE clone. These data will enrich our knowledge of triterpenoid saponin biosynthesis and provide a theoretical foundation for molecular research on H. helix.

Identification of critical genes associated with lignin biosynthesis in radish (Raphanus sativus L.) by de novo transcriptome sequencing

Radish is an important root vegetable crop with high nutritional, economic, and medicinal value. Lignin is an important secondary metabolite possessing a great effect on plant growth and product quality. To date, lignin biosynthesis-related genes have been identified in some important plant species. However, little information on characterization of critical genes involved in plant lignin biosynthesis is available in radish. In this study, a total of 71,148 transcripts sequences were obtained from radish root, of which 66 assembled unigenes and ten candidate genes were identified to be involved in lignin monolignol biosynthesis. Full-length cDNA sequences of seven randomly selected genes were isolated and sequenced from radish root, and the assembled unigenes covered more than 80% of their corresponding cDNA sequences. Moreover, the lignin content gradually accumulated in leaf during the developmental stages, and it increased from pre-cortex to cortex splitting stage, followed by a decrease at thickening stage and then increased at mature stage in root. RT-qPCR analysis revealed that all these genes except RsF5H exhibited relatively low expression level in root at thickening stage. The expression profiles of Rs4CL5, RsCCoAOMT1, and RsCOMT genes were consistent with the changes of root lignin content, implying that these candidate genes may play important roles in lignin formation in radish root. These findings would provide valuable information for identification of lignin biosynthesis-related genes and facilitate dissection of molecular mechanism underlying lignin biosynthesis in radish and other root vegetable crops.

Comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri identifies potential genes related to triterpenoid saponin biosynthesis

Background

Bacopa monnieri commonly known as Brahmi is utilized in Ayurveda to improve memory and many other human health benefits. Bacosides enriched standardized extract of Bacopa monnieri is being marketed as a memory enhancing agent. In spite of its well known pharmacological properties it is not much studied in terms of transcripts involved in biosynthetic pathway and its regulation that controls the secondary metabolic pathway in this plant. The aim of this study was to identify the potential transcripts and provide a framework of identified transcripts involved in bacosides production through transcriptome assembly.

Results

We performed comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri in two independent biological replicate and obtained 22.48 million and 22.0 million high quality processed reads in shoot and root respectively. After de novo assembly and quantitative assessment total 26,412 genes got annotated in root and 18,500 genes annotated in shoot sample. Quality of raw reads was determined by using SeqQC-V2.2. Assembled sequences were annotated using BLASTX against public database such as NR or UniProt. Searching against the KEGG pathway database indicated that 37,918 unigenes from root and 35,130 unigenes from shoot were mapped to 133 KEGG pathways. Based on the DGE data we found that most of the transcript related to CYP450s and UDP-glucosyltransferases were specifically upregulated in shoot tissue as compared to root tissue. Finally, we have selected 43 transcripts related to secondary metabolism including transcription factor families which are differentially expressed in shoot and root tissues were validated by qRT-PCR and their expression level were monitored after MeJA treatment and wounding for 1, 3 and 5 h.

Conclusions

This study not only represents the first de novo transcriptome analysis of Bacopa monnieri but also provides information about the identification, expression and differential tissues specific distribution of transcripts related to triterpenoid sapogenin which is one of the most important pharmacologically active secondary metabolite present in Bacopa monnieri. The identified transcripts in this study will establish a foundation for future studies related to carrying out the metabolic engineering for increasing the bacosides biosynthesis and its regulation for human health benefits.

Electronic supplementary material

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

De novo assembly and annotation of the Zhe-Maidong (Ophiopogon japonicus (L.f.) Ker-Gawl) transcriptome in different growth stages

Zhe-Maidong (Ophiopogon japonicus (L.f.) Ker-Gawl) is a traditional medicinal herb in the family Liliaceae that has significant pharmacological effects on immunity and cardiovascular disease. In this study, three different growth stages of Zhe-Maidong were investigated using RNA-seq, and a total of 16.4 Gb of raw data was obtained. After filtering and assembling, 96,738 unigenes with an average length of 605.3 bp were ultimately generated. A total of 77,300 unigenes were annotated using information from five databases, including the NT, NR, SwissProt, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases. Additionally, the mechanisms of flavonoid, saponin and polysaccharide biosynthesis and of accumulation at different stages of tuber development were also characterized. From the first to third years, the contents of flavonoids, saponins and polysaccharides all increased, whereas the expression levels of related genes decreased in the flavonoid and saponin pathways and first increased and then decreased in the polysaccharide pathway. The results of this study provide the most comprehensive expressed sequence resource for Zhe-Maidong and will expand the available O. japonicus gene library and facilitate further genome-wide research and analyses of this species.

Effect of in vitro drought stress on phenolic acids, flavonols, saponins, and antioxidant activity in Agave salmiana

This work evaluated the effect of in vitro drought stress on morphological characteristics, phenolic compounds, flavonols, saponin content, and antioxidant activity in plantlets of Agave salmiana. Drought stress was induced with polyethylene glycol (PEG) at 0, 10, 20 and 30% w/v in Murashige and Skoog solid medium. The determination of specific flavonols and saponins was achieved via HPLC-DAD and HPLC-ELSD, respectively. Compared with the control, plants grown in 30% PEG showed a change in the width of the leaves and a different color, showing less clarity and more darkening (L = 21.18, b = 14.27) and also had the lowest flavonol content, but the highest total saponin content (tigogenin glycoside, 163 mg of protodioscin equivalents/g dw) and the highest antioxidant activity. Total phenolic compounds did not significantly differ between treatments. Agave salmiana plants cultured in vitro increased their saponin content and antioxidant activity in response to drought stress induced via PEG.

Next-generation sequencing of representational difference analysis products for identification of genes involved in diosgenin biosynthesis in fenugreek (Trigonella foenum-graecum)

Representational difference analysis of cDNA was performed and differential products were sequenced and annotated. Candidate genes involved in biosynthesis of diosgenin in fenugreek were identified. Detailed mechanism of diosgenin synthesis was proposed. Fenugreek (Trigonella foenum-graecum L.) is a valuable medicinal and crop plant. It belongs to Fabaceae family and has a unique potential to synthesize valuable steroidal saponins, e.g., diosgenin. Elicitation (methyl jasmonate) and precursor feeding (cholesterol and squalene) were used to enhance the content of sterols and steroidal sapogenins in in vitro grown plants for representational difference analysis of cDNA (cDNA-RDA). To identify candidate genes involved in diosgenin biosynthesis, differential, factor-specific libraries were subject to the next-generation sequencing. Approximately 9.9 million reads were obtained, trimmed, and assembled into 31,491 unigenes with an average length of 291 bp. Then, functional annotation and gene ontogeny enrichment analysis was performed by aligning all-unigenes with public databases. Within the transcripts related to sterol and steroidal saponin biosynthesis, we discovered novel candidate genes of diosgenin biosynthesis and validated their expression using quantitative RT-PCR analysis. Based on these findings, we supported the idea that diosgenin is biosynthesized from cycloartenol via cholesterol. This is the first report on the next-generation sequencing of cDNA-RDA products. Analysis of the transcriptomes enriched in low copy sequences contributed substantially to our understanding of the biochemical pathways of steroid synthesis in fenugreek.

Spatial transcriptome analysis provides insights of key gene(s) involved in steroidal saponin biosynthesis in medicinally important herb Trillium govanianum

Trillium govanianum, an endangered medicinal herb native to the Himalaya, is less studied at the molecular level due to the non-availability of genomic resources. To facilitate the basic understanding of the key genes and regulatory mechanism of pharmaceutically important biosynthesis pathways, first spatial transcriptome sequencing of T. govanianum was performed. 151,622,376 (~11.5 Gb) high quality reads obtained using paired-end Illumina sequencing were de novo assembled into 69,174 transcripts. Functional annotation with multiple public databases identified array of genes involved in steroidal saponin biosynthesis and other secondary metabolite pathways including brassinosteroid, carotenoid, diterpenoid, flavonoid, phenylpropanoid, steroid and terpenoid backbone biosynthesis, and important TF families (bHLH, MYB related, NAC, FAR1, bZIP, B3 and WRKY). Differentially expressed large number of transcripts, together with CYPs and UGTs suggests involvement of these candidates in tissue specific expression. Combined transcriptome and expression analysis revealed that leaf and fruit tissues are the main site of steroidal saponin biosynthesis. In conclusion, comprehensive genomic dataset created in the current study will serve as a resource for identification of potential candidates for genetic manipulation of targeted bioactive metabolites and also contribute for development of functionally relevant molecular marker resource to expedite molecular breeding and conservation efforts in T. govanianum.

Transcriptome Analysis of Salicylic Acid Treatment in Rehmannia glutinosa Hairy Roots Using RNA-seq Technique for Identification of Genes Involved in Acteoside Biosynthesis

Rehmannia glutinosa is a common bulk medicinal material that has been widely used in China due to its active ingredients. Acteoside, one of the ingredients, has antioxidant, antinephritic, anti-inflammatory, hepatoprotective, immunomodulatory, and neuroprotective effects, is usually selected as a quality-control component for R. glutinosa herb in the Chinese Pharmacopeia. The acteoside biosynthesis pathway in R. glutinosa has not yet been clearly established. Herein, we describe the establishment of a genetic transformation system for R. glutinosa mediated by Agrobacterium rhizogenes. We screened the optimal elicitors that markedly increased acteoside accumulation in R. glutinosa hairy roots. We found that acteoside accumulation dramatically increased with the addition of salicylic acid (SA); the optimal SA dose was 25 μmol/L for hairy roots. RNA-seq was applied to analyze the transcriptomic changes in hairy roots treated with SA for 24 h in comparison with an untreated control. A total of 3,716, 4,018, and 2,715 differentially expressed transcripts (DETs) were identified in 0 h-vs.-12 h, 0 h-vs.-24 h, and 12 h-vs.-24 h libraries, respectively. KEGG pathway-based analysis revealed that 127 DETs were enriched in "phenylpropanoid biosynthesis." Of 219 putative unigenes involved in acteoside biosynthesis, 54 were found to be up-regulated at at least one of the time points after SA treatment. Selected candidate genes were analyzed by quantitative real-time PCR (qRT-PCR) in hairy roots with SA, methyl jasmonate (MeJA), AgNO₃ (Ag⁺), and putrescine (Put) treatment. All genes investigated were up-regulated by SA treatment, and most candidate genes were weakly increased by MeJA to some degree. Furthermore, transcription abundance of eight candidate genes in tuberous roots of the high-acteoside-content (HA) cultivar QH were higher than those of the low-acteoside-content (LA) cultivar Wen 85-5. These results will pave the way for understanding the molecular basis of acteoside biosynthesis in R. glutinosa, and can serve as a basis for future validation studies.

Integrated analysis of transcriptome and metabolites reveals an essential role of metabolic flux in starch accumulation under nitrogen starvation in duckweed

BACKGROUND

Duckweed is considered a promising source of energy due to its high starch content and rapid growth rate. Starch accumulation in duckweed involves complex processes that depend on the balanced expression of genes controlled by various environmental and endogenous factors. Previous studies showed that nitrogen starvation induces a global stress response and results in the accumulation of starch in duckweed. However, relatively little is known about the mechanisms underlying the regulation of starch accumulation under conditions of nitrogen starvation.

RESULTS

In this study, we used next-generation sequencing technology to examine the transcriptome responses of Lemna aequinoctialis 6000 at three stages (0, 3, and 7 days) during nitrogen starvation in the presence of exogenously applied sucrose. Overall, 2522, 628, and 1832 differentially expressed unigenes (DEGs) were discovered for the treated and control samples. Clustering and enrichment analysis of DEGs revealed several biological processes occurring under nitrogen starvation. Genes involved in nitrogen metabolism showed the earliest responses to nitrogen starvation, whereas genes involved in carbohydrate biosynthesis were responded subsequently. The expression of genes encoding nitrate reductase, glutamine synthetase, and glutamate synthase was down-regulated under nitrogen starvation. The expression of unigenes encoding enzymes involved in gluconeogenesis was up-regulated, while the majority of unigenes involved in glycolysis were down-regulated. The metabolite results showed that more ADP-Glc was accumulated and lower levels of UDP-Glc were accumulated under nitrogen starvation, the activity of AGPase was significantly increased while the activity of UGPase was dramatically decreased. These changes in metabolite levels under nitrogen starvation are roughly consistent with the gene expression changes in the transcriptome.

CONCLUSIONS

Based on these results, it can be concluded that the increase of ADP-glucose and starch contents under nitrogen starvation is a consequence of increased output from the gluconeogenesis and TCA pathways, accompanied with the reduction of lipids and pectin biosynthesis. The results provide novel insights into the underlying mechanisms of starch accumulation during nitrogen starvation, which provide a foundation for the improvement of advanced bioethanol production in duckweed.

Transcriptome and Metabolite analysis reveal candidate genes of the cardiac glycoside biosynthetic pathway from Calotropis procera

Calotropis procera is a medicinal plant of immense importance due to its pharmaceutical active components, especially cardiac glycosides (CG). As genomic resources for this plant are limited, the genes involved in CG biosynthetic pathway remain largely unknown till date. Our study on stage and tissue specific metabolite accumulation showed that CG's were maximally accumulated in stems of 3 month old seedlings. De novo transcriptome sequencing of same was done using high throughput Illumina HiSeq platform generating 44074 unigenes with average mean length of 1785 base pair. Around 66.6% of unigenes were annotated by using various public databases and 5324 unigenes showed significant match in the KEGG database involved in 133 different pathways of plant metabolism. Further KEGG analysis resulted in identification of 336 unigenes involved in cardenolide biosynthesis. Tissue specific expression analysis of 30 putative transcripts involved in terpenoid, steroid and cardenolide pathways showed a positive correlation between metabolite and transcript accumulation. Wound stress elevated CG levels as well the levels of the putative transcripts involved in its biosynthetic pathways. This result further validated the involvement of identified transcripts in CGs biosynthesis. The identified transcripts will lay a substantial foundation for further research on metabolic engineering and regulation of cardiac glycosides biosynthesis pathway genes.

De Novo Transcriptome Assembly and Annotation of the Leaves and Callus of Cyclocarya Paliurus (Bata1) Iljinskaja

Cyclocarya Paliurus (Bata1) Iljinskaja contains various bioactive secondary metabolites especially in leaves, such as triterpenes, flavonoids, polysaccharides and alkaloids, and its leaves are widely used as an hyperglycemic tea in China. In the present paper, we sequenced the transcriptome of the leaves and callus of Cyclocarya Paliurus using Illumina Hiseq 4000 platform. After sequencing and de novo assembly, a total of 65,654 unigenes were generated with an N50 length of 1,244bp. Among them, 35,041 (53.37%) unigenes were annotated in NCBI Non-Redundant database, 19,453 (29.63%) unigenes were classified into Gene Ontology (GO) database, and 7,259 (11.06%) unigenes were assigned to Clusters of Orthologous Group (COG) categories. Furthermore, 11,697 (17.81%) unigenes were mapped onto 335 pathways in Kyoto Encyclopedia of Genes and Genomes (KEGG), among which 1,312 unigenes were identified to be involved in biosynthesis of secondary metabolites. In addition, a total of 11,247 putative simple sequence repeats (SSRs) were detected. This transcriptome dataset provides a comprehensive sequence resource for gene expression profiling, genetic diversity, evolution and further molecular genetics research on Cyclocarya Paliurus.

Genome wide transcriptome profiling reveals differential gene expression in secondary metabolite pathway of Cymbopogon winterianus

Advances in transcriptome sequencing provide fast, cost-effective and reliable approach to generate large expression datasets especially suitable for non-model species to identify putative genes, key pathway and regulatory mechanism. Citronella (Cymbopogon winterianus) is an aromatic medicinal grass used for anti-tumoral, antibacterial, anti-fungal, antiviral, detoxifying and natural insect repellent properties. Despite of having number of utilities, the genes involved in terpenes biosynthetic pathway is not yet clearly elucidated. The present study is a pioneering attempt to generate an exhaustive molecular information of secondary metabolite pathway and to increase genomic resources in Citronella. Using high-throughput RNA-Seq technology, root and leaf transcriptome was analysed at an unprecedented depth (11.7 Gb). Targeted searches identified majority of the genes associated with metabolic pathway and other natural product pathway viz. antibiotics synthesis along with many novel genes. Terpenoid biosynthesis genes comparative expression results were validated for 15 unigenes by RT-PCR and qRT-PCR. Thus the coverage of these transcriptome is comprehensive enough to discover all known genes of major metabolic pathways. This transcriptome dataset can serve as important public information for gene expression, genomics and function genomics studies in Citronella and shall act as a benchmark for future improvement of the crop.

Transcriptomic Analysis of Grapevine (cv. Summer Black) Leaf, Using the Illumina Platform

Proceeding to illumina sequencing, determining RNA integrity numbers for poly RNA were separated from each of the four developmental stages of cv. Summer Black leaves by using Illumina HiSeq^™ 2000. The sums of 272,941,656 reads were generated from vitis vinifera leaf at four different developmental stages, with more than 27 billion nucleotides of the sequence data. At each growth stage, RNA samples were indexed through unique nucleic acid identifiers and sequenced. KEGG annotation results depicted that the highest number of transcripts in 2,963 (2Avs4A) followed by 1Avs4A (2,920), and 3Avs4A (2,294) out of 15,614 (71%) transcripts were recorded. In comparison, a total of 1,532 transcripts were annotated in GOs, including Cellular component, with the highest number in “Cell part” 251 out of 353 transcripts (71.1%), followed by intracellular organelle 163 out of 353 transcripts (46.2%), while in molecular function and metabolic process 375 out of 525 (71.4%) transcripts, multicellular organism process 40 out of 525 (7.6%) transcripts in biological process were most common in 1Avs2A. While in case of 1Avs3A, cell part 476 out of 662 transcripts (71.9%), and membrane-bounded organelle 263 out of 662 transcripts (39.7%) were recorded in Cellular component. In the grapevine transcriptome, during the initial stages of leaf development 1Avs2A showed single transcript was down-regulated and none of them were up-regulated. While in comparison of 1A to 3A showed one up-regulated (photosystem II reaction center protein C) and one down regulated (conserved gene of unknown function) transcripts, during the hormone regulating pathway namely SAUR-like auxin-responsive protein family having 2 up-regulated and 7 down-regulated transcripts, phytochrome-associated protein showed 1 up-regulated and 9 down-regulated transcripts, whereas genes associated with the Leucine-rich repeat protein kinase family protein showed 7 up-regulated and 1 down-regulated transcript, meanwhile Auxin Resistant 2 has single up-regulated transcript in second developmental stage, although 3 were down-regulated at lateral growth stages (3A and 4A). In the present study, 489 secondary metabolic pathways related genes were identified during leaf growth, which mainly includes alkaloid (40), anthocyanins (21), Diterpenoid (144), Monoterpenoid (90) and Flavonoids (93). Quantitative real-time PCR was applied to validate 10 differentially expressed transcripts patterns from flower, leaf and fruit metabolic pathways at different growth stages.

De novo sequencing and assembly analysis of transcriptome in the Sodom apple (Calotropis gigantea)

Background

The Sodom apple (Calotropis gigantea), a member of the Asclepiadaceae family, is a large evergreen shrub native to continental Asia and northern Africa. As an important medicinal shrub and a fiber resource plant, there is an urgent need for developing molecular markers to facilitate breeding and genetic improvement of varieties.

Results

In this study, using the Illumina high throughput sequencing technique we obtained about 45 million paired end sequencing reads, De novo assembled and generated a total of 133,634 transcripts with a mean of 1837.47 bp in length. Based on protein homology searches against available databases, a total of 21,851 unigenes were functionally annotated. In particular, many transcripts that encode for putative proteins involved in fiber and secondary metabolite biosynthesis were identified and analyzed. Key fiber genes identified were validated experimentally through Real-Time PCR technique. Various transcription factors involved in regulating plant response to abiotic stress were also identified. In addition, based on the unigene sequences assembled, 11,623 microsatellites loci were detected, which provide very useful resources for developing microsatellite molecular markers.

Conclusion

This study is the first report on transcriptome information in the Calotropis species and provides rich gene transcript resources for conducting further studies on understanding the molecular basis of fiber and secondary metabolite biosynthesis, serving the genetic improvement and resource utilization in Calotropis plants.

Electronic supplementary material

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

Andrographis paniculata transcriptome provides molecular insights into tissue-specific accumulation of medicinal diterpenes

Background

Kalmegh (Andrographis paniculata) has been widely exploited in traditional medicine for the treatment of infectious diseases and health disorders. Ent-labdane-related diterpene (ent-LRD) specialized (i.e., secondary) metabolites of kalmegh such as andrographolide, neoandrographolide and 14-deoxy-11,12-didehydroandrographolide, are known for variety of pharmacological activities. However, due to the lack of genomic and transcriptomic information, underlying molecular basis of ent-LRDs biosynthesis has remained largely unknown. To identify candidate genes of the ent-LRD biosynthetic pathway, we performed comparative transcriptome analysis using leaf and root tissues that differentially accumulate ent-LRDs.

Results

De novo assembly of Illumina HiSeq2000 platform-generated paired-end sequencing reads resulted into 69,011 leaf and 64,244 root transcripts which were assembled into a total of 84,628 unique transcripts. Annotation of these transcripts to the Uniprot, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Carbohydrate-Active Enzymes (CAZy) databases identified candidate transcripts of the ent-LRD biosynthetic pathway. These included transcripts that encode enzymes of the plastidial 2C-methyl-D-erythritol-4-phosphate pathway which provides C5 isoprenoid precursors for the ent-LRDs biosynthesis, geranylgeranyl diphosphate synthase, class II diterpene synthase (diTPS), cytochrome P450 monooxygenase and glycosyltransferase. Three class II diTPSs (ApCPS1, ApCPS2 and ApCPS3) that showed distinct tissue-specific expression profiles and are phylogenetically related to the dicotyledon ent-copalyl diphosphate synthases, are identified. ApCPS1, ApCPS2 and ApCPS3 encode for 832-, 817- and 797- amino acids proteins of 55–63 % identity, respectively. Spatio-temporal patterns of transcripts and ent-LRDs accumulation are consistent with the involvement of ApCPS1 in general (i.e., primary) metabolism for the biosynthesis of phytohormone gibberellin, ApCPS2 in leaf specialized ent-LRDs biosynthesis and ApCPS3 in root diterpene biosynthesis. Moreover, simple sequence repeats (SSRs) that might assist in genotyping and developing specific chemotypes were identified in transcripts of the specialized metabolic pathways, including ent-LRDs.

Conclusions

Comparative analysis of root and leaf transcriptomes disclosed novel genes of the ent-LRD biosynthetic pathway, including three class II diTPSs that showed discrete spatio-temporal expression patterns; thus, suggesting their participation into distinct diterpene metabolic pathways of kalmegh. Overall, these results will be useful in understanding molecular basis of the medicinal ent-LRDs biosynthesis and developing breeding strategies for improving their yields.

Electronic supplementary material

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

Next Generation Sequencing and Transcriptome Analysis Predicts Biosynthetic Pathway of Sennosides from Senna (Cassia angustifolia Vahl.), a Non-Model Plant with Potent Laxative Properties

Senna (Cassia angustifolia Vahl.) is a world's natural laxative medicinal plant. Laxative properties are due to sennosides (anthraquinone glycosides) natural products. However, little genetic information is available for this species, especially concerning the biosynthetic pathways of sennosides. We present here the transcriptome sequencing of young and mature leaf tissue of Cassia angustifolia using Illumina MiSeq platform that resulted in a total of 6.34 Gb of raw nucleotide sequence. The sequence assembly resulted in 42230 and 37174 transcripts with an average length of 1119 bp and 1467 bp for young and mature leaf, respectively. The transcripts were annotated using NCBI BLAST with 'green plant database (txid 33090)', Swiss Prot, Kyoto Encylcopedia of Genes & Genomes (KEGG), Cluster of Orthologous Gene (COG) and Gene Ontology (GO). Out of the total transcripts, 40138 (95.0%) and 36349 (97.7%) from young and mature leaf, respectively, were annotated by BLASTX against green plant database of NCBI. We used InterProscan to see protein similarity at domain level, a total of 34031 (young leaf) and 32077 (mature leaf) transcripts were annotated against the Pfam domains. All transcripts from young and mature leaf were assigned to 191 KEGG pathways. There were 166 and 159 CDS, respectively, from young and mature leaf involved in metabolism of terpenoids and polyketides. Many CDS encoding enzymes leading to biosynthesis of sennosides were identified. A total of 10,763 CDS differentially expressing in both young and mature leaf libraries of which 2,343 (21.7%) CDS were up-regulated in young compared to mature leaf. Several differentially expressed genes found functionally associated with sennoside biosynthesis. CDS encoding for many CYPs and TF families were identified having probable roles in metabolism of primary as well as secondary metabolites. We developed SSR markers for molecular breeding of senna. We have identified a set of putative genes involved in various secondary metabolite pathways, especially those related to the synthesis of sennosides which will serve as an important platform for public information about gene expression, genomics, and functional genomics in senna.