- Research article
- Open Access
Comparative transcriptome analysis of root, stem, and leaf tissues of Entada phaseoloides reveals potential genes involved in triterpenoid saponin biosynthesis
BMC Genomics volume 21, Article number: 639 (2020)
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.
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.
The dataset reported here will facilitate the research about the functional genomics of triterpenoid saponin biosynthesis and genetic engineering of Entada phaseoloides.
Entada phaseoloides (L.) Merr. is a liana belonging to Fabaceae family. It grows in Southern China and other tropical countries. The stem of Entada phaseoloides is popularly used in traditional medicine because of its significant pharmacological activities [1,2,3]. The stem of Entada phaseoloides, also called “Guo Gang Long,” produces curative effects that dispel wind and dampness and exhibits remarkable anti-inflammatory activity. Its main bioactive ingredients are triterpenoid saponins compounds . Various types of triterpene saponins have been isolated from E. phaseoloides. The representative saponins of Entada phaseoloides are oleanane-type triterpene saponins which contain seven sugar chains.
The mevalonic acid (MVA) pathway is an important metabolic pathway in plants [4, 5]. Triterpenoid saponins comprise six isoprene units and are derived from a C-30 hydrocarbon precursor, squalene. Squalene is synthesized from isopentenyl diphosphate (IPP) via the MVA pathway. Subsequently, squalene epoxidase (SQE) catalyzes the conversion of squalene to 2,3-oxidosqualene. The diversifying step in triterpenoid backbone biosynthesis is the cyclization of 2,3-oxidosqualene catalyzed by a class of oxidosqualene cyclases (OSCs) [6, 7]. Cytochrome P450 monooxygenases (CYP450s) and UDP-glycosyltransferases (UGTs) govern the hydroxylation, oxidation, and glycosylation steps, yielding triterpenoid saponins [8,9,10]. However, the key genes related to triterpenoid saponin biosynthesis in Entada phaseoloides have not been identified.
High-throughput sequencing analysis is a useful method to clarify the molecular mechanism of plant secondary metabolism [11, 12]. Recently, transcriptome assay with next-generation sequencing has been extensively used to explore the novel genes underlying active-ingredient biosynthesis pathways in medicinal plants. Some include the excavation of genes encoding enzymes that catalyze distinct steps related to the biosynthetic pathway of ginsenosides in Panax ginseng , triterpenoid saponin biosynthesis in Bacopa monnieri , artemisinin in Artemisia annua [15, 16], flavonoid biosynthesis in safflower , glycyrrhizin in Glycyrrhiza uralensis , rubber in Parthenium argentatum , cardiac glycoside in Calotropis procera , terpenoid in Cinnamomum camphora , cannabinoids in Cannabis sativa , withanolide in Withania somnifera , picrosides in Picrorhiza kurrooa , paclitaxel in Taxus chinensis  and steroidal saponins in Asparagus racemosus .
Because the synthesis and accumulation of specific metabolites in different tissues depends on the age of the plant, and are greatly affected by the different developmental stages. It was found that the content of triterpenes accumulated in the leaves of P. ginseng are higher in the early growth stage, while the content of triterpenes in the roots of old plants are higher . Comparative transcriptome analysis including root and leaf tissues to excavate transcripts related to saponin biosynthesis is already reported for many plants, such as Hedera helix , Panax notoginseng  and Asparagus racemosus . The maximum triterpenoid saponin content was identified in stem. So, in this study, comparative transcriptome analysis of root, stem, and leaf tissues of Entada phaseoloides was performed to identify genes related to triterpenoid saponin. We obtained thousands of putative genes, including a series of genes related to triterpene saponin biosynthesis. Moreover, different expression patterns of CYP450s and UGTs in the three tissues were analyzed. This work was established to functionally research the genes related to triterpene saponin biosynthesis and provide more information about this species.
Illumina sequencing and de novo assembly
To characterize the transcriptomes of Entada phaseoloides, we sequenced nine cDNA libraries prepared from the root, stem, and leaf tissues with three biological repeats by using the Illumina Hiseq 2500 platform. A total of 53.26 Gb clean data were obtained after removing adaptors, poly-A tails, and primer sequences, short (< 50 bp), and low-quality sequences. A total of 57–60 million for each tissue were generated (Additional file 1). The high-quality reads were assembled using the Trinity program  and the TGI clustering tool (TGICL)  to remove redundant sequences. Finally, 116,910 unigenes were identified, with an average N50 length of 1218 bp (Additional file 1). The correlation indices between repeated samples were > 0.9 (Additional file 2), indicating that the Illumina sequencing results are credible.
All assembled unigenes were searched against the Non-redundant (Nr), Uniprot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Pfam, Gene Ontology (GO), and Clusters of Orthologous Groups (COG) databases using the BLASTx program with E-value <1e-5. Among the 116,910 sequences, 42,191 (36.1%), 41,228 (35.3%), 28,126 (24.1%), 26,874 (23.0%), 15,119 (13.0%) and 11,812 (10.1%) unigenes showed significant similarity to known proteins in NR, Uniprot, GO, Pfam, KEGG and COG database, respectively. The result of BLASTX with different databases and their annotation were listed in Additional file 3. Based on the Nr database, the E-value distribution indicated that 70.67% of the matched unigenes ranged from 1e-5 to 1e-100 (Fig. 1a). For the similarity distribution, 46.43% unigenes exhibited a similarity of above 80%, whereas 50.79% of the unigenes showed a similarity of 40–80% (Fig. 1b). Furthermore, 14.63% of the Entada phaseoloides unigenes shared high similarity with the genes of Cicer arietinum, 13.46% similarity with Cajanus cajan, 12.84% similarity with Glycine max, and 9.07% similarity with Medicago truncatula (Fig. 1c).
GO analysis included three main domains that describe biological processes, cellular components, and molecular functions. When GO was used to classify gene functions, 28,126 unigenes were assigned to 60 functional categories (Additional files 4 and 5). Within the biological process domain, the three most enriched categories were “biosynthetic process,” “cellular nitrogen compound metabolic process,” and “response to stress.” In the cellular component domain, the three most matched categories were “cellular component.” “nucleus,” and “protein complex.” In the molecular function domain, the three most common categories were “ion binding,” “molecular function,” and “kinase activity.”
To better understand the functions of specific metabolic pathways in Entada phaseoloides, we mapped the annotated unigenes to the reference biological pathways in the KEGG database. A total of 15,119 unigenes (13.0%) could be assigned to five main categories and 31 sub-categories (Fig. 2, Additional file 6). These enzymes feature assigned functions in 28 secondary metabolic pathways in KEGG (Table 1). Among these unigenes, 147 encode key enzymes are related to the pathways for terpenoid biosynthesis, including the synthesis of the terpenoid backbone (63 unigenes), monoterpenoids (7 unigenes), diterpenoids (18 unigenes), sesquiterpenoids and triterpenoids (16 unigenes), and other terpenoid-quinone complexes (43 unigenes). Fifty unigenes are involved in alkaloid biosynthesis, including isoquinoline alkaloid (24 unigenes) and tropane, piperidine, and pyridine alkaloid biosynthesis (26 unigenes). Exactly 210 unigenes were associated with the flavonoid biosynthesis pathway, including the phenylpropanoid (162 unigenes), flavonoid (36 unigenes), flavone and flavonol (7 unigenes), and isoflavonoid (5 unigenes) biosynthesis pathways. Unigenes involved in these pathways should be further identified to understand their functions in the biosynthesis of active ingredients in leguminous plants.
Differentially expressed gene (DEG) analysis
The clean reads were mapped back onto the assembled unigenes by using the alignment via Burrows–Wheeler aligner (BWA) program to analyze the DEGs among different tissues . The Fragments per Kilobase Million (FPKM) value was calculated for each unigene in each tissue of Entada phaseoloides. The DEGs were identified (Additional file 7) using FDR ≤ 0.001 and |log2Ratio| ≥ 1 . The lowest number of DEGs was observed between the stem and leaf tissues, and the highest was noted between the root and leaf. Furthermore, the DEGs in one tissue were studied and compared with those in the other two tissues. The stem contained the largest number of highly expressed unigenes, having 8962 unigenes more abundant in the stem. Figure 3 shows the other expression differences among various tissues.
KEGG enrichment analyses were performed with the DEGs among the root, stem, and leaf tissues to investigate the genes regulating the distribution of triterpenoid saponin. These DEGs were evidently enriched in specific pathways. Meanwhile, the top 20 significant pathways were analyzed based on the FDR ≤ 0.01. Between the leaf and stem, plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, terpenoid backbone biosynthesis, and phosphatidyllnositol signaling system showed significant enrichment (Fig. 4a). Between the leaf and root, plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, ubiquinone and other terpenoid-quinone biosynthesis, and cyanoamino acid metabolism pathways showed visible differential expression (Fig. 4b). Between the stem and root, plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, cyanoamino acid metabolism, and terpenoid backbone biosynthesis showed significant enrichment (Fig. 4c). In addition to the common pathways of primary metabolism, enriched secondary metabolic pathways, including terpenoid and phenylpropanoid biosynthesis, were also found between different tissues, indicating the possible distinct distribution of secondary metabolites in different tissues.
Putative genes involved in triterpenoid saponin backbone biosynthesis
Triterpenes are synthetized from a five-carbon isoprene unit through the cytosolic MVA pathway. Triterpenoid saponins are composed of six isoprene units and are derived from the C-30 hydrocarbon precursor, squalene. Squalene is synthesized from isopentenyl diphosphate (IPP) via the MVA pathway. All genes encoding the enzymes associated with the upstream regions of triterpenoid biosynthesis were successfully detected in the Entada phaseoloides transcriptome. Their expression value was monitored in three biological replicates along with their mean values (Table 2, Fig. 5). Most unigenes related to MVA pathway were specifically upregulated in the stem tissue. Hydroxymethylglutaryl-CoA reductase showed the highest expression, which is the rate limiting step MVA pathway for saponin biosynthesis.
The diversifying step in triterpenoid backbone biosynthesis is the cyclization of 2,3-oxidosqualene catalyzed by a class of OSCs. The major saponins in Entada phaseoloides are oleanane-type triterpenoid saponins derived from β-amyrin. The Illumina sequencing of Entada phaseoloides revealed 21 OSC sequences, among which eight unigenes were putative β-amyrin synthases. A full-length OSC sequence (EpBAS) with high identity to β-amyrin synthase was obtained (Additional file 8). The EpBAS cDNA included a 2289 bp full open reading frame fragment. The deduced amino acid sequence of EpBAS (762 amino acids) shared 89.37 and 89.34% similarity with β-amyrin synthase in Abrus precatorius (ApBAS) and GiBAS in Glycyrrhiza inflata (Fig. 6), respectively. The relatively high similarities of the EpBAS protein with other β-amyrin synthases suggest that this gene encodes β-amyrin synthase in Entada phaseoloides.
CYP450s and UGTs
Earlier studies suggested that CYP450s and UGTs may account for the biosynthesis and accumulation of triterpene saponins in specific organs . Tissue-specific transcriptome analysis of Entada phaseoloides suggests that the enzymes involved in triterpenoid saponin backbone are present in all the three tissues. Based on DEG analysis using transcriptome data, there is the possibility of further modifications such as oxidation and glycosylation using CYP450s and UGTs occur in the stem. Although the enzymes related to precursor biosynthesis are also present in root and leaf tissues which suggest the involvement of all the three tissues in the metabolic pathway. However, the metabolic analysis has indicated that it is the stem which mostly contains higher triterpenoid saponin content and utilized widely for its excellent pharmacological activity. In this study, in total of 326 CYP450s and 148 UGTs were found. Among the DEGs, 26 CYP450s and 17 UGTs were upregulated in the stem compared with the root and leaf tissues (Fig. 7a and b).
qRT-PCR validation of candidate genes involved in triterpenoid saponin biosynthesis
To verify the expression profiles obtained from Illumina sequencing, we performed qRT-PCR on nine selected genes related to triterpene saponin biosynthesis (Fig. 8). Consistent with the Illumina data, most of these genes showed strong expression levels in the stem compared with the root and leaf, and acetyl-CoA acetyltransferase, hydroxymethylglutaryl-CoA synthase, hydroxymethylglutaryl-CoA reductase and SQE genes were expressed abundantly. The expression fold changes were also close to the RNA-seq results. qRT-PCR results indicate that the RNA-seq data in this studty were reliable.
Entada phaseoloides is an important traditional medicinal plant with various pharmaceutical activities. Although this plant is pharmacologically important, its genomic or transcriptomic information is highly limited. In NCBI, only 38 protein sequences are accessible for Entada phaseoloides. We revealed the comparative transcriptome analysis of the root, stem, and leaf tissues of Entada phaseoloides. The dataset reported here is useful in understanding the biosynthetic pathway of pharmacodynamic triterpenoid saponin and genetic engineering of this species.
In this study, a total of 53.26 Gb clean data were generated from nine RNA-seq libraries of the root, stem, and leaf. De novo assembly acquired 116,910 unigenes, with an average N50 length of 1218 bp, which is similar to that of previously reported non-model plants, such as Raphanus sativus  and Isodon Amethystoides . The best match for each unigene search against the Nr and KEGG databases was of help to assign GO functional annotation under biological process, cellular component, and molecular function categories. The varied GO assignments to unigenes represented the possible assortment of genes in the Entada phaseoloides transcriptome. Several unigenes mapped onto KEGG are related to distinct secondary metabolic pathways. Most unmatched unigenes are short sequence proteins with no domain, untranslated regions, non-coding RNA or assembly mistakes. In support of the annotation, all the unigenes encoding enzymes related to the upstream regions of the MVA pathway for saponin biosynthesis from acetyl-CoA to squalene were found.
SQE enzymes catalyze the oxidation of squalene to 2,3-oxidosqualene. In our transcriptomic analysis, sequences encoding SQE represented the highest number (16) of unigenes associated with the MVA pathway. Single copies of SQE were identified in mouse and yeast, and the destruction of SQE in these species is lethal . However, two or more copies of SQE are usually found in plants. Hwang et al.  examined 17 SQE sequences in Eleutherococcus senticosus. In Arabidopsis thaliana, six SQE enzymes have been identified, and three of them encode functional SQEs . The expression of PgSQE1 regulates the biosynthesis of ginsenoside in Panax ginseng . Thus, SQE is possibly an important enzyme in the saponin biosynthetic pathway. The SQE enzyme responsible for the saponin biosynthesis in the 16 SQE sequences in Entada phaseoloides remains to be identified.
The cyclization of 2,3-oxidosqualene is a branch point of saponin synthesis. The major saponins in the stem of Entada phaseoloides are oleanane-type triterpenoids. Oleanolic acid sapogenin was derived from β-amyrin after hydroxylation by CYP450s and glycosylation by UGTs [41, 42]. A total of 8 β-amyrin synthase, 326 CYP450, and 148 UGT sequences were observed in our transcriptome. The high expression of β-amyrin synthase, an important enzyme related to triterpenoid sapogenin biosynthesis at later stages, further reveals the high concentration of sapogenins in the stem of Entada phaseoloides. Similar tissue-specific concentrations of triterpenoid sapogenins have already been reported in other plants [43,44,45]. Moreover, 26 CYP450s and 17 UGTs were found to be upregulated in the stem. Further characterization of these candidate enzymes is needed to confirm the pathway of triterpenoid saponin biosynthesis in Entada phaseoloides.
In the present study, the comparative transcriptome analysis of root, stem and leaf tissues of Entada phaseoloides was performed to investigate the putative genes involved in triterpenoid saponin biosynthetic pathway of an important medicinal plant. The differential expression pattern of pathway genes suggest tissue-specific synthesis. The identified data will help the further discovery and functional genomics and transcriptomics analysis of Entada phaseoloides.
Three-year-old healthy wild-type Entada phaseoloides plants were collected from the experimental farm of South China Botanical Garden, Guangzhou City, Guangdong Province, P.R. China, in May 2019. After cleaning with ultrapure water, the roots, stems, and leaves were collected separately, immediately frozen in liquid nitrogen, and stored at − 80 °C.
RNA extraction, cDNA synthesis, and sequencing
Total RNA from approximately 1.0 g of each tissue was extracted using TRIzol (Invitrogen, Canada) following the manufacturer’s instructions. Three replicates were employed for each experiment. mRNA was isolated from total RNA by using Oligo (dT) magnetic beads. By mixing with fragmentation buffer, the mRNA was broken into short fragments.
The short fragments were purified and resolved with EB buffer. After end repair and single base “A”addition, adapters were ligated to the cDNA molecules. To select suitable cDNA fragments for PCR amplification, we purified the sample library with the AMPure XP system (Beckman Coulter, USA). Finally, PCR products were purified, and library quality was were determined using an Agilent Bioanalyzer 2100 system (Agilent Technologies, USA) and a Qubit 3.0 fluorometer (Invitrogen, USA). Each cDNA library was sequenced in a single lane of the Illumina Hiseq 2500 platform.
Data filtering and de novo assembly
The raw reads were first filtered to exclude the reads containing adaptors or with ambiguous nucleotides (‘N’). Next, the low-quality reads having more than 20% Q < 20 bases were also trimmed. The yielded high-quality clean reads were used to the develop sequence assembly by using the Trinity software. After the removal of redundant Trinity-generated sequences by using the TGICL, clusters and unigenes were finally obtained.
Functional annotation and classification
All assembled unigenes were annotated by BLASTx analysis against the Nr (http://www.ncbi.nlm.nih.gov/), UniProt (http://www.uniprot.org/downloads), Pfam (http://pfam.xfam.org/), COG (http://www.ncbi.nlm.nih.gov/COG/) databases with an E-value <1e-5. Only the top hit results were extracted for each unigene. GO (http://www.geneontology.org) terms were functionally classified based on Nr annotations by using the Blast2go program (http://www.Blast2 go.de/). KEGG (http://www.genome.jp/kegg/) was used to draw metabolic maps. The KEGG analysis results included KEGG orthology (KO) numbers and enzyme commission (EC) numbers.
Clean reads were mapped back onto the assembled unigenes by using the BWA program. The FPKM value was calculated for each unigene in each tissue of Entada phaseoloides. The expression difference was analyzed by Fisher’s exact test, and the FDR for each gene were obtained. DEGs were required to have thresholds of FDR ≤ 0.001 and |log2Ratio| ≥ 1. KEGG pathways were then reconstructed on DEGs.
Verification of gene expression by using qRT-PCR
Nine genes related to triterpene saponin biosynthesis were selected for validation by qRT-PCR. The primers for qRT-PCR analysis are listed in Additional file 9. All reactions were performed on the CFX96 real-time PCR system (Bio-Rad, USA) with a SYBR® Premix Ex Taq™ kit (Takara, China). The Actin gene was used as an internal control (Additional file 10). Each qRT-PCR experiment was performed with three biological repeats. The relative gene expression was calculated using the 2−ΔΔCT method.
Availability of data and materials
Final sequences obtained from all the three tissues were submitted to the SRA database of NCBI with accession number PRJNA597694.
Cytochrome P450 monooxygenase
The TGI clustering tool
Kyoto Encyclopedia of Genes and Genomes
Clusters of Orthologous Groups
Differentially expressed gene
Quantitative reverse transcription-polymerase chain reaction
Xiong H, Ding X, Yang XZ, Yang GZ, Mei ZN. Triterpene Saponins from the stems of Entada phaseoloides. Planta Med. 2014;80:710–8.
Dong Y, Shi H, Yang H, Peng Y, Wang M, Li X. Antioxidant phenolic compounds from the stems of Entada phaseoloides. Chem Biodivers. 2012;9:68–79.
Xiong H, Zheng YN, Yang GZ, Wang HX, Mei ZN. Triterpene saponins with anti-inflammatory activity from the stems of Entada phaseoloides. Fitoterapia. 2015;103:33–45.
Haralampidis K, Trojanowska M, Osbourn AE. Biosynthesis of triterpenoid saponins in plants. Adv Biochem Eng Biotechnol. 2002;75:31–49.
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.
Abe I, Rohmer M, Prestwich GD. Enzymatic cyclization of squalene and oxidosqualene to sterols and triterpenes. Chem Rev. 1993;93:2189–206.
Xu R, Fazio GC, Matsuda SPT. On the origins of triterpenoid skeletal diversity. Phytochemistry. 2004;65:261–91.
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(Suppl 5):S5.
Zhao YJ, Li C. Biosynthesis of plant triterpenoid saponins in microbial cell factories. Agr Food Chem. 2018;66:12155–65.
Seki H, Tamura K, Muranaka T. P450s and UGTs: key players in the structural diversity of triterpenoid saponins. Plant Cell Physiol. 2015;56:1463–71.
Wang Z, Gerstein M, Snyder M. RNA-Seq a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63.
Belair CD, Hu T, Chu B, Freimer JW, Cooperberg MR, Blelloch RH. High-throughput, efficient, and unbiased capture of small RNAs from low-input samples for sequencing. Sci Rep. 2019;9:2262.
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.
Jeena GS, Fatima S, Tripathi P, Upadhyay S, Shukla RK. Comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri identifies potential genes related to triterpenoid saponin biosynthesis. BMC Genomics. 2017;18:1–15.
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.
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.
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.
Chen J, Tang XH, Ren CX, Wei B, Wu YY, Wu QH, et al. Full-length transcriptome sequences and the identification of putative genes for flavonoid biosynthesis in safflower. BMC Genomics. 2018;19:548.
Stonebloom SH, Scheller HV. Transcriptome analysis of rubber biosynthesis in guayule (Parthenium argentatum gray). BMC Plant Biol. 2019;19:71.
Pandey A, Swarnkar V, Pandey T, Srivastava P, Kanojiya MDK, et al. Transcriptome and metabolite analysis reveal candidate genes of the cardiac glycoside biosynthetic pathway from Calotropis procera. Sci Rep. 2016;6:34464.
Chen CH, Zheng YJ, Zhong YD, Wu YF, Li ZT, Xu LA, et al. Transcriptome analysis and identification of genes related to terpenoid biosynthesis in Cinnamomum camphora. BMC Genomics. 2018;19:550.
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.
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.
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.
Liao WF, Zhao SY, Zhang M, Dong KG, Chen Y, Fu CH, et al. Transcriptome assembly and systematic identification of novel cytochrome P450s in Taxus chinensis. Front Plant Sci. 2017;8:1468.
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.
Kim YJ, Zhang D, Yang DC. Biosynthesis and biotechnological production ofginsenosides. Biotechnol Adv. 2015;33:717–35.
Sun HP, Li F, Xu ZJ, Sun ML, Cong HQ, Qiao F, et al. De novo leaf and root transcriptome analysis to identify putative genes involved in triterpenoid saponins biosynthesis in Hedera helix L. PLoS One. 2017;12:e0182243.
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.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat Biotechnol. 2011;29:644–52.
Perteal G, Huang XQ, Liang F, Antonescu V, Sultanal R, Karamycheval S, et al. TIGR gene indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003;19:651–2.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25:1754–60.
Wang LK, Feng ZX, Wang X, Wang XW, Zhang XG. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26:136–8.
Yendo AC, de Costa F, Gosmann G, Fett-Neto AG. Production of plant bioactive triterpenoid saponins: elicitation strategies and target genes to improve yields. Mol Biotechnol. 2010;46:94–104.
Zhang L, Jia H, Yin Y, Wu G, Xia H, Wang X, et al. Transcriptome analysis of leaf tissue of Raphanus sativus by RNA sequencing. PLoS One. 2013;8:e80350.
Zhao F, Sun M, Zhang W, Jiang C, Teng J, Sheng W, et al. Comparative transcriptome analysis of roots, stems and leaves of Isodon amethystoides reveals candidate genes involved in Wangzaozins biosynthesis. BMC Plant Biol. 2018;18:272.
Landl KM, Klösch B, Turnowsky F. ERG1, encoding squalene epoxidase, is located on the right arm of chromosome VII of Saccharomyces cerevisiae. Yeast. 1996;12:609–13.
Hwang HS, Lee H, Choi YE. Transcriptomic analysis of Siberian ginseng (Eleutherococcus senticosus) to discover genes involved in saponin biosynthesis. BMC Genomics. 2015;16:180.
Rasbery JM, Shan H, LeClair RJ, Norman M, Matsuda SP, Bartel B. Arabidopsis thaliana squalene epoxidase 1 is essential for root and seed development. J Biol Chem. 2007;282:17002–13.
Han JY, In JG, Kwon YS, Choi YE. Regulation of ginsenoside and phytosterol biosynthesis by RNA interferences of squalene epoxidase gene in Panax ginseng. Phytochemistry. 2010;71:36–46.
Fukushima EO, Seki H, Ohyama K, Ono E, Umemoto N, Mizutani M, et al. CYP716A subfamily members are multifunctional oxidases in triterpenoid biosynthesis. Plant Cell Physiol. 2011;52:2050–61.
Han JY, Kim MJ, Ban YW, Hwang HS, Choi YE. The involvement of β-amyrin 28-oxidase (CYP716A52v2) in oleanane-type ginsenoside biosynthesis in Panax ginseng. Plant Cell Physiol. 2013;54:2034–46.
Perez SL, Scossa F, Proost S, Bitocchi E, Papa R, Tohge T, et al. Multi-tissue integration of transcriptomic and specialized metabolite profiling provides tools for assessing the common bean (Phaseolus vulgaris) metabolome. Plant J. 2019;97:1132–53.
Wu Q, Ma X, Zhang K, Feng X. Identification of reference genes for tissue-specific gene expression in Panax notoginseng using quantitative real-time PCR. Biotechnol Lett. 2015;37:197–204.
Alonso-Serra J, Safronov O, Lim KJ, Fraser-Miller SJ, Blokhina OB, Campilho A, et al. Tissue-specific study across the stem reveals the chemistry and transcriptome dynamics of birch bark. New Phytol. 2019;222:1816–31.
The authors acknowledge Prof. Zulin Ning (South China Botanical Garden) for providing the root, stem, and leaf tissues of Entada phaseoloides.
This work was financially supported by the Hubei Provincial Natural Science Foundation of China (Grant no. 2019CFB265), the Research and Innovation Initiatives of WHPU (Grant no. 2019Y03), and the National Natural Science Foundation of China (Grant no. 31700100). The Funding bodies were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of transcriptome sequencing and assembly results.
Correlation indices between different samples.
Functional annotation of Entada phaseoloides unigenes.
Frequencies of unigenes matching GO terms.
GO enrichment of unigenes.
Unigenes for KEGG analysis.
DEGs in different tissues.
Sequence of β-amyrin synthase (EpBAS) gene in Entada phaseoloides.
List of qRT-PCR primer sequences.
Melting curves of reference gene Actin for qRT-PCR amplification.
About this article
Cite this article
Liao, W., Mei, Z., Miao, L. et al. Comparative transcriptome analysis of root, stem, and leaf tissues of Entada phaseoloides reveals potential genes involved in triterpenoid saponin biosynthesis. BMC Genomics 21, 639 (2020). https://doi.org/10.1186/s12864-020-07056-1
- Entada phaseoloides
- Triterpenoid saponins
- Secondary metabolites