Function of AP2/ERF Transcription Factors Involved in the Regulation of Specialized Metabolism in Ophiorrhiza pumila Revealed by Transcriptomics and Metabolomics

Metabolite-based genome-wide association studies enable the dissection of the genetic bases of flavonoids, betaine and spermidine in wolfberry (Lycium)

Wolfberry is a plant with medicinal and food values. However, its bioactive ingredients and the corresponding genetic bases have not been determined. Here, we de novo generated a chromosome-level genome assembly for wolfberry, yielding a genome sequence of ~1.77 Gb with contig N50 of 50.55 Mb and 39 224 predicted gene models. A variation map, using 307 re-sequenced accessions, was called based on this genome assembly. Furthermore, the fruit metabolome of these accessions was profiled using 563 annotated metabolites, which separated Lycium barbarum L. and non-L. barbarum L. The flavonoids, coumarins, alkaloids and nicotinic acid contents were higher in the former than in the latter. A metabolite-based genome-wide association study mapped 156 164 significant single nucleotide polymorphisms corresponding to 340 metabolites. This included 19 219 unique lead single nucleotide polymorphisms in 1517 significant association loci, of which three metabolites, flavonoids, betaine and spermidine, were highlighted. Two candidate genes, LbUGT (evm.TU.chr07.2692) and LbCHS (evm.TU.chr07.2738), with non-synonymous mutations, were associated with the flavonoids content. LbCHS is a structural gene that interacts with a nearby MYB transcription factor (evm.TU.chr07.2726) both in L. barbarum and L. ruthenicum. Thus, these three genes might be involved in the biosynthesis/metabolism of flavonoids. LbSSADH (evm.TU.chr09.627) was identified as possibly participating in betaine biosynthesis/metabolism. Four lycibarbarspermidines (E-G and O) were identified, and only the lycibarbarspermidines O content was higher in L. barbarum varieties than in non-L. barbarum varieties. The evm.TU.chr07.2680 gene associated with lycibarbarspermidines O was annotated as an acetyl-CoA-benzylalcohol acetyltransferase, suggesting that it is a candidate gene for spermidine biosynthesis. These results provide novel insights into the specific metabolite profile of non-L. barbarum L. and the genetic bases of flavonoids, betaine and spermidine biosynthesis/metabolism.

The Ginkgo biloba microRNA160-ERF4 module participates in terpene trilactone biosynthesis

Terpene trilactones (TTLs) are important secondary metabolites in ginkgo (Ginkgo biloba); however, their biosynthesis gene regulatory network remains unclear. Here, we isolated a G. biloba ethylene response factor 4 (GbERF4) involved in TTL synthesis. Overexpression of GbERF4 in tobacco (Nicotiana tabacum) significantly increased terpenoid content and upregulated the expression of key enzyme genes (3-hydroxy-3-methylglutaryl-CoA reductase [HMGR], 3-hydroxy-3-methylglutaryl-CoA synthase [HMGS], 1-deoxy-D-xylulose-5-phosphate reductoisomerase [DXR], 1-deoxy-D-xylulose-5-phosphate synthase [DXS], acetyl-CoA C-acetyltransferase [AACT], and geranylgeranyl diphosphate synthase [GGPPS]) in the terpenoid pathway in tobacco, suggesting that GbERF4 functions in regulating the synthesis of terpenoids. The expression pattern analysis and previous microRNA (miRNA) sequencing showed that gb-miR160 negatively regulates the biosynthesis of TTLs. Transgenic experiments showed that overexpression of gb-miR160 could significantly inhibit the accumulation of terpenoids in tobacco. Targeted inhibition and dual-luciferase reporter assays confirmed that gb-miR160 targets and negatively regulates GbERF4. Transient overexpression of GbERF4 increased TTL content in G. biloba, and further transcriptome analysis revealed that DXS, HMGS, CYPs, and transcription factor genes were upregulated. In addition, yeast 1-hybrid and dual-luciferase reporter assays showed that GbERF4 could bind to the promoters of the HMGS1, AACT1, DXS1, levopimaradiene synthase (LPS2), and GGPPS2 genes in the TTL biosynthesis pathway and activate their expression. In summary, this study investigated the molecular mechanism of the gb-miR160-GbERF4 regulatory module in regulating the biosynthesis of TTLs. It provides information for enriching the understanding of the regulatory network of TTL biosynthesis and offers important gene resources for the genetic improvement of G. biloba with high contents of TTLs.

CaERF1- mediated ABA signal positively regulates camptothecin biosynthesis by activating the iridoid pathway in Camptotheca acuminata

Camptotheca acuminata is one of the primary sources of camptothecin (CPT), which is widely used in the treatment of human malignancies because of its inhibitory activity against DNA topoisomerase I. Although several transcription factors have been identified for regulating CPT biosynthesis in other species, such as Ophiorrhiza pumila, the specific regulatory components controlling CPT biosynthesis in C. acuminata have yet to be definitively determined. In this study, CaERF1, an DREB subfamily of the APETALA2/ethylene response factors (AP2ERFs), was identified in C. acuminata. The transient overexpression and silencing of CaERF1 in C. acuminata leaves confirmed that it positively regulates the accumulation of CPT by inducing the expression of CaCYC1 and CaG8O in the iridoid pathway. Results of transient transcriptional activity assay and yeast one-hybrid assays have showed that CaERF1 transcriptionally activates the expression of CaCYC1 and CaG8O by binding to RAA and CEI elements in the promoter regions of these two genes. Furthermore, the expression of CaCYC1 and CaG8O in CaERF1-silenced leaves was less sensitive to ABA treatment, indicating that CaERF1 is a crucial component involved in ABA-regulated CPT biosynthesis in C. acuminata.

Molecular regulation of the key specialized metabolism pathways in medicinal plants

The basis of modern pharmacology is the human ability to exploit the production of specialized metabolites from medical plants, for example, terpenoids, alkaloids, and phenolic acids. However, in most cases, the availability of these valuable compounds is limited by cellular or organelle barriers or spatio-temporal accumulation patterns within different plant tissues. Transcription factors (TFs) regulate biosynthesis of these specialized metabolites by tightly controlling the expression of biosynthetic genes. Cutting-edge technologies and/or combining multiple strategies and approaches have been applied to elucidate the role of TFs. In this review, we focus on recent progress in the transcription regulation mechanism of representative high-value products and describe the transcriptional regulatory network, and future perspectives are discussed, which will help develop high-yield plant resources.

Plant disease resistance outputs regulated by AP2/ERF transcription factor family

Plants have evolved a complex and elaborate signaling network to respond appropriately to the pathogen invasion by regulating expression of defensive genes through certain transcription factors. The APETALA2/ethylene response factor (AP2/ERF) family members have been determined as key regulators in growth, development, and stress responses in plants. Moreover, a growing body of evidence has demonstrated the critical roles of AP2/ERFs in plant disease resistance. In this review, we describe recent advances for the function of AP2/ERFs in defense responses against microbial pathogens. We summarize that AP2/ERFs are involved in plant disease resistance by acting downstream of mitogen activated protein kinase (MAPK) cascades, and regulating expression of genes associated with hormonal signaling pathways, biosynthesis of secondary metabolites, and formation of physical barriers in an MAPK-dependent or -independent manner. The present review provides a multidimensional perspective on the functions of AP2/ERFs in plant disease resistance, which will facilitate the understanding and future investigation on the roles of AP2/ERFs in plant immunity.

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.

Phylogenomics of transcriptionally active AP2/ERF and bHLH transcription factors and study of their promoter regions in Nothapodytes nimmoniana (J.Graham) Mabb

Nothapodytes nimmoniana is a medicinally important plant producing anticancer monoterpene indole alkaloid (MIA), camptothecin (CPT). The CPT is synthesised through the strictosidine intermediate following the MIA pathway; however, transcriptional regulation of CPT pathway is still elusive in N. nimmoniana. Biosynthesis of MIA is regulated by various transcription factors (TFs) belonging to AP2/ERF, bHLH, MYB, and WRKY families. The present study identified transcriptionally active full-length 105 AP2/ERF and 68 bHLH family TFs from the N. nimmoniana. AP2/ERF TFs were divided into three subfamilies along with a soloist, while bHLH TFs were divided into 10 subfamilies according to their phylogenetic similarities. Three group IXa ERFs, Nn-ERF22, Nn-ERF29, and Nn-ERF41, one subfamily IVa TF Nn-bHLH7, and three subfamilies IIIe Nn-bHLH33, Nn-bHLH51, and Nn-bHLH52 clustered with the TFs regulating alkaloid biosynthesis in Catharanthus roseus, tomato, tobacco, and Artemisia annua. Expression of these TFs in N. nimmoniana was higher in roots, which is a primary CPT accumulating tissue. Moreover, genome skimming approach was used to reconstruct the promoter regions of candidate ERF genes to identify the cis-regulatory elements. The presence of G-boxes and other jasmonic acid-responsive elements in the promoter suggests the regulation of ERFs by bHLHs. The present study effectively generated and used genomics resource for characterisation of regulatory TFs from non-model medicinal plant.

Comprehensive analysis of OpHD-ZIP transcription factors related to the regulation of camptothecin biosynthesis in Ophiorrhiza pumila

Ophiorrhiza pumila, as a folk herb belonging to the Rubiaceae family, has become a potential source of camptothecin (CPT), which is a monoterpenoid indole alkaloid with good antitumor property. However, the camptothecin content in this herb is low, and is far from meeting the increasing clinical demand. Understanding the transcriptional regulation of camptothecin biosynthesis provides an effective strategy for improvement of camptothecin yield. Previous studies have demonstrated several transcription factors that are related to camptothecin biosynthesis, while the functions of HD-ZIP members in O. pumila have not been investigated yet. In this study, 32 OpHD-ZIP transcription factor members were genome-wide identified. Phylogenetic tree showed that these OpHD-ZIP proteins are divided into four subfamilies. Based on the transcriptome data, nine OpHD-ZIP genes were shown to be predominantly expressed in O. pumila roots, which were in line with the camptothecin biosynthetic genes. Co-expression analysis showed that OpHD-ZIP7 and OpHD-ZIP20 were potentially related to the modulation of camptothecin biosynthesis. Dual-luciferase reporter assays (Dual-LUC) showed that both OpHD-ZIP7 and OpHD-ZIP20 could activate the expression of camptothecin biosynthetic genes OpIO and OpTDC. In conclusion, this study offered the promising data for exploring the roles of OpHD-ZIP transcription factors in regulating camptothecin biosynthesis.

A liquid chromatography-mass spectrometry-based metabolomics strategy to explore plant metabolic diversity

Plants are expert chemists producing millions of metabolites, only a fraction of which are known to date. Plant metabolomics explores the rationale for highly diverse metabolites evolved and synthesized by plants. Over two-thirds of modern medicines are somehow inspired and/or derived from plants, making the identification of phytochemicals a means of discovering new medicines to challenge existing and emerging diseases. This chapter introduces our established liquid chromatography-tandem mass spectrometry-based untargeted metabolomics approach centered around discovering specialized metabolites (so-called secondary metabolites) across broad lineages of nonmodel plant species. Detecting hundreds to thousands of metabolite peaks, including assigning chemical identity, makes metabolomics data generation and analysis a very complex process. Various mass spectrometry techniques are currently being developed to approach the comprehensive metabolome. Among them, untargeted metabolomics can provide new biological insights by simultaneously and unbiasedly measuring and analyzing all detected metabolites. We have provided a hands-on modular account for untargeted plant metabolomics, from preparing plant biological samples to data analysis and processing using ultra-performance liquid chromatography-quadrupole time-of-flight mass spectrometry. The methods described here offer a foundation and expert opinion on plant metabolome analysis.

OpNAC1 transcription factor regulates the biosynthesis of the anticancer drug camptothecin by targeting loganic acid O-methyltransferase in Ophiorrhiza pumila

Camptothecin (CPT) is an anticancer pentacyclic quinoline alkaloid widely used to treat cancer patients worldwide. However, the biosynthetic pathway and transcriptional regulation of camptothecin are largely unknown. Ophiorrhiza pumila, the herbaceous plant from the Rubiaceae family, has emerged as a model plant for studying camptothecin biosynthesis and regulation. In this study, a high-quality reference genome of O. pumila with estimated size of ~456.90 Mb was reported, and the accumulation level of camptothecin in roots was higher than that in stems and leaves. Based on its spatial distribution in the plant, we examined gene functions and expression by combining genomics with transcriptomic analysis. Two loganic acid O-methyltransferase (OpLAMTs) were identified in strictosidine-producing plant O. pumila, and enzyme catalysis assays showed that OpLAMT1 and not OpLAMT2 could convert loganic acid into loganin. Further knock-out of OpLAMT1 expression led to the elimination of loganin and camptothecin accumulation in O. pumila hairy roots. Four key residues were identified in OpLAMT1 protein crucial for the catalytic activity of loganic acid to loganin. By co-expression network, we identified a NAC transcription factor, OpNAC1, as a candidate gene for regulating camptothecin biosynthesis. Transgenic hairy roots and biochemical assays demonstrated that OpNAC1 suppressed OpLAMT1 expression. Here, we reported on two camptothecin metabolic engineering strategies paving the road for industrial-scale production of camptothecin in CPT-producing plants.

Profiles of Metabolic Genes in Uncaria rhynchophylla and Characterization of the Critical Enzyme Involved in the Biosynthesis of Bioactive Compounds-(iso)Rhynchophylline

Rhynchophylline (RIN) and isorhynchophylline (IRN), two of the representative types of indole alkaloids, showed the unique spiroindole structures produced in Uncaria rhynchophylla. As the bioactive constituent of U. rhynchophylla, IRN has recently drawn extensive attention toward antihypertensive and neuroprotective activities. Despite their medicinal importance and unique chemical structure, the biosynthetic pathways of plant spiroindole alkaloids are still largely unknown. In this study, we used U. rhynchophylla, extensively used in traditional Chinese medicine (TCM), a widely cultivated plant of the Uncaria genus, to investigate the biosynthetic genes and characterize the functional enzymes in the spiroindole alkaloids. We aim to use the transcriptome platform to analyse the tissue-specific gene expression in spiroindole alkaloids-producing tissues, including root, bud, stem bark and leaf. The critical genes involved in the biosynthesis of precursors and scaffold formation of spiroindole alkaloids were discovered and characterized. The datasets from this work provide an essential resource for further investigating metabolic pathways in U. rhynchophylla and facilitate novel functional enzyme characterization and a good biopharming approach to spiroindole alkaloids.

Interplay of transcription factors orchestrating the biosynthesis of plant alkaloids

Plants produce a range of secondary metabolites primarily as defence molecules. A plant has to invest considerable energy to synthesise alkaloids, and sometimes they are even toxic to themselves. Hence, the biosynthesis of alkaloids is a spatiotemporally regulated process under quantitative feedback regulation which is accomplished by the signal reception, transcriptional/translational regulation, transport, storage and accumulation. The transcription factors (TFs) initiate the biosynthesis of alkaloids after appropriate cues. The present study recapitulates last decade understanding of the role of TFs in alkaloid biosynthesis. The present review discusses TF families, viz. AP2/ERF, bHLH, WRKY, MYB involved in the biosynthesis of various types of alkaloids. It also highlights the role of the jasmonic acid cascade and post-translational modifications of TF proteins. A thorough understanding of TFs will help us to decide a strategy to facilitate successful pathway manipulation and in vitro production.

Function of Cajal Bodies in Nuclear RNA Retention in A. thaliana Leaves Subjected to Hypoxia

Retention of RNA in the nucleus precisely regulates the time and rate of translation and controls transcriptional bursts that can generate profound variability in mRNA levels among identical cells in tissues. In this study, we investigated the function of Cajal bodies (CBs) in RNA retention in A. thaliana leaf nuclei during hypoxia stress was investigated. It was observed that in ncb-1 mutants with a complete absence of CBs, the accumulation of poly(A⁺) RNA in the leaf nuclei was lower than that in wt under stress. Moreover, unlike in root cells, CBs store less RNA, and RNA retention in the nuclei is much less intense. Our results reveal that the function of CBs in the accumulation of RNA in nuclei under stress depends on the plant organ. Additionally, in ncb-1, retention of introns of mRNA RPB1 (largest subunit of RNA polymerase II) mRNA was observed. However, this isoform is highly accumulated in the nucleus. It thus follows that intron retention in transcripts is more important than CBs for the accumulation of RNA in nuclei. Accumulated mRNAs with introns in the nucleus could escape transcript degradation by NMD (nonsense-mediated mRNA decay). From non-fully spliced mRNAs in ncb-1 nuclei, whose levels increase during hypoxia, introns are removed during reoxygenation. Then, the mRNA is transferred to the cytoplasm, and the RPB1 protein is translated. Despite the accumulation of isoforms in nuclei with retention of introns in reoxygenation, ncb-1 coped much worse with long hypoxia, and manifested faster yellowing and shrinkage of leaves.

Transcriptome revealing the dual regulatory mechanism of ethylene on the rhynchophylline and isorhynchophylline in Uncaria rhynchophylla

Rhynchophylline (RIN) and isorhynchophylline (IRN) are extracted from Uncaria rhynchophylla, which are used to treat Alzheimer's disease. However, the massive accumulation of RIN and IRN in U. rhynchophylla requires exogenous stimulation. Ethylene is a potential stimulant for RIN and IRN biosynthesis, but there is no study on the role of ethylene in RIN or IRN synthesis. This study investigated the regulation of ethylene in RIN and IRN biosynthesis in U. rhynchophylla. An increase in the content of RIN and IRN was observed that could be attributed to the release of ethylene from 18 mM ethephon, while ethylene released from 36 mM ethephon reduced the content of RIN and IRN. The transcriptome and weighted gene co-expression network analysis indicated the up-regulation of seven key enzyme genes related to the RIN/IRN biosynthesis pathway and starch/sucrose metabolism pathway favored RIN/IRN synthesis. In comparison, the down-regulation of these seven key enzyme genes contributed to the reduction of RIN/IRN. Moreover, the inhibition of photosynthesis is associated with a reduction in RIN/IRN. Photosynthesis was restrained owing to the down-regulation of Lhcb1 and Lhcb6 after 36 mM ethephon treatment and further prevented supply of primary metabolites (such as α-D-glucose) for RIN/IRN synthesis. However, uninterrupted photosynthesis ensured a normal supply of primary metabolites at 18 mM ethephon treatment. AP2/ERF1, bHLH1, and bHLH2 may positively regulate the RIN/IRN accumulation, while NAC1 may play a negative regulatory role. Our results construct the potential bidirectional model for ethylene regulation on RIN/IRN synthesis and provide novel insight into the ethylene-mediated regulation of the metabolism of terpenoid indole alkaloids.

Genome-wide survey of the GATA gene family in camptothecin-producing plant Ophiorrhiza pumila

Background

Ophiorrhiza pumila (Rubiaceae) is capable of producing camptothecin (CPT), one monoterpene indole alkaloid extensively employed in the treatment of multiple cancers. Transcription factors (TFs) GATA are a group of transcription regulators involved in plant development and metabolism, and show the feature of binding to the GATA motif within the promoters of target genes. However, GATA TFs have not been characterized in O. pumila.

Result

In this study, a total of 18 GATA genes classified into four subfamilies were identified, which randomly distributed on 11 chromosomes of O. pumila. Synteny analysis of GATA genes between O. pumila and other plant species such as Arabidopsis thaliana, Oryza sativa, Glycine max, Solanum lycopersicum, Vitis vinifera, and Catharanthus roseus genomes were analyzed. Tissue expression pattern revealed that OpGATA1 and OpGATA18 were found to be correlated with ASA, MK, CPR and GPPS, which were highly expressed in leaves. OpGATA7, showed high expression in roots as most of the CPT biosynthetic pathway genes did, suggesting that these OpGATAs may be potential candidates regulating CPT biosynthesis in O. pumila.

Conclusions

In this study, we systematically analyzed the OpGATA TFs, and provided insights into the involvement of OpGATA TFs from O. pumila in CPT biosynthesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08484-x.

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.

Genome-Wide Identification and Comparative Analysis of the Teosinte Branched 1/Cycloidea/Proliferating Cell Factors 1/2 Transcription Factors Related to Anti-cancer Drug Camptothecin Biosynthesis in Ophiorrhiza pumila

Ophiorrhiza pumila (O. pumila; Op) is a medicinal herbaceous plant, which can accumulate camptothecin (CPT). CPT and its derivatives are widely used as chemotherapeutic drugs for treating malignant tumors. Its biosynthesis pathway has been attracted significant attention. Teosinte branched 1/cycloidea/proliferating cell factors 1/2 (TCP) transcription factors (TFs) regulate a variety of physiological processes, while TCP TFs are involved in the regulation of CPT biosynthesis remain unclear. In this study, a systematic analysis of the TCP TFs family in O. pumila was performed. A total of 16 O. pumila TCP (OpTCP) genes were identified and categorized into two subgroups based on their phylogenetic relationships with those in Arabidopsis thaliana. Tissue-specific expression patterns revealed that nine OpTCP genes showed the highest expression levels in leaves, while the other seven OpTCPs showed a higher expression level in the stems. Co-expression, phylogeny analysis, and dual-luciferase (Dual-LUC) assay revealed that OpTCP15 potentially plays important role in CPT and its precursor biosynthesis. In addition, the subcellular localization experiment of candidate OpTCP genes showed that they are all localized in the nucleus. Our study lays a foundation for further functional characterization of the candidate OpTCP genes involved in CPT biosynthesis regulation and provides new strategies for increasing CPT production.

Strictosidine synthase, an indispensable enzyme involved in the biosynthesis of terpenoid indole and β-carboline alkaloids

Terpenoid indole (TIAs) and β-carboline alkaloids (BCAs), such as suppressant reserpine, vasodilatory yohimbine, and antimalarial quinine, are natural compounds derived from strictosidine. These compounds can exert powerful pharmacological effects but be obtained from limited source in nature. the whole biosynthetic pathway of TIAs and BCAs, The Pictet-Spengler reaction catalyzed by strictosidine synthase (STR; EC: 4.3.3.2) is the rate-limiting step. Therefore, it is necessary to investigate their biosynthesis pathways, especially the role of STR, and related findings will support the biosynthetic generation of natural and unnatural compounds. This review summarizes the latest studies concerning the function of STR in TIA and BCA biosynthesis, and illustrates the compounds derived from strictosidine. The substrate specificity of STR based on its structure is also summarized. Proteins that contain six-bladed four-stranded β-propeller folds in many organisms, other than plants, are listed. The presence of these folds may lead to similar functions among organisms. The expression of STR gene can greatly influence the production of many compounds. STR is mainly applied to product various valuable drugs in plant cell suspension culture and biosynthesis in other carriers.

Current advances of endophytes as a platform for production of anti-cancer drug camptothecin

Camptothecin (CPT), a well-known monoterpenoid indole alkaloid with broad-spectrum anti-cancer activity, is produced from plants and endophytes. In view of the limitations of plants as sources of camptothecin in productivity and efficiency, endophytes serve as the fast growth, high cost-effectiveness, good reproducibility, and feasible genetic manipulation, so they have the potential to meet the huge market demand of the pharmaceutical industry. In this review, we summarized the isolation, identification and fermentation of CPT-producing endophytes, as well as the biosynthesis, extraction and detection of camptothecin from endophytes. Among them, we put emphasis on increasing the production of camptothecin in endophytes through different strategies such as changing the proportion of carbon, nitrogen and phosphate source, adding the precursors, elicitors or adsorbent resin, utilizing co-culture fermentation or fermenter culture. However, cell subculture and metabolic reprogramming affect the expression of camptothecin biosynthetic genes in CPT-producing endophytes, which poses a challenge to the industrial production of camptothecin. Therefore, it will be useful to gain insights through the review of these researches and provide alternative approaches to develop economical, eco-friendly and reliable natural products.

Expression Profiles of Alkaloid-Related Genes across the Organs of Narrow-Leafed Lupin (Lupinus angustifolius L.) and in Response to Anthracnose Infection

The main restraint obstructing the wider adoption of lupins as protein crops is the presence of bitter and toxic quinolizidine alkaloids (QAs), whose contents might increase under exposure to stressful environmental conditions. A poor understanding of how QAs accumulate hinders the breeding of sweet varieties. Here, we characterize the expression profiles of QA-related genes, along with the alkaloid content, in various organs of sweet and bitter narrow-leafed lupin (NLL, Lupinus angustifolius L.). Special attention is paid to the RAP2-7 transcription factor, a candidate regulator of the QA pathway. We demonstrate the upregulation of RAP2-7 and other QA-related genes, across the aerial organs of a bitter cultivar and the significant correlations between their expression levels, thus supporting the role of RAP2-7 as an important regulatory gene in NLL. Moreover, we showed that the initial steps of QA synthesis might occur independently in all aerial plant organs sharing common regulatory mechanisms. Nonetheless, other regulatory steps might be involved in RAP2-7-triggered QA accumulation, given its expression pattern in leaves. Finally, the examination of QA-related gene expression in plants infected with Colletotrichum lupini evidenced no connection between QA synthesis and anthracnose resistance, in contrast to the important role of polyamines during plant–pathogen interactions.

Genome-wide identification and analysis of AP2/ERF transcription factors related to camptothecin biosynthesis in Camptotheca acuminata

Camptotheca acuminata produces camptothecin (CPT), a monoterpene indole alkaloid (MIA) that is widely used in the treatment of lung, colorectal, cervical, and ovarian cancers. Its biosynthesis pathway has attracted significant attention, but the regulation of CPT biosynthesis by the APETALA2/ethylene-responsive factor (AP2/ERF) transcription factors (TFs) remains unclear. In this study, a systematic analysis of the AP2/ERF TFs family in C. acuminata was performed, including phylogeny, gene structure, conserved motifs, and gene expression profiles in different tissues and organs (immature bark, cotyledons, young flower, immature fruit, mature fruit, mature leaf, roots, upper stem, and lower stem) of C. acuminata. A total of 198 AP2/ERF genes were identified and divided into five relatively conserved subfamilies, including AP2 (26 genes), DREB (61 genes), ERF (92 genes), RAV (18 genes), and Soloist (one gene). The combination of gene expression patterns in different C. acuminata tissues and organs, the phylogenetic tree, the co-expression analysis with biosynthetic genes, and the analysis of promoter sequences of key enzymes genes involved in CPT biosynthesis pathways revealed that eight AP2/ERF TFs in C. acuminata might be involved in CPT synthesis regulation, which exhibit relatively high expression levels in the upper stem or immature bark. Among these, four genes (CacAP2/ERF123, CacAP2/ERF125, CacAP2/ERF126, and CacAP2/ERF127) belong to the ERF-B2 subgroup; two genes (CacAP2/ERF149 and CacAP2/ERF152) belong to the ERF-B3 subgroup; and two more genes (CacAP2/ERF095 and CacAP2/ERF096) belong to the DREB-A6 subgroup. These results provide a foundation for future functional characterization of the AP2/ERF genes to enhance the biosynthesis of CPT compounds of C. acuminata.

Chromosome-level genome assembly of Ophiorrhiza pumila reveals the evolution of camptothecin biosynthesis

Plant genomes remain highly fragmented and are often characterized by hundreds to thousands of assembly gaps. Here, we report chromosome-level reference and phased genome assembly of Ophiorrhiza pumila, a camptothecin-producing medicinal plant, through an ordered multi-scaffolding and experimental validation approach. With 21 assembly gaps and a contig N50 of 18.49 Mb, Ophiorrhiza genome is one of the most complete plant genomes assembled to date. We also report 273 nitrogen-containing metabolites, including diverse monoterpene indole alkaloids (MIAs). A comparative genomics approach identifies strictosidine biogenesis as the origin of MIA evolution. The emergence of strictosidine biosynthesis-catalyzing enzymes precede downstream enzymes' evolution post γ whole-genome triplication, which occurred approximately 110 Mya in O. pumila, and before the whole-genome duplication in Camptotheca acuminata identified here. Combining comparative genome analysis, multi-omics analysis, and metabolic gene-cluster analysis, we propose a working model for MIA evolution, and a pangenome for MIA biosynthesis, which will help in establishing a sustainable supply of camptothecin.

ERF Gene Clusters: Working Together to Regulate Metabolism

Plants produce structurally diverse specialized metabolites, including bioactive alkaloids and terpenoids, in response to biotic and abiotic environmental stresses. The APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) family of transcription factors (TFs) play key roles in regulating biosynthesis of specialized metabolites. Increasing genomic and functional evidence shows that a subset of the ERF genes occurs in clusters on the chromosomes. These jasmonate-responsive ERF TF gene clusters control the biosynthesis of many important metabolites, from natural products, such as nicotine and steroidal glycoalkaloids (SGAs), to pharmaceuticals, such as artemisinin, vinblastine, and vincristine. Here, we review the function, regulation, and evolution of ERF clusters and highlight recent advances in understanding the distinct roles of clustered ERF genes and their possible application in metabolic engineering.

The transcription factor OpWRKY2 positively regulates the biosynthesis of the anticancer drug camptothecin in Ophiorrhiza pumila

The limited bioavailability of plant-derived natural products with anticancer activity poses major challenges to the pharmaceutical industry. An example of this is camptothecin, a monoterpene indole alkaloid with potent anticancer activity that is extracted at very low concentrations from woody plants. Recently, camptothecin biosynthesis has been shown to become biotechnologically amenable in hairy-root systems of the natural producer Ophiorrhiza pumila. Here, time-course expression and metabolite analyses were performed to identify novel transcriptional regulators of camptothecin biosynthesis in O. pumila. It is shown here that camptothecin production increased over cultivation time and that the expression pattern of the WRKY transcription factor encoding gene OpWRKY2 is closely correlated with camptothecin accumulation. Overexpression of OpWRKY2 led to a more than three-fold increase in camptothecin levels. Accordingly, silencing of OpWRKY2 correlated with decreased camptothecin levels in the plant. Further detailed molecular characterization by electrophoretic mobility shift, yeast one-hybrid and dual-luciferase assays showed that OpWRKY2 directly binds and activates the central camptothecin pathway gene OpTDC. Taken together, the results of this study demonstrate that OpWRKY2 acts as a direct positive regulator of camptothecin biosynthesis. As such, a feasible strategy for the over-accumulation of camptothecin in a biotechnologically amenable system is presented.

The transcription factor OpWRKY2 positively regulates the biosynthesis of the anticancer drug camptothecin in Ophiorrhiza pumila

The limited bioavailability of plant-derived natural products with anticancer activity poses major challenges to the pharmaceutical industry. An example of this is camptothecin, a monoterpene indole alkaloid with potent anticancer activity that is extracted at very low concentrations from woody plants. Recently, camptothecin biosynthesis has been shown to become biotechnologically amenable in hairy-root systems of the natural producer Ophiorrhiza pumila. Here, time-course expression and metabolite analyses were performed to identify novel transcriptional regulators of camptothecin biosynthesis in O. pumila. It is shown here that camptothecin production increased over cultivation time and that the expression pattern of the WRKY transcription factor encoding gene OpWRKY2 is closely correlated with camptothecin accumulation. Overexpression of OpWRKY2 led to a more than three-fold increase in camptothecin levels. Accordingly, silencing of OpWRKY2 correlated with decreased camptothecin levels in the plant. Further detailed molecular characterization by electrophoretic mobility shift, yeast one-hybrid and dual-luciferase assays showed that OpWRKY2 directly binds and activates the central camptothecin pathway gene OpTDC. Taken together, the results of this study demonstrate that OpWRKY2 acts as a direct positive regulator of camptothecin biosynthesis. As such, a feasible strategy for the over-accumulation of camptothecin in a biotechnologically amenable system is presented.

Solanum steroidal glycoalkaloids: structural diversity, biological activities, and biosynthesis

Chemical structures of typical Solanum steroidal glycoalkaloids from eggplant, tomato, and potato.

Genome-wide identification of AP2/ERF transcription factor-encoding genes in California poppy (Eschscholzia californica) and their expression profiles in response to methyl jasmonate

With respect to the biosynthesis of plant alkaloids, that of benzylisoquinoline alkaloids (BIAs) has been the most investigated at the molecular level. Previous investigations have shown that the biosynthesis of BIAs is comprehensively regulated by WRKY and bHLH transcription factors, while promoter analyses of biosynthesis enzyme-encoding genes have also implicated the involvement of members of the APETALA2/ethylene responsive factor (AP2/ERF) superfamily. To investigate the physiological roles of AP2/ERF transcription factors in BIA biosynthesis, 134 AP2/ERF genes were annotated using the draft genome sequence data of Eschscholzia californica (California poppy) together with transcriptomic data. Phylogenetic analysis revealed that these genes could be classified into 20 AP2, 5 RAV, 47 DREB, 60 ERF and 2 Soloist family members. Gene structure, conserved motif and orthologous analyses were also carried out. Gene expression profiling via RNA sequencing in response to methyl jasmonate (MeJA) indicated that approximately 20 EcAP2/ERF genes, including 10 group IX genes, were upregulated by MeJA, with an increase in the expression of the transcription factor-encoding gene EcbHLH1 and the biosynthesis enzyme-encoding genes Ec6OMT and EcCYP719A5. Further quantitative RT-PCR confirmed the MeJA responsiveness of the EcAP2/ERF genes, i.e., the increased expression of 9 group IX, 2 group X and 2 group III ERF subfamily genes. Transactivation activity of group IX EcAP2/ERFs was also confirmed by a luciferase reporter assay in conjunction with the promoters of the Ec6OMT and EcCYP719A5 genes. The physiological roles of AP2/ERF genes in BIA biosynthesis and their evolution in the regulation of alkaloid biosynthesis are discussed.

Biotechnological Exploration of Transformed Root Culture for Value-Added Products

Medicinal plants produce valuable secondary metabolites with anticancer, analgesic, anticholinergic or other activities, but low metabolite levels and limited available tissue restrict metabolite yields. Transformed root cultures, also called hairy roots, provide a feasible approach for producing valuable secondary metabolites. Various strategies have been used to enhance secondary metabolite production in hairy roots, including increasing substrate availability, regulating key biosynthetic genes, multigene engineering, combining genetic engineering and elicitation, using transcription factors (TFs), and introducing new genes. In this review, we focus on recent developments in hairy roots from medicinal plants, techniques to boost production of desired secondary metabolites, and the development of new technologies to study these metabolites. We also discuss recent trends, emerging applications, and future perspectives.

Metabolome and Transcriptome Analysis of Hexaploid Solidago canadensis Roots Reveals its Invasive Capacity Related to Polyploidy

Polyploid plants are more often invasive species than their diploid counterparts. As the invasiveness of a species is often linked to its production of allelopathic compounds, we hypothesize that differences in invasive ability between cytotypes may be due to their different ability to synthesize allelopathic metabolites. We test this using two cytotypes of Solidago canadensis as the model and use integrated metabolome and transcriptome data to resolve the question. Metabolome analysis identified 122 metabolites about flavonoids, phenylpropanoids and terpenoids, of which 57 were differentially accumulated between the two cytotypes. Transcriptome analysis showed that many differentially expressed genes (DEGs) were enriched in ‘biosynthesis of secondary metabolites’, ‘plant hormone signal transduction’, and ‘MAPK signaling’, covering most steps of plant allelopathic metabolite synthesis. Importantly, the differentially accumulated flavonoids, phenylpropanoids and terpenoids were closely correlated with related DEGs. Furthermore, 30 miRNAs were found to be negatively associated with putative targets, and they were thought to be involved in target gene expression regulation. These miRNAs probably play a vital role in the regulation of metabolite synthesis in hexaploid S. canadensis. The two cytotypes of S. canadensis differ in the allelopathic metabolite synthesis and this difference is associated with regulation of expression of a range of genes. These results suggest that changes in gene expression may underlying the increased invasive potential of the polyploidy.

Mutually Regulated AP2/ERF Gene Clusters Modulate Biosynthesis of Specialized Metabolites in Plants

APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) gene clusters regulate the biosynthesis of diverse specialized metabolites, including steroidal glycoalkaloids in tomato (Solanum lycopersicum) and potato (Solanum tuberosum), nicotine in tobacco (Nicotiana tabacum), and pharmaceutically valuable terpenoid indole alkaloids in Madagascar periwinkle (Catharanthus roseus). However, the regulatory relationships between individual AP2/ERF genes within the cluster remain unexplored. We uncovered intracluster regulation of the C. roseus AP2/ERF regulatory circuit, which consists of ORCA3, ORCA4, and ORCA5 ORCA3 and ORCA5 activate ORCA4 by directly binding to a GC-rich motif in the ORCA4 promoter. ORCA5 regulates its own expression through a positive autoregulatory loop and indirectly activates ORCA3 In determining the functional conservation of AP2/ERF clusters in other plant species, we found that GC-rich motifs are present in the promoters of analogous AP2/ERF clusters in tobacco, tomato, and potato. Intracluster regulation is evident within the tobacco NICOTINE2 (NIC2) ERF cluster. Moreover, overexpression of ORCA5 in tobacco and of NIC2 ERF189 in C. roseus hairy roots activates nicotine and terpenoid indole alkaloid pathway genes, respectively, suggesting that the AP2/ERFs are functionally equivalent and are likely to be interchangeable. Elucidation of the intracluster and mutual regulation of transcription factor gene clusters advances our understanding of the underlying molecular mechanism governing regulatory gene clusters in plants.

A cheminformatics approach to characterize metabolomes in stable-isotope-labeled organisms

We report a computational approach (implemented in MS-DIAL 3.0; http://prime.psc.riken.jp/) for metabolite structure characterization using fully ¹³C-labeled and non-labeled plants and LC-MS/MS. Our approach facilitates carbon number determination and metabolite classification for unknown molecules. Applying our method to 31 tissues from 12 plant species, we assigned 1,092 structures and 344 formulae to 3,604 carbon-determined metabolite ions, 69 of which were found to represent structures currently not listed in metabolome databases.

Transcriptome-derived investigation of biosynthesis of quinolizidine alkaloids in narrow-leafed lupin (Lupinus angustifolius L.) highlights candidate genes linked to iucundus locus

Unravelling the biosynthetic pathway of quinolizidine alkaloids (QAs), regarded as antinutritional compounds of narrow-leafed lupin (NLL) seeds, is fundamental to best exploit NLL as food or feed. We investigated 12 candidate genes connected to QA biosynthesis, selecting them by transcriptomic and genomic approaches, from the landscape of genes differentially expressed in leaves of the high- and low-alkaloid NLL accessions. Linkage analysis enabled the assessment of the location of the candidate genes in relation to iucundus, a major locus of unknown identity, that confers reduced QA content in seeds. The key finding was the identification of APETALA2/ethylene response transcription factor, RAP2-7, cosegregating with the iucundus locus and located within a region with highly significant QTLs that affect QA composition. We additionally identified a 4-hydroxy-tetrahydrodipicolinate synthase (DHDPS) gene involved in L-lysine biosynthesis as being closely linked to iucundus. The distributed location of other remaining candidates (including previously known QA genes) across different linkage groups, also indirectly supports the transcription factor as a possible regulator of lupin alkaloid biosynthesis. Our findings provide crucial insight into QA biosynthesis in NLL. Additionally, we evaluated and selected appropriate reference genes for qRT-PCRs to analyse the expression levels of QA genes in NLL.

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.

Analysis of long-non-coding RNAs associated with ethylene in tomato

Long-Non-Coding RNAs (LncRNAs) are a class of non-coding endogenous RNAs contributing to numerous biological processes. LeERF1 is a tomato ethylene response factor (ERF) near the end of the ethylene signal transduction pathway. To identify lncRNAs in tomato and elucidate their roles in ethylene signaling, deep sequencing was deployed in over-expression and repression LeERF1 transgenic and control tomato fruits. A total of 397 lncRNAs were identified, including 169 tomato lncRNAs that had not previously been identified. Among these, 12 were differentially expressed between the transgenic and control tomato fruits. Numerous lncRNA target genes were identified including many associated with ethylene signaling including auxin response factors and auxin-induced proteins, F-box proteins, ERFs and MADS-box proteins. In addition, two lncRNAs were found to be the precursor of three miRNAs and four lncRNAs could be targeted by five miRNAs. We propose a regulatory model highlighting the relationships between lncRNAs and their targets involved in ethylene signal transduction which establishes a foundation for addressing the role of LncRNAs in ethylene response.

[Development of Plant Metabolomics and Medicinal Plant Genomics]

 A variety of chemicals produced by plants, often referred to as 'phytochemicals', have been used as medicines, food, fuels and industrial raw materials. Recent advances in the study of genomics and metabolomics in plant science have accelerated our understanding of the mechanisms, regulation and evolution of the biosynthesis of specialized plant products. We can now address such questions as how the metabolomic diversity of plants is originated at the levels of genome, and how we should apply this knowledge to drug discovery, industry and agriculture. Our research group has focused on metabolomics-based functional genomics over the last 15 years and we have developed a new research area called 'Phytochemical Genomics'. In this review, the development of a research platform for plant metabolomics is discussed first, to provide a better understanding of the chemical diversity of plants. Then, representative applications of metabolomics to functional genomics in a model plant, Arabidopsis thaliana, are described. The extension of integrated multi-omics analyses to non-model specialized plants, e.g., medicinal plants, is presented, including the identification of novel genes, metabolites and networks for the biosynthesis of flavonoids, alkaloids, sulfur-containing metabolites and terpenoids. Further, functional genomics studies on a variety of medicinal plants is presented. I also discuss future trends in pharmacognosy and related sciences.

Transcriptome analysis in tissue sectors with contrasting crocins accumulation provides novel insights into apocarotenoid biosynthesis and regulation during chromoplast biogenesis

Crocins, the red soluble apocarotenoids of saffron, accumulate in the flowers of Crocus species in a developmental and tissue-specific manner. In Crocus sieberi, crocins accumulate in stigmas but also in a distinct yellow tepal sector, which we demonstrate contains chromoplast converted from amyloplasts. Secondary metabolites were analysed by LC-DAD-HRMS, revealing the progressive accumulation of crocetin and crocins in the yellow sector, which were also localized in situ by Raman microspectroscopy. To understand the underlying mechanisms of crocin biosynthesis, we sequenced the C. sieberi tepal transcriptome of two differentially pigmented sectors (yellow and white) at two developmental stages (6 and 8) by Illumina sequencing. A total of 154 million high-quality reads were generated and assembled into 248,099 transcripts. Differentially expressed gene analysis resulted in the identification of several potential candidate genes involved in crocin metabolism and regulation. The results provide a first profile of the molecular events related to the dynamics of crocetin and crocin accumulation during tepal development, and present new information concerning apocarotenoid biosynthesis regulators and their accumulation in Crocus. Further, reveals genes that were previously unknown to affect crocin formation, which could be used to improve crocin accumulation in Crocus plants and the commercial quality of saffron spice.

A Novel AP2/ERF Transcription Factor CR1 Regulates the Accumulation of Vindoline and Serpentine in Catharanthus roseus

As one type of the most important alkaloids in the world, terpenoid indole alkaloids (TIAs) show a wide range of pharmaceutical activities that are beneficial for clinical treatments. Catharanthus roseus produces approximately 130 identified TIAs and is considered to be a model plant to study TIA biosynthesis. In order to increase the production of high medical value metabolites whose yields are extremely low in C. roseus, genetic engineering combined with transcriptional regulation has been applied in recent years. By using bioinformatics which is based on RNA sequencing (RNA-seq) data from methyl jasmonate (MeJA)-treated C. roseus as well as phylogenetic analysis, the present work aims to screen candidate genes that may be involved in the regulation of TIA biosynthesis, resulting in a novel AP2/ERF transcription factor, CR1 (Catharanthus roseus 1). Subsequently, virus-induced gene silencing (VIGS) of CR1 was carried out to identify the involvement of CR1 in the accumulations of several TIAs and quantitative real-time PCR (qRT-PCR) was then applied to detect the expression levels of 7 genes in the related biosynthetic pathway in silenced plants. The results show that all the 7 genes were upregulated in CR1-silenced plants. Furthermore, metabolite analyses indicate that silencing CR1 could increase the accumulations of vindoline and serpentine in C. roseus. These results suggest a novel negative regulator which may be involved in the TIAs biosynthetic pathway.

Integrated omics analysis of specialized metabolism in medicinal plants

Medicinal plants are a rich source of highly diverse specialized metabolites with important pharmacological properties. Until recently, plant biologists were limited in their ability to explore the biosynthetic pathways of these metabolites, mainly due to the scarcity of plant genomics resources. However, recent advances in high-throughput large-scale analytical methods have enabled plant biologists to discover biosynthetic pathways for important plant-based medicinal metabolites. The reduced cost of generating omics datasets and the development of computational tools for their analysis and integration have led to the elucidation of biosynthetic pathways of several bioactive metabolites of plant origin. These discoveries have inspired synthetic biology approaches to develop microbial systems to produce bioactive metabolites originating from plants, an alternative sustainable source of medicinally important chemicals. Since the demand for medicinal compounds are increasing with the world's population, understanding the complete biosynthesis of specialized metabolites becomes important to identify or develop reliable sources in the future. Here, we review the contributions of major omics approaches and their integration to our understanding of the biosynthetic pathways of bioactive metabolites. We briefly discuss different approaches for integrating omics datasets to extract biologically relevant knowledge and the application of omics datasets in the construction and reconstruction of metabolic models.