Identification of a WRKY protein as a transcriptional regulator of benzylisoquinoline alkaloid biosynthesis in Coptis japonica

Transcription Factor CmWRKY13 Regulates Cucurbitacin B Biosynthesis Leading to Bitterness in Oriental Melon Fruit (Cucumis melo var. Makuwa Makino)

Bitterness, caused by cucurbitacin B (CuB), is one of the important traits that affects melon fruit quality and consumer acceptance. Therefore, the detailed mechanism behind the regulation of CuB biosynthesis on melon fruit needs to be further explored. This study investigated CuB content and transcriptomes of "YMR" melon fruit treated by 5 and 20 mg L-1 CPPU. The content of CuB reaches its peak in 5 days and then decreases. WGCNA identified the WRKY transcription factor (TF), CmWRKY13, coexpressed with CuB biosynthetic genes (Cm180, Cm170, Cm160, and CmACT). Yeast one-hybrid, dual-luciferase, and transient gene expression assays were conducted and suggested that the nucleus-localized CmWRKY13 transactivated the promoters of CuB biosynthetic genes and participated in the regulation of CuB biosynthesis. Furthermore, CmWRKY13 could interact with CmBt, the fruit bitterness-specific TF, which synergistically activated CuB biosynthetic gene expression. These findings provide a novel mechanistic insight for CuB biosynthesis and regulation in melon cultivation.

Transcriptome analysis of the genes and regulators involving in vitamin E biosynthesis in Elaeagnus mollis diels

Elaeagnus mollis is an important newly developing woody oil plant species and the vitamin E (VitE) content in its kernel oil is relatively high. In the present study, the VitE component content and functional genes involving in VitE biosynthesis in E. mollis kernel at different developmental stage were investigated. The VitE content increased with kernel development, reaching up to ~ 7.96 mg/g oil in kernel mature stage. The content of tocopherol was much higher than that of tocotrienol and γ-tocopherol became the dominant component. E. mollis kernel extracts had relatively strong antioxidant capacity. We identified 17 genes (16 VTEs and 1 homogentisic acid geranylgeranyl transferase (HGGT)) directly involving in VitE biosynthesis in RNA-Seq data. Phylogenetic and qRT-PCR results indicated that the annotation and reliability of the RNA-Seq were accurate. Transient overexpression of EmVTE3 and EmWRKY13 in tobacoo leaves increased and decreased the VitE content to 192.18 and 118.29 µg/g, respectively. Weighted gene co-expression analysis elucidated that the blue module showed significant correlation with tocopherol content. Co-expression network analysis revealed that 2-methyl-6-phytobenzoquinone methyltransferase (MPBQ-MT/VTE3) played a vital role and EmWRKY13 may be a key negative regulator in E. mollis VitE biosynthesis. This study not only revealed the traditional VitE biosynthesis pathway in E. mollis, but also set a solid foundation for future genetic breeding of this species.

RNA-Seq and WGCNA Analyses Reveal Key Regulatory Modules and Genes for Salt Tolerance in Cotton

The problem of soil salinization has seriously hindered agricultural development. Cotton is a pioneering salinity-tolerant crop, so harvesting its key salinity-tolerant genes is important for improving crop salt tolerance. In this study, we analyzed changes in the transcriptome expression profiles of the salt-tolerant cultivar Lu Mian 28 (LM) and the salt-sensitive cultivar Zhong Mian Suo 12 (ZMS) after applying salt stress, and we constructed weighted gene co-expression networks (WGCNA). The results indicated that photosynthesis, amino acid biosynthesis, membrane lipid remodeling, autophagy, and ROS scavenging are key pathways in the salt stress response. Plant-pathogen interactions, plant hormone signal transduction, the mitogen-activated protein kinase (MAPK) signaling pathway, and carotenoid biosynthesis are the regulatory networks associated with these metabolic pathways that confer cotton salt tolerance. The gene-weighted co-expression network was used to screen four modules closely related to traits, identifying 114 transcription factors, including WRKYs, ERFs, NACs, bHLHs, bZIPs, and MYBs, and 11 hub genes. This study provides a reference for acquiring salt-tolerant cotton and abundant genetic resources for molecular breeding.

Salicylic acid modulates secondary metabolism and enhanced colchicine accumulation in long yellow daylily (Hemerocallis citrina)

Salicylic acid (SA) is an essential phytoregulator that is widely used to promote the synthesis of high-value nutraceuticals in plants. However, its application in daylily, an ornamental plant highly valued in traditional Chinese medicine, has not been reported. Herein, we investigated the exogenous SA-induced physiological, transcriptional and biochemical changes in long yellow daylily (LYD). We found that 2 mg/L foliar SA treatment significantly improved LYD plant growth and yield. Transcriptome sequencing and differentially expressed genes (DEGs) analysis revealed that the phenylpropanoid biosynthesis, isoquinoline alkaloid biosynthesis, sulfur metabolism, plant hormone signal transduction and tyrosine metabolism were significantly induced in SA-treated leaves. Many transcription factors and antioxidant system-related DEGs were induced under the SA treatment. Biochemical analyses showed that the leaf contents of soluble sugar, soluble protein (Cpr), ascorbic acid (AsA) and colchicine were significantly increased by 15.15% (from 30.16 ± 1.301 to 34.73 ± 0.861 mg/g), 19.54% (from 60.3 ± 2.227 to 72.08 ± 1.617 mg/g), 30.45% (from 190.1 ± 4.56 to 247.98 ± 11.652 μg/g) and 73.05% (from 3.08 ± 0.157 to 5.33 ± 0.462 μg/g), respectively, under the SA treatment. Furthermore, we identified 15 potential candidate genes for enhancing the growth, production and phytochemical content of LYD. Our results provide support for the bioaccumulation of colchicine in yellow daylily and valuable resources for biotechnological-assisted production of this important nutraceutical in Hemerocallis spp.

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.

De novo sequencing, assembly, and characterization of Asparagus racemosus transcriptome and analysis of expression profile of genes involved in the flavonoid biosynthesis pathway

Asparagus racemosus is known for its diverse content of secondary metabolites, i.e., saponins, alkaloids, and a wide range of flavonoids. Flavonoids, including phenols and polyphenols, have a significant role in plant physiology and are synthesized in several tissues. Despite the diverse role of flavonoids, genetic information is limited for flavonoid biosynthesis pathways in A. racemosus. The current study explores full-scale functional genomics information of A. racemosus by de novo transcriptome sequencing using Illumina paired-end sequencing technology to elucidate the genes involved in flavonoid biosynthesis pathways. The de novo assembly of high-quality paired-end reads resulted in ∼2.3 million high-quality reads with a pooled transcript of 45,647 comprising ∼76 Mb transcriptome with a mean length (bp) of 1,674 and N50 of 1,868bp. Furthermore, the coding sequence (CDS) prediction analysis from 45,647 pooled transcripts resulted in 45,444 CDS with a total length and mean length of 76,398,686 and 1,674, respectively. The Gene Ontology (GO) analysis resulted in a high number of CDSs assigned to 25,342 GO terms, which grouped the predicted CDS into three main domains, i.e., Biological Process (19,550), Molecular Function (19,873), and Cellular Component (14,577). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database was used to categorize 6,353 CDS into 25 distinct biological pathway categories, in which the majority of mapped CDS were shown to be related to translation (645), followed by signal transduction (532), carbohydrate metabolism (524), folding, sorting, and degradation (522). Among these, only ∼64 and 14 CDSs were found to be involved in the phenylpropanoid and flavonoid biosynthesis pathways, respectively. Quantitative Real-time PCR was used to check the expression profile of fourteen potential flavonoid biosynthesis pathway genes. The qRT-PCR analysis result matches the transcriptome sequence data validating the Illumina sequence results. Moreover, a large number of genes associated with the flavonoids biosynthesis pathway were found to be upregulated under the induction of methyl jasmonate. The present-day study on transcriptome sequence data of A. racemosus can be utilized for characterizing genes involved in flavonoid biosynthesis pathways and for functional genomics analysis in A. racemosus using the reverse genetics approach (CRISPR/Cas9 technology).

Comparative Transcriptome Profiling Reveals Two WRKY Transcription Factors Positively Regulating Polysaccharide Biosynthesis in Polygonatum cyrtonema

Polygonatum cyrtonema (P. cyrtonema) is a valuable rhizome-propagating traditional Chinese medical herb. Polysaccharides (PCPs) are the major bioactive constituents in P. cyrtonema. However, the molecular basis of PCP biosynthesis in P. cyrtonema remains unknown. In this study, we measured the PCP contents of 11 wild P. cyrtonema germplasms. The results showed that PCP content was the highest in Lishui Qingyuan (LSQY, 11.84%) and the lowest in Hangzhou Lin'an (HZLA, 7.18%). We next analyzed the transcriptome profiles of LSQY and HZLA. Through a qRT-PCR analysis of five differential expression genes from the PCP biosynthesis pathway, phosphomannomutase, UDP-glucose 4-epimerase (galE), and GDP-mannose 4,6-dehydratase were determined as the key enzymes. A protein of a key gene, galE1, was localized in the chloroplast. The PCP content in the transiently overexpressed galE1 tobacco leaves was higher than in the wild type. Moreover, luciferase and Y1H assays indicated that PcWRKY31 and PcWRKY34 could activate galE1 by binding to its promoter. Our research uncovers the novel regulatory mechanism of PCP biosynthesis in P. cyrtonema and is critical to molecular-assisted breeding.

Advances in the Research on Plant WRKY Transcription Factors Responsive to External Stresses

The WRKY transcription factors are a class of transcriptional regulators that are ubiquitous in plants, wherein they play key roles in various physiological activities, including responses to stress. Specifically, WRKY transcription factors mediate plant responses to biotic and abiotic stresses through the binding of their conserved domain to the W-box element of the target gene promoter and the subsequent activation or inhibition of transcription (self-regulation or cross-regulation). In this review, the progress in the research on the regulatory effects of WRKY transcription factors on plant responses to external stresses is summarized, with a particular focus on the structural characteristics, classifications, biological functions, effects on plant secondary metabolism, regulatory networks, and other aspects of WRKY transcription factors. Future research and prospects in this field are also proposed.

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.

California poppy (Eschscholzia californica), the Papaveraceae golden girl model organism for evodevo and specialized metabolism

California poppy or golden poppy (Eschscholzia californica) is the iconic state flower of California, with native ranges from Northern California to Southwestern Mexico. It grows well as an ornamental plant in Mediterranean climates, but it might be invasive in many parts of the world. California poppy was also highly prized by Native Americans for its medicinal value, mainly due to its various specialized metabolites, especially benzylisoquinoline alkaloids (BIAs). As a member of the Ranunculales, the sister lineage of core eudicots it occupies an interesting phylogenetic position. California poppy has a short-lived life cycle but can be maintained as a perennial. It has a comparatively simple floral and vegetative morphology. Several genetic resources, including options for genetic manipulation and a draft genome sequence have been established already with many more to come. Efficient cell and tissue culture protocols are established to study secondary metabolite biosynthesis and its regulation. Here, we review the use of California poppy as a model organism for plant genetics, with particular emphasis on the evolution of development and BIA biosynthesis. In the future, California poppy may serve as a model organism to combine two formerly separated lines of research: the regulation of morphogenesis and the regulation of secondary metabolism. This can provide insights into how these two integral aspects of plant biology interact with each other.

PeCLH2 Gene Positively Regulate Salt Tolerance in Transgenic Populus alba × Populus glandulosa

Salt is an important environmental stress factor, which seriously affects the growth, development and distribution of plants. Chlorophyllase plays an important role in stress response. Nevertheless, little is known about the physiological and molecular mechanism of chlorophyll (Chlase, CLH) genes in plants. We cloned PeCLH2 from Populus euphratica and found that PeCLH2 was differentially expressed in different tissues, especially in the leaves of P. euphratica. To further study the role of PeCLH2 in salt tolerance, PeCLH2 overexpression and RNA interference transgenic lines were established in Populus alba × Populus glandulosa and used for salt stress treatment and physiologic indexes studies. Overexpressing lines significantly improved tolerance to salt treatment and reduced reactive oxygen species production. RNA interference lines showed the opposite. Transcriptome analysis was performed on leaves of control and transgenic lines under normal growth conditions and salt stress to predict genes regulated during salt stress. This provides a basis for elucidating the molecular regulation mechanism of PeCLH2 in response to salt stress and improving the tolerance of poplar under salt stress.

Transcriptome and metabolome profiling of the medicinal plant Veratrum mengtzeanum reveal key components of the alkaloid biosynthesis

Veratrum mengtzeanum is the main ingredient for Chinese folk medicine known as "Pimacao" due to its unique alkaloids. A diverse class of plant-specific metabolites having key pharmacological activities. There are limited studies on alkaloid synthesis and its metabolic pathways in plants. To elucidate the alkaloid pathway and identify novel biosynthetic enzymes and compounds in V. mengtzeanum, transcriptome and metabolome profiling has been conducted in leaves and roots. The transcriptome of V. mengtzeanum leaves and roots yielded 190,161 unigenes, of which 33,942 genes expressed differentially (DEGs) in both tissues. Three enriched regulatory pathways (isoquinoline alkaloid biosynthesis, indole alkaloid biosynthesis and tropane, piperidine and pyridine alkaloid biosynthesis) and a considerable number of genes such as AED3-like, A4U43, 21 kDa protein-like, 3-O-glycotransferase 2-like, AtDIR19, MST4, CASP-like protein 1D1 were discovered in association with the biosynthesis of alkaloids in leaves and roots. Some transcription factor families, i.e., AP2/ERF, GRAS, NAC, bHLH, MYB-related, C3H, FARI, WRKY, HB-HD-ZIP, C2H2, and bZIP were also found to have a prominent role in regulating the synthesis of alkaloids and steroidal alkaloids in the leaves and roots of V. mengtzeanum. The metabolome analysis revealed 74 significantly accumulated metabolites, with 55 differentially accumulated in leaves compared to root tissues. Out of 74 metabolites, 18 alkaloids were highly accumulated in the roots. A novel alkaloid compound viz; 3-Vanilloylygadenine was discovered in root samples. Conjoint analysis of transcriptome and metabolome studies has also highlighted potential genes involved in regulation and transport of alkaloid compounds. Here, we have presented a comprehensive metabolic and transcriptome profiling of V. mengtzeanum tissues. In earlier reports, only the roots were reported as a rich source of alkaloid biosynthesis, but the current findings revealed both leaves and roots as significant manufacturing factories for alkaloid biosynthesis.

A lesion-mimic mutant of Catharanthus roseus accumulates the opioid agonist, akuammicine

Catharanthus roseus is a medicinal plant that produces an abundance of monoterpenoid indole alkaloids (MIAs), notably including the anticancer compounds vinblastine and vincristine. While the canonical pathway leading to these drugs has been resolved, the regulatory and catalytic mechanisms controlling many lateral branches of MIA biosynthesis remain largely unknown. Here, we describe an ethyl methanesulfonate (EMS) C. roseus mutant (M₂-117523) that accumulates high levels of MIAs. The mutant exhibited stunted growth, partially chlorotic leaves, with deficiencies in chlorophyll biosynthesis, and a lesion-mimic phenotype. The lesions were sporadic and spontaneous, appearing after the first true bifoliate and continuing throughout development. The lesions are also the site of high concentrations of akuammicine, a minor constituent of wild type C. roseus leaves. In addition to akuammicine, the lesions were enriched in 25 other MIAs, resulting, in part, from a higher metabolic flux through the pathway. The unique metabolic shift was associated with significant upregulation of biosynthetic and regulatory genes involved in the MIA pathway, including the transcription factors WRKY1, CrMYC2, and ORCA2, and the biosynthetic genes STR, GO, and Redox1. Following the lesion-mimic mutant (LMM) phenotype, the accumulation of akuammicine is jasmonate (JA)-inducible, suggesting a role in plant defence response. Akuammicine is medicinally significant, as a weak opioid agonist, with a preference for the κ-opioid receptor, and a potential anti-diabetic. Further study of akuammicine biosynthesis and regulation can guide plant and heterologous engineering for medicinal uses.

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.

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.

Transcriptomic analysis of genes related to alkaloid biosynthesis and the regulation mechanism under precursor and methyl jasmonate treatment in Dendrobium officinale

Dendrobium officinale is both a traditional herbal medicine and a plant of high ornamental and medicinal value. Alkaloids, especially terpenoid indole alkaloids (TIAs), with pharmacological activities are present in the tissues of D. officinale. A number of genes involved in alkaloid biosynthetic pathways have been identified. However, the regulatory mechanisms underlying the precursor and methyl jasmonate (MeJA)-induced accumulation of alkaloids in D. officinale are poorly understood. In this study, we collected D. officinale protocorm-like bodies (PLBs) and treated them with TIA precursors (tryptophan and secologanin) and MeJA for 0 (T0), 4 (T4) and 24 h (T24); we also established control samples (C4 and C24). Then, we measured the total alkaloid content of the PLBs and performed transcriptome sequencing using the Illumina HiSeq 2,500 system. The total alkaloid content increased significantly after 4 h of treatment. Go and KEGG analysis suggested that genes from the TIA, isoquinoline alkaloid, tropane alkaloid and jasmonate (JA) biosynthetic pathways were significantly enriched. Weighted gene coexpression network analysis (WGCNA) uncovered brown module related to alkaloid content. Six and seven genes related to alkaloid and JA bisosynthetic pathways, respectively, might encode the key enzymes involved in alkaloid biosynthesis of D. officinale. Moreover, 13 transcription factors (TFs), which mostly belong to AP2/ERF, WRKY, and MYB gene families, were predicted to regulate alkaloid biosynthesis. Our data provide insight for studying the regulatory mechanism underlying TIA precursor and MeJA-induced accumulation of three types of alkaloids in D. officinale.

Transcriptome Analysis of Stephania tetrandra and Characterization of Norcoclaurine-6-O-Methyltransferase Involved in Benzylisoquinoline Alkaloid Biosynthesis

Stephania tetrandra (S. Moore) is a source of traditional Chinese medicine that is widely used to treat rheumatism, rheumatoid arthritis, edema, and hypertension. Benzylisoquinoline alkaloids (BIAs) are the main bioactive compounds. However, the current understanding of the biosynthesis of BIAs in S. tetrandra is poor. Metabolite and transcriptomic analyses of the stem, leaf, xylem, and epidermis of S. tetrandra were performed to identify candidate genes associated with BIAs biosynthesis. According to the metabolite analysis, the majority of the BIAs accumulated in the root, especially in the epidermis. Transcriptome sequencing revealed a total of 113,338 unigenes that were generated by de novo assembly. Among them, 79,638 unigenes were successfully annotated, and 42 candidate structural genes associated with 15 steps of BIA biosynthesis identified. Additionally, a new (S)-norcoclaurine-6-O-methyltransferase (6OMT) gene was identified in S. tetrandra, named St6OMT2. Recombinant St6OMT2 catalyzed (S)-norcoclaurine methylation to form (S)-coclaurine in vitro. Maximum activity of St6OMT2 was determined at 30°C and pH 6.0 in NaAc-HAc buffer. Its half-life at 50°C was 22 min with the K_m and k_cat of 28.2 μM and 1.5 s^-1, respectively. Our results provide crucial transcriptome information for S. tetrandra, shedding light on the understanding of BIAs biosynthesis and further gene functional characterization.

Capsicum Leaves under Stress: Using Multi-Omics Analysis to Detect Abiotic Stress Network of Secondary Metabolism in Two Species

The plant kingdom contains an enormous diversity of bioactive compounds which regulate plant growth and defends against biotic and abiotic stress. Some of these compounds, like flavonoids, have properties which are health supporting and relevant for industrial use. Many of these valuable compounds are synthesized in various pepper (Capsicum sp.) tissues. Further, a huge amount of biomass residual remains from pepper production after harvest, which provides an important opportunity to extract these metabolites and optimize the utilization of crops. Moreover, abiotic stresses induce the synthesis of such metabolites as a defense mechanism. Two different Capsicum species were therefore exposed to chilling temperature (24/18 ℃ vs. 18/12 ℃), to salinity (200 mM NaCl), or a combination thereof for 1, 7 and 14 days to investigate the effect of these stresses on the metabolome and transcriptome profiles of their leaves. Both profiles in both species responded to all stresses with an increase over time. All stresses resulted in repression of photosynthesis genes. Stress involving chilling temperature induced secondary metabolism whereas stresses involving salt repressed cell wall modification and solute transport. The metabolome analysis annotated putatively many health stimulating flavonoids (apigetrin, rutin, kaempferol, luteolin and quercetin) in the Capsicum biomass residuals, which were induced in response to salinity, chilling temperature or a combination thereof, and supported by related structural genes of the secondary metabolism in the network analysis.

Jasmonate-Responsive Transcription Factors NnWRKY70a and NnWRKY70b Positively Regulate Benzylisoquinoline Alkaloid Biosynthesis in Lotus (Nelumbo nucifera)

Lotus (Nelumbo nucifera) is a large aquatic plant that accumulates pharmacologically significant benzylisoquinoline alkaloids (BIAs). However, little is known about their biosynthesis and regulation. Here, we show that the two group III WRKY transcription factors (TFs), NnWRKY70a and NnWRKY70b, positively regulate the BIA biosynthesis in lotus. Both NnWRKY70s are jasmonic acid (JA) responsive, with their expression profiles highly correlated to the BIA concentration and BIA pathway gene expression. A dual-luciferase assay showed that NnWRKY70a could transactivate the NnTYDC promoter, whereas NnWRKY70b could activate promoters of the three BIA structural genes, including NnTYDC, NnCYP80G, and Nn7OMT. In addition, the transient overexpression of NnWRKY70a and NnWRKY70b in lotus petals significantly elevated the BIA alkaloid concentrations. Notably, NnWRKY70b seems to be a stronger BIA biosynthesis regulator, because it dramatically induced more BIA structural gene expressions and BIA accumulation than NnWRKY70a. A yeast two-hybrid assay further revealed that NnWRKY70b physically interacted with NnJAZ1 and two other group III WRKY TFs (NnWRKY53b and NnWRKY70a), suggesting that it may cooperate with the other group III WRKYs to adjust the lotus BIA biosynthesis via the JA-signaling pathway. To illustrate the mechanism underlying NnWRKY70b-mediated BIA regulation in the lotus, a simplified model is proposed. Our study provides useful insights into the regulatory roles of WRKY TFs in the biosynthesis of secondary metabolites.

Identification of WRKY transcription factors involved in regulating the biosynthesis of the anti-cancer drug camptothecin in Ophiorrhiza pumila

Camptothecin is a chemotherapeutic drug widely used to treat various cancers. Ophiorrhiza pumila is an ideal plant model for the study of camptothecin production, with various advantages for studying camptothecin biosynthesis and regulation. The DNA-binding WRKY transcription factors have a key regulatory role in secondary metabolite biosynthesis in plants. However, little is currently known about their involvement in camptothecin biosynthesis in O. pumila. We identified 46 OpWRKY genes unevenly distributed on the 11 chromosomes of O. pumila. Phylogenetic and multiple sequence alignment analyses divided the OpWRKY proteins into three subfamilies. Based on spatial expression and co-expression, we targeted the candidate gene OpWRKY6. Overexpression of OpWRKY6 significantly reduced the accumulation of camptothecin compared with the control. Conversely, camptothecin accumulation increased in OpWRKY6 knockout lines. Further biochemical assays showed that OpWRKY6 negatively regulates camptothecin biosynthesis from both the iridoid and shikimate pathways by directly downregulating the gene expression of OpGES, Op10HGO, Op7DLH, and OpTDC. Our data provide direct evidence for the involvement of WRKYs in the regulation of camptothecin biosynthesis and offer valuable information for enriching the production of camptothecin in plant systems.

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.

Scientometric Analysis of Medicinal and Edible Plant Coptis

Objective: A scientometric analysis to obtain knowledge mapping of Coptis revealed the current research situation, knowledge base and research hotspots in Coptis research.

Methods:Coptis-related documents published from 1987 to 2020 were selected through the Web of Science Core Collection. CiteSpace, VOSviewer and Microsoft Excel were used to construct knowledge maps of the Coptis research field.

Results: A total of 367 documents and their references were analyzed. These papers were primarily published in mainland China (214), followed by Japan (57) and South Korea (52), and they each formed respective cooperation networks. The document co-citation analysis suggested that the identification of Coptis Salisb. species, the production of alkaloids, and the mechanisms of action of these alkaloids formed the knowledge bases in this field. A keyword analysis further revealed that the research hotspots were primarily concentrated in three fields of research involving berberine, Coptis chinensis Franch, and Coptis japonica (Thunb) Makino. Oxidative stress, rat plasma (for the determination of plasma alkaloid contents), and Alzheimer’s disease are recent research hotspots associated with Coptis.

Conclusion:Coptis research was mainly distributed in three countries: China, Japan, and South Korea. Researchers were concerned with the identification of Coptis species, the production of Coptis alkaloids, and the efficacy and pharmacological mechanism of the constituent alkaloids. In addition, the anti-oxidative stress, pharmacokinetics, and Alzheimer’s disease treatment of Coptis are new hotspots in this field. This study provides a reference for Coptis researchers.

Genome-wide transcriptional analysis unveils the molecular basis of organ-specific expression of isosteroidal alkaloids biosynthesis in critically endangered Fritillaria roylei Hook

Fritillaria roylei Hook. is a critically endangered high altitude Himalayan medicinal plant species with rich source of pharmaceutically active structurally diverse steroidal alkaloids. Nevertheless, except few marker compounds, the chemistry of the plant remains unexplored. Therefore, in the current study, transcriptome sequencing efforts were made to elucidate isosteroidal alkaloids biosynthesis by creating first organ-specific genomic resource using bulb, stem, and leaf tissues derived from natural populations of Indian Himalayan region. Overall, 349.9 million high quality paired-end reads obtained using NovaSeq 6000 platform were assembled (de novo) into 82,848 unigenes and 31,061 isoforms. Functional annotation and organ specific differential expression (DE) analysis identified 2488 significant DE transcripts, grouped into three potential sub-clusters (sub-cluster I: 728 transcripts; sub-cluster II: 446 transcripts and Sub-cluster III: 1314 transcripts). Subsequently, pathway enrichment (GO, KEGG) and protein-protein network analysis revealed significantly higher enrichment of phenyl-propanoid and steroid backbone including terpenoid, sesquiterpenoid and triterpenoid biosynthesis in bulb. Additionally, upregulated expression of cytochrome P450, UDP-dependent Glucuronosyltransferase families and key transcription factor families (FAR1, bHLH, GRAS, C2H2, TCP and MYB) suggests 'bulb' as a primary site of MVA mediated isosteroidal alkaloids biosynthesis. The comprehensive elucidation of molecular insights in this study is a first step towards the understanding of isosteroidal alkaloid biosynthesis pathway in F. roylei. Furthermore, key genes and regulators identified here can facilitate metabolic engineering of potential bioactive compounds at industrial scale.

Dual function of VvWRKY18 transcription factor in the β-aminobutyric acid-activated priming defense in grapes

Induction of phytoalexin production after invading pathogens is recognized as an essential aspect of the plant-induced resistance. The WRKY family includes plant-specific transcriptional factors associated with plant defense responses, but the comprehensive mechanisms are poorly understood. Here, we attempted to elaborate the regulatory function of VvWRKY18 from the group IIa of WRKY transcription factor (TF) from Vitis vinifera, in the regulation of β-aminobutyric acid (BABA)-activated stilbene phytoalexins biosynthesis and PATHOGENESIS-RELATED (PR) genes expressions in grapes. BABA at 10 mmol L^-1 triggered a priming protection in grapes and conferred a potentiation of the expression levels of VvWRKY18, VvNPR1, and several salicylic acid (SA)-responsive genes, which was accompanied by enhanced stilbene production upon Botrytis cinerea infection. In addition, a physical interaction between VvWRKY18 and the regulatory protein VvNPR1 was detected in vivo and in vitro by yeast-2-hybrid (Y2H), pull-down and co-immunoprecipitation assay (Co-IP) assays. Furthermore, yeast-1-hybrid (Y1H) and dual-luciferase reporter (DLR) assays indicated that VvWRKY18 activated the transcription of STILBENE SYNTHASE (STS) genes, including VvSTS1 and VvSTS2, by directly binding the W-box elements within the specific promoters and resultantly enhancing stilbene phytoalexins biosynthesis. Further investigation demonstrated that heterologous expression of VvWRKY18 elevated the transcriptions of STS and PR genes, thus contributing to potentiating the defense of transgenic Arabidopsis thaliana plants and resultantly inhibiting B. cinerea invasion. Hence, VvWRKY18 serves as a singular effector involved in the synthesis of stilbene phytoalexins in grapes and its interaction with VvNPR1 provided DNA binding ability required for VvNPR1 to initiate systemic acquired resistance (SAR) defense.

Transcriptional Factor-Mediated Regulation of Active Component Biosynthesis in Medicinal Plants

Plants produce thousands of chemically diverse secondary metabolites, many of which have valuable pharmaceutical properties. There is much interest in the synthesis of these pharmaceuticallyvaluable compounds, including the key enzymes and the transcription factors involved. The function and regulatory mechanism of transcription factors in biotic and abiotic stresses have been studied in depth. However, their regulatory roles in the biosynthesis of bioactive compounds, especially in medicinal plants, have only begun. Here, we review what is currently known about how transcription factors contribute to the synthesis of bioactive compounds (alkaloids, terpenoids, flavonoids, and phenolic acids) in medicinal plants. Recent progress has been made in the cloning and characterization of transcription factors in medicinal plants on the genome scale. So far, several large transcription factors have been identified in MYB, WRKY, bHLH, ZIP, AP2/ERF transcription factors. These transcription factors have been predicted to regulate bioactive compound production. These transcription factors positively or negatively regulate the expression of multiple genes encoding key enzymes, and thereby control the metabolic flow through the biosynthetic pathway. Although the research addressing this niche topic is in its infancy, significant progress has been made, and advances in high-throughput sequencing technology are expected to accelerate the discovery of key regulatory transcription factors in medicinal plants. This review is likely to be useful for those interested in the synthesis of pharmaceutically- valuable plant compounds, especially those aiming to breed or engineer plants that produce greater yields of these compounds.

Genome-Wide Profiling of WRKY Genes Involved in Benzylisoquinoline Alkaloid Biosynthesis in California Poppy (Eschscholzia californica)

Transcription factors of the WRKY family play pivotal roles in plant defense responses, including the biosynthesis of specialized metabolites. Based on the previous findings of WRKY proteins regulating benzylisoquinoline alkaloid (BIA) biosynthesis, such as CjWRKY1-a regulator of berberine biosynthesis in Coptis japonica-and PsWRKY1-a regulator of morphine biosynthesis in Papaver somniferum-we performed genome-wide characterization of the WRKY gene family in Eschscholzia californica (California poppy), which produces various BIAs. Fifty WRKY genes were identified by homology search and classified into three groups based on phylogenetic, gene structure, and conserved motif analyses. RNA sequencing showed that several EcWRKY genes transiently responded to methyl jasmonate, a known alkaloid inducer, and the expression patterns of these EcWRKY genes were rather similar to those of BIA biosynthetic enzyme genes. Furthermore, tissue expression profiling suggested the involvement of a few subgroup IIc EcWRKYs in the regulation of BIA biosynthesis. Transactivation analysis using luciferase reporter genes harboring the promoters of biosynthetic enzyme genes indicated little activity of subgroup IIc EcWRKYs, suggesting that the transcriptional network of BIA biosynthesis constitutes multiple members. Finally, we investigated the coexpression patterns of EcWRKYs with some transporter genes and discussed the diversified functions of WRKY genes based on a previous finding that CjWRKY1 overexpression in California poppy cells enhanced BIA secretion into the medium.

Analysis of the Coptis chinensis genome reveals the diversification of protoberberine-type alkaloids

Chinese goldthread (Coptis chinensis Franch.), a member of the Ranunculales, represents an important early-diverging eudicot lineage with diverse medicinal applications. Here, we present a high-quality chromosome-scale genome assembly and annotation of C. chinensis. Phylogenetic and comparative genomic analyses reveal the phylogenetic placement of this species and identify a single round of ancient whole-genome duplication (WGD) shared by the Ranunculaceae. We characterize genes involved in the biosynthesis of protoberberine-type alkaloids in C. chinensis. In particular, local genomic tandem duplications contribute to member amplification of a Ranunculales clade-specific gene family of the cytochrome P450 (CYP) 719. The functional versatility of a key CYP719 gene that encodes the (S)-canadine synthase enzyme involved in the berberine biosynthesis pathway may play critical roles in the diversification of other berberine-related alkaloids in C. chinensis. Our study provides insights into the genomic landscape of early-diverging eudicots and provides a valuable model genome for genetic and applied studies of Ranunculales.

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.

Genome-Wide Identification of WRKY Genes in Artemisia annua: Characterization of a Putative Ortholog of AtWRKY40

Artemisia annua L. is well-known as the plant source of artemisinin, a sesquiterpene lactone with effective antimalarial activity. Here, a putative ortholog of the Arabidopsis thaliana WRKY40 transcription factor (TF) was isolated via reverse transcription-polymerase chain reaction and rapid amplification of cDNA ends in A. annua and named AaWRKY40. A putative nuclear localization domain was identified in silico and experimentally confirmed by using protoplasts of A. annua transiently transformed with AaWRKY40-GFP. A genome-wide analysis identified 122 WRKY genes in A. annua, and a manually curated database was obtained. The deduced proteins were categorized into the major WRKY groups, with group IIa containing eight WRKY members including AaWRKY40. Protein motifs, gene structure, and promoter regions of group IIa WRKY TFs of A. annua were characterized. The promoter region of AaWRKY group IIa genes contained several abiotic stress cis-acting regulatory elements, among which a highly conserved W-box motif was identified. Expression analysis of AaWRKY40 compared to AaWRKY1 in A. annua cell cultures treated with methyl jasmonate known to enhance artemisinin production, suggested a possible involvement of AaWRKY40 in terpenoid metabolism. Further investigation is necessary to study the role of AaWRKY40 and possible interactions with other TFs in A. annua.

Comparative Transcriptome Analysis Reveals Candidate Genes Involved in Isoquinoline Alkaloid Biosynthesis in Stephania tetrandra

The roots of Stephania tetrandra are used as a traditional Chinese medicine. Isoquinoline alkaloids are considered to be the most important and effective components in this herb, but little is known about the molecular mechanism underlying their biosynthesis. In this context, this study aimed to reveal candidate genes related to isoquinoline alkaloid biosynthesis in S. tetrandra. Determination of tetrandrine and fangchinoline in the roots and leaves of S. tetrandra by HPLC showed that the roots had much higher contents of the two isoquinoline alkaloids than the leaves. Thus, a comparative transcriptome analysis of the two tissues was performed to uncover candidate genes involved in isoquinoline alkaloid biosynthesis. A total of 71 674 unigenes was obtained and 31 994 of these were assigned putative functions based on BLAST searches against 6 annotation databases. Among the 79 isoquinoline alkaloid-related unigenes, 51 were differentially expressed, with 42 and 9 genes upregulated and downregulated, respectively, when the roots were compared with the leaves. The upregulated differentially expressed genes were consistent with isoquinoline alkaloid accumulation in roots and thus were deemed key candidate genes for isoquinoline alkaloid biosynthesis in the roots. Moreover, the expression profiles of 10 isoquinoline alkaloid-related differentially expressed genes between roots and leaves were validated by quantitative real-time polymerase chain reaction, which indicated that our transcriptome and gene expression profiles were reliable. This study not only provides a valuable genomic resource for S. tetrandra but also proposes candidate genes involved in isoquinoline alkaloid biosynthesis and transcription factors related to the regulation of isoquinoline alkaloid biosynthesis. The results lay a foundation for further studies on isoquinoline alkaloid biosynthesis in this medicinal plant.

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.

Full-length transcriptome analysis of Coptis deltoidea and identification of putative genes involved in benzylisoquinoline alkaloids biosynthesis based on combined sequencing platforms

The study carry out comprehensive transcriptome analysis of C. deltoidea and exploration of BIAs biosynthesis and accumulation based on UHPLC-MS/MS and combined sequencing platforms. Coptis deltoidea is an important medicinal plant with a long history of medicinal use, which is rich in benzylisoquinoline alkaloids (BIAs). In this study, Ultra performance liquid chromatography-electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) and combined sequencing platforms were performed for exploration of BIAs biosynthesis, accumulation and comprehensive transcriptome analysis of C. deltoidea. By metabolism profiling, the accumulation of ten BIAs was analyzed using UHPLC-MS/MS and different contents were observed in different organs. From transcriptome sequencing result, we applied single-molecule real-time (SMRT) sequencing to C. deltoidea and generated a total of 75,438 full-length transcripts. We proposed the candidate biosynthetic pathway of tyrosine, precursor of BIAs, and identified 64 full length-transcripts encoding enzymes putatively involved in BIAs biosynthesis. RNA-Seq data indicated that the majority of genes exhibited relatively high expression level in roots. Transport of BIAs was also important for their accumulation. Here, 9 ABC transporters and 2 MATE transporters highly homologous to known alkaloid transporters related with BIAs transport in roots and rhizomes were identified. These findings based on the combined sequencing platforms provide valuable genetic information for C. deltoidea and the results of transcriptome combined with metabolome analysis can help us better understand BIAs biosynthesis and transport in this medicinal plant. The information will be critical for further characterization of C. deltoidea transcriptome and molecular-assisted breeding for this medicinal plant with scarce resources.

Fusarium oxysporum infection activates the plastidial branch of the terpenoid biosynthesis pathway in flax, leading to increased ABA synthesis

Upregulation of the terpenoid pathway and increased ABA content in flax upon Fusarium infection leads to activation of the early plant's response (PR genes, cell wall remodeling, and redox status). Plants have developed a number of defense strategies against the adverse effects of fungi such as Fusarium oxysporum. One such defense is the production of antioxidant secondary metabolites, which fall into two main groups: the phenylpropanoids and the terpenoids. While functions and biosynthesis of phenylpropanoids have been extensively studied, very little is known about the genes controlling the terpenoid synthesis pathway in flax. They can serve as antioxidants, but are also substrates for a plethora of different compounds, including those of regulatory functions, like ABA. ABA's function during pathogen attack remains obscure and often depends on the specific plant-pathogen interactions. In our study we showed that in flax the non-mevalonate pathway is strongly activated in the early hours of pathogen infection and that there is a redirection of metabolites towards ABA synthesis. The elevated synthesis of ABA correlates with flax resistance to F. oxysporum, thus we suggest ABA to be a positive regulator of the plant's early response to the infection.

Comprehensive Analysis and Functional Studies of WRKY Transcription Factors in Nelumbo nucifera

The WRKY family is one of the largest transcription factor (TF) families in plants and plays central roles in modulating plant stress responses and developmental processes, as well as secondary metabolic regulations. Lotus (Nelumbo nucifera) is an aquatic crop that has significant food, ornamental and pharmacological values. Here, we performed an overview analysis of WRKY TF family members in lotus, and studied their functions in environmental adaptation and regulation of lotus benzylisoquinoline alkaloid (BIA) biosynthesis. A total of 65 WRKY genes were identified in the lotus genome and they were well clustered in a similar pattern with their Arabidopsis homologs in seven groups (designated I, IIa-IIe, and III), although no lotus WRKY was clustered in the group IIIa. Most lotus WRKYs were functionally paired, which was attributed to the recently occurred whole genome duplication in lotus. In addition, lotus WRKYs were regulated dramatically by salicilic acid (SA), jasmonic acid (JA), and submergence treatments, and two lotus WRKYs, NnWRKY40a and NnWRKY40b, were significantly induced by JA and promoted lotus BIA biosynthesis through activating BIA biosynthetic genes. The investigation of WRKY TFs for this basal eudicot reveals new insights into the evolution of the WRKY family, and provides fundamental information for their functional studies and lotus breeding.

Hybrid sequencing of the Gynostemma pentaphyllum transcriptome provides new insights into gypenoside biosynthesis

Background

Gypenosides are a group of triterpene saponins from Gynostemma pentaphyllum that are the same as or very similar to ginsenosides from the Panax species. Several enzymes involved in ginsenoside biosynthesis have been characterized, which provide important clues for elucidating the gypenoside biosynthetic pathway. We suppose that gypenosides and ginsenosides may have a similar biosynthetic mechanism and that the corresponding enzymes in the two pathways may have considerable similarity in their sequences. To further understand gypenoside biosynthesis, we sequenced the G. pentaphyllum transcriptome with a hybrid sequencing-based strategy and then determined the candidate genes involved in this pathway using phylogenetic tree construction and gene expression analysis.

Results

Following the PacBio standard analysis pipeline, 66,046 polished consensus sequences were obtained, while Illumina data were assembled into 140,601 unigenes with Trinity software. Then, these output sequences from the two analytical routes were merged. After removing redundant data with CD-HIT software, a total of 140,157 final unigenes were obtained. After functional annotation, five 2,3-oxidosqualene cyclase genes, 145 cytochrome P450 genes and 254 UDP-glycosyltransferase genes were selected for the screening of genes involved in gypenoside biosynthesis. Using phylogenetic analysis, several genes were divided into the same subfamilies or closely related evolutionary branches with characterized enzymes involved in ginsenoside biosynthesis. Using real-time PCR technology, their expression patterns were investigated in different tissues and at different times after methyl jasmonate induction. Since the genes in the same biosynthetic pathway are generally coexpressed, we speculated that GpOSC1, GpCYP89, and GpUGT35 were the leading candidates for gypenoside biosynthesis. In addition, six GpWRKYs and one GpbHLH might play a possible role in regulating gypenoside biosynthesis.

Conclusions

We developed a hybrid sequencing strategy to obtain longer length transcriptomes with increased accuracy, which will greatly contribute to downstream gene screening and characterization, thus improving our ability to elucidate secondary metabolite biosynthetic pathways. With this strategy, we found several candidate genes that may be involved in gypenoside biosynthesis, which laid an important foundation for the elucidation of this biosynthetic pathway, thus greatly contributing to further research in metabolic regulation, synthetic biology and molecular breeding in this species.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6000-y) contains supplementary material, which is available to authorized users.

Characterization of LuWRKY36, a flax transcription factor promoting secoisolariciresinol biosynthesis in response to Fusarium oxysporum elicitors in Linum usitatissimum L. hairy roots

The involvement of a WRKY transcription factor in the regulation of lignan biosynthesis in flax using a hairy root system is described. Secoisolariciresinol is the main flax lignan synthesized by action of LuPLR1 (pinoresinol-lariciresinol reductase 1). LuPLR1 gene promoter deletion experiments have revealed a promoter region containing W boxes potentially responsible for the response to Fusarium oxysporum. W boxes are bound by WRKY transcription factors that play a role in the response to stress. A candidate WRKY transcription factor, LuWRKY36, was isolated from both abscisic acid and Fusarium elicitor-treated flax cell cDNA libraries. This transcription factors contains two WRKY DNA-binding domains and is a homolog of AtWRKY33. Different approaches confirmed LuWRKY36 binding to a W box located in the LuPLR1 promoter occurring through a unique direct interaction mediated by its N-terminal WRKY domain. Our results propose that the positive regulator action of LuWRKY36 on the LuPLR1 gene regulation and lignan biosynthesis in response to biotic stress is positively mediated by abscisic acid and inhibited by ethylene. Additionally, we demonstrate a differential Fusarium elicitor response in susceptible and resistant flax cultivars, seen as a faster and stronger LuPLR1 gene expression response accompanied with higher secoisolariciresinol accumulation in HR of the resistant cultivar.

Comparative transcriptome analysis revealing the potential mechanism of seed germination stimulated by exogenous gibberellin in Fraxinus hupehensis

BACKGROUND

Fraxinus hupehensis is an endangered tree species that is endemic to in China; the species has very high commercial value because of its intricate shape and potential to improve and protect the environment. Its seeds show very low germination rates in natural conditions. Preliminary experiments indicated that gibberellin (GA₃) effectively stimulated the seed germination of F. hupehensis. However, little is known about the physiological and molecular mechanisms underlying the effect of GA₃ on F. hupehensis seed germination.

RESULTS

We compared dormant seeds (CK group) and germinated seeds after treatment with water (W group) and GA₃ (G group) in terms of seed vigor and several other physiological indicators related to germination, hormone content, and transcriptomics. Results showed that GA₃ treatment increases seed vigor, energy requirements, and trans-Zetain (ZT) and GA₃ contents but decreases sugar and abscisic acid (ABA) contents. A total of 116,932 unigenes were obtained from F. hupehensis transcriptome. RNA-seq analysis identified 31,856, 33,188 and 2056 differentially expressed genes (DEGs) between the W and CK groups, the G and CK groups, and the G and W groups, respectively. Up-regulation of eight selected DEGs of the glycolytic pathway accelerated the oxidative decomposition of sugar to release energy for germination. Up-regulated genes involved in ZT (two genes) and GA₃ (one gene) biosynthesis, ABA degradation pathway (one gene), and ABA signal transduction (two genes) may contribute to seed germination. Two down-regulated genes associated with GA₃ signal transduction were also observed in the G group. GA₃-regulated genes may alter hormone levels to facilitate germination. Candidate transcription factors played important roles in GA₃-promoted F. hupehensis seed germination, and Quantitative Real-time PCR (qRT-PCR) analysis verified the expression patterns of these genes.

CONCLUSION

Exogenous GA₃ increased the germination rate, vigor, and water absorption rate of F. hupehensis seeds. Our results provide novel insights into the transcriptional regulation mechanism of effect of exogenous GA₃ on F. hupehensis seed germination. The transcriptome data generated in this study may be used for further molecular research on this unique species.

Comparative analysis of the root and leaf transcriptomes in Chelidonium majus L

Chelidonium majus is a traditional medicinal plant, which commonly known as a rich resource for the major benzylisoquinoline alkaloids (BIAs), including morphine, sanguinarine, and berberine. To understand the biosynthesis of C. majus BIAs, we performed de novo transcriptome sequencing of its leaf and root tissues using Illumina technology. Following comprehensive evaluation of de novo transcriptome assemblies produced with five programs including Trinity, Bridger, BinPacker, IDBA-tran, and Velvet/Oases using a series of k-mer sizes (from 25 to 91), BinPacker was found to produce the best assembly using a k-mer of 25. This study reports the results of differential gene expression (DGE), functional annotation, gene ontology (GO) analysis, classification of transcription factor (TF)s, and SSR and miRNA discovery. Our DGE analysis identified 6,028 transcripts that were up-regulated in the leaf, and 4,722 transcripts that were up-regulated in the root. Further investigations showed that most of the genes involved in the BIA biosynthetic pathway are significantly expressed in the root compared to the leaf. GO analysis showed that the predominant GO domain is "cellular component", while TF analysis found bHLH to be the most highly represented TF family. Our study further identified 10 SSRs, out of a total of 39,841, that showed linkage to five unigenes encoding enzymes in the BIA pathway, and 10 conserved miRNAs that were previously not detected in this plant. The comprehensive transcriptome information presented herein provides a foundation for further explorations on study of the molecular mechanisms of BIA synthesis in C. majus.

Transcriptome analysis reveals important candidate genes involved in grain-size formation at the stage of grain enlargement in common wheat cultivar "Bainong 4199"

Grain-size is one of the yield components, and the first 14 days after pollination (DAP) is a crucial stage for wheat grain-size formation. To understand the mechanism of grain-size formation at the whole gene expression level and to identify the candidate genes related to grain pattern formation, cDNA libraries from immature grains of 5 DAP and 14 DAP were constructed. According to transcriptome analysis, a total of 12,555 new genes and 9,358 differentially expressed genes (DEGs) were obtained. In DEGs, 2,876, 3,357 and 3,125 genes were located on A, B and D subgenome respectively. 9,937 (79.15%) new genes and 9,059 (96.80%) DEGs were successfully annotated. For DEGs, 4,453 were up-regulated and 4,905 were down-regulated at 14 DAP. The Gene Ontology (GO) database indicated that most of the grain-size-related genes were in the same cluster. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database analysis showed that 130, 129 and 20 DEGs were respectively involved in starch and sucrose metabolism, plant hormone signal transduction and ubiquitin-mediated proteolysis. Expression levels of 8 randomly selected genes were confirmed by qRT-PCR, which was consistent with the transcriptome data. The present database will help us understand the molecular mechanisms underlying early grain development and provide the foundation for increasing grain-size and yield in wheat breeding programs.

VvWRKY8 represses stilbene synthase genes through direct interaction with VvMYB14 to control resveratrol biosynthesis in grapevine

Resveratrol (Res) is a stilbenoid, a group of plant phenolic metabolites derived from stilbene that possess activities against pests, pathogens, and abiotic stresses. Only a few species, including grapevine (Vitis), synthesize and accumulate Res. Although stilbene synthases (STSs) have been isolated and characterized in several species, the gene regulatory mechanisms underlying stilbene biosynthesis are still largely unknown. Here, we characterize a grapevine WRKY transcription factor, VvWRKY8, that regulates the Res biosynthetic pathway. Transient and stable overexpression of VvWRKY8 in grapevine results in decreased expression of VvSTS15/21 and VvMYB14, as well as in a reduction of Res accumulation. VvWRKY8 does not bind to or activate the promoters of VvMYB14 and VvSTS15/21; however, it physically interacts with VvMYB14 proteins through their N-terminal domains to prevent them from binding to the VvSTS15/21 promoter. Application of exogenous Res results in the stimulation of VvWRKY8 expression and in a decrease of VvMYB14 and VvSTS15/21 expression in grapevine suspension cells, and in the activation of the VvWRKY8 promoter in tobacco leaves. These results demonstrate that VvWRKY8 represses VvSTS15/21 expression and Res biosynthesis through interaction with VvMYB14. In this context, the VvMYB14-VvSTS15/21-Res-VvWRKY8 regulatory loop may be an important mechanism for the fine-tuning of Res biosynthesis in grapevine.

Transcriptome analysis provides insight into prickle development and its link to defense and secondary metabolism in Solanum viarum Dunal

Prickles are epidermal outgrowth found on the aerial surface of several terrestrial plants. Microscopic studies on prickles of S. viarum Dunal indicated a crucial role of glandular trichomes (GTs) in their development. A spontaneously obtained prickleless mutant showed normal epidermal GTs, but its downstream developmental process to prickle was perturbed. Thus, prickleless mutant offers an ideal opportunity to unveil molecular regulators working downstream to GTs in the prickle formation. Differential transcriptome analysis of epidermis of prickly and prickleless mutant revealed that expression of several defense regulators like ethylene, salicylic acid, PR-proteins, etc. were significantly down-regulated in prickleless mutant, provide an important link between defense and prickle development. It was also noteworthy that the expression of few essential development related TFs like MADS-box, R2R3-MYB, REM, DRL1, were also down-regulated in the stem, petioles, and leaves of prickleless mutant indicating their potential role in prickle development. Interestingly, the gene expression of terpenoid, steroid, flavonoid, glucosinolate, and lignin biosynthesis pathways were up-regulated in prickleless mutant. The biochemical and qRT-PCR analysis also confirmed metabolite elevation. These results indicated that the loss of prickle was compensated by elevated secondary metabolism in the prickleless mutant which played important role in the biotic and abiotic stress management.

Benzylisoquinoline Alkaloids Biosynthesis in Sacred Lotus

Sacred lotus (Nelumbo nucifera Gaertn.) is an ancient aquatic plant used throughout Asia for its nutritional and medicinal properties. Benzylisoquinoline alkaloids (BIAs), mostly within the aporphine and bisbenzylisoquinoline structural categories, are among the main bioactive constituents in the plant. The alkaloids of sacred lotus exhibit promising anti-cancer, anti-arrhythmic, anti-HIV, and anti-malarial properties. Despite their pharmacological significance, BIA metabolism in this non-model plant has not been extensively investigated. In this review, we examine the diversity of BIAs in sacred lotus, with an emphasis on the distinctive stereochemistry of alkaloids found in this species. Additionally, we discuss our current understanding of the biosynthetic genes and enzymes involved in the formation of 1-benzylisoquinoline, aporphine, and bisbenzylisoquinoline alkaloids in the plant. We conclude that a comprehensive functional characterization of alkaloid biosynthetic enzymes using both in vitro and in vivo methods is required to advance our limited knowledge of BIA metabolism in the sacred lotus.

A WRKY transcription factor, PcWRKY33, from Polygonum cuspidatum reduces salt tolerance in transgenic Arabidopsis thaliana

PcWRKY33 is a transcription factor which can reduce salt tolerance by decreasing the expression of stress-related genes and increasing the cellular levels of reactive oxygen species (ROS). WRKY transcription factors play important roles in the regulation of biotic and abiotic stresses. Here, we report a group I WRKY gene from Polygonum cuspidatum, PcWRKY33, that encodes a nucleoprotein, which specifically binds to the W-box in the promoter of target genes to regulate their expression. The results from qPCR and promoter analysis show that expression of PcWRKY33 can be induced by various abiotic stresses, including NaCl and plant hormones. Overexpression of PcWRKY33 in Arabidopsis thaliana reduced tolerance to salt stress. More specifically, several physiological parameters (such as root length, seed germination rate, seedling survival rate, and chlorophyll concentration) of the transgenic lines were significantly lower than those of the wild type under salt stress. In addition, following exposure to salt stress, transgenic plants showed decreased expression of stress-related genes, a weakened ability to maintain Na⁺/K⁺ homeostasis, decreased activities of reactive oxygen species- (ROS-) scavenging enzymes, and increased accumulation of ROS. Taken together, these results suggest that PcWRKY33 negatively regulates the salt tolerance in at least two ways: by down-regulating the induction of stress-related genes and by increasing the level of cellular ROS. In sum, our results indicate that PcWRKY33 is a group I WRKY transcription factor involved in abiotic stress regulation.

Comprehensive metabolic and transcriptomic profiling of various tissues provide insights for saponin biosynthesis in the medicinally important Asparagus racemosus

Asparagus racemosus (Shatavari), belongs to the family Asparagaceae and is known as a “curer of hundred diseases” since ancient time. This plant has been exploited as a food supplement to enhance immune system and regarded as a highly valued medicinal plant in Ayurvedic medicine system for the treatment of various ailments such as gastric ulcers, dyspepsia, cardiovascular diseases, neurodegenerative diseases, cancer, as a galactogogue and against several other diseases. In depth metabolic fingerprinting of various parts of the plant led to the identification of 13 monoterpenoids exclusively present in roots. LC-MS profiling led to the identification of a significant number of steroidal saponins (33). However, we have also identified 16 triterpene saponins for the first time in A. racemosus. In order to understand the molecular basis of biosynthesis of major components, transcriptome sequencing from three different tissues (root, leaf and fruit) was carried out. Functional annotation of A. racemosus transcriptome resulted in the identification of 153 transcripts involved in steroidal saponin biosynthesis, 45 transcripts in triterpene saponin biosynthesis, 44 transcripts in monoterpenoid biosynthesis and 79 transcripts in flavonoid biosynthesis. These findings will pave the way for better understanding of the molecular basis of steroidal saponin, triterpene saponin, monoterpenoids and flavonoid biosynthesis in A. racemosus.

Combinatorial Regulation of Stilbene Synthase Genes by WRKY and MYB Transcription Factors in Grapevine (Vitis vinifera L.)

Stilbene synthase (STS) is the key enzyme leading to the biosynthesis of resveratrol. Recently we reported two R2R3-MYB transcription factor (TF) genes that regulate the stilbene biosynthetic pathway in grapevine: VviMYB14 and VviMYB15. These genes are strongly co-expressed with STS genes under a range of stress and developmental conditions, in agreement with the specific activation of STS promoters by these TFs. Genome-wide gene co-expression analysis using two separate transcriptome compendia based on microarray and RNA sequencing data revealed that WRKY TFs were the top TF family correlated with STS genes. On the basis of correlation frequency, four WRKY genes, namely VviWRKY03, VviWRKY24, VviWRKY43 and VviWRKY53, were further shortlisted and functionally validated. Expression analyses under both unstressed and stressed conditions, together with promoter-luciferase reporter assays, suggested different hierarchies for these TFs in the regulation of the stilbene biosynthetic pathway. In particular, VviWRKY24 seems to act as a singular effector in the activation of the VviSTS29 promoter, while VviWRKY03 acts through a combinatorial effect with VviMYB14, suggesting that these two regulators may interact at the protein level as previously reported in other species.

The effects of promoter variations of the N-Methylcanadine 1-Hydroxylase (CYP82Y1) gene on the noscapine production in opium poppy

Noscapine is an antitumor alkaloid produced in opium poppy (Papaver somniferum) and some members of the Papaveraceae family. It has been primarily used for its antitussive effects; more recently, its anticancer properties were shown. Herein, we detected an SSR embedded in the promoter region of the CYP82Y1 gene, which was found to be the first committed-step enzyme in the noscapine biosynthesis pathway, using the MISA program. Some collected ecotypes of P. somniferum were investigated for understanding of SSRs role in the regulation of gene expression and metabolite content. Quantitative PCR showed that a variation in the motif repeat number (either a decrease or increase) down-regulated the expression of the CYP82Y1 gene. Furthermore, the analysis of noscapine content suggested that a variation in the promoter region influence noscapine amount. Moreover, P. bracteatum was analyzed in both transcript and metabolite levels, and illustrated much less expression and metabolite level in comparison to P. somniferum. By exploiting the transcriptome data from the eight genera of the Papaveraceae family, we found that noscapine biosynthesis genes are present in P. bracteatum and are not shared in other genera of the Papaveraceae family. This results may explain production of a confined metabolite within a genus.

Investigation of benzylisoquinoline alkaloid biosynthetic pathway and its transcriptional regulation in lotus

Lotus predominantly accumulates benzylisoquinoline alkaloids (BIAs), but their biosynthesis and regulation remain unclear. Here, we investigated structural and regulatory genes involved in BIA accumulation in lotus. Two clustered CYP80 genes were identified to be responsible for the biosynthesis of bis-BIAs and aporphine-type BIAs, respectively, and their tissue-specific expression causes divergence in alkaloid component between leaf and embryo. In contrast with the common (S)-reticuline precursor for most BIAs, aporphine alkaloids in lotus leaf may result from the (S)-N-methylcoclaurine precursor. Structural diversity of BIA alkaloids in the leaf is attributed to enzymatic modifications, including intramolecular C-C phenol coupling on ring A and methylation and demethylation at certain positions. Additionally, most BIA biosynthetic pathway genes show higher levels of expression in the leaf of high-BIA cultivar compared with low-BIA cultivar, suggesting transcriptional regulation of BIA accumulation in lotus. Five transcription factors, including three MYBs, one ethylene-responsive factor, and one basic helix-loop-helix (bHLH), were identified to be candidate regulators of BIA biosynthesis in lotus. Our study reveals a BIA biosynthetic pathway and its transcriptional regulation in lotus, which will enable a deeper understanding of BIA biosynthesis in plants.