The Basic Helix-Loop-Helix Transcription Factor GubHLH3 Positively Regulates Soyasaponin Biosynthetic Genes in Glycyrrhiza uralensis

Glycyrrhizin Production in Licorice Hairy Roots Based on Metabolic Redirection of Triterpenoid Biosynthetic Pathway by Genome Editing

Glycyrrhizin, a type of the triterpenoid saponin, is a major active ingredient contained in the roots of the medicinal plant licorice (Glycyrrhiza uralensis, G. glabra and G. inflata), and is used worldwide in diverse applications, such as herbal medicines and sweeteners. The growing demand for licorice threatens wild resources and therefore a sustainable method of supplying glycyrrhizin is required. With the goal of establishing an alternative glycyrrhizin supply method not dependent on wild plants, we attempted to produce glycyrrhizin using hairy root culture. We tried to promote glycyrrhizin production by blocking competing pathways using CRISPR/Cas9-based gene editing. CYP93E3 CYP72A566 double-knockout (KO) and CYP93E3 CYP72A566 CYP716A179 LUS1 quadruple-KO variants were generated, and a substantial amount of glycyrrhizin accumulation was confirmed in both types of hairy root. Furthermore, we evaluated the potential for promoting further glycyrrhizin production by simultaneous CYP93E3 CYP72A566 double-KO and CYP88D6-overexpression. This strategy resulted in a 3-fold increase (∼1.4 mg/g) in glycyrrhizin accumulation in double-KO/CYP88D6-overexpression hairy roots, on average, compared with that of double-KO hairy roots. These findings demonstrate that the combination of blocking competing pathways and overexpression of the biosynthetic gene is important for enhancing glycyrrhizin production in G. uralensis hairy roots. Our findings provide the foundation for sustainable glycyrrhizin production using hairy root culture. Given the widespread use of genome editing technology in hairy roots, this combined with gene knockout and overexpression could be widely applied to the production of valuable substances contained in various plant roots.

Disruption of a licorice cellulose synthase-derived glycosyltransferase gene demonstrates its in planta role in soyasaponin biosynthesis

KEY MESSAGE

CRISPR-Cas9-mediated disruption of a licorice cellulose synthase-derived glycosyltransferase gene, GuCSyGT, demonstrated the in planta role of GuCSyGT as the enzyme catalyzing 3-O-glucuronosylation of triterpenoid aglycones in soyasaponin biosynthesis. Triterpenoid glycosides (saponins) are a large, structurally diverse group of specialized metabolites in plants, including the sweet saponin glycyrrhizin produced by licorice (Glycyrrhiza uralensis) and soyasaponins that occur widely in legumes, with various bioactivities. The triterpenoid saponin biosynthetic pathway involves the glycosylation of triterpenoid sapogenins (the non-sugar part of triterpenoid saponins) by glycosyltransferases (GTs), leading to diverse saponin structures. Previously, we identified a cellulose synthase-derived GT (CSyGT), as a newly discovered class of triterpenoid GT from G. uralensis. GuCSyGT expressed in yeast, which could transfer the sugar glucuronic acid to the C3 position of glycyrrhetinic acid and soyasapogenol B, which are the sapogenins of glycyrrhizin and soyasaponin I, respectively. This suggested that GuCSyGT is involved in the biosynthesis of glycyrrhizin and soyasaponin I. However, the in planta role of GuCSyGT in saponin biosynthesis remains unclear. In this study, we generated GuCSyGT-disrupted licorice hairy roots using CRISPR-Cas9-mediated genome editing and analyzed the saponin content. This revealed that soyasaponin I was completely absent in GuCSyGT-disrupted lines, demonstrating the in planta role of GuCSyGT in saponin biosynthesis.

Contemporary understanding of transcription factor regulation of terpenoid biosynthesis in plants

KEY MESSAGE

This review summarized how TFs function independently or in response to environmental factors to regulate terpenoid biosynthesis via fine-tuning the expression of rate-limiting enzymes. Terpenoids are derived from various species and sources. They are essential for interacting with the environment and defense mechanisms, such as antimicrobial, antifungal, antiviral, and antiparasitic properties. Almost all terpenoids have high medicinal value and economic performance. Recently, the control of enzyme genes on terpenoid biosynthesis has received a great deal of attention, but transcriptional factors regulatory network on terpenoid biosynthesis and accumulation has yet to get a thorough review. Transcription factors function as activators or suppressors independently or in response to environmental stimuli, fine-tuning terpenoid accumulation through regulating rate-limiting enzyme expression. This study investigates the advancements in transcription factors related to terpenoid biosynthesis and systematically summarizes previous works on the specific mechanisms of transcription factors that regulate terpenoid biosynthesis via hormone signal-transcription regulatory networks in plants. This will help us to better comprehend the regulatory network of terpenoid biosynthesis and build the groundwork for terpenoid development and effective utilization.

Transcriptomics and metabolomics reveal the primary and secondary metabolism changes in Glycyrrhiza uralensis with different forms of nitrogen utilization

The roots and rhizomes of Glycyrrhiza uralensis Fisch. represent the oldest and most frequently used herbal medicines in Eastern and Western countries. However, the quality of cultivated G. uralensis has not been adequate to meet the market demand, thereby exerting increased pressure on wild G. uralensis populations. Nitrogen, vital for plant growth, potentially influences the bioactive constituents of plants. Yet, more information is needed regarding the effect of different forms of nitrogen on G. uralensis. G. uralensis seedlings were exposed to a modified Hoagland nutrient solution (HNS), varying concentrations of nitrate (KNO3), or ammonium (NH4)2SO4. We subsequently obtained the roots of G. uralensis for physiology, transcriptomics, and metabolomics analyses. Our results indicated that medium-level ammonium nitrogen was more effective in promoting G. uralensis growth compared to nitrate nitrogen. However, low-level nitrate nitrogen distinctly accelerated the accumulation of flavonoid ingredients. Illumina sequencing of cDNA libraries prepared from four groups-treated independently with low/medium NH4 + or NO3 - identified 364, 96, 103, and 64 differentially expressed genes (DEGs) in each group. Our investigation revealed a general molecular and physiological metabolism stimulation under exclusive NH4 + or NO3 - conditions. This included nitrogen absorption and assimilation, glycolysis, Tricarboxylic acid (TCA) cycle, flavonoid, and triterpenoid metabolism. By creating and combining putative biosynthesis networks of nitrogen metabolism, flavonoids, and triterpenoids with related structural DEGs, we observed a positive correlation between the expression trend of DEGs and flavonoid accumulation. Notably, treatments with low-level NH4 + or medium-level NO3 - positively improved primary metabolism, including amino acids, TCA cycle, and glycolysis metabolism. Meanwhile, low-level NH4 + and NO3 - treatment positively regulated secondary metabolism, especially the biosynthesis of flavonoids in G. uralensis. Our study lays the foundation for a comprehensive analysis of molecular responses to varied nitrogen forms in G. uralensis, which should help understand the relationships between responsive genes and subsequent metabolic reactions. Furthermore, our results provide new insights into the fundamental mechanisms underlying the treatment of G. uralensis and other Glycyrrhiza plants with different nitrogen forms.

Recent advances in triterpenoid pathway elucidation and engineering

Triterpenoids are among the most assorted class of specialized metabolites found in all the taxa of living organisms. Triterpenoids are the leading active ingredients sourced from plant species and are utilized in pharmaceutical and cosmetic industries. The triterpenoid precursor 2,3-oxidosqualene, which is biosynthesized via the mevalonate (MVA) pathway is structurally diversified by the oxidosqualene cyclases (OSCs) and other scaffold-decorating enzymes such as cytochrome P450 monooxygenases (P450s), UDP-glycosyltransferases (UGTs) and acyltransferases (ATs). A majority of the bioactive triterpenoids are harvested from the native hosts using the traditional methods of extraction and occasionally semi-synthesized. These methods of supply are time-consuming and do not often align with sustainability goals. Recent advancements in metabolic engineering and synthetic biology have shown prospects for the green routes of triterpenoid pathway reconstruction in heterologous hosts such as Escherichia coli, Saccharomyces cerevisiae and Nicotiana benthamiana, which appear to be quite promising and might lead to the development of alternative source of triterpenoids. The present review describes the biotechnological strategies used to elucidate complex biosynthetic pathways and to understand their regulation and also discusses how the advances in triterpenoid pathway engineering might aid in the scale-up of triterpenoid production in engineered hosts.

Construction of a UDP-Arabinose Regeneration System for Efficient Arabinosylation of Pentacyclic Triterpenoids

Glycosylation is an important method of modifying natural products and is usually catalyzed by uridine 5'-diphosphate (UDP)-glycosyltransferase. UDP-β-l-arabinose (UDP-Ara) confers specific functions to natural products such as pentacyclic triterpenoids. However, UDP-arabinosyltransferase with high regioselectivity toward pentacyclic triterpenoids has rarely been reported. In addition, UDP-Ara is mainly biosynthesized from UDP-α-d-glucose (UDP-Glc) through several reaction steps, resulting in the high cost of UDP-Ara. Herein, UGT99D1 was systematically characterized for specifically transferring one moiety of arabinose to the C-3 position of typical pentacyclic triterpenoids. Subsequently, 15 enzymes from plants, mammals, and microorganisms were characterized, and a four-enzyme cascade comprising sucrose synthase, UDP-Glc dehydrogenase, UDP-α-d-glucuronic acid decarboxylase, and UDP-Glc 4-epimerase was constructed to transform sucrose into UDP-Ara with UDP recycling. This system was demonstrated to efficiently produce the arabinosylated derivative (Ara-BA) of typical pentacyclic triterpenoid betulinic acid (BA). Finally, the in vitro cytotoxicity tests indicated that Ara-BA showed much higher anticancer activities than BA. The established arabinosylation platform shows the potential to enhance the pharmacological activity of natural products.

Genome-wide identification of SrbHLH transcription factors highlights its potential role in rebaudioside A (RA) biosynthesis in Stevia rebaudiana

Stevia rebaudiana Bertoni is a valuable medicinal plant and an essential source of natural sweetener, steviol glycosides (SGs), with rebaudioside A (RA) being one of the main components of SGs. bHLH transcription factors play a crucial role in plant development and secondary metabolism. In this study, 159 SrbHLH genes were identified from the S. rebaudiana genome, and each gene was named based on its chromosome location. The SrbHLH proteins were then clustered into 18 subfamilies through phylogenetic analysis. The analysis of conserved motifs and gene structure further supported the classification of the SrbHLH family. Chromosomal location and gene duplication events of SrbHLH genes were also studied. Moreover, based on the RNA-Seq data of different tissues of S. rebaudiana, 28 SrbHLHs were co-expressed with structural genes involved in RA biosynthesis. The expression pattern of candidate SrbHLH genes were confirmed by qPCR. Finally, dual luciferase reporter assays (DLAs) and subcellular localization analysis verified SrbHLH22, SrbHLH111, SrbHLH126, SrbHLH142, and SrbHLH152 are critical regulators of RA biosynthesis. This study provides new insights into the function of SrbHLHs in regulating SGs biosynthesis and lays the foundation for future applications of SrbHLH genes in molecular breeding of S. rebaudiana.

OeBAS and CYP716C67 catalyze the biosynthesis of health-beneficial triterpenoids in olive (Olea europaea) fruits

The bioactive properties of olive (Olea europaea) fruits and olive oil are largely attributed to terpenoid compounds, including diverse triterpenoids such as oleanolic, maslinic and ursolic acids, erythrodiol, and uvaol. They have applications in the agri-food, cosmetics, and pharmaceutical industries. Some key steps involved in the biosynthesis of these compounds are still unknown. Genome mining, biochemical analysis, and trait association studies have been used to identify major gene candidates controlling triterpenoid content of olive fruits. Here, we identify and functionally characterize an oxidosqualene cyclase (OeBAS) required for the production of the major triterpene scaffold β-amyrin, the precursor of erythrodiol, oleanolic and maslinic acids, and a cytochrome P450 (CYP716C67) that mediates 2α oxidation of the oleanane- and ursane-type triterpene scaffolds to produce maslinic and corosolic acids, respectively. To confirm the enzymatic functions of the entire pathway, we have reconstituted the olive biosynthetic pathway for oleanane- and ursane-type triterpenoids in the heterologous host, Nicotiana benthamiana. Finally, we have identified genetic markers associated with oleanolic and maslinic acid fruit content on the chromosomes carrying the OeBAS and CYP716C67 genes. Our results shed light on the biosynthesis of olive triterpenoids and provide new gene targets for germplasm screening and breeding for high triterpenoid content.

Transcription Factors ZjMYB39 and ZjMYB4 Regulate Farnesyl Diphosphate Synthase- and Squalene Synthase-Mediated Triterpenoid Biosynthesis in Jujube

Jujube (Ziziphus jujuba Mill.) is rich in valuable bioactive triterpenoids. However, the regulatory mechanism underlying triterpenoid biosynthesis in jujube remains poorly studied. Here, we characterized the triterpenoid content in wild jujube and cultivated jujube. The triterpenoid content was higher in wild jujube than in cultivated jujube, triterpenoids were most abundant in young leaves, buds, and later stages of development. The transcriptome analysis and correlation analysis showed that differentially expressed genes (DEGs) were enriched in the terpenoid synthesis pathways, and triterpenoids content was strongly correlated with farnesyl diphosphate synthase (ZjFPS), squalene synthase (ZjSQS), and transcription factors ZjMYB39 and ZjMYB4 expression. Gene overexpression and silencing analysis indicated that ZjFPS and ZjSQS were key genes in triterpenoid biosynthesis and transcription factors ZjMYB39 and ZjMYB4 regulated triterpenoid biosynthesis. Subcellular localization experiments showed that ZjFPS and ZjSQS were localized to the nucleus and endoplasmic reticulum and ZjMYB39 and ZjMYB4 were localized to the nucleus. Yeast one-hybrid, glucuronidase activity, and dual-luciferase activity assays suggested that ZjMYB39 and ZjMYB4 regulate triterpenoid biosynthesis by directly binding and activating the promoters of ZjFPS and ZjSQS. These findings provide insights into the underlying regulatory network of triterpenoids metabolism in jujube and lay theoretical and practical foundation for molecular breeding.

Methyl Jasmonate- and Salicylic Acid-Induced Transcription Factor ZjWRKY18 Regulates Triterpenoid Accumulation and Salt Stress Tolerance in Jujube

Triterpenoids are important, pharmacologically active substances in jujube (Ziziphus jujuba Mill.), and play an important role in the plant's resistance to abiotic stress. However, regulation of their biosynthesis, and the underlying mechanism of their balance with stress resistance, remain poorly understood. In this study, we screened and functionally characterized the ZjWRKY18 transcription factor, which is associated with triterpenoid accumulation. The transcription factor is induced by methyl jasmonate and salicylic acid, and its activity was observed by gene overexpression and silencing experiments, combined with analyses of transcripts and metabolites. ZjWRKY18 gene silencing decreased the transcription of triterpenoid synthesis pathway genes and the corresponding triterpenoid content. Overexpression of the gene promoted the biosynthesis of jujube triterpenoids, as well as triterpenoids in tobacco and Arabidopsis thaliana. In addition, ZjWRKY18 binds to W-box sequences to activate promoters of 3-hydroxy-3-methyl glutaryl coenzyme A reductase and farnesyl pyrophosphate synthase, suggesting that ZjWRKY18 positively regulates the triterpenoid synthesis pathway. Overexpression of ZjWRKY18 also increased tolerance to salt stress in tobacco and Arabidopsis thaliana. These results highlight the potential use of ZjWRKY18 to improve triterpenoid biosynthesis and salt stress tolerance in plants, and provide a strong basis for metabolic engineering to improve the content of triterpenoids and breeding of jujube varieties that are resistant to stress.

Integrative analysis of microRNAs and mRNAs reveals the regulatory networks of triterpenoid saponin metabolism in Soapberry (Sapindus mukorossi Gaertn.)

Triterpenoid saponin are important secondary metabolites and bioactive constituents of soapberry (Sapindus mukorossi Gaertn.) and are widely used in medicine and toiletry products. However, little is known about the roles of miRNAs in the regulation of triterpenoid saponin biosynthesis in soapberry. In this study, a total of 3036 miRNAs were identified, of which 1372 miRNAs were differentially expressed at different stages of pericarp development. Important KEGG pathways, such as terpenoid backbone biosynthesis, sesquiterpenoid and triterpenoid biosynthesis, and basal transcription factors were highlighted, as well the roles of some key miRNAs, such as ath-miR5021, han-miR3630-3p, and ppe-miR858, which may play important roles in regulating triterpenoid saponin biosynthesis. In addition, 58 miRNAs might participate in saponin biosynthesis pathways by predicting the targets of those miRNAs to 53 saponin biosynthesis structural genes. And 75 miRNAs were identified to potentially play vital role in saponin accumulation by targeting transcript factor genes, bHLH, bZIP, ERF, MYB, and WRKY, respectively, which are candidate regulatory genes in the pathway of saponin biosynthesis. The results of weighted gene coexpression network analysis (WGCNA) suggested that two saponin-specific miRNA modules and 10 hub miRNAs may participate in saponin biosynthesis. Furthermore, multiple miRNA-mRNA regulatory networks potentially involved in saponin biosynthesis were generated, e.g., ath-miR5021-SmIDI2/SmGPS5/SmbAS1/SmCYP71D-3/SmUGT74G-2, han-miR3630-3p-SmCYP71A-14/SmbHLH54/SmMYB135/SmWRKY32, and ppe-miR858-SmMYB5/SmMYB32. qRT-PCR analysis validated the expression patterns of nine miRNAs and 12 corresponding target genes. This study represents the first comprehensive analysis of miRNAs in soapberry and lays the foundation for further understanding of miRNA-based regulation in triterpenoid saponin biosynthesis.

Genome-Wide Identification of SnRK1 Catalytic α Subunit and FLZ Proteins in Glycyrrhiza inflata Bat. Highlights Their Potential Roles in Licorice Growth and Abiotic Stress Responses

Sucrose non-fermenting-1-related protein kinase-1 (SnRK1) and its scaffolding proteins, FCS-like zinc finger proteins (FLZs), are well conserved in land plants and involved in various processes of plant growth and stress responses. Glycyrrhiza inflata Bat. is a widely used licorice species with strong abiotic stress resistance, in which terpenoids and flavonoids are the major bioactive components. Here, we identified 2 SnRK1 catalytic α subunit encoding genes (GiSnRK1α1 and GiSnRK1α2) and 21 FLZ genes in G. inflata. Polygenetic analysis showed that the 21 GiFLZs could be divided into three groups. A total of 10 representative GiFLZ proteins interact with GiSnRK1α1, and they display overlapped subcellular localization (mainly in the nucleus and the cytoplasm) when transiently expressed in Nicotiana benthamiana leaf cells. Coinciding with the existence of various phytohormone-responsive and stress-responsive cis-regulatory elements in the GiSnRK1α and GiFLZ gene promoters, GiFLZs are actively responsive to methyl jasmonic acid (MeJA) and abscisic acid (ABA) treatments, and several GiFLZs and GiSnRK1α1 are regulated by drought and saline-alkaline stresses. Interestingly, GiSnRK1α and 20 of 21 GiFLZs (except for GiFLZ2) show higher expression in the roots than in the leaves. These data provide comprehensive information on the SnRK1 catalytic α subunit and the FLZ proteins in licorice for future functional characterization.

Saponin Biosynthesis in Pulses

Pulses are a group of leguminous crops that are harvested solely for their dry seeds. As the demand for plant-based proteins grows, pulses are becoming important food crops worldwide. In addition to being a rich source of nutrients, pulses also contain saponins that are traditionally considered anti-nutrients, and impart bitterness and astringency. Saponins are plant secondary metabolites with great structural and functional diversity. Given their diverse functional properties and biological activities, both undesirable and beneficial, saponins have received growing attention. It can be expected that redirecting metabolic fluxes to control the saponin levels and produce desired saponins would be an effective approach to improve the nutritional and sensory quality of the pulses. However, little effort has been made toward understanding saponin biosynthesis in pulses, and, thus there exist sizable knowledge gaps regarding its pathway and regulatory network. In this paper, we summarize the research progress made on saponin biosynthesis in pulses. Additionally, phylogenetic relationships of putative biosynthetic enzymes among multiple pulse species provide a glimpse of the evolutionary routes and functional diversification of saponin biosynthetic enzymes. The review will help us to advance our understanding of saponin biosynthesis and aid in the development of molecular and biotechnological tools for the systematic optimization of metabolic fluxes, in order to produce the desired saponins in pulses.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

Genetic divergence in transcriptional regulators of defense metabolism: insight into plant domestication and improvement

A number of mutational changes in transcriptional regulators of defense metabolism have occurred during plant domestication and improvement. Plant domestication and improvement entail genetic changes that underlie divergence in development and metabolism, providing a tremendous model of biological evolution. Plant metabolism produces numerous specialized alkaloids, terpenoids, phenolics, and cyanogenic glucosides with indispensable roles in defense against herbivory and microbial infection. Many compounds toxic or deterrent to predators have been eliminated through domestication and breeding. Series of genes involved in defense metabolism are coordinately regulated by transcription factors that specifically recognize cis-regulatory elements in promoter regions of downstream target genes. Recent developments in DNA sequencing technologies and genomic approaches have facilitated studies of the metabolic and genetic changes in chemical defense that have occurred via human-mediated selection, many of which result from mutations in transcriptional regulators of defense metabolism. In this article, we review such examples in almond (Prunus dulcis), cucumber (Cucumis sativus), pepper (Capsicum spp.), potato (Solanum tuberosum), quinoa (Chenopodium quinoa), sorghum (Sorghum bicolor), and related species and discuss insights into the evolution and regulation of metabolic pathways for specialized defense compounds.

A collection of inducible transcription factor-glucocorticoid receptor fusion lines for functional analyses in Arabidopsis thaliana

Arabidopsis possesses approximately 2000 transcription factors (TFs) in its genome. They play pivotal roles in various biological processes but analysis of their function has been hampered by the overlapping nature of their activities. To uncover clues to their function, we generated inducible TF lines using glucocorticoid receptor (GR) fusion techniques in Arabidopsis. These TF-GR lines each express one of 1255 TFs as a fusion with the GR gene. An average 14 lines of T₂ transgenic TF-GR lines were generated for each TF to monitor their function. To evaluate these transcription lines, we induced the TF-GR lines of phytochrome-interacting factor 4, which controls photomorphogenesis, with synthetic glucocorticoid dexamethasone. These phytochrome-interacting factor 4-GR lines showed the phenotype described in a previous report. We performed screening of the other TF-GR lines for TFs involved in light signaling under blue and far-red light conditions and identified 13 novel TF candidates. Among these, we found two lines showing higher anthocyanin accumulation under light conditions and we examined the regulating genes. These results indicate that the TF-GR lines can be used to dissect functionally redundant genes in plants and demonstrate that the TF-GR line collection can be used as an effective tool for functional analysis of TFs.

Chromosome-level genome of Entada phaseoloides provides insights into genome evolution and triterpenoid saponins biosynthesis

As a medicinal herbal plant, Entada phaseoloides has high levels of secondary metabolites, particularly triterpenoid saponins, which are important resources for scientific research and medical applications. However, the lack of a reference genome for this genus has limited research on its evolution and utilization of its medicinal potential. In this study, we report a chromosome-scale genome assembly for E. phaseoloides using Illumina, Nanopore long reads, and high-throughput chromosome conformation capture technology. The assembled reference genome is 456.18 Mb (scaffold N50=30.9 Mb; contig N50=6.34 Mb) with 95.71 % of the sequences anchored onto 14 pseudochromosomes. E. phaseoloides was estimated to diverge from the Leguminosae lineage at approximately 72.0 million years ago. With the integration of transcriptomic and metabolomic data, gene expression patterns and metabolite profiling of E. phaseoloides were determined in different tissues. The pattern of gene expression and metabolic profile of the kernel were distinct from those of other tissues. Furthermore, the evolution of certain gene families involved in the biosynthesis of triterpenoid saponins and terpenes was analyzed and offer new insights into the formation of these two metabolites. Four CYP genes, one UGT gene and related transcription factors were identified as candidate genes contributing to regulation of triterpenoid saponins biosynthesis. As the first high-quality assembled reference genome in the genus Entada, it will not only provide new information for the evolutionary study of this genus and conservation biology of E. phaseoloides but also lay a foundation for the formation and utilization of secondary metabolites in medicinal plants.

De novo assembly and analysis of Polygonatum cyrtonema Hua and identification of genes involved in polysaccharide and saponin biosynthesis

Background

The investigation of molecular mechanisms involved in polysaccharides and saponin metabolism is critical for genetic engineering of Polygonatum cyrtonema Hua to raise major active ingredient content. Up to now, the transcript sequences are available for different tissues of P. cyrtonema, a wide range scanning about temporal transcript at different ages’ rhizomes was still absent in P. cyrtonema.

Results

Transcriptome sequencing for rhizomes at different ages was performed. Sixty-two thousand six hundred thirty-five unigenes were generated by assembling transcripts from all samples. A total of 89 unigenes encoding key enzymes involved in polysaccharide biosynthesis and 56 unigenes encoding key enzymes involved in saponin biosynthesis. The content of total polysaccharide and total saponin was positively correlated with the expression patterns of mannose-6-phosphate isomerase (MPI), GDP-L-fucose synthase (TSTA3), UDP-apiose/xylose synthase (AXS), UDP-glucose 6-dehydrogenase (UGDH), Hydroxymethylglutaryl CoA synthase (HMGS), Mevalonate kinase (MVK), 2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase (ispF), (E)-4-hydroxy-3-methylbut-2-enyl-diphosphate synthase (ispG), 4-hydroxy-3-methylbut-2-enyl diphosphate reductase (ispH), Farnesyl diphosphate synthase (FPPS). Finally, a number of key genes were selected and quantitative real-time PCR were performed to validate the transcriptome analysis results.

Conclusions

These results create the link between polysaccharides and saponin biosynthesis and gene expression, provide insight for underlying key active substances, and reveal novel candidate genes including TFs that are worth further exploration for their functions and values.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08421-y.

Integrated Analysis of Basic Helix Loop Helix Transcription Factor Family and Targeted Terpenoids Reveals Candidate AarbHLH Genes Involved in Terpenoid Biosynthesis in Artemisia argyi

Artemisia argyi is a valuable traditional medicinal plant in Asia. The essential oil from its leaves is rich in terpenoids and has been used to enhance health and well-being. In China, the market scale of industries related to A. argyi has attained tens of billions of Chinese Yuan. The basic helix-loop-helix (bHLH) family is one of the largest transcription factors families in plants that plays crucial roles in diverse biological processes and is an essential regulatory component of terpenoid biosynthesis. However, the bHLH TFs and their regulatory roles in A. argyi remain unknown. Here, 53 AarbHLH genes were identified from the transcriptome of A. argyi and were classified into 15 subfamilies based on the classification of bHLH proteins in Arabidopsis thaliana. The MEME analysis showed that the conserved motif 1 and motif 2 constituted the most conserved bHLH domain and distributed in most AarbHLH proteins. Additionally, integrated analysis of the expression profiles of AarbHLH genes and the contents of targeted terpenoids in different tissues group and JA-treated group were performed. Eleven up-regulated AarbHLHs and one down-regulated AarbHLH were screened as candidate genes that may participate in the regulation of terpenoid biosynthesis (TPS-AarbHLHs). Correlation analysis between gene expression and terpenoid contents indicated that the gene expression of these 12 TPS-AarbHLHs was significantly correlated with the content changes of 1,8-cineole or β-caryophyllene. Protein-protein interaction networks further illustrated that these TPS-AarbHLHs might be involved in terpenoid biosynthesis in A. argyi. This finding provides a basis to further investigate the regulation mechanism of AarbHLH genes in terpenoid biosynthesis, and will be helpful to improve the quality of A. argyi.

Comprehensive identification of bHLH transcription factors in Litsea cubeba reveals candidate gene involved in the monoterpene biosynthesis pathway

Litsea cubeba (Lour.) Person, an economically important aromatic plant producing essential oils, has lemon-like fragrance and 96.44-98.44% monoterpene contents. bHLH transcription factor plays an important role in plant secondary metabolism and terpene biosynthesis. In this study, we used bioinformatics to identify bHLH transcription factors in L. cubeba, 173 bHLH genes were identified from L. cubeba and divided these into 26 subfamilies based on phylogenetic analysis. The majority of bHLHs in each subfamily shared comparable structures and motifs. While LcbHLHs were unevenly distributed across 12 chromosomes, 10 tandem repeats were discovered. Expression profiles of bHLH genes in different tissues demonstrated that LcbHLH78 is a potential candidate gene for regulating monoterpene biosynthesis. LcbHLH78 and the terpene synthase LcTPS42 showed comparable expression patterns in various tissues and fruit development stages of L. cubeba. Subcellular localization analysis revealed that LcbHLH78 protein localizes to the nucleus, consistent with a transcription factor function. Importantly, transient overexpression of LcbHLH78 increased geraniol and linalol contents. Our research demonstrates that LcbHLH78 enhances terpenoid biosynthesis. This finding will be beneficial for improving the quality of L. cubeba and provides helpful insights for further research into the control mechanism of LcbHLH genes over terpenoid biosynthesis.

Bridging the Gaps between Plant and Human Health: A Systematic Review of Soyasaponins

Saponins, prominent secondary plant metabolites, are recognized for their roles in plant defense and medicinal benefits. Soyasaponins, commonly derived from legumes, are a class of triterpenoid saponins that demonstrate significant potential for plant and human health applications. Previous research and reviews largely emphasize human health effects of soyasaponins. However, the biological effects of soyasaponins and their implications for plants in the context of human health have not been well-discussed. This review provides comprehensive discussions on the biological roles of soyasaponins in plant defense and rhizosphere microbial interactions; biosynthetic regulation and compound production; immunological effects and potential for therapeutics; and soyasaponin acquisition attributed to processing effects, bioavailability, and biotransformation processes based on recent soyasaponin research. Given the multifaceted biological effects elicited by soyasaponins, further research warrants an integrated approach to understand molecular mechanisms of regulations in their production as well as their applications in plant and human health.

Transcription Factors in Alkaloid Engineering

Plants produce a large variety of low-molecular-weight and specialized secondary compounds. Among them, nitrogen-containing alkaloids are the most biologically active and are often used in the pharmaceutical industry. Although alkaloid chemistry has been intensively investigated, characterization of alkaloid biosynthesis, including biosynthetic enzyme genes and their regulation, especially the transcription factors involved, has been relatively delayed, since only a limited number of plant species produce these specific types of alkaloids in a tissue/cell-specific or developmental-specific manner. Recent advances in molecular biology technologies, such as RNA sequencing, co-expression analysis of transcripts and metabolites, and functional characterization of genes using recombinant technology and cutting-edge technology for metabolite identification, have enabled a more detailed characterization of alkaloid pathways. Thus, transcriptional regulation of alkaloid biosynthesis by transcription factors, such as basic helix-loop-helix (bHLH), APETALA2/ethylene-responsive factor (AP2/ERF), and WRKY, is well elucidated. In addition, jasmonate signaling, an important cue in alkaloid biosynthesis, and its cascade, interaction of transcription factors, and post-transcriptional regulation are also characterized and show cell/tissue-specific or developmental regulation. Furthermore, current sequencing technology provides more information on the genome structure of alkaloid-producing plants with large and complex genomes, for genome-wide characterization. Based on the latest information, we discuss the application of transcription factors in alkaloid engineering.

Biosynthesis of Soyasapogenol B by Engineered Saccharomyces cerevisiae

Soyasapogenol B is an oleanane-type pentacyclic triterpene that has various applications in food and healthcare and has a higher biological activity than soyasaponin. Saccharomyces cerevisiae is a potential platform for terpenoid production with mature genetic tools for metabolic pathway manipulation. In this study, we developed a biosynthesis method to produce soyasapogenol B. First, we expressed β-amyrin synthase derived from Glycyrrhiza glabra in S. cerevisiae to generate β-amyrin, as the precursor of soyasapogenol B. Several different types of promoters were then used to regulate the expression of key genes in the mevalonate pathway (MVA), and this subsequently increased the yield of β-amyrin to 17.6 mg/L, 25-fold more than that produced in the original strain L01 (0.68 mg/L). Then, using the β-amyrin-producing strain, we expressed soyasapogenol B synthases from Medicago truncatula (CYP93E2 and CYP72A61V2) and from G. glabra (CYP93E3 and CYP72A566). Soyasapogenol B yields were then optimized by using soyasapogenol B synthases and cytochrome P450 reductase from G. glabra. The most effective soyasapogenol B production strain was used for fermentation, and the yield of soyasapogenol B reached 2.9 mg/L in flask and 8.36 mg/L in a 5-L bioreactor with fed glucose and ethanol. This study demonstrated the heterologous synthesis of soyasapogenol B in S. cerevisiae using the combined expression of CYP93E3 and CYP72A566 in the synthesis pathway, which significantly increased the production of soyasapogenol B and provides a reference method for the biosynthesis of other triterpenes.

A post processing strategy to score and rank the annotation confidence of saponins in natural products by integrating MS2 spectral similarity and fragment interpretation

Tandem mass spectrometry-spectra-based annotation in natural products challenges a lot because of ambiguous structural characterization. It still lacks an efficiency method to score and rank the annotation confidence. Herein, we develop a novel approach to rank the annotation confidences of saponins. Annotations were accomplished according to fragmentation patterns. The corresponding diagnostic fragments and their abundances were recorded. Average abundances were taken as a reference spectrum, and the cosine similarity score (CSS) was calculated to measure how well the spectral matched. According to CSS values, statistic description for confidence levels can be effectively provided. Next, the fragment interpretation score (FIS) was proposed to investigate the deviators' characteristic fragmentation. FIS offset the effect from the deviators' unique fragments. Suspicious annotations involving low CSS and high FIS, may derived from the MS² spectral background interferences or co-elution. Annotations with low CSS and FIS rank as low confidences, as these annotations need more attention. Using this method, novel saccharide sequences, specific fragmentation preferences, undistinguished precursors, even new structures can also be well traced. By proposed new scoring system, confidence evaluations can be ranked, resulting in significantly enhanced annotation reliability.

Characterization of UDP-glucose dehydrogenase isoforms in the medicinal legume Glycyrrhiza uralensis

Uridine 5'-diphosphate (UDP)-glucose dehydrogenase (UGD) produces UDP-glucuronic acid from UDP-glucose as a precursor of plant cell wall polysaccharides. UDP-glucuronic acid is also a sugar donor for the glycosylation of various plant specialized metabolites. Nevertheless, the roles of UGDs in plant specialized metabolism remain poorly understood. Glycyrrhiza species (licorice), which are medicinal legumes, biosynthesize triterpenoid saponins, soyasaponins and glycyrrhizin, commonly glucuronosylated at the C-3 position of the triterpenoid scaffold. Often, several different UGD isoforms are present in plants. To gain insight into potential functional differences among UGD isoforms in triterpenoid saponin biosynthesis in relation to cell wall component biosynthesis, we identified and characterized Glycyrrhiza uralensis UGDs (GuUGDs), which were discovered to comprise five isoforms, four of which (GuUGD1-4) showed UGD activity in vitro. GuUGD1-4 had different biochemical properties, including their affinity for UDP-glucose, catalytic constant, and sensitivity to feedback inhibitors. GuUGD2 had the highest catalytic constant and highest gene expression level among the GuUGDs, suggesting that it is the major isoform contributing to the transition from UDP-glucose to UDP-glucuronic acid in planta. To evaluate the contribution of GuUGD isoforms to saponin biosynthesis, we compared the expression patterns of GuUGDs with those of saponin biosynthetic genes in methyl jasmonate (MeJA)-treated cultured stolons. GuUGD1-4 showed delayed responses to MeJA compared to those of saponin biosynthetic genes, suggesting that MeJA-responsive expression of GuUGDs compensates for the decreased UDP-glucuronic acid pool due to consumption during saponin biosynthesis.

Full-length transcriptome sequences by a combination of sequencing platforms applied to isoflavonoid and triterpenoid saponin biosynthesis of Astragalus mongholicus Bunge

BACKGROUND

Astragalus mongholicus Bunge is an important medicinal plant used in traditional Chinese medicine. It is rich in isoflavonoids and triterpenoid saponins. Although these active constituents of A. mongholicus have been discovered for a long time, the genetic basis of isoflavonoid and triterpenoid saponin biosynthesis in this plant is virtually unknown because of the lack of a reference genome. Here, we used a combination of next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing to identify genes involved in the biosynthetic pathway of secondary metabolites in A. mongholicus.

RESULTS

In this study, NGS, SMRT sequencing, and targeted compound analysis were combined to investigate the association between isoflavonoid and triterpenoid saponin content, and specific gene expression in the root, stem, and leaves of A. mongholicus. Overall, 643,812 CCS reads were generated, yielding 121,107 non-redundant transcript isoforms with an N50 value of 2124 bp. Based on these highly accurate transcripts, 104,756 (86.50%) transcripts were successfully annotated by any of the seven databases (NR, NT, Swissprot, KEGG, KOG, Pfam and GO). Levels of four isoflavonoids and four astragalosides (triterpenoid saponins) were determined. Forty-four differentially expressed genes (DEGs) involved in isoflavonoid biosynthesis and 44 DEGs from 16 gene families that encode enzymes involved in triterpenoid saponin biosynthesis were identified. Transcription factors (TFs) associated with isoflavonoid and triterpenoid saponin biosynthesis, including 72 MYBs, 53 bHLHs, 64 AP2-EREBPs, and 11 bZIPs, were also identified. The above transcripts showed different expression trends in different plant organs.

CONCLUSIONS

This study provides important genetic information on the A. mongholicus genes that are essential for isoflavonoid and triterpenoid saponin biosynthesis, and provides a basis for developing the medicinal value of this plant.

Spatio-temporal selection of reference genes in the two congeneric species of Glycyrrhiza

Glycyrrhiza, a genus of perennial medicinal herbs, has been traditionally used to treat human diseases, including respiratory disorders. Functional analysis of genes involved in the synthesis, accumulation, and degradation of bioactive compounds in these medicinal plants requires accurate measurement of their expression profiles. Reverse transcription quantitative real-time PCR (RT-qPCR) is a primary tool, which requires stably expressed reference genes to serve as the internal references to normalize the target gene expression. In this study, the stability of 14 candidate reference genes from the two congeneric species G. uralensis and G. inflata, including ACT, CAC, CYP, DNAJ, DREB, EF1, RAN, TIF1, TUB, UBC2, ABCC2, COPS3, CS, R3HDM2, were evaluated across different tissues and throughout various developmental stages. More importantly, we investigated the impact of interactions between tissue and developmental stage on the performance of candidate reference genes. Four algorithms, including geNorm, NormFinder, BestKeeper, and Delta Ct, were used to analyze the expression stability and RefFinder, a comprehensive software, provided the final recommendation. Based on previous research and our preliminary data, we hypothesized that internal references for spatio-temporal gene expression are different from the reference genes suited for individual factors. In G. uralensis, the top three most stable reference genes across different tissues were R3HDM2, CAC and TUB, while CAC, CYP and ABCC2 were most suited for different developmental stages. CAC is the only candidate recommended for both biotic factors, which is reflected in the stability ranking for the spatio (tissue)-temporal (developmental stage) interactions (CAC, R3HDM2 and DNAJ). Similarly, in G. inflata, COPS3, R3HDM2 and DREB were selected for tissues, while RAN, COPS3 and CS were recommended for developmental stages. For the tissue-developmental stage interactions, COPS3, DREB and ABCC2 were the most suited reference genes. In both species, only one of the top three candidates was shared between the individual factors and their interactions, specifically, CAC in G. uralensis and COPS3 in G. inflata, which supports our overarching hypothesis. In summary, spatio-temporal selection of reference genes not only lays the foundation for functional genomics research in Glycyrrhiza, but also facilitates these traditional medicinal herbs to reach/maximize their pharmaceutical potential.

The Heat Shock Protein 40-Type Chaperone MASH Supports the Endoplasmic Reticulum-Associated Degradation E3 Ubiquitin Ligase MAKIBISHI1 in Medicago truncatula

Jasmonates (JA) are oxylipin-derived phytohormones that trigger the production of specialized metabolites that often serve in defense against biotic stresses. In Medicago truncatula, a JA-induced endoplasmic reticulum-associated degradation (ERAD)-type machinery manages the production of bioactive triterpenes and thereby secures correct plant metabolism, growth, and development. This machinery involves the conserved RING membrane-anchor (RMA)-type E3 ubiquitin ligase MAKIBISHI1 (MKB1). Here, we discovered two additional members of this protein control apparatus via a yeast-based protein-protein interaction screen and characterized their function. First, a cognate E2 ubiquitin-conjugating enzyme was identified that interacts with MKB1 to deliver activated ubiquitin and to mediate its ubiquitination activity. Second, we identified a heat shock protein 40 (HSP40) that interacts with MKB1 to support its activity and was therefore designated MKB1-supporting HSP40 (MASH). MASH expression was found to be co-regulated with that of MKB1. The presence of MASH is critical for MKB1 and ERAD functioning because the dramatic morphological, transcriptional, and metabolic phenotype of MKB1 knock-down M. truncatula hairy roots was phenocopied by silencing of MASH. Interaction was also observed between the Arabidopsis thaliana (Arabidopsis) homologs of MASH and MKB1, suggesting that MASH represents an essential and plant-specific component of this vital and conserved eukaryotic protein quality control machinery.

Insights into the diversification of subclade IVa bHLH transcription factors in Fabaceae

BACKGROUND

Fabaceae plants appear to contain larger numbers of subclade IVa basic-helix-loop-helix (bHLH) transcription factors than other plant families, and some members of this subclade have been identified as saponin biosynthesis regulators. We aimed to systematically elucidate the diversification of this subclade and obtain insights into the evolutionary history of saponin biosynthesis regulation in Fabaceae.

RESULTS

In this study, we collected sequences of subclade IVa bHLH proteins from 40 species, including fabids and other plants, and found greater numbers of subclade IVa bHLHs in Fabaceae. We confirmed conservation of the bHLH domain, C-terminal ACT-like domain, and exon-intron organisation among almost all subclade IVa members in model legumes, supporting the results of our classification. Phylogenetic tree-based classification of subclade IVa revealed the presence of three different groups. Interestingly, most Fabaceae subclade IVa bHLHs fell into group 1, which contained all legume saponin biosynthesis regulators identified to date. These observations support the co-occurrence and Fabaceae-specific diversification of saponin biosynthesis regulators. Comparing the expression of orthologous genes in Glycine max, Medicago truncatula, and Lotus japonicus, orthologues of MtTSAR1 (the first identified soyasaponin biosynthesis regulatory transcription factor) were not expressed in the same tissues, suggesting that group 1 members have gained different expression patterns and contributions to saponin biosynthesis during their duplication and divergence. On the other hand, groups 2 and 3 possessed fewer members, and their phylogenetic relationships and expression patterns were highly conserved, indicating that their activities may be conserved across Fabaceae.

CONCLUSIONS

This study suggests subdivision and diversification of subclade IVa bHLHs in Fabaceae plants. The results will be useful for candidate selection of unidentified saponin biosynthesis regulators. Furthermore, the functions of groups 2 and 3 members are interesting targets for clarifying the evolution of subclade IVa bHLH transcription factors in Fabaceae.

Current advances in environmental stimuli regulating the glycyrrhizic acid biosynthesis pathway

Glycyrrhizic acid, the main active ingredient of licorice, has good antibacterial, anti-tumor, anti-viral, anti-inflammatory, and immunostimulatory activities. However, the content of glycyrrhizic acid fluctuates greatly in different licorice cultivars, and production depends on plant sources, which greatly limits its development and applications. Therefore, increasing glycyrrhizic acid content has become a research priority. In recent years, regulation of the glycyrrhizic acid biosynthesis pathway has been analyzed, the downstream synthesis pathway in licorice has been fully investigated, some key genes have been cloned, polymorphisms have been studied, and the content of glycyrrhizic acid was shown to be regulated by environmental stimuli. This work has provided a basis for studying the regulation mechanism of the glycyrrhizic acid synthesis pathway. This review summarizes and discusses relevant research to provide a current understanding of the glycyrrhizic acid synthesis pathway and its regulation in licorice.

The Arabidopsis transcription factor NAI1 activates the NAI2 promoter by binding to the G-box motifs

Brassicaceae plants, including Arabidopsis thaliana, develop endoplasmic reticulum (ER)-derived structures called ER bodies, which are involved in chemical defense against herbivores. NAI1 is a basic helix-loop-helix (bHLH) type transcription factor that regulates two downstream genes, NAI2 and BGLU23, that are responsible for the ER body formation and function. Here, we examined the transcription factor function of NAI1, and found that NAI1 binds to the promoter region of NAI2 and activates the NAI2 promoter. The recombinant NAI1 protein recognizes the canonical and non-canonical G-box motifs in the NAI2 promoter. Furthermore, we examined the DNA binding activity of NAI1 toward several E-box motifs in the NAI2 and BGLU23 promoters and found that NAI1 binds to a DNA fragment that includes an E-box motif from the BGLU23 promoter. Subcellular localization of NAI1 was evident in the nucleus, which is consistent with its transcription factor function. Transient expression experiments in Nicotiana benthamiana leaves showed that GFP-NAI1 protein activated the NAI2 promoter by binding to the two G-boxes of the promoter. Disruption of the G-boxes abolished the NAI1-dependent activation of the NAI2 promoter. These results indicate that NAI1 has a DNA binding activity in a motif-dependent manner and suggest that NAI1 regulates NAI2 and BGLU23 gene expressions through binding to these DNA motifs in their promoters.

A cellulose synthase-derived enzyme catalyses 3-O-glucuronosylation in saponin biosynthesis

Triterpenoid saponins are specialised metabolites distributed widely in the plant kingdom that consist of one or more sugar moieties attached to triterpenoid aglycones. Despite the widely accepted view that glycosylation is catalysed by UDP-dependent glycosyltransferase (UGT), the UGT which catalyses the transfer of the conserved glucuronic acid moiety at the C-3 position of glycyrrhizin and various soyasaponins has not been determined. Here, we report that a cellulose synthase superfamily-derived glycosyltransferase (CSyGT) catalyses 3-O-glucuronosylation of triterpenoid aglycones. Gene co-expression analyses of three legume species (Glycyrrhiza uralensis, Glycine max, and Lotus japonicus) reveal the involvement of CSyGTs in saponin biosynthesis, and we characterise CSyGTs in vivo using Saccharomyces cerevisiae. CSyGT mutants of L. japonicus do not accumulate soyasaponin, but the ectopic expression of endoplasmic reticulum membrane-localised CSyGTs in a L. japonicus mutant background successfully complement soyasaponin biosynthesis. Finally, we produced glycyrrhizin de novo in yeast, paving the way for sustainable production of high-value saponins.

Advances in biosynthesis of triterpenoid saponins in medicinal plants

In recent years, biosynthesis of triterpenoid saponins in medicinal plants has been widely studied because of their active ingredients with diverse pharmacological activities. Various oxidosqualene cyclases, cytochrome P450 monooxygenases, uridine diphosphate glucuronosyltransferases, and transcription factors related to triterpenoid saponins biosynthesis have been explored and identified. In the biosynthesis of triterpenoid saponins, the progress of gene mining by omics-based sequencing, gene screening, gene function verification, catalyzing mechanism of key enzymes and gene regulation are summarized and discussed. By the progress of the biosynthesis pathway of triterpenoid saponins, the large-scale production of some triterpenoid saponins and aglycones has been achieved through plant tissue culture, transgenic plants and engineered yeast cells. However, the complex biosynthetic pathway and structural diversity limit the biosynthesis of triterpenoid saponins in different system. Special focus can further be placed on the systematic botany information of medicinal plants obtained from omics large dataset, and triterpenoid saponins produced by synthetic biology strategies, gene mutations and gene editing technology.

The MYB transcription factor CiMYB42 regulates limonoids biosynthesis in citrus

BACKGROUND

Limonoids are major bioactive compounds that are produced by the triterpenoid metabolic pathway. The detailed biochemical process of limonoid biosynthesis and the mechanism of its molecular regulation remain elusive. The identification of transcription factors that regulate limonoid biosynthetic pathways is very important for understanding the underlying regulatory mechanisms. This information could also provide tools for manipulating biosynthesis genes to modulate limonoid production.

RESULTS

In this study, the CiMYB42 transcription factor was isolated to identify its role in limonoid biosynthesis. Multiple alignment analysis and phylogenetic analysis demonstrated that CiMYB42 is a typical R2R3MYB transcription factor that shares high similarity of its amino acid sequence with AtMYB42. Limonoids contents were higher in Citrus sinensis and Citrus grandis than in other species. Limonoid accumulation during leaf development also showed diverse trends in different genotypes. The expression of CiMYB42 was significantly related to the limonoid content and the expression of CiOSC in some citrus accessions. The overexpression of CiMYB42 in sweet orange resulted in significant accumulation of limonin, whereas the downregulation of CiMYB42 by RNAi resulted in a dwarf phenotype and less nomilin accumulation. Furthermore, the results of a yeast one-hybrid assay and EMSA indicated that CiMYB42 binds exclusively to the TTGTTG sequence (type II MYB core) in the promoter of CiOSC. Together, these results suggest that CiMYB42 positively regulates limonoid biosynthesis by regulating the expression of CiOSC by binding to the TTGTTG sequence (type II MYB core) of its promoter.

CONCLUSIONS

CiMYB42 is an important transcription activator involved in limonoid biosynthesis that regulates the expression of CiOSC by binding to the TTGTTG sequence (type II MYB core).

Lotus japonicus Triterpenoid Profile and Characterization of the CYP716A51 and LjCYP93E1 Genes Involved in Their Biosynthesis In Planta

Lotus japonicus is an important model legume plant in several fields of research, such as secondary (specialized) metabolism and symbiotic nodulation. This plant accumulates triterpenoids; however, less information regarding its composition, content and biosynthesis is available compared with Medicago truncatula and Glycine max. In this study, we analyzed the triterpenoid content and composition of L. japonicus. Lotus japonicus accumulated C-28-oxidized triterpenoids (ursolic, betulinic and oleanolic acids) and soyasapogenols (soyasapogenol B, A and E) in a tissue-dependent manner. We identified an oxidosqualene cyclase (OSC) and two cytochrome P450 enzymes (P450s) involved in triterpenoid biosynthesis using a yeast heterologous expression system. OSC9 was the first enzyme derived from L. japonicus that showed α-amyrin (a precursor of ursolic acid)-producing activity. CYP716A51 showed triterpenoid C-28 oxidation activity. LjCYP93E1 converted β-amyrin into 24-hydroxy-β-amyrin, a metabolic intermediate of soyasapogenols. The involvement of the identified genes in triterpenoid biosynthesis in L. japonicus plants was evaluated by quantitative real-time PCR analysis. Furthermore, gene loss-of-function analysis of CYP716A51 and LjCYP93E1 was conducted. The cyp716a51-mutant L. japonicus hairy roots generated by the genome-editing technique produced no C-28 oxidized triterpenoids. Likewise, the complete abolition of soyasapogenols and soyasaponin I was observed in mutant plants harboring Lotus retrotransposon 1 (LORE1) in LjCYP93E1. These results indicate that the activities of these P450 enzymes are essential for triterpenoid biosynthesis in L. japonicus. This study increases our understanding of triterpenoid biosynthesis in leguminous plants and provides information that will facilitate further studies of the physiological functions of triterpenoids using L. japonicus.

High-Throughput Screening and Characterization of a High-Density Soybean Mutant Library Elucidate the Biosynthesis Pathway of Triterpenoid Saponins

Triterpenes (C30) constitute one of the diverse class of natural products with potential applications in food, cosmetic and pharmaceutical industries. Soyasaponins are oleanane-type triterpenoids widespread among legumes and particularly abundant in soybean seeds. They have associated with various pharmacological implications and undesirable taste properties of soybean-based food products. Uncovering the biosynthetic genes of soyasaponins will provide new opportunities to control the pathway for human benefits. However, the pathway of soyasaponin biosynthesis has not been fully elucidated in part because of a paucity of natural mutants. Here, we applied a structured high-density soybean mutant library for the forward genetic screening of triterpenoid biosynthesis. The seed soyasaponin polymorphism in the mutant library was evaluated using a high-throughput thin-layer chromatography and liquid chromatography tandem mass spectrometry analysis. This screening identified 35 mutants (3.85% of 909 mutant lines) with seven unusual soyasaponin phenotypes (Categories 1-7), which was greater than the number of natural mutants reported previously (22 mutants, 0.18% of ∼12,428 accessions). Nine unique intermediates of soyasaponin biosynthesis were identified and their chemical structures were estimated based on their MS/MS fragment patterns. Based on published information, 19 mutants could be associated with loss of function of four individual soyasaponin biosynthesis genes identified through expressed sequence tag mining or positional cloning, whereas the remaining 16 mutants were novel and may facilitate discovery of the unknown biosynthetic genes of soyasaponins. Our approach and library may help to identify new phenotype materials and causative genes associated with specialized metabolite production and other traits.

The Recruitment Model of Metabolic Evolution: Jasmonate-Responsive Transcription Factors and a Conceptual Model for the Evolution of Metabolic Pathways

Plants produce a vast array of structurally diverse specialized metabolites with various biological activities, including medicinal alkaloids and terpenoids, from relatively simple precursors through a series of enzymatic steps. Massive metabolic flow through these pathways usually depends on the transcriptional coordination of a large set of metabolic, transport, and regulatory genes known as a regulon. The coexpression of genes involved in certain metabolic pathways in a wide range of developmental and environmental contexts has been investigated through transcriptomic analysis, which has been successfully exploited to mine the genes involved in various metabolic processes. Transcription factors are DNA-binding proteins that recognize relatively short sequences known as cis-regulatory elements residing in the promoter regions of target genes. Transcription factors have positive or negative effects on gene transcription mediated by RNA polymerase II. Evolutionarily conserved transcription factors of the APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) and basic helix-loop-helix (bHLH) families have been identified as jasmonate (JA)-responsive transcriptional regulators of unrelated specialized pathways in distinct plant lineages. Here, I review the current knowledge and propose a conceptual model for the evolution of metabolic pathways, termed "recruitment model of metabolic evolution." According to this model, structural genes are repeatedly recruited into regulons under the control of conserved transcription factors through the generation of cognate cis-regulatory elements in the promoters of these genes. This leads to the adjustment of catalytic activities that improve metabolic flow through newly established passages.