Transcriptome Sequencing of Gynostemma pentaphyllum to Identify Genes and Enzymes Involved in Triterpenoid Biosynthesis

Genome-wide characterization of the bHLH gene family in Gynostemma pentaphyllum reveals its potential role in the regulation of gypenoside biosynthesis

BACKGROUND

Gynostemma pentaphyllum, an ancient Chinese herbal medicine, serves as a natural source of gypenosides with significant medicinal properties. Basic helix-loop-helix (bHLH) transcription factors play pivotal roles in numerous biological processes, especially in the regulation of secondary metabolism in plants. However, the characteristics and functions of the bHLH genes in G. pentaphyllum remain unexplored, and their regulatory role in gypenoside biosynthesis remains poorly elucidated.

RESULTS

This study identified a total of 111 bHLH members in G. pentaphyllum (GpbHLHs), categorizing them into 26 subgroups based on shared conserved motif compositions and gene structures. Collinearity analysis illustrated that segmental duplications predominately lead to the evolution of GpbHLHs, with most duplicated GpbHLH gene pairs undergoing purifying selection. Among the nine gypenoside-related GpbHLH genes, two GpbHLHs (GpbHLH15 and GpbHLH58) were selected for further investigation based on co-expression analysis and functional prediction. The expression of these two selected GpbHLHs was dramatically induced by methyl jasmonate, and their nuclear localization was confirmed. Furthermore, yeast one-hybrid and dual-luciferase assays demonstrated that GpbHLH15 and GpbHLH58 could bind to the promoters of the gypenoside biosynthesis pathway genes, such as GpFPS1, GpSS1, and GpOSC1, and activate their promoter activity to varying degrees.

CONCLUSIONS

In conclusion, our findings provide a detailed analysis of the bHLH family and valuable insights into the potential use of GpbHLHs to enhance the accumulation of gypenosides in G. pentaphyllum.

Diploid chromosome-level reference genome and population genomic analyses provide insights into Gypenoside biosynthesis and demographic evolution of Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum (Thunb.) Makino is a perennial creeping herbaceous plant in the family Cucurbitaceae, which has great medicinal value and commercial potential, but urgent conservation efforts are needed due to the gradual decreases and fragmented distribution of its wild populations. Here, we report the high-quality diploid chromosome-level genome of G. pentaphyllum obtained using a combination of next-generation sequencing short reads, Nanopore long reads, and Hi-C sequencing technologies. The genome is anchored to 11 pseudo-chromosomes with a total size of 608.95 Mb and 26 588 predicted genes. Comparative genomic analyses indicate that G. pentaphyllum is estimated to have diverged from Momordica charantia 60.7 million years ago, with no recent whole-genome duplication event. Genomic population analyses based on genotyping-by-sequencing and ecological niche analyses indicated low genetic diversity but a strong population structure within the species, which could classify 32 G. pentaphyllum populations into three geographical groups shaped jointly by geographic and climate factors. Furthermore, comparative transcriptome analyses showed that the genes encoding enzyme involved in gypenoside biosynthesis had higher expression levels in the leaves and tendrils. Overall, the findings obtained in this study provide an effective molecular basis for further studies of demographic genetics, ecological adaption, and systematic evolution in Cucurbitaceae species, as well as contributing to molecular breeding, and the biosynthesis and biotransformation of gypenoside.

Integrated Metabolomic and Transcriptomic Analysis and Identification of Dammarenediol-II Synthase Involved in Saponin Biosynthesis in Gynostemma longipes

Gynostemma longipes contains an abundance of dammarane-type ginsenosides and gypenosides that exhibit extensive pharmacological activities. Increasing attention has been paid to the elucidation of cytochrome P450 monooxygenases (CYPs) and UDP-dependent glycosyltransferases (UGTs) that participate downstream of ginsenoside biosynthesis in the Panax genus. However, information on oxidosqualene cyclases (OSCs), the upstream genes responsible for the biosynthesis of different skeletons of ginsenoside and gypenosides, is rarely reported. Here, an integrative study of the metabolome and the transcriptome in the leaf, stolon, and rattan was conducted and the function of GlOSC1 was demonstrated. In total, 46 triterpenes were detected and found to be highly abundant in the stolon, whereas gene expression analysis indicated that the upstream OSC genes responsible for saponin skeleton biosynthesis were highly expressed in the leaf. These findings indicated that the saponin skeletons were mainly biosynthesized in the leaf by OSCs, and subsequently transferred to the stolon via CYPs and UGTs biosynthesis to form various ginsenoside and gypenosides. Additionally, a new dammarane-II synthase (DDS), GlOSC1, was identified by bioinformatics analysis, yeast expression assay, and enzyme assays. The results of the liquid chromatography-mass spectrometry (LC-MS) analysis proved that GlOSC1 could catalyze 2,3-oxidosqualene to form dammarenediol-II via cyclization. This work uncovered the biosynthetic mechanism of dammarenediol-II, an important starting substrate for ginsenoside and gypenosides biosynthesis, and may achieve the increased yield of valuable ginsenosides and gypenosides produced under excess substrate in a yeast cell factory through synthetic biology strategy.

Gene excavation and expression analysis of CYP and UGT related to the post modifying stage of gypenoside biosynthesis in Gynostemma pentaphyllum (Thunb.) Makino by comprehensive analysis of RNA and proteome sequencing

Previous studies have revealed that gypenosides produced from Gynostemma pentaphyllum (Thunb.) Makino are mainly dammarane-type triterpenoid saponins with diverse structures and important biological activities, but the mechanism of diversity for gypenoside biosynthesis is still unclear. In this study, a combination of isobaric tags for relative and absolute quantification (iTRAQ) proteome analysis and RNA sequencing transcriptome analysis was performed to identify the proteins and genes related to gypenoside biosynthesis. A total of 3925 proteins were identified by proteomic sequencing, of which 2537 were quantified. Seventeen cytochrome P450 (CYP) and 11 uridine 5’-diphospho-glucuronosyltransferase (UDP-glucuronosyltransferase, UGT) candidate genes involved in the side chain synthesis and modification of gypenosides were found. Seven putative CYPs (CYP71B19, CYP77A3, CYP86A7, CYP86A8, CYP89A2, CYP90A1, CYP94A1) and five putative UGTs (UGT73B4, UGT76B1, UGT74F2, UGT91C1 and UGT91A1) were selected as candidate structural modifiers of triterpenoid saponins, which were cloned for gene expression analysis. Comprehensive analysis of RNA sequencing and proteome sequencing showed that some CYPs and UGTs were found at both the transcription and translation levels. In this study, an expression analysis of 7 CYPs and 5 UGTs that contributed to gypenoside biosynthesis and distribution in G. pentaphyllum was performed, providing consistent results that will inspire more future research on vital genes/proteins involved in gypenoside biosynthesis.

Identification of three groups of ginsenoside biosynthetic UDP-glycosyltransferases from Gynostemma pentaphyllum

Ginsenosides are glycosylated dammarene-type triterpenes that have been identified in distantly related Panax ginseng and Gynostemma pentaphyllum. The phylogenetic relatedness of the ginsenoside biosynthetic genes in the two species was previously unknown. The final steps of ginsenoside biosynthesis are the glycosylations of hydroxylated triterpenes, protopanaxadiol (PPD) and protopanaxatriol (PPT), and their glycosylated forms by UDP-glycosyltransferases (UGTs). Ginsenoside biosynthetic UGTs have been identified in Panax but not in Gynostemma. Through a biochemical screening of Gynostemma UGTs (GpUGTs), we herein identified three groups of ginsenoside biosynthetic GpUGTs. These groups comprise: two GpUGTs that belong to the UGT71 family and glucosylate the C20-OH positions of PPD- and PPT-type ginsenosides; one GpUGT that belongs to the UGT74 family and glucosylates the C3-OH position of PPD-type ginsenosides; and two GpUGTs that belong to the UGT94 family and add a glucose to the C3-O-glucosides of PPD-type ginsenosides. These GpUGTs belong to the same UGT families as the ginsenoside biosynthetic Panax UGTs (PgUGTs). However, GpUGTs and PgUGTs belong to different subfamilies. Furthermore, cucumber UGTs orthologous to GpUGTs do not glucosylate ginsenosides. These results collectively suggest that, during evolution, P. ginseng and G. pentaphyllum independently opted to use the same UGT families to synthesize ginsenosides.

NAA and 6-BA promote accumulation of oleanolic acid by JA regulation in Achyranthes bidentata Bl

Application of plant growth regulators has become one of the most important means of improving yield and quality of medicinal plants. To understand the molecular basis of phytohormone-regulated oleanolic acid metabolism, RNA-seq was used to analyze global gene expression in Achyranthes bidentata treated with 2.0 mg/L 1-naphthaleneacetic acid (NAA) and 1.0 mg/L 6-benzyladenine (6-BA). Compared with untreated controls, the expression levels of 20,896 genes were significantly altered with phytohormone treatment. We found that 13071 (62.5%) unigenes were up-regulated, and a lot of differentially expressed genes involved in hormone or terpenoid biosynthesis, or transcription factors were significantly up-regulated. These results suggest that oleanolic acid biosynthesis induced by NAA and 6-BA occurs due to the expression of key genes involved in jasmonic acid signal transduction. This study is the first to analyze the production and hormonal regulation of medicinal A. bidentata metabolites at the molecular level. The results herein contribute to a better understanding of the regulation of oleanane-type triterpenoid saponins accumulation and define strategies to improve the yield of these useful metabolites.

Positive Selection of Squalene Synthase in Cucurbitaceae Plants

Triterpenoid saponins are secondary metabolites synthesized through isoprenoid pathways in plants. Cucurbitaceae represent an important plant family in which many species contain cucurbitacins as secondary metabolites synthesized through isoprenoid and triterpenoid pathways. Squalene synthase (SQS) is required for the biosynthesis of isoprenoids, but the forces driving the evolution of SQS remain undetermined. In this study, 10 SQS cDNA sequences cloned from 10 species of Cucurbitaceae and 49 sequences of SQS downloaded from GenBank and UniProt databases were analyzed in a phylogenetic framework to identify the evolutionary forces for functional divergence. Through phylogenetic construction and positive selection analysis, we found that SQS sequences are under positive selection. The sites of positive selection map to functional and transmembrane domains. 180L, 189S, 194S, 196S, 265I, 289P, 389P, 390T, 407S, 408A, 410R, and 414N were identified as sites of positive selection that are important during terpenoid synthesis and map to transmembrane domains. 196S and 407S are phosphorylated and influence SQS catalysis and triterpenoid accumulation. These results reveal that positive selection is an important evolutionary force for SQS in plants. This provides new information into the molecular evolution of SQS within the Cucurbitaceae family.

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.