Comparative glandular trichome transcriptome-based gene characterization reveals reasons for differential (-)-menthol biosynthesis in Mentha species

Isoprenyl diphosphate synthases of terpenoid biosynthesis in rose-scented geranium (Pelargonium graveolens)

The essential oil of Pelargonium graveolens (rose-scented geranium), an important aromatic plant, comprising mainly mono- and sesqui-terpenes, has applications in food and cosmetic industries. This study reports the characterization of isoprenyl disphosphate synthases (IDSs) involved in P. graveolens terpene biosynthesis. The six identified PgIDSs belonged to different classes of IDSs, comprising homomeric geranyl diphosphate synthases (GPPSs; PgGPPS1 and PgGPPS2), the large subunit of heteromeric GPPS or geranylgeranyl diphosphate synthases (GGPPSs; PgGGPPS), the small subunit of heteromeric GPPS (PgGPPS.SSUI and PgGPPS.SSUII), and farnesyl diphosphate synthases (FPPS; PgFPPS).All IDSs exhibited maximal expression in glandular trichomes (GTs), the site of aroma formation, and their expression except PgGPPS.SSUII was induced upon treatment with MeJA. Functional characterization of recombinant proteins revealed that PgGPPS1, PgGGPPS and PgFPPS were active enzymes producing GPP, GGPP/GPP, and FPP respectively, whereas both PgGPPS.SSUs and PgGPPS2 were inactive. Co-expression of PgGGPPS (that exhibited bifunctional G(G)PPS activity) with PgGPPS.SSUs in bacterial expression system showed lack of interaction between the two proteins, however, PgGGPPS interacted with a phylogenetically distant Antirrhinum majus GPPS.SSU. Further, transient expression of AmGPPS.SSU in P. graveolens leaf led to a significant increase in monoterpene levels. These findings provide insight into the types of IDSs and their role in providing precursors for different terpenoid components of P. graveolens essential oil.

Top and Side Lighting Induce Morphophysiological Improvements in Korean Ginseng Sprouts (Panax ginseng C.A. Meyer) Grown from One-Year-Old Roots

Nowadays, not only the roots, but also leaves and flowers of ginseng are increasingly popular ingredients in supplements for healthcare products and traditional medicine. The cultivation of the shade-loving crop, ginseng, is very demanding in terms of the light environment. Along with the intensity and duration, light direction is another important factor in regulating plant morphophysiology. In the current study, three lighting directions-top (T), side (S), or top + side (TS)-with an intensity of 30 ± 5 μmol·m-2·s-1 photosynthetic photon flux density (PPFD) were employed. Generally, compared with the single T lighting, the composite lighting direction, TS, was more effective in shaping the ginseng with improved characteristics, including shortened, thick shoots; enlarged, thick leaves; more leaf trichomes; earlier flower bud formation; and enhanced photosynthesis. The single S light resulted in the worst growth parameters and strongly inhibited the flower bud formation, leading to the latest flower bud observation. Additionally, the S lighting acted as a positive factor in increasing the leaf thickness and number of trichomes on the leaf adaxial surface. However, the participation of the T lighting weakened these traits. Overall, the TS lighting was the optimal direction for improving the growth and development traits in ginseng. This preliminary research may provide new ideas and orientations in ginseng cultivation lodging resistance and improving the supply of ginseng roots, leaves, and flowers to the market.

SDG26 Is Involved in Trichome Control in Arabidopsis thaliana: Affecting Phytohormones and Adjusting Accumulation of H3K27me3 on Genes Related to Trichome Growth and Development

Plant trichomes formed by specialized epidermal cells play a role in protecting plants from biotic and abiotic stresses and can also influence the economic and ornamental value of plant products. Therefore, further studies on the molecular mechanisms of plant trichome growth and development are important for understanding trichome formation and agricultural production. SET Domain Group 26 (SDG26) is a histone lysine methyltransferase. Currently, the molecular mechanism by which SDG26 regulates the growth and development of Arabidopsis leaf trichomes is still unclear. We found that the mutant of Arabidopsis (sdg26) possessed more trichomes on its rosette leaves compared to the wild type (Col-0), and the trichome density per unit area of sdg26 is significantly higher than that of Col-0. The content of cytokinins and jasmonic acid was higher in sdg26 than in Col-0, while the content of salicylic acid was lower in sdg26 than in Col-0, which is conducive to trichome growth. By measuring the expression levels of trichome-related genes, we found that the expression of genes that positively regulate trichome growth and development were up-regulated, while the negatively regulated genes were down-regulated in sdg26. Through chromatin immunoprecipitation sequencing (ChIP-seq) analysis, we found that SDG26 can directly regulate the expression of genes related to trichome growth and development such as ZFP1, ZFP5, ZFP6, GL3, MYB23, MYC1, TT8, GL1, GIS2, IPT1, IPT3, and IPT5 by increasing the accumulation of H3K27me3 on these genes, which further affects the growth and development of trichomes. This study reveals the mechanism by which SDG26 affects the growth and development of trichomes through histone methylation. The current study provides a theoretical basis for studying the molecular mechanism of histone methylation in regulating leaf trichome growth and development and perhaps guiding the development of new crop varieties.

Transcriptome analysis of transcription factors and enzymes involved in monoterpenoid biosynthesis in different chemotypes of Mentha haplocalyx Briq

Background

The main active ingredients of Mentha haplocalyx Briq. essential oils are monoterpenes. According to the component of essential oils, M. haplocalyx can be divided into different chemotypes. Chemotype variation is widespread in Mentha plants but its formation mechanism is unclear.

Methods

We selected the stable chemotype l-menthol, pulegone, and carvone of M. haplocalyx for transcriptome sequencing. To further investigate the variation of chemotypes, we analyzed the correlation between differential transcription factors (TFs) and key enzymes.

Results

Fourteen unigenes related to monoterpenoid biosynthesis were identified, among which (+)-pulegone reductase (PR) and (-)-menthol dehydrogenase (MD) were significantly upregulated in l-menthol chemotype and (-)-limonene 6-hydroxylase was significantly upregulated in carvone chemotype. In addition, 2,599 TFs from 66 families were identified from transcriptome data and the differential TFs included 113 TFs from 34 families. The families of bHLH, bZIP, AP2/ERF, MYB, and WRKY were highly correlated with the key enzymes PR, MD, and (-)-limonene 3-hydroxylase (L3OH) in different M. haplocalyx chemotypes (r > 0.85). The results indicate that these TFs regulate the variation of different chemotypes by regulating the expression patterns of PR, MD, and L3OH. The results of this study provide a basis for revealing the molecular mechanism of the formation of different chemotypes and offer strategies for effective breeding and metabolic engineering of different chemotypes in M. haplocalyx.

Glandular trichome specificity of menthol biosynthesis pathway gene promoters from Mentha × piperita

Several cis-elements including Myb-binding motifs together confer glandular trichome specificity as revealed from heterologous expression and analysis of menthol biosynthesis pathway gene promoters. Glandular Trichomes (GTs) are result of division of epidermal cells that produce diverse metabolites. Species of mint family are important for their essential oil containing many high-value terpenoids, biosynthesized and stored in these GTs. Hence, GTs constitute attractive targets for metabolic engineering and GT-specific promoters are important. In this investigation, the upstream regions of the Mentha × piperita menthol biosynthetic pathway genes (-)-limonene synthase, (-)-P450 limonene-3- hydroxylase, (-)-trans-isopiperitenol dehydrogenase, (-)-Isopiperitenone reductase, ( +)-Pulegone reductase, (-)-Menthone reductase/ (-)-Menthol dehydrogenase and a branched pathway gene ( +)-menthofuran synthase were isolated and characterized. These fragments, fused to β-glucuronidase (GUS) reporter gene of pBI101 binary vector, are able to drive high level gene expression in transgenic tobacco trichomes with strong signals in GTs, except for (-)-Isopiperitenone reductase. The GT-enriched tissue from transformed plants were analysed for GUS enzyme activity and RNA expression which correlates the GUS staining. To characterize the cis-elements responsible for GT-specific expression, a series of 5' deletion constructs for MpPLS and MpPMFS were cloned and analysed in stable transgenic tobacco lines. The specificity of trichome expression was located to -  797 to-  598 bp sequence for (-)-limonene synthase and-  629 to -   530 bp for ( +)-menthofuran synthase promoters containing specific Myb-binding motifs in addition to other unique motifs described for developmental regulation without any defined pattern. All other pathway promoters also recruits specific but different Myb factors as indicated by this analysis.

Altered Developmental and Metabolic Gene Expression in Basil Interspecific Hybrids

To understand the altered developmental changes and associated gene expression in inter-genomic combinations, a study was planned in two diverse yet closely related species of Ocimum, targeting their hybrid F1 and amphidiploids. The existing developmental variations between F1 and amphidiploids was analyzed through phenotypical and anatomical assessments. The absence of 8330 transcripts of F1 in amphidiploids and the exclusive presence of two transcripts related to WNK lysine-deficient protein kinase and geranylgeranyl transferase type-2 subunit beta 1-like proteins in amphidiploids provided a set of genes to compare the suppressed and activated functions between F1 and amphidiploids. The estimation of eugenol and methyleugenol, flavonoid, lignin and chlorophyll content was correlated with the average FPKM and differential gene expression values and further validated through qRT-PCR. Differentially expressed genes of stomatal patterning and development explained the higher density of stomata in F1 and the larger size of stomata in amphidiploids. Gene expression study of several transcription factors putatively involved in the growth and developmental processes of plants clearly amalgamates the transcriptome data linking the phenotypic differences in F1 and amphidiploids. This investigation describes the influence of interspecific hybridization on genes and transcription factors leading to developmental changes and alleviation of intergenomic instability in amphidiploids.

Molecular Mechanisms of Plant Trichome Development

Plant trichomes, protrusions formed from specialized aboveground epidermal cells, provide protection against various biotic and abiotic stresses. Trichomes can be unicellular, bicellular or multicellular, with multiple branches or no branches at all. Unicellular trichomes are generally not secretory, whereas multicellular trichomes include both secretory and non-secretory hairs. The secretory trichomes release secondary metabolites such as artemisinin, which is valuable as an antimalarial agent. Cotton trichomes, also known as cotton fibers, are an important natural product for the textile industry. In recent years, much progress has been made in unraveling the molecular mechanisms of trichome formation in Arabidopsis thaliana, Gossypium hirsutum, Oryza sativa, Cucumis sativus, Solanum lycopersicum, Nicotiana tabacum, and Artemisia annua. Here, we review current knowledge of the molecular mechanisms underlying fate determination and initiation, elongation, and maturation of unicellular, bicellular and multicellular trichomes in several representative plants. We emphasize the regulatory roles of plant hormones, transcription factors, the cell cycle and epigenetic modifications in different stages of trichome development. Finally, we identify the obstacles and key points for future research on plant trichome development, and speculated the development relationship between the salt glands of halophytes and the trichomes of non-halophytes, which provides a reference for future studying the development of plant epidermal cells.

Transcriptomic and functional analysis provides molecular insights into multicellular trichome development

Trichomes, the hair-like structures located on aerial parts of most vascular plants, are associated with a wide array of biological processes and affect the economic value of certain species. The processes involved in unicellular trichome formation have been well-studied in Arabidopsis (Arabidopsis thaliana). However, our understanding of the morphological changes and the underlying molecular processes involved in multicellular trichome development is limited. Here, we studied the dynamic developmental processes involved in glandular and nonglandular multicellular trichome formation in cucumber (Cucumis sativus L.) and divided these processes into five sequential stages. To gain insights into the underlying mechanisms of multicellular trichome formation, we performed a time-course transcriptome analysis using RNA-sequencing analysis. A total of 711 multicellular trichome-related genes were screened and a model for multicellular trichome formation was developed. The transcriptome and co-expression datasets were validated by reverse transcription-quantitative PCR and in situ hybridization. In addition, virus-induced gene silencing analysis revealed that CsHOMEOBOX3 (CsHOX3) and CsbHLH1 are involved in nonglandular trichome elongation and glandular trichome formation, respectively, which corresponds with the transcriptome data. This study presents a transcriptome atlas that provides insights into the molecular processes involved in multicellular trichome formation in cucumber and can be an important resource for future functional studies.

Transcriptome mining of genes in Zanthoxylum armatum revealed ZaMYB86 as a negative regulator of prickly development

Zanthoxylum armatum DC. is an important economic tree species. Prickle is a type of trichome with special morphology, and there are a lot of prickles on the leaves of Z. armatum, which seriously restricts the development of Z. armatum industry. In this study, the leaves of Z. armatum cv. Zhuye (ZY) and its budding variety 'Rongchangwuci' (WC) (A less prickly mutant variety) at different developmental stages were used as materials, and the transcriptome sequencing data were analyzed. A total of 96,931 differentially expressed genes (DEGs) were identified among the samples, among which 1560 were candidate DEGs that might be involved in hormone metabolism. The contents of JA, auxin and CK phytohormones in ZY leaves were significantly higher than those in WC leaves. Combined with weighted gene co-expression network analysis, eight genes (MYC, IAA, ARF, CRE/AHK, PP2C, ARR-A, AOS and LOX) were identified, including 25 transcripts, which might affect the metabolism of the three hormones and indirectly participate in the formation of prickles. Combining with the proteins successfully reported in other plants to regulate trichome formation, ZaMYB86, a transcription factor of R2R3 MYB family, was identified through local Blast and phylogenetic tree analysis, which might regulate prickle formation of Z. armatum. Overexpression of ZaMYB86 in mutant A. thaliana resulted in the reduction of trichomes in A. thaliana leaves, which further verified that ZaMYB86 was involved in the formation of pickles. Yeast two-hybrid results showed that ZaMYB86 interacted with ZaMYB5. Furthermore, ZaMYB5 was highly homologous to AtMYB5, a transcription factor that regulated trichomes development, in MYB DNA binding domain. Taken together, these results indicated that ZaMYB86 and ZaMYB5 act together to regulate the formation of prickles in Z. armatum. Our findings provided a new perspective for revealing the molecular mechanism of prickly formation.

Glucosylation of (±)-Menthol by Uridine-Diphosphate-Sugar Dependent Glucosyltransferases from Plants

Menthol is a cyclic monoterpene alcohol of the essential oils of plants of the genus Mentha, which is in demand by various industries due to its diverse sensorial and physiological properties. However, its poor water solubility and its toxic effect limit possible applications. Glycosylation offers a solution as the binding of a sugar residue to small molecules increases their water solubility and stability, renders aroma components odorless and modifies bioactivity. In order to identify plant enzymes that catalyze this reaction, a glycosyltransferase library containing 57 uridine diphosphate sugar-dependent enzymes (UGTs) was screened with (±)-menthol. The identity of the products was confirmed by mass spectrometry and nuclear magnetic resonance spectroscopy. Five enzymes were able to form (±)-menthyl-β-d-glucopyranoside in whole-cell biotransformations: UGT93Y1, UGT93Y2, UGT85K11, UGT72B27 and UGT73B24. In vitro enzyme activity assays revealed highest catalytic activity for UGT93Y1 (7.6 nkat/mg) from Camellia sinensis towards menthol and its isomeric forms. Although UGT93Y2 shares 70% sequence identity with UGT93Y1, it was less efficient. Of the five enzymes, UGT93Y1 stood out because of its high in vivo and in vitro biotransformation rate. The identification of novel menthol glycosyltransferases from the tea plant opens new perspectives for the biotechnological production of menthyl glucoside.

Transcriptome analysis of Clinopodium gracile (Benth.) Matsum and identification of genes related to Triterpenoid Saponin biosynthesis

BACKGROUND

Clinopodium gracile (Benth.) Matsum (C. gracile) is an annual herb with pharmacological properties effective in the treatment of various diseases, including hepatic carcinoma. Triterpenoid saponins are crucial bioactive compounds in C. gracile. However, the molecular understanding of the triterpenoid saponin biosynthesis pathway remains unclear.

RESULTS

In this study, we performed RNA sequencing (RNA-Seq) analysis of the flowers, leaves, roots, and stems of C. gracile plants using the BGISEQ-500 platform. The assembly of transcripts from all four types of tissues generated 128,856 unigenes, of which 99,020 were mapped to several public databases for functional annotation. Differentially expressed genes (DEGs) were identified via the comparison of gene expression levels between leaves and other tissues (flowers, roots, and stems). Multiple genes encoding pivotal enzymes, such as squalene synthase (SS), or transcription factors (TFs) related to triterpenoid saponin biosynthesis were identified and further analyzed. The expression levels of unigenes encoding important enzymes were verified by quantitative real-time PCR (qRT-PCR). Different chemical constituents of triterpenoid saponins were identified by Ultra-Performance Liquid Chromatography coupled with quadrupole time-of-flight mass spectrometry (UPLC/Q-TOF-MS).

CONCLUSIONS

Our results greatly extend the public transcriptome dataset of C. gracile and provide valuable information for the identification of candidate genes involved in the biosynthesis of triterpenoid saponins and other important secondary metabolites.

Updates on molecular mechanisms in the development of branched trichome in Arabidopsis and nonbranched in cotton

Trichomes are specialized epidermal cells and a vital plant organ that protect plants from various harms and provide valuable resources for plant development and use. Some key genes related to trichomes have been identified in the model plant Arabidopsis thaliana through glabrous mutants and gene cloning, and the hub MYB-bHLH-WD40, consisting of several factors including GLABRA1 (GL1), GL3, TRANSPARENT TESTA GLABRA1 (TTG1), and ENHANCER OF GLABRA3 (EGL3), has been established. Subsequently, some upstream transcription factors, phytohormones and epigenetic modification factors have also been studied in depth. In cotton, a very important fibre and oil crop globally, in addition to the key MYB-like factors, more important regulators and potential molecular mechanisms (e.g. epigenetic modifiers, distinct metabolic pathways) are being exploited during different fibre developmental stages. This occurs due to increased cotton research, resulting in the discovery of more complex regulation mechanisms from the allotetraploid genome of cotton. In addition, some conservative as well as specific mediators are involved in trichome development in other species. This study summarizes molecular mechanisms in trichome development and provides a detailed comparison of the similarities and differences between Arabidopsis and cotton, analyses the possible reasons for the discrepancy in identification of regulators, and raises future questions and foci for understanding trichome development in more detail.

Ocimum metabolomics in response to abiotic stresses: Cold, flood, drought and salinity

Ocimum tenuiflorum is a widely used medicinal plant since ancient times and still continues to be irreplaceable due to its properties. The plant has been explored chemically and pharmacologically, however, the molecular studies have been started lately. In an attempt to get a comprehensive overview of the abiotic stress response in O. tenuiflorum, de novo transcriptome sequencing of plant leaves under the cold, drought, flood and salinity stresses was carried out. A comparative differential gene expression (DGE) study was carried out between the common transcripts in each stress with respect to the control. KEGG pathway analysis and gene ontology (GO) enrichment studies exhibited several modifications in metabolic pathways as the result of four abiotic stresses. Besides this, a comparative metabolite profiling of stress and control samples was performed. Among the cold, drought, flood and salinity stresses, the plant was most susceptible to the cold stress. Severe treatments of all these abiotic stresses also decreased eugenol which is the main secondary metabolite present in the O. tenuiflorum plant. This investigation presents a comprehensive analysis of the abiotic stress effects in O. tenuiflorum. Current study provides an insight to the status of pathway genes’ expression that help synthesizing economically valuable phenylpropanoids and terpenoids related to the adaptation of the plant. This study identified several putative abiotic stress tolerant genes which can be utilized to either breed stress tolerant O. tenuiflorum through pyramiding or generating transgenic plants.

Effect of Heat Stress on Yield, Monoterpene Content and Antibacterial Activity of Essential Oils of Mentha x piperita var. Mitcham and Mentha arvensis var. piperascens

Heat stress affects the yield of medicinal plants and can reduce biomass and/or metabolite production. In order to evaluate the effect of heat-induced stress on the essential oil production in Mentha x piperita L. var. Mitcham (Mitcham mint) and Mentha arvensis var. piperascens Malinv. ex L. H. Bailey (Japanese mint), we studied the chemical composition of the oils of the two mint species under different heat shock stresses in growth chambers. The antibacterial activity of the essential oils was also evaluated; microscopic observation (fluorescence and electron transmission) was used to assess the effect of the tested samples on bacterial growth. The results obtained shed light on the mint essential oils composition and biological activity in relation to heat stress.

Regulation of trichome development in tobacco by JcZFP8, a C2H2 zinc finger protein gene from Jatropha curcas L

Trichomes are epidermal outgrowths of plant tissues that can secrete or store large quantities of secondary metabolites, which contribute to plant defense responses against stress. The use of bioengineering methods for regulating the development of trichomes and metabolism is a widely researched topic. In the present study, we demonstrate that JcZFP8, a C2H2 zinc finger protein gene from Jatropha curcas L., can regulate trichome development in transgenic tobacco. To understand the underlying mechanisms, we performed transcriptome profiling of overexpression JcZFP8 transgenic plants and wild-type tobacco. Based on the analysis of differentially expressed genes, we determined that genes of the plant hormone signal transduction pathway was significantly enriched, suggesting that these pathways were modulated in the transgenic plants. In addition, the transcript levels of the known trichome-related genes in Arabidopsis were not significantly changed, whereas CycB2 and MYB genes were differentially expressed in the transgenic plants. Despite tobacco and Arabidopsis have different types of trichomes, all the pathways were associated with C2H2 zinc finger protein genes. Our findings help us to understand the regulation of multicellular trichome formation and suggest a new metabolic engineering method for the improvement of plants.

De novo transcriptomic analysis of leaf and fruit tissue of Cornus officinalis using Illumina platform

Cornus officinalis is one of the most widely used medicinal plants in China and other East Asian countries to cure diseases such as liver, kidney, cardiovascular diseases and frequent urination for thousands of years. It is a Level 3 protected species, and is one of the 42 national key protected wild species of animals and plants in China. However, the genetics and molecular biology of C. officinalis are poorly understood, which has hindered research on the molecular mechanism of its metabolism and utilization. Hence, enriching its genomic data and information is very important. In recent years, the fast-growing technology of next generation sequencing has provided an effective path to gain genomic information from nonmodel species. This study is the first to explore the leaf and fruit tissue transcriptome of C. officinalis using the Illumina HiSeq 4000 platform. A total of 57,954,134 and 60,971,652 clean reads from leaf and fruit were acquired, respectively (GenBank number SRP115440). The pooled reads from all two libraries were assembled into 56,392 unigenes with an average length 856 bp. Among these, 41,146 unigenes matched with sequences in the NCBI nonredundant protein database. The Gene Ontology database assigned 24,336 unigenes with biological process (83.26%), cellular components (53.58%), and molecular function (83.93%). In addition, 10,808 unigenes were assigned a KOG functional classification by the KOG database. Searching against the KEGG pathway database indicated that 18,435 unigenes were mapped to 371 KEGG pathways. Moreover, the edgeR database identified 4,585 significant differentially expressed genes (DEGs), of which 1,392 were up-regulated and 3,193 were down-regulated in fruit tissue compared with leaf tissue. Finally, we explored 581 transcription factors with 50 transcription factor gene families. Most DEGs and transcription factors were related to terpene biosynthesis and secondary metabolic regulation. This study not only represented the first de novo transcriptomic analysis of C. officinalis but also provided fundamental information on its genes and biosynthetic pathway. These findings will help us explore the molecular metabolism mechanism of terpene biosynthesis in C. officinalis.

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.