The jasmonate-responsive AP2/ERF transcription factors AaERF1 and AaERF2 positively regulate artemisinin biosynthesis in Artemisia annua L

Multilayered regulation of secondary metabolism in medicinal plants

Medicinal plants represent a huge reservoir of secondary metabolites (SMs), substances with significant pharmaceutical and industrial potential. However, obtaining secondary metabolites remains a challenge due to their low-yield accumulation in medicinal plants; moreover, these secondary metabolites are produced through tightly coordinated pathways involving many spatiotemporally and environmentally regulated steps. The first regulatory layer involves a complex network of transcription factors; a second, more recently discovered layer of complexity in the regulation of SMs is epigenetic modification, such as DNA methylation, histone modification and small RNA-based mechanisms, which can jointly or separately influence secondary metabolites by regulating gene expression. Here, we summarize the findings in the fields of genetic and epigenetic regulation with a special emphasis on SMs in medicinal plants, providing a new perspective on the multiple layers of regulation of gene expression.

Transcriptome and Metabolomics Integrated Analysis Reveals MdMYB94 Associated with Esters Biosynthesis in Apple (Malus × domestica)

Volatile esters are major aromas contributing to the organoleptic quality of apple fruit. However, the molecular mechanisms underlying the regulation of volatile ester biosynthesis in apple remain elusive. This study investigated the volatile profiles and transcriptomes of 'Qinguan' (QG) apple fruit during development and/or postharvest storage. Although the constitution of volatiles varied widely between the peel and flesh, the volatile profiles of the peel and flesh of ripening QG fruit were dominated by volatile esters. WGCNA results suggested that 19 genes belonging to ester biosynthesis pathways and 11 hub transcription factor genes potentially participated in the biosynthesis and regulation of esters. To figure out key regulators of ester biosynthesis, correlation network analysis, dual-luciferase assays, and yeast one-hybrid assay were conducted and suggested that MdMYB94 trans-activated the MdAAT2 promoter and participated in the regulation of ester biosynthesis. This study provides a framework for understanding ester biosynthesis and regulation in apple.

RNA sequencing in Artemisia annua L explored the genetic and metabolic responses to hardly soluble aluminum phosphate treatment

Artemisia annua L. is a medicinal plant valued for its ability to produce artemisinin, a molecule used to treat malaria. Plant nutrients, especially phosphorus (P), can potentially influence plant biomass and secondary metabolite production. Our work aimed to explore the genetic and metabolic response of A. annua to hardly soluble aluminum phosphate (AlPO₄, AlP), using soluble monopotassium phosphate (KH₂PO₄, KP) as a control. Liquid chromatography-mass spectrometry (LC-MS) was used to analyze artemisinin. RNA sequencing, gene ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were applied to analyze the differentially expressed genes (DEGs) under poor P conditions. Results showed a significant reduction in plant growth parameters, such as plant height, stem diameter, number of leaves, leaf areas, and total biomass of A. annua. Conversely, LC-MS analysis revealed a significant increase in artemisinin concentration under the AlP compared to the KP. Transcriptome analysis revealed 762 differentially expressed genes (DEGs) between the AlP and the KP. GH3, SAUR, CRE1, and PYL, all involved in plant hormone signal transduction, showed differential expression. Furthermore, despite the downregulation of HMGR in the artemisinin biosynthesis pathway, the majority of genes (ACAT, FPS, CYP71AV1, and ALDH1) were upregulated, resulting in increased artemisinin accumulation in the AlP. In addition, 12 transcription factors, including GATA and MYB, were upregulated in response to AlP, confirming their importance in regulating artemisinin biosynthesis. Overall, our findings could contribute to a better understanding the parallel transcriptional regulation of plant hormone transduction and artemisinin biosynthesis in A. annua L. in response to hardly soluble phosphorus fertilizer.

JA-regulated AaGSW1-AaYABBY5/AaWRKY9 complex regulates artemisinin biosynthesis in Artemisia annua

Artemisinin, a sesquiterpene lactone from A. annua, is an essential therapeutic against malaria. YABBY family transcription factor; AaYABBY5 is an activator of AaCYP71AV1 (cytochrome P450-dependent hydroxylase) and AaDBR2 (double bond reductase 2); however, the protein-protein interactions of AaYABBY5, as well as the mechanism of its regulation, are not elucidated before. AaWRKY9 protein is a positive regulator of artemisinin biosynthesis that activates AaGSW1 (Glandular trichome specific WRKY1) and AaDBR2 (double bond reductase 2), respectively. In this study, YABBY-WRKY interactions are revealed to indirectly regulate artemisinin production. AaYABBY5 significantly increased the activity of the luciferase (LUC) gene fused to the promoter of AaGSW1. Towards the molecular basis of this regulation, AaYABBY5 interaction with AaWRKY9 protein was found. The combined effectors AaYABBY5 + AaWRKY9 showed synergistic effects toward the activities of AaGSW1, and AaDBR2 promoters, respectively. In AaYABBY5 over-expression plants, the expression of GSW1 was found significantly increase when compared to that of AaYABBY5 antisense or control plants. Secondly, AaGSW1 was seen as an upstream activator of AaYABBY5. Thirdly, it was found that AaJAZ8, a transcriptional repressor of jasmonates signaling, interacted with AaYABBY5 and attenuated its activity. Co-expression of AaYABBY5 and antiAaJAZ8 in A. annua increased the activity of AaYABBY5 towards artemisinin biosynthesis. For the first time, the current study provided the molecular basis of regulation of artemisinin biosynthesis through YABBY-WRKY interactions and its regulation through AaJAZ8. This knowledge provides AaYABBY5 overexpression plants as a powerful genetic resource for artemisinin biosynthesis.

Transcriptional regulatory network of high-value active ingredients in medicinal plants

High-value active ingredients in medicinal plants have attracted research attention because of their benefits for human health, such as the antimalarial artemisinin, anticardiovascular disease tanshinones, and anticancer Taxol and vinblastine. Here, we review how hormones and environmental factors promote the accumulation of active ingredients, thereby providing a strategy to produce high-value drugs at a low cost. Focusing on major hormone signaling events and environmental factors, we review the transcriptional regulatory network mediating biosynthesis of representative active ingredients. In this network, many transcription factors (TFs) simultaneously control multiple synthase genes; thus, understanding the molecular mechanisms affecting transcriptional regulation of active ingredients will be crucial to developing new breeding possibilities.

AaWIN1, an AP2/ERF protein, positively regulates glandular secretory trichome initiation in Artemisia annua

Exploring the genetic network of glandular trichomes and manipulating genes relevant to secondary metabolite biosynthesis are of great importance and value. Artemisinin, a key antimalarial drug ingredient, is synthesized and stored in glandular secretory trichomes (GSTs) in Artemisia annua. WIN/SHN proteins, a clade of AP2/ERF family, are known as regulators for cuticle biosynthesis. However, their function in glandular trichome development is less unknown. In this study, we identified a WIN/SHN gene from A. annua and named it as AaWIN1. AaWIN1 was predominantly expressed in buds, flowers and trichomes, and encoded a nuclear-localized protein. Overexpressing AaWIN1 in A. annua significantly increased the density of GST as well as the artemisinin content. Furthermore, AaGSW2 was reported to play an important role in promoting GST initiation, and the expression of AaGSW2 was induced in AaWIN1-overexpression lines. AaMIXTA1, a MYB protein positively regulating trichome initiation and cuticle biosynthesis, was confirmed to interact with AaWIN1. In addition, the ectopic expression of AaWIN1 resulted in slender and curled leaves, fewer trichomes, and rising expressions of cuticle biosynthesis genes in Arabidopsis thaliana. Taken together, based on phenotype observations, content measurements and gene expression detections, AaWIN1 was considered as a positive regulator for GST initiation in A. annua.

LcERF134 increases the production of monoterpenes by activating the terpene biosynthesis pathway in Litsea cubeba

Litsea cubeba, an aromatic species of the Lauraceae family, produces a diverse array of monoterpenes. The biosynthesis of monoterpenes is regulated by transcriptional factors (TFs), such as APETALA2/ethylene response factor (AP2/ERF). However, the regulatory mechanisms that control the AP2/ERF gene responsible for the biosynthesis of monoterpenes in L. cubeba have yet to be elucidated. Here, we identified an AP2/ERF gene, LcERF134, as an activator for the accumulation of citral and other monoterpenes. The expression level of LcERF134 was consistent with terpene synthase LcTPS42 in the pericarp. The transient overexpression of LcERF134 significantly increased monoterpene production in L. cubeba as well as the expression of rate-limiting genes involved in the monoterpene biosynthesis pathway. Furthermore, yeast one-hybrid, dual-luciferase and electrophoretic mobility shift assays demonstrated that LcERF134 activated the monoterpene biosynthesis pathway by directly binding to the GCC-box elements of the LcTPS42 and LcGPPS.SSU1 promoters. However, the overexpression of LcERF134 in tomatoes had no impact on the synthesis of monoterpenes, thus indicating that LcERF134 is a species-specific TF. Our research demonstrated that LcERF134 significantly increased the biosynthesis of monoterpenes by inducing the expression of LcTPS42 and LcGPPS.SSU1, thus offering insight into how to enhance the flavor of L. cubeba essential oil.

Signal Molecules Regulate the Synthesis of Secondary Metabolites in the Interaction between Endophytes and Medicinal Plants

Signaling molecules act as the links and bridges between endophytes and host plants. The recognition of endophytes and host plants, the regulation of host plant growth and development, and the synthesis of secondary metabolites are not separated by the participation of signaling molecules. In this review, we summarized the types and characteristics of signaling molecules in medicinal plants and the recent processes in intracellular conduction and multi-molecular crosstalk of signaling molecules during interactions between endophytic bacteria and medicinal plants. In addition, we overviewed the molecular mechanism of signals in medical metabolite accumulation and regulation. This work provides a reference for using endophytic bacteria and medicinal plants to synthesize pharmaceutical active ingredients in a bioreactor.

CmWRKY41 activates CmHMGR2 and CmFPPS2 to positively regulate sesquiterpenes synthesis in Chrysanthemum morifolium

Chrysanthemum morifolium is one of the most significant multipurpose crops with ornamental, medicinal, and edible value. Terpenoids, an essentials component of volatile oils, are abundant in chrysanthemum. However, the transcriptional regulation of terpenoid biosynthesis in chrysanthemums remains unclear. In the present investigation, we identified CmWRKY41, whose expression pattern is similar to that of terpenoid content in chrysanthemum floral scent, as a candidate gene that may promote terpenoid biosynthesis in chrysanthemum. Two structural genes 3-hydroxy-3-methylglutaryl-CoA reductase 2 (CmHMGR2) and farnesyl pyrophosphate synthase 2 (CmFPPS2), play key role in terpene biosynthesis in chrysanthemum. CmWRKY41 can directly bind to the promoters of CmHMGR2 or CmFPPS2 through GTGACA or CTGACG elements and activate its expression to promote sesquiterpene biosynthesis. In summary, these results indicate that CmWRKY41 targets CmHMGR2 and CmFPPS2 to positively regulate sesquiterpene biosynthesis in chrysanthemums. This study preliminarily revealed the molecular mechanism of terpenoid biosynthesis in chrysanthemum while enriching the secondary metabolism regulatory network.

RNA-Sequencing Reveals the Involvement of Sesquiterpene Biosynthesis Genes and Transcription Factors during an Early Response to Mechanical Wounding of Aquilaria sinensis

Plants respond to wounding by reprogramming the expression of genes involved in secondary metabolism. Aquilaria trees produce many bioactive secondary metabolites in response to wounding, but the regulatory mechanism of agarwood formation in the early response to mechanical wounding has remained unclear. To gain insights into the process of transcriptome changes and to determine the regulatory networks of Aquilaria sinensis to an early response (15 days) to mechanical wounding, we collected A. sinensis samples from the untreated (Asc1) and treated (Asf1) xylem tissues and performed RNA sequencing (RNA-seq). This generated 49,102,523 (Asc1) and 45,180,981 (Asf1) clean reads, which corresponded to 18,927 (Asc1) and 19,258 (Asf1) genes, respectively. A total of 1596 differentially expressed genes (DEGs) were detected in Asf1 vs. Asc1 (|log2 (fold change)| ≥ 1, Padj ≤ 0.05), of which 1088 were up-regulated and 508 genes were down-regulated. GO and KEGG enrichment analysis of DEGs showed that flavonoid biosynthesis, phenylpropanoid biosynthesis, and sesquiterpenoid and triterpenoid biosynthesis pathways might play important roles in wound-induced agarwood formation. Based on the transcription factor (TF)-gene regulatory network analysis, we inferred that the bHLH TF family could regulate all DEGs encoding for farnesyl diphosphate synthase, sesquiterpene synthase, and 1-deoxy-D-xylulose-5-phosphate synthase (DXS), which contribute to the biosynthesis and accumulation of agarwood sesquiterpenes. This study provides insight into the molecular mechanism regulating agarwood formation in A. sinensis, and will be helpful in selecting candidate genes for improving the yield and quality of agarwood.

MYC2: A Master Switch for Plant Physiological Processes and Specialized Metabolite Synthesis

The jasmonic acid (JA) signaling pathway plays important roles in plant defenses, development, and the synthesis of specialized metabolites synthesis. Transcription factor MYC2 is a major regulator of the JA signaling pathway and is involved in the regulation of plant physiological processes and specialized metabolite synthesis. Based on our understanding of the mechanism underlying the regulation of specialized metabolite synthesis in plants by the transcription factor MYC2, the use of synthetic biology approaches to design MYC2-driven chassis cells for the synthesis of specialized metabolites with high medicinal value, such as paclitaxel, vincristine, and artemisinin, seems to be a promising strategy. In this review, the regulatory role of MYC2 in JA signal transduction of plants to biotic and abiotic stresses, plant growth, development and specialized metabolite synthesis is described in detail, which will provide valuable reference for the use of MYC2 molecular switches to regulate plant specialized metabolite biosynthesis.

AabHLH113 integrates jasmonic acid and abscisic acid signaling to positively regulate artemisinin biosynthesis in Artemisia annua

Artemisinin, a sesquiterpene lactone isolated from Artemisia annua, is in huge market demand due to its efficient antimalarial action, especially after the COVID-19 pandemic. Many researchers have elucidated that phytohormones jasmonic acid (JA) and abscisic acid (ABA) positively regulate artemisinin biosynthesis via types of transcription factors (TFs). However, the crosstalk between JA and ABA in regulating artemisinin biosynthesis remains unclear. Here, we identified a novel ABA- and JA-induced bHLH TF, AabHLH113, which positively regulated artemisinin biosynthesis by directly binding to the promoters of artemisinin biosynthetic genes, DBR2 and ALDH1. The contents of artemisinin and dihydroartemisinic acid increased by 1.71- to 2.06-fold and 1.47- to 2.23-fold, respectively, in AabHLH1113 overexpressed A. annua, whereas they decreased by 14-36% and 26-53%, respectively, in RNAi-AabHLH113 plants. Furthermore, we demonstrated that AabZIP1 and AabHLH112, which, respectively, participate in ABA and JA signaling pathway to regulate artemisinin biosynthesis, directly bind to and activate the promoter of AabHLH113. Collectively, we revealed a complex network in which AabHLH113 plays a key interrelational role to integrate ABA- and JA-mediated regulation of artemisinin biosynthesis.

Making small molecules in plants: A chassis for synthetic biology-based production of plant natural products

Plant natural products have been extensively exploited in food, medicine, flavor, cosmetic, renewable fuel, and other industrial sectors. Synthetic biology has recently emerged as a promising means for the cost-effective and sustainable production of natural products. Compared with engineering microbes for the production of plant natural products, the potential of plants as chassis for producing these compounds is underestimated, largely due to challenges encountered in engineering plants. Knowledge in plant engineering is instrumental for enabling the effective and efficient production of valuable phytochemicals in plants, and also paves the way for a more sustainable future agriculture. In this manuscript, we briefly recap the biosynthesis of plant natural products, focusing primarily on industrially important terpenoids, alkaloids, and phenylpropanoids. We further summarize the plant hosts and strategies that have been used to engineer the production of natural products. The challenges and opportunities of using plant synthetic biology to achieve rapid and scalable production of high-value plant natural products are also discussed.

Biosynthesis of α-Bisabolol by Farnesyl Diphosphate Synthase and α-Bisabolol Synthase and Their Related Transcription Factors in Matricaria recutita L

The essential oil of German chamomile (Matricaria recutita L.) is widely used in food, cosmetics, and the pharmaceutical industry. α-Bisabolol is the main active substance in German chamomile. Farnesyl diphosphate synthase (FPS) and α-bisabolol synthase (BBS) are key enzymes related to the α-bisabolol biosynthesis pathway. However, little is known about the α-bisabolol biosynthesis pathway in German chamomile, especially the transcription factors (TFs) related to the regulation of α-bisabolol synthesis. In this study, we identified MrFPS and MrBBS and investigated their functions by prokaryotic expression and expression in hairy root cells of German chamomile. The results suggest that MrFPS is the key enzyme in the production of sesquiterpenoids, and MrBBS catalyzes the reaction that produces α-bisabolol. Subcellular localization analysis showed that both MrFPS and MrBBS proteins were located in the cytosol. The expression levels of both MrFPS and MrBBS were highest in the extension period of ray florets. Furthermore, we cloned and analyzed the promoters of MrFPS and MrBBS. A large number of cis-acting elements related to light responsiveness, hormone response elements, and cis-regulatory elements that serve as putative binding sites for specific TFs in response to various biotic and abiotic stresses were identified. We identified and studied TFs related to MrFPS and MrBBS, including WRKY, AP2, and MYB. Our findings reveal the biosynthesis and regulation of α-bisabolol in German chamomile and provide novel insights for the production of α-bisabolol using synthetic biology methods.

Genome-wide analysis of AP2/ERF superfamily in Isatis indigotica

OBJECTIVE

AP2/ERF (APETALA2/ethylene-responsive factor) superfamily is one of the largest gene families in plants and has been reported to participate in various biological processes, such as the regulation of biosynthesis of active lignan. However, few studies have investigated the genome-wide role of the AP2/ERF superfamily in Isatis indigotica. This study establishes a complete picture of the AP2/ERF superfamily in I. indigotica and contributes valuable information for further functional characterization of IiAP2/ERF genes and supports further metabolic engineering.

METHODS

To identify the IiAP2/ERF superfamily genes, the AP2/ERF sequences from Arabidopsis thaliana and Brassica rapa were used as query sequences in the basic local alignment search tool. Bioinformatic analyses were conducted to investigate the protein structure, motif composition, chromosome location, phylogenetic relationship, and interaction network of the IiAP2/ERF superfamily genes. The accuracy of omics data was verified by quantitative polymerase chain reaction and heatmap analyses.

RESULTS

One hundred and twenty-six putative IiAP2/ERF genes in total were identified from the I. indigotica genome database in this study. By sequence alignment and phylogenetic analysis, the IiAP2/ERF genes were classified into 5 groups including AP2, ERF, DREB (dehydration-responsive element-binding factor), Soloist and RAV (related to abscisic acid insensitive 3/viviparous 1) subfamilies. Among which, 122 members were unevenly distributed across seven chromosomes. Sequence alignment showed that I. indigotica and A. thaliana had 30 pairs of orthologous genes, and we constructed their interaction network. The comprehensive analysis of gene expression pattern in different tissues suggested that these genes may play a significant role in organ growth and development of I. indigotica. Members that may regulate lignan biosynthesis in roots were also preliminarily identified. Ribonucleic acid sequencing analysis revealed that the expression of 76 IiAP2/ERF genes were up- or down-regulated under salt or drought treatment, among which, 33 IiAP2/ERF genes were regulated by both stresses.

CONCLUSION

This study undertook a genome-wide characterization of the AP2/ERF superfamily in I. indigotica, providing valuable information for further functional characterization of IiAP2/ERF genes and discovery of genetic targets for metabolic engineering.

Identification of key genes involved in secondary metabolite biosynthesis in Digitalis purpurea

The medicinal plant Digitalis purpurea produces cardiac glycosides that are useful in the pharmaceutical industry. These bioactive compounds are in high demand due to ethnobotany's application to therapeutic procedures. Recent studies have investigated the role of integrative analysis of multi-omics data in understanding cellular metabolic status through systems metabolic engineering approach, as well as its application to genetically engineering metabolic pathways. In spite of numerous omics experiments, most molecular mechanisms involved in metabolic pathways biosynthesis in D. purpurea remain unclear. Using R Package Weighted Gene Co-expression Network Analysis, co-expression analysis was performed on the transcriptome and metabolome data. As a result of our study, we identified transcription factors, transcriptional regulators, protein kinases, transporters, non-coding RNAs, and hub genes that are involved in the production of secondary metabolites. Since jasmonates are involved in the biosynthesis of cardiac glycosides, the candidate genes for Scarecrow-Like Protein 14 (SCL14), Delta24-sterol reductase (DWF1), HYDRA1 (HYD1), and Jasmonate-ZIM domain3 (JAZ3) were validated under methyl jasmonate treatment (MeJA, 100 μM). Despite early induction of JAZ3, which affected downstream genes, it was dramatically suppressed after 48 hours. SCL14, which targets DWF1, and HYD1, which induces cholesterol and cardiac glycoside biosynthesis, were both promoted. The correlation between key genes and main metabolites and validation of expression patterns provide a unique insight into the biosynthesis mechanisms of cardiac glycosides in D. purpurea.

PpERF1b-like enhances lignin synthesis in pear (Pyrus pyrifolia) 'hard-end' fruit

The hard-end is a disorder of pear fruit, however, the mechanisms underlying its development remain unknown. In this study, we found that the hard-end fruit contained a higher transcript abundance level of ethylene-response factor 1b-like (PpERF1b-like) and released more ethylene compared to normal pear. In the ethephon treated normal fruit, flesh tissues accumulated more lignin together with elevated expression of PpERF1b-like. Overexpressing PpERF1b-like transiently in fruit and stably in callus increased lignin accumulation and the expression of lignin biosynthesis genes; the opposite results were observed in fruit showing repressed expression of PpERF1b-like. These results confirmed the role of PpERF1b-like in promoting hard-end formation through promoting lignin synthesis. This study provided valuable information for further clarifying the regulation of hard-end formation in pear.

Transcriptome profiling of transcription factors in Ganoderma lucidum in response to methyl jasmonate

Ganoderma lucidum is a traditional Chinese medicine and its major active ingredients are ganoderma triterpenoids (GTs). To screen for transcription factors (TFs) that involved in the biosynthetic pathway of GTs in G. lucidum, the chemical composition in mycelia, primordium and fruiting body were analyzed, and the transcriptomes of mycelia induced by methyl jasmonate (MeJA) were analyzed. In addition, the expression level data of MeJA-responsive TFs in mycelia, primordia and fruiting body were downloaded from the database, and the correlation analysis was carried out between their expression profiles and the content of total triterpenoids. The results showed that a total of 89 components were identified, and the content of total triterpenoids was the highest in primordium, followed by fruiting body and mycelia. There were 103 differentially expressed TFs that response to MeJA-induction including 95 upregulated and 8 downregulated genes. These TFs were classified into 22 families including C2H2 (15), TFII-related (12), HTH (9), fungal (8), bZIP (6), HMG (5), DADS (2), etc. Correlation analysis showed that the expression level of GL23559 (MADS), GL26472 (HTH), and GL31187 (HMG) showed a positive correlation with the GTs content, respectively. While the expression level of GL25628 (fungal) and GL26980 (PHD) showed a negative correlation with the GTs content, respectively. Furthermore, the over expression of the Glmhr1 gene (GL25628) in Pichia pastoris GS115 indicated that it might be a negative regulator of GT biosynthesis through decreasing the production of lanosterol. This study provided useful information for a better understanding of the regulation of TFs involved in GT biosynthesis and fungal growth in G. lucidum.

Identification and Characterization of Transcription Factors Involved in Geraniol Biosynthesis in Rosa chinensis

Fragrance is an important characteristic of rose flowers and is largely determined by the terpenes. Rose has a unique NUDX1 (NUDIX HYDROLASES 1)-dependent monoterpene geraniol biosynthesis pathway, but little is known about its transcriptional regulation. In this study, we characterized two China rose (Rosa chinensis) materials from the 'Old Blush' variety with contrasting aromas. We profiled the volatile metabolome of both materials, and the results revealed that geraniol was the main component that distinguishes the aroma of these two materials. We performed a comparative transcriptome analysis of the two rose materials, from which we identified the hydrolase RcNUDX1 as a key factor affecting geraniol content, as well as 17 transcription factor genes co-expressed with RcNUDX1. We also determined that the transcription factor RcWRKY70 binds to four W-box motifs in the promoter of RcNUDX1, repressing RcNUDX1 expression, based on yeast one-hybrid and transient dual-luciferase assays. These results provide important information concerning the transcriptional regulatory framework underlying the control of geraniol production in rose.

From Plant to Yeast-Advances in Biosynthesis of Artemisinin

Malaria is a life-threatening disease. Artemisinin-based combination therapy (ACT) is the preferred choice for malaria treatment recommended by the World Health Organization. At present, the main source of artemisinin is extracted from Artemisia annua; however, the artemisinin content in A. annua is only 0.1-1%, which cannot meet global demand. Meanwhile, the chemical synthesis of artemisinin has disadvantages such as complicated steps, high cost and low yield. Therefore, the application of the synthetic biology approach to produce artemisinin in vivo has magnificent prospects. In this review, the biosynthesis pathway of artemisinin was summarized. Then we discussed the advances in the heterologous biosynthesis of artemisinin using microorganisms (Escherichia coli and Saccharomyces cerevisiae) as chassis cells. With yeast as the cell factory, the production of artemisinin was transferred from plant to yeast. Through the optimization of the fermentation process, the yield of artemisinic acid reached 25 g/L, thereby producing the semi-synthesis of artemisinin. Moreover, we reviewed the genetic engineering in A. annua to improve the artemisinin content, which included overexpressing artemisinin biosynthesis pathway genes, blocking key genes in competitive pathways, and regulating the expression of transcription factors related to artemisinin biosynthesis. Finally, the research progress of artemisinin production in other plants (Nicotiana, Physcomitrella, etc.) was discussed. The current advances in artemisinin biosynthesis may help lay the foundation for the remarkable up-regulation of artemisinin production in A. annua through gene editing or molecular design breeding in the future.

Genome-wide transcriptome analysis reveals the regulatory network governing terpene trilactones biosynthesis in Ginkgo biloba

Ginkgo biloba L. is currently the only remaining gymnosperm of the Ginkgoaceae Ginkgo genus, and its history can be traced back to the Carboniferous 200 million years ago. Terpene trilactones (TTLs) are one of the main active ingredients in G. biloba, including ginkgolides and bilobalide. They have a good curative effect on cardiovascular and cerebrovascular diseases because of their special antagonistic effect on platelet-activating factors. Therefore, it is necessary to deeply mine genes related to TTLs and to analyze their transcriptional regulation mechanism, which will hold vitally important scientific and practical significance for quality improvement and regulation of G. biloba. In this study, we performed RNA-Seq on the root, stem, immature leaf, mature leaf, microstrobilus, ovulate strobilus, immature fruit and mature fruit of G. biloba. The TTL regulatory network of G. biloba in different organs was revealed by different transcriptomic analysis strategies. Weighted gene co-expression network analysis (WGCNA) revealed that the five modules were closely correlated with organs. The 12 transcription factors, 5 structural genes and 24 Cytochrome P450 (CYP450) were identified as candidate regulators for TTL accumulation by WGCNA and cytoscape visualization. Finally, 6 APETALA2/ethylene response factors, 2 CYP450s and bHLH were inferred to regulate the metabolism of TTLs by correlation analysis. This study is the comprehensive in authenticating transcription factors, structural genes and CYP450 involved in TTL biosynthesis, thereby providing molecular evidence for revealing the comprehensive regulatory network involved in TTL metabolism in G. biloba.

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.

An Integrative Transcriptional Network Revealed Spatial Molecular Interplay Underlying Alantolactone and Inulin Biosynthesis in Inula racemosa Hook f

Inula racemosa Hook. f. (Pushkarmula), a perennial Himalayan herb known for its aromatic and phytopharmaceutical attributes, is not yet explored at genomic/transcriptomic scale. In this study, efforts were made to unveil the global transcriptional atlas underlying organ-specific specialized metabolite biosynthesis by integrating RNA-Seq analysis of 433 million sequenced reads with the phytochemical analysis of leaf, stem, and root tissues. Overall, 7242 of 83,772 assembled nonredundant unigenes were identified exhibiting spatial expression in leaf (3761), root (2748), and stem (733). Subsequently, integration of the predicted transcriptional interactome network of 2541 unigenes (71,841 edges) with gene ontology and KEGG pathway enrichment analysis revealed isoprenoid, terpenoid, diterpenoid, and gibberellin biosynthesis with antimicrobial activities in root tissue. Interestingly, the root-specific expression of germacrene-mediated alantolactone biosynthesis (GAS, GAO, G8H, IPP, DMAP, and KAO) and antimicrobial activities (BZR1, DEFL, LTP) well-supported with both quantitative expression profiling and phytochemical accumulation of alantolactones (726.08 μg/10 mg) and isoalantolactones (988.59 μg/10 mg), which suggests “roots” as the site of alantolactone biosynthesis. A significant interaction of leaf-specific carbohydrate metabolism with root-specific inulin biosynthesis indicates source (leaf) to sink (root) regulation of inulin. Our findings comprehensively demonstrate the source-sink transcriptional regulation of alantolactone and inulin biosynthesis, which can be further extended for upscaling the targeted specialized metabolites. Nevertheless, the genomic resource created in this study can also be utilized for development of genome-wide functionally relevant molecular markers to expedite the breeding strategies for genetic improvement of I. racemosa.

Overexpression of PagERF072 from Poplar Improves Salt Tolerance

Extreme environments, especially drought and high salt conditions, seriously affect plant growth and development. Ethylene-responsive factor (ERF) transcription factors play an important role in salt stress response. In this study, a significantly upregulated ERF gene was identified in 84K (Populus alba × P. glandulosa), which was named PagERF072. PagERF072 was confirmed to be a nuclear-localized protein. The results of yeast two-hybrid (Y2H) assay showed that PagERF072 protein exhibited no self-activating activity, and yeast one-hybrid (Y1H) demonstrated that PagERF072 could specifically bind to GCC-box element. Under salt stress, the transgenic poplar lines overexpressing PagERF072 showed improved salt tolerance. The activities of peroxidase (POD), superoxide dismutase (SOD) and catalase (CAT) in transgenic poplars were significantly increased relative to those of wild-type (WT) plants, whereas malondialdehyde (MDA) content showed an opposite trend. In addition, reactive oxygen species (ROS) was significantly reduced, and the expression levels of POD- and SOD-related genes were significantly increased in transgenic poplars under salt stress compared with WT. All results indicate that overexpression of the PagERF072 gene can improve the salt tolerance of transgenic poplars.

AabHLH112, a bHLH transcription factor, positively regulates sesquiterpenes biosynthesis in Artemisia annua

The bHLH transcription factors play important roles in the regulation of plant growth, development, and secondary metabolism. β-Caryophyllene, epi-cedrol, and β-farnesene, three kinds of sesquiterpenes mainly found in plants, are widely used as spice in the food industry and biological pesticides in agricultural production. Furthermore, they also have a significant value in the pharmaceutical industry. However, there is currently a lack of knowledge on the function of bHLH family TFs in β-caryophyllene, epi-cedrol, and β-farnesene biosynthesis. Here, we found that AabHLH112 transcription factor had a novel function to positively regulate β-carophyllene, epi-cedrol, and β-farnesene biosynthesis in Artemisia annua. Exogenous MeJA enhanced the expression of AabHLH112 and genes of β-caryophyllene synthase (CPS), epi-cedrol synthase (ECS), and β-farnesene synthase (BFS), as well as sesquiterpenes content. Dual-LUC assay showed the activation of AaCPS, AaECS, and AaBFS promoters were enhanced by AabHLH112. Yeast one-hybrid assay showed AabHLH112 could bind to the G-box (CANNTG) cis-element in promoters of both AaCPS and AaECS. In addition, overexpression of AabHLH112 in A. annua significantly elevated the expression levels of AaCPS, AaECS, and AaBFS as well as the contents of β-caryophyllene, epi-cedrol, and β-farnesene, while suppressing AabHLH112 expression by RNAi reduced the expression of the three genes and the contents of the three sesquiterpenes. These results suggested that AabHLH112 is a positive regulator of β-caryophyllene, epi-cedrol, and β-farnesene biosynthesis in A. annua.

MADS-box gene AaSEP4 promotes artemisinin biosynthesis in Artemisia annua

The plant Artemisia annua is well known for its production of artemisinin, a sesquiterpene lactone that is an effective antimalarial compound. Although remarkable progress has been made toward understanding artemisinin biosynthesis, the effect of MADS-box family transcription factors on artemisinin biosynthesis is still poorly understood. In this study, we identified a MADS transcription factor, AaSEP4, that was predominantly expressed in trichome. AaSEP4 acts as a nuclear-localized transcriptional activator activating the expression of AaGSW1 (GLANDULAR TRICHOME-SPECIFIC WRKY1). Dual-luciferase and Yeast one-hybrid assays revealed that AaSEP4 directly bound to the CArG motif in the promoter region of AaGSW1. Overexpression of AaSEP4 in A. annua significantly induced the expression of AaGSW1 and four artemisinin biosynthesis genes, including amorpha-4,11-diene synthase (ADS), cytochrome P450 monooxygenase (CYP71AV1), double-bond reductase 2 (DBR2) and aldehyde dehydrogenase 1 (ALDH1). Furthermore, the results of high-performance liquid chromatography (HPLC) showed that the artemisinin content was significantly increased in the AaSEP4-overexpressed plants. In addition, RT-qPCR results showed that AaSEP4 was induced by methyl jasmonic acid (MeJA) treatment. Taken together, these results explicitly demonstrate that AaSEP4 is a positive regulator of artemisinin biosynthesis, which can be used in the development of high-artemisinin yielding A. annua varieties.

Variation on a theme: the structures and biosynthesis of specialized fatty acid natural products in plants

Plants are able to construct lineage-specific natural products from a wide array of their core metabolic pathways. Considerable progress has been made toward documenting and understanding, for example, phenylpropanoid natural products derived from phosphoenolpyruvate via the shikimate pathway, terpenoid compounds built using isopentyl pyrophosphate, and alkaloids generated by the extensive modification of amino acids. By comparison, natural products derived from fatty acids have received little attention, except for unusual fatty acids in seed oils and jasmonate-like oxylipins. However, scattered but numerous reports show that plants are able to generate many structurally diverse compounds from fatty acids, including some with highly elaborate and unique structural features that have novel bioproduct functionalities. Furthermore, although recent work has shed light on multiple new fatty acid natural product biosynthesis pathways and products in diverse plant species, these discoveries have not been reviewed. The aims of this work, therefore, are to (i) review and systematize our current knowledge of the structures and biosynthesis of fatty acid-derived natural products that are not seed oils or jasmonate-type oxylipins, specifically, polyacetylenic, very-long-chain, and aromatic fatty acid-derived natural products, and (ii) suggest priorities for future investigative steps that will bring our knowledge of fatty acid-derived natural products closer to the levels of knowledge that we have attained for other phytochemical classes.

Cold stress regulates accumulation of flavonoids and terpenoids in plants by phytohormone, transcription process, functional enzyme, and epigenetics

Plants make different defense mechanisms in response to different environmental stresses. One common way is to produce secondary metabolites. Temperature is the main environmental factor that regulates plant secondary metabolites, especially flavonoids and terpenoids. Stress caused by temperature decreasing to 4-10 °C is conducive to the accumulation of flavonoids and terpenoids. However, the accumulation mechanism under cold stress still lacks a systematic explanation. In this review, we summarize three aspects of cold stress promoting the accumulation of flavonoids and terpenoids in plants, that is, by affecting (1) the content of endogenous plant hormones, especially jasmonic acid and abscisic acid; (2) the expression level and activity of important transcription factors, such as bHLH and MYB families. This aspect also includes post-translational modification of transcription factors caused by cold stress; (3) key enzyme genes expression and activity in the biosynthesis pathway, in addition, the rate-limiting enzyme and glycosyltransferases genes are responsive to cold stress. The systematic understanding of cold stress regulates flavonoids, and terpenoids will contribute to the future research of genetic engineering breeding, metabolism regulation, glycosyltransferases mining, and plant synthetic biology.

Genome-wide Investigation and Expression Analysis of the AP2/ERF Family for Selection of Agarwood Related Genes in Aquilaria sinensis (Lour.) Gilg

Aquilaria sinensis is an important non-timber tree species for producing high-value agarwood, which is widely used as a traditional medicine and incense. Agarwood is the product of Aquilaria trees in response to injury and fungal infection. AP2/ERF transcription factors play important roles in plant stress responses and metabolite biosynthesis. In this study, 119 AsAP2/ERF genes were identified from the A. sinensis genome and divided into ERF, AP2, RAV and Soloist subfamilies. Their conserved motif, gene structure, chromosomal localization, and subcellular localization were characterized. A stress/defense-related ERF-associated amphiphilic repression (EAR) motif and an EDLL motif were identified. Moreover, 11 genes that were highly expressed in the agarwood layer in response to whole-tree agarwood induction technique (Agar-Wit) treatment were chosen, and their expression levels in response to MeJA, SA or salt treatment were further analyzed using qRT-PCR. Among the 11 genes, eight belonged to subgroup B-3. All 11 genes were significantly upregulated under salt treatment, while eight genes were significantly induced by both MeJA and SA. In addition, the gene clusters containing these upregulated genes on chromosomes were observed. The results obtained from this research not only provide useful information for understanding the functions of AP2/ERF genes in A. sinensis but also identify candidate genes and gene clusters to dissect their regulatory roles in agarwood formation for future research.

Exploration of the Mechanisms of Differential Indole Alkaloid Biosynthesis in Dedifferentiated and Cambial Meristematic Cells of Catharanthus roseus Using Transcriptome Sequencing

Catharanthus roseus produces terpenoid indole alkaloids (TIAs) of high medicinal importance. The current research focuses on finding an efficient production system such as cell suspension cultures for high TIA concentrations. Catharanthus roseus cambial meristematic cells (CMCs) offer multiple advantages over dedifferentiated cells (DDCs) regarding growth, homogeneity, and shear resistance. Our lab has established a CMC culture system induced by C. roseus cambium. We determined the concentrations of TIAs in CMCs and DDCs. CMCs produced significantly higher concentrations of total alkaloids, vindoline, vinblastine, catharanthine, and ajmalicine as compared to DDCs. We then performed Illumina HiSeq transcriptome sequencing of CMCs and DDCs and explored the differential transcriptomic signatures. Of the 96,004 unigenes, 9,564 were differentially expressed between the 2 cell suspension types. These differentially expressed genes (DEGs) were enriched in 137 KEGG pathways. Most importantly, genes from the indole alkaloid biosynthesis and the upstream pathways i.e., tryptophan metabolism, monoterpenoid biosynthesis, tropane, piperidine, and pyridine alkaloid biosynthesis, and terpenoid backbone biosynthesis showed differential transcriptomic signatures. Remarkably, the expression of genes associated with plant hormone biosynthesis, signaling, and MAPK signaling pathways was relatable to the different TIA concentrations in CMCs and DDCs. These results put forward multiple target genes, transcription factors, and regulators to develop a large-scale TIA production system using C. roseus CMCs.

Transcriptome analysis reveals regulation mechanism of methyl jasmonate-induced terpenes biosynthesis in Curcuma wenyujin

Curcuma wenyujin is the source plant of three traditional Chinese medicines, which have been widely used in clinical treatment over 1000 years. The content of terpenes, the major medicinal active ingredients, is relatively low in this plant. Studies have shown that MeJA can promote terpenes biosynthesis in plants. However, the mechanism underlying the effect of MeJA in C. wenyujin remains unclear. In this work, the transcriptome of C. wenyujin leaves with MeJA treatment was analyzed to elucidate the regulation mechanism of MeJA-mediated terpene biosynthesis. Based on the RNA-seq data, 7,246 unigenes were differentially expressed with MeJA treatment. Expression pattern clustering of DEGs revealed that unigenes, related to JA biosynthesis and signal transduction, responded to exogenous MeJA stimulation on the early stage and maintained throughout the process. Subsequently, unigenes related to terpene biosynthesis pathway showed a significant up-regulation with 6 h treatment. The analysis results suggested that MeJA induced the expression of JA biosynthesis genes (such as LOXs, AOSs, AOCs, OPRs, and MFPs) and JA signal transduction core genes (JAZs and MYCs) to activate JA signaling pathway. Meanwhile, downstream JA-responsive genes presented up-regulated expression levels such as AACT, HMGSs, HMGRs, DXSs, DXRs, MCTs, HDSs, and HDRs, thus promoting terpenes biosynthesis. The transcriptional expressions of these genes were validated by qRT-PCR. In addition, six CwTPS genes in response to MeJA were identified. With MeJA treatment, the expression levels of CwTPSs were increased as well as those of the transcription factors MYB, NAC, bZIP, WRKY, AP2/ERF, and HLH. These TFs might potentially regulate terpenes biosynthesis. These results provide insights for regulation mechanism of terpenes biosynthesis.

PatDREB Transcription Factor Activates Patchoulol Synthase Gene Promoter and Positively Regulates Jasmonate-Induced Patchoulol Biosynthesis

The production of patchoulol in the patchouli (Pogostemon cablin) plant determines its application value, as it is the principal active sesquiterpene of essential oil extracted from this plant. Here, the promoter of patchoulol synthase gene (PatPTS_pro) was isolated and found to be methyl jasmonate (MeJA)-induced. A nucleus-localized AP2/ERF transcription factor PatDREB was identified as a transcription activator binding to PatPTS_pro, regulating patchoulol biosynthesis through modulating the gene expression. PatDREB also interacts with jasmonate ZIM-domain 4 (JAZ4). Furthermore, PatDREB could physically interact with the MYB-related transcription factor PatSWC4 and synergistically facilitate patchoulol biosynthesis. However, the transcriptional activation activity of the PatDREB-PatSWC4 complex could be inhibited by PatJAZ4, and JA could reverse this interference. Overall, we demonstrated the positive roles of PatDREB and the PatDREB-PatSWC4 complex in regulating patchoulol production, which advance our understanding of the regulatory network of patchoulol biosynthesis.

The Promoter Analysis of VvPR1 Gene: A Candidate Gene Identified through Transcriptional Profiling of Methyl Jasmonate Treated Grapevine (Vitis vinifera L.)

Methyl jasmonate (MeJA) plays a vital role in plant disease resistance and also induces the expression of disease resistance genes in plants. In this study, a transcriptome analysis was performed on grapevine leaves after 12, 24 and 48 h of MeJA-100 μM treatment. A total of 1242 differentially expressed genes (DEGs) were identified from the transcriptome data, and the analysis of the DEGs showed that genes related to phytohormone signal transduction, jasmonic acid-mediated defense, Mitogen-activated protein kinase (MAPK), and flavonoid biosynthetic pathways were upregulated. As Pathogenesis-related gene 1 (PR1) is an important marker gene in plant defense also upregulated by MeJA treatment in RNA-seq data, the VvPR1 gene was selected for a promoter analysis with β-glucuronidase (GUS) through transient expression in tobacco leaves against abiotic stress. The results showed that the region from −1837 bp to −558 bp of the VvPR1 promoter is the key region in response to hormone and wound stress. In this study, we extended the available knowledge about induced defense by MeJA in a grapevine species that is susceptible to different diseases and identified the molecular mechanisms by which this defense might be mediated.

The Transcription Factor MdERF78 Is Involved in ALA-Induced Anthocyanin Accumulation in Apples

As a friendly plant growth regulator to the environment, 5-aminolevulinic acid (ALA) has been widely used in plant production, such as fruit coloration, stress resistance, and so on. Previous studies have identified some genes that have a function in the anthocyanin accumulation induced by ALA. However, the regulatory mechanism has not been well revealed. In the current study, we proposed that an ALA-responsive transcription factor, MdERF78, regulated anthocyanin accumulation. MdERF78, overexpressed in apple peels or calli, resulted in a significant increase of anthocyanins, while MdERF78 interference had an opposite trend. Furthermore, the anthocyanin accumulation induced by MdERF78 overexpression was enhanced by exogenous ALA treatment, suggesting that MdERF78 was involved in the ALA-induced anthocyanin accumulation. Yeast one-hybrid and dual luciferase reporter assays revealed that MdERF78 bound to the promoters of MdF3H and MdANS directly and activated their expressions. Additionally, MdERF78 interacted with MdMYB1 and enhanced the transcriptional activity of MdMYB1 to its target gene promoters. Based on these, it can be concluded that MdERF78 has a positive function in ALA-induced anthocyanin accumulation via the MdERF78-MdF3Hpro/MdANSpro and MdERF78-MdMYB1-MdDFRpro/MdUFGTpro/MdGSTF12pro regulatory network. These findings provide new insights into the regulatory mechanism of ALA-promoted anthocyanin accumulation.

Overexpression of a novel heat-inducible ethylene-responsive factor gene LlERF110 from Lilium longiflorum decreases thermotolerance

AP2/ERF (APETALA2/ethylene-responsive factor) family transcription factors are involved in various plant-specific processes, especially in plant development and response to abiotic stress. However, their roles in thermotolerance are still largely unknown. In the current study, we identified a heat-inducible ERF member LlERF110 from Lilium longiflorum that was rapidly induced by high temperature. Its protein was localized in the nucleus, and transcriptional activation activity was observed in yeast and plant cells. In addition, LlERF110 was able to bind to GCC- and CGG-elements, but not to DRE-elements. Overexpression of LlERF110 conferred delayed bolting and bushy phenotype, with decreased thermotolerance accompanied by a disrupted ROS (reactive oxygen species) homeostasis in transgenic plants. The accumulation of LlERF110 may activate certain repressors related to heat stress response (HSR) and indirectly damage the normal expression of heat stress (HS)-protective genes such as AtHSFA2, which consequently leads to reduced thermotolerance. Our results implied that LlERF110 might function as a heat-inducible gene but may hinder the establishment of thermotolerance.

Cloning and homologous characterization of geranylgeranyl pyrophosphate synthase (GGPPS) from Withania somnifera revealed alterations in metabolic flux towards gibberellic acid biosynthesis

Overexpression of a novel geranylgeranyl pyrophosphate synthase gene (WsGGPPS) in planta resulted in increased levels of gibberellic acid and decrease in withanolide content. Withania somnifera (L.) Dunal, the herb from family Solanaceae is one of the most treasured medicinal plant used in traditional medicinal systems owing to its unique stockpile of pharmaceutically active secondary metabolites. Phytochemical and pharmacological studies in this plant were well established, but the genes affecting the regulation of biosynthesis of major metabolites were not well elucidated. In this study cloning and functional characterization of a key enzyme in terpenoid biosynthetic pathway viz. geranylgeranyl pyrophosphate synthase (EC 2.5.1.29) gene from Withania somnifera was performed. The full length WsGGPPS gene contained 1,104 base pairs that encode a polypeptide of 365 amino acids. The quantitative expression analysis suggested that WsGGPPS transcripts were expressed maximally in flower tissues followed by berry tissues. The expression levels of WsGGPPS were found to be regulated by methyl jasmonate (MeJA) and salicylic acid (SA). Amino acid sequence alignment and phylogenetic studies suggested that WsGGPPS had close similarities with GGPPS of Solanum tuberosum and Solanum pennellii. The structural analysis provided basic information about three dimensional features and physicochemical parameters of WsGGPPS protein. Overexpression of WsGGPPS in planta for its functional characterization suggested that the WsGGPPS was involved in gibberellic acid biosynthesis.

Comparative transcriptomic analysis of genes in the triterpene saponin biosynthesis pathway in leaves and roots of Ardisia kteniophylla A. DC., a plant used in traditional Chinese medicine

Ardisia kteniophylla (Primulaceae) is highly valued in traditional medicine due to its production of the pharmacologically active secondary metabolites, especially triterpenoid saponins in its roots. Although A. kteniophylla is very important in traditional medicine, the genetic basis for its production of triterpenoid saponins remains largely unknown. Therefore, we sequenced transcriptomes of A. kteniophylla to identify putative genes involved in production of triterpenoid saponins in both leaves and roots, and we used the transcriptomes to compare expression levels of these genes between the two organ systems. The production of triterpenoid saponins in plants is usually induced through hormonal signaling on account of the presence of pests. Thus, we treated plants with the hormones salicylic acid (SA) and methyl jasmonate (MeJA) and used quantitative real‐time PCR (qRT‐PCR) to investigate expression levels of genes involved in triterpenoid saponin biosynthesis. In total, we obtained transcriptomes for leaf and root tissues representing 52,454 unigenes. Compared with the leaf transcriptome, we found that 6092 unigenes were upregulated in the root, especially enzymes involved in the direct synthesis of triterpenoid saponins, while 6001 genes appeared downregulated, including those involved in precursory steps in the triterpenoid saponin biosynthesis pathway. Our results from qRT‐PCR indicate that genes within the upstream parts of the triterpenoid saponin biosynthesis pathway may be upregulated under exposure to the applied hormones, but downstream genes are downregulated. This suggests possible conflicting effects of SA and MeJA in promoting the production of secondary metabolites on the one hand, and, on the other, limiting plant growth processes to devote energy to combating pests. We also performed an analysis of transcription factors (TFs) and found 997 unique transcripts belonging to 16 TF families. Our data may help to facilitate future work on triterpene saponins biosynthesis in A. kteniophylla with potential pharmacological and molecular breeding applications.

Integrated metabolite profiling and transcriptome analysis reveals tissue-specific regulation of terpenoid biosynthesis in Artemisia argyi

Artemisia argyi L. is a widely distributed medicinal plant in China. The major bioactive substances of essential oils extracted from leaves are terpenoids. Although many researches have studied the pharmacological effects of the essential oils, the tissue-specific accumulation of terpenoid biosynthesis and the regulatory networks in A. argyi are poorly understood. This study conducted an integrated metabolomic and transcriptomic analysis of roots, stems, and leaves to investigate the tissue-specific regulatory network of terpenoid biosynthesis in A. argyi. We identified 77 unigenes putatively involved in terpenoid backbone biosynthesis. Three rate-determining enzyme genes (DXS, DXR, and HDR) of the methylerythritol phosphate pathway were predominantly expressed in leaves, and strongly co-expressed with eight transcription factors (2 MYBs, 4 WRKYs, and 2 AP2s). An metabolite-transcript correlation analysis revealed 26 putative cytochrome P450s related to terpenoid metabolism in leaves. These results provide a foundation for the future metabolic engineering of useful terpenoids in A. argyi.

The jasmonate-induced bHLH gene SlJIG functions in terpene biosynthesis and resistance to insects and fungus

Jasmonic acid (JA) is a key regulator of plant defense responses. Although the transcription factor MYC2, the master regulator of the JA signaling pathway, orchestrates a hierarchical transcriptional cascade that regulates the JA responses, only a few transcriptional regulators involved in this cascade have been described. Here, we identified the basic helix-loop-helix (bHLH) transcription factor gene in tomato (Solanum lycopersicum), METHYL JASMONATE (MeJA)-INDUCED GENE (SlJIG), the expression of which was strongly induced by MeJA treatment. Genetic and molecular biology experiments revealed that SlJIG is a direct target of MYC2. SlJIG knockout plants generated by gene editing had lower terpene contents than the wild type from the lower expression of TERPENE SYNTHASE (TPS) genes, rendering them more appealing to cotton bollworm (Helicoverpa armigera). Moreover, SlJIG knockouts exhibited weaker JA-mediated induction of TPSs, suggesting that SlJIG may participate in JA-induced terpene biosynthesis. Knocking out SlJIG also resulted in attenuated expression of JA-responsive defense genes, which may contribute to the observed lower resistance to cotton bollworm and to the fungus Botrytis cinerea. We conclude that SlJIG is a direct target of MYC2, forms a MYC2-SlJIG module, and functions in terpene biosynthesis and resistance against cotton bollworm and B. cinerea.

Comparative methylomics and chromatin accessibility analysis in uncovers regulation of genic transcription and mechanisms of key floral scent production

Linalool and ionone are two important aromatic components in sweet osmanthus petals, and the regulatory mechanisms that produce these two components remain unclear. In this study, we employed whole-genome methylation sequencing and ATAC-seq technology to analyze the genomic DNA methylation status and chromatin accessibility of the sweet osmanthus cultivars 'Zaohuang' and 'Chenghong Dangui'. Results showed that the promoter region of TPS2, a key gene in the linalool synthesis pathway, was less methylated in 'Chenghong Dangui' than in 'Zaohuang'. The chromatin was more accessible in 'Chenghong Dangui' than in 'Zaohuang', which resulted in a much stronger expression of this gene in 'Chenghong Dangui' than in 'Zaohuang'. This eventually led to a high quantity of linalool and its oxides in the petals of 'Chenghong Dangui', but there were lower levels present in the petals of 'Zaohuang'. These results suggest that DNA methylation and chromatin accessibility play major roles in linalool synthesis in sweet osmanthus. The methylation level of the promoter region of CCD4, a key gene for ionone synthesis, was higher in 'Zaohuang' than in 'Chenghong Dangui'. The chromatin accessibility was lower in 'Zaohuang' than in 'Chenghong Dangui', although the expression of this gene was significantly higher in 'Zaohuang' than in 'Chenghong Dangui'. ChIP-seq analysis and a series of experiments showed that the differential expression of CCD4 and CCD1 in the two cultivars may predominantly be the result of regulation by ERF2 and other transcription factors. However, a 183-bp deletion involving the CCD4 promoter region in 'Chenghong Dangui' may be the main reason for the low expression of this gene in its petals. This study provides an important theoretical basis for improving selective breeding of key floral fragrance components in sweet osmanthus.

LcERF19, an AP2/ERF transcription factor from Litsea cubeba, positively regulates geranial and neral biosynthesis

The APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) transcription factors (TFs) are involved in the regulation of specialized terpenoid biosynthesis. However, the AP2/ERF TFs in Litsea cubeba have not been characterized and their role in the biosynthesis of terpenoids is unknown. Here, 174 LcAP2/ERF TFs were identified in L. cubeba and categorized into four subfamilies: 27 AP2, 7 RAV, 1 Soloist, and 139 ERF. Transcriptomic and qRT-PCR assays both showed that the expression levels of LcERF19 were similar to that of terpene synthase LcTPS42 in the pericarp, which is related to the synthesis of geranial and neral in L. cubeba. LcERF19 was further shown to encode a nuclear-localized protein and its expression was strongly induced by jasmonate. Yeast one-hybrid and dual-luciferase assays showed that LcERF19 associated with GCC box elements of the LcTPS42 promoter and promoted its activity. Transient overexpression of LcERF19 in L. cubeba and overexpression of LcERF19 in tomato resulted in a significant increase in geranial and neral. Our findings show that LcERF19 enhances geranial and neral biosynthesis through activation of LcTPS42 expression, which provides a strategy to improve the flavor of tomato and other fruits.

Environmental and Genetic Factors Involved in Plant Protection-Associated Secondary Metabolite Biosynthesis Pathways

Plant specialized metabolites (PSMs) play essential roles in the adaptation to harsh environments and function in plant defense responses. PSMs act as key components of defense-related signaling pathways and trigger the extensive expression of defense-related genes. In addition, PSMs serve as antioxidants, participating in the scavenging of rapidly rising reactive oxygen species, and as chelators, participating in the chelation of toxins under stress conditions. PSMs include nitrogen-containing chemical compounds, terpenoids/isoprenoids, and phenolics. Each category of secondary metabolites has a specific biosynthetic pathway, including precursors, intermediates, and end products. The basic biosynthetic pathways of representative PSMs are summarized, providing potential target enzymes of stress-mediated regulation and responses. Multiple metabolic pathways share the same origin, and the common enzymes are frequently to be the targets of metabolic regulation. Most biosynthetic pathways are controlled by different environmental and genetic factors. Here, we summarized the effects of environmental factors, including abiotic and biotic stresses, on PSM biosynthesis in various plants. We also discuss the positive and negative transcription factors involved in various PSM biosynthetic pathways. The potential target genes of the stress-related transcription factors were also summarized. We further found that the downstream targets of these Transcription factors (TFs) are frequently enriched in the synthesis pathway of precursors, suggesting an effective role of precursors in enhancing of terminal products. The present review provides valuable insights regarding screening targets and regulators involved in PSM-mediated plant protection in non-model plants.

Linalool synthesis related PpTPS1 and PpTPS3 are activated by transcription factor PpERF61 whose expression is associated with DNA methylation during peach fruit ripening

The monoterpene linalool is a major contributor to flavor of multiple fruit species. Although great progress has been made in identifying genes related to linalool formation, transcriptional regulation for the pathway remains largely unknown. As a super transcription factor family, roles of AP2/ERF in regulating linalool production have not been elucidated. Peach linalool is catalyzed by terpene synthases PpTPS1 and PpTPS3. Here, we observed that expression of PpERF61 correlated with these two PpTPSs during fruit ripening by transcriptome co-expression analysis. Dual-luciferase assay and EMSA results indicated that PpERF61 activated the PpTPS1 and PpTPS3 transcription by binding to the DRE/CRT motif in their promoters. Transient overexpressing PpERF61 in peach fruit significantly increased PpTPS1 and PpTPS3 expression and linalool content. Further study revealed significant correlation between PpERF61 transcripts and linalool contents across 30 peach cultivars. Besides transcriptional regulation, accumulated linalool was associated with DNA demethylation of PpERF61 during peach fruit ripening. In addition, interactions between PpERF61 and PpbHLH1 were evaluated, indicating these two transcription factors were associated with linalool production during peach fruit ripening. Overall, our results revealed a new insight into the regulation of linalool synthesis in fruit.

Metabolome and transcriptome analyses identify the plant immunity systems that facilitate sesquiterpene and lignan biosynthesis in Syringa pinnatifolia Hemsl

BACKGROUND

Syringa pinnatifolia Hemsl. is a shrub belonging to the Oleaceae family. The peeled woody stems and roots of S. pinnatifolia are used in Chinese traditional medicine. This plant has been used for centuries, and modern pharmacological research has revealed its medicinal value. However, the wild populations of S. pinnatifolia have been decreasing, and it has been listed as an endangered plant in China. To elucidate the molecular mechanism leading to the synthesis of the major components of S. pinnatifolia for its further development and sustainable use, this study compared peeled stems and twigs at the metabolic and molecular levels.

RESULTS

Peeled stems with the purple substance visible (SSP) and peeled twigs without the purple substance (TSP) were compared at different levels. Microscopic observation showed resin-like fillers in SSP and wood fiber cell walls approximately 1.0 μm thicker than those in TSP (wood fiber cell thickness approximately 2.7 μm). In addition, 104 volatile organic compounds and 870 non-volatile metabolites were detected in the non-targeted and widely-targeted metabolome analyses, respectively. Among the 76 differentially accumulated metabolites (DAMs) detected, 62 were up-accumulated in SSP. Most of these DAMs were terpenes, of which 90% were identified as sesquiterpenes in the volatile organic compound analysis. In the analysis of the non-volatile metabolites, 21 differentially accumulated lignans were identified, of which 18, including five subtypes, were accumulated in SSP. RNA sequencing revealed 4,421 upregulated differentially expressed genes (DEGs) and 5,522 downregulated DEGs in SSP compared with TSP, as well as 33,452 genes that were not differentially expressed. Analysis of the DEGs suggested that sesquiterpenes and lignans were mostly biosynthesized via the mevalonate and phenylpropanoid pathways, respectively. Additionally, in SSP, the enriched Gene Ontology terms included response to biotic stimulus and defense response, while the enriched Kyoto Encyclopedia of Genes and Genomes pathways included plant-pathogen interaction and many other pathways related to plant immunity.

CONCLUSIONS

This study provides metabolome and transcriptome information for S. pinnatifolia, suggesting that biotic stimuli, including pathogens, are potential and valuable approaches to promoting the biosynthesis of the metabolites linked to the medicinal properties of this plant.

The Diesel Tree Sindora glabra Genome Provides Insights Into the Evolution of Oleoresin Biosynthesis

Sindora glabra is an economically important tree that produces abundant oleoresin in the trunk. Here, we present a high-quality chromosome-scale assembly of S. glabra genome by combining Illumina HiSeq, Pacific Biosciences sequencing, and Hi-C technologies. The size of S. glabra genome was 1.11 Gb, with a contig N50 of 1.27 Mb and 31,944 predicted genes. This is the first sequenced genome of the subfamily Caesalpinioideae. As a sister taxon to Papilionoideae, S. glabra underwent an ancient genome triplication shared by core eudicots and further whole-genome duplication shared by early-legume in the last 73.3 million years. S. glabra harbors specific genes and expanded genes largely involved in stress responses and biosynthesis of secondary metabolites. Moreover, 59 terpene backbone biosynthesis genes and 64 terpene synthase genes were identified, which together with co-expressed transcription factors could contribute to the diversity and specificity of terpene compounds and high terpene content in S. glabra stem. In addition, 63 disease resistance NBS-LRR genes were found to be unique in S. glabra genome and their expression levels were correlated with the accumulation of terpene profiles, suggesting potential defense function of terpenes in S. glabra. These together provide new resources for understanding genome evolution and oleoresin production.

The Multifaceted Roles of MYC2 in Plants: Toward Transcriptional Reprogramming and Stress Tolerance by Jasmonate Signaling

Environmental stress is one of the major restrictions on plant development and foodstuff production. The adaptive response in plants largely occurs through an intricate signaling system, which is crucial for regulating the stress-responsive genes. Myelocytomatosis (MYC) transcription factors are the fundamental regulators of the jasmonate (JA) signaling branch that participates in plant development and multiple stresses. By binding to the cis-acting elements of a large number of stress-responsive genes, JA-responsive transcription factors activate the stress-resistant defense genes. The mechanism of stress responses concerns myriad regulatory processes at the physiological and molecular levels. Discovering stress-related regulatory factors is of great value in disclosing the response mechanisms of plants to biotic or abiotic stress, which could guide the genetic improvement of plant resistance. This review summarizes recent researches in various aspects of MYC2-mediated JA signaling and emphasizes MYC2 involvement in plant growth and stress response.

Changes in Polar Metabolites Content during Natural and Methyl-Jasmonate-Promoted Senescence of Ginkgo biloba Leaves

The present study clarified changes in the contents of polar metabolites (amino acids, organic acids, saccharides, cyclitols, and phosphoric acid) in leaf senescence in Ginkgo biloba with or without the application of methyl jasmonate (JA-Me) in comparison with those in naturally senescent leaf blades and petioles. The contents of most amino acids and citric and malic acids were significantly higher in abaxially, and that of myo-inositol was lower in abaxially JA-Me-treated leaves than in adaxially JA-Me-treated and naturally senescent leaves. The levels of succinic and fumaric acids in leaves treated adaxially substantially high, but not in naturally senescent leaves. In contrast, sucrose, glucose, and fructose contents were much lower in leaf blades and petioles treated abaxially with JA-Me than those treated adaxially. The levels of these saccharides were also lower compared with those in naturally senescent leaves. Shikimic acid and quinic acid were present at high levels in leaf blades and petioles of G. biloba. In leaves naturally senescent, their levels were higher compared to green leaves. The shikimic acid content was also higher in the organs of naturally yellow leaves than in those treated with JA-Me. These results strongly suggest that JA-Me applied abaxially significantly enhanced processes of primary metabolism during senescence of G. biloba compared with those applied adaxially. The changes in polar metabolites in relation to natural senescence were also discussed.

Genome-wide investigation of the AP2/ERF gene family in ginger: evolution and expression profiling during development and abiotic stresses

BACKGROUND

AP2/ERF transcription factors (TFs) constitute one of the largest TF families in plants, which play crucial roles in plant metabolism, growth, and development as well as biotic and abiotic stresses responses. Although the AP2/ERF family has been thoroughly identified in many plant species and several AP2/ERF TFs have been functionally characterized, little is known about this family in ginger (Zingiber officinale Roscoe), an important affinal drug and diet vegetable. Recent completion of the ginger genome sequencing provides an opportunity to investigate the expression profiles of AP2/ERF genes in ginger on a genome-wide basis.

RESULTS

A total of 163 AP2/ERF genes were obtained in the Z.officinale genome and renamed according to the chromosomal distribution of the ZoAP2/ERF genes. Phylogenetic analysis divided them into three subfamilies, of which 35 belonged to the AP2 subfamily, 120 to ERF, three to RAV, and five to Sololist, respectively, which is in accordance with the number of conserved domains and gene structure analysis. A total of 10 motifs were detected in ZoAP2/ERF genes, and some of the unique motifs were found to be important for the function of ZoAP2/ERF genes. The chromosomal localization, gene structure, and conserved protein motif analyses, as well as the characterization of gene duplication events provided deep insight into the evolutionary features of these ZoAP2/ERF genes. The expression profiles derived from the RNA-seq data and quantitative reserve transcription (qRT-PCR) analysis of ZoAP2/ERFs during development and responses to abiotic stresses were investigated in ginger.

CONCLUSION

A comprehensive analysis of the AP2/ERF gene expression patterns in various tissues by RNA-seq and qRT-PCR showed that they played an important role in the growth and development of ginger, and genes that might regulate rhizome and flower development were preliminary identified. In additionally, the ZoAP2/ERF family genes that responded to abiotic stresses were also identified. This study is the first time to identify the ZoAP2/ERF family, which contributes to research on evolutionary characteristics and better understanding the molecular basis for development and abiotic stress response, as well as further functional characterization of ZoAP2/ERF genes with an aim of ginger crop improvement.

AaHY5 ChIP-seq based on transient expression system reveals the role of AaWRKY14 in artemisinin biosynthetic gene regulation

ChIP-seq (Chromatin immunoprecipitation with sequencing) is the gold standard for determining genome-wide in vivo transcription factor binding sites, the first step for targets prediction and network construction. For non-model plants, it is challenging to perform ChIP-seq due to the difficulty in generating stable transgenic plants. AaHY5 is a positive regulator in artemisinin biosynthesis, whose detailed mode of action remains elusive. Here, we established a protoplast transformation procedure for Artemisia annua by optimizing different conditions in protoplast isolation and transfection. We then performed AaHY5 ChIP-seq based on the established transient expression system. Combining RNA-seq data for various tissues, we identified four transcription factors (one MYB and three WRKY family members) in AaHY5 targets that potentially regulated artemisinin biosynthesis. The three WRKY transcription factors could be induced by light and the overexpression of AaHY5 and upregulate two artemisinin biosynthetic genes, ADS and CYP71AV1. Furthermore, AaWRKY14 showed transcriptional activation activity on artemisinin biosynthetic gene CYP71AV1. Together, AaWRKY14 was identified as a potential transcription factor linking AaHY5 and the artemisinin biosynthetic gene regulation.