Comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri identifies potential genes related to triterpenoid saponin biosynthesis
- Gajendra Singh Jeena1,
- Shahnoor Fatima†1,
- Pragya Tripathi†1,
- Swati Upadhyay†1 and
- Rakesh Kumar Shukla1Email author
© The Author(s). 2017
Received: 4 March 2017
Accepted: 15 June 2017
Published: 28 June 2017
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.
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.
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.
Bacopa monnieri (L.) Wettst. Commonly known as “Brahmi” is a medicinal herb of great significance. It has been used worldwide due to its broad use of pharmaceutically important triterpenoid saponins mainly bacosides . It is a perennial, creeping herb basically found in the wetlands of Australia, southern India, Europe, Asia, Africa, and South America. Bacopa is a genus of 70–100 aquatic plants of family Scrophulariaceae (Plantaginaceae) having succulent, oblong and 4–6 mm thick leaves arranged oppositely on the stem.
In the conventional system of Indian medicine, Bacopa is well treated as a drug to augment intelligence and memory function and combat the effects of mental stress . It is also used for gastrointestinal infections, rejuvenation, skin disorders, epilepsy, pyrexia and analgesia . Apart from that Bacopa monnieri extract is familiar to possess potent antioxidant and anticancer properties . Bacpa monnieri has many efficacious compounds in addition to alkaloids, flavonoids, betulic acid, stigmasterol, beta-sitosterol and saponins (bacoside A, bacoside B, bacopasaponin C, bacopaside I, bacopaside II, bacopaside X, bacopaside N2). Due to its affluent active components, the plant has been largely utilized as a nootropic digestive attention and for improving respiratory functions . Earlier studies also revealed that bacosides from Bacopa monnieri leaf extract is used to treat memory functions and impairment . Evaluation and extraction of bacosides from areal part of the plant using HPLC chromatographic technique have been performed and reported that highest concentration of bacoside A is present in stolon (9.54 mg/g dry wt) followed by leaves (4.73 mg/g dry wt) and roots .
In bacteria and higher eukaryotes MVA pathway is one of the important metabolic pathways. The product of MEP and MVA pathway such as 3,3-dimethylallyl diphosphate (DMAPP) and isopentenyl diphosphate (IPP) act as the important mediator in several secondary metabolites production in plants [8, 9]. Several enzymes take part in structural alteration of these intermediates by undergoing substitution, oxidation, and glycosylation producing various glycosylated triterpenoids. Squalene synthase which is an important governing step in sterol and triterpenoid biosynthesis catalyzes the conversion of farnesyl pyrophosphate (FPP) molecules to squalene. Squalene epoxidase is involved in the epoxidation of squalene to form 2,3-oxidosqualene which is the branch point in the biosynthesis of steroidal and terpenoid sapogenin. Cyclization of 2,3-oxidosqualene catalyzed by a class of oxidosqualene cyclases (OSCs) to form distinct secondary metabolite backbones. Later cytochrome P450-dependent monooxygenases (CYP450) and glycosyltransferases (GTs) govern oxidation, hydroxylation, and glycosylation step to yield triterpenoid saponins and sterols with distinct structures and biological activities .
Emerging high-throughput sequencing technologies has changed the pace of DNA sequencing in plants and animals which perform comprehensive profiling of RNA expression [11, 12] other than model plant organisms were constrained molecular genetics applications have been implemented. RNA sequencing administers complete transcriptome expression profiles of selected plant cells or tissues, thereby allowing novel approaches to analyze functional genomics, in addition more precise measurement of transcript level and their isoforms than other methods. Transcriptome assay using NGS sequencing has been extensively used to identify genes encoding enzymes ramified in distinct steps of biosynthetic pathways in medicinal plants. Some include the description of genes encoding enzyme catalyzed metabolic steps mediated in the biosynthetic pathway of artemisinin in Artemisia annua [13, 14], withanolides in Withania somnifera [15, 16], cannabinoids in Cannabis sativa , ginsenosides in Panax ginseng , glycyrrhizin in Glycyrrhiza uralensis , picrosides in Picrorhiza kurrooa , cardiac glycoside in Calotropis procera , menthol content in Mentha species , medecinal diterpene in Andrographis paniculata  and biosynthesis of steroidal saponins in Asparagus racemosus .
As the synthesis and accumulation of specialized plant metabolites in different tissues depends on the age of the plant and is also highly influenced by the different developmental stages. It was reported earlier that in P. ginseng and P. quinquefolius, leaves accumulate higher content of triterpenoids during earlier growth stages whereas older plants have a higher content of triterpenoids in root . Comparative transcriptome analysis including root and shoot tissues to identify transcripts involved in saponin biosynthesis is already reported for many plant species like Asparagus racemosus , Panax notoginseng  Gynostemma pentaphyllum . The maximum bacoside content was identified in stolon followed by leaf and usually the aerial portion extract is used for checking its efficacy in memory impairment responses. So in this study comparative shoot versus root paired-end transcriptome sequencing was performed for Bacopa monnieri to identify genes involved in triterpenoid saponin. Efforts were also made to identify potential transcription factors, CYP450s and UGTs that play a key role in regulation and diversification of secondary metabolites. Paired-end sequencing strategy extends the mapped fragment length to 200–500 bp so it is expected to be useful in obtaining longer reads, allows for detecting alternate splice junctions, deletions, insertions, and is useful for de novo transcriptome assembly . Transcripts generated after assembly were functionally elucidated in different gene ontology, some of the secondary metabolic pathway related transcripts were selected for further validation. The fold change induction of differentially expressed transcripts in the shoot and root tissue was evaluated using DGE (Digital gene expression) profile, followed by q-RT PCR analysis. The transcriptome data generated in this study will provide new insight into the enzyme transcripts involved in the secondary metabolic pathway of this important medicinal plant.
Transcriptome sequencing and de novo assembly
Summary of HI-Seq Illumina Transcriptome assembly of Bacopa monnieri root and shoot tissue after two independent biological replicate
Number of transcripts identified
Maximum Contig Length
Minimum Contig Length
Average Contig Length
999.7 ± 826.2
931.5 ± 793.6
Median Contig Length
Total Contigs Length
Total Number of Non-ATGC Characters
Percentage of Non-ATGC Characters
Contigs > = 100 bp
Contigs > = 200 bp
Contigs > = 500 bp
Contigs > = 1 Kbp
Contigs > = 10 Kbp
Functional characterization and GO analysis
Mapping and enrichment of KEGG pathway
The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway provides better insight about pathway maps revealing the knowledge on the molecular interaction and reaction networks. To better understand about biological pathways operating in Bacopa monnieri, a BLASTx search against the KEGG protein database was performed on the assembled unigenes with an E-value threshold of 10–5. Pathway analysis was done by using KAAS Server. Arabidopsis thaliana were considered as reference organisms for pathway identification. Among 85,579 annotated unigenes in root and 74,624 unigenes in the shoot, 37,918 and 35,130 unigenes were mapped to 133 KEGG pathways in root and shoot respectively. Based upon the number of unigenes assigned, carbohydrate metabolism, translation, lipid metabolism, amino acid metabolism and signal transduction are the dominant metabolic pathways operating in shoot and root. Among secondary metabolic pathways, the most represented category was phenylpropanoid biosynthesis (365 unigenes), terpenoid backbone biosynthesis (340 unigenes), carotenoid biosynthesis (125 unigenes), diterpenoid biosynthesis (77 unigenes) and zeatin biosynthesis (57 unigenes) in root followed by terpenoid backbone biosynthesis (358 unigenes), phenylpropanoid biosynthesis (286 unigenes), carotenoid biosynthesis (126 unigenes), diterpenoid biosynthesis (52 unigenes) and zeatin biosynthesis (51 unigenes) in shoot.
Differential expression analysis of unigenes
Summary of DGE data obtained after two independent biological replicate and its statistical analysis
Total Unigenes available
No. of unigenes in each category for Root and Shoot samples
No. of unigenes expressed in both samples
No. of unigenes expressed only in Root sample
No. of unigenes expressed only in Shoot sample
No. of P-significant unigenes
No. of Q-significant unigenes
Software used for carrying out DGE
Sample taken as control for DGE
Sample taken as Treated for DGE
Mevalonic acid pathway and transcripts related to Saponin biosynthesis
List of transcripts found in DGE data related to Steroidal and Terpenoid sapogenins biosynthesis obtained after two independent biological replicates along with their mean values
Expression value (Root)
Expression value (shoot)
Acetyl-CoA C acetyltransferase [EC:220.127.116.11]
Hydroxymethyl gutryl CoA synthase [EC:18.104.22.168]
Hydroxymethyl gutryl CoA reductase [EC:22.214.171.124]
Mevalonate kinase [EC:126.96.36.199]
Phosphomevalonate kinase [EC: 188.8.131.52]
Diphosphomevalonate decarboxylase [EC: 184.108.40.206]
Geranylgeranyl diphosphate synthase [EC:220.127.116.11]
Isopentenyl pyrophosphate [EC:18.104.22.168]
Farnesyl diphosphate synthase [EC:22.214.171.124.126.96.36.199]
Squalene synthase [EC:188.8.131.52]
Squalene monooxygenase [EC:184.108.40.206]
Cycloartenol synthase [EC:220.127.116.11]
Sterol 24-methyltransferase [EC:18.104.22.168]
Cycloeucalenol cycloisomerase [EC:22.214.171.124]
Obtusifoliol 14α demethylase [EC:126.96.36.199]
Cholestenol delta isomerise [EC:188.8.131.52]
Beta-amyrin 28 oxidase [EC:184.108.40.206]
Cytochrome P450 [EC:1.14.]
1-deoxy-D-xylulose 5-phosphate reductoisomerase [EC:220.127.116.117]
1-deoxy-D-xylulose-5-phosphate synthase [EC:18.104.22.168]
Study of metabolic pathway genes
Putative transcripts involved in Phenylpropanoid biosynthesis pathway
Most of the known transcripts associated with the phenylpropanoid and monolignoid pathway were found to be clearly expressing in root in comparison with shoot tissue of Bacopa monnieri. In our transcriptome data, we found that 286 unigenes in the shoot and 365 unigenes in root tissue were specifically involved in phenylpropanoid biosynthetic pathway. Chalcone synthase (24 unigenes, EC 22.214.171.124) which is an important primary enzyme in the phenylpropanoid pathway was found to be highly expressed in root tissue leads to the production of various secondary metabolites being flavonoids, isoflavonoids, and anthocyanins. The expression level of cinnamoyl alcohol dehydrogenase (CAD) was detected and was found upregulated in root tissue which is assumed to be involved in monolignoid biosynthesis.
Transcription factor families
Functionally annotated transcripts validation by using qRT-PCR analysis
Analysis of secondary metabolic pathway related transcripts by quantitative real-time PCR (q RT-PCR) after wounding and methyl jasmonate (MeJA) treatment
MeJA acts as an important elicitor in a wide spectrum of signaling pathways and also regulates the activity of distinct enzymes involved in the secondary metabolic pathways. MeJA has been proved to elicit the production of several compounds such as alkaloids, terpenoid and phenolic phytoalexins, coumarins, and taxanes in many plant species . Wounding is a common damage that occurs to plants as a result of abiotic and biotic stress. It is hypothesized that plants may have evolved several mechanisms in response to wounding. We have selected 43 transcripts including transcription factor families from shoot and root related to secondary metabolite biosynthesis specifically to triterpenoid and phenylpropanoid pathway and also monitor their relative expression after 1, 3 and 5 h of exogenous application of MeJA treatment and wounding in tissues respectively. All the transcripts and transcription factor families showed increased accumulation after 3 h of MeJA treatment (Additional file 8A, B). In Arabidopsis, MeJA was found to activate both the general and downstream aspects of the phenylpropanoid biosynthesis pathway . To be steady with Arabidopsis, the expression of most phenylpropanoid biosynthesis genes in Bacopa monnieri was upregulated by MeJA in different degree (Fig. 4b).
B. monnieri is an important medicinal plant with the number of pharmacological uses. Although it is pharmacologically important but the genome or transcriptome information is very confined. In NCBI only 68 ESTs, 68 nucleotide sequences and 35 protein sequences are accessible from Bacopa monnieri. The present manuscript has studied the comparative de novo transcriptome analysis of shoot and root tissue of Bacopa monnieri. The identified transcript from shoot and root will further be useful to understand the biosynthetic pathway of pharmacologically useful secondary metabolite present in different tissues of B. monnieri.
De novo comparative shoot and root transcriptome of Bacopa monnieri generates more than 44.68 million high-quality reads which were further used for assembly of data. The mean length of unigenes in the resulting assembly was more than 1000 bp which was found to be longer than the result showed in earlier studies [33, 34]. We reflect the pace of high-quality sequencing by bringing about more than 100% high-quality bases for both shoot and root tissues. Assembly process could be hampered by the presence of low-quality reads including adapters resulting in misassemblies or degraded contigs. Consequently, low-quality bases including adapter sequences were removed before continuing with further assembly. Our results proposed that the de novo transcriptome assembly from Bacopa monnieri were assembled correctly, which was further monitored by the high fraction of unigenes coordinated with known proteins along with higher N50 value.
The best match for individual unigene search against the NCBI non-redundant database and KEGG was promoted to assign functional GO annotation under biological process, cellular component groups and molecular function. The abundance of varied GO assignments to unigenes featured the assortment of genes possibly represented in Bacopa monnieri shoot and root transcriptome data. Unigenes were mapped onto KEGG, we had found that several unigenes are involved in distinct secondary metabolites biosynthesis. Unigenes without the match are most likely related to untranslated regions, noncoding RNA and short sequence protein without any domain or assembly mistakes. Moreover, the sequences encoding a large range of transcripts characterized under the cellular component category which points out the requirement of various transcripts for maintenance of the cellular structure. In support of the annotation, it reveals that the transcripts encoding a large number of enzymes involved in triterpenoid backbone biosynthesis specifically occur in roots.
Bacosides are the main constituent of Bacopa monnieri that was chemically a triterpenoid saponin and the pharmacological properties are mainly attributed to these triterpenoid saponin compounds present in the plant extract . Three new dammarane-type triterpenoid saponins i.e. bacopasaponins A, B and C alongwith pseudojujubogenin (bacopasaponin D) were identified by spectroscopic and chemical transformation methods . In our transcriptome, we have found transcripts related to both steroidal and terpenoid saponin which suggests the presence of both types of saponins in this plant. Transcripts related to bacosides biosynthesis like β amyrin synthase, CYP450 monooxygenases, GTs predominates in root tissue as compared to shoot but expression analysis of differentially expressed transcripts related to CYP450 and UGTs showed that they are highly expressed in the shoot as compared to root tissue. This implies that later stage modification of saponin backbone that involves several hydroxylation and glycosylation steps takes place in shoot tissue. In addition to saponin, several other transcripts involved in Phenylpropanoid biosynthesis were also found in our transcriptome data which suggests the presence of some flavonoids, lignin and phenylpropanoids compounds that was not reported earlier in Bacopa monnieri.
MeJA acts as an elicitor and plays critical roles in various plants metabolic processes by enhancing the expression of genes associated with secondary metabolic pathways. On the other hand wounding generally, induces the pathways genes related to secondary metabolites involved in defense response . Transcript of Cinnamyl alcohol dehydrogenase (CAD) from Bacopa monnieri was showing induced expression in response to MeJA and wounding in shoot tissue which is involved in monolignoid biosynthesis. Similarly wounding in Silanum elaeagnifolium showed significant increased stress response Caryophyllene, Geranyl linalool and two terpenes transcripts .
Our result showed differential expression of transcripts involved in secondary metabolism in both root and shoot of Bacopa monnieri. Transcripts of Cinnamoyl CoA dehydrogenase, Caffeoyl-CoA methyltransferase, HMG-CoA reductase, Squalene monooxygenase and Isoflavone 2-hydroxylase of Bacopa monnieri showed induced expression after MeJA treatment in shoot specific manner. In root caffeoyl-CoA-O-methyltransferase, diphospho mevalonate decarboxylase, β-amyrin synthase, HMG- CoA-synthase showed induced expression in a root-specific manner. Identified and selected transcripts related to saponin and other secondary metabolic pathway related transcripts were showing better response towards MeJA treatment in shoot tissues. The increased expression of an important enzyme of triterpenoid backbone biosynthesis after MeJA treatment in shoots such as HMGCR, 5-Phospho mevalonate kinase, Mevalonate 5-pyrophosphate decarboxylase, Isopentenyl 5-pyrophosphate isomerase, Squalene synthase and Squalene epoxidase further showed that primary reactions of isoprenoid backbone biosynthesis take place in the shoot. The enhanced expression of β-amyrin synthase in response to MeJA and wounding in root tissue of Bacopa monnieri which is an important enzyme involved in triterpenoid sapogenin biosynthesis at later stages further reveals that triterpenoid saponin biosynthesis might be enhanced when the plant recognizes certain elicitor under stressed conditions. Similar stimulatory effects of MeJA on the biosynthetic pathway of other triterpenoid saponin compounds have already been reported . It was also showed that triterpenoid biosynthetic enzyme β-amyrin synthase increased up to 50-fold by introducing MeJA to cell suspension cultures of M. truncatula . Illumina paired-end RNA-seq for de novo reconstruction of P. minus leaf transcriptome to identify differentially expressed genes in response to MeJA elicitation .
Transcription factors interact with the promoter region of the downstream gene and regulate its expression. These are the site of single step manipulation in the complex metabolic pathway as they can regulate simultaneously a number of steps in a pathway. During transcriptome analysis, we have identified the specific transcription factors which are expressing in either shoot or root. The real-time expression analysis of these transcription factors in response to MeJA and wounding has further identified NAC and bHLH transcription factor showing the better response in comparison with other families such as WRKY and AP2/ERF protein. Similarly, dynamic changes in the expression pattern were observed in different transcription factor families in Populus hopeiensis after ABA treatment . Also, the transcriptome based tissue-specific expression of transcription factor genes provides important information for understanding the development and transcriptional regulation of the paper mulberry . The transcriptional regulatory responses of these transcripts after MeJA treatment and wounding would be useful further to identify transcripts responsive to elicitation. The above information can be utilized as an effective component through genetic engineering for enhancing the metabolite content.
In the present study, the comparative transcriptome analysis of shoot and root tissues of Bacopa monnieri was performed to know about the putative transcripts involved in the secondary metabolite biosynthetic pathway of an important medicinal herb. The differential expression pattern of pathway genes and transcription factors in the shoot and root tissues of Bacopa monnieri suggest tissue-specific synthesis, accumulation followed by modification of these metabolites might occur by different signaling cascades. This study provides a useful resource of biosynthetic pathway transcripts identified from this important Ayurvedic medicinal plant with human health benefits. The identified transcripts will pave the foundation to metabolically engineer this plant with pharmacologically increased metabolite content.
Sampling and RNA extraction
For transcriptome analysis wild type Bacopa monnieri (CIM-JAGRITI) plants were collected from the green house of CSIR-CIMAP (Lucknow). Samples of the shoot and root tissues were collected from the 1 month old plant and immediately chilled in liquid nitrogen and stored at −80 °C. Total RNA from approximately 100 mg of frozen tissue of both shoot and root was extracted using TRIzol (Invitrogen) following the manufacturer’s recommendations.
cDNA library preparation and sequencing
Illumina multiplexed Sequencing mRNA library preparation was carried out at Genotypic Technology’s Genomics Facility using Sure-select strand-specific RNA library preparation kit following manufacture’s protocol (Cat # 5500–0116). Briefly, the mRNA was fragmented for 4 min at elevated temperature (94oC) in the presence of divalent cations and first strand cDNA is synthesized. Strand specificity is maintained by the addition of actinomycin D. Second strand cDNA was synthesized using second strand synthesis mix. The cDNA was cleaned by AgencourtAmpure XP SPRI beads (Beckman-Coulter). After end repair and addition of “A” base, adapters were ligated to the cDNA molecules. After ligation SPRI cleanup was performed. The library was amplified and indexed using 10 cycles of PCR for improvement of adapter-ligated fragments. Finally, the concentration and size of the cDNA library were determined with a Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, US) and an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).
De novo assembly and bioinformatics analysis of sequences
After quality determination, sequencing was performed by using paired-end Illumina Hi-sequencing 2000. Quality control programme i.e. SeqQC-V2.2 was used to determine the aspects of sequence reads. Reads with adaptor and vector contamination was discarded. Finally, the reads with more than 20% Q < 20 bases were also removed. The high quality concatenated reads of shoots and roots were used to develop de novo assembly by De Bruijn graph algorithm based software tools. Reads were massed into contigs at distinct k-mer value using Velvet (https://www.ebi.ac.uk/~zerbino/velvet/). After these resulting contigs were merged into transcript isoforms using Oases_v0.2.08 (https://www.ebi.ac.uk/~zerbino/oases/), which was specially designed for the de novo assembly of transcripts using short reads. Different hash lengths were used to select the best de novo assembly. Construction of transcript isoforms was done by using paired-end information. Different specifications of assembly such as N50 length, average contig length, the number of contigs were upgraded to retrieve the best-assembled data with high coverage length. Combined transcripts from the multiple k-mer assemblies were run through the CD-HIT (CD-HIT v4.5.4 http://www.bioinformatics.org/cd-hit/) program to remove redundant transcripts sharing 100% identity.
The putative function of unigenes obtained in transcriptome was subjected to BLASTx analysis against the non-redundant protein sequences (nr) of PDB, UniProt, PIR, all annotated protein sequences of Arabidopsis (available at the Arabidopsis information resource) and Oryza sativa. The results of only the best hit were extracted and the hits with an E-value <1E-05 were considered to be significant. Uniprot, Swiss-Prot, and TrEmBL databases were utilized for homology search. Gene ontology (GO) annotations were made by using Blast2GO program that was based on Nr annotations of NCBI and unigenes were categorized into three different GO terms i.e. molecular function, cellular component, and biological process. After assigning Enzyme Commission numbers (ECs) to unigenes from Fast Annotator, KEGG (Kyoto Encyclopedia of Genes and Genomes) was used to determine metabolic pathways related with them. Bidirectional best hit information was achieved, related to the metabolic pathway with the help of KAAS (KEGG Automatic Annotation Server http://www.genome.jp/tools/kaas/) by using reference database of model plant Arabidopsis thaliana (thale cress) and Oryza sativa japonica (Japanese rice).
Identification and analysis of differentially expressed genes (DEGs)
The definite algorithm was developed to identify genes differentially transcribed in the shoot and root tissue of Bacopa monnieri. Both leaf and root samples were combined and clustered using CD-hit. Reads of both the samples were used against unigenes using Bowtie-2 to obtained read count profile. DESeq was used to identify the differential gene expression. RPKM values were adjusted with RNA and total read number to normalize transcript molar concentration and length. To compare the RPKM of DEGs we use false discovery rate (FDR), a statistical method to test the correction for comparisons. We obtained the p-value using the differential gene expression test. FDR adjustment was used to determine the p-value threshold in multiple tests and analyses. When a single gene is represented by two or more unigenes, the longest unigenes among them was selected to check the relative expression level and its coverage. For genes with missing values in a specific sample, the RPKMs were adjusted to 0.001. To identify differentially expressed genes, more precise criteria, with smaller FDR and bigger fold-change values, were used. Functional correlation can be observed within genes with similar expression patterns. For each of these analyses, a P value <0.05 was required for differences to be considered statistically significant.
Transcription factor analysis
Transcripts coding for transcription factors in the transcriptome of Bacopa monnieri shoot and root tissues were anticipated by searching against all transcription factor protein sequences present in PlantTFcat (http://plantgrn.noble.org/PlantTFcat/about.gy) using BLASTX with an E-value cutoff 1e−6 and domain analysis within transcription factor was performed by using InterPro.
Quantitative real-time PCR (qRT-PCR) analysis
Twenty-three transcripts and 20 transcription factors which are either highly up-regulated or tissue-specific in shoot and root with inherent aspects in secondary metabolite biosynthesis were selected for validation by qRT-PCR (list of primers for qRT-PCR analysis is given in Additional file 10: Table S1). For control sample, 2-month-old seedlings of Bacopa monnieri plants were harvested, shoot and root tissues were separated, quickly wiped, frozen in liquid nitrogen and stored at −80 °C for further use. Total RNA was isolated from the shoot and root tissue separately by GITC method. To study the effect of MeJA and wounding on the expression pattern of randomly selected 23 potential transcripts related to secondary metabolite biosynthesis, and 20 transcripts related to transcription factors were selected. Green house grown 1-month-old seedlings of Bacopa monnieri were used for MeJA treatment and wounding. MeJA (250 μM) solution in DMSO and Triton-X was prepared and used for spraying over aerial parts of healthy plants. After spraying the solution the seedlings were covered by perforated transparent plastic bags to maintain proper transpiration. Wounding treatment was given to the plantlets by puncturing the whole plant on leaves and stem. Samples of MeJA and Wound treated plants were collected after 1, 3 and 5 h of treatment separately and washed with MQ to remove soil and other contaminants. Total RNA from the shoot and root samples were extracted separately by GITC method. The RNA samples were dissolved in MQ and examined using the NanoDrop ND-8000. cDNA synthesis was done by using 2 μg of RNA from the shoot and root tissues of treated as well as control plant with high capacity cDNA Reverse Transcription Kit (Applied Biosystems). Diluted cDNA samples were utilized as a template for qRT-PCR analysis. The reaction was carried out on a 7500 FAST Real-Time PCR System (Applied Biosystems, USA) using the RealMasterMix (SYBR Green, Takara). Actin and ubiquitin were used as endogenous control for normalization of expression levels of unigenes. Relative gene expression levels were calculated using the 2-ΔΔCt.
The authors acknowledge the national gene bank of Central Institute of Medicinal and Aromatic Plants (CIMAP) Lucknow for providing Bacopa monnieri seedlings. RKS acknowledge Director CSIR-CIMAP for providing necessary facilities. Gajendra, Swati, acknowledge CSIR and UGC for fellowship.
Partly from INSA young scientist project (GAP-274) from Indian National Science Academy, INDIA.
Availability of data and materials
Final sequences obtained from both the tissues were submitted in the SRA database of NCBI with accession number SAMN04216962.
GJ has performed the transcriptome analysis and did the real time analysis to validate them, SF and PT has also performed the real time analysis with GJ. SU has helped GJ in transcriptome analysis. GJ has prepared the manuscript and RKS has edited the manuscript and supervised the work. All authors read and approved the final manuscript.
Biotechnology Division, Central Institute of Medicinal and Aromatic Plants (CSIR-CIMAP), Kukrail Picnic Spot Road, P.O. CIMAP, Lucknow 226015, India.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Darokar MP, Khanuja Suman PS, Shasany AK, Kumar S. Low levels of genetic diversity detected by RAPD analysis in geographically distinct accessions of Bacopa monnieri. Genet Resourc Crop Evol. 2001;48:555–8.View ArticleGoogle Scholar
- Sivaramakrishna C, Rao CV, Trimurtulu G. Triterpenoid glycosides from Bacopa monnieri. Phytochemistry. 2005;66:2719–28.View ArticlePubMedGoogle Scholar
- Aguiar S, Borowski T. Neuropharmacological review of the nootropic herb Bacopa monnieri. Rejuvenation Res. 2013;16:313–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Elangovan V, Govindaswamy S, Ramamoorthy N, Balasubramanian K. In vitro studies on the anticancer activity of B. monniera. Fitoterapia. 1995;66:211–5.Google Scholar
- Russo A, Borrelli F. Bacopa monnieri, a reputed nootropic plant: an overview. Phytomedicine. 2005;12:305–17.View ArticlePubMedGoogle Scholar
- Hota SK, Barhwal K, Baitharu I, Prasad D, Singh SB, Ilavazhagan G. Bacopa monniera leaf extract ameliorates hypobaric hypoxia induced spatial memory impairment. Neurobiol Dis. 2009;34(1):23–39.View ArticlePubMedGoogle Scholar
- Naik PM, Manohar SH, Praveen N, Upadhya V, Murthy NH. Evaluation of Bacoside a content in different accessions and various organs of Bacopa monnieri (L.) Wettst. J Herbs Spices Med Plants. 2012;18:387–95.View ArticleGoogle Scholar
- Verpoorte R, Van der Heijden R, Moreno PRH. Biosynthesis of terpenoid indole alkaloids in Catharanthus roseus cells In: Cordell, GA. editor. The alkaloids. 1997;Vol.49:San Diego: Academic Press;221–299.Google Scholar
- Contin A, Collu G, Van der Heijden R, Verpoorte R. The effects of phenobarbital and ketoconazole on the alkaloid biosynthesis in Catharanthus roseus cell suspension cultures. Plant Physiol Biochem. 1999;37:139–44.View ArticleGoogle Scholar
- Niu Y, Luo H, Sun C, Yang TJ, Dong L, Huang L, et al. Expression profiling of the triterpene saponin biosynthesis genes FPS, SS, SE, and DS in the medicinal plant Panax notoginseng. Gene. 2014;533:295–303.View ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M. RNA-Seq a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Ozsolak F, Milos PM. RNA sequencing: advances, challenges and opportunities. Nat Rev Genet. 2011;12:87–98.View ArticlePubMedGoogle Scholar
- Soetaert SS, Van Neste CM, Vandewoestyne ML, Head SR, Goossens A, Van Nieuwerburgh FC, et al. Differential transcriptome analysis of glandular and filamentous trichomes in Artemisia annua. BMC Plant Biol. 2013;13:220.View ArticlePubMedPubMed CentralGoogle Scholar
- Nair P, Misra A, Singh A, Shukla AK, Gupta MM, Gupta AK, et al. Differentially expressed genes during contrasting growth stages of Artemisia annua for artemisinin content. PLoS One. 2013;8:e60375.View ArticlePubMedPubMed CentralGoogle Scholar
- Senthil K, Wasnik NG, Kim YJ, Yang DC. Generation and analysis of expressed sequence tags from leaf and root of Withania somnifera (Ashwgandha). Mol Biol Rep. 2010;37:893–902.View ArticlePubMedGoogle Scholar
- Gupta P, Goel R, Pathak S, Srivastava A, Singh SP, Sangwan RS, et al. De novo assembly, functional annotation and comparative analysis of Withania somnifera leaf and root transcriptomes to identify putative genes involved in the withanolides biosynthesis. PLoS One. 2013;8:e62714.View ArticlePubMedPubMed CentralGoogle Scholar
- Gagne SJ, Stout JM, Liu E, Boubakir Z, Clark SM, Page JE. Identification of olivetolic acid cyclase from Cannabis sativa reveals a unique catalytic route to plant polyketides. Proc Natl Acad Sci. 2012;109:12811–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Li C, Zhu Y, Guo X, Sun C, Luo H, Song J, et al. Transcriptome analysis reveals ginsenosides biosynthetic genes, microRNAs and simple sequence repeats in Panax ginseng. BMC Genomics. 2013;14:245.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramilowski JA, Sawai S, Seki H, Mochida K, Yoshida T, Sakurai T, et al. Glycyrrhiza uralensis transcriptome landscape and study of phytochemicals. Plant Cell Physiol. 2013;54:697–710.View ArticlePubMedGoogle Scholar
- Gahlan P, Singh HR, Shankar R, Sharma N, Kumari A, Chawla V, et al. De novo sequencing and characterization of Picrorhiza kurrooa transcriptome at two temperatures showed major transcriptome adjustments. BMC Genomics. 2012;13:126.View ArticlePubMedPubMed CentralGoogle Scholar
- Pandey A, Swarnkar V, Pandey T, Srivastava P, Kanojiya, Mishra DK, et al. Transcriptome and metabolite analysis reveal candidate genes of the cardiac glycoside biosynthetic pathway from Calotropis procera. Sci Rep. 2016;6:34464.View ArticlePubMedPubMed CentralGoogle Scholar
- Akhtar MQ, Qamar N, Yadav P, Kulkarni P, Kumar A, Shasany AK. Comparative glandular trichome transcriptome based gene characterization reveals reasons for differential (−)-menthol biosynthesis in Mentha species. Physiol Plant. 2017;160(2):128–41.View ArticlePubMedGoogle Scholar
- Garg A, Agrawal L, Mishra RC, Sharma S, Gosh S. Andrographis paniculata Transcriptome provides molecular insights into tissue-specific accumulation of medicinal diterpenes. BMC Genomics. 2015;16:659.View ArticlePubMedPubMed CentralGoogle Scholar
- Upadhyay S, Phukan UJ, Mishra S, Shukla RK. De novo leaf and root transcriptome analysis identified novel genes involved in steroidal sapogenin biosynthesis in Asparagus racemosus. BMC Genomics. 2014;15:746.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim YJ, Zhang D, Yang DC. Biosynthesis and biotechnological production of ginsenosides. Biotechnology advances. 2015;33(6):717–35.Google Scholar
- Liu MH, Yang BR, Cheung WF, Yang KY, Zhou HF, Kwok JS, et al. Transcriptome analysis of leaves, roots and flowers of Panax notoginseng identifies genes involved in ginsenoside and alkaloid biosynthesis. BMC Genomics. 2015;16:265.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Q, Ma C, Qian J, Lan X, Chao N, Sun J, et al. Transcriptome sequencing of Gynostemma pentaphyllum to identify genes and enzymes involved in Triterpenoid biosynthesis. Int J Genomics. 2016;2016:7840914.View ArticlePubMedPubMed CentralGoogle Scholar
- Gan RC, Chen TW, Wu TH, Huang PJ, Lee CC, Yeh YM, et al. PARRoT- a homology-based strategy to quantify and compare RNA-sequencing from non-model organisms. BMC Bioinformatics. 2016;17:513.View ArticlePubMedPubMed CentralGoogle Scholar
- Kumar S, Kalra S, Kumar S, Kaur J, Singh K. Differentially expressed transcripts from leaf and root tissue of Chlorophytum borivilianum: a plant with high medicinal value. Gene. 2012;511:79–87.View ArticlePubMedGoogle Scholar
- Kalra S, Kumar S, Lakhanpal N, Kaur J, Singh K. Characterization of Squalene synthase Gene from Chlorophytum borivilianum (Sant. And Fernand.). Mol Biotechnol. 2013;54:944–53.View ArticlePubMedGoogle Scholar
- Van der Fits L, Memelink J. ORCA3, a jasmonate responsive transcriptional regulator of plant primary and secondary metabolism. Science. 2000;289:295–7.View ArticlePubMedGoogle Scholar
- Pauwels L, Morreel K, De Witte E, Lammertyn F, Van Montagu M, Boerjan W, et al. Mapping methyl jasmonate mediated transcriptional reprogramming of metabolism and cell cycle progression in cultured Arabidopsis cells. Proc Natl Acad Sci. 2008;105:1380–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Sui C, Zhang J, Wei J, Chen S, Li Y, Xu J, et al. Transcriptome analysis of Bupleurum chinense focusing on genes involved in the biosynthesis of saikosaponins. BMC Genomics. 2011;12:539.View ArticlePubMedPubMed CentralGoogle Scholar
- Luo H, Sun C, Sun Y, Wu Q, Li Y, Song J, et al. Analysis of the transcriptome of Panax notoginseng root uncovers putative triterpene saponin-biosynthetic genes and genetic markers. BMC Genomics. 2011;12:S5.View ArticlePubMedPubMed CentralGoogle Scholar
- Garai S, Mahato SB, Ohtani K, Yamasaki K. Dammarane-type triterpenoid saponins from Bacopa monniera. Phytochemistry. 1996;42:815–20.View ArticlePubMedGoogle Scholar
- Cheong YH, Chang HS, Gupta R, Wang H, Zhu T. Transcriptional profiling reveals novel interactions between wounding, pathogen, abiotic stress, and hormonal responses in Arabidopsis. Plant Physiol. 2002;129:661–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsaballa A, Nikolaidis A, Trikka F, Ignea C, Kampranis SC, Makris AM, et al. Use of the de novo transcriptome analysis of silver leaf nightshade (Solanum elaeagnifolium) to identify gene expression changes associated with wounding and terpene biosynthesis. BMC Genomics. 2015;16:504.View ArticlePubMedPubMed CentralGoogle Scholar
- Namdeo AG. Plant cell elicitation for production of secondary metabolites: a review. Pharmacogn Rev. 2007;1:69–79.Google Scholar
- Suzuki H, Reddy MS, Naoumkina M, Aziz N, May GD, Huhman DV, et al. Methyl jasmonate and yeast elicitor induce differential transcriptional and metabolic reprogramming in cell suspension cultures of the model legume Medicago truncatula. Planta. 2005;220:696–707.View ArticlePubMedGoogle Scholar
- Rahnamaie-Tajadod R, Loke KK, Goh HH, Noor NM. Differential Gene expression analysis in Polygonum minus leaf upon 24 h of methyl Jasmonate elicitation. Front Plant Sci. 2017;8:109.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Z, Ji L, Wang J, Jin J, Yang X, Rao P, et al. Dynamic changes in the transcriptome of Populus hopeiensis in response to abscisic acid. Sci Rep. 2017;7:42708.View ArticlePubMedPubMed CentralGoogle Scholar
- Peng X, Wang Y, He R, Zhao M, Shen S. Global transcriptomics identification and analysis of transcriptional factors in different tissues of the paper mulberry. BMC Plant Biol. 2014;14:194.View ArticlePubMed CentralGoogle Scholar